MATLAB--第七章.ppt

附：蝴蝶效应 蝴蝶效应这一说法来源，是气象学家洛仑兹制作了一个电脑程序，模拟气候的变化，并用图像来表示。最后他发现，图像是很敏感地依赖于初值，而且十分像一只蝴蝶张开的双翅，因而他形象的将这一图形以“蝴蝶扇动翅膀”的方式进行阐释。蝴蝶效应是指对初始条件敏感性的一种依赖现象。输入端微小的差别会迅速放大到输出端。蝴蝶效应在经济生活中比比皆是：中国宣布发射导弹，港台100亿美元流向美国。今天的蝴蝶效应”或者“广义的蝴蝶效应”已不限于当初洛仑兹的蝴蝶效应仅对天气预报而言，而是一切复杂系统对初值极为敏感性的代名词或同义语。 核心理念：看似微不足道的细小变化，却能以某种方式对社会产生微妙的影响，甚至影响整个社会系统的正常运行。关注细节，防微杜渐，注重关联，控制全局。 * */19 * 微分方程与计算机模拟 常微分方程数值求解方法 追击问题的计算机仿真 有阻力抛射曲线实验 思考题与练习题 ? ? ? ? 例7.1 以 1994 年我国人口为12亿为初值，求解马尔萨斯模型常微分方程 N(t)表示人口数量,取人口变化率r =0.015,微分方程 function z=fun1(t,N) z=0.015*N; ode23(fun1,[1994,2020],12) [T,N]=ode23(fun1,[1994,2020],12) 命令窗口 ? 编辑窗口 ? 例2. Logistic模型 创建微分方程右端函数: function z=fun2(t,u) z=u.*(1-u); ode23(fun2,[0,6],1.8) ode23(fun2,[0,6],0.2) 在命令窗口求数值解 [x,y]=meshgrid(0:.25:6,0:.05:2); k=y.*(1-y); d=sqrt(1+k.^2); px=1./d;py=k./d; quiver(x,y,px,py),hold on u=dsolve(Du=u*(1-u),u(0)=.2); v=dsolve(Dv=v*(1-v),v(0)=1.8); ezplot(u,[0,6]) ezplot(v,[0,6]) 根据微分方程右端函数 f(x，y)= y(1 – y)，区域 D = { (x, y) | 0≤ x ≤ 6, 0≤ y ≤ 2} 内未知函数的导数值，确定解函数曲线的切线对应单位向量，绘制向量场。 常微分方程组初值问题 一阶常微分方程组初值问题数值求解方法 [T,y] = ode23( F ,Tspan,y0) 其中, F是函数文件, 表示 微分方程右端函数 Tspan = [t0 Tfinal] —— 求解区域; y0 —— 初始条件 注: 函数F(t,y) 必须返回列向量. 数值解 y 的每一行对应于列向量T中的每一行数据 捕食者与被捕食者问题 海岛上有狐狸和野兔,当野兔数量增多时，狐狸捕食野兔导致狐群数量增长;大量兔子被捕食使狐群进入饥饿状态其数量下降;狐群数量下降导致兔子被捕食机会减少,兔群数量回升。微分方程模型如下 计算 x(t)，y(t) 当t∈[0，20]时的数据。绘图并分析捕食者和被捕食者的数量变化规律。 x(0)= 100 y(0)=20 创建MATLAB的函数文件 function z=fox(t,y) z(1,:)=y(1)-0.015*y(1).*y(2); z(2,:)=-y(2)+0.01*y(1).*y(2); Y0=[100,20]; [t,Y]=ode23(fox,[0,20],Y0); x=Y(:,1);y=Y(:,2); figure(1),plot(t,x,b,t,y,r) figure(2),plot(x,y) 求微分方程数值解并绘解函数图形 ------兔子数量; ------狐狸数量 兔-狐数量 变化相位图 “蝴蝶效应”来源于洛伦兹一次讲演。模型如下 求微分方程数值解, 并绘出解函数在Y-X平面的投影曲线 取 ?=8/3，?=10，?=28。 x(0)=0，y(0)=0，z(0)=0.01。 t∈[0，80]， 记向量 [y1，y2，y3] = [x，y，z]，创建MATLAB函数文件如下 function z=flo(t,y) z(1,:)=-8*y(1)/3+y(2).*y(3); z(2,:)=-10*(y(2)-y(3)); z(3,:)=-y(1).*y(2)+28*y(2)-y(3); 用MATLAB命令求解并绘出Y-X平面的投影图 y0=[0;0;0.01]; [x,y]=ode23(flo,[0, 80],y0); figure(1),plot(y(:,2),y(:,1)) figure(2),comet3(y(:,1),y(: