Rossler系统的分叉图Matlab实现

function u =rosser(t,x)

global  c

a=0.2; b=0.2; 振动资讯Mg8U%Q6K.h9T!K"u*i

x0=[0.5,0.4,0.3]';振动资讯%} pM'^7I:b

u =[-x(2)-x(3), x(1)+a*x(2), b+x(3)*(x(1)-c)]';振动资讯Q ?eV4y5G

 振动资讯9Q#piE*l

clear all

global c

zhang=[];振动资讯4x#nP#KX#\b5v

M=[2.5:0.001:6.5];振动资讯8Ic9S_b*tf_

counter=1;

for counter=1:length(M)

       c=M(counter);振动资讯#p jCp_!G9iZ3s

    i=2;振动资讯-G([9p%L(Mg#xm

xmax=0;振动资讯:n
|iHs&^*R
E

xmaxold=0;振动资讯3p+[u/[ J

frmdata=[];振动资讯 x"mu1tA^x0?W.Z

error=0;振动资讯jD oX6N

tspan=[0 300];振动资讯Dt_WQ6d3M

var=1;

y0=[0.5;0.3;0.2];振动资讯x&O W`#x;~

[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)振动资讯Tm6J/x#q/{C

            xmax=x(i,var);振动资讯?:dtwN
PbU

            if xmaxold ~= 0

                frmdata=[frmdata ; xmax xmaxold];

            end

                xmaxold=xmax;振动资讯h
j8W1Q,@8f

        end振动资讯m
F:u,J4H%{D*q:|]'C

        i=i+1;

end振动资讯*[[H'X:Y}C"Pr_

r= length(frmdata)-20:length(frmdata);振动资讯|W{j*fh!t:B?1n

fradata1(1,r)=frmdata(r,1);振动资讯n_)?-TSv&K3A8A,[!R

zhang=[zhang;fradata1];

end

plot(M,zhang,'k.','markersize',1);

        xlabel(sprintf('c'));

ylabel(sprintf('x'))振动资讯;^r?;cn@Y-j

 振动资讯X.r] lc5v)@)uj

Rossler Bifurcation versus the control parameter C