计算方法-9.5线性多步法课件.ppt

《计算方法-9.5线性多步法课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2021/8/819.5线性多步法单步法计算时只用到前一步的结果，因此只要给定初值，计算就可以进行下去。但是Euler等单步法的精度都较低，龙格-库塔方法虽然可以得到较高的精度，但这类算法为了提高精度，需要增加一些非节点处的函数值的计算，在每一步都需要先预报这些非节点上的斜率值，计算量比较大。考虑到计算yi+1之前已得出一系列节点上的斜率值，能否利用这些已知值来减少计算量呢？这就是线性多步法的设计思想，可以在计算量增加不多的情况下获得较高的精度。用已知的若干节点处的y及y‘值的线性组合来近似y(

2、xn+1)。线性多步法通式可写为：2021/8/82当10时，为隐式公式；1=0则为显式公式。基于数值积分的构造法将在上积分，得到只要近似地算出右边的积分，则可通过近似y(xn+1)。而选用不同近似式Ir，可得到不同的计算公式。2021/8/83构造线性多步法的主要方法：数值积分法和泰勒展开法。对积分式分别采用矩形公式和梯形公式可得到欧拉公式和改进欧拉公式，截断误差分别为O(h2)和O(h3)。为此，我们自然可以想到，若用更高次的插值多项式来代替f(x,y)，则所得公式的精度会更高。

3、这就是基于数值积分方法构造线性多步法的起源思想。2021/8/842021/8/85若积分用节点作为积分点，则有积分系数这是显式格式，q+1阶r+1步格式。局部截断误差2021/8/86例：建立r=1，q=2的显式格式r=1，积分区间为q=2，显式格式，积分节点为所以2021/8/87同样，若以为积分节点，可以构造r+1步q+1阶隐格式2021/8/88例：建立r=2，q=2的隐格式r=2，积分区间为q=2，隐式格式，积分节点为所以2021/8/89它的截断误差较显式格式小，通常也具有更好的稳定

4、性。2021/8/810Adams方法是线性多步法的一个代表，它是利用插值多项式进行积分得出来的，这样构造线性多步法的方法称为数值求积法，它是构造线性多步法的一种途径，另外还有Taylor法。(1)显式Adams方法2021/8/811r=0，积分区间为q=1，显式格式，积分节点为从简单情况入手Adams公式－－r=0时候的多步法2021/8/812所以二阶显式Adams方法类似方法可通过增加节点得到更高精度的三阶显式Adams公式r=0，积分区间为q=2，显式格式，积分节点为所以2021/8

5、/813三阶显式Adams方法2021/8/814同理再增加节点得到四阶显式Adams公式r=0，积分区间为q=3，显式格式，积分节点为所以2021/8/815四阶显式Adams方法2021/8/816其对应的局部截断误差为2021/8/817注：一般有，其中Bq与yn+1计算公式中fn,…,fnq各项的系数均可查表得到。10123qfnfn1fn2fn3…Bq…………………常用的是q=3的4阶阿当姆斯显式公式(2)隐式Adams方法2021/8/819隐式格式表明构造定积分的近似公式中

6、包含了节点xn+1。类似显式公式的推导过程，可得到不同精度的隐式Adams公式r=0，积分区间为q=1，显式格式，积分节点为二阶隐式Adams方法r=0，积分区间为q=2，显式格式，积分节点为三阶隐式Adams方法r=0，积分区间为q=2，显式格式，积分节点为四阶隐式Adams方法2021/8/820利用q+1个节点上的被积函数值fn+1,fn,…,fnq+1构造q阶牛顿前插多项式。与显式多项式完全类似地可得到一系列隐式公式，并有，其中与fn+1,fn,…,fnq+1的系数亦可查表得到。10

7、123qfn+1fnfn1fn2…Bq…………………~常用的是q=3的4阶阿当姆斯隐式公式较同阶显式稳定2021/8/822基于Taylor展开的构造法（待定系数法）将上式中的右端各项yn1,…,ynr;fn+1,fn1,…,fnr分别在xn点作泰勒展开，与精确解y(xn+1)在xn点的泰勒展开作比较。通过令同类项系数相等，得到足以确定待定系数a0,…,ar;1,0,…,r的等式，则可构造出各阶线性多步法的公式。例：推导最高阶的二步线性多步法二步显式多步法为左端展开后相同项

8、系数相等，得到2021/8/823解得因此得阶数最高的二步显式线性多步法为其局部截断误差为是三阶方法2021/8/824对应的二步隐式线性多步法为二步隐式多步法为五个可选参数是四阶方法2021/8/825由泰勒展开推导Adams公式确定式中待定系数，使公式具有四阶精度。由泰勒展开得2021/8/826这里有7个未知量，5个方程，若令求解该方程即得四步四阶显式Adams公式：2021/8/827Adams预估校正公式由于隐式Adams公式需要用迭代法进行求解，比较麻烦，仿照欧拉预估校正公式，常把阿