13讲：ch5-1基本迭代法.pptVIP

Gauss-Seidel迭代法 即： Gauss-Seidel迭代法 Gauss-Seidel迭代法 则修正项为： 用高斯-赛德尔迭代法解方程组 例5.1 Gauss-Seidel迭代法举例 如此继续下去，可通过程序观察结果 Gauss迭代法的特点 Gauss迭代法的特点 例5.2 解 Gauss-Seidel 迭代法举例 Gauss_Seidel迭代法一般程序 A={{8,-3,2},{4,11,-1},{6,3,12}}; MatrixForm[%]; b={20,33,36}; I1=IdentityMatrix[3]; DD=DiagonalMatrix[A[[1,1]],A[[2,2]],A[[3,3]]]; L=-{{0,0,0},{A[[2,1]],0,0},{A[[3,1]],A[[3,2]],0}}; DDLN=Inverse[DD-L]; G=I1-DDLN.A; MatrixForm[%]; R=Max[Abs[Eigenvalues[N[G]]]]-1 f=DDLN.b; xa={0,0,0}; Do[xb=G.xa+f//N;z=Max[Abs[xb-xa]]//N;Print[k, ,xb, ,z];xa=xb,{k,1,12}]; LinearSolve[A,b] Gauss-Seidel迭代法 Jacobi与Gau