哪位大神帮我看看我的分岔图的程序哪出错了啊

我恨李雅普诺夫

分岔图程序如下：

1. clear;clc;clf
   
2. global f
   
3. range=linspace(5,20,100);
   
4. for f=range
   
5. y0=[0.001 0.01 0.001 0.010]';
   
6. period=(2*pi)/f;
   
7. step=period/100;  
   
8.    tspan=[0:step:100*period];
   
9.   [t,x]=ode45('Driver_vehicle1',tspan,y0);
   
10.   plot(f,x(5000:100:end,1),'k.','markersize',10);
    
11. hold on
    
12. end

复制代码

跑出来的图根本就不是分岔图啊，我测试了几个f=5，8，10三种情况，庞加莱映射图都是几个点，但我跑的分岔图每个f对应的只是一个单一点，跪求大神解决啊。。。

我恨李雅普诺夫

求助啊。。。我用这段程序跑其他典型混沌能跑出来分岔图，但是我的自己的微分方程就是出不来啊。。。。。

gghhjj

没头没尾的一个问题，没办法回答

系统没给，庞加莱映射图如何获取没给，庞加莱映射图结果什么样没有

kuri2012

根据，以前的帖子，你好像，没取到周期点

gghhjj

kuri2012 发表于 2014-5-13 20:09
根据，以前的帖子，你好像，没取到周期点

1. plot(f,x(5000:100:end,1),'k.','markersize',10);

复制代码

是取得周期点

我恨李雅普诺夫

gghhjj 发表于 2014-5-17 06:37
是取得周期点

我取的是周期点啊，100个点是一个周期。。

s13081104004

我恨李雅普诺夫 发表于 2014-6-30 20:34
我取的是周期点啊，100个点是一个周期。。

那有可能你的系统没有出现分叉，只有稳定的单周期

mumianhua

试着把初始条件改变一下  试试看   

huizi_nice

我恨李雅普诺夫 发表于 2014-2-21 22:48
求助啊。。。我用这段程序跑其他典型混沌能跑出来分岔图，但是我的自己的微分方程就是出不来啊。。。。。

如果用你自己的方程能画出来Poincare映射图可出不来分岔图那就改变初值和参数试试

cuocheng

{:{44}:}