章线性方程组的迭代解法.docVIP

第14章 线性方程组的迭代解法 Jocobi迭代法 Jacobi迭代算法 1. 输入变量个数n、初值向量x（0）、迭代精度eps、系数矩阵A、常数项b 和迭代最大次数nmax 2 For i=1,2,…,n 如果|aii|eps1,则输出“迭代失败”提示并终止 3. Bj ( E-D-1A 4. gj ( D-1b 5. For k=1,2,…,nmax x (Bj.x0+ gj 如果||x-x0||eps ,输出解向量x ，终止；否则x0( x 6. 如果||x-x0||eps ,输出迭代失败，终止。 Jacobi 迭代法程序 Clear[a,b,x]; nmax=500; n=Input[“线性方程组阶数n=”]; a=Input[系数矩阵A=]； b=Input[常数项b=]； x0=Input[输入迭代初值向量x0]； eps1=0.000001; eps=Input[输入精度控制eps=]； Do[If[Abs[a[[i,i]]]eps1,t1=1;Return[],t1=0],{i,1,n}]; If[t1==1,Print[Jacobi迭代法失效], d=DiagonalMatrix[Table[a[[i,i]],{i,1,n}]]; d1=Inverse[d]; bj=IdentityMatrix[