大侠们帮看一下，为什么运行不了

xxzz · 发表于 2009-4-13 16:44

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

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

x

function xdot=fun2_4(t,x);
xdot1(1)=36*x(1)-36*x(2);
xdot1(2)=-x(1)*x(3)-3*x(2);
xdot1(3)=x(1)*x(2)-20*x(3);
xdot=xdot1';

clear
x0=[30 0 -37];
tspan=[0,30];
[t,x]=ode45('fun2_4',tspan,x0);
plot(t,x(:,1),'k',t,x(:,2),'k--',t,x(:,3),'k:');
xlabel('time(tian)t0=0,tf=30');
legend('x1','x2','x3');
grid;

ChaChing · 发表于 2009-4-13 22:35

建议楼主看下本版规则! 给齐完整格式：出错代码和出错提示!

xxzz · 发表于 2009-4-14 15:19

Error in ==> ode45 at 467
idx = oldnout+1:nout;
Error in ==> pp at 4
[t,x]=ode45('fun2_4',tspan,x0);

[ 本帖最后由 ChaChing 于 2009-4-14 23:23 编辑 ]

ChaChing · 发表于 2009-4-14 23:21

LZ的程序试跑过, 好像当机似的, 偶是使用Ctr-C强迫切断!
猜测会不会LZ这个不收敛?
待高手确认

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

[编程技巧] 大侠们帮看一下，为什么运行不了

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

各位大侠帮看一下，混沌系统的程序，怎么运行不了