NA2_4线性方程组的迭代法分解.ppt

2.4 解线性方程组的迭代法 /* Iterative Methods for Solving Linear Systems */;2.4.1 迭代法的一般形式;将上式写为迭代过程 这种迭代过程称为逐次逼近法，B称为迭代矩阵。 ;问题： 是否是方程组Ax=b的解？;逐次逼近法收敛的条件;要检验一个矩阵的谱半径小于1比较困难，所以我们希望用别的办法判断收敛性 ;定理4;①;;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 */;;下面从另一个角度来说明;;注：这二种方法都存在收敛性问题。在讨论收敛性之前我们先来讲一些预备知识和有关的定理 ;二、 可约矩阵与不可约矩阵;三、 有关性质;定理9 设A是不可约对角占优矩阵， 那么A是非奇异矩阵.;2.4.5 Jacobi迭代法和Gauss-Seidel迭代法的收敛性;定理11 如果A是按行(列)严格对角占优的矩阵，那么Jacobi和G－S迭代法都收敛;注意的问题;举例;例1：分别Jacobi及G-S迭代法解下列线性方程组。初值均取(0，0，0)T;%用雅可比法解例1 a=[9,-1,-1;-1,8,0;-1,0,9]; D=-(a-triu(a)-tril(a)); % D=diag(diag(a)); L=-(tril(a)-D); U=-(triu(a)-D); xo=[0;0;0];bo=[7;7;8]; ep=0.0001;dx=1;k=0; while dxep k=k+1; x=D\(L+U)*xo+D\bo; dx=abs(norm(x)-norm(xo)); xo=x; end k,x;%用雅可比法解例1 a=[9,-1,-1;-1,8,0;-1,0,9]; bo=[7;7;8]; ep=0.0001;dx=1;k=0; 结果： k=7, 迭代如下： 0.7778 0.8750 0.8889 0.9738 0.9722 0.9753 0.9942 0.9967 0.9971 0.9993 0.9993 0.9994 0.9998 0.9999 0.9999 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000;从计算结果可以看到:如果两种迭代法都收 敛,那么Jacob迭代法慢于G-S迭代法.这 个结论具有一般意义.