資源描述：
《數(shù)值計(jì)算方法 解線性方程組的數(shù)值解法》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、3.2矩陣的三角分解法我們知道對矩陣進(jìn)行一次初等變換，就相當(dāng)于用相應(yīng)的初等矩陣去左乘原來的矩陣。因此我們這個觀點(diǎn)來考察Gauss消元法用矩陣乘法來表示，即可得到求解線性方程組的另一種直接法：矩陣的三角分解。3.2.1Gauss消元法的矩陣形式3.2.2Doolittle分解Doolittle分解若矩陣A有分解：A=LU，其中L為單位下三角陣，U為上三角陣，則稱該分解為Doolittle分解，可以證明，當(dāng)A的各階順序主子式均不為零時，Doolittle分解可以實(shí)現(xiàn)并且唯一。A的各階順序主子式均不為零，即Doolittle分解Doolittle分解Doolittle
2、分解Doolittle分解Doolittle分解Doolittle分解例題例題例題例題例題Doolittle分解3.2.3對稱矩陣的Cholesky分解在應(yīng)用數(shù)學(xué)中，線性方程組大多數(shù)的系數(shù)矩陣為對稱正定這一性質(zhì)，因此利用對稱正定矩陣的三角分解式求解對稱正定方程組的一種有效方法，且分解過程無需選主元，有良好的數(shù)值穩(wěn)定性。對稱矩陣的Cholesky分解A對稱：AT=AA正定：A的各階順序主子式均大于零。即對稱矩陣的Cholesky分解由Doolittle分解，A有唯一分解對稱矩陣的Cholesky分解定理3.2.4設(shè)A為對稱正定矩陣，則存在唯一分解A=LDLT，其中
3、L為單位下三角陣，D=diag(d1,d2,…,dn)且di>0(i=1,…,n)對稱矩陣的Cholesky分解證明：對稱矩陣的Cholesky分解對稱矩陣的Cholesky分解對稱矩陣的Cholesky分解推論：設(shè)A為對稱正定矩陣，則存在唯一分解其中L為具有主對角元素為正數(shù)的下三角矩陣。對稱矩陣的Cholesky分解證明：Cholesky分解的求法Cholesky分解的求法Cholesky分解的求法Cholesky分解法Cholesky分解法缺點(diǎn)及優(yōu)點(diǎn)優(yōu)點(diǎn)：可以減少存儲單元。缺點(diǎn)：存在開方運(yùn)算，可能會出現(xiàn)根號下負(fù)數(shù)。改進(jìn)Cholesky分解法改進(jìn)的choles
4、ky分解A=LDLT改進(jìn)的cholesky分解改進(jìn)的cholesky分解改進(jìn)的cholesky分解算法改進(jìn)的cholesky分解算法例題例題例題例題A=LDLT分解，既適合于解對稱正定方程組，也適合求解A為對稱，而各階順序主子式不為零的方程組而對A=LLT只適合于對稱正定方程組3.2.4三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法其計(jì)算工作量為5n-4次乘除法。工作量小，其實(shí)現(xiàn)的條件為qi不為零。有以下定理可得證三對角矩陣求解的充分性條件。解三對角矩陣線性方程組的追趕法程序框圖3.3矩陣
5、求逆矩陣求逆矩陣求逆為使求逆過程不斷提高求解精度，因此增加選主元工作，最常用的是選列主元求逆。因此增加一個數(shù)組Z(n)，記錄選主元的交換號，最后在消元工作完成后，根據(jù)Z(n)對A中的元素進(jìn)行相應(yīng)的列交換，得到A-1Gauss—Jordan原地求逆法算法（原地求逆法）例題例題例題例題例題例題
