数值分析－第七章常微分方程数值解.ppt

常微分方程表示方法 在微分方程中, 自变量的个数只有一个, 称为 常微分方程. 自变量的个数为两个或两个以上的微分方程 叫偏微分方程。 常微分方程解法回顾 7.5 亚当姆斯方法 7.5.1 亚当姆斯格式 龙格-库塔方法是一类重要算法，但这类算法在每一步都需要先预报几个点上的斜率值，计算量比较大。考虑到计算yi+1之前已得出一系列节点上 的斜率值，能否利用这些已知值来减少计算量呢？ 这就是亚当姆斯（Adams）方法的设计思想。 设用xi,xi+1两点的斜率值加权平均作为区间 上的平均斜率，有计算格式 （7.21） 选取参数λ，使格式（7.21）具有二阶精度。 将 在xi处Taylor展开 代入计算格式(7.21)化简,并假设, 因此有 与y(xi+1)在xi处的Taylor展开式 相比较, 需取 才使格式(7.21)具有二阶精度。这样导出的计算格式 称之为二阶亚当姆斯格式。类似地可以导出三阶亚当姆斯格式。 和四阶亚当姆斯格式。 (7.22) 这里和下面均记 上述几种亚当姆斯格式都是显式的，算法比较简单，但用节点 的斜率值来预报区间 上的平均斜率是个外推过程，效果不够理想。为了进一步改善精度，变外推为内插，即增加节点xi+1的斜率值来得出 上的平均斜率。譬如考察形如 （7.23） 的隐式格式,设(7.23)右端的 Taylor展开有 可见要使格式(7.23)具有二阶精度,需令 , 这样就可构造二阶隐式亚当姆斯格式 其实是梯形格式。类似可导出三阶隐式亚当姆斯格式 和四阶隐式亚当姆斯格式 （7.24） 7.5.2 亚当姆斯预报-校正格式 参照改进的欧拉格式的构造方法，以四阶亚当姆斯为例，将显式（7.22）和隐式（7.24）相结合，用显式公式做预报，再用隐式公式做校正，可构成亚当姆斯预报-校正格式 （7.25） 预报： 校正： ( 7.5 ) (7.5)式的右端含有未知的yi+1,它是一个关于yi+1的函数方程,这类数值方法称为隐式方法。相反地,欧拉法是关于yi+1的一个直接的计算公式， 这类数值方法称为显式方法。 例7.3 用梯形公式求下面初值问题的解 在x=0.01上的值y(0.01) y= y y (0) = 1 解：取h=0.01, 由梯形公式得 y=ex e0.01=1.010050167 7.3.3 两步欧拉公式 对方程 的两端在区间上 积分得 ( 7.6 ) 改用中矩形公式计算其积分项，即 代入上式,并用yi近似代替式中y(xi)即可得到两步欧拉公式 ( 7.7 ) 前面介绍过的数值方法,无论是欧拉方法,还是梯形方法，它们都是单步法,其特点是在计算yi+1时只用到前一步的信息yi;可是公式(7.7)中除了yi外,还用到更前一步的信息yi-1,即调用了前两步的信息,故称其为两步欧拉公式 7.3.4． 欧拉法的局部截断误差 衡量求解公式好坏的一个主要标准是求解公式的精度, 因此引入局部截断误差和阶数的概念。 定义7.1 在yi准确的前提下, 即 时, 用数值方法计算yi+1的误差 , 称为该数值方法计算时yi+1的局部截断误差。 对于欧拉公式，假定 ，则有 而将真解y(x)在xi处按二阶泰勒展开 因此有 定义7.2 数值方法的局部截断误差为 ,则称这种数值方法的阶数是P。步长(h1) 越小,P越高, 则局部截断误差越小,计算精度越高。欧拉公式的局部截断误差为 , 欧拉方法仅为一阶方法。 两步欧拉公式比欧拉公式精度也是高一个数值方法,设 , 前两步准确,则两步欧拉公式 把y(xi-1)在xi处展开成Taylor级数，即 将 在xn处按二阶泰勒展开 略去余项 用近似值yn代替y(xn)有 由 再将y(xi+1)在xi处进行泰勒展开 (7.8) (7.9) 所以, 由(7.8)和(7.9)可得两步欧拉公式的局部截断