9.5 线性多步法.ppt

9.5 线性多步法 9.5.1 线性多步法的一般公式 9.5.2 阿当姆斯显式与隐式公式 9.5.3 米尔尼方法与辛普森方法 9.5.4 汉明方法 9.5.5 预测-校正方法 9.5.6 构造多步法公式的注记和例 利用米尔尼公式（5.11）和汉明公式（5.15）相匹配， 并利用截断误差（5.12），（5.16）改进计算结果，可类 似地建立四阶修正米尔尼-汉明预测-校正格式PMECME）： P： M： E： C： M： E： 将例6的初值问题用修正的米尔尼-汉明预测-校 正公式计算 及 ，初值 仍用已算出的精 确解， 例7 即 给出计算结果及误差. 解 其中 根据修正的米尔尼-汉明预测-校正公式可得 误差 从结果看，此方法的误差比四阶阿当姆斯隐式法和四阶汉 明方法小，这与理论分析一致. 构造多步法公式有基于数值积分和泰勒展开两种途径， 而用泰勒展开则可构造任意多步法公式，其做法是根 据多步法公式的形式，直接在 处做泰勒展开. 确定多步法（5.1）的系数 及 时不 必套用系数公式（5.4），因为多步法公式不一定如（5.1） 的形式， 数值积分方法只能用于可转化为等价的积分方程的微分方程， 它是有局限性的. 而且套用公式容易记错. 在逐步推进的求解过程中，计算 之前事实上已经 求出了一系列的近似值 如果充分利用前面多步的信息来预测 ，则可以期望获得较高的精度. 这就是构造所谓线性多步法的基本思想. 本节主要介绍基于泰勒展开的构造方法. 构造多步法的主要途径是基于数值积分方法和基于泰 勒展开方法，前者可直接由方程 两端积分后 利用插值求积公式得到. 一般的线性多步法公式可表示为 （5.1） 其中 为 的近似， 计算时需先给出前面 个近似值 ,再由 （5.1）逐次求出 . 如果计算 时，除了使用 的值，还用到 的值，则称此方法为线性多步法. 为常数. 及 不全为零， 则称为线性 步法. 若 如果 ，称（5.1）为显式 步法，这时 可 直接由（5.1）算出. 如果 ，则（5.1）称为隐式 步法. 求解时与梯形法（2.7）相同，要用迭代法方可算出 . （2.7） 设 是初值问题(1.1)，(1.2)的准确解，线性多步法（5.1）在 上的局部截断误差为 定义8 （5.2） （5.1）中系数 及 可根据方法的局部截断误差及阶确定，其定义如下. 若 ，则称方法（5.1）是 阶的， 则称方法（5.1）与方程（1.1）是相容的. （5.1） 由定义8，对 在 处做泰勒展开，由于 代入（5.2）, （5.3） （5.2） 得 其中 （5.4） 若在公式（5.1）中选择系数 及 ，使它满足 由定义可知,此时所构造的多步法是 阶的. （5.5） 称右端第一项为局部截断误差主项， 称为误差常数. 根据相容性定义， , 即 故方法（5.1）与微分方程（1.1）相容的充分必要条件是 （5.6）成立. 且 （5.6） 由（5.4）得 当 时，若 ，则由（5.6）可求得 此时公式（5.1）为 即为欧拉法. 从（5.4）可求得 ，故方法为一阶精度， 这和第2节给出的定义及结果是一致的. 且局部截断误差为 对 , 若 ，方法为隐式公式. 为了确定系数 ，可由 解得 于是得到公式 即为梯形法. 由（5.4）可求得 , 故 ，所以梯形法 是3阶方法， 其局部截断误差主项是 . 这与第2节中的讨论也是一致的. 对 的多步法公式都可利用（5.4）确定系数 ， 并由（5.5）给出局部截断误差. （5.5） 考虑形如 （5.7） 的 步法，称为阿当姆斯（Adams)方法. 为显式方法， 为隐式方