ode45系统二阶

时间:2015-03-02 20:03:34

标签: system ode differential-equations

我已经编写了这个程序,但它无法正常工作

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)

1 个答案:

答案 0 :(得分:1)

您的系统尺寸为6,初始值为尺寸3,这不起作用,尺寸必须相同。

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