weifen3.m

rossler混沌系统图像

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[t,Y1] = ode45(@rossler2,[0 150],[0 1 1],options);%rigid or rossler2 or duffing2  2.4,4,4.23,4.3,5,5.7
[t,Y2] = ode45(@rossler4,[0 150],[0 1 1],options);
[t,Y3] = ode45(@rossler423,[0 150],[0 1 1],options);
[t,Y4] = ode45(@rossler43,[0 150],[0 1 1],options);
[t,Y5] = ode45(@rossler5,[0 150],[0 1 1],options);
[t,Y6] = ode45(@rossler57,[0 150],[0 1 1],options);
subplot(3,2,1),plot(Y1(:,1),Y1(:,2)), xlabel('x'),ylabel('y'),zlabel('z'),title('rossler  tspan=150  u=2.4');
subplot(3,2,2),plot(Y2(:,1),Y2(:,2));  xlabel('x'),ylabel('y'),title('rossler   u=4'); 
subplot(3,2,3),plot(Y3(:,1),Y3(:,2)); xlabel('x'),ylabel('y'),title('rossler   u=4.23');
subplot(3,2,4),plot(Y4(:,1),Y4(:,2)); xlabel('x'),ylabel('y'),title('rossler  u=4.3 ');
subplot(3,2,5),plot(Y4(:,1),Y4(:,2)); xlabel('x'),ylabel('y'),title('rossler   u=5');
subplot(3,2,6),plot(Y4(:,1),Y4(:,2)); xlabel('x'),ylabel('y'),title('rossler   u=5.7');