rosslerdis.asv

用来计算rossler吸引子的程序

function Rosslerdis
clear
clc,clf,
[t,y]=ode45('Rossler',[0 :0.05:2000],[-1,0,1]);
%表示在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  

X = y(:,1);
Y = y(:,2);
Z = y(:,3);

subplot(2,2,4);
%plot(y(:,2),y(:,3))
plot3(X,Y,Z);
view([20,42]);
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('三维关系图')
%save 'd:\rossler.dat' y -ascii;