常微分方程数值解.pptVIP

*第94页，共151页，星期日，2025年，2月5日例6用四阶阿当姆斯显式和隐式方法解初值问题取步长.解本题.从四阶阿当姆斯隐式公式得到从四阶阿当姆斯显式公式得到*第95页，共151页，星期日，2025年，2月5日由此可直接解出而不用迭代，得到计算结果见表9-8，其中显式方法中的初值及隐式方法中的均用准确解得到，对一般方程，可用四阶R-K方法计算初始近似.(表略，见教材.)*第96页，共151页，星期日，2025年，2月5日从以上例子看到同阶的阿当姆斯方法，隐式方法要比显式方法误差小.这可以从两种方法的局部截断误差主项的系数大小得到解释.这里分别为及.*第97页，共151页，星期日，2025年，2月5日9.5.3米尔尼方法与辛普森方法考虑另一个的显式公式其中为待定常数.由（5.4）可知，再令得到可根据使公式的阶尽可能高这一条件来确定其数值.*第98页，共151页，星期日，2025年，2月5日解此方程组得于是得到四步显式公式（5.11）称为米尔尼(Milne）方法.由于，故方法为4阶，其局部截断误差为（5.12）米尔尼方法也可以通过方程两端积分*第99页，共151页，星期日，2025年，2月5日得到.右端积分通过辛普森求积公式就有（5.13）称为辛普森方法.（5.14）它是隐式二步四阶方法，其局部截断误差为若将方程从到积分，可得*第100页，共151页，星期日，2025年，2月5日9.5.4汉明方法辛普森公式是二步方法中阶数最高的，但它的稳定性差，为了改善稳定性，考察另一类三步法公式其中系数及为常数.若希望导出的公式是四阶的，则系数中至少有一个自由参数.若取，则可得到辛普森公式.若取，仍利用泰勒展开，由（5.4），令*第101页，共151页，星期日，2025年，2月5日解此方程组得于是有（5.15）则可得到*第102页，共151页，星期日，2025年，2月5日称为汉明(Hamming)方法.由于，故方法是四阶的，且局部截断误差（5.16）*第103页，共151页，星期日，2025年，2月5日9.5.5预测-校正方法对于隐式的线性多步法，计算时要进行迭代，计算量较大.为了避免进行迭代，通常采用显式公式给出的一个初始近似，记为，称为预测(predictor).接着计算的值(evaluation)，再用隐式公式计算，称为校正（corrector）.在（2.13）中用欧拉法做预测，再用梯形法校正，得到改进欧拉法，它就是一个二阶预测-校正方法.（2.13）*第104页，共151页，星期日，2025年，2月5日例如用四阶的阿当姆斯显式方法做预测，再用四阶阿当姆斯隐式公式做校正，得到以下格式：一般情况下，预测公式与校正公式都取同阶的显式方法与隐式方法相匹配.预测P：求值E：校正C：求值E：此公式称为阿当姆斯四阶预测-校正格式（PECE).*第105页，共151页，星期日，2025年，2月5日依据四阶阿当姆斯公式的截断误差，对于PECE的预测步P有对校正步C有两式相减得于是有下列事后误差估计*第106页，共151页，星期日，2025年，2月5日定理4阶方法(4.1)与初值问题(1.1)，(1.2)相容的充分必要条件是由定理3可知单步法（4.1）收敛的充分必要条件是(4.1)是相容的.以上讨论表明阶方法（4.1）当时与(1.1),(1.2)相容，反之相容方法至少是1阶的.*第62页，共151页，星期日，2025年，2月5日9.4.2绝对稳定性与绝对稳定域