数学实行常微分方程数值解.docxVIP

数学实验报告 实验 4 常微分方程数值解 实验 4 常微分方程数值解 分 1 黄浩 2011011743 一、 实验目的 掌握用 MATLAB软件求微分方程初值问题数值解的方法； 通过实例学习用微分方程模型解决简化的实际问题； 了解欧拉方法和龙格 - 库塔方法的基本思想和计算公式，及稳定性等概念。 二、 实验内容 《数学实验》第一版（问题 2） 问题叙述： 小型火箭初始重量为 1400kg，其中包括 1080kg 燃料。火箭竖直向上发射时燃料燃烧率为 18kg/s ，由此产生 32000N的推力，火箭引擎在燃烧用尽时关 闭。设火箭上升时空气阻力正比于速度的平方，比例系数为 0.4kg/m ，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 模型转换及实验过程： （一）从发射到引擎关闭 设火箭总质量为 m，上升高度为 h，瞬时速度为 v，瞬时加速度为 a，由燃料燃烧时间 t=60s ，可列如下的方程组： ??(??) = ??0 ? ????= 1400 ? 18?? ? (??) = ??(??) 其中 ??∈[0,60] ??(??) = ??(??) ?? ? ?? 2 推 阻 32000 ? 0.4??(??) ? 9.8 ??( )?= ? ??= 1400 ? 18?? ??(??) 因此，上述方程为二元常微分方程组，选择 t 为自变量， h 和 v 为因变量 进行分析。 初值条件： ?(0) = 0 , ??(0) = 0 对上述模型，使用 ode45() 函数求数值解（程序见四 .1 、四 .2 ），结果如下： 1 数学实验报告 实验 4 常微分方程数值解 x h(t) v(t) a(t) 0.00 0.00 0.00 13.0571 1.00 6.57 13.19 13.3045 2.00 26.44 26.58 13.4533 3.00 59.76 40.06 13.4972 4.00 106.57 53.54 13.4331 5.00 166.79 66.89 13.2613 6.00 240.27 80.02 12.9853 7.00 326.72 92.83 12.6122 8.00 425.79 105.22 12.1520 9.00 536.99 117.11 11.6169 10.00 659.80 128.43 11.0213 11.00 793.63 139.14 10.3800 12.00 937.85 149.18 9.7083 13.00 1091.79 158.55 9.0209 14.00 1254.71 167.23 8.3309 15.00 1425.93 175.22 7.6502 16.00 1604.83 182.55 6.9901 17.00 1790.78 189.22 6.3593 18.00 1983.13 195.27 5.7646 19.00 2181.24 200.75 5.2095 20.00 2384.47 205.70 4.6946 21.00 2592.36 210.18 4.2220 22.00 2804.52 214.19 3.7943 23.00 3020.56 217.79 3.4120 24.00 3240.08 221.01 3.0730 25.00 3462.65 223.92 2.7726 26.00 3687.88 226.56 2.5044 27.00 3915.58 228.97 2.2677 28.00 4145.60 231.14 2.0633 29.00 4377.76 233.11 1.8898 30.00 4611.86 234.91 1.7433 31.00 4847.68 236.57 1.6178 32.00 5085.02 238.14 1.5062 33.00 5323.85 239.61 1.4095 34.00 5564.11 240.99 1.3293 35.00 5805.77 242.28 1.2650 36.00 6048.72 243.50 1.2139 37.00 6292.87 244.68 1.1708 38.00 6538.11 245.83 1.1303 2 数学实验报告 实验 4 常微分方程数值解 39.00 6784.48 246.96 1.0947 40.00 7031.96 248.05 1.0663 41.00 7280.54 249.10 1.0456 42.00 7530.19 250.12 1.0308 43.00 7780.85 251.14 1.0178 44.00 8032.49 252.