數(shù)值分析5用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第1頁
數(shù)值分析5用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第2頁
數(shù)值分析5用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第3頁
數(shù)值分析5用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第4頁
數(shù)值分析5用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第5頁
已閱讀5頁,還剩7頁未讀 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

1、程組Ax二B的標準程序,并求下列方程組的解可取初始向量X(o)=(O,0,Of;迭代終止條件|x(k+1,-x(k)|二10e-68-11210-111-5X11x24x3.3521-142.2-310.X1x2x3-二1220.3.Jacobi迭代法:流程圖判斷b中的最大值有沒有比 謀差大1給X賦M值進行迭代求出x,弱到100次還沒到,警告不收斂程序clear;clc;A斗&?:U;2J0,01;:L,:L,?5;b=l;4;3;e=le-6;x0=0;0;01;n=length(A);x=zeros(n/l);k=0;r=max(abs(b);whilerefori=l:nd=A(iJ);

2、ifabs(d)100warningf不收斂,);endendx=xO;程序結(jié)果(1)第1次迭代;與上次計算結(jié)果的3人:0.6000000.12500.4000-0.6000第2次迭代;與上次計算結(jié)果的距善:0.1250000.25000.4350-0.4950第3次迭代;與上次計茸結(jié)果的距八:0.0355000.24120.3995-0.4630第4次迭代:與上次計算結(jié)果的距S:0.0088500.23280.3981-0.4718第5次迭代;與上次計算結(jié)果的3人:0.0025720.23370.4006-0.4738策6次迭代:與上次計算結(jié)果的距八:0.0006990.23430.4006

3、-0.4731第:次迭代:與上次計算結(jié)果的距離;0.00Q1840.23420.4005-0.4730第8次迭代;勻上次計算結(jié)果的距離:0.0000540.23420.4005-0.4731第9次迭代:與上次計算結(jié)果的距離:0.0000130.23420.4005-0.4731第10次迭代:與上次計算結(jié)果的距離:0.0000040.23420.4005-0.4731第11次迭代;與上次計算結(jié)果的距離;0?0000010.23420.4005-0.4731(2)第i次迭代:與上次計算結(jié)果的距離:5.000000-2.40005.00000.3000第2次迭代:與上次計尊結(jié)果的3人:2.06000

4、04?46004.25002.2800第3次迭代:與上次計算結(jié)果的距離:1?505000-4.55602.74502.4670第4次迭代:與上次計算結(jié)累的距離;0.564600-3.99142.62752.0347第5次迭代;與上次計算結(jié)果的距離:0?357300-3.85792.98481.8865第6次迭代:與上次計算結(jié)果的距高:0?113286-3.97123.09221.9670第13次迭代:與上次計算結(jié)果的距八:0.001023-4.00022.99922.0002第14次迭代:與上次計算結(jié)果的距離;0?000626-3.99972.99981.9998第15次迭代:與上次計算結(jié)果的

5、距離:0.000341-3.99993.00021.9999第16次迭代:與上次計茸結(jié)果的距離:0.000155-4.00003.00012.0000第訂次迭代;與上次計算結(jié)果的距離:o.000106?4.00003.00002.0000第18次迭代:與上次計算結(jié)果的距離:0.000041-4.00003.00032.0000第 12 次迭代;第;次迭代:與上次計算結(jié)果的距離:0.068570-4.03033.02372.0219第8次迭代:與上次計算結(jié)果的匪離;0.042216-4.01392.98152.0132第9次迭代:與上次計算結(jié)果的距離:0.018637-3.99522.99001

6、.9972第10次迭代:與上次計算結(jié)果的范離:0.0L2637-3.99543.00261.9960第口次迭代:與上次計算結(jié)果的徒離:0.004819-4.00023.00311.9999與上次計磬結(jié)果的施禹:0.003121-00123.00002.0010第19次迭代:與上次計算結(jié)果的昵離:3.000027-4.00003.00002.0000第如次迭代;與上次i+算結(jié)果的柜胡:3.000009-4.00003.00002.0000第21次迭代:與上次i+算結(jié)果的距韶:1000006-4.00003.00002.0000金周22次送代:與上次計算結(jié)果的捱離:1000003?4.00003.

7、00002.0000第23次迭代:與上次計算結(jié)果的距離:3.OOOOOL-4.00003.00002.0000與上次計算結(jié)果的拒離:3.OOOOOL-4.00003.00002.000016加次逑代:Gauss-Seide迭代法:SJ?clear;clc;%A=8r14;240,01;14z-5;%b=l;4;3;A=5,24;-l/4,2;2,-340;b=-12;20;3;m=size(A);ifm(l)r=m(2)errorC矩陣A不是方陣J;endn=length(b);%初始化N=0,%迭代次數(shù)L=zeros(n)八6分解A=D+L+U,D是對角陣,L是下三角陣,U是上三角陣U=zeros(n);D=zeros(n);G=zeros(n)A6G=-inv(D+L)*Ud=zeros(n/l)56d=inv(D+L)w,bx=zeros(n4);fori=l:n%初始化L和Uforj=l:nifijU(iJ)=A(iJ);endendendfori=l:n%初始化DD(iJ)=A(iJ);endG=-inv(D+L)*UA6初始化Gd=(D+L)

溫馨提示

  • 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)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論