[求助]快疯了。。高人帮忙看一下啊。
解这个式子。。w=pi,t取60到240。。搞的程序如下function xdot=lorenzeq(t,x)
xdot=[0.7*x(2)+10*(0.1-(x(2)^2));
-x+0.25*sin(pi*t)]保存;
t_final=240; x0=;
=ode45('lorenzeq',,x0);plot(t,x)
但是好象不行。。。。到底要怎么样才能解出来啊。。。郁闷
回复:(ko1982)[求助]快疯了。。高人帮忙看一下啊。...
function xdot=lorenzeq(t,x)xdot=[0.7*x(2)+10*(0.1-(x(2)^2));
-x(1)+0.25*sin(pi*t)]保存;
t_final=240; x0=;
=ode45('lorenzeq',,x0);plot(t,x)
另外你仔细看一下你的第一个方程,给得好像不对,应该是
0.7*x(2)+10*x(1)*(0.1-(x(2)^2));
页:
[1]