计算方法件第八章常微分方程初值问题的数值解法.ppt

复习泰勒公式 由于二级方法 ，方程也可写为 要使Rn尽量小，应首先让h，h2项的系数为零，即 上述方程所确定的解都能使二级 R-K方法成为一个二阶方法。 此时Rn为 为使Rn最小，应令 从而有 局部截断误差为： 故有二阶R-K方法为： 该方法又称为二阶Heun（休恩）方法。 Euler预估-校正格式 若取 Remark1：我们可以构造无穷多个二级R-K方法，这些方法的截断误差均为O(h3),即都是二阶方法。其中二阶Heun方法是截断误差项数最少，且允许f 任意变化的情况下截断误差最小的二阶方法。 Remark2：二级R-K方法不可能达到三阶 Remark3：同样可构造其他阶的R-K方法，它们都有无穷多组解，且三级R-K方法阶数不超过3，四级R-K方法阶数不超过4。 Remark4：更高阶的方法由于计算量较大，一般不再采用。 标准（经典）四级四阶R-K公式 关于R-K方法计算量的讨论 二阶R-K方法需计算两个函数值，四阶R-K方法需计算四个函数值，但精度要比二阶方法高出两阶。因此，要达到同样的精度，用低阶方法需步长取得比较小，但若用高阶方法则可以将步长取得大一些，从而降低计算量。 四阶经典R－K方法的稳定性条件是 关于R-K方法稳定性的讨论 二阶R－K方法的稳定性条件是 线性多步法的基本思想：如果充分利用前面多步的信息来预测yn+