各位老兄帮忙看看这个Poincare怎么画成这个样子

jiangxiaowei123 · 发表于 2008-12-8 22:06

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

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

x

function dy=jxw(t,y);
global det
global w
global r
dy=[y(2);-det*y(2)-(y(1))^3+r*cos(w*t)];

function poincare
global det
global r
global w
det=0.5;w=1;
r=7.5;
y0=[0.001;0];
tspan=[0:0.1*2*pi:70000];
[t,y]=ode45('jxw',tspan,y0);
m=zeros(10000,1);
n=zeros(10000,1);
for i=1:10000
m(i,1)=y(10*i,1);
n(i,1)=y(10*i,2);
end
plot(m,n,'k.','markersize',1);

jiangxiaowei123 · 发表于 2008-12-8 22:08

:@( 本人愚钝，实在是没有什么办法了，来寻求高手帮忙
谢谢啦

baolun · 发表于 2008-12-9 23:58

你取的系数有问题吧，我画了下，最后稳态的是个1周期运动，poincare截面上是个点.

jiangxiaowei123 · 发表于 2008-12-10 13:23

是初始值问题，我换了个初始值[3，4]，做出来了，是这个结果。

无水1324 · 发表于 2008-12-11 20:16

恩
初值会导致这样的结果的，楼主也可以把你做的贴出来欣赏一下
另外，你这个初值是自己猜测的吗？

wxl789 · 发表于 2008-12-12 10:15

对。
楼主分享一下参数调试的经验。。。

我在求动力学系统混沌的时候调了整整5天。。。。。。快要抓狂了。:lol

无水1324 · 发表于 2008-12-12 20:41

呵呵，是啊，这种都没有什么规律可循的，都是碰运气

wxl789 · 发表于 2008-12-12 22:17

呵呵
觉得这个确实是看运气了:@P

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

[分形与混沌] 各位老兄帮忙看看这个Poincare怎么画成这个样子

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

请各位老兄帮我看看啊，

回复地板 jiangxiaowei123 的帖子

回复 5楼无水1324 的帖子

回复 6楼 wxl789 的帖子

回复 7楼无水1324 的帖子

[分形与混沌] 各位老兄帮忙看看这个Poincare怎么画成这个样子

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

请各位老兄帮我看看啊，

回复 地板 jiangxiaowei123 的帖子

回复 5楼 无水1324 的帖子

回复 6楼 wxl789 的帖子

回复 7楼 无水1324 的帖子

回复地板 jiangxiaowei123 的帖子

回复 5楼无水1324 的帖子

回复 7楼无水1324 的帖子