数值分析--第9章_常微分方程初值问题数值解法-(精品课件).pptVIP

类似可得三阶及四阶R-K方法的E(hλ)分别为 由|1+hλ|1可得到相应的绝对稳定域. 当λ为实数时，则得绝对稳定区间，它们分别为 三阶显式R-K方法： 四阶显式R-K方法： 　　从以上讨论可知显式R-K方法的绝对稳定域均为有限域，都对步长h有限制. 如果h不在所给的绝对稳定区间内，方法就不稳定. 例4 分别取h=0.1及h=0.2，用经典的四阶R-K方法(3.1)计算初值问题 解 本例λ=-20,hλ分别为-2及-4，前者在绝对稳定区间内，后者则不在，用四阶R-K方法计算其误差见下表 xn 0.2 0.4 0.6 0.8 1.0 h=0.1 h=0.2 0.93×10-1 4.98 0.12×10-1 25.0 0.14×10-2 125.0 0.15×10-3 625.0 0.17×10-4 3125.0 　　从以上结果看到，如果步长h不满足绝对稳定条件，误差增长很快. 　　对隐式单步法,可以同样讨论方法的绝对稳定性,例如对后退欧拉法，用他解模型方程可得 故 由|E(hλ)|1可得绝对稳定域为|1-hλ|1，这是以(1,0)为圆心，1为半径的单位圆外部. 故方法的绝对稳定区间为-∞hλ0. 当λ0时，则0h∞，即对任何步长均为稳定的. 对隐式梯形法，它用于解模型方程(4.8)得 故 对Re(λ)0有|E(hλ)|1，故绝对稳定域为μ=hλ的左半平面，绝对稳定区间为-∞hλ0，即0h∞时隐式梯形法均是稳定的. 9.5 线性多步法 在逐步推进的求解过程中，计算yn+1之前事实上已经求出了一系列的近似值y0,y1,?,yn，如果充分利用前面多步的信息来预测yn+1，则可以期望会获得较高的精度. 这就是构造所得线性多步法的基本思想. 构造多步法的主要途径基于数值积分方法和基于泰勒展开方法，前者可直接由方程(1.1)两端积分后利用插值求积公式得到. 本节主要介绍基于泰勒展开的构造方法. 9.5.1 线性多步法的一般公式 如果计算yn+k时，除用yn+k-1的值，还要用到yn+i (i=0,1,?,k-2)的值，则称此方法为线性多步法. 一般的线性多步法公式可表示为 其中yn+1为y(xn+1)的近似，fn+i=f(xn+i, yn+i), 这里xn+i=xn+ih，αi, βi为常数， α0及β0不全为零，则称(5.1)为线性k步法，计算时需先给出前面k个近似值y0,y1,?,yk-1，再由(5.1)逐次求出yk,yk+1,?. 如果βk=0，则(5.1)称为显式k步法，这时yn+k可直接由(5.1)算出；如果βk≠0, 则(5.1)称为隐式k步法,求解时与梯形法(2.7)相同, 要用迭代法方可算出yn+k. (5.1)中系数αi及βi可根据方法的局部截断误差及阶确定，其定义为 定义7 设y(x)是初值问题(1.1), (1.2)的准确解，线性多步法(5.1)在xn+k上局部截断误差为 若Tn+k=O(hp+1)，则称方法(5.1)是p阶的，p≥1则称方法(5.1)与方程(1.1)是相容的. 　　由定义7，对Tn+k在xn处泰勒展开，由于 代入(5.2)得 其中 　　若在公式(5.1)中选择系数αi及βi，使它满足 由定义可知此时所构造的多步法是p阶的，且 称右端第一项为局部截断误差主项, cp+1称为误差常数. 　　根据相容性定义, p≥1，即c0=c1=0，由(5.4)得 故线性多步法(5.1)与微分方程(1.1)相容的充分必要条件是(5.6)成立. 　　显然，当k=1时，若β1=0，则由(5.6)可求得 α0=1，β0=1. 此时公式(5.1)为 即为欧拉法. 从(5.4)可求得c2=1/2≠0，故方法为1阶精度，且局部截断误差为 这和第2节给出的定义及结果是一致的. 　　对k=1，若β1≠0，此时方法为隐式公式，为了确定系数α0,β0,β1，可由c0=c1=c2=0解得α0=1, β0=β1=1/2.于是得到公式 即为梯形公式. 　　由(5.4)可求得c2=-1/12，故p=2，所以梯形法是二阶方法，其局部截断误差主项是 这与第2节中讨论是一致的. 　　对k≥2的多步法公式都可利用(5.4) 确定系数αi,βi,并由(5.5)给出局部截断误差，下面只就若干常用的多步法导出具体公式.　 9.5.2 阿当姆斯显式与隐式公式 　　p362—自学. 9.5.3 米尔尼方法与辛普森方法 　　p366—自学. 9.5.4 汉明方法 　　p367—自学. 9.5.5 预测-校正方法 　　p368—自学. 9.5.6 构造多步