第八章常微分方程的数值解法精选.ppt

第八章 常微分方程的数值解法 第二节 欧拉方法 第二节 龙格—库塔法 在上一节中，我们得到了一些求微分方程近似解的数值方法，这些方法的局部截断误差较大，精度较低，我们希望得到有更高阶精度的方法。 一阶龙格 —库塔方法 如果以y(x)在xi处的斜率作为y(x)在 [xi，xi+1]上的平均斜率k*，即 二阶龙格 —库塔方法 在[xi，xi+1]上取两点xi，xi+p（0 p≤1）的斜率值k1，k2的线性组合λ1k1+λ2k2作为 k*的近似值(λ1、λ2为待定常数)，此公式 一般形式可写成 这就是二阶龙格 —库塔法公式。 三阶龙格 —库塔公式 为了提高精度，考虑在[xi，xi+1]上 取三点xi，xi+p，xi+q的斜率值分别为k1，k2，k3，将k1，k2，k3的线性组合作为平均斜率k*的近似值，其中 k1 k2 k* 公式中含有三个参数λ1，λ2和p，如果我们适当选取参数的值，可以使公式的局部截断误差为O(h3)。 对k1和k2作泰勒展开 代入（8-7）得 （*） 又 y (x)在xi处的二阶泰勒展开式为 当x = xi+1时， ，有 （**） （**） 比较（*）与（**）的系数即可发现， 要使公式（7-7）的局部截断