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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

解这个式子。。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=[0;1];
[t,x]=ode45('lorenzeq',[0,t_final],x0);plot(t,x)
但是好象不行。。。。到底要怎么样才能解出来啊。。。郁闷

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

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=[0;1];
[t,x]=ode45('lorenzeq',[0,t_final],x0);plot(t,x)

另外你仔细看一下你的第一个方程，给得好像不对，应该是
0.7*x(2)+10*x(1)*(0.1-(x(2)^2));

账号		自动登录	找回密码
密码			我要加入

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（ko1982）[求助]快疯了。。高人帮忙看一下啊。...