[科普中國]-列主元消去法

針對問題

高斯消去法從第k步到第k+1步的消元過程，必須滿足條件。而這個元素即被稱為第k步的主元(素)。顯然，高斯消去法是按方程排列的自然順序產(chǎn)生主元的，這樣，一旦出現(xiàn)計算就歸于失敗，
而且即使，但若其絕對值很小，也將會因用它作除數(shù)，引起其他元素的數(shù)量級及舍人誤差急劇增大，導(dǎo)致最終計算結(jié)果不可靠。為了避免在高斯消去法應(yīng)用中可能出現(xiàn)的這類問題，就發(fā)展形成了列主元、全主元等多種消去法。這些方法的基本點在于對高斯消去法的過程作某些技術(shù)性修改，全面或局部地選取絕對值最大的元素為主元素，從而構(gòu)成了相應(yīng)的主元(素)消去法。列主元(素)消去法以處理簡單、相對計算量小的特點，在各類主元消去法中得到最為廣泛的應(yīng)用。1

基本思想列主元消去法的基本思想是：在進(jìn)行第 步消元時，從第k列的 及其以下的各元素中選取絕對值最大的元素，然后通過行變換將它交換到主元素 的位置上，再進(jìn)行消元。2

過程概述高斯消去法中第k-1步消元得到的結(jié)果可由分塊矩陣記為

式中，b1、b2為對應(yīng)于右端常數(shù)列的兩個子塊，而

與高斯消去法不同的是，列主元消去法在第k步消元之前，先在它的第k列主對角元 及其下方的所有元素中(亦即Akk的第一列元素中)選出絕對值最大的元素 作為這一列的主元，即

對式(1)，作初等行變換

這樣，就把 換到主對角元位置上。經(jīng)過選列主元與行交換之后，再如高斯消去法一樣作行的消元變換。上述過程從第一步消元開始執(zhí)行，即k=1，2，…，n一1，這就構(gòu)成了列****主元消去法。綜上可見，列主元消去法僅需矩陣行的交換，因此不發(fā)生未知數(shù)排列次序的調(diào)換。而全主元消去法則要花費相當(dāng)多的時間選取主元，且還需記錄系數(shù)矩陣列交換信息。1

算法本算法將消元結(jié)果覆蓋A，乘數(shù)mik沖掉 方程組解x存放在b內(nèi)。1

1) 對于k=1，2，…，n-1，依次做到第5)步；

2) 按列選主元，即確定ik，使

3) 如果 則輸出無解信息，停機(jī)；

4) 如果ik=k，則轉(zhuǎn)向5)，否則換行

5) 消元過程

6) 如果 則停機(jī)；

7) 回代求解

8) 輸出結(jié)果x。

功能將矩陣A化為上三角矩陣，然后回代求解線性方程組Ax=b。

MATLAB程序function z= Gausselimpiv (A,b, ep)

[m，n]=size(A);

if m~=n

disp(’輸入錯誤，系數(shù)矩證陣只能是方陣')

end

if n~=length (b)

disp('輸入錯誤，常數(shù)項的個數(shù)應(yīng)與方程的個數(shù)相同')

end

if nargin==2

ep= eps ;

end

for k=1:n-1

p=A(k,k);I=k;

for i=k:n

if abs (A (i, k)) > abs (p)

end

end

if p