第三章线性代数(11.18第二部分)预案.pptVIP

二、Gauss-Seidel迭代法及其收敛条件 1、迭代公式 其中 cij= -aij / aii (i=1,...,n, j=1,....,n, j≠i) di=bi/aii -----分量形式 则Gauss-Seidel迭代可写成: ( D-L ) x (k+1) = U x (k) + b x (k+1) = ( D-L )-1 U x (k) + (D-L) –1 b 记 G = (D-L)-1 U (迭代矩阵） f = (D-L)-1 b 则Gauss-Seidel迭代的矩阵形式为： x (k+1) = Gx (k) + f 2.矩阵形式 若 A分解为 A= D - L – U, D, L,U 含义同上 例4：将例1改写成Gauss-seidel迭代 8x1-3x2+2x3=20 4x1+11x2-x3=33 6x1+3x2+12x3=36 jacobi迭代： Gauss-seidel迭代 例5 用赛德尔迭代法解方程组 解 将原方程组写成等价形式并构造赛德 尔迭代公式 开始迭代……. 对比上面简单迭代, 显然Gauss-seidel收敛较快 ! 下表是jacobi迭代(简单迭代)的结果 3.收敛性 充分条件1 :方程组