差分格式的稳定性与收敛性9.pdf

差分格式的稳定性与收敛性 第一部分：介绍了差分格式稳定性与收敛性的概念； 第二部分：给出了一个二阶常微分方程边值问题的一种差分格式； 第三部分：对上述问题的差分解的稳定性与收敛性作了讨论； 第四部分：小结. 1 基本概念 所谓稳定性问题是指在数值计算过程中产生的误差的积累 和传播是否受到控制.在应用差分格式求近似解的过程中，由于 我们是按节点逐次递推进行，所以误差的传播是不可避免的， 如果差分格式能有效的控制误差的传播，使它对于计算结果不 会产生严重的影响，或者说差分方程的解对于边值和右端具有 某种连续相依的性质，就叫做差分格式的稳定性. h 差分格式的收敛性是指在步长 足够小的情况下，由它所 u 确定的差分解 能够以任意指定的精度逼近微分方程边值问题 m u(x ) {u } 的精确解 .下面给出收敛性的精确定义：设 是差分格式 m m 定义的差分解，如果当h 0 并且um x 时，有um u(x) 0 ，则 称此格式是收敛的. 2 差分方程的建立 对于二阶边值问题  Lu u q (x )u f (x ),a x b,  (1) u(a) ,u(b) ,  其中 、f (x ) C a ,b ,q (x ) 0. q(x)   N x a mh,m 0,1,,N , 将区间 分成 等份，记分点为 a,b   m 这里步长h b a .利用泰勒公式，得 N 1 2 [(u (xm1 ) 2u (xm ) u (xm1 )] u (xm ) Rm (2) h 1 h2 其中 R  u ( 4 )( ), (x ,x ) (3) m m m m1 m1 12 把式(2)代入式(1)中的微分方程，有 1 L u (x )  [(u (x ) 2u (x ) u (x )] q (x )u (x ) h m h2 m1 m m1 m m f (x m ) R m (4) R x 略去余项 ，便得到(1)式中的微分方程在内部节点 的差分 m