声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1470|回复: 3

[稳定性与分岔] 分叉图说明了什么运动

[复制链接]
发表于 2011-9-26 10:00 | 显示全部楼层 |阅读模式

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

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

x
分叉图
      

分叉图.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2011-9-26 10:01 | 显示全部楼层
希望高手多多帮助。
发表于 2011-9-27 11:50 | 显示全部楼层
把程序贴出来看看,我也遇到这样的问题,还没解决
不过这样应该是不对的 参数应该有问题
 楼主| 发表于 2011-9-27 19:42 | 显示全部楼层
方程如下:
function dy=accele(t,y,r1)
dy=zeros(2,1);
f=0.8; a=0.001; K=2000; L=0.1;
r0=3; w=3; r=r0+r1*sin(w*t);

k11=(r^2-1-(f^2+a*r)*(pi)^2)*(pi)^2;
k12=3*K^2*(pi)^4/8;
k13=3*L*(pi)^4/4;

dy(1)=y(2);
dy(2)=k11*y(1)-k12*(y(1))^3-k13*(y(1))^2*y(2);
分叉图程序:clear
A=0:0.01:3;
options=odeset('RelTol',1e-7);

hold on
for j=1:length(A)
    r1=A(j);
    tt=2*pi/3;
    [t,y]=ode45(@accele,[0:tt/100:80*tt],[0,0],options,r1);
    plot(A(j),y(8000:100:end,1),'k')
end
不知道哪里出了问题,还请多多帮助,谢谢。

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-19 08:16 , Processed in 0.070424 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表