追击问题.pptVIP

追击问题 我缉私雷达发现，距离c处有一走私船正以匀速a沿直线行驶，缉私舰立即以最大速度（匀速b）追赶。 若用雷达进行跟踪，保持船的瞬时速度方向始终指向走私船， 缉私舰的运动轨迹是怎样的？是否能够追上走私船？ 如果能追上，需要用多长时间？ 应用、思考和练习(追击问题，如果雷达失效) 当缉私舰雷达发现d处有一走私船后，雷达突然损坏 若假定走私船作匀速直线运动（但不知方向），且缉私舰艇速度v大于走私船速度a， 则缉私舰应采用什么样的航行路线，不管走私船从哪个方向逃跑，都能追捕上它？ * * 1：数学建模 为找出x与y的关系，我们设法消除变量t。对第一个式子关于x求导。 由于： （式子中的－是由于s随x的增大而减小） 从而可以得到二阶微分方程的初值问题 其中： 即为缉私艇追击走私船的数学模型。 2 求解析解 前面给出的二阶微分方程初值问题是属于可降阶方程，故令： 可得微分方程： 这又是一个一阶微分方程的初值问题，此方程 的解依赖于参数r 1）当r1（即走私船的速度a小于缉私艇的速度b）时，方程的解为： 即为缉私艇的追赶的路线函数 当x＝0时，缉私艇追赶上走私船，此时走私船走过的距离： 追赶时间： 如图:c=3km, a=0.4km/min, 分别取b=0.6,0.8,1.2km/min时，缉私艇追赶路线图形。 追赶时间分别为： 9.0000 5.0000 2.8125 2）：当r=1时 即走私船的速度a等于缉私艇的速度b时，方程的解为： 式中当x（0）时，y（inf）说明这种情况下，缉私艇不可能追上走私船。 3）当r1时 此时走私船的速度a大于缉私艇的速度b 方程的解为： 式中当x(0)时，y（inf）这种情况下缉私艇不可能追上走私船。 3 用MATLAB软件求解 对降阶的微分方程 可用MATLAB求解 syms x y r c; dsolve(Dy=1/2*((x/c)^r-(c/x)^r), y(c)=0,x) ans = 1/2*(x/c)^(r+1)*c/(r+1)+ 1/2/(-1+r)*x*(c/x)^r-c*r/(r^2-1) 4 用MATLAB软件求数值解 我们可以根据MATLAB提供的数值求解命令求初值问题的数值解，多数情况要求我们利用MATLAB软件自己编程计算问题的数值解，这里我们使用二三阶龙格库塔算法和四五阶龙格库塔算法求解。 取 c=3, a=0.4, b=0.8, r=a/b=0.5 function y=zx(t,y) y=0.5*((t/3))^0.5-(3/t)^0.5; ode23(zx,3,0.0005,0) %注意该命令中终值0.0005可以选择的 追击路线图 ode45(‘zx,3,0.0005,0) 追击路线图 5 计算机仿真算法 一般情况下建立动态微分方程是比较困难的，我们可以用计算机仿真法对系统进行分析研究。 所谓计算机仿真就是利用计算机对实际动态系统的结构和行为进行编程，模拟和计算，由此预测系统的行为效果。 下面我们对缉私艇追击走私船的过程进行计算机仿真。 实际过程： 一步步的进行模拟，动态过程。 取时间步长 则在时刻 缉私艇的位置： 从而： 计算过程中，直线段代替了曲线段，当时间段划分较细时，计算近似准确。 仿真算法： 第一步: 设置时间步长 速度a,b及初始位置 第二步: 计算缉私艇在时刻 坐标 对走私船在时刻 坐标 第三步: 计算缉私艇与走私船这两个动点之间 的距离： 如果 小于事先给定的距离，则认为缉私 艇已经追上了走私船，退出循环，否则让时间产生一个步长，进入下一循环。 第四步: 当从述循环退出后，将点列 与 分别绘制 两条曲线即为缉私艇和走私船走过的轨迹。 实验过程: 取 c=3, a=0.4, b=0.8, r=a/b=0.5 c=3;a=0.4/60;b=0.8/60; jstxb=[];jstyb=[];zscxb=[];zscyb=[]; d=0.01;dt=2;t=0; jstx=3;jsty=0;zscx=0;zscy=0; while (sqrt((jstx-zscx)^2+ … (jsty-zscy)^2)=d) t=t+dt; jstx=jstx-b*dt*jstx/ … sqrt(jstx^2+(a*t-jsty)^2); jstxb=[jstxb,jstx]; jsty=jsty+b*dt*(a*t-jsty)/ sqrt(jstx^2+(a*t-jsty)^2); jstyb=[jstyb,jsty]; zscy=a*t; zscyb=[zscyb,zscy]; end zscxb=zeros(length(zscyb)); plot(jstxb,jstyb,o,