数学实验—常微分方程数值解.docx

. .. 实验4 常微分方程数值解 分1 黄浩 2011011743 实验目的 掌握用MATLAB软件求微分方程初值问题数值解的方法； 通过实例学习用微分方程模型解决简化的实际问题； 了解欧拉方法和龙格-库塔方法的基本思想和计算公式，及稳定性等概念。 实验内容 《数学实验》第一版（问题2） 问题叙述： 小型火箭初始重量为1400kg，其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃烧用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 模型转换及实验过程： （一）从发射到引擎关闭 设火箭总质量为m，上升高度为h，瞬时速度为v，瞬时加速度为a，由燃料燃烧时间t=60s，可列如下的方程组： mt ht=v vt a 因此，上述方程为二元常微分方程组，选择t为自变量，h和v为因变量进行分析。 初值条件：h 对上述模型，使用ode45()函数求数值解（程序见四.1、四.2），结果如下： 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 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.15 1.0024 45.00 8285.12 253.16 0.9876 46.00 8538.75 254.15 0.9763 47.00 8793.39 255.12 0.9696 48.00 9049.01 256.07 0.9663 49.00 9305.58 257.03 0.9624 50.00 9563.08 257.99 0.9527 51.00 9821.52 258.95 0.9412 52.00 10080.93 259.90 0.9337 53.00 10341.30