第六讲常微分方程初值问题的数值解法2合编.ppt

* §6.3 Runge-Kutta法 考虑改进Euler法 通常写成如下的预报—校正公式 预报—校正公式是由梯形公式和Euler公式复合而成， 它具有二阶代数精度。 一、Runge-Kutta方法的基本思想 对于常微分方程的初值问题 的解 即 引入记号 则 改进的欧拉公式计算时进一步修改为预报—校正公式 推广上面的格式 其中 这就是龙格—库塔方法的思想。 如何确定其中的待定系数？ 通过Taylor展开式求上述式中的待定系数。 二、Runge-Kutta公式的推导 此时，有 利用二元函数的Taylor展开形式 再利用 得到 比较上面的两个式子，得 满足上式的解有无穷多，每一公式都是二阶精度。统称为二阶龙格—库塔方法。 若取 则正好是预报—校正公式。 若取 则得到 同样的推导方法可以构造出三阶、四阶的Runge-Kutta 公式。 一般常用的三阶Runge-Kutta公式为 局部截断误差为 ，三阶精度 实际中更常用的是下列四阶Runge-Kutta公式，又称为经典的龙格—库塔公式。 例1 用经典的龙格-库塔法计算上节例题，取步长 解： 从而 利用上面公式可计算出 欧拉法 预报-校正法 准确解 1 1.183216 1.341641 1.483240 1.6124