ode45系统二阶答案

【问题标题】：ode45 system second orderode45系统二阶
【发布时间】：2015-03-02 20:03:33
【问题描述】：

我已经写了这个程序，但它不工作

function xprime = exh(t,z)
xprime = [z(2);
          30*z(3) - 40*z(1); 
          z(4);
          30*z(1)+30*z(5)-60*z(3);
          z(6);
          30*z(3)-40*z(5)];

tspan=[0,200];
Z0=[0.05 0.04 0.03];
[t,z] = ode45('exh',tspan,Z0);
plot(t,z)

【问题讨论】：

你能更详细地描述一下什么不起作用吗？
它给了我一个错误
odearguments 中的错误（第 87 行） f0 = feval(ode,t0,y0,args{:}); % ODE15I 将 args{1} 设置为 yp0。 ode45 错误（第 113 行）[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

标签： system ode differential-equations

【解决方案1】：

你的系统是6维的，你的初始值是3维的，这个不行，维数必须一致。

换句话说，对于 2 阶系统，您不仅需要初始位置，还需要初始速度。

【讨论】：