混沌表达式：
x(n+1)=ax(n)+by(n)z(n)
y(n+1)=cy(n)+dx(n)z(n)
z(n+1)=ex(n)y(n)+kz(n)+mx(n)w(n)
w(n)=ny(n)

生成混沌：

function s=odone(t,y)
a=8;
b=-1;
c=-40;
d=1;
e=2;
m=1;
n=-2;
k=-14;
s=[a*y(1)+b*y(2)*y(3);
      c*y(2)+d*y(1)*y(3);
      e*y(1)*y(2)+k*y(3)+m*y(1)*y(4);
      n*y(2)];
end

绘制吸引子曲线：

clc;clear;
[t,b]=ode45(@odone,[-100,100],[0.35;0.21;0.43;0.89]);
figure
plot(b(:,1),b(:,2));
title('x-y plane');
xlabel('x');
ylabel('y');
figure
plot(b(:,1),b(:,3));
title('x-z plane');
xlabel('x');
ylabel('z');
figure
plot(b(:,1),b(:,4));
title('x-w plane');
xlabel('x');
ylabel('w');
figure
plot3(b(:,1),b(:,2),b(:,3));
title('x-y-z plane');
xlabel('x');
ylabel('y');
zlabel('z');

结果：

以上内容仅作为个人学习笔记，无商业行为。