基于matlab的Lorenz系统的仿真可研报告.doc

MATLAB 课 程 期 末 作 业 报告完成的是大作业第七题 7. Simulink仿真在高等数学课程中的应用 宋沛儒 基于MATLAB/Simulink对Lorenz系统仿真研究 1.引言 1963年Lorenz通过观察大量大气现象并进行数值实验和理论思考，得到了一系列混沌运动的基本特征，提出了第一个奇异吸引子—Lorenz吸引子[1] ，Lorenz通过计算机模拟一个由三阶微分方程描述的天气模型时发现，在某些条件下同一个系统可以表现出非周期的无规则行为。Lorenz揭示了一系列混沌运动的基本特征，成为后人研究混沌理论的基石和起点，具有非常重要的意义。Lorenz系统方程如下： （1），，b，c为实数 本人利用了matlab，对Lorenz系统进行了研究，了对其的认知。 matlab求解orenz系统 首先创建文件“Lorenz.m”Lorenz方程，固定=10，b=2.6667，c=30,程序如下： 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 然后利用45（Runge-Kutta算法）求解并绘制图形，初值取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)) 运行后，得如下波形： ）Lorenz混沌吸引子在x-z上的投影，b）在y-z上的投影，c）在x-y上的投影，d）Lorenz混沌吸引子。都类似于“8”形。 . Lorenz系统对初值的敏感性固定参数=10，b=2.6667，c=30时，为混沌系统，初值具有敏感性初值很小的差异会引起系统大变化。例如上例中初值x=z=0.1y=0.11，绘制此时混沌吸引子在x-z上的投影，并与x=y=z=0.1同一张图比较。（为x=y=z=0.1投影用蓝色，初值为x=z=0.1y=0.11时投影用红色： clf x0=[0.1,0.1,0.1]; [t,x]=ode45(ex_lorenz,[0,100],x0); plot(x(:,1),x(:,3)) hold on x0=[0.1,0.11,0.1]; [t,x]=ode45(ex_lorenz,[0,100],x0); plot(x(:,1),x(:,3),r*) 得到图形如下： y仅变化图中红色与蓝色不重合出明显。证明Lorenz系敏感性 4.matlab对Lorenz系统的仿真 [1]可知在上述方程组（1）中，令，当c1时，系统有三个平衡点：， ，。当c=1时，系统在原点失去稳定。当c1时，原点是唯一的平衡点并且是汇点。 利用matlab的Simulink功能，搭建Lorenz系统模型，探讨参数对系统影响。仿真模型如图： 根据文献[1]的分析， 当参数0c1时，只有一个稳定平衡点O（0,0,0）。取初值为x=y=z=2,参数c=0.5，仿真停止时间取为50，运行仿真。得到x、y、z的相图以及x-z，y-z，x-y的图形依次如下所示： 可见，系统很快地趋向并稳定在O（0,0,0），验证了前面所述。 当c1时，系统有三个平衡点：原点O(0,0,0)和S+，S-。c13.926时，不稳定流形最终螺旋地趋于与之同侧的平衡点S+或S-；当c=13.926时，不稳定流形刚好无限趋于原点O，即出现同宿轨；当c13.926时，不稳定流形将绕到另一侧，最终趋于与之异侧的S+或S-。可见，c是一个同宿分岔点。因此，取初值x=y=z=2，c=8，仿真停止时间为50，运行仿真，得到x、y、z的相图以及x-z，y-z，x-y的图形依次如下所示： 可以看到，系统趋于与之同侧的平衡点S+或S-。 取初值x=y=z=2，c=18，仿真停止时间为50，运行仿真，得到x、y、z的相图以及x-z，y-z，x-y的图形依次如下所示： 可以