Question

我对朱莉亚的前途语言并不陌生，希望它能加速我僵硬的常微分方程。这是东西： 1）必须使用尺寸为400x400的Matlab之外的质量，阻尼，刚度矩阵，以矩阵形式定义方程。实现了二阶ODE的公共状态空间表示法。 2）除了线性动力学，还有非线性力作用，这取决于其某些状态。这些力必须在ode函数中定义。

但是，由于初始条件的影响，状态变量根本不会改变，尽管应该改变。这是一个用于较小原型的示例代码，具有较小的矩阵：

#Load packages
using LinearAlgebra
using OrdinaryDiffEq
using DifferentialEquations
using Plots

# Define constant matrices (here generated for the example)
const M=Matrix{Float64}(I, 4, 4) #  Mass matrix
const C=zeros(Float64, 4, 4) #  Damping matrix
const K=[10.0 0.0 0.0 0.0;  0.0 7.0 0.0 0.0; 0.0 0.0 6.0 0.0;0.0 0.0 5.0 0.0] #  Stiffness matrix

x0 = [0.0;0.0;0.0;0.0; 1.0; 1.0; 1.0; 1.0] #  Initial conditions
tspan = (0.,1.0) #  Simulation time span

#Define the underlying equation
function FourDOFoscillator(xdot,x,p,t)
  xdot=[-inv(M)*C -inv(M)*K; Matrix{Float64}(I, 4, 4) zeros(Float64, 4, 4)]*x
end

#Pass to Solvers
prob = ODEProblem(FourDOFoscillator,x0,tspan)
sol = solve(prob,alg_hints=[:stiff],reltol=1e-8,abstol=1e-8)
plot(sol)

我想念什么？谢谢

槟榔

Answer 1

您不会改变输出，而是创建一个新数组。如果您执行xdot.=，它将起作用。

#Load packages
using LinearAlgebra
using OrdinaryDiffEq
using DifferentialEquations
using Plots

# Define constant matrices (here generated for the example)
const M=Matrix{Float64}(I, 4, 4) #  Mass matrix
const C=zeros(Float64, 4, 4) #  Damping matrix
const K=[10.0 0.0 0.0 0.0;  0.0 7.0 0.0 0.0; 0.0 0.0 6.0 0.0;0.0 0.0 5.0 0.0] #  Stiffness matrix

x0 = [0.0;0.0;0.0;0.0; 1.0; 1.0; 1.0; 1.0] #  Initial conditions
tspan = (0.,1.0) #  Simulation time span

#Define the underlying equation
function FourDOFoscillator(xdot,x,p,t)
  xdot.=[-inv(M)*C -inv(M)*K; Matrix{Float64}(I, 4, 4) zeros(Float64, 4, 4)]*x
end

#Pass to Solvers
prob = ODEProblem(FourDOFoscillator,x0,tspan)
sol = solve(prob,alg_hints=[:stiff],reltol=1e-8,abstol=1e-8)
plot(sol)

Julia -matrix形式的ODE可能吗？

1 个答案: