数学实验第二次作业.docx

3.问题：小型火箭初始质量为1400kg，其中包括1080kg燃料，火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度及火箭到达最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。模型：设速度为v，根据牛顿第二定律，可得微分方程在0t60时，dv/dt=(32000-0.4*v^2-9.8*(-18*t+1400))/(-18*t+1400)在引擎关闭，火箭上升至最高点之前dv/dt=-(18*v^2+9.8*320)/320;计算方法：编写dv/dt的函数function dv=rocket1(t,v)dv=(32000-0.4*v^2-9.8*(-18*t+1400))/(-18*t+1400);function dv=rocket2(t,v)dv=-(18*v^2+9.8*320)/320;用龙格-库塔方法求这两个常微分方程，再利用梯形公式求出v对t的积分，得到火箭上升的高度。并输出60s时的速度、加速度和高度,以及总的上升高度。ts1=0:0.1:60;v0=0;[t1,v1]=ode45(@rocket1,ts1,v0);H1=trapz(t1,v1)dv=(32000-0.4*v1(601)^2-9.8*(-18*t1(601)+1400))/(-18*t1(601)+1400)ts2=60:0.1:71.3;v0=v1(601)[t2,v2]=ode45(@rocket2,ts2,v0);H2=trapz(t2,v2)H=H1+H2dv1=(32000-0.4*v1.^2-9.8.*(-18.*t1+1400))./(-18.*t1+1400);dv2=-(18.*v2.^2+9.8.*320)./320;dv=[dv1;dv2];plot(t,dv),grid,[t]=[t1;t2];[v]=[v1;v2];pause,plot(t,v),grid,h(1)=0;for i=2:length(t);h(i)=trapz(t(1:i),v(1:i));endpause,plot(t,h),grid,实验结果：t=60s时，v=267.2612407732609m/sa=0.914984734734975m/s2H=12.18976913272247km后一段上升高度为0.9255331614101061km总上升高度为13111530229413258km显然此时的加速度为-9.8m/s2高度图像如下：速度图像如下：加速度图像如下：实验结果分析与讨论：从图上可知，60s之后，速度、加速度发生突变，速度在较短的时间里变为0。6.题目：一只小船渡过宽为d的河流，目标是起点A正对着的另一岸B点。已知河水流速v1与船在静水中的速度v2之比为k。（图略）建立描述小船航线的数学模型，求其解析解；设d=100m,v1=1m/s,v2=2m/s,用数值解求渡河所需时间、任意时刻小船的位置及航行曲线，作图，并与解析解比较。若流速v1=0,0.5,1.5,2(m/s),结果将如何。模型：以水流方向为x轴，A到B的方向为y轴，小船在任意时刻的坐标为（x,y,t）;由图容易得到速度与位移之间的关系：初始条件为（0，-d）两式相除得令推出：推出：故有有代入初值得计算方法：首先编写船过河的微分方程，存入boat.mfunction dx=boat(t,x)d=100;v1=1;v2=2;s=sqrt(x(1)^2+(d-x(2))^2);dx=[-v2*x(1)/s+v1;v2*(d-x(2))/s];再编写解析解方程：function x=boat2(y,k) x=0.5*(-0.01)^(-k)*y.^(1-k)-0.5*(-0.01)^k*y.^(1+k); ts=0:0.01:100;x0=[0,-100];option=odeset(reltol,1e-6,abstol,1e-9);[t,x]=ode15s(@boat,ts,x0,option,1,2);%这里需使用解刚性方程得龙格—库塔公式计算plot(t,x),gridgtext(x(t))gtext(y(t))pauseplot(x(:,1),x(:,2)),grid,y=0:-0.01:-100; k=0.5;x=boat2(y,0.5); pause,plot(x,y),grid;计算结果：数值解x,y-t图像：解析解x,y-t图像：x-y图像：当t=66.66s时，x=0.0069734m，y= -1.9451e-006m，可认为小船已经达到目的地B点。结果分析与讨论：可以发祥数值解、解析解的差距不是很大。将