第八章节常微分方程的数值解法幻灯片.pptVIP

k1 k2 k* 公式中含有三个参数λ1，λ2和p，如果我们适当选取参数的值，可以使公式的局部截断误差为O(h3)。 对k1和k2作泰勒展开 代入（8-7）得 （*） 又 y (x)在xi处的二阶泰勒展开式为 当x = xi+1时， ，有 （**） （**） 比较（*）与（**）的系数即可发现， 要使公式（7-7）的局部截断误差满足 ，即要求公式具有二阶 精度只要下列条件成立即可。 （8-8） 满足条件（8-8）的一簇公式统称为 二阶龙格—库塔公式。 特别的，当 塔公式就成为改进欧拉公式。 时，龙格-库 　　改进欧拉公式就是以y(x)在xi和xi+1处 的斜率k1和k2的算术平均值作为y(x)在 [xi，xi+1]上的平均斜率k*来进行计算的。 若取 时，龙格-库塔 公式就称为变形的欧拉公式，其形式为 （8-9） 此处的 就是欧拉方法预报出的中点 处的近似解；而 等于中点的斜率值 则近似 ； 所以公式可以看作用中点斜率近似 代替平均斜率k*，因此，公式（8-9）也 称作中点公式。 三阶龙格 —库塔公式 为了提高精度，考虑在[xi，xi+1]上 取三点xi，xi+p，xi+q的斜率值分别为k1，k2，k3，将k1，k2，k3的线性组合作为平均斜率k*的近似值，其中 xi+q=xi+qh （0p≤q1） 这时的计算公式为 其中 为了预报点xi+q的斜率k3，一种很自然的 想法是用欧拉法预报，即取 但是，这样做效率比较差。因为在区间 [xi，xi+q]内已有两个斜率可以使用，所以把k1， k2的线性组合作为[xi，xi+q]上平均斜率的近似 值，当然比用欧拉法预报y(x i+q)精度要好。由 此，得到y(xi+q)的预报值。 于是可取 从而得到三点计算公式的形式 （8-10） 下面列出其中的一种形式 类似前面的推导，利用泰勒展开适当选 择参数λ1，λ2和p，q，r，s的值使上述公式 具有三阶精度，即局部截断误差为O(h4)。这 类公式统称为三阶龙格-库塔公式。 （8-11） 高阶龙格-库塔法 为了进一步提高精度，在[xi，xi+1]上可取多个点，预报相应点的斜率值，对这些斜率值加权平均作为平均斜率值。利用泰勒展开，比较相应系数，从而确定在尽可能高的精度下有关参数应满足的条件。 从理论上讲，可以构造任意高阶的龙格 -库塔公式，但实践证明，高于四阶的龙格-库塔公式，不但计算量大，而且精确度并不一定提高。在实际计算中，四阶龙格-库塔 公式是精度及计算量较理想的公式。 经典的四阶龙格-库塔公式为 （8-12） 龙格-库塔法的步长选择 上面所讨论的龙格-库塔方法是定 步长的，单从每一步来说，步长 h 越 小，局部截断误差就越小；但随着步 长的缩小，不但引起计算量的增加， 而且也有可能引起舍入误差的严重积 累；而步长h太大又不能达到预期的精度要求，所以怎样选择合适的步长h，这在实际计算中是很重要的。 我们以四阶龙格-库塔公式为例，从 节点xi出发，先以某个h为步长求出一个 近似值，记为yi+1(h)，然后将步长折半， 即取h/2为步长从xi计算两步到xi+1，再求 得一个近似值yi+1(h/2) 。 通过计算| yi+1(h/2) -yi+1(h) |ε是否成 立，成立则步长选择h/2，否则继续将区 间分成h/4，再次重复上述步骤至满足精 度要求。 例 使用高阶R-K方法计算初值问题 解： (1) 使用三阶R-K方法 其余结果如下: (2) 如果使用四阶R-K方法 i xi k1 k2 k3 yi 1.0000 0.1000 1.0000 1.1025 1.2555 1.1111 2.0000 0.2000 1.2345 1.3755 1.5945 1.2499 3.0000 0.3000 1.5624 1.7637 2.0922 1.4284 4.0000 0.4000 2.0404 2.3423 2.8658 1.6664 5.0000 0.5000 2.7768 3.2587 4.1634