微分方程模型0.ppt

微分方程微分方程 模模 型型 马马 戈戈导导 弹弹 追追 踪踪 问问 题题 设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰发射导弹，导弹头始终对准乙舰.如果乙舰以最大的速度v0(是常数)沿平行于y轴的直线行驶，导弹的速度是5v0，求导弹运行的曲线方程.又乙舰行驶多远时，导弹将它击中？解法一解法一（解析法）导导 弹弹 追追 踪踪 问问 题题由(1),(2)消去t整理得模型:导导 弹弹 追追 踪踪 问问 题题解法二解法二(数值解)令y1=y,y2=y1’，将方程（3）化为一阶微分方程组导导 弹弹 追追 踪踪 问问 题题1.建立m-文件eq1.m function dy=eq1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5*sqrt(1+y(1)^2)/(1-x); 2. 取x0=0，xf=0.9999，建立主程序ff6.m如下： x0=0，xf=0.9999 [x,y]=ode15s('eq1',[x0 xf],[0 0]); plot(x,y(:,1),’b.') hold on y=0:0.01:2; plot(1,y,’b*') 结论结论: 导弹大致在（导弹大致在（1，，0.2）处击中乙舰）处击中乙舰导导 弹弹 追追 踪踪 问问 题题解法三解法三(建立参数方程求数值解) 设时刻t乙舰的坐标为(X(t)，Y(t))，导弹的坐标为(x(t)，y(t)).导导 弹弹 追追 踪踪 问问 题题3．因乙舰以速度v0沿直线x=1运动，设v0=1，则w=5，X=1，Y=t导导 弹弹 追追 踪踪 问问 题题4. 解导弹运动轨迹的参数方程建立m-文件eq2.m如下： function dy=eq2(t,y) dy=zeros(2,1); dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2); dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2); 取t0=0，tf=2，建立主程序chase2.m如下： [t,y]=ode45('eq2',[0 2],[0 0]); Y=0:0.01:2; plot(1,Y,'-')， hold on plot(y(:,1),y(:,2),'*')导导 弹弹 追追 踪踪 问问 题题5. 结果见下图1导弹大致在（1，0.2）处击中乙舰，与前面的结论一致.图1图2 在chase2.m中，按二分法逐步修改tf，即分别取tf=1，0.5,0.25,…,直到tf=0.21时，得图2.结论：时刻结论：时刻t=0.21时，导弹在（时，导弹在（1，，0.21）处击中乙舰。

处击中乙舰导导 弹弹 追追 踪踪 问问 题题慢跑者与狗慢跑者与狗 一个慢跑者在平面上沿椭圆以恒定的速率v=1跑步,设椭圆方程为: x=10+20cost, y=20+5sint. 突然有一只狗攻击他. 这只狗从原点出发,以恒定速率w跑向慢跑者,狗的运动方向始终指向慢跑者.分别求出w=20,w=5时狗的运动轨迹.1. 模型建立 设时刻t慢跑者的坐标为(X(t)，Y(t))，狗的坐标为(x(t)，y(t)). 则X=10+20cost, Y=20+15sint, 狗从(0,0)出发,与导弹追踪问题类似，建立狗的运动轨迹的参数方程:慢跑者与狗慢跑者与狗2. 模型求解function dy=eq3(t,y) dy=zeros(2,1); dy(1)=20*(10+20*cos(t)-y(1))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); dy(2)=20*(20+15*sin(t)-y(2))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2);(1)w=20时时, 慢跑者与狗慢跑者与狗建立m-文件eq3.m如下:取t0=0，tf=10，建立主程序chase3.m如下： t0=0;tf=10; [t,y]=ode45('eq3',[t0 tf],[0 0]); T=0:0.1:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,'-') hold on plot(y(:,1),y(:,2),'*') 在chase3.m，不断修改tf的值,分别取tf=5, 2.5, 3.5,…,至3.15时,狗刚好追上慢跑者.慢跑者与狗慢跑者与狗 function dy=eq4(t,y) dy=zeros(2,1); dy(1)=5*(10+20*cos(t)-y(1))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); dy(2)=5*(20+15*sin(t)-y(2))/sqrt ((10+20*cos(t)- y(1))^2+(20+15*sin(t)-y(2))^2);(2) w=5时时建立m-文件eq4.m如下:慢跑者与狗慢跑者与狗取t0=0，tf=10，建立主程序chase4.m如下： t0=0;tf=10; [t,y]=ode45('eq4',[t0 tf],[0 0]); T=0:0.1:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,'-') hold on plot(y(:,1),y(:,2),'*') 在chase3.m，不断修改tf的值,分别取tf=20, 40, 80,…,可以看出,狗永远追不上慢跑者.慢跑者与狗慢跑者与狗地中海鲨鱼问题地中海鲨鱼问题 意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究，他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中，发现鲨鱼等的比例有明显增加（见下表），而供其捕食的食用鱼的百分比却明显下降.显然战争使捕鱼量下降，食用鱼增加，鲨鱼等也随之增加，但为何鲨鱼的比例大幅增加呢？ 地中海鲨鱼问题地中海鲨鱼问题 他无法解释这个现象，于是求助于著名的意大利数学家V.Volterra，希望建立一个食饵—捕食系统的数学模型，定量地回答这个问题.相关数据表：地中海鲨鱼问题地中海鲨鱼问题 该 模型反映了在没有人工捕获的自然环境中食饵与捕食者之间的制约关系，没有考虑食饵和捕食者自身的阻滞作用，是Volterra提出的最简单的模型.地中海鲨鱼问题地中海鲨鱼问题3．模型建立与求解．模型建立与求解 模型（一） 不考虑人工捕获  针对一组具体的数据用Matlab软件进行计算.地中海鲨鱼问题地中海鲨鱼问题首先，建立m-文件shier.m如下： function dx=shier(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-0.1*x(2)); dx(2)=x(2)*(-0.5+0.02*x(1));其次，建立主程序shark.m如下： [t,x]=ode45('shier',[0 15],[25 2]); plot(t,x(:,1),'-',t,x(:,2),'*') plot(x(:,1),x(:,2))地中海鲨鱼问题地中海鲨鱼问题求求 解解 结结 果果： 左图反映了x1（t）与x2（t）的关系。

可以猜测： x1（t）与x2（t）都是周期函数模型（二） 考虑人工捕获 设表示捕获能力的系数为e，相当于食饵的自然增长率由r1 降为r1-e，捕食者的死亡率由r2 增为 r2+e 地中海鲨鱼问题地中海鲨鱼问题地中海鲨鱼问题地中海鲨鱼问题 设战前捕获能力系数e=0.3, 战争中降为e=0.1, 则战前与战争中的模型分别为:模模 型型 求求 解解:1、分别用m-文件shier1.m和shier2.m定义上述两个方程2、建立主程序shark1.m, 求解两个方程，并画出两种情况下鲨鱼数在鱼类总数中所占比例 x2(t)/[x1(t)+x2(t)] 实线为战前的鲨鱼比例，“*”线为战争中的鲨鱼比例结论：战争中鲨鱼的比例比战前高！结论：战争中鲨鱼的比例比战前高！ 再见！。