第五章 第二節(jié) Gauss-Seidel迭代法_第1頁
第五章 第二節(jié) Gauss-Seidel迭代法_第2頁
第五章 第二節(jié) Gauss-Seidel迭代法_第3頁
第五章 第二節(jié) Gauss-Seidel迭代法_第4頁
第五章 第二節(jié) Gauss-Seidel迭代法_第5頁
已閱讀5頁,還剩14頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

第二節(jié)迭代法一、迭代格式二迭代法的收斂性

三、小結(jié)一、迭代格式考慮方程組(2.1)即以后,把它們代入(2.1)中第一個(gè)方程算出對方程組(2.1)作迭代時(shí),取定初始近似值顯然,迭代格式收斂的話,則比更接近于的第一個(gè)分量所以在計(jì)算時(shí),我們不再像迭代法那樣以代入(2.1)中第二式的右邊,而是把新算出的及代入該式右邊,得到即計(jì)算下一個(gè)分量時(shí),要用到剛算出的新分量。這樣或許能收到更好的效果。按這樣方式建立的迭代格式稱為迭代格式,其一般形式為(2.2)

用矩陣表示就是

(2.3)

其中,

由(2.3)式可知,因存在,所以迭代格式(2.3)也可表示為

(2.4)我們稱為迭代法的迭代矩陣。

由(2.4)式可見,對方程組作迭代,等價(jià)于對方程組

(2.5)作迭代。

二迭代法的收斂性

定理3對于任意右端向量和初始向量,迭代法收斂的充要條件是

其中由于對方程組作簡單迭代是一回事,故由定理1有作迭代同對方程組即特征方程的根的絕對值小于1。

而由于類似于定理2,我們還可以給出如下收斂的充分條件。

(2.6)所以在實(shí)際問題中,只需求出方程的根。

定理4對于任意右端向量初始向量,迭代法收斂的充分條件是

由此定理可知,條件(1)或(2)被滿足時(shí),則迭代法與迭代法都收斂。

可以證明,當(dāng)條件(2)被滿足時(shí),迭代法比迭代法收斂得快些。

例4

分別用和迭代法解方程組解由于,故迭代法迭代法都收斂。取,首先采用迭代格式,計(jì)算求得與其精確解相比,其誤差為

再利用迭代格式,計(jì)算求得

其誤差為

從此例可以看出,當(dāng)充分條件(2)被滿足時(shí),迭代法確實(shí)比迭代法收斂快些。

然而,迭代法并不總比迭代法好。有時(shí)迭代法還比迭代法收斂得慢些,有時(shí)甚至在迭代法收斂時(shí),它卻不收斂。例5設(shè)方程組的系數(shù)矩陣為

試證明迭代法收斂,而迭代法不收斂。

證明顯然,迭代法的迭代矩陣為

因?yàn)橛啥ɡ?可知,迭代法收斂。

,則有

令又,其中

則有

故迭代法不收斂。

類似地方法,可以證明,若系數(shù)矩陣為

時(shí),迭代法不收斂,而迭代法收斂。

這個(gè)問題留給同學(xué)做練習(xí)。下面我們給出判斷Gauss-Seidel迭代法收斂的其它方法。定理5設(shè)方程組為(1)如果矩陣是對稱正定的,則迭代法收斂。

(2)如果是嚴(yán)格對角占優(yōu)的,則迭代法收斂。

證略。

由上可見,迭代法的使用與問題的特點(diǎn)有著密切的關(guān)系,使用時(shí)應(yīng)

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論