14章线性方程组及迭代解法.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[n]-d1.a; gj=d1.b; Do[x=bj.x0+gj; err=Max[Abs[x-x0]]; Print[x=,x//N, i=,i, err=,err//N]; If[N[err]eps,Break[],x0=x], {i,1,nmax}]; If[err=eps,Print[迭代失败 ]] ] 说明：本程序用于求线性方程组Ax=b的解。程序执行后，先通过键盘输入线性方程组阶数n、系数矩阵A、常数项b、迭代初值向量x0和输入精度控制eps，程序即可给出每次迭代的次数和对应的迭代向量序列x(K),其中最后输出的结果即为所求的根。如果迭代超出500次还没有求出满足精度的根则输出迭代失败提示，如果出现主对角线元素aii=0给出Jacobi迭代法失效提示。 程序中变量说明： x0：存放初始向量和迭代过程中的向量x（k） x: 存放迭代过程中的向量x（k+1） nmax:存放迭代允许的最大次数 err：存放误差||x-x0||( t1：临时变量 注：迭代最大次数可以修改为其他数字。 例题与实验 用Jacobi 迭代法解如下线性方程组 5x1+2x2+3x3= -12 -x1+4x2+2x3= 20 2x1-3x2+10x3= 3 要求误差||x(k+1)-x(k)||(10-4，并用取不同初值的方法实验观察迭代收敛的情况。 解：执行Jacobi迭代法程序后在输入的4个窗口中按提示分别输入 3，{{5, 2, 1}, {-1, 4, 2}, {2, -3, 10}}，{-12, 20, 3}，{0, 0, 0}，0.0001 每次输入后用鼠标点击窗口的“OK”按扭，得如下输出结果。 x={-2.4, 5., 0.3} i=1 err=5. x={-4.46, 4.25, 2.28} i=2 err=2.06 x={-4.556, 2.745, 2.467} i=3 err=1.505 x={-3.9914, 2.6275, 2.0347} i=4 err=0.5646 x={-3.85794, 2.9848, 1.88653} i=5 err=0.3573 x={-3.97123, 3.09225, 1.96703} i=6 err=0.113286 x={-4.03031, 3.02368, 2.02192} i=7 err=0.0685705 x={-4.01386, 2.98146, 2.01316} i=8 err=0.042216 x={-3.99522, 2.98995, 1.99721} i=9 err=0.0186374 x={-3.99542, 3.00259, 1.99603} i=10 err=0.0126367 x={-4.00024, 3.00313, 1.99986} i=11 err=0.0048186 x={-4.00122, 3.00001, 2.00099} i=12 err=0x={-4.0002, 2.9992, 2.00025}