第3章03龙格-库塔方法.pptVIP

* * * * * * * * * * * * * * * * * * * * * * 第三章 常微分方程的差分方法 3.1 欧拉方法 3.2 改进的欧拉方法 3.3 龙格-库塔方法 3.4 亚当姆斯方法 3.5 收敛性与稳定性 3.6 方程组和高阶方程的情形 3.7 边值问题 3.3 龙格-库塔(Runge-Kutta)方法 3.3.1 龙格-库塔方法的设计思想 3.3.2 二阶龙格-库塔方法 3.3.3 三阶龙格-库塔方法 3.3.4 四阶龙格-库塔方法 3.3.5 变步长的龙格-库塔方法 3.3 龙格-库塔（Runge-Kutta)法 龙格-库塔方法的设计思想 龙格-库塔方法的设计思想 龙格-库塔方法的设计思想 3.3.4 四阶龙格-库塔方法 仿照上述的讨论，用四个不同点上的函数值的线性组合,可导出一族四阶龙格－库塔格式: 四阶龙格－库塔公式每步要四次计算函数值，具有四阶精度，局部截断误差是O(h5) . 四阶龙格－库塔格式 下列经典格式是其中常用的一种: 经典龙格－库塔公式每步要四次计算函数值，具有 四阶精度，局部截断误差是O(h5) . 四阶经典龙格－库塔法的框图 框图 3.3.5 变步长的龙格-库塔方法 对每一步: 步长越小,截断误差越小; 但步长越小,完成的步数越多,计算量越大,可能造成舍入误差的严重积累. 选择步长的原则: 如何衡量和检验计算结果的精度; 如何依据所判定的精度来调整步长. 四阶经典龙格-库塔方法步长折半前后的局部截断误差 可以通过检查步长折半前后两次计算结果的偏差 来判断所选取的步长是否合适. 对于给定的精度ε,如果δ ε,则反复将步长折半进行计算,直到δ ε为止,这时取步长折半后的新值作为结果; 如果δ ε,则反复将步长加倍进行计算,直到δε为止,这时取步长加倍前的老值作为结果. 变步长方法表面上看为了选择步长,每一步的计算量增加了,但总体考虑往往是合算的. 例1 例2 用经典Runge-Kutta方法求解 取 h = 0.1计算结果见下表 解: ? Euler方法 预估－校正方法 经典R-K方法 ? ? ? ? ? ? 0.0 1.000000 0.0 1.000000 0.0 10.0 0.1 1.000000 4.8×10－3 1.005000 1.6×10－4 18.2×10－8 0.2 1.010000 8.7×10－3 1.019025 2.9×10－4 11.5×10－7 0.3 1.029000 1.2×10－2 1.041218 4.0×10－4 12.0×10－7 0.4 1.056100 1.4×10－2 1.070800 4.8×10－4 12.4×10－7 0.5 1.090490 1.6×10－2 1.107076 5.5×10－4 12.7×10－7 0.6 1.131441 1.7×10－2 1.149404 5.9×10－4 12.9×10－7 0.7 1.178297 1.8×10－2 1.197210 6.2×10－4 13.1×10－7 0.8 1.230467 1.9×10－2 1.249975 6.5×10－4 13.2×10－7 0.9 1.287420 1.9×10－2 1.307228 6.6×10－4 13.3×10－7 1.0 1.348678 1.9×10－2 1.368514 6.6×10－4 1.3678977 3.3×10－7 整体截断误差 （四级4阶） （二级2阶） （一级1阶） 说明:（1）通过比较可以发现在相同步长下经典Runge-Kutta （2）经典Runge-Kutta方法的计算量是Euler方法的4倍,预估 注:（2）中不与梯形方法比较，因梯形方法与这三种不同类。 方法的结果比Euler方法、梯形方法、预估－校正Euler方法好的 多，即更准确，误差较小。 的步长取 ，预估－校正Euler方法取 ，它们的计算量大致 －校正方法的2倍。 相同，此时，经典Runge-Kutta方法仍比Euler方法、预估－校正 若经典Runge-Kutta方法步长取h，Euler方法 Euler方法好的多。 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *