线性多步法详解演示文稿

线性多步法详解演示文稿当前第1页\共有20页\编于星期二\13点（优选）线性多步法当前第2页\共有20页\编于星期二\13点.

一、Adams外插法线性多步法的实质是用若干节点处的函数值及导数值的线性组合来逼近y(xi+1)的值．由于在计算yi+1时需要用到其前k+1个值：yi,yi−1,yi−2,,yi−k，故多步法又称为k+1步法当−1=0时，此k+1步法是显式的当−10时，此k+1步法是隐式的原则上，一切多步法都可用Taylor级数展开的方法来导出(即确定其中的系数j,j)，但有些多步法也可用数值积分法来构造．下面讨论实际中较为常用的几种线性多步法当前第3页\共有20页\编于星期二\13点(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外插法当前第4页\共有20页\编于星期二\13点外插公式略去积分余项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)式为外插公式当前第5页\共有20页\编于星期二\13点考虑到插值点x靠近区间[xi−k,xi]的最后一个节点xi，采用Newton向后插值公式，于是有其中，而▽为向后差分算子,即所以称为Adams外插公式(亦有称Adams-Bashforth公式)．不难看出，这是一类k+1步显式方法.当前第6页\共有20页\编于星期二\13点容易计算出bjj012345bj11/25/123/825/72095/288(2)Adams外插公式的局部截断误差注意到，若假定yi−j=y(xi−j),j=0,1,,k精确成立，则有．因此，Adams外插公式的局部截断误差为其中当前第7页\共有20页\编于星期二\13点Adams外插公式是一类k+1步k+1阶的显式方法.几个常用的Adams外插公式如下①单步法(k=0)②二步法(k=1)当前第8页\共有20页\编于星期二\13点③三步法(k=2)④四步法(k=3)当前第9页\共有20页\编于星期二\13点现在以k+2个节点x−k,xi−k+1,,xi,xi+1作为插值节点,作函数f(x,y(x))的k+1次插值多项式pk+1(x)，从而有其中R(x)为插值余项去掉上式中的积分余项，得二、Adams内插法(1)公式的推导当前第10页\共有20页\编于星期二\13点同样，采用Newton向后插值公式，并重复Adams外插公式的推导过程，可得其中dj的前几个值j012345dj1-1/2-1/12-1/24-19/720-3/160由于插值点x

现在是落在积分区间[xi

,xi+1]之内，所以称为Adams内插公式(Adams-Moutton公式)．这是一类k+1步隐式方法当前第11页\共有20页\编于星期二\13点(2)局部截断误差其中而i是介于xi−k与xi+1之间的某个值.所以,Adams内插公式是一类k+1步k+2阶的隐式方法.当前第12页\共有20页\编于星期二\13点几个常用的Adams内插值公式①单步法(k=0)②二步法(k=1)③三步法(k=2)当前第13页\共有20页\编于星期二\13点优点步数相同的Adams内插公式比外插公式在精度上要高一阶，而阶数相同的内插公式的截断误差也比外插值公式的截断误差小许多缺点内插法是隐式的，求解用迭代法，因而计算量较大对初值问题分别用四步四阶Adams外插法和三步四阶Adams内插法求解取步长h=0.1，即N=10例1解当前第14页\共有20页\编于星期二\13点将f(x,y)=2x+y,h=0.1,xi

=0.1i

代入，得四步四阶Adams外插法的公式为当前第15页\共有20页\编于星期二\13点三步四阶Adams内插法的公式为将f(x,y)=2x+y,h=0.1,xi

=0.1i

代入，得本例可以解出yi+1

使其成为显式当前第16页\共有20页\编于星期二\13点本例的精确解为y(x)=3ex−2x−2.利用此精确解求出y1=y(x1),y2=y(x2),y3=y(x3)，并用y0,y1,y2,y3

作为上述外插公式的