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

《单步法的收敛性和稳定性讲稿.ppt》由会员分享，可在线阅读，更多相关《单步法的收敛性和稳定性讲稿.ppt（42页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、单步法的收敛性和稳定性第一页，讲稿共四十二页哦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条件

2、第六页，讲稿共四十二页哦 若若某某算算法法对于于任任意意固固定定的的 x=xi=x0+i h，当当 h0(同同时 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方法收敛.第八页

3、，讲稿共四十二页哦 例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)式的每一步计算都是准确的，即不考虑计算中的舍入误差然而这一假定是不切合实际的，用(35)式进行实际数值计算时，每一步都不可避免地含有

4、舍入误差稳定性就是讨论计算过程中的舍入误差对最终结果的影响！第十一页，讲稿共四十二页哦定义4 如果一种数值方法在节点xi的值y i有大小为i的扰动，而由这个扰动引起以后各节点上值y i(j i)的偏差j均满足|j|i|，则称该数值方法是绝对稳定的定的.第十二页，讲稿共四十二页哦 考虑一般的单步法(35)式若值y i有一个扰动i，那么用(35)式计算，得到的值yi+1就会产生一个偏差i+1若记则可将yi+1视为单步法公式的准确结果用(43)式减去(35)式，得第十三页，讲稿共四十二页哦或由此可知，单步法(35)式绝对稳定的条件是由于增量函数与微分方程的右端f 有关，从而给考察单步法的稳定性带来了

5、困难为了简化讨论，通常是用试验方程y=y (为复常数)来检验数值方法的稳定性！第十四页，讲稿共四十二页哦（1）首先考察）首先考察Euler方法的方法的稳定性定性此时增量函数(x,y,h)=f(x,y)=y，因而有 因此对于试验方程(46)，Euler方法稳定的条件是|1+|1(47)由于可以是复数，故在h的平面上，(47)式表示以点-1为中心的单位圆及其内部区域这个区域称为Euler方法的方法的绝对稳定区域定区域第十五页，讲稿共四十二页哦（2）讨论改改进Euler方法的方法的稳定性定性此时增量函数改进Euler方法的稳定性条件为第十六页，讲稿共四十二页哦（3）经典典Runge-Kutta方法的

6、方法的稳定性定性此时增量函数第十七页，讲稿共四十二页哦代入后得于是有由此得出经典Runge-Kutta方法的稳定性条件为第十八页，讲稿共四十二页哦如果仅限于讨论是实数的情形，则上述几个单步法的稳定性条件可分别简化为 Euler法稳定性条件：2 h 0，改进Euler法稳定性条件：2 h 0，经典Runge-Kutta法稳定性条件：2.785 h 0.由上面的讨论可以看到，如果方法的绝对稳定区域或区间是有限的，那么，步长h的选取要受绝对稳定性的约束.第十九页，讲稿共四十二页哦例7 对初值问题取h=0.1和0.2，用经典Runge-Kutta方法求解.解本例中=20,h分别为-2和-4前者属绝对稳

7、定区间-2.785,0，后者不属此区间问题的准确解为y=e20 x计算结果误差见表8-6第二十页，讲稿共四十二页哦隐式单步法的稳定性讨论.（1）考察向后Euler方法.对于试验方程(46)，其向后Euler法的公式为 yi+1=y i+h yi+1解出yi+1，有 ，从而得到误差(扰动)公式为由此得到绝对稳定的条件为其绝对稳定区域是以1为半径、以1为中心的圆外部，如图8-5所示.第二十一页，讲稿共四十二页哦（2）讨论梯形法的稳定性对于试验方程(46)，相应的梯形法公式为解出yi+1，有由此得出绝对稳定的条件为其绝对稳定区域为Re(h)0 的整个复平面当为实数时，其绝对稳定区间为h 0.第二十二

8、页，讲稿共四十二页哦 从以上的分析讨论可以看到，隐式方法的稳定性比显式方法好，这也是隐式方法的主要优点！第二十三页，讲稿共四十二页哦 8.1.4 线性多步法 求解初值问题的数值方法都是“步进式”的，即求解过程从初值y0开始，顺着节点的排列次序，一步一步地向前推进所以，在计算yi+1 时，前面的i+1 个值y0,y1,y i 都是已知的如果在计算yi+1 时能充分利用这些已有的信息，而不是像单步法中那样，只用其前一步的值y i，则可望构造出精度高，但计算量小的求解公式线性多步法就是基于这一思想发展起来的，其计算公式可表示为其中，而j,j 都是常数.第二十四页，讲稿共四十二页哦 线性多步法(50)

9、式的实质是用若干节点处的函数值及导数值的线性组合来逼近y(xi+1)的值由于在计算yi+1时需要用到其前k+1个值：y i,yi1,yi2,y i k，故多步法(50)式又称为k+1步法，且当1=0时，此k+1步法是显式的，而在1 0时，此k+1步法是隐式的.原则上，一切形如(50)式的多步法都可用Taylor级数展开的方法来导出(即确定其中的系数j,j)，但有些多步法也可用数值积分法来构造下面讨论实际中较为常用的几种线性多步法.一、Adams外插法第二十五页，讲稿共四十二页哦 1公式的推导 对方程y0=f(x,y)的两边从xi到xi+1积分，得为了近似计算(51)式中的积分，我们以xi k,

10、xik+1,xi1,xi 为插值节点，作函数f(x,y(x)的k 次插值多项式p k(x)，从而有f(x,y(x)=p k(x)+R(x),其中，R(x)为插值余项将上式代入(51)式，得第二十六页，讲稿共四十二页哦略去积分余项R(x)d x，并用y i代替y(xi)，可得到计算公略去积分余项R(x)d x，并用y i代替y(xi)，可得到计算公式 注意，这里的 与p k(x)不同，它是将插值多项式p k(x)中用到的函数值f(xi,y(xi)均以近似值f i=f(xi,y i)代替后所得到的表达式.由于(52)式的积分是在区间xi,xi+1上进行的，而 的插值区间却为xi k,xi，即插值点

11、x xi,xi+1，位于插值区间之外，故称(53)式为外插公式.第二十七页，讲稿共四十二页哦考虑到插值点x靠近区间xi k,xi 的最后一个节点xi，我们采用Newton向后插值公式，于是有其中 ，而为向后差分算子,即 ，等.第二十八页，讲稿共四十二页哦将(54)式代入(53)式并作变量代换，得其中(55)式称为Adams外插公式(亦有称Adams-Bashforth公式)不难看出，这是一类k+1步显式方法.第二十九页，讲稿共四十二页哦由(56)式容易计算出b j，它的前几个值见表8-7.2Adams外插公式的局部截断误差 注意到，若假定y i j=y(xi j),j=0,1,k精确成立，则有

12、 因此，Adams外插公式(55)的局部截断误差为第三十页，讲稿共四十二页哦其中而i是介于xi k与xi+1之间的某个值.由(57)式可知，Adams外插公式(55)是一类k+1步k+1阶的显式方法.几个常用的Adams外插公式如下：单步法(k=0)二步法(k=1)第三十一页，讲稿共四十二页哦 三步法(k=2)四步法(k=3)第三十二页，讲稿共四十二页哦 二、Adams内插法 现在以k+2个节点x k,xik+1,xi,xi+1作为插值节点，作函数f(x,y(x)的k+1次插值多项式pk+1(x)，从而有其中R(x)为插值余项去掉上式中的积分余项，得第三十三页，讲稿共四十二页哦 同样，采用Ne

13、wton向后插值公式，并重复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)第三十六页，讲稿共四十二页哦从上面所列的

14、公式可以看到，步数相同的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内插法比同阶的外插法精确.第四十一页，讲稿共四十二页哦第四十二页，讲稿共四十二页哦