2015-03-02 57 views
-2

我有写这个程序,但它不工作ODE45系统二阶

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) 
+0

系统的初始位置,而且初期的速度,你可以描述哪些项目没有更详细一点的工作? – TofferJ 2015-03-02 20:26:12

+0

它给了我一个错误 – 2015-03-03 06:15:00

+0

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,... – 2015-03-03 06:15:09

回答

1

你的系统有尺寸6,你的初始值的维数3,这可不行,尺寸必须是相同的。

换句话说,你需要的不仅是订单2.