第9讲常微分方程数值解.ppt

第9章 常微分方程初值问题数值解法 9.1 引 言 9.2 简单的数值方法 9.2.2 梯形方法 9.2.3 改进欧拉公式 9.2.4 单步法的局部截断误差与阶 9.3 龙格-库塔方法 9.3.3 三阶与四阶显式R-K方法 谢 谢！ 9.3.4 变步长的龙格-库塔方法 9.4 单步法的收敛性与稳定性 9.4.2 绝对稳定性与绝对稳定域 9.5 线性多步法 9.5.1 线性多步法的一般公式 9.5.2 阿当姆斯显式与隐式公式 9.5.3 米尔尼方法与辛普森方法 9.5.4 汉明方法 9.5.5 预测-校正方法 9.5.6 构造多步法公式的注记和例 9.7 一阶方程组与刚性方程 9.7.2 化高阶方程为一阶方程组 9.7.3 刚性方程组 由此可直接解出 而不用迭代，得到 计算结果见表9-8，其中显式方法中的初值 及 隐式方法中的 均用准确解 得到， 对一般方程，可用四阶R-K方法计算初始近似. (表略，见教材. ) 从以上例子看到同阶的阿当姆斯方法，隐式方法要比 显式方法误差小. 这可以从两种方法的局部截断误差主项 的系数大小得到解释. 这里 分别为 及 . 考虑另一个 的显式公式 其中 为待定常数. 由（5.4）可知 ，再令 得到 可根据使公式的阶尽可能高这一条件来确定其数值. 解此方程组得 于是得到四步显式公式 （5.11） 称为米尔尼(Milne）方法. 由于 ，故方法为4阶，其局部截断误差为 （5.12） 米尔尼方法也可以通过方程 两端积分 得到. 右端积分通过辛普森求积公式就有 （5.13） 称为辛普森方法. （5.14） 它是隐式二步四阶方法，其局部截断误差为 若将方程 从 到 积分，可得 辛普森公式是二步方法中阶数最高的，但它的稳定性差，为了改善稳定性，考察另一类三步法公式 其中系数 及 为常数. 若希望导出的公式是四阶的，则系数中至少有一个自由参数. 若取 ，则可得到辛普森公式. 若取 ，仍利用泰勒展开，由（5.4），令 解此方程组得 于是有 （5.15） 则可得到 称为汉明(Hamming)方法. 由于 ，故方法是四阶的，且局部截断误差 （5.16） 对于隐式的线性多步法，计算时要进行迭代，计算量 较大. 为了避免进行迭代，通常采用显式公式给出 的一个初始近似，记为 ，称为预测(predictor). 接着计算 的值(evaluation)，再用隐式公式计算 ，称为校正（corrector）. 在（2.13）中用欧拉法做预测，再用梯形法校正，得 到改进欧拉法，它就是一个二阶预测-校正方法. （2.13） 例如用四阶的阿当姆斯显式方法做预测，再用四阶阿 当姆斯隐式公式做校正，得到以下格式： 一般情况下，预测公式与校正公式都取同阶的显式方 法与隐式方法相匹配. 预测P： 求值E： 校正C： 求值E： 此公式称为阿当姆斯四阶预测-校正格式（PECE). 依据四阶阿当姆斯公式的截断误差，对于PECE的预测步P有 对校正步C有 两式相减得 于是有下列事后误差估计 容易看出 比 更好. （5.17） 但在 的表达式中 是未知的，因此计算时用 上一步代替，从而构造一种修正预测-校正格式(PMECME): P： M： E： C： M： E： 在PMECME格式中已将（5.17）的 及 分别改为 及 . 利用米尔尼公式（5.11）和汉明公式（5.15）相匹配， 并利用截断误差（5.12），（5.16）改进计算结果，可类 似地建立四阶修正米尔尼-汉明预测-校正格式PMECME）： P： M： E：