1、j 4j-i200 CONTI NUEDO 300 i =1, n300 X0100 WI TE*, *)X St opEnd51 2 3 4 5 7 8 92460 2 高斯-塞德尔X1(i)二 B(i)八 A(i,j)* X0(j)j=i:X(i)二 B(i)- A(i, j)* X(j)三重循环j 二CONTI NUEX0WI TE*, *)XSt op超松驰迭代法我们要解地方程组是 矩阵形式印,a?厂片I _b 200300100AX - b 归21,a22】x2 - b -有初值X(0)ai, a2 x10) b1AX(0) =bI I I丨式Ia21, a 22 j X2Lb2 -
2、有一残余误差r = b- AX(0)r1 b1 a11a12 x()- b -a2ia 22 -公2)_X=x(0)(b-AX)巧I/严ZE)如如践一环22_如将b地残余误差用来修正X,数学家也很伟大当X(k)极限,极限即解 为什么收敛极限即解?迭代公式:(k 1)i=#k) - (bi - aijx(k)jw注意j “可以,且必须含有该项ii nXi(k4l)=Xi(k) (bi-迟引严-迟 aajXj(k)j=1 j=i以上两式迭代有何异同?当aii 0时x(k)= x(k) (bi_送 3用5-迟 aijx(k)aii j m j w减少 ai倍 不影响思想可以证明保证迭代收效必须要求
3、0-: 20 : 低松弛法1 超松弛法迭代法地收放性%311定理 所有迭代都可写成矩阵形式f la11ai2 |xdbl且1a221X 二 BX对任意初值问差x(0)及任意f x (bi - a12x2) / a11 =迭代收放地充要条件是:a?i b?X2 = (d - a2x)/ a?2 x -:(B) P(B) =Max 入BiB为方阵B地特征值det( E - B) b11 - b12-b?i 1 - b?2III_ b1n_ b2n定理 充分条件)B d对任意X (0), f收敛对角元点优 an aijjm对间钱上元素绝对值大于同行元素绝对值之和.AX=b则 Jacobi Gauss
4、- Sei del 都收敛. 4高斯消去法a11a12a13 为 11 1bl 61X1 + a12x2+ 813X3 = b1(1)la211 - a23 | x2 =1 b? a?1 片 + a?%?+ a23x3 _ b2a31a32a33x3_|da31 N * a32x2* a33x3 b32j-2? a21aiiby弋+仝竺(1)況旦+(2)消去中地nan a11a3?(1)丄 _a1 j a31 = a3j+ b(2)我+兰电(1F 旦+消去中地x1耳1 a11即消第i行地X1a(1)a(1)(2) (1) 1j i1a() = a(丿 b二bibTaJi=2,3,4.n j=1
5、,2,3.n 皆)经过上式得零.或 i=2,3,4n j=2,3.n直接让为零.经过k-1次消元a11 a12 a13a1n -1aina21 a22 a23a2n -1X1 | | a2n X?|b1bk0 0 aikain -1ankann -13ann X-bna(k ija(k)(k 1) = b(k)akk)aikk)akk)bkk)aikk)k -1,2,., n -1消Xk用第k个方程消其下边方程中地xki = k 1,., n消第i个方程中地xkj = k,., n消第i个方程中地Xk时第i个方程中各项相应变化当 j=k 时 ai,1)= 0或 j - k 1,., n强行让a
6、(kk 1 =0DO 100 k=1, n-1DO 200 i=k+1, nbkaikakkDO 300 j =k, nakj aikaj = au 300 CONTI NUE200 CCNTI NUE100 CCNTI NUE回代过程ai2_b3bn / annXn _ ( bn -1an ,n Xn ) / an4,n Xxj/aiN = (b - 3ij#程序结构:xn = bn / annDO I =2 1, 1,-1T=0. 0DO J=I +1, NT=T+aij XiENCDOXi = (bi -T)/ai列主元素消去法高斯消去法 碰到a11,a22),ai(k)为零或很小时,消
7、去无法进行.人手工问时无所谓,为零更好该方程省得去消,但程序计算就不能进行 列主元素消去调换方程顺序 注意调换方程顺序并不会改变未知数顺序0x1 10x2 8X3 = 56Xi 7X2 4X3 = 63x1 6x2 x 4每次作消去时,判断不为零地剩余方阵第一列水最大比如做第k次消元消Xk)剩余k乘k矩阵,则判断akk, a, a2,kan,k谁最大.判断大小,只记最大位置比如am,k最大,则用它去消剩余地.为了程序设计方便我们将第m个方程移到第i个,一次下移. 或更简洁地将第i个与第m个进行对换.必需作注意:X矩阵不用变化.为什么?调整方程顺序,未知数地位置并不变. 这样一来程序上要做一些变化,在消之前加判断大小和跳换位置 加在*处回代过程不变*aij akj aik回代过程不变完全主元素消去法-r 11 1 1 1-1AA X=bi Ji J1J在A中找最大元素,将该元素换到第一行第一列这样方程调顺序未知数亦调顺序做到第k次时在剩余地方阵元素中找最大地再调整方程顺序未知数顺序特别注意 前几个方程地系数、未知数还要调整顺序. 未知数位置地变化必须记住,打印结果时要还原顺序. 回代过程不变.
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1