《单步法的收敛性和稳定性.ppt》由会员分享，可在线阅读，更多相关《单步法的收敛性和稳定性.ppt（42页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、单步法的收敛性和稳定性现在学习的是第1页，共42页1、局部截断误差在在假假设 yi=y(xi)，即即第第 i 步步计算算是是精精确确的的前前提提下下，考考虑的的截截断断误差差 Ri=y(xi+1)yi+1 称称为局局部部截截断断误差差/*local truncation error*/。假定“yi=y(xi)”称为局部化假定局部化假定现在学习的是第2页，共42页2、整体截断误差现在学习的是第3页，共42页3、局部截断误差与整体截断误差的关系现在学习的是第4页，共42页再利用引理1就可得到现在学习的是第5页，共42页结论：若单步法的局部截断误差为O（hP+1)整体截断误差为O（hP)条件：满足L
2、ipschitz条件现在学习的是第6页，共42页 若若某某算算法法对于于任任意意固固定定的的 x=xi=x0+i h，当当 h0(同同时 i )时有有 yi y(xi)，则称称该算算法法是是收收敛的。的。单步法的收敛性定义结论1：收敛整体截断误差Ei 0结论2：只要单步法(35)式是高于零阶的方法，判断单步法(35)式的收敛性就归结为验证其增量函数(x,y,h)是否满足对y的Lipschitz条件现在学习的是第7页，共42页例5 Euler方法是收敛的.证明：由于Euler方法是一阶方法，且其增量函数(x,y,h)=f(x,y)而初值问题是要求函数f(x,y)对y满足Lipschitz条件的，
3、故Euler方法收敛.现在学习的是第8页，共42页 例6 改进Euler方法是收敛的.证明 改进Euler方法是二阶方法，其增量函数为下面证明，当f(x,y)满足对y的Lipschitz条件时，(42)式中的(x,y,h)也满足对y的Lipschitz条件.现在学习的是第9页，共42页由(42)式有假定h h0(h0为定数)，并记 ，则有即(x,y,h)满足对y的Lipschitz条件，故改进的Euler方法是收敛的.现在学习的是第10页，共42页二、单步法的稳定性 上面讨论单步法的收敛性，是假定(35)式的每一步计算都是准确的，即不考虑计算中的舍入误差然而这一假定是不切合实际的，用(35)式
4、进行实际数值计算时，每一步都不可避免地含有舍入误差稳定性就是讨论计算过程中的舍入误差对最终结果的影响！现在学习的是第11页，共42页定义4 如果一种数值方法在节点xi的值y i有大小为i的扰动，而由这个扰动引起以后各节点上值y i(j i)的偏差j均满足|j|i|，则称该数值方法是绝对稳定的定的.现在学习的是第12页，共42页 考虑一般的单步法(35)式若值y i有一个扰动i，那么用(35)式计算，得到的值yi+1就会产生一个偏差i+1若记则可将yi+1视为单步法公式的准确结果用(43)式减去(35)式，得现在学习的是第13页，共42页或由此可知，单步法(35)式绝对稳定的条件是由于增量函数与
5、微分方程的右端f 有关，从而给考察单步法的稳定性带来了困难为了简化讨论，通常是用试验方程y=y (为复常数)来检验数值方法的稳定性！现在学习的是第14页，共42页（1）首先考察）首先考察Euler方法的方法的稳定性定性此时增量函数(x,y,h)=f(x,y)=y，因而有 因此对于试验方程(46)，Euler方法稳定的条件是|1+|1(47)由于可以是复数，故在h的平面上，(47)式表示以点-1为中心的单位圆及其内部区域这个区域称为Euler方法的方法的绝对稳定区域定区域现在学习的是第15页，共42页（2）讨论改改进Euler方法的方法的稳定性定性此时增量函数改进Euler方法的稳定性条件为现在
6、学习的是第16页，共42页（3）经典典Runge-Kutta方法的方法的稳定性定性此时增量函数现在学习的是第17页，共42页代入后得于是有由此得出经典Runge-Kutta方法的稳定性条件为现在学习的是第18页，共42页如果仅限于讨论是实数的情形，则上述几个单步法的稳定性条件可分别简化为 Euler法稳定性条件：2 h 0，改进Euler法稳定性条件：2 h 0，经典Runge-Kutta法稳定性条件：2.785 h 0.由上面的讨论可以看到，如果方法的绝对稳定区域或区间是有限的，那么，步长h的选取要受绝对稳定性的约束.现在学习的是第19页，共42页例7 对初值问题取h=0.1和0.2，用经典
7、Runge-Kutta方法求解.解本例中=20,h分别为-2和-4前者属绝对稳定区间-2.785,0，后者不属此区间问题的准确解为y=e20 x计算结果误差见表8-6现在学习的是第20页，共42页隐式单步法的稳定性讨论.（1）考察向后Euler方法.对于试验方程(46)，其向后Euler法的公式为 yi+1=y i+h yi+1解出yi+1，有 ，从而得到误差(扰动)公式为由此得到绝对稳定的条件为其绝对稳定区域是以1为半径、以1为中心的圆外部，如图8-5所示.现在学习的是第21页，共42页（2）讨论梯形法的稳定性对于试验方程(46)，相应的梯形法公式为解出yi+1，有由此得出绝对稳定的条件为其
8、绝对稳定区域为Re(h)0 的整个复平面当为实数时，其绝对稳定区间为h 0.现在学习的是第22页，共42页 从以上的分析讨论可以看到，隐式方法的稳定性比显式方法好，这也是隐式方法的主要优点！现在学习的是第23页，共42页 8.1.4 线性多步法 求解初值问题的数值方法都是“步进式”的，即求解过程从初值y0开始，顺着节点的排列次序，一步一步地向前推进所以，在计算yi+1 时，前面的i+1 个值y0,y1,y i 都是已知的如果在计算yi+1 时能充分利用这些已有的信息，而不是像单步法中那样，只用其前一步的值y i，则可望构造出精度高，但计算量小的求解公式线性多步法就是基于这一思想发展起来的，其计
9、算公式可表示为其中，而j,j 都是常数.现在学习的是第24页，共42页 线性多步法(50)式的实质是用若干节点处的函数值及导数值的线性组合来逼近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外插法现在学习的是第25页，共42页 1公式的推导 对方程y0
10、=f(x,y)的两边从xi到xi+1积分，得为了近似计算(51)式中的积分，我们以xi k,xik+1,xi1,xi 为插值节点，作函数f(x,y(x)的k 次插值多项式p k(x)，从而有f(x,y(x)=p k(x)+R(x),其中，R(x)为插值余项将上式代入(51)式，得现在学习的是第26页，共42页略去积分余项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)代替后所得到的表达式.
11、由于(52)式的积分是在区间xi,xi+1上进行的，而 的插值区间却为xi k,xi，即插值点x xi,xi+1，位于插值区间之外，故称(53)式为外插公式.现在学习的是第27页，共42页考虑到插值点x靠近区间xi k,xi 的最后一个节点xi，我们采用Newton向后插值公式，于是有其中 ，而为向后差分算子,即 ，等.现在学习的是第28页，共42页将(54)式代入(53)式并作变量代换，得其中(55)式称为Adams外插公式(亦有称Adams-Bashforth公式)不难看出，这是一类k+1步显式方法.现在学习的是第29页，共42页由(56)式容易计算出b j，它的前几个值见表8-7.2Ad
12、ams外插公式的局部截断误差 注意到，若假定y i j=y(xi j),j=0,1,k精确成立，则有 因此，Adams外插公式(55)的局部截断误差为现在学习的是第30页，共42页其中而i是介于xi k与xi+1之间的某个值.由(57)式可知，Adams外插公式(55)是一类k+1步k+1阶的显式方法.几个常用的Adams外插公式如下：单步法(k=0)二步法(k=1)现在学习的是第31页，共42页 三步法(k=2)四步法(k=3)现在学习的是第32页，共42页 二、Adams内插法 现在以k+2个节点x k,xik+1,xi,xi+1作为插值节点，作函数f(x,y(x)的k+1次插值多项式pk
13、+1(x)，从而有其中R(x)为插值余项去掉上式中的积分余项，得现在学习的是第33页，共42页 同样，采用Newton向后插值公式，并重复Adams外插公式的推导过程，可得其中表8-8列出了d j的前几个值.现在学习的是第34页，共42页 由于插值点x 现在是落在积分区间xi,xi+1之内，所以称公式(60)为Adams内插公式(亦称Adams-Moutton公式)容易看出，这是一类k+1步隐式方法.由(59)式可知，公式(60)的局部截断误差为其中而i是介于xi k与xi+1之间的某个值所以，Adams内插公式(60)是一类k+1步k+2阶的隐式方法.现在学习的是第35页，共42页几个常用的
14、Adams内插值公式如下：单步法(k=0)二步法(k=1)三步法(k=2)现在学习的是第36页，共42页从上面所列的公式可以看到，步数相同的Adams内插公式比外插公式在精度上要高一阶，而阶数相同的内插公式的截断误差也比外插值公式的截断误差小许多，这是内插法的优点但内插法是隐式的，求解用迭代法，因而计算量较大，这是它的缺点.例8 对初值问题分别用四步四阶Adams外插法和三步四阶Adams内插法求解.解取步长h=0.1，即N=10.现在学习的是第37页，共42页四步四阶Adams外插法的公式为将 f(x,y)=2x+y，h=0.1，xi=0.1i 代入，得现在学习的是第38页，共42页三步四阶Adams内插法的公式为将f(x,y)=2x+y，h=0.1,xi=0.1i 代入，得现在学习的是第39页，共42页 本例可以解出yi+1 使其成为显式 本例的精确解为y(x)=3ex 2x 2利用此精确解求出y1=y(x1),y2=y(x2),y3=y(x3)，并用y0,y1,y2,y3 作为上述外插公式的起步值，计算结果见表8-9.现在学习的是第40页，共42页 从表8-9可以看到，Adams内插法比同阶的外插法精确.现在学习的是第41页，共42页现在学习的是第42页，共42页