陈氏超混沌系统相图

陈氏超混沌系统相图（该系统由陈关荣提出，故命名：Chen’s hyperchaos）

建立函数：

function dydt=HyChen(t,y) 
global a b c d r 
a=35; b=3; c=12; d=7; r=0.58;%input('r='); 
dydt=[a*y(2,:)-a*y(1,:)+y(4,:) 
d*y(1,:)+c*y(2,:)-y(1,:)*y(3,:) 
y(1,:)*y(2,:)-b*y(3,:) 
y(2,:)*y(3,:)+r*y(4,:) 
];      

建立文件：

y0=[-4 -3 5 -6];
x0=0; %起始时间点 
xf=200;%终止时间点 
[t,ym]=ode45('HyChen',x0:0.001:xf,y0); 
figure(1) 
subplot(221) 
plot(ym(:,1),ym(:,2)); 
grid on
xlabel('y1');ylabel('y2'); 
subplot(222) 
plot(ym(:,1),ym(:,3)); 
grid on
xlabel('y1');ylabel('y3'); 
subplot(223) 
plot(ym(:,1),ym(:,4)); 
grid on
xlabel('y1');ylabel('y4'); 
subplot(224) 
plot(ym(:,2),ym(:,3)); 
grid on
xlabel('y2');ylabel('y3'); 
figure(2) 
subplot(221) 
plot(ym(:,2),ym(:,4)); 
grid on
xlabel('y2');ylabel('y4'); 
subplot(222) 
plot(ym(:,3),ym(:,4)); 
grid on
xlabel('y3');ylabel('y4'); 
subplot(223) 
plot3(ym(:,1),ym(:,2),ym(:,3)); 
grid on
xlabel('y1');ylabel('y2');zlabel('y3') 
subplot(224) 
plot3(ym(:,1),ym(:,2),ym(:,4)); 
grid on
xlabel('y1');ylabel('y2');zlabel('y4')      

显示图像：