第四章 第二節(jié) 矩陣三角分解法

第二節(jié)矩陣三角分解

法一、分解法二、平方根法

1、

平方根法

2、分解法

三、追趕法設矩陣A存在LU分解，即有比較兩端的第一行，有

一、分解法比較兩端的第一列，有

再比較兩端第二行的其余元素得

再比較兩端第二列的其余元素得

從上面的計算過程，不難歸納出一般的計算公式

再對計算

同時，由，按上邊順序依次求得：這個分解過程叫做分解，用這種方法求解方程組所需要的計算量和用

消元法計算需要的計算量基本相同，但這種方法把

對系數(shù)矩陣的計算和對右端項的計算分開了，這就使我們在計算系數(shù)矩陣相同而右端項不同的一系列方程組時變得特別方便。

與此類似，如果A存在LU分解，則有

其中，是以為對角元素的對角矩陣，是單位上三角矩陣，是下三角矩陣，和分解一樣,我們可以依次得矩陣和的元素和，從而將矩陣A分解為這樣的分解稱為分解。例2

用分解求解方程組和解容易驗證系數(shù)矩陣A的順序主子式全不為零，故A存在唯一的分解，由得再由可計算出

對第一個方程組，由和可求得

再由和，求得解對第二個方程組，由和得

由和得

二、平方根法

求解對稱正定方程組，在許多工程實際計算問題中經(jīng)常遇到，由于其本身的特點，使用前述的解法是不利的，應該采用適合其特點的解法。這里介紹的分解法，其運算量和存儲量較LU分解法均節(jié)省一半左右，是目前在計算機上解這類問題最有效的方法之一。設A為對稱正定矩陣，則A的各階主子式均大于零，于是有如下三角分解其中為下三角矩陣，為單位上三角矩陣。

式中為單位下三角矩陣，為非奇異的對角矩陣，記為于是成為因為，于是有由于的對角元素皆為非零元素，則也有分解式

這里，和同為下三角矩陣，和同為單位上三角矩陣，于是據(jù)分解式的唯一性，有，故若用表示對角矩陣，則有從而式又可以表為

其中已不是單位下三角矩陣了。

我們稱和為對稱正定矩陣A的分解，下面據(jù)這兩種分解式來討論方程組的計算方法，從分解式出發(fā)而得到的計算方法稱為平方根法；從分解式出發(fā)而得到的計算方法稱為分解法（改進的平方根法），這兩種方法統(tǒng)稱為分解法

。其中

和分解一樣，比較式兩端對應元素，可以依次算出，今設從第1列到第k-1列的元素均已算出，現(xiàn)在來計算L的第k列元素由關系式，有

1、

平方根法

設對稱正定陣A具有形如的分解：于是當時有

當時，有

對，重復n次的計算后，便求得矩陣L。這樣，求解對稱正定方程組可化為解三角形方程組及，它們的求解公式為：由式可以看出這說明L的元素其絕對值不會變得很大，所以分解過程中各元素是完全可以控制的，舍入誤差的增長也就是可控的。因而上述平方根法的計算過程是穩(wěn)定的。但在使用這一方法時，要完成n次開方運算，這樣將增加不少運算量，為了避免開方運算，可直接使用分解式進行計算。其中L是單位下三角矩陣，D是對角矩陣且對角元素為正數(shù)，比較式兩端的元素，有于是及L的第一列元素已算出，又于是又可算出和L的第二列元素：2、分解法

設對稱正定陣A具有形如的分解：設及L的第列元素已算出，由關系式便可得到計算矩陣D和L的計算公式其中是D的對角元素，在計算時求和號中每個元素都由三個數(shù)相乘得到，這就增加了計算量，為此引進輔助量，并將改寫為求得矩陣L和D以后，解方程組可分三步完成：1）解，得；2）解，得，即；

3）解，得。例3用分解求方程組的解。解容易驗證系數(shù)矩陣A的順序主子式都大于零，故A為對稱正定矩陣，現(xiàn)在按計算D和L：這樣就得到依次解方程組：解得解得解得三、追趕法求三次樣條插值函數(shù)及用差分法求解二階常微分方程的邊值問題等，都需要求解所謂對角占優(yōu)的三對角方程組，即

(4.1)其中由于在實際問題中遇到的這種方程組的階數(shù)都較高，所以零元占優(yōu)勢，為節(jié)省存貯和計算量，不直接用前述的（或）分解公式，而是根據(jù)矩陣的特殊形狀，推導出更為簡捷的三角分解計算和求解公式。由矩陣的特殊形狀可假設的分解具有如下形式：

(4.2)

比較兩端對應元素，得到如下計算公式：

容易驗證，分解過程（4.3）是數(shù)值穩(wěn)定的。

這樣，解方程組（4.1）可分為兩步求解和，計算公式是求解公式（4.4）和（4.5）直接用及表示，