第9章微分方程初值问题的数值解法-2教程.ppt

* 五、线性多步法 第九章 微分方程初值问题的数值解法 1. 线性多步法的一般理论 线性多步法的基本思想: 如果充分利用前面多步的信息来预测 yn+k , 则可以期望获得较高的精度. 记 y(xn)的近似值为 yn , xn = x0+ nh , 并记 fn = f (xn , yn ), 则 k 步线性多步法的一般形式为: 其中 均为常数, 不全为零. 由于上式两端可同乘一常数, 故一般可设 , 则上式可写成: 或写成 注: k=1时, 即为一种单步法. 如k=1, α0= -1, β0=1 ,β1=0 就是Euler方法; k=1, α0= -1, β0= β1=1/2 就是梯形方法. 其中 任意给出, s = 0,1,2,…,迭代到满足给定的精度要求为止. 对于隐式情形的公式( ), 由于f (x , y)一般是非线性函数, 很难求得 yn+k 的显式表达式, 故常用迭代法求解: 可以证明, 当f (x , y)关于 y 满足Lipschitz条件(或 )时, 若 , 则为隐式方法; 若 , 则为显式方法. 上式称为k 步线性法. 只要