第九章-常微分方程数值解法2.pptVIP

* 一、Runge-Kutta方法的基本思想 §3 龙格-库塔（Runge-Kutta）方法 显式单步法的一般形式： R-K方法是利用一些点的线性组合构造增量函数， 使得相应方法的局部截断误差的阶数尽可能高。 ?二阶Runge-Kutta方法 确定参数 ，使得 与 在点 的Taylor展开式有尽可能多的相同项。 比较两式的相同项得 方程组有无穷多解 若取其一组解 则得到改进的Euler公式（二阶方法） 若取其另一组解 则得到二阶的Heun（休恩）公式（见教材）。 二、显式Runge-Kutta方法及其稳定性 和 设 是一个正整数，代表使用函数值 的个数， 是一些特定的权因子（均为 实数），则称下列方法（公式） 为初值问题 的m级显式Runge–Kutta公式， 其中 类似前面的处理方法，可以得到四级方法：m =4 局部截断误差 最常用的一种四阶方法：经典显式Runge-Kutta公式 解： 例2：用经典的四阶Runge-Kutta 方法求解下列初值问题 。 经典的四阶Runge-Kutta公式： 注：?对于显式N级R-K方法，最多只能得到N阶方法。 ?上述方法的缺陷：计算非常复杂。 ?可通过积分方法确定参数。 例2：确定如下三级三阶显式Runge-Kutta公式中的参数： 解： 对微分方程 两边积分得 采用Simpson公式计算上式右端积分项 可设参数 则有 选择剩余参数，使得 取 取 利用Taylor展开式 代入 当 时， 例3：求经典四阶的R-K方法的绝对稳定域。 解： 其绝对稳定域为 三、隐式Runge-Kutta方法 m级隐式R–K方法的一般形式 其中系数的确定方法同显式R–K方法完全类似 (1)一级二阶的隐式中点方法： (2)二级四阶的隐式R-K方法： N级隐式R-K法 可以达到2N阶 缺陷：需要求解 非线性方程（组） 一、k步线性多步法 §4 线性多步法与预估-校正格式 /*Linear Mutistep Method and Predictor-Corrector Format*/ 所谓的线性多步法，指的是某一步解的公式不仅与前一步的值有关，而且与前面若干步解的值有关的方法。 对初值问题 两边积分得 将 换为节点 取节点 ，构造 的q+1个点的 Lagrange插值多项式： ?多步显式公式 其中 记 若函数值 已知，则得r+1步显式方法 如 时，可得二步显式阿达姆斯（Adams）格式 其中 ? Adams显式公式的局部截断误差： 由Lagrange插值余项知 其中 （第二积分中值定理） q+1阶方法 取节点 ，构造 的q+1个点的 Lagrange插值多项式： ?多步隐式公式 其中 记 则得到r+1步q+1阶的隐式方法 如 时，可得二步隐式阿达姆斯（Adams）格式 梯形公式 ?常用的一种预报-校正公式： 四阶Adams预报-校正公式： （显式） （隐式） 初始迭代值由4阶R-K方法计算 例4：用Adams预报-校正公式 求解下列初值问题 。 解： Adams预报-校正公式： 1.7320507198 1.6733199993 1.6124515364 1.5491933804 1.4832398242 1.4142138334 1.3416413571 Adams预-校法 1.7320508075 1.0 1.6733200530 0.9 1.6124515496 0.8 1.5491933384 0.7 1.4832396974 0.6 1.4142135623 0.5 1.3416407864 0.4 1.2649110640 1.264912 0.3 1.1832159566 1.183217 0.2 1.0954451153 1.095446 0.1 1.0000000000 1 0 精确解 R-K方法 一、单步法的收敛性 /*Convergence of Onestep Met