一.利用MATLAB数学软件对上面微分方程求解,进行数值模.PDFVIP

Lorenz 混沌系统 美国气象学家洛伦兹于1963 年在大气科学杂志上提出第一个表现奇异吸引子的动力学系统。该混沌系统模型 可以用下列微分方程组描述： dx 10(x y ) dt dy 28x y xz dt dz 8 xy  z dt 3 一． 利用MATLAB 数学软件对上面微分方程求解，进行数值模拟。 编程表示Lorenz 方程，其中y （1）表示x，y （2）表示y，y （3）表示z ，程序如下： function dy=Lorenz(t,y) dy=zeros(3,1); dy(1)=10*(-y(1)+y(2)); dy(2)=28*y(1)-y(2)-y(1)*y(3); dy(3)=y(1)*y(2)-8*y(3)/3; end 用ode45 命令求解，程序如下 [t,y]=ode45(Lorenz,[0,30],[12,2,9]); figure(1) plot(t,y(:,1)) figure(2) plot(t,y(:,2)) figure(3) plot(t,y(:,3)) figure(4) plot3(y(:,1),y(:,2),y(:,3)); view([20,42]) 取时间为0 到30 秒，设初始值为（12,2,9） 求解结果,可以得到上式的结果图像： y(1)与时间的关系 y(2)与时间的关系 y(3)与时间的关系 Y 的三个分量与时间的关系 动态绘制吸引子，程序如下： clc; clear; fig=figure; aviobj=avifile(C:\Users\Administrator\Dynamical Systems\Lorenz); for i = 1:200 [t,y]=ode45(Lorenz,[0,0.1*i],[12,2,9]); clf axis([-20,20,-25,25,10,50]) view([20,42]); hold on plot3(y(:,1),y(:,2),y(:,3)); MOV=getframe(fig); aviobj=addframe(aviobj,MOV); end close(fig) aviobj=close(aviobj) 具体动态过程见视频。 二．Lorenz 吸引子 可以通过图像得到Lorenz 吸引子的图像，如下： 吸引子图像 x-y 平面吸引子图像 x-z 平面吸引子图像 y-z 平面吸引子图像 三．倍周期 通过系数的调试可以得到Lorenz 混沌的单倍周期，2 倍周期和多倍周期，如下： 单倍周期 2 倍周期 多倍周期 四．蝴蝶效应 Lorenz 方程的解对初始值有着高度敏感性，为验证这一