研究生数值剖析12高斯 赛德尔Gauss Seidel迭代法.pptVIP

研究雅可比迭代法，我们发现在逐个求 ( 1) k X ? 的分量时，当计算到 ( 1) k i x ? 时 , 分量 ( 1) ( 1) 1 1 , , k k i x x ? ? ? 都已经求得，而仍用旧分量 ( ) ( ) 1 1 , , k k i x x ? 计算 ( 1) k i x ? 。由于新计算出的分量比旧分量准确些， ( 1) ( 1) 1 1 , , k k i x x ? ? ? 求出，马上就用新分量 ( 1) ( 1) 1 1 , , k k i x x ? ? ? 代替雅可比迭代法中 ( ) ( ) 1 1 , , k k i x x ? 来求 ( 1) k i x ? , 这就是高斯 - 赛德尔 (Gauss-Seidel) 迭代法。 2 高斯 - 赛德尔（ Gauss-Seidel ）迭代法 因此设想一旦新分量 高斯 - 赛德尔迭代公式如下： ( 1) ( ) ( ) ( ) 1 12 2 13 3 1 1 11 ( 1) ( 1) ( ) ( ) 2 21 1 23 3 2 2 11 ( 1) ( 1) ( 1) ( 1) ( ) ( ) 1 1 2 2 , 1 1 , 1 1 ( 1) ( 1 1 1 ( ) 1 ( ) 1 ( ) 1 ( k k k k n n k k k k n n k k k k k k i i i i i i i i i in n i ii k n n nn x a x a x a x b a x a x a x a x b a x a x a x a x a x a x b a x a x a ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1) ( 1) ( 1) 2 2 , 1 1 ) k k k n n n n n a x a x b ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （ 5 ） 其矩阵表示形式为 ( 1) 1 ( 1) ( ) ( ) k k k X D LX UX b ? ? ? ? ? ? 现将 ( 1) k X ? 显式化，由 ( 1) ( ) ( ) k k D L X UX b ? ? ? ? 得 ( 1) 1 ( ) 1 ( ) ( ) k k X D L UX D L b ? ? ? ? ? ? ? 令 1 ( ) G B D L U ? ? ? （ 称为高斯 - 赛德尔（ Gauss-Seidel ）迭代矩阵）， 1 ( ) G f D L b ? ? ? 则得 ( 1) ( ) k k G G X B X f ? ? ? 为高斯 - 赛德尔迭代法的矩阵表示形式。 1 ( ) G B D L U ? ? ? 0 G I B ? ? ? 1 ( ) 0 I D L U ? ? ? ? ? 1 ( ) ( ) 0 D L D L U ? ? ? ? ? ? 1 ( ) 0 D L ? ? ? ( ) 0 D L U ? ? ? ? 上式左端为将系数矩阵 A 的对角线及对角线 以下元素同乘以 λ 后所得新矩阵的行列式。 我们用定理 2 来判断高斯 - 赛德尔迭代公式是否 收敛，需要考虑高斯 - 赛德尔迭代矩阵 的特征方程 即 将上式写成 由于 所以 例 9 用高斯 - 赛德尔迭代法解方程组 1 2 3 1 2 3 1 2 3 10 2 3 2 10 15 2 5 10 x x x x x x x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 解： 相应的高斯 - 赛德尔迭代公式为 ( 1) ( ) ( ) 1 2 3 ( 1) ( 1) ( ) 2 1 3 ( 1) ( 1) ( 1) 3 1 2 0.2 0.1 0.3 0.2 0.1 1.5 0.2 0.4 2 k k k k k k k k k x x x x x x x x x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 取迭代初值 (0) (0) (0) (0) 1 2 3 ( , , ) (0,0,0) T T X x x x ? ? 按此迭代公式进行迭代，计算结果为 k ( ) 1 k x ( ) 2 k x ( ) 3 k x 0 1 2 3 4 5 0 0.3 0.8804 0.9843 0.9978 0.9997 0