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

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

function dy=jxw(t,y);
global det
global w
global r
dy=;

functionpoincare
global det
global r
global w
det=0.5;w=1;
r=7.5;
y0=;
tspan=;
=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

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

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

回复 地板 jiangxiaowei123 的帖子

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

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

回复 5楼 无水1324 的帖子

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

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

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

回复 6楼 wxl789 的帖子

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

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

回复 7楼 无水1324 的帖子

呵呵
觉得这个确实是看运气了:@P
页: [1]
查看完整版本: 各位老兄帮忙看看这个Poincare怎么画成这个样子