西安交通大学matlab数学实验课件6.pptVIP

* * * 实验6 缉私艇追击走私船 1 实验目的 (1)提高学生根据实际问题建立微分方程模型的能力； (2)学习求微分方程解析解与数值解的方法; (3)学习简单的计算机仿真方法。 海上边防缉私艇发现距c公里处有一走私船正以匀速a沿直线行驶,缉私艇立即以最大速度b追赶,在雷达的引导下，缉私艇的方向始终指向走私船。问:缉私艇何时追赶上走私船?并求出缉私艇追赶的路线。 x y c o 2 实验问题 3 建立模型 x c o y 走私船初始位置在点(0,0)， 行驶方向为y轴正方向， 缉私艇的初始位置在点(c,0)， 缉私艇行驶的历程为s 。 在时刻t： 缉私艇到达点 走私船的位置到达点 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阶导数。 完整的调用格式是: dsolve(eqn1,eqn2, ...) 其中‘eqn1’,‘eqn2’, ...是输入宗量，包括三部分：微分方程、初始条件、指定变量，若不指定变量，则默认小写字母t为独立变量。 （2）用MATLAB软件求解析解 例 求微分方程 的通解。 指令：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) （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) 此时缉私艇的位置坐标是（0.0005，1.9601） 执行下面的命令: ode45(zx,3,0.0005,0) 若想看图中“o”点的坐标可执行下面的命令： [t,y]=ode45(zx,3,0.0005,0) plot(t,y) 此时缉私艇的位置坐标是（0.0005，1.9675 ） 若想看图中“o”点的坐标可执行下面的命令： [t,y]=ode23(zx,3,0.0005,0) plot(t,y) ode23(function,domain,initial condition) MATLAB软件仿真法 当建立动态系统的微分方程模型很困难时， 我们可以用计算机仿真法对系统进行分析研究。 所谓计算机仿真就是利用计算机对实际动态系统的结构和行为进行编程、模拟和计算，以此来预测系统的行为效果。 MATLAB软件仿真法 x c o y 走私船初始位在点(0,0)，方向为y轴正方向， 缉私艇的初始位在点(c,0)， 追赶方向可用方向余弦表示为： 走私船的位置： 缉私艇的位置： o 缉私艇的位置： 则 时间步长 仿真算法： 第二步： 计算动点缉私艇D在时刻 时的坐标 计算走私船R在时刻 时的坐标 第一步：设置时间步长 , 速度a, b及初始位置 o 第三步：计算缉私艇与走私船这两个动点之间的距离： 根据事先给定的距离，判断缉私艇是否已经追上了走私船，从而判断退出循环还是让时间产生一个步长，返回到第二步继续进入下一次循环； 第四步：当从上述循环退出后，由点列 和 可分别绘制成两条曲线 即为缉私艇和走私船走过的轨迹曲线。 程序zhuixiantu.m c=3; a=0.4/60; b=0.8/60; jstxb=[];jstyb=[];zscxb=[];zscyb=[]; d=0.01;dt=2;t=0; jstx=c;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)^