數(shù)值分析課后習(xí)題_第1頁
數(shù)值分析課后習(xí)題_第2頁
數(shù)值分析課后習(xí)題_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡(jiǎn)介

1、6.1方程式5x12x3=-12-x1x2x3=202 x1-3x2x3=3(a )利用雅可比迭代法、高斯-塞德爾迭代法,考察求解該方程式的收斂性(b )如果用雅可比迭代法和高斯-塞德爾迭代法求解該方程式,則要求迭代在xk 1-xk10-4時(shí)結(jié)束。解: (a )系數(shù)矩陣在行的嚴(yán)格對(duì)角上占優(yōu)勢(shì),因此雅可比法和高斯塞德爾法收斂。(b )雅可比法的迭代格式為:x1(k1)=-25x2k-15x3k-125 x2(k1)=14x1k-12 x3k5x3(k1)=-15x1k 310 x2k 310取x (0)=(1,1,1 ) t,重復(fù)18次要求精度x (18 )=(-3.99999964,2.999

2、99739,1.999999 ) t高斯-塞德爾的迭代形式是x1(k1)=-25x2k-15x3k-125 x2(k1)=14x1k1- 12 x3k5x3(k1)=-15x1k 1310 x2k 1310取x (0)=(1,1,1 ) t,代替8次要求精度。x (8)=(-4.000036,2.9999,2.000003 ) t建立兩個(gè)方程式(a ) x 10.4 x 20.4 x3=10.4 x1 x 20.8 x3=20.4 x 10.8 x2 x3=3(b ) x12x2- 2x3=1x1 x2 x3=12 x1 x2 x3=1試著考察解決這個(gè)方程式的雅科比迭代法和高斯賽德爾迭代法的收

3、斂性。解(a )雅可比法的迭代矩陣bj=d-1 llu=0- 0.4-0.4-0.40-0.8-0.4-0.80I-BJ=(-0.8)(2 0.8-0.32 )由于(BJ)=1.09282031,所以雅可比迭代法不收斂。高斯塞德爾迭代矩陣Bs=(D-L)-1U=0-2202-3002I-Bs=(-2)2,(Bs)=21高斯-賽德爾迭代法不收斂。limkAk=A的命中條件證明對(duì)于任何向量x都是limkAkx=Ax。證明所需的條件是limkAk=A,limkxijk=aijAk-A0(k) .對(duì)于任意的xAK-ax0 (k)即AkxAx,limkAkx=Ax充分的條件對(duì)于任意的xRn設(shè)定為AkxA

4、x(k)xi=(0,0,1,0,0 ) t (I=1,2,n )Akxi=(x1ik,x2ik,xnik)TAxi (k)Axi=(a1i、a2i、ani)TPSPS (j=1,2,n; I=1、2、n )即AkA,limkAk=A假設(shè)Ax=b,其中a是對(duì)稱正定的,試試問該方程式的雅可比迭代法是否必定收斂的練習(xí)問題2(a )方程式。解并不一定是一定的. 這是因?yàn)楣庾V半徑(BJ )不一定小于1關(guān)于練習(xí)題2(a )、a對(duì)稱,另外因?yàn)?=11,2=0.840,3=a=0.2960,所以a是正常的,但是雅可比迭代法沒有收斂。用SOR法解方程式(分別取緩和因子=1.03、=1、=1.1 )。4x1-x2

5、=1,-x1 4x2-x3=4-x2 4x3=-3。為了準(zhǔn)確地解x *=12,1,-12T .需要在x*-x(k)510-6處結(jié)束迭代,并針對(duì)每個(gè)值確定迭代次數(shù)。解決SOR方法的迭代形式如下x1(k1)=x1k(14-x1k 14 x2k ) x2k1=x2k(114 x1k1- x2k 14 x3k ) x3(k1)=x3k(-3414 x2k1- x3k )如果=1.03并且初始值x (0)=(0,0,0 ) t,則重復(fù)5次要求高精度X5=(0.5000035,0.9999989,-0.5000003)T用SOR法解方程式(=0.9 )5 x 12 x2x3=-12-x1x2x3=202 x1-3x2x3=3在xk 1-xk10-4處要求反復(fù)結(jié)束。SOR迭代的形式如下x1 (k1 )=x1k(-125-x1k-25x2k-15x3k ) x2 k1=x2k(514 x 1k1-x2k-12x3k ) x3 (k1 )=x3k(310-15x1k 1310 x2

溫馨提示

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

最新文檔

評(píng)論

0/150

提交評(píng)論