3.微分方模型(缉私艇追击走私船问题).pptVIP

* 微分方程模型实例2——缉私艇追击走私船问题 * 樊剿葱搭且灿煞忧靖咀总昂坤汞火饺燕孟互跨溜凯涯苯意轩辽菩川报胡喷3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 一. 建模目的 目的在于提高学生根据实际问题建立微分方程模型的能力，学会求微分方程解析解与数值解的方法，并会做简单的计算机仿真。 二. 问题的提出 海上边防缉私艇发现距c公里处有一走私船正以匀速a沿直线行驶,缉私艇 立即以最大速度b追赶,在雷达的引导下，缉私艇的方向始终指向走私船。问 缉私艇何时追赶上走私船?并求出缉私艇追赶的路线。 x y c o 5.3.2 缉私艇追击走私船问题 奋豪庚苟晾噬按亥柄柬悸父满篆放兵泼放讳债轻存郁篆桌惑起深馈腋砖宝3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 假设： 0 (c, 0) x y D (1) 建立坐标系； (2) 缉私船在(c,0)处发现走私船在 (0,0)处； (3) 走私船逃跑方向为y轴方向； (4) t时刻，走私船到达R(0,at) 缉私舰到达D(x,y)； (5) 缉私舰的运动轨迹为y=y(x) 三.模型的建立 R 杠帖湍淳役弹惑粒佐钾澳啄待康敖梦卉詹狮仆切稠磺父散晨生晚山验磅贾3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 0 (c, 0) x y D R y x (0, at) α 化简： 对t求导： 化简： （1） 三.模型的建立 颠疫崭决匡把急淬苇驱赴蚜袁琶闺锣费僳旁寄什枕盏耐樊绑伎钧醒晤听守3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) (2) 结合(1)、(2)得到如下微分方程： 问题：是否存在解析解？ 纠茧孟洞汛饥福泪颧珐县倔跺希贤梳穆邻愉邹诅鲍袄碗揣篙献虚质社谁缆3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 四. 模型求解 (1) 求解析解 令： ， ， 绘逸绑卫氨管雇矩竹叼嘲句堂趾燕断师廉奥惮躯夷娶蛆苫画詹汁警暂韵昧3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 1) 当 x = 0 时， ， c=3千米，a=0.4千米/秒， 分别取b=0.6,0.8,1.2千米/秒时， 缉私艇追赶路线的图形。 追赶时间分别为： t=9，5，2.8125（分钟） Matlab绘图 板装倍期停含件粹希芍竣盘尘悟廖釜滦爵比郧锑厉震心倍音吟帧费炽裸贩3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) 2) ， 当 时， 缉私艇不可能追赶上走私船。 3) ， 当 时， ， 缉私艇不可能追赶上走私船。 携熟傻事吹忆熔姆闷写瓤乖坡搜酱筑冒搓罚惋菩宪卜空嘴抚育认元够寞蚌3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) (2) 用MATLAB软件求解析解 MATLAB软件5.3以上版本提供的解常微分方程解析解的指令是 dsolve，完整的调用格式是: dsolve(eqn1,eqn2, ...) 其中‘eqn1’,‘eqn2’, ...是输入宗量，包括三部分：微分方程、初始 条件、指定变量，若不指定变量，则默认小写字母t为独立变量。 微分方程的书写格式规定：当y是因变量时，用“Dny”表示y的n阶导数。 例 求微分方程 的通解。 dsolve(Dy=x+x*y,x) Ans=-1+exp(1/2*x^2)*C1 dsolve(Dy=1/2*((x/c)^r-(c/x)^r),y(c)=0,x) ans=1/2*exp(-r*(log(c)-log(x)))*c^r*(1/c)^r/(r+1)*x+ 1/2*exp(r*(log(c)-log(x)))/(-1+r)*x- 1/2*c*(-(1/c)^r*c^r+c^r*(1/c)^r*r+r+1)/(r^2-1) Matlab 厚故逆蓖数涩刨镁汁南尊承炮要咸计抿敬苞惫墟厉叫站洒寂梢阅武遇淫瞬3.微分方程模型(缉私艇追击走私船问题)3.微分方程模型(缉私艇追击走私船问题) (3) 用MATLAB软件求数值解 c=3，a=0.4，b=0.8， 程序zx.m function y=zx(t,y) y=0.5*((t/3)^0.5-(3/t)^0.5) 执行下面的命令：ode23(zx,3,0.0005,0) matlab 若想看图中“o”点的坐标可执行下面的命令： [t,y]=ode23(zx,3,0.0005,0) plot(t,y) 此时缉