[数学]数值方法.pptVIP

7.5 亚当姆斯方法 7.5.1 亚当姆斯格式 龙格-库塔方法是一类重要算法，但这类算法在每一步都需要先预报几个点上的斜率值，计算量比较大。考虑到计算yi+1之前已得出一系列节点上 的斜率值，能否利用这些已知值来减少计算量呢？ 这就是亚当姆斯（Adams）方法的设计思想。 设用xi,xi+1两点的斜率值加权平均作为区间 上的平均斜率，有计算格式 （7.21） 选取参数λ，使格式（7.21）具有二阶精度。 将 在xi处Taylor展开 代入计算格式(7.21)化简,并假设, 因此有 与y(xi+1)在xi处的Taylor展开式 相比较, 需取 才使格式(7.21)具有二阶精度。这样导出的计算格式 称之为二阶亚当姆斯格式。类似地可以导出三阶亚当姆斯格式。 和四阶亚当姆斯格式。 (7.22) 这里和下面均记 上述几种亚当姆斯格式都是显式的，算法比较简单，但用节点 的斜率值来预报区间 上的平均斜率是个外推过程，效果不够理想。为了进一步改善精度，变外推为内插，即增加节点xi+1的斜率值来得出 上的平均斜率。譬如考察形如 （7.23） 的隐式格式,设(7.23)右端的 Taylor展开有 可见要使格式(7.23)具有二阶精度,需令 , 这样就可构造二阶隐式亚当姆斯格式 其实是梯形格式。类似可导出三阶隐式亚当姆斯格式 和四阶隐式亚当姆斯格式 （7.24） 7.5.2 亚当姆斯预报-校正格式 参照改进的欧拉格式的构造方法，以四阶亚当姆斯为例，将显式（7.22）和隐式（7.24）相结合，用显式公式做预报，再用隐式公式做校正，可构成亚当姆斯预报-校正格式 （7.25） 预报： 校正： 这种预报-校正格式是四步法，它在计算yi+1时不但用到前一步的信息 ，而且要用到再前面三步的信息 ，因此它不能自行启动。在实际计算时，可借助于某种单步法，譬如四阶龙格—库塔法提供开始值 。 例7.5 取步长h=0.1,用亚当姆斯预报-校正公式求解 初值问题 的数值解。 解: 用四阶龙格-库塔公式求出发值 ,计算得： 表中的 ,yi和y(xi)分别为预报值、校正值和准确解( ),以比较计算结果的精度。 再使用亚当姆斯预报-校正公式(7.25),见教材P166列表算得其余的计算结果 7.6 算法的稳定性及收敛性 7.6.1 稳定性 稳定性在微分方程的数值解法中是一个非常重要的问题。因为微分方程初值问题的数值方法是用差分格式进行计算的，而在差分方程的求解过程中，存在着各种计算误差，这些计算误差如舍入误差等引起的扰动，在传播过程中，可能会大量积累，对计算结果的准确性将产生影响。这就涉及到算法稳定性问题。 当在某节点上xi的yi值有大小为δ的扰动时，如果在其后的各节点 上的值yi产生的偏差都不大于δ，则称这种方法是稳定的。 稳定性不仅与算法有关，而且与方程中函数f(x,y)也有关，讨论起来比较复杂。为简单起见，通常只针对模型方程 来讨论。一般方程若局部线性化，也可化为上述形式。模型方程相对比较简单，若一个数值方法对模型方程是稳定的，并不能保证该方法对任何方程都稳定，但若某方法对模型方程都不稳定，也就很难用于其他方程的求解。 先考察显式Euler方法的稳定性。模型方程 的Euler公式为 将上式反复递推后，可得 或 式中 要使yi有界，其充要条件为 即 由于 ，故有 （7.26） 可见，如欲保证算法的稳定，显式Euler格式的步长h的选取要受到式（7.26）的限制。 的绝对值越大，则限制的h值就越小。 用隐式Euler格式，对模型方程 的计算公式为，可化为 由于 ,则恒有 ,故恒有 因此，隐式Euler格式是绝对稳定的（无条件稳定）的（对任何h0）。 7.6.2 收敛性 常微分方程初值问题的求解，是将微分方程转化为差分方程来求解，并用计算值yi来近似替代y(xi)，这种近似替代是否合理，还须看分割区间 的长度h越来越小时，即 时， 是否成立。若成立，则称该方法是收敛的，否则称为不收敛