陈氏超混沌系统相图

陈氏超混沌系统相图（该系统由陈关荣提出，故命名：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')

显示图像：