MATLAB缉私艇追击走私船 .ppt

取c=3千米，a=0.4千米/分钟，b=0.8千米/分钟， c=3; a=0.4/60; b=0.8/60; d=0.01;dt=2;t=0; jstx=c;jsty=0;zscx=0;zscy=0; hold on axis([0,3,0,2]) while (sqrt((jstx-zscx)^2+(jsty-zscy)^2)d) t=t+dt; jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2); jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2); zscy=a*t; plot(jstx,jsty,r+,zscx,zscy,b*) pause(0.2) end jstx,jsty,zscx,zscy,t 程序zhuixiantu.m 历时：296秒， 艇：（0.000041，1.97） 船：（0，1.97） 显示船与艇行进路线程序 c=3; a=0.4/60; b=0.8/60; d=0.01;dt=2;t=0; jstx=c;jsty=0;zscx=0;zscy=0; while (sqrt((jstx-zscx)^2+(jsty-zscy)^2)d) pause(0.01) hold on axis([0,3,0,2]) t=t+dt; jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2); jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2); zscy=a*t; plot(jstx,jsty,rO,zscx,zscy, b*) end * * * * * * * * * * * * * 实验3 缉私艇追击走私船 3.1 实验目的 (1)学会用MATLAB软件求解微分方程的初值问题; (2)了解微分方程数值解的思想,掌握微分方程数值解的方法; (3)学会根据实际问题建立简单微分方程数学模型,提高解决问题能力; (4)了解简单的计算机仿真和数据模拟的基本方法. 3.2 实验问题 海上边防缉私艇发现距c公里处有一走私船正以匀速a沿直线行驶,缉私艇 立即以最大速度b追赶,在雷达的引导下，缉私艇的方向始终指向走私船。问 缉私艇何时追赶上走私船?并求出缉私艇追赶的路线。 x y c o 2.3 建立模型 x c o y 走私船初始位置在点(0,0)， 行驶方向为y轴正方向， 缉私艇的初始位置在点(c,0)， 缉私艇行驶的路程为s 。 在时刻t： 缉私艇到达点 走私船的位置到达点 2.4 模型求解 (1) 求解析解 令： ， ， 1) , 当 x = 0 时， ， c=3千米，a=0.4千米/分， 分别取b=0.6,0.8,1.2千米/分时， 缉私艇追赶路线的图形。 追赶时间分别为： t=9，5，2.8125（分钟） 2) 当 时， ， 缉私艇不可能追赶上走私船。 3） ， ， 当 时， ， 缉私艇不可能追赶上走私船。 （2）用MATLAB软件求解析解 MATLAB软件5.3以上版本提供的解常微分方程解析解的指令是 dsolve，完整的调用格式是: 微分方程的书写格式规定：当y是因变量时，用“Dny”表示y的n阶导数。 例1 求微分方程 的通解。 解 输入命令: y=dsolve(Dy=x+x*y,x) Ans=-1+exp(1/2*x^2)*C1 dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’) 解 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x) 结 果 为 : y =3e-2xsin（5x） 解 输入命令 ： [x,y,z]=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t); x=simple(x) % 将x化简 y=simple(y) z=simple(z) 结 果 为：x = (c1-c2+c3+c2e -3t-c3e-3t)e2t y = -c1e-4t+c2e-4t+c2e-3t-c3e-3t+c1-c2+c3)e2t z = (-c1e-4t+c2e-4t+c1-c2+c3)e2t 在实际生产和科研中