Chapter33_线性方程的组的迭代法.pptVIP

3.3 解线性方程组的迭代法 /* Iterative Methods for Solving Linear Systems */;2.2.1 逐次逼近法(迭代格式的构造);将上式写为迭代过程 这种迭代过程称为逐次逼近法，B称为迭代矩阵。 ;问题： 是否是方程组Ax=b的解？;逐次逼近法收敛的条件;要检验一个矩阵的谱半径小于1比较困难，所以我们希望用别的办法判断收敛性 ;定理2.2.4;①;;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;Jacobi 迭代的分量形式：; Algorithm: Jacobi Iterative Method Solve given an initial approximation . Input: the number of equations and unknowns n; the matrix entries a[ ][ ]; the entries b[ ]; the initial approximation X0[ ]; tolerance TOL; maximum number of iterations Mmax. Output: approximate solution X[ ] or a message of failure. Step 1 Set k = 1; Step 2 While ( k ? Mmax) do steps 3-6 Step 3 For i = 1, …, n Set ; /* compute xk */ Step 4 If then Output (X[ ]); STOP; /* successful */ Step 5 For i = 1, …, n Set X 0[ ] = X [ ]; /* update X0 */ Step 6 Set k ++; Step 7 Output (Maximum number of iterations exceeded); STOP. /* unsuccessful */;;下面从另一个角度来说明;;注：这二种方法都存在收敛性问题。在讨论收敛性之前我们先来讲一些预备知识和有关的定理 ;二、 可约矩阵与不可约矩阵;三、 有关性质;定理2.2.9 设A是不可约对角占优矩阵， 那么A是非奇异矩阵.;2.2.5 Jacobi迭代法和Gauss-Seidel迭代法的收敛性;对于BJ=I-D-1A,有;的特征值为2－?i;定理（2.2.12，2.2.13，2.2.14 ）如果A是按行(列)严格对角占优的矩阵，那么Jacobi和G－S迭代法都收敛;定理2.2.13 A按行严格对角占优，则Gauss-Seidel迭代收敛。 证 设 是 任一特征值，x是相应特征向量。设 若 则;注意的问题;举例;系数矩阵A是正定矩阵，因此用 Gauss-Seidel法收敛