4_YQJ_Ch2_4方组迭代Hh.pptVIP

§2.4 解线性方程组的迭代法 ( Iterative Methods for Solving Linear Systems );解线性方程组的迭代法 ; 在用直接法解线性方程组时，要对系数矩阵不断变换。如果方程组的阶数很高，则运算量将会很大并且大量占用计算机资源，因此对线性方程组 要求找寻更经济、适用的数值解法。;如果能将线性方程组(1)变换为;这种方式就称为迭代法，以上过程称为迭代过程. 迭代法产生一个序列： 如果其极限存在，即 则称迭代法收敛，否则称为发散.;例 用迭代法求解方程组，误差不超过1e-4.;解 现将方程改写;;;依此类推，得方程组满足精度的解为x12：;2.4.1 迭代法的一般形式及其收敛性;将（2.23）式写成迭代格式 上迭代过程称为迭代法，G 称为迭代矩阵。;问题：x* 是否为方程组 Ax = b 的解？;证 依据范数的等价定义，可设;性质 谱半径的性质 （1）（特征值上界） 即A 的谱半径不超过A 的任何一种范数。 （2）如果A 为对称矩阵，则 ;迭代法收敛的条件;注:要检验一个矩阵的谱半径小于1，比较困难，所以，我们希望用别的办法判断收敛性. ;定理2.9 (充分条件) 如果矩阵的某种范数||G|| 1, 则迭代收敛，且有下列误差估计：;(1);写成矩阵形式：;令;2.4.2 Jacobi 迭代法 ( Jacobi Iterative Methods );依此类推，线性方程组Ax = b可化为;;;;;由定理2.8, 有 定理2.10 Jacobi 迭代法收敛的充要条件是?(GJ) 1. 由定理2.9, 得到 定理2.11 若 || GJ || 1，则 Jacobi 迭代法收敛. (充分条件);注：为了得到更方便的收敛判别方法，先来讲一些预备知识。 ;例：;引理2.1 若 n 阶方阵 A 是主对角按行(或按列)严格占优 阵，则 A 非奇异。;定理2.12 若方程组 Ax=b 的系数阵 A 是主对角按行(或按列)严格占优阵，则用 Jacobi 迭代法求解必收敛。;当 A 是按列严格占优阵时，有;det(m I -GJ) = det(m I + D-1( L + U )] ) = det{mD-1[D+m-1 ( L + U )]} = mn ?detD-1 ? det [D+m-1 ( L + U )] (*);例6 证明：用 Jacobi 迭代法求解下方程组必收敛, 并求解, 要求 || xk- xk-1|| ∞≤10-5.;因 ,故 J-迭代法收敛.;例7 证明： 用 J-迭代法求解下方程组不收敛.;其特征方程为 ;例8 判别,下列方程组使用 J-迭代法求解是否收敛?; 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; ( suc