缉私艇追走私船模型实验解析.pptxVIP

缉私艇追走私船模型实验解析;2、了解微分方程数值解思想，掌握两种简单的微分方程数值解方法。;微分方程数值解(理论). 用MATLAB软件求微分方程解析解、数值解. 微分方程模型实验:缉私艇追赶走私船. 问题的进一步扩展.;1. 微分方程数值解;数值解方法: 欧拉方法、欧拉两步法、改进的欧拉方法;所以，过初始点A0(x0,y0)，做切线;这样依次类推得到上述微分方程的解在x=xn+1处 函数值的近似计算公式，;特别地，如果节点之间是等分的，则有;将上两式算术平均即得梯形公式：;2）欧拉两步法;3）改进的欧拉公式;步骤2 对预测值再用梯形公式校正一次得yn+1，称为校正值;再解释：;2. 用MATLAB软件求微分方程的解析解;3. 编程计算微分方程数值解;（1）直接在命令窗执行命令 dsolve('Dy=3/x*y+x^3*(exp(x)+cos(x))-2*x', 'y(pi)=(exp(pi)+2/pi)*pi^3','x') ans= x^3*exp(x)+x^3*sin(x)+2*x^2 ;szy_eu=[]; y=(exp(pi)+2/pi)*pi^3; for x=pi:0.1:2*pi y=y+0.1*(3/x*y+x^3*(exp(x)+cos(x))-2*x); szy_eu=[szy_eu,y]； end szy_eu t=pi:0.1:2*pi; f=t.^3.*exp(t)+t.^3.*sin(t)+2*t.^2; plot(t,szy_eu,'r-','linewidth',2) hold on plot(t,f,'b*--');第18页/共41页;x=pi:0.1:2*pi; h=0.1; y=(exp(pi)+2/pi)*pi^3; szy_imeu=[y]; for i=1:length(x)-1 jzj=y+h*(3/x(i)*y+x(i)^3*(exp(x(i))+cos(x(i)))-2*x(i)); yy=y+h/2*(3/x(i)*y+x(i)^3*(exp(x(i))+cos(x(i)))- 2*x(i)+3/x(i+1)*jzj+x(i+1)^3*(exp(x(i+1))+cos(x(i+1)))-2*x(i+1)); szy_imeu=[szy_imeu,yy]; y=yy; end szy_imeu plot(x,szy_imeu,'r-','linewidth',2) f=x.^3.*exp(x)+x.^3.*sin(x)+2*x.^2; hold on plot(t,f,'b*--');第20页/共41页;4. 用MATLAB软件求微分方程数值解;例3;5. 模型实验: 缉私艇追赶走私船;(1) 建立模型;可得缉私艇追击走私船过程的数学模型：;(2) 模型求解;--zhao1203.m;取c=3千米，a=0.4千米/分钟，分别取b=0.6,0.9, 1.2千米/分钟时，缉私艇追赶路线的图形为;第29页/共41页;定义m文件函数： function y=zx(t,y) y=0.5*((t/3)^0.5-(3/t)^0.5) ;6 动态仿真;走私船初始位在点(0,0)，方向为y轴正方向，;时间步长;仿真算法： ;第三步：计算缉私艇与走私船这两个动点之间的距离： ;取c=3千米，a=0.4千米/分钟，b=0.8千米/分钟，r=a/b=0.5 显示船与艇行进路线程序;结果分析;上机任务;1. 目的 2. 内容（对每一???）： 模型（对应用题）；算法设计；计算结果；结果分析；附程序（必要时加说明语句）。 3. 收获和建议。 于第17周周五（6月15日）之前提交实验报告！ ;THANK YOUR ATTENDING;谢谢观看！