A=-L-UDB8.2Jacobi迭代法与Gauss-Seidel迭代法Jacobi.pptVIP

Jacobi 迭代的分量形式： 定理 设线性方程组x=Bx+g有惟一解，那么逐次逼近法对任意初始向量X０收敛的充分必要条件是迭代矩阵B的谱半径 ? ( B ) 1 定理 如果A是按行(列)严格对角占优的矩阵，那么 Jacobi和G－S迭代法都收敛 定理 A按行严格对角占优，则Gauss-Seidel迭代收敛。 证 设 是 任一特征值，x是相应特征向量。设 若 则 注 （2）Jacobi迭代法和Gauss-Seidel迭代法的收敛性没有必然的联系： 用Jacobi迭代法求解不收敛，但用 Gauss-Seidel法收敛 系数矩阵A是正定矩阵，因此用 Gauss-Seidel法收敛 * 将 改写为 等价形式 建立迭代 　 从初值 出发，得到序列 如何建立迭代格式？　 收敛速度？ 向量序列的收敛条件？ 误差估计？ A = -L -U D B Jacobi 迭代法 BG-S 作A的另一个分裂： 其迭代格式的矩阵形式为 Gauss-Seidel 迭代法 写成分量形式： … … … … 定理 设A是有正对角元的n阶对称矩阵，那么 Jacobi迭代法收敛的充分必要条件是 A和2D-A同为正定矩阵. 证明： 略 定理 设A是不可约对角占优矩阵， 那么Jacobi迭代法与G－S迭代法都收敛. 定理 设A是n阶正定矩阵，那么，G－S迭代法收敛. 即当Gauss-Seidel法收敛时，Jacobi法可能不收敛； 而Jacobi法收敛时， Gauss-Seidel法也可能不收敛。 （1）Jacobi迭代法和Gauss-Seidel迭代法的 迭代矩阵不同： BJ =D-1(L+U), B G-S = (D-L) -1U 用Jacobi迭代法求解收敛， 但用 Gauss-Seidel法不收敛。 BJ的特征值为0，0，0， 而BG－S的特征值为0，2，2 不是正定矩阵，因此用 Jacobi迭代法不收敛 A是有正对角元的 n阶对称矩阵