(完整版)基于matlab的Lorenz系统仿真研究.docVIP

基于 Matlab 的 Lorenz 系统仿真研究 摘要：本文利用 matlab 这一数学工具对 Lorenz 系统进行了研究。首先使用 matlab 分析求解 Lorenz 方程，利用 matlab 的绘图功能，直观地观察了 Lorenz 混沌吸引子的三维图形，并简单观察了 Lorenz 混沌系统对初值的敏感性； 然后对 Lorenz 系统进行仿真，比较分析在不同参数下的 Lorenz 系统仿真 结果；最后验证了通过添加反馈控制的方式， 可以使 Lorenz 方程不稳定的 平衡点成为稳定的平衡点。 关键词： Lorenz 系统；matlab ；混沌系统 1. 引言 Lorenz方程是由美国著名的气象学家 Lorenz在1963年为研究气候变化，通过 对对流实验的研究，建立的三个确定性一阶非线性微分方程。这三个方程是混沌 领域的经典方程， Lorenz系统也是第一个表现奇怪吸引子的连续动力系统，具有 着举足轻重的作用。 Lorenz方程的表达式如下： { dx dt = σ(y- x) dy dt = ( μ- z) x - y dz = -bz + xy dt 其中，σ、μ、 b为正实常数。 本文利用matlab 这一数学工具，对Lorenz系统进行了研究，得到了仿真结果， 加深了对 Lorenz系统的认识。 2.matlab 求解Lorenz方程并绘图 首先建立 m文件“Lorenz.m”来定义 Lorenz 方程，固定σ=10，μ=30，b=8/3， 程序如下所示： function dx=Lorenz(t,x) dx=[-10*(x(1)-x(2));30*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-2.6667*x(3)]; end 然后利用 ode45 命令来求解 Lorenz 方程并绘制图形，初值取 x=y=z=0.1 。程 序如下所示： clf x0=[0.1,0.1,0.1]; [t,x]=ode45(Lorenz,[0,100],x0); subplot(2,2,1) plot(x(:,1),x(:,3)) title((a)) subplot(2,2,2) plot(x(:,2),x(:,3)) title((b)) subplot(2,2,3) plot(x(:,1),x(:,2)) title((c)) subplot(2,2,4) plot3(x(:,1),x(:,2),x(:,3)) title((d)) 运行上述程序，可得到如下波形： (a) (b) 60 60 40 40 20 20 0 -20 -10 0 10 20 30 0 -40 -20 0 20 40 (d) (c) 40 20 50 40 30 20 0 10 -20 40 -40 -20 -10 0 10 20 30 20 0 -20 -40 -20 0 20 其中，图（a）为 Lorenz 混沌吸引子在 x-z 平面上的投影，图（ b）为Lorenz 混沌吸引子在 y-z 平面上的投影，图（ c）为 Lorenz 混沌吸引子在 x-y 平面上的 投影，图（ d）为 Lorenz 混沌吸引子的三维图。可以看到，混沌吸引子在各平面 上的投影类似于横写的“ 8”字形。 由于参数σ =10，μ=30，b=8/3 时为混沌系统，对初值具有敏感性，初值很 小的差异会引起系统行为的显著改变。因此，将初值改为 x=z=0.1,y=0.11, 绘制 此时混沌吸引子在 x-z 平面上的投影，并与初值为 x=y=z=0.1 时混沌吸引子在 x-z 平面上的投影放在同一张图中比较。为了区别两者，初值为 x=y=z=0.1 时混沌吸 引子在 x-z 平面上的投影用蓝色， 初值改为 x=z=0.1,y=0.11 时混沌吸引子在 x-z 平面上的投影用红色。程序如下所示： clf x0=[0.1,0.1,0.1]; [t,x]=ode45(Lorenz,[0,100],x0); plot(x(:,1),x(:,3)) hold on x0=[0.1,0.1,0.1]; [t,x]=ode45(Lorenz,[0,100],x0); x0=[0.1,0.11,0.1]; [t,x]=ode45(Lorenz,[0,100],x0); plot(x(:,1),x(:,3),r*) 得到的图形如下所示： 60 50 40 30 20 10 0 -20 -15 -10 -5 0 5 10 15 20 25 可以看到， 虽然初值只有 0.01 的改变，红色与蓝色图形明显不重合， 这证明 了系统的敏感性。 3.matlab 对 Lorenz 系统仿真 首先利用 matlab 的 Simulink