2015-03-02 103 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.