混沌系统分岔图2页.docx

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 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求解