Chapter单步法的收敛性和稳定性.ppt

③ 三步法(k=2) ④ 四步法(k=3) 二、Adams内插法 现在以k +2个节点x ?k , xi?k+1, ? , xi , xi+1作为插值节点，作函数f (x, y (x) )的k + 1次插值多项式pk+1(x)，从而有 其中R (x)为插值余项．去掉上式中的积分余项，得 同样，采用Newton向后插值公式，并重复Adams外插公式的推导过程，可得 其中 表8-8列出了d j的前几个值. 由于插值点x 现在是落在积分区间[xi , xi+1]之内，所以称公式(60)为Adams内插公式(亦称Adams-Moutton公式)．容易看出，这是一类k + 1步隐式方法. 由(59)式可知，公式(60)的局部截断误差为 其中 而?i是介于xi ?k与xi+1之间的某个值．所以，Adams内插公式(60)是一类k + 1步k + 2阶的隐式方法. 几个常用的Adams内插值公式如下： ① 单步法(k=0) ② 二步法(k=1) ③ 三步法(k=2) 从上面所列的公式可以看到，步数相同的Adams内插公式比外插公式在精度上要高一阶，而阶数相同的内插公式的截断误差也比外插值公式的截断误差小许多，这是内插法的优点．但内插法是隐式的，求解用迭代法，因而计算量较大，这是它的缺点. 例8 对初值问题 分别用四步四阶Adams外插法和三步四阶Adams内插法求解. 解取步长h = 0.1，即N = 10. 四步四阶Adams外插法的公式为 将 f (x, y) = 2x + y，h = 0.1， xi = 0.1i 代入，得 三步四阶Adams内插法的公式为 将f (x, y) = 2x + y， h = 0.1, xi = 0.1i 代入，得 本例可以解出yi+1 使其成为显式 本例的精确解为y (x) = 3ex ?2x ?2．利用此精确解求出y1 = y(x1), y2 = y(x2), y3 = y(x3)，并用y0, y1, y2, y3 作为上述外插公式的起步值，计算结果见表8-9. 从表8-9可以看到，Adams内插法比同阶的外插法精确. 在讨论收敛性之前，先介绍局部截断误差、整体截断误差的定义及其他们之间的关系 8.1.3 单步法的收敛性和稳定性 (Convergency and Stability) 一、单步法的收敛性 求解初值问题 的一般显式单步法可以写成如下形式： 1、局部截断误差 定义 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差 /* local truncation error */。 假定“yi = y(xi)”称为局部化假定 2、整体截断误差 3、局部截断误差与整体截断误差的关系 再利用引理1就可得到 结论：若单步法的局部截断误差为O（hP+1) 整体截断误差为O（hP) 条件：满足Lipschitz条件 定义 　　　 若某算法对于任意固定的 x = xi = x0 + i h，当 h?0 ( 同时 i ? ?) 时有 yi ? y( xi )，则称该算法是收敛的。 单步法的收敛性定义 结论1：收敛 整体截断误差Ei 0 结论2：只要单步法(35)式是高于零阶的方法， 判断单步法(35)式的收敛性就归结为验证其增量 函数?(x, y, h)是否满足对y的Lipschitz条件 例5 Euler方法是收敛的. 证明： 由于Euler方法是一阶方法，且其增量函数?(x, y, h) =f (x, y)．而初值问题是要求函数f (x, y)对y满足Lipschitz条件的， 故Euler方法收敛. 例6 改进Euler方法是收敛的. 证明 改进Euler方法是二阶方法，其增量函数为 下面证明，当f (x, y)满足对y的Lipschitz条件时，(42)式中的 ?(x, y, h)也满足对y的Lipschitz条件. 由(42)式有 假定h ? h0(h0为定数)，并记 ，则有 即?(x, y, h)满足对y的Lipschitz条件，故改进的Euler方法是收敛的. 二、单步法的稳定性 上面讨论单步法的收敛性，是假定(35)式的每一步计算都是准确的