系统仿真与数值计算综合实验例题..docVIP

1. 小型火箭初始质量为900千克，其中包括600千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生30000牛顿的恒定推力。当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数为0.4（千克/米）。重力加速度取9.8米/秒2. 建立火箭升空过程的数学模型（微分方程）； 求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时间和高度。 解： 建立模型： ——时刻的火箭高度； =30000（牛顿）——火箭推力，当时间40秒时，T=0； ——火箭飞行过程中的质量，40秒时，千克 =900（千克）——火箭初始质量； =600（千克）——燃料质量； =15（千克/秒）——燃料的燃烧速率； =0.4（千克/米）——空气阻力系数； =9.8（米/秒）——重力加速度 由能量守恒定律，可得到火箭飞行过程的方程： 这是一个初值问题，初始条件为 设，则问题化为求下列微分方程组的初值问题： 关闭引擎时（秒），所求的是此时火箭的高度；速度； 加速度，及火箭到最高点的时间和高度。具体的Matlab程序如下： 首先建立微分方程的的m文件： function y=huojian(t,x) k=0.4;g=9.8;m0=900;T=30000;m=m0-15*t; if t40 T=0; m=300; end y=[x(2),-(k/m)*x(2)^2+T/m-g]; 主程序： %feixing.m k=0.4;g=9.8;m0=900;T=30000; x0=[0,0]; ts=0:1:55; [t,x]=ode23(huojian,ts,x0); [t,x(:,1)] %------ a=[t,x]; x40=a(41,2) %燃料用尽时的高度 v40=a(41,3) %燃料用尽时的速度 a40=-(k/300)*v40^2+T/300-g %燃料用尽时的加速度 %------- xmax=max(x(:,1)) %火箭到达最高点的高度 subplot(2,1,1),plot(t,x(:,1)),title(altitude) subplot(2,1,2),plot(t,x(:,2)),title(speed) 运行结果为： 1.0e+003 * 0 0 0.0010 0.0118 0.0020 0.0475 0.0030 0.1067 0.0040 0.1889 0.0050 0.2927 0.0060 0.4168 0.0070 0.5592 0.0080 0.7180 0.0090 0.8914 0.0100 1.0770 0.0380 7.8051 0.0390 8.0631 0.0400 8.3224 0.0410 8.5400 0.0420 8.7004 0.0430 8.8242 0.0440 8.9218 0.0450 8.9994 0.0460 9.0607 0.0470 9.1083 0.0480 9.1437 0.0490 9.1681 0.0500 9.1822 0.0510 9.1864 0.0520 9.1807 0.0530 9.1651 0.0540 9.1390 0.0550 9.1018 x40 =8322.4 v40 = 254.1728 a40 = 4.0616 xmax =9186.4 关闭引擎时（秒），此时火箭的高度h==8322.4米，速度v==254.1728米/秒，加速度为a== 4.0616米/秒，火箭到最高点的时间=51秒，高度=9186.4米。 2．小型火箭初始质量为1200千克，其中包括900千克燃料。火箭竖直向上发射时燃料以15千克/秒的速率燃烧掉，由此产生40000牛顿的恒定推力．当燃料用尽时引擎关闭。设火箭上升的整个过程中，空气阻力与速度平方成正比，比例系数记作k，火箭升空过程的数学模型为 其中为火箭在时刻t的高度，m=1200-15t为火箭在时刻t的质量，T（=30000牛顿）为推力，g (=9.8米/秒