数学建模实验三Lorenz模型与食饵模型

1、数学建模实验三 Lorenz模型与食饵模型一、实验目的1、学习用Mathematica求常微分方程的解析解和数值解，并进行定性分析；2、学习用MATLAB求常微分方程的解析解和数值解，并进行定性分析。二、实验材料2.1问题图3.3.1是著名的洛仑兹(E.N.Lorenz)混沌吸引子，洛仑兹吸引子已成为混沌理论的徽标，好比行星轨道图代表着哥白尼、开普勒理论一样。洛仑兹是学数学出身的，1948年起在美国麻省理工学院（MIT）作动力气象学博士后工作，1963年他在大气科学杂志上发表的论文确定性非周期流是混沌研究史上光辉的著作。以前科学家们不自觉地认为微分方程的解只有那么几类：1)发散轨道；2)不动点

2、；3)极限环 ；4)极限环面。除此以外，大概没有新的运动类型了，这是人们的一种主观猜测，谁也没有给出证明。事实上这种想法是非常错误的。1963年美国麻省理工学院气象科学家洛仑兹给出一个具体模型，就是著名的Lorenz模型，清楚地展示了一种新型运动体制：混沌运动，轨道既不收敛到极限环上也不跑掉。而今Lorenz 模型在科学与工程计算中经常运用的问题。例如，数据加密中。我们能否绘制出洛仑兹吸引子呢？图3.3.1 洛仑兹(E.N.Lorenz)混沌吸引子假设狐狸和兔子共同生活在同一个有限区域内，有足够多的食物供兔子享用，而狐狸仅以兔子为食物.x为兔子数量,y表狐狸数量。假定在没有狐狸的情况下，兔子增

3、长率为400。如果没有兔子，狐狸将被饿死，死亡率为90。狐狸与兔子相互作用的关系是，狐狸的存在使兔子受到威胁，且狐狸越多兔子增长受到阻碍越大，设增长的减小与狐狸总数成正比，比例系数为0.02。而兔子的存在又为狐狸提供食物，设狐狸在单位时间的死亡率的减少与兔子的数量成正比，设比例系数为0.001。建立数学模型，并说明这个简单的生态系统是如何变化的。2.2预备知识1、求解常微分方程的Euler折线法求初值问题 （12.1）在区间上的数值解，并在区间插入了结点。由导数的定义，即微商。（右端称为差商）从而可在每个结点上用差商来近似替代导数，将微分方程转化为代数方程组（此处的代数方程组常称为差分方程）

4、,加上初值条件则可确定一组解。求解这一差分方程即可得到微分方程初值问题的数值解。变形上述方程有 ,记,，从而，则有 这就是求解微分方程初值问题的欧拉(Euler)折线法。之所以称为欧拉折线法是因为：就几何角度而言，所求得的近似解是初值问题精确解的折线逼近，而且此折线的起点是初值条件所对应的点。2、微分方程的Mathematica求解（1）求解命令 有两个命令：DSolve 与NDSolve。命令格式分别为 DSolve方程，y，x NDSolve方程，y，x，xl，x2。 其中方程必须为微分方程及相应初始条件，x，xl，x2说明要给出数值解的范围为区间x1，x2。（2）使用的注意事项 方程中的

5、函数应写成完整形式yx,以表明y是x的函数； 方程应写成=的形式； 重复使用时，应随时清除要涉及变量的以前定义，方法是Cleary； 使用NDSolve时，所加初始条件的个数应等于微分方程的阶数，同时方程中也不含其它参数，否则给不出正确结果。（3）解的表示形式Mathematica给出的微分方程的解是以纯函数（或数学中的算子）定义的形式给出的，例如：DSolvey'x+ 3*yx=2x,y,x的结果是3、微分方程的MATLAB求解（1）求解析解命令dsolve；（2）求数值解命令ODE或 Simulink。2.3建立模型问题（1）的洛仑兹吸引子可以用下面的微分方程得到，著名的Loren

6、z 模型的状态方程可表示为若令 且初值为，e 为一个小常数，假设。求微分方程的数值解，并绘制出时间曲线与相空间曲线。问题（2）是著名的食饵模型，数学模型为 2.4练习题1、求解微分方程的通解。求解的Mathematica命令为： DSolvey'x+2*x*yx= x*E(-x2),y,x 或者 DSolveDyx,x+2*x*yx= x*E(-x2),y,x2、求微分方程在初始条件下的特解。应给出的命令为：DSolvex*y'x+ yx-Ex=0,y1=2E,y,x3、求在初始条件下的特解，并画出解的图形。要求分别求解析解与数值解并作比较。清除要涉及变量的命令为：Clearx

7、,y求解析解的命令为：sc=DSolve(x2-1)y'x+2x*yx-Cosx=0,y0=1,y,x 画解析解图像的命令为：y=y/.sc1g1=Plotyx,x,0,1,PlotStyle->RGBColor1,0,0注：也可将画图范围变为Plotyx,x,0,4求数值解的命令为：sn=NDSolve(x2-1)y'x+2x*yx-Cosx=0,y0=1, y,x,0,1 画数值解图像的命令为： y=y/.sn1g2=Plotyx,x,0,1 比较解析解图像与数值解图像的命令为： Showg1,g24、求微分方程组 在初始条件,下的解，并画出解函数的图形。求解微分方程

8、组的命令为：Clearx,y,txy=DSolvex't+5*xt+yt=Et,y't-xt-3*yt=0,x0=1,y0=0,x,y,t画解的相位图的命令为：y=y/.xy1;x=x/.xy1;ParametricPlotxt,yt,t,0,3,PlotRange->-10,2,0,5注：图中反应出y随x的变化关系。三、实验准备 认真阅读实验目的与实验材料后要正确地解读实验，在此基础上制定实验计划（修改、补充或编写程序，提出实验思路，明确实验步骤），为上机实验做好准备。四、实验思路提示4.1实验步骤1、求解问题（2）中的食饵模型的微分方程组，并画出解的图形和相位图。（1

9、）以x=800，y=100为初始值，计算x（t）,y（t），当t0，14时的数据。绘出解的图形，并分析捕食者和被捕食者的数量变化规律。可以先用下面的命令求解析解：Clearx,y,txy=DSolvex't=4*xt-0.02*xt*yt, y't=-0.9*yt+0.001*xt*yt,x0=800, y0=100,x,y,t注：可以发现不能求出解析解。修改代码如下，可以求数值解：Clearx,y,txy=NDSolvex't=4*xt-0.02*xt*yt, y't=-0.9*yt+0.001*xt*yt,x0=800, y0=100,x,y,t,0,14绘

10、出解的图形：y=y/.xy1;x=x/.xy1;Plotxt,yt,t,0,14,PlotStyle->RGBColor0,0,1,RGBColor1,0,0 图3.3.2 捕食者和被捕食者的数量变化（2）以x为横坐标，y为纵坐标绘制相位图。根据图形分析被捕食者数量增加(减少)对捕食者数量的影响。绘制相位图的命令：ParametricPlotxt,yt,t,0,14 图3.3.3 相位图2、用MATLAB求解问题（1）中Lorenz 模型的微分方程。（1）打开MATLAB的编辑器；（2）在编辑器中用下面的几个语句描述微分方程，并将其保存在lorenzeq.m的m文件中：f unction

11、 xdot = lorenzeq(t,x)xdot=-8/3*x(1)+x(2)*x(3);      -10*x(2)+10*x(3);      -x(1)*x(2)+28*x(2)-x(3);（3）新建命令文件：t_final=100; x0=0;0;1e-10;   t,x=ode45('lorenzeq',0,t_final,x0);   plot(t,x),   figure; plot3(x(:,1),x(:,2),x(:,3); axis(10 40 -20 20 -20 20);绘制出时间曲线与相空间曲线，如下图所示。图3.3.4时间曲线与相空间曲线4.2思考问题1、运用Mathematica求解Lorenz 模型的微分方程组，从而了解系统状态是如何变化的。2、求解以下问题（广告的效用）： 某公司生产一种耐用消费品，产品一上市，该公司