线性方程组的矩阵求解算法.docxVIP

线性方程组的矩阵求解算法 摘要 线性方程组的矩阵求解算法，只需在约当消元法的基础上,再对方程组的 增广矩阵的行最简形进行行(列)删除和增加行，交换行等运算即可得到方程组的 解，并且这种方法既可求解有唯一解的方程组.因而算法简单，易于实现. 关键词 线性方程组;解向量;解法；约当消元法 1矩阵求解算法 设冇线性方程组二b,其增广矩阵忑如)算法的步骤如下： 第一步:利用约当消元法，把增广矩阵勿化为行最简形,设行最简形为B闷心). 若/⑴r(A)= r,则方程组无解;否则设r(A)= R,并执行以F步骤； 第二步:删除B屮的所冇零行和每一行第一个非零元素(这个非零元素一定 是1)所在的列，得到矩阵Drx(n_r+1),并记录每行的第一个非零元所在的列标，放在 一维数组f(l,L 如第，行的第一个非零元在第丿列，则= -1 L(D\ -1 L 第三步:构造矩阵，其中F = F丿 V vz — w 7(/t-r)x(n-r4-l) 第四步:对矩阵H中的行作交换运算:把H中的第，行(z = r,r-l,L 1,即从第厂 行开始直到第一行)依次与其下一行交换，使之成为第f(i)行,交换运算结果后的 矩阵记为G ,则G屮的前72-厂个斤维列向量即为方程组的一个基础解系,最后一 列向量即为方程组的一个特解； 第五步：写岀方程组的通解. 2算法证明 先证一个特殊情形，增广矩阵久的行最简形矩阵B的左上角为一r阶的单位 矩阵，即第7?行的第一个非零元的列标为匚即r(0 = z(l/r),所以设B为 10L0 1 0 L 0 1 L B = L L L 0 0 L L L L 0 0 L o S L ° C2,r+1 L L L L 1 J+i L L L L 0 0 L C2n 〃2 L L Sn山 L L 由上述算法可得H为 / C|,F+1 Cl,r+\ L 5 C2,r+1 C2,r+2 L L L L L L H = ^r,r+l Cr,r+2 L 5 d,. -1 0 L 0 0 0 -1 L 0 0 L L L L L o 0 L -1 0丿 由于r(z) = z(lzr),故从H得到G时，H屮的行不需交换位置，即G = H. 那么矩阵〃的增广矩阵的线性方程组为 西=%一6+代+]一L 5丿“ x2=d2-c2 r+[xr+l-L c2nxn, “ LLLLLLLL ￡=d「_c“]_L 一％ 令 e = / \ Cl,r+1 C2,r+1 M cr, r +1 -1 , 勺= / 、 Cl,r+2 C2,r+2 M cr, r 4- 2 0 ,J , ocn_r = / 、 5 C2n M 5 0 仏、 d2 M dr 0 0 -1 0 0 M M M M o -1 \T丿 (0丿 可以验证a,a2X 是方程组（1）所对应的齐次线性方程组的解，〃是方程组（1） 的特解，乂冷也丄 % 的后〃 的特解，乂冷也丄 % 的后 〃-厂个分量构成的向量组 ?线性 无先把它扩充成维向量组后也线性无关，所以弘也丄—线性无关,又因为 r（A）= r,所以方程组⑴的基础解系中有n-r（A）= n-r个向量，因此 少,勺丄即为方程组（1）的基础解系，特殊情形得证. 对于行最简形矩阵B为一般情形吋，可以通过若干次列交换把它变形为上述 特殊情形，但是，列交换将会导致最后结果中对应未知数的次序混乱，即在进行第 ，列与第丿?列的交换后，最后结果中齐与??次序也就被交换了，因此，在这过程中, 必须记住所进行的一切列交换，以便在最后结果屮恢复，但若使用木矩阵求解算 法，则可避免上述麻烦，为了叙述方便，还是只证一种特殊情形. （1 C|2 0 L 0 6+2 L 0 0 1 L 0 C2,r+2 L ￡ L L L O L L L L L 设 B = 0 0 L L 1 Cr,r+2 L 5 dr 0 0 L L 0 0 L 0 0 L L L L L L L L L 0 0 L L 0 0 L 0 0丿 即/(1)= 1,应）十 1(2 i 「）， 则 丫 Cd Cl,r+2 L 5 d\ n - 0 ^2,r+2 L C、2” d2 Lx — L L L L L 0 \ 1+2 L 5 ￡  / C12 0 L Cl,r+2 C2,r+2 L L L L C2n L L / -1 0 Cl,r+2 0 C2,r+2 L L L 0 0 d. ■ H = 0 Cr,r+2 L ^rn dr , G = L L L L L -1 0 L 0 0 0 Cr,r+2 L 5 dr 0 -1 L 0 0 0 -1 L 0 0 L L L L L L L L L L 0 0 L -1 0丿 o 0 L -1 0丿 现在证明G的mn-r个列向量是B所对应的方程的基础