Lorenz混沌基于FPGA的Lorenz混沌系统verilog实现

Posted 2022-04-25 fpga&matlab

tags:

篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了Lorenz混沌基于FPGA的Lorenz混沌系统verilog实现相关的知识，希望对你有一定的参考价值。

1.软件版本

matlab2021a，quartusii

2.本算法理论知识

Lorenz在1963年洛伦兹在研究大气对流现象时，发现了第一个结构简单的三维自治混沌系统，这就是著名的“蝴蝶效应”模型，其数学模型为：

当r≤1时，对所有的x，y，z都成立，等式仅在r<1，x=y=z=0和r=1，x=y，z=0时才成立。

当r>1时，系统开始出现不稳定性，当r增加至rc时，有一个亚临界的Hopf分岔。

当r>rc时，在吸引子附近相邻的轨道平均来说是以指数分离的，所以两条开始在一起非常接近的轨道很快丧失所有的相关性，呈现混沌状。

3.核心代码


% Dx=(25a+10)(y-x)
% Dy=(28-35a)x-xz+(29a-1)y
% Dz=xy-(a+8)z/3
% 当0<=a<0.8，广义Lorenz系统
% 当a=0.8，广义吕系统
% 当0.8<a<=1，广义陈系统
clear;clc


 x=1.2;y=1.3;z=1.6;
%x=1 ;y=1 ;z=1 ;

dt=0.005;
a=10.0;c=28.0;b=2.666667;
figure;
for i=1:10000
    newx=x+a*(y-x)*dt;
    newy=y+(c*x-y-x*z)*dt;
    newz=z+(x*y-b*z)*dt;
    if i>1000
 %   plot(x,y)
    % plot(x,z)
    % plot(y,z)
     plot3(x,y,z)
    hold on
    end
    x=newx;y=newy;z=newz;
end

figure;
for i=1:10000
    newx=x+a*(y-x)*dt;
    newy=y+(c*x-y-x*z)*dt;
    newz=z+(x*y-b*z)*dt;
    if i>1000
    plot(x,y)
    % plot(x,z)
    % plot(y,z)
     %plot3(x,y,z)
    hold on
    end
    x=newx;y=newy;z=newz;
end


figure;
for i=1:10000
    newx=x+a*(y-x)*dt;
    newy=y+(c*x-y-x*z)*dt;
    newz=z+(x*y-b*z)*dt;
    if i>1000
   % plot(x,y)
     plot(x,z)
    % plot(y,z)
     %plot3(x,y,z)
    hold on
    end
    x=newx;y=newy;z=newz;
end


figure;
for i=1:10000
    newx=x+a*(y-x)*dt;
    newy=y+(c*x-y-x*z)*dt;
    newz=z+(x*y-b*z)*dt;
    if i>1000
    %plot(x,y)
    % plot(x,z)
     plot(y,z)
     %plot3(x,y,z)
    hold on
    end
    x=newx;y=newy;z=newz;
end