數(shù)值分析第5章1-3節(jié)(編輯修改稿)
2025-02-15 11:24 本頁面
【文章內(nèi)容簡介】 ???將 約化為 等價(jià)的三角形方程組 (). bAx ?則可通過高斯消去法 (a) 消元計(jì)算 )1,2,1( ?? nk ?????????? ??? ),1(/ )()( nkiaam kkkkikik?),1,()()()1( nkjiamaa kkjikkijkij ??????).,1()()()1( nkibmbb kkikkiki ??????（ ） .)()2(2)1(121)()2(2)2(22)1(1)1(12)1(11???????????????????????????????????????????nnnnnnnnbbbxxxaaaaaa???????且計(jì)算公式如下： 29 (b) 回代計(jì)算 ????? ? ,/ )()( nnnnn abx).1,2,1(/)( )(1)()( ?????? ???nniaxabx iiinijjiijiii (2) 如果 為非奇異矩陣，則可通過高斯消去法 (及交 換兩行的初等變換 )將方程組 約化為 (). AbAx ?（ ） .)()2(2)1(121)()2(2)2(22)1(1)1(12)1(11???????????????????????????????????????????nnnnnnnnbbbxxxaaaaaa???????30 算法 1(高斯算法 ) 對(duì)于 sk ,2,1 ?? (1) 如果 則計(jì)算停止 ,0?kka (2) 對(duì)于 mki ,1 ??? (a) kkikikik aama /?? (b) 對(duì)于 nkj ,1 ???.* kjikijij amaa ?? 本算法用高斯方法將 約化為上梯形，且 覆蓋 ，乘數(shù) 覆蓋 . ),2,1( sk ?? AA)(kA ikm ika),1m i n (),1(R nmsmnm ???? ?A設(shè) 如果 0)( ?kkka31 當(dāng) 時(shí)，總共大約需要 次乘法運(yùn)算 . nm? 3/3n 數(shù) 稱為約化的 主元素 . )(kkka 算法 2(回代算法 ) 上三角陣， 設(shè) 其中 為非奇異 ,bUx ? nn?? RU本算法計(jì)算 的解 . bUx ? 對(duì)于 1,?ni? (1) ii bx ? 算法 1第 步需要作 次除法， 次乘法，因此，本算法 (從第 1步到第 步消元計(jì)算總的計(jì)算量 )大約需要 次乘法 (對(duì)相當(dāng)大的 ). k km? ))(( knkm ??smn ssnms ??? 2/)(3/ 23 s32 (2) 對(duì)于 nij ,1 ???jijii xuxx *?? (3) iiii uxx /? 這個(gè)算法需要 乘除法運(yùn)算 . 2/)1( ?nn 高斯消去法對(duì)于某些簡單的矩陣可能會(huì)失敗， .01 10 ?????????A由此，需要對(duì)算法 1進(jìn)行修改， 例如 ).,2,1(0)( ??? ka kkk在什么條件下才能保證 A首先需要研究原來的矩陣 33 定理 6 約化的主元素 的充要條件 ),2,1(0)( kia iii ???是矩陣 的順序主子式 A).,2,1(0 kiD i ???,0111 ?? aD即 （ ） ).,2,1(01111kiaaaaDiiiii ???????? 證明 顯然，當(dāng) 時(shí)，定理 6成立 . 1?k 現(xiàn)設(shè)定理 6充分性對(duì) 是成立的，求證定理 6充分性對(duì) 亦成立 . 1?kk首先利用歸納法證明定理 6的充分性 . 34 設(shè) ),2,1(0 kiD i ???),1,2,1( ?? ki ?可用高斯消去法將 約化到 ， )1(A )(kA,)()()()()2(2)2(2)2(22)1(1)1(1)1(12)1(11)()1(??????????????????????knnknkkknkkknknkkaaaaaaaaaaa???????????AA且有 0)( ?iiia于是由歸納法假設(shè)有 即 ,0 )2(22)1(11)2(22)1(12)1(112 aaaaaD ??35 .)()2(22)1(11)()1(1)1(11kkkkkkkk aaaaaaD ?????? （ ） 由設(shè) ),2,1(0 kiD i ???定理 6充分性對(duì) 亦成立 . k 顯然，由假設(shè) ),2,1(0)( kia iii ???利用 ()式， ,0)( ?kkka則有 利用 ()式亦可 ).,2,1(0 kiD i ???推出 36 推論 ????? ? ,)1(11 1Da如果 的順序主子式 A ),1,2,1(0 ??? nk ?kD 則 ).,3,2(/)( nka kkk ??? ? 1kk DD37 于是對(duì) ()施行第一次消元后化為 ()， 矩陣的三角分解 下面借助矩陣?yán)碚撨M(jìn)一步對(duì)消去法作些分析，從而建 立高斯消去法與矩陣因式分解的關(guān)系 . 設(shè) ()的系數(shù)矩陣 的各順序主子式均不為零 . nn?? RA由于對(duì) 施行行的初等變換相當(dāng)于用初等矩陣左乘 ， A A這時(shí) 化為 )1(A,)2(A化為 )1(b ，)2(b, )2()1(1)2()1(1 bbLAAL ??即 38 其中 .1111131211????????????????????nmmm??L 一般第 步消元， k 化為 , )(kA )1( ?kA化為 ， )(kb )1( ?kb, )1()()1()( ?? ?? kkkkkk bbLAAL相當(dāng)于 其中 39 .1111,1????????????????????????knkkkmm???L重復(fù)這過程，最后得到 ??????? 。)()1(121 nn AALLL ? （ ） .)()1(121 nn bbLLL ?? ? 記上三角矩陣 為 ，由 ()得到 )(nA U40 ,1 11211 LUULLLA ?? ? ??? n?其中 1 11211 ? ???? nLLLL ?為單位下三角矩陣 . 這就是說，高斯消去法實(shí)質(zhì)上產(chǎn)生了一個(gè)將 分解為 兩個(gè)三角矩陣相乘的因式分解，于是得到如下重要定理，
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1