|
楼主 |
发表于 2006-8-31 09:30
|
显示全部楼层
非常感谢您!
子空间迭代法,不太懂这个方法,能说的具体点不?
我在matlab中自己这样做的,可是出问题了function y=vdp_eq1(t,x,flag,k,c)
y=[eye(3)*x(2);
-k*x(1)-c*x(2)];
这是我编写的函数,x0=[0;0;0;0;0;0];k=[1 2 3;4 2 5;6 7 2];c=[3 5 6;4 4 7;3 5 7];[t1,y1]=ode45('vdp_eq1',[0,20],x0,h_opt,k,c),在matlab空间下运行这句话,就出现了这样的错误,??? Error using ==> funfun\private\odearguments
VDP_EQ1 must return a column vector.
Error in ==> ode45 at 171
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...
能给指点一下这是怎么回事啊? |
|