ode45系统二阶

Question

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

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)

Answer 1

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

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