高斯消去法和因子表分解法

《高斯消去法和因子表分解法》由会员分享，可在线阅读，更多相关《高斯消去法和因子表分解法（7页珍藏版）》请在装配图网上搜索。

1、高斯消去法目前，电力网络方程主要用高斯消去法求解。计算机在电力系统应用的初期， 曾经因为内存容量的限制采用过迭代法求解电力网络的线性方程式组。迭代法的 致命缺点是存在收敛性问题。因此，自从稀疏技术成功地在电力系统应用之后， 迭代法几乎完全为高斯消去法所代替。高斯消去法求解线性方程式组由消去运算和回代运算两部分组成。消去运算 又叫前代运算，可以按行进行，也可以按列进行。同样，回代运算可以按行进行， 也可以按列进行。通常采用“消去运算按列进行，回代运算按行进行”的方式较 多。设有n阶线性方程组AX = B .其中矩阵A和向量B的元素可以是实数或复 数。由于消去运算只对A和B进行，因此，为了算法叙述

2、方便，把B作为第n+1列附在A之后，形成nx (n +1)阶增广矩阵：A = LabLa11a21a .12a.22 a1n a2nb -1 b2=a11a21a12a22. a1n. a2n a ,1 n+1a ,2 n+1an1a.n2annbnan1an2. annan,n+1为了方便讨论，上式中用a替代了 b ( j= 1,2,., n ) 。j ,n+1j按列消去的运算步骤如下： 第一步，消去第一列。首先，把增广矩阵A的第一行规格化为1-6 )1 a (1) a (1) a (1)12131,n+1式中：aa (1)二 s (j = 2,3,., n +1)1 j a11然后，用式(

3、1-6)所表示的行消去A的第一列对角线以下各元素a ,a，,a， 2131n1结果使A的第2到第n行其他元素化为ai (1) = a - a a (1)(j = 2,3,., n +1; i = 2,3,., n)ji ji 1 1 j式中：上标(1)表示该元素第一次运算的结果。这时矩阵A变为A : 1=AB =11a (1). a (1)a (1)121n1,n+1a (1). a (1)a (1)222n2,n+1. .a (1). a (1)a (1)n2nnn,n+1与之对应的方程组是AX二B，它与AX = B同解。矩阵未标出的元素为零,11下同。第二步,消去第二列。首先，把增广矩阵A

4、的第二行规格化为10 1 a (2) a (2)(1-7)232,n+l式中：a (1)a (2)= 2(j 二 3,4,., n +1)2 j a (1)22然后，用式(1-7)所表示的行消去A的第二列对角线以下各元素1a (1),a (1) ,., a (1) ，3242n 2结果使A的第3到n行其他元素化为1a (2)二 a a a (2)(j = 3,4,.,n +1;i = 3,4,.,n) ijiji 22 j式中：上标(2)表示该元素第二次运算的结果。这时矩阵A变为石:121a (1)a (1). a (1)a(1)12131n1,n+11a (2). a (2)a(2)232n

5、2,n+1A = A B =a (2). a (2)a(2)2 2 2333n3,n+1a (2). a (2)a(2)-n3nnn,n+1一般地，在消去第k列时要做以下的运算:a (k1)a (k) = k(j = k +1,.,n +1)(1-8)kja (k1)kka (k)二 a (k-1) a (k-1)a (k)(j = k +1,., n +1; i = k +1,., n)( 1-9)ijijikkj经过对矩阵A的N次消去运算，即k从1依次取到N按式(1-8),(1-9)运 算，使矩阵 A 对角线以下的元素全部化为零，从而得到增广矩阵1 a121A = AB =nn na (1

6、)13a (2)231a(1)a(1)1n1, n +1a(2)a(2)2n2,n+1a(3)a(3)3n3,n+11 a (n)n,n+11-10)与之对应的方程组是A X二B，即nnx + a x + a (1)x +. + a (1)x = a (1)1 1221331nn1,n +1x +a(2)x +.+a (2)x =a (2)2 2332 nn2,n +11-11)x +.+a(3)x =a (3)3 3nn3,n +1x = a(n)nn ,n +1它与原方程组AX = B同解。现在来讨论按行回代过程。对于方程组( 1-11)，回代运算自下而上进行 首先由第n个方程可知x =

7、a( n)nn ,n+1然后将x代入第n-1个方程，解出 nx = a( -1) a(n-1)xn1n1,n+1n1,nn再将 x 和 x 代入第 n-2 个方程，可解出 x 。一般地，把已求出的n 1nn2x ,x ,., x 代入第 i 个方程，即可求出i +1 i + 2nx = a - a (i)x (i=n,2,1)(1-12)ii , n +1i j jj=i+1式(1-12)就是按行回代的一般公式。因子表和三角分解法在实际计算中，经常遇到这种情况：对于方程组需要多次求解，每次仅改变 其常数项B,而系数矩阵A是不变的。这时，为了提高计算速度，可以利用因 子表求解。因子表可以理解为高

8、斯消去法解线性方程组的过程中对常数项 B 全部运算 的一种记录表格。如前所述，高斯消去法分为消去过程和回代过程。回代过程的 运算由对系数矩阵进行消去运算后得到的上三角矩阵元素确定，见式( 1-10)。 为了对常数项进行消去运算(又叫前代过程)，还必须记录消去过程运算所需要 的运算因子。消去过程中的运算又分为规格化运算和消去运算，以按列消去过程为列，由式(1-8)、(1-9)可知，消去过程中对常数项B中的第i个元素b (即 ia )的运算包括i ,n+1(i=1,2,.,n)1-13)b (k ) i=b Qik-1) - a (k-1)b (k ) ikk(k=1,2,. ,i-1)将上式中的

9、运算因子 a ,a (1)ai-2 及 a(i-1)逐彳丁放在i1i 2i ,i -1ii(1-10)的上三角矩阵元素合在一起，就得到了因子表aa (1)a (1)a (1) .a (1)111213141naa (1)a (2)a (2) .a (2)212223242naa (1)a (2)a (3) .a (3)313233343naa (1)a (2)a (3) .a (4)414243444naa (1)a (2)a (3) .a( n-1)n1n2n3n4nn其中下三角元素用来对常数项 B 进行消去(前代)运算进行消去回代运算。因子表也可以表示为如下形式：duuu .u111213

10、141nlduu .u212223242nlldu .u313233343nllld .u414243444nllll .dn1n2n3n4nn式中d=iia (i -1) iiu=ija (i ) ij( ij)l=a (i -1)(ji)三角部分，和式，上三角元素用来 b (i) b (i)二ii a (i-1)ii1-14 )1-15)ij ij不难看出，因子表中下三角部分的元素就是系数矩阵在消去过程中曾用以进 行运算的元素，因此只要把它们保留在原来的位置，并把对角元素取倒数就可以 得到因子表的下三角部分。而因子表中上三角部分的元素就是系数矩阵在消去过 程完成后的结果。对于方程组，需要多次求解，每次仅改变其常数项B而系数矩阵A是不变 的情况，应首先对其系数矩阵 A 进行消去运算，形成因子表。有了因子表，就 可以对不同的常数项 B 求解。这时，可以直接应用因子表中的元素，用下面的 公式代替式(1-13)、(1-14)，进行消去运算：b(i) = b (i) / diiiib(k) = b(k-i)一 l b (k)(i=k+l,n)iiik k用以下公式代替式(1-12 )进行回代运算:x = b (n)nnx = b (订 - u xi ii j jj=i+1参考文献王锡凡主编，方万良，杜正春.现代电力系统分析M.北京:科学出版社,2003.