9.3 龙格-库塔方法.ppt

9.3 龙格-库塔方法 9.3.1 显式龙格-库塔法的一般形式 9.3.2 二阶显式R-K方法 9.3.3 三阶与四阶显式R-K方法 9.3.4 变步长的龙格-库塔方法 上节给出了显式单步法的表达式 其局部截断误差为 对欧拉法 ，即方法为 阶. （3.1） 若用改进欧拉法，它可表示为 此时增量函数 （3.2） 与欧拉法的 相比，增加了计算一个 右函数 的值，可望 . 若要使得到的公式阶数 更大， 就必须包含更多的 值. （3.3） 从方程 等价的积分形式（2.4），即 若要使公式阶数提高，就必须使右端积分的数值求积公式 精度提高，必然要增加求积节点. 为此可将（3.3）的右端用求积公式表示为 点数 越多，精度越高， 上式右端相当于增量函数 ， 为得到便于计算的显式方法，可类似于改进欧拉法，将公式表示为 （3.4） 其中 （3.5） 这里 均为常数. (3.4)和(3.5)称为 级显式龙格-库塔（Runge-Kutta）法， 简称R-K方法. 当 时，就是欧拉法， 此时方法的阶为 . 当 时，改进欧拉法(3.1)，(3.2)也是其中的一种. 下面将证明阶 . 要使公式(3.4),(3.5)具有更高的阶 ,就要增加点数 . 下面就 推导R-K方法. （3.4） （3.5） 对 的R-K方法，计算公式如下 （3.6） 这里 均为待定常数. 希望适当选取这些系数，使公式阶数 尽量高. 根据局部截断误差的定义，（3.6）的局部截断误差为 （3.7） 这里 . 为得到 的阶 ，要将上式各项在 处做泰 勒展开， 由于 是二元函数，故要用到二元泰勒展开， 其中 各项展开式为 （3.8） 将以上结果代入局部截断误差公式则有 要使公式（3.6）具有 阶，必须使 即 （3.9）的解是不惟一的. 令 ，则得 这样得到的公式称为二阶R-K方法， 如取 ，则 这就是改进欧拉法（3.1）. （3.9） 若取 则 . 称为中点公式， （3.10）也可表示为 得计算公式 （3.10） 相当于数值积分的中矩形公式. 的R-K公式(3.6)的局部误差不可能提高到 . 把 多展开一项，从（3.8）的 看到展开式中 的项是不能通过选择参数消掉的. 实际上要使 的项为零，需增加3个方程，要确定4个 参数 ，这是不可能的. 故 的显式R-K方法的阶只能是 ，而不能得 到三阶公式. 要得到三阶显式R-K方法，必须 . （3.11） 其中 及 均为待定参数. 此时（3.4）,（3.5）的公式表示为 公式（3.11）的局部截断误差为 只要将 按二元函数泰勒展开，使 ， 可得待定参数满足方程 （3.12） 这是8个未知数6个方程的方程组，解也不是惟一的. 所以 这是一簇公式. 满足条件（3.12）的公式（3.11）统称为三阶R-K公式. 一个常见的公式为 此公式称为库塔三阶方法. 继续上述过程，经过较复杂的数学演算，可以导出各 种四阶龙格-库塔公式，下列经典公式是其中常用的一个： 可以证明其截断误差为 . 四阶龙格-库塔方法的每一步需要计算4次函数值 ， （3.13） 例3 解 设取步长 ，从 到 用四阶龙格- 库塔方法求解初值问题 这里 ，经典的四阶龙格-库塔公式为 表9-3列出了计算结果，同时列出了相应的精确解. 比较例3和例2的计算结果，显然龙格-库塔方法的精度高. 但由于这里放大了步长 ，所以表9-3和表9-2 所耗费的计算量几乎相同. 龙格-库塔