资源描述：
《混沌系统分岔图.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、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
2、-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)holdonclearXm
3、axendxlabel('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];endendendfigurep
4、lot(Z,’.’,’markersize’,1)title(’chen映射分岔图’)xlabel(’r’)，ylabel(’y’)
