binran115 发表于 2011-7-8 20:04

大家看下下面的分岔图

本帖最后由 binran115 于 2011-7-8 20:21 编辑

下面程序是我画的chen电路延迟系统的分岔图，结果这么也感觉不对呢，谁能指点一下呢，谢谢大家了！思想：chen系统是一个自治系统，所以画的方法我也是按自治系统的方法来画，采取找局部极大值方法选取poincare截面进行绘制。
chen电路延迟系统动态方程如下：
主程序如下（有人可能下不到，放上来大家给指导一下）：
clc
clear
global K
x0=rand(1,4);
range=;
for K=range
   sol = dde23(@chen_ddex1de,,,);
   %step=0.1;
   %x=0:step:200;
   x=linspace(0,100,11681);
   kk=deval(sol,x);
   data=kk(2,:);
   n=length(data)
   N=round(n/2)
   a=data(N-2);
   b=data(N-1);
for i=N:n
    if a <= b & b >= data(i)
      plot(K,b,'k.','markersize',1);
      hold on;
    else
    end
    a=b;
    b=data(i);
end
end




binran115 发表于 2011-7-9 18:03

没人知道这个么这弄吗，！

meiyongyuandeze 发表于 2011-7-10 20:22

参数取的太稀疏了吧！

lihaitao123 发表于 2011-7-11 10:05

回复 3 # meiyongyuandeze 的帖子

取细一点，间距小了。怎么看还是混沌那种感觉

mustangyhz 发表于 2011-7-27 09:44

回复 1 # binran115 的帖子

我用禹思敏的书上P176上的程序也是得到类似的图，你搞定了没有？
W1=[];
W2=[];
global a;
for a=8:0.002:9;
    =sim('chua',300);
    data=X(:,1);
    n=length(T);
    N=round(n/2);
    a1=data(N-2);
    b1=data(N-1);
    for i=N:n
      if(b1>=a1)&(b1>=data(i));
            plot(a,b1);
            hold on;
            axis on;
            W1=;
            W2=;
      else
      end
      a1=b1;
      b1=data(i);
    end
end
   
   

dwgboy 发表于 2011-9-10 10:57

路过不要错过，{:{27}:}
页: [1]
查看完整版本: 大家看下下面的分岔图