混沌系统分岔图

CHEN系统求解分岔图一：functiondy=chen(t,y)globalr;dy=zeros(3,1);dy(1)=35*(y(2)-y(1));dy(2)=-7*y(1)+28*y(2)-y(1)*y(3);dy(3)=y(1)*y(2)-r*y(3)end对r的分岔图：function[Xmax]=getmax(y)a=length(y);j=1;fori=(a-1)/2:ab=(y(i,1)-y(i-2,1))/2;c=(y(i,1)+y(i-2,1))/2-y(i-1,1);ify(i-2,1)=y(i-1,1)&y(i-1,1)=y(i,1)&c==0Xmax(j)=y(i-1,1);elseify(i-2,1)=y(i-1,1)&y(i-1,1)=y(i,1)Xmax(j)=y(i-1,1)-b^2/(4*c);j=j+1;endend求解：clearallglobalrt0=[0100];y0=[0.015,0.758,1.83];forr=1:5.5:1000[t,y]=ode45('FUN',t0,y0);[Xmax]=getmax(y(:,1));plot(r,Xmax,'b','markersize',1)holdonclearXmaxendxlabel('r')ylabel('Xmax')求解分岔图二：clc;clearglobalr；Z=[];forr=linspace(0,0.01,100)options=odeset(’RelTol’,1e-6,’AbsTol’,[1e-41e-41e-5]);[T,Y]=ode45(’FUN’,[0,50],[0.0150.7581.83],options);n=length(Y);fork=round(n/2):nifabs(Y(k,1))1Z=[Z,r+abs(Y(k,2))*i];endendendfigureplot(Z,’.’,’markersize’,1)title(’chen映射分岔图’)xlabel(’r’)，ylabel(’y’)