chendis.m

用于生成陈吸引子的程序

function chendis
clc,clf
[t,y]=ode45('chen_attractor',[0 30],[0,1,0]);
%表示在0－30秒内求解，在零时刻设y(1)为12，y(2)为2，y(3)为9
for j=1:3
      subplot(2,2,j);
      plot(t,y(:,j),'r')
      
end  

%plot(t,y(:,1))       %显示y(1)，即x与时间的关系图
%pause                %暂停
%plot(t,y(:,2))       %显示y(2)，即y与时间的关系图
%pause
%plot(t,y(:,3))       %显示y(3)，即z与时间的关系图
%pause
%--------------------
subplot(2,2,4);
plot3(y(:,1),y(:,2),y(:,3));
view([10,30]);
grid on
subplot(2,2,1);
title('y(1)，即x与时间的关系图')
subplot(2,2,2);
title('y(2)，即y与时间的关系图')
subplot(2,2,3);
title('y(3)，即z与时间的关系图')
subplot(2,2,4);
title('三维关系图')