声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 788|回复: 0

[综合讨论] 请问如下的分岔图错在那?

[复制链接]
发表于 2008-4-6 20:17 | 显示全部楼层 |阅读模式

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

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

x
在曲面y(2)=-0.5*y(1)+0.00004作poincare映射关于f与y(1)的分岔图,坐标f是0.00035到0.00055,y(1)是-0.03到 -0.01。程序如下:
global f;
rang=0.00035:0.00001:0.00055;
hold on
for f=rang
    [t,y]=ode45(@Xiang10,[0,10000],[-0.031,0.01,0.01,0.02]);
    if abs(y(2)+0.5*y(1)-0.0004<0.001)
        plot(f(j),y(:,1)),hold on
    end
end
其中
function dx=Xiang10(t,y);
global f;
w1=1.0;w2=0.515;
s1=0.03;
s2=0.004;
c1=0.005;
c2=0.005;f=0.000497;
dy=[-1/2*c1*y(1)-s2*y(2)+3*w2*w2/(4*w1)*y(3)*y(4);
-1/2*c1*y(2)+s2*y(1)-3*w2*w2/(8*w1)*(y(3)*y(3)-y(4)*y(4))+f/(2*w1);
-1/2*c2*y(3)-1/2*(s2-s1)*y(4)-1/4*(2*w1-w2)*(y(1)*y(4)-y(2)*y(3));
-1/2*c2*y(4)+1/2*(s2-s1)*y(3)-1/4*(2*w1-w2)*(y(1)*y(3)+y(2)*y(4))];

[ 本帖最后由 eight 于 2008-4-7 16:09 编辑 ]
回复
分享到:

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-11-18 08:16 , Processed in 0.092721 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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