Chapter5_3_线性多步法.pptVIP

§5 线性多步法/*Linear multistep method*/ . 一、Adams外插法 * 一、Adams外插法 二、Adams内插法 三、Taylor级数法 求解初值问题的数值方法都是“步进式”的，即求 解过程从初值y0开始，顺着节点的排列次序，一 步一步地向前推进．所以，在计算yi+1 时，前面 的i + 1 个值y0, y1, ?, y i 都是已知的．如果在计算 yi+1 时能充分利用这些已有的信息，而不是像单 步法中那样，只用其前一步的值yi，则可望构造 出精度高，但计算量小的求解公式．线性多步法 就是基于这一思想发展起来的，其计算公式可表 示为 其中 线性多步法的实质是用若干节点处的函数值及导 数值的线性组合来逼近y(xi+1)的值． 由于在计算yi+1时需要用到其前k + 1个值：yi , yi?1, yi?2,?, y i ?k，故多步法又称为k + 1步法 当??1= 0时，此k + 1步法是显式的 当??1? 0时，此k + 1步法是隐式的 原则上，一切多步法都可用Taylor级数展开的方 法来导出(即确定其中的系数?j ,?j )，但有些多步 法也可用数值积分法来构造． 下面讨论实际中较为常用的几种线性多步法 (1)公式的推导 对方程y’ = f (x, y)的两边从xi到xi+1积分，得 为了近似计算式中的积分，以xi?k , xi?k+1, ?, xi?1, xi 为插值节点，作函数f (x, y (x)) 的k 次插值多项 式pk (x)，从而有 f (x, y (x) ) = pk (x) + R (x), 其中，R (x)为插值余项 一、Adams外插法 外插公式 略去积分余项R(x)dx，并用yi代替y (xi )，可得到 计算公式 注意，这里的 与pk(x)不同，它是将插值多 项式pk(x)中用到的函数值f(xi,y(xi))均以近似值 f i = f(xi,yi)代替后所得到的表达式. 由于积分是在区间[xi , xi+1]上进行的，而 的 插值区间却为[xi?k,xi]，即插值点x ? [xi , xi+1]， 位于插值区间之外，故称(53)式为外插公式 考虑到插值点x靠近区间[xi?k,xi]的最后一个节点 xi，采用Newton向后插值公式，于是有 其中 ， 而▽为向后差分算子,即 所以 称为Adams外插公式(亦有称Adams-Bashforth公式)． 不难看出，这是一类k + 1步显式方法. 容易计算出bj 95/288 25/720 3/8 5/12 1/2 1 bj 5 4 3 2 1 0 j (2)Adams外插公式的局部截断误差 注意到，若假定yi?j=y(xi?j) , j= 0,1,?,k精确成立， 则有 ．因此，Adams外插公式的 局部截断误差为 其中 Adams外插公式是一类k + 1步k +1阶的显式方法. 几个常用的Adams外插公式如下 ① 单步法(k=0) ② 二步法(k=1) ③ 三步法(k=2) ④ 四步法(k=3) 现在以k +2个节点x?k , xi?k+1, ? , xi , xi+1作为插值 节点,作函数f(x,y(x))的k+1次插值多项式pk+1(x)， 从而有 其中R (x)为插值余项 去掉上式中的积分余项，得 二、Adams内插法 (1)公式的推导 同样，采用Newton向后插值公式，并重复Adams 外插公式的推导过程，可得 其中 d j的前几个值 -3/160 -19/720 -1/24 -1/12 -1/2 1 dj 5 4 3 2 1 0 j 由于插值点x 现在是落在积分区间[xi , xi+1]之内， 所以称为Adams内插公式(Adams-Moutton公式)． 这是一类k + 1步隐式方法 (2)局部截断误差 其中 而?i是介于xi?k与xi+1之间的某个值. 所以, Adams 内插公式是一类k + 1步k + 2阶的隐式方法. 几个常用的Adams内插值公式 ① 单步法(k=0) ② 二步法(k=1) ③ 三步法(k=2) 优点 步数相同的Adams内插公式比外插公式在精 度上要高一阶，而阶数相同的内插公式的截 断误差也比外插值公式的截断误差小许多 缺点 内插法是隐式的，求解用迭代法，因而计算 量较大 对初值问题 分别用四步四阶Adams外插法和三步四阶Adams 内插法求解 取步长h = 0.1，即N = 10 例1 解 将 f (x, y) = 2x + y, h = 0.1, xi = 0.1i 代入，得 四步四阶Adams外