微分方程应用案例

案例

导弹追踪问题

设位于坐标原点的甲舰向位于x轴上点A(1，0)处的舰发射导弹，导弹头始终对准了舰.如果了,舰以最大的速度Vo(是常数)沿平行于y轴的直线行驶，导弹的速度是5vo，求导弹运行的曲线方程.又乙舰行驶多远时，导弹将它击中 ?

解法

解法一( 解析法
假设导弹在t时刻的位置为P(x(t)，y(t))，乙舰位于Q(,V)由于导弹头始终对准乙舰，故此时直线PQ就是导弹的轨迹曲线弧OP在点P处的切线

p%t_y即有1-X即(1)vot = (1-x)y'+y
又根据题意，弧OP的长度为 AQ 的5倍,即 1+ ydx = 5v,t (2)

由(1),(2)消去t整理得模型
(1-x)y"=(3初值条件为: y(0)=0 y'(0) =0解即为导弹的运行轨迹
2当X=1时y=即当乙舰航行到点(1,处时被导弹击中被击中时间为:t = -若 v-1，则在 -0.21 处被击中24vo

导弹轨迹程序 :
clear
X-0:0.01:1;
y--5*(1-x).^(4/5)/8+5*(1-x).^(6/5)/12+5/24

plot(x,y,

解法二( 数值解
令yi-y;y,-y，，将方程(3) 化为一阶微分方程组。y= Y(-y"= vI+yh=VI+Y/0-3)

1.建立m-文件fun.m
function dy=fun(x,y)dy=zeros (2,1) ;dy(1)=y (2) .dy(2)=1/5*sgrt(1