实验四_常微分方程数值解new

《实验四_常微分方程数值解new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验四常微分方程数值解材料系材932009011976邓陟一．实验目的1.用MATLAB软件掌握求微分方程数值解的方法，并对结果作初步分析。2.通过实例学习用微分方程模型解决简化的实际问题。二．实验内容1.小型火箭初始质量为1400kg，其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭达到最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。可以根据牛顿第二定律Fma列

2、出火箭不同飞行时期满足的运动方程。1)燃料未用尽前(0t60)火箭飞行时受到燃料燃烧产生的推力、自身的重力以及空气阻力（向下）。列出微分方程dhvdtdvFkv2pga=dtmutmut00hv00,002其中F为燃料燃烧产生的推力，g为重力加速度（取9.8m/s），k为空气阻力p比例系数，m为火箭初始质量，u为燃料燃烧率。0编写M文件rocket1.mfunctiondx=rocket1(t,x)Fp=32000;g=9.8;k=0.4;m0=1400;u=18;dx=[x(2);Fp/(m0-u*t)-g-k*x(2

3、)^2/(m0-u*t)];编写程序ts=0:0.1:60;x0=[0,0];[t,x]=ode45(@rocket1,ts,x0)得到引擎关闭瞬间火箭的高度h6012190m，速度v60267.26m/s，2带入微分方程算得加速度a0.92m/s。2)燃料用尽瞬间至火箭达到最高点火箭飞行时受到自身的重力以及空气阻力（向下）。列出微分方程dhvdt2dvkvga=dtmhv6012190,60267.262其中g为重力加速度（取9.8m/s），k为空气阻力比例系数，m为燃料用尽后的火箭质量。编写M文件rocket2

4、.mfunctiondx=rocket2(t,x)g=9.8;k=0.，4;m=320;dx=[x(2);-g-k*x(2)^2/m];编写程序ts=60:0.1:75;x0=[12190,267.26];[t,x]=ode45(@rocket2,ts,x0)算得在t71.3s时v0，此时火箭的高度h71.313116m，加速度2ag9.8m/s。3)最高点变加速下落到最终匀速下落至地面火箭飞行时受到自身的重力以及空气阻力（向上）。列出微分方程dhvdt2dvkvga=dtmhv71.313116,71.302

5、其中g为重力加速度（取9.8m/s），k为空气阻力比例系数，m为燃料用尽后的火箭质量。为画出火箭高度、速度、加速度随时间变化的图形，须描述火箭整个的运动过程。编写M文件rocket.mfunctiondx=rocket(t,x)Fp=32000;g=9.8;k=0.4;m0=1400;u=18;m=320;ift<=60dx=[x(2);Fp/(m0-u*t)-g-k*x(2)^2/(m0-u*t)];endift>60&t<=71.3dx=[x(2);-g-k*x(2)^2/m];endift>71.3dx=[x(2);-g+k*x(2)^2/m];end编写程序

6、ts=0:0.1:225.2;x0=[0,0];[t,x]=ode45(@rocket,ts,x0);fori=1:2253ift(i)<=60a(i)=(32000-0.4*x(i,2).^2)/(1400-18*t(i))-9.8;endift(i)>60&t(i)<=71.3a(i)=-0.4*x(i,2).^2/320-9.8;endift(i)>71.3a(i)=0.4*x(i,2).^2/320-9.8;endendsubplot(3,1,1),plot(t,x(:,1));subplot(3,1,2),plot(t,x(:,2));subplot(3,

7、1,3),plot(t,a);高度、速度、加速度随时间变化的图形2.一只小船渡过宽为d的河流，目标是起点A正对着的另一岸B点。已知河水流速v1与船在静水中的速度v之比为k。2(1)建立描述小船航线的数学模型，求其解析解；(2)设d100m，v1m/s，v2m/s，用数值解法求渡河所需时间、任意时刻12小船的位置及航行路线，作图，并与解析解比较；(3)若流速v0,0.5,1.5,2(m/s)，结果将如何。11)建立数学模型以B点为坐标系原点建立直角坐标系，A点坐标为0,d。设小船航线上任意一点的坐标为xy,，v与v夹角为，列出微分方程21dx