数学建模综合实验缉私艇走私艇.docVIP

一.实验目的 本综合实验旨在考察及训练学生对微分方程建模及Matlb编程的灵活运用。通过本实验了解数学建模的基本思想，并熟练掌握用数学软件解决数学问题的方法。提高学生的综合能力。 二.实验内容 1、已知微分方程组 满足初始条件． 求上述微分方程组初值问题的特解(解析解)，并画出解函数的图形． 分别用ode23、ode45求上述微分方程组初值问题数值解(近似解)，求解区间为．利用画图比较两种求解器之间的差异． 2、分别用 Euler 折线法和四阶 Runge-Kutta 法求解微分方程初值问题 的数值解(步长h取0.1)，求解范围为区间[0,3] ． 3、海防某部缉私艇上的雷达发现正东方向15海里处有一艘走私船正以20海里/小时的速度向正北方向行驶，缉私艇立即以40 海里/小时的速度前往拦截。用雷达进行跟踪时，可保持缉私艇的速度方向始终指向走私船。建立任意时刻缉私艇的位置和缉私艇航线的数学模型，确定缉私艇追上走私船的位置，求出追上的时间，画出航线图形，并通过改变速度等参数进行讨论。 三. 实验方案（程序设计说明） 第1题：使用 dsolve 函数、ode23、ode45求解器编程求解； 第2题：利用 Euler 折线法和四阶 Runge-Kutta 法的递推公式编程求解； 第3题：实验方案如下： (一)建立模型 以时刻缉私艇位置为原点，正东方向为正轴方向，正北方向为正y轴方向建立直角坐标系，则缉私艇与走私船的初始距离，设缉私艇行驶的路程为，缉私艇航线任一点切线与轴正向夹角为， 缉私艇：速度，初始位置， 时刻位置 走私艇：速度，初始位置， 时刻位置。 由题意有方程： (1) 方程(1)即为缉私艇任意时刻位置的数学模型。 ，两边对求导并化简得： (2) 而，，故 (3) 令，由(2)(3)得： (4) 方程(4)即为缉私艇航线的数学模型。 (二)模型求解 令，，代入(4)， ，积分并化简得： ，取两边倒数有： ，两式相加并化简： ，积分并化简得： ，其中。 上式即为缉私艇的航线方程。 由题意，缉私艇追上走私艇时，有： ，， 按本题条件：，，，得：，， 即缉私艇在半小时后于点恰好追上走私艇。 若取，，得：， 若取，，得：， 程序及运行结果见后文。 四.实验步骤或程序（经调试正确的源程序 第1题： M函数文件fun.m function xx=fun(t,x) xx=[-x(1)-x(2); x(2)-x(1)]; 程序文件prog1.m clear;syms x y t fprintf(特解为:\n) [x,y]=dsolve(Dx+x+y=0,Dy+x-y=0,x(0)=1,y(0)=0,t) figure(1);ezplot(x,y,[0,2]); clear;figure(2); y0=[1;0]; [t,xx]=ode45(fun,[0,2],y0); subplot(1,2,1);plot(xx(:,1),xx(:,2),r-,LineWidth,2); title(ode45) hold on clear;y0=[1;0]; [t,xx]=ode23(fun,[0,2],y0); subplot(1,2,2);plot(xx(:,1),xx(:,2),b-,LineWidth,3); title(ode23) 第2题： Euler折线法：prog2_1.m clear f=sym(y-exp(x)*cos(x)); a=0; b=3; h=0.1; n=(b-a)/h+1; x=0; y=1; szj=[x,y]; for i=1:n-1 y=y+h*subs(f,{x,y},{x,y}); x=x+h; szj=[szj;x,y]; end szj plot(szj(:,1),szj(:,2),or-) Runge-Kutta 法：prog2_2.m clear; f=sym(y-exp(x)*cos(x)); a=0; b=3; h=0.1; n=(b-a)/h+1; x=0; y=1; szj=[x,y]; for i=1:n-1 l1=subs(f,{x,y},{x,y}); l2=subs(f,{x,y},{x+h/2,y+l1*h/2}); l3=subs(f,{x,y},{x+h/2,y+l2*h/2}); l4=subs(f,{x,y},{x+h,y