单步法的收敛性和稳定性-8.5线性多步法及其预测-校正格式1.pptVIP

例分别用改进的欧拉格式和四阶龙格—库塔格式解初值问题（取步长h=0.2）节点改进欧拉法四阶龙格—库塔法准确解01110.21.1866671.1832291.1832160.41.3483121.3416671.3416410.61.4937041.4832811.4832400.81.6278611.6125141.61245211.7542051.7321421.732051（注：已指出过准确解)x四阶二阶真解四阶误差二阶误差0.01.0000001.0000001.0000000.00000.0000000.11.1048291.1024501.1048291.60E-72.38E-30.21.2185971.2115071.2185973.40E-77.09E-30.31.3401411.3257661.3401415.48E-71.44E-20.41.4681751.4436711.4681757.69E-72.45E-20.51.6012781.5635061.6012799.95E-73.78E-20.61.7378801.6833741.7378811.20E-65.45E-20.71.8762461.8011791.8762471.42E-67.51E-20.82.0144571.9146032.0144591.68E-69.99E-20.92.1503952.0210862.1503971.96E-61.29E-11.02.2817162.1178002.2817182.32E-61.64E-1例关于Runge-Kutta方法1.4.4误差控制和法对任意给定，可以求出最大的，使得格式(1.41)的整体截断误差小于。(1.41)逼近初值问题(1.1)，(1.2)解的理想单步格式：应该具有性质：其局部截断误差是，其中下面我们以欧拉格式为例，给出一种局部误差估计的技术。而梯形法的局部截断误差是。如果所以那么(1.42)因此可以作为欧拉法的局部误差的近似值。假设有二个逼近初值问题(1.1)，(1.2)的单步格式，其中一个格式的局部截断误差是，而另一个差分格式的局部截断误差是。下面我们给出如何利用格式局部截断误差的估计值来确定用于控制整体截断误差的最佳步长。类似于对欧拉法的分析，我们有(1.43)(1.44)这里是，所以存在常数，使得由式(1.43)和(1.44)得(1.45)用代替截断误差中的，其中是一个大于零的有界常数。用表示这一截断误差。由式(1.43)和(1.44)可得(1.46)为了使，我们选择常数满足即(1.47)使用不等式(1.47)进行误差控制的一种常用的方法是由在1970年提出的法。Runge-Kutta-Fehlberg方法Fehlberg设计了一个精巧的嵌套方法如下：Runge-Kutta-Fehlberg方法Fehlberg给出的四阶、五阶公式RKF4(5)如下：Runge-Kutta-Fehlberg方法七阶、八阶RKF7(8)Runge-Kutta-Fehlberg方法七阶、八阶RKF7(8)线性多步法常微分方程初值问题的数值解法中，除了Runge-Kutta型公式等单步法之外，还有另一种类型的解法，即某一步的公式不仅与前一步解的值有关，而且与前若干步解的值有关，利用前面多步的信息预测下一步的值，这就是多步法的基本思想，可以期望获得较高的精度。构造多步法有多种途径，下面先讨论基于数值积分的方法。1.5线性多步法前面利用欧拉法或梯形法求未知函数在的近似值，基本思想是在积分表达式中被积函数，用水平直线或连接，两点的直线代替。然而为了近似中的