第3章 解線性方程組的數(shù)值解法2-lu分解法

第3章 解線性方程組的數(shù)值解法2-lu分解法

ID:19537453

大?。?00.50 KB

頁數(shù):64頁

時間:2018-10-03

第3章 解線性方程組的數(shù)值解法2-lu分解法_第1頁
第3章 解線性方程組的數(shù)值解法2-lu分解法_第2頁
第3章 解線性方程組的數(shù)值解法2-lu分解法_第3頁
第3章 解線性方程組的數(shù)值解法2-lu分解法_第4頁
第3章 解線性方程組的數(shù)值解法2-lu分解法_第5頁
資源描述:

《第3章 解線性方程組的數(shù)值解法2-lu分解法》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、3.2矩陣的三角分解法我們知道對矩陣進行一次初等變換,就相當于用相應的初等矩陣去左乘原來的矩陣。因此我們這個觀點來考察Gauss消元法用矩陣乘法來表示,即可得到求解線性方程組的另一種直接法:矩陣的三角分解。3.2.1Gauss消元法的矩陣形式3.2.2Doolittle分解Doolittle分解若矩陣A有分解:A=LU,其中L為單位下三角陣,U為上三角陣,則稱該分解為Doolittle分解,可以證明,當A的各階順序主子式均不為零時,Doolittle分解可以實現(xiàn)并且唯一。A的各階順序主子式均不為零,即Doolittle分解Doolittle分解

2、Doolittle分解Doolittle分解Doolittle分解Doolittle分解例題例題例題例題例題Doolittle分解3.2.3對稱矩陣的Cholesky分解在應用數(shù)學中,線性方程組大多數(shù)的系數(shù)矩陣為對稱正定這一性質(zhì),因此利用對稱正定矩陣的三角分解式求解對稱正定方程組的一種有效方法,且分解過程無需選主元,有良好的數(shù)值穩(wěn)定性。對稱矩陣的Cholesky分解A對稱:AT=AA正定:A的各階順序主子式均大于零。即對稱矩陣的Cholesky分解由Doolittle分解,A有唯一分解對稱矩陣的Cholesky分解定理3.2.4設A為對稱正定矩

3、陣,則存在唯一分解A=LDLT,其中L為單位下三角陣,D=diag(d1,d2,…,dn)且di>0(i=1,…,n)對稱矩陣的Cholesky分解證明:對稱矩陣的Cholesky分解對稱矩陣的Cholesky分解對稱矩陣的Cholesky分解推論:設A為對稱正定矩陣,則存在唯一分解其中L為具有主對角元素為正數(shù)的下三角矩陣。對稱矩陣的Cholesky分解證明:Cholesky分解的求法Cholesky分解的求法Cholesky分解的求法Cholesky分解法Cholesky分解法缺點及優(yōu)點優(yōu)點:可以減少存儲單元。缺點:存在開方運算,可能會出現(xiàn)根

4、號下負數(shù)。改進Cholesky分解法改進的cholesky分解A=LDLT改進的cholesky分解改進的cholesky分解改進的cholesky分解算法改進的cholesky分解算法例題例題例題例題A=LDLT分解,既適合于解對稱正定方程組,也適合求解A為對稱,而各階順序主子式不為零的方程組而對A=LLT只適合于對稱正定方程組3.2.4三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法三對角方程組求解的追趕法其計算工作量為5n-4次乘除法。工作量小,其實現(xiàn)的條件為qi不為零。有以下定理可得證三對

5、角矩陣求解的充分性條件。解三對角矩陣線性方程組的追趕法程序框圖3.3矩陣求逆矩陣求逆矩陣求逆為使求逆過程不斷提高求解精度,因此增加選主元工作,最常用的是選列主元求逆。因此增加一個數(shù)組Z(n),記錄選主元的交換號,最后在消元工作完成后,根據(jù)Z(n)對A中的元素進行相應的列交換,得到A-1Gauss—Jordan原地求逆法算法(原地求逆法)例題例題例題例題例題例題

當前文檔最多預覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當前文檔最多預覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學公式或PPT動畫的文件,查看預覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負責整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時聯(lián)系客服。
3. 下載前請仔細閱讀文檔內(nèi)容,確認文檔內(nèi)容符合您的需求后進行下載,若出現(xiàn)內(nèi)容與標題不符可向本站投訴處理。
4. 下載文檔時可能由于網(wǎng)絡波動等原因無法下載或下載錯誤,付費完成后未能成功下載的用戶請聯(lián)系客服處理。