第10章_常微分方程数值解..pptVIP

* 结束 8.8.2 高阶常微分方程 以两阶常微方程为例： 则可令z=y′,化为一阶方程组求解： 完全套用8.6.1的各解法即可解决. * 结束 当μh处于某区间(μ0, 0)误差传递不扩大，则称区间(μ0, 0)为绝对稳定区间。显然μ0越小,方法越稳定,h的选择余地就越大。例如：欧拉法的绝对稳定区间为(-2,0),四阶经典R-K法为(-2.785,0) * 结束 * 结束 一般讲,隐式方法比显式方法稳定性好.若选择某公式计算时发现误差扩大现象,应减少步长h,以期获得稳定的结果. 下面列出已介绍过的方法的稳定区间 * 结束 欧拉法 (-2,0) 欧拉隐式法 (-∞,0) 欧拉中点法 — 梯形法 (-∞,0) 梯形公式预估校正法 (-2,0) 二阶R-K法 (-2,0) 三阶R-K法 (-2.51,0) 四阶RK法 (-2.78,0) 米尔尼公式 — 哈明公式 (-2.6,0) 阿达姆斯外推公式 (-0.3,0) 阿达姆斯内插公式 (-3,0) 辛普森公式 — * 结束 §8.5 线性多步法 单步法只利用前一步的结果,只要给出初值,就能开始计算.但也因为它只利用前一步的值,为了提高精度就要计算一些非结点处的函数值,增加了计算量.R-K法就是通过这一途径提高精度的.下面介绍的线性多步法,在求yi+1时,不仅用到yi的值, 还用到前若干步的yi-1,…,yi-k的值,这些值都是已知的,因此可在计算量增加不多的情况下提高精度. 8.5.1 用待定系数法构造线性多步法 线性多步法的一般形式是： 或写为： * 结束 其中αj,βj(j=0,1,…,k)都是实常数,且αk≠0, ｜α0｜+｜β0｜≠0, fi+j= f (xi+j,yi+j), j=0,1,…,k, 由(8.15)可看出要计算yi+k,要利用它前面的k个值yi,yi+1,…,yi+k-1,又因为(8.15)关于yi+j和 f i+j都是线性组合,所以这一类方法都称为线性k步法.欧拉法,隐式欧拉法和梯形法都是线性一步法,欧拉中点公式是线性二步法. 把y(x+jh)和y’(x+jh)作Taylor展开： * 结束 其中 * 结束 若选择αj,βj,使C0=C1=…=Cp=0,Cp+1≠0,则 将x=xi代入上式,设αk=1,并注意到y?(xi+jh)=f(xi+jh,y(xi+jh)),可推出 即 * 结束 设yi=y(xi),yi+1=y(xi+1),…,yi+k-1=y(xi+k-1),记左端为 yi+k,并舍去最后两项,(8.17)变为： 就是一种p阶的线性k步方法,Cp+1hp+1y(p+1)(xi)称为局部截断误差的主项.当βk=0时,是显式方法,当βk≠0时是隐式方法. 下面构造几个实用的线性多步法公式 例4 形如： yi+4= -α0yi+h(β1fi+1+β2fi+2+β3fi+3)的线性4步法公式,试确定α0,β1,β2,β3并求其局部截断误差主项. * 结束 解: 由(8.15)知α1=α2=α3=0,α4=1,β0=β4=0,因为有4个待定系数,由(8.16)写出前4个方程： 解之,得α0=-1,β1=β3=8/3,β2=-4/3 故所求的公式为 * 结束 将α0,β1,β2,β3代入C4 再求C5 公式(8.19)称为米尔尼(Milne)公式,它的局部截断误差为： 它也可以写成 * 结束 例5 试确定下列公式的系数和局部截断误差 yi+1=a yi+b yi-2+h (c fi+1+d fi+e fi-1) 解 按(8.15)和(8.16)可知k=3,α0= - b,α1=0,α2= - a,α3=1, β0=0,β1=e,β2=d,β3=c,有五个未知参数,写出前五个方程： * 结束 解之得：a=9/8, b= - 1/8, c=3/8, d=6/8, e= - 3/8代入C5 此公式称为哈明(Hamming)公式,写为： 局部截断误差记为： 它是一个四阶隐式三步法. 在(8.15)中若α0=α1=…=αk-2=0,αk=1,则称此类方法为阿达姆斯(Adams)类方法.当βk=0时,称显式,否则为隐式.这类方法也可以用以上待定系数法求出这类方法的通式为： * 结束 例6 求(8.24)中的隐式三步法 解:k=3,α0=α1=0,α2=-1,α3=1 * 结束 解之,β0=1/24,β1=-5/24,β2=19/24,β3=