ko1982 发表于 2006-4-25 17:33

[求助]快疯了。。高人帮忙看一下啊。

解这个式子。。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)
但是好象不行。。。。到底要怎么样才能解出来啊。。。郁闷

happy 发表于 2006-4-25 21:03

回复:(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]
查看完整版本: [求助]快疯了。。高人帮忙看一下啊。