声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1709|回复: 5

[混合编程] 这样的分岔图怎么解释

[复制链接]
发表于 2015-11-16 11:16 | 显示全部楼层 |阅读模式

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

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

x
请问大神,图片中的分岔图怎么解释啊,说明了什么
QQ图片20151116111740.jpg
回复
分享到:

使用道具 举报

发表于 2015-11-17 08:48 | 显示全部楼层
从图上看感觉你的计算没有收敛或者没有去除必要的瞬态相
猜测该系统的非线性演化过程为
周期运动->hopf分岔->周期运动->hopf分岔
或者是
周期运动->倍周期分岔->周期运动->hopf分岔

评分

1

查看全部评分

 楼主| 发表于 2015-11-18 09:35 | 显示全部楼层
谢谢,这样之后是正确的吗?
QQ图片20151118093630.png

点评

还是差点意思  详情 回复 发表于 2015-11-18 09:46
发表于 2015-11-18 09:46 | 显示全部楼层
qingting11 发表于 2015-11-18 09:35
谢谢,这样之后是正确的吗?

还是差点意思,可以考虑再多迭代些周期,前面去掉的计算结果再多一些
 楼主| 发表于 2015-11-18 16:23 | 显示全部楼层
能具体点吗,现贴上程序,指点一下:
clear;
clc;
syms  w0 T

a=2000:5:10000;
w0=pi/6;
T=2*pi/w0;

for j=1:length(a)
    [t,y]=ode45('kutalunwensancifang',[0:T/200:50*T],[0 0.1],[],a(j),w0);
    plot(a(j),y(3000:200:end,1),'b','marker','.','markersize',5)
end
xlabel('动载荷p1/N');
ylabel('位移x/m');
发表于 2015-11-18 17:32 | 显示全部楼层
for j=1:length(a)
    [t,y]=ode45('kutalunwensancifang',[0:T/200:200*T],[0 0.1],[],a(j),w0);
    plot(a(j),y(20000:200:end,1),'b','marker','.','markersize',5)

end

评分

1

查看全部评分

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

本版积分规则

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

GMT+8, 2024-11-5 15:51 , Processed in 0.081715 second(s), 27 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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