《混沌系统分岔图.docx》由会员分享,可在线阅读,更多相关《混沌系统分岔图.docx(2页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、精品文档,仅供学习与交流,如有侵权请联系网站删除CHEN系统求解分岔图一:function dy=chen(t,y)global r;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; for i=(a-1)/2:a b=(y(i,1)-y(i-2,1)/2; c=(y(i,1)+y(i-2,1)/2-y(i-1,1); if y(i-2,1)=y(i,1)&c=0
2、 Xmax(j)=y(i-1,1); elseif y(i-2,1)=y(i,1) Xmax(j)=y(i-1,1)-b2/(4*c); j=j+1; end end求解:clear all global r t0=0 100; y0=0.015,0.758,1.83; for r=1:5.5:1000 t,y=ode45(FUN,t0,y0); Xmax=getmax(y(:,1); plot(r,Xmax,b,markersize,1) hold on clear Xmax end xlabel(r) ylabel(Xmax)求解分岔图二:clc;clearglobal r;Z=;for r=linspace(0,0.01,100)options=odeset(RelTol,1e-6,AbsTol,1e-4 1e-4 1e-5);T,Y=ode45(FUN,0,50,0.015 0.758 1.83,options);n=length(Y);for k=round(n/2):nif abs(Y(k,1)1Z=Z,r+abs(Y(k,2)*i;endendendfigureplot(Z,.,markersize,1)title(chen映射分岔图)xlabel(r),ylabel(y)【精品文档】第 2 页