Chapter8_2_单步法收敛性与稳定性.pptVIP

在讨论收敛性之前，先介绍局部截断误差、整体截断误差的定义及其他们之间的关系; 1、局部截断误差;2、整体截断误差;3、局部截断误差与整体截断误差的关系;再利用引理1就可得到;结论：若单步法的局部截断误差为O（hP+1) 整体截断误差为O（hP);定义; 例5 Euler方法是收敛的. 证明： 由于Euler方法是一阶方法，且其增量函数?(x, y, h) =f (x, y)．而初值问题是要求函数f (x, y)对y满足Lipschitz条件的， 故Euler方法收敛. ; 例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)式． 若值y i有一个扰动?i，那么用(35)式计算，得到的值yi+1就会产生一个偏差?i+1．若记 ;或 由此可知，单步法(35)式绝对稳定的条件是 ;（1）首先考察Euler方法的稳定性． 此时增量函数? (x, y,h) =f (x, y)=y，因而有 ;（2）讨论改进Euler方法的稳定性．此时增量函数 ;（3）经典Runge-Kutta方法的稳定性． 此时增量函数;代入后得 于是有 ;如果仅限于讨论?是实数的情形，则上述几个单步法的稳定性条件可分别简化为 Euler法稳定性条件：?2 ? ? h? 0， 改进Euler法稳定性条件：?2 ? ? h ? 0， 经典Runge-Kutta法稳定性条件：?2.785 ? ? h ? 0.;例7 对初值问题 取h = 0.1和0.2，用经典Runge-Kutta方法求解.;隐式单步法的稳定性讨论. （1）考察向后Euler方法. 对于试验方程(46)，其向后Euler法的公式为 yi+1 = y i + h? yi+1 解出yi+1，有 ，从而得到误差(扰动)公式为 由此得到绝对稳定的条件为 其绝对稳定区域是以1为半径、以1为中心的圆外部，如图8-5所示.; （2）讨论梯形法的稳定性． 对于试验方程(46)，相应的梯形法公式为 解出yi+1，有 ; 从以上的分析讨论可以看到，隐式方法的稳定性比显式方法好，这也是隐式方法的主要优点！ ; 8.1.4 线性多步法 求解初值问题的数值方法都是“步进式”的，即求解过程从初值y0开始，顺着节点的排列次序，一步一步地向前推进．所以，在计算yi+1 时，前面的i + 1 个值y0, y1, ?, y i 都是已知的．如果在计算yi+1 时能充分利用这些已有的信息，而不是像单步法中那样，只用其前一步的值y i，则可望构造出精度高，但计算量小的求解公式．线性多步法就是基于这一思想发展起来的，其计算公式可表示为 其中， ，而?j ,?j 都是常数.; 线性多步法(50)式的实质是用若干节点处的函数值及导数值的线性组合来逼近y(xi+1)的值．由于在计算yi+1时需要用到其前k + 1个值：y i , yi?1, yi?2,?, 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 , xi?k+1,?,xi?1, 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代