第九章 常微分方程初值问题数值解法-数值的分析.ppt

* 预估校正法 对于隐式线性多步法, 计算时要进行迭代, 计算量大, 通常 采用显式方法先给出yn+k 的初始近似, 记为y p n+k 称为P(Predictor) 预估, 再计算 fn+k 的函数值E(Evaluation) , 继而以隐式公式计 算yn+k 称为C(Corrector)校正及M(Modifier)改进 阿达姆斯四阶预测-校正公式的PECE模式 * 阿达姆斯公式的PMECME模式 由四阶阿达姆斯公式的截断误差,对于PECE有 二式相减, * 由此得事后误差估计： 于是加上两个改进步(M),最后再求一次函数值(E)供下一次预测用,就形成阿达姆斯法的PMECME模式 * 修正米尔尼-汉明预测-校正法(PMECME) 利用米尔尼公式(5.11)进行预测,再用哈明公式(5.15)进行校正,中间利用两公式局部截断误差的主项系数的比例进行改进两次,就得到PMECME模式,此模式称修正米尔尼-汉明预测-校正法 有时侯并不一定要按系数公式(5.4)确定多步法(5.1)的 系数 αi , βi ,可直接用泰勒展开进行构造. 例 解初值问题的显式二步法 试确定系数αi , βi使方法阶数尽量的高,并求局部截断误差 解: * * 为使方法阶数尽量高,令 解得: 此时为三阶公式, 且所求局部误差为: 例 证明