第六章線性方程組迭代法

第六章線性方程組迭代法第一頁，共二十八頁，編輯于2023年，星期五對(duì)方程組做等價(jià)變換如：令，則則，我們可以構(gòu)造序列若同時(shí)：所以，序列收斂與初值的選取無關(guān)第二頁，共二十八頁，編輯于2023年，星期五定義：（收斂矩陣）定理：即：矩陣B為收斂矩陣當(dāng)且僅當(dāng)B的譜半徑<1由知，若有某種范數(shù)則迭代收斂.第三頁，共二十八頁，編輯于2023年，星期五1.Jacobi迭代法第四頁，共二十八頁，編輯于2023年，星期五格式很簡(jiǎn)單：第五頁，共二十八頁，編輯于2023年，星期五2.Gauss－Seidel迭代法在Jacobi迭代中，使用最新計(jì)算出的分量值第六頁，共二十八頁，編輯于2023年，星期五

迭代矩陣記A=-L-UD第七頁，共二十八頁，編輯于2023年，星期五易知，Jacobi迭代有第八頁，共二十八頁，編輯于2023年，星期五

迭代矩陣第九頁，共二十八頁，編輯于2023年，星期五JacobiiterationGauss-Seideliteration計(jì)算x(k+1)時(shí)需要x(k)的所有分量，因此需開兩組存儲(chǔ)單元分別存放x(k)和x(k+1)計(jì)算xi(k+1)時(shí)只需要x(k)的i+1~n個(gè)分量，因此x(k+1)的前i個(gè)分量可存貯在x(k)的前i個(gè)分量所占的存儲(chǔ)單元，無需開兩組存儲(chǔ)單元.第十頁，共二十八頁，編輯于2023年，星期五迭代公式:例用Gauss-seidel迭代法解方程組Ax=b計(jì)算結(jié)果:第十一頁，共二十八頁，編輯于2023年，星期五3逐次超松弛迭代法(SOR)記則可以看作在前一步上加一個(gè)修正量。若在修正量前乘以一個(gè)因子，有對(duì)Gauss－Seidel迭代格式整理得引入松弛因子第十二頁，共二十八頁，編輯于2023年，星期五寫成分量形式，有第十三頁，共二十八頁，編輯于2023年，星期五迭代矩陣第十四頁，共二十八頁，編輯于2023年，星期五

SOR方法收斂的快慢與松弛因子的選擇有密切關(guān)系.但是如何選取最佳松弛因子,即選取=*,使(B)達(dá)到最小,是一個(gè)尚未很好解決的問題.實(shí)際上可采用試算的方法來確定較好的松弛因子.經(jīng)驗(yàn)上可取1.4<<1.6.第十五頁，共二十八頁，編輯于2023年，星期五第十六頁，共二十八頁，編輯于2023年，星期五第十七頁，共二十八頁，編輯于2023年，星期五第十八頁，共二十八頁，編輯于2023年，星期五4迭代法的收斂性定義

設(shè)有矩陣序列及，如果則稱收斂于，記為第十九頁，共二十八頁，編輯于2023年，星期五一些關(guān)于收斂的定義及定理定理定理設(shè)，則其中為的譜半徑。定理(迭代法基本定理）設(shè)有方程組對(duì)于任意初始向量及任意，解此方程組的迭代法收斂的充要條件是第二十頁，共二十八頁，編輯于2023年，星期五定義稱為迭代法的收斂速度.定理(迭代法收斂的充分條件)如果方程組的迭代公式為，且迭代矩陣的某一種范數(shù)，則1）迭代法收斂，即對(duì)任取，有

2）3）實(shí)際計(jì)算中,通常利用作為控制迭代的終止條件.不過要注意,當(dāng)時(shí),較大,盡管已非常小,但誤差向量的?？赡芎艽?迭代法收斂將是緩慢的.第二十一頁，共二十八頁，編輯于2023年，星期五特別的,Jacobi迭代法收斂G-S迭代法收斂SOR迭代法收斂第二十二頁，共二十八頁，編輯于2023年，星期五

定理若SOR方法收斂,則0<<2.

證設(shè)SOR方法收斂,則(B)<1,所以|det(B)|=|12…n|<1而det(B)=det[(D-L)-1((1-)D+U)]

=det[(E-D-1L)-1]det[(1-)E+D-1U)]

=(1-)n于是|1-|<1,或0<<2第二十三頁，共二十八頁，編輯于2023年，星期五

定理設(shè)A是對(duì)稱正定矩陣,0<<2時(shí),則解方程組

Ax=b的SOR方法收斂.

第二十四頁，共二十八頁，編輯于2023年，星期五注意的問題（2）Jacobi迭代法和Gauss-Seidel迭代法的收斂性沒有必然的聯(lián)系：即當(dāng)Gauss-Seidel法收斂時(shí)，Jacobi法可能不收斂；而Jacobi法收斂時(shí)，Gauss-Seidel法也可能不收斂。（1）Jacobi迭代法和Gauss-Seidel迭代法的迭代矩陣不同：BJ=D-1(L+U),BG-S=(D-L)-1U第二十五頁，共二十八頁，編輯于2023年，星期五用Jacobi迭代法求解不收斂，但用G-S法收斂。用Jacobi迭代法求解收斂，但用G-S法不收斂。BJ的特征值為0，0，0，而BG－S的特征值為0，2，2第二十六頁，共二十八頁，編輯于2023年，星期五系數(shù)矩陣A是正定矩陣，因此用Gauss-Seidel法收斂不是正定矩陣，因此用Jacobi迭代法不收斂A是有正對(duì)角元的n階對(duì)稱矩陣第二十七頁，共二十八頁，編輯于2023年，星