大家看下下面的分岔图

binran115

下面程序是我画的chen电路延迟系统的分岔图，结果这么也感觉不对呢，谁能指点一下呢，谢谢大家了！思想：chen系统是一个自治系统，所以画的方法我也是按自治系统的方法来画，采取找局部极大值方法选取poincare截面进行绘制。
chen电路延迟系统动态方程如下：
主程序如下（有人可能下不到，放上来大家给指导一下）：
clc
clear
global K
x0=rand(1,4);
range=[0:0.1:1.3];
for K=range
     sol = dde23(@chen_ddex1de,[0.3],[0.1,0.1,0.1],[0,500]);
     %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

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

meiyongyuandeze

参数取的太稀疏了吧！

lihaitao123

回复 3 # meiyongyuandeze 的帖子

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

mustangyhz

回复 1 # binran115 的帖子

我用禹思敏的书上P176上的程序也是得到类似的图，你搞定了没有？
W1=[];
W2=[];
global a;
for a=8:0.002:9;
    [T,X]=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=[W1;a];
            W2=[W2;b1];
        else
        end
        a1=b1;
        b1=data(i);
    end
end
   
   

dwgboy

路过不要错过，