antonio1111 发表于 2012-5-13 10:55

matlab求解混沌动力学方程

为什么我用matlab求解混沌动力学方程组的时候，matlab一直响应busy,我仔细看过了方程，没有敲错，求各位大神相助！代码如下：
function dy=der(t,y)
dy=zeros(2,1);
dy(1)=4.938*10^2*y(2)*(y(1)+1)-1.2*y(1)*(y(1)+1);
dy(2)=1.458*0.4/3*(1+y(1)*(t-2.7))-1.458*y(2)-0.486*y(1)-y(2)*(y(1)+1)+2.43*10^(-3)*y(1)*(y(1)+1);
>>=ode45(@der,[-100 100],[-0.1;-2]);
plot(y(:,1),y(:,2))

ChaChing 发表于 2012-5-14 01:30

matlab一直响应busy
为何不先将范围取小些!?

antonio1111 发表于 2012-5-14 09:55

回复 2 # ChaChing 的帖子

取了[-0.5,0.5]后可以运行了，但是图像不太对，取小了后不是混沌状态了，纠结中。。。

antonio1111 发表于 2012-5-15 18:34

各位有做过混沌的吗？出来的图像不对啊，帮下忙吧。

antonio1111 发表于 2012-5-15 20:46

方程中有一个延迟时间t,可以直接写吗？如同上面所写。然后再用ode45求解，可是出来的图像一直往里收敛，不是闭合的，怎么办啊？
页: [1]
查看完整版本: matlab求解混沌动力学方程