bif_rosser.m

rossler分岔代码

clear all
global c
zhang=[];
M=[2.5:0.01:6.5];
counter=1;
for counter=1:length(M)
    counter
    c=M(counter);
    i=2;
    xmax=0;
    xmaxold=0;
    frmdata=[];
    error=0;%此句的作用?
    tspan=[0 200];
    var=1;
    y0=[0.5;0.3;0.2];
    [t,x]=ode45('rosser',tspan,y0);
    while i<size(x,1)
        if x(i-1,var)<x(i,var)&x(i+1,var)<=x(i,var)
            xmax=x(i,var);
            if xmaxold~=0
                frmdata=[frmdata;xmax xmaxold];
            end
              xmaxold=xmax;
            end
            i=i+1;
        end
        r=length(frmdata)-20:length(frmdata);
%           r=1:length(frmdata);
        fradata1(1,r)=frmdata(r,1);
        zhang=[zhang;fradata1];
    end
    plot(M,zhang,'k.','markersize',1);
%      plot(M,fradata1,'k.','markersize',1);
    xlabel(sprintf('c'));
    ylabel(sprintf('x'))