jacobi迭代法实验MATLAB程序数值分析.docVIP

12 1220183311136112x(0)=z 12 12 20183311136112 x(0)= z r(o)、 Al ⑴、 Al x(々)= A\ 彳) 彳) 8x, - 3x2 + 2x3 = 20, 4; +1 lx2 - x3 = 33， 6x, + 3x2 +1 2x3 = 36. 得近似解。精确解为x*=[3, 2，1]’ 解：对方程进行移项就得 8%, - 3x 8%, - 3x2 + 2x3 = 20, 4%, +1 lx2 - =33,= 6x, +3x2 +12x3 = 36. Xj = — (3x’ — 2x3 + 20), -^2 = (_4^i + x3 + 33), — (-6a*| - 3x? + 36). 记为Ax=b,或写为x=Box+f,其中 _2 i 丄 n 0 取初始值x(G =(0,0,0)\代入原方程组可得=(2.5,3,3)\再将 把它代入可得x2).反复利用这个汁算过程，得到一向量序列和一般的计算公式(迭 代公式) x^+,)= (3x^-2%^+ 20)/8, +l)=(-4x，V)+33)/11， x^+,) = (-6》⑷-3x[k) +36)/12. 简写为氏;vu)+/j表示迭代次数a = 0，l，2，...).迭代到第10次有 x(10) = (3.000032,1.999838,0.9998813)y; (io)0.000187 (io) 0.000187 (f(10)=x(10)-x： oo 从此例看出，由迭代法产生的向量序列x(k)逐步逼近方程组的精确解X*. 6.1常用迭代法 定义1 ( i)对于给定的方程组用公式xu+u = Bx(k} +/，逐步代入 求近似解的方法称为迭代法(或称为一阶定常迭代法,这里5 与k无关). (ii)如果lim，存在(记为/),称此迭代法收敛，显然x*就是方 程组的解，否则称此迭代法发散. 迭代法的流程图为： 为初始向量，A⑼，_4°)，…，xf); r是判断条件，即|x, -xQ| ￡时停止运行 k是循环次数。 X = Sx + /中带入初始值,然后赋给X , ①Jacobi迭代法 从第一个方程解出XI,第二个方程解出X2,…，记成 bnx2 + Z?13x3 + …+ bXnxn +g, x2=Z?2lXi+ 么3心+…+ b2nxn +g2 Xn= bnXXX + bn2X2 + bn3X3 + ??? + bn,n-lXn +Sn 用矩阵写法即x=Bx+g, B的对角元皆零，可拆成 B=L+U L是B下三角部分，C/是B上三角部分Jacobi迭代法如下述. 任取初始近似x(0)，对*=1，2,…计算 r (⑽. xi 一x2(A-+1)= b2\X r (⑽. xi 一 x2 (A-+1) = b2\Xl b^x-,{k} + H—.+ + Z?23x3(A) + ??? + b2^nk)+S2 (A+l) (k) (A) = bn,x,{k} +bn2x^K, +bn3x；Kf + + bn,n-lXn (A) 直至II I 预定的精度.用矩阵记号，即任取初始近似x(())，对仁1，2,…计 算 x(k+])=Bx{k)+g, 直至II I ?通常对迭代法限最大迭代次数也是必要的).Jacobi迭代法 的流程图为： 开始 在以上的流程图屮，先读入数据，即先输入系数矩阵A，常数向量6,初始值，停 止条件和最大循环次数。图屮％是，在我们迭代公式中的^+|)。k是循环次数， N是最大循环次数。 例2. 利用Jacobi方法求方程组 10x, -x2 - 2x3 = 7.2 —Xj +10xo = 8.3 -%| — + 5x^ — 4.2 的近似解。 解把原方程改为 x, = 0. lx2 +0.2x3 +0.72 x2 = O.lx, +0.2x3 +0.83 = 0.2x, +0.2x2 +0.84 任取初始近似x(0),对hl，2,…计算 +|) = 0. lx?) + 0.2x^ + 0.72 4+|) =0.1;)+0.2;^)+0.83 x广=0.2^^ + 0.2x^ +0.84 直至I || 预定的精度. 此即Jacobi迭代法.计算结果如下表. k v (k) v (k) 0 0 0 0 1 0.720000000000 0.830000000000 0.840000000000 2 0.971000000000 1.070000000000 1.150000000000 3 1.057000000000 1.157100000000 1.248200000000 4 1.085350000000 1.185340000000 1.282820000000 5 1.095098000