经典Lorenz混沌映射代码实现（MATLAB）

最新推荐文章于 2024-09-09 11:41:27 发布

颦颦的账号

最新推荐文章于 2024-09-09 11:41:27 发布

阅读量3.6k

点赞数 6

分类专栏：信息安全文章标签： matlab 开发语言图像处理

本文链接：https://blog.csdn.net/cj1628342983/article/details/122043730

版权

本文展示了如何在MATLAB中实现经典的Lorenz混沌映射，通过代码生成了x-y、x-z、y-z平面图以及三维的x-y-z图，详细解析了映射过程。

摘要由CSDN通过智能技术生成

创建一个.m文件，复制以下内容

h=0.002;n=50000;t=800;a=10;b=28;c=8/3;
% 生成1×n的全零矩阵
xn=zeros(1,n);yn=zeros(1,n);zn=zeros(1,n);
%设定初始值
x0=1.1;y0=2.2;z0=3.3;
% 龙格库塔法进行离散化,获得3个长50000的序列
for i=1:n+t
k11=-a*(x0-y0);k12=-a*((x0+k11*h/2)-y0);
k13=-a*((x0+k12*h/2)-y0);k14=-a*((x0+k13*h)-y0);
x1=x0+(k11+k12+k13+k14)*h/6;
 
k21=b*x1-y0-x1*z0;k22=b*x1-(y0+k21*h/2)-x1*z0;
k23=b*x1-(y0+k22*h/2)-x1*z0;k24=b*x1-(y0+k23*h/2)-x1*z0;
y1=y0+(k21+k22+k23+k24)*h/6;
 
k31=x1*y1-c*z0;k32=x1*y1-c*(z0+k31*h/2);
k33=x1*y1-c*(z0+k32*h/2);k34=x1*y1-c*(z0+k33*h/2);
z1=z0+(k31+k32+k33+k34)*h/6;
   
x0=x1;y0=y1;z0=z1;
if i>t
    xn(i-t)=x1; yn(i-t)=y1; zn(i-t)=z1;
end 
end
%x-y平面
figure(1);plot(xn,yn);
%x-z平面
figure(2);plot(xn,zn);
%y-z平面
figure(3);plot(yn,