biaozhun.m

此程序用于画混沌系统的相图及其庞加莱截面

function biaozhun
for w=1.5:0.1:2.5
%w=80;
[T,Y]=ode45(@chen,[0:0.01:100],[0.9 0.7 0.02 0.3 0.23],[],w);
%figure,plot(Y(4000:end,1),Y(4000:end,3));
figure,plot(Y(4000:end,1),Y(4000:end,3),'black');
%subplot(2,2,2),plot(Y(:,1),Y(:,2));
%subplot(2,2,3),plot(Y(:,2),Y(:,3));
end

function dy=chen(t,y,w)
dy=zeros(5,1);
%b=15;       
%c=35;    
%d=2.5;
%w=12;
dy(1)=10*(y(2)-y(1));
dy(2)=50*y(1)-2*y(1)*y(3)+y(4);
dy(3)=4*y(1)*y(1)-4*y(3);
dy(4)=sin(y(5))*y(2)*y(3);
dy(5)=w;