數(shù)值分析用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第1頁
數(shù)值分析用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第2頁
數(shù)值分析用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第3頁
數(shù)值分析用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第4頁
數(shù)值分析用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組_第5頁
免費(fèi)預(yù)覽已結(jié)束,剩余2頁可下載查看

付費(fèi)下載

下載本文檔

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

文檔簡介

1、作業(yè)六:分別編寫用Jacobi迭代法和Gauss-Seidel迭代法求解線性方程組Ax二B的標(biāo)準(zhǔn)程序,并求下列方程組的解。可取初始向量X(°)=(0,0,Of;迭代終止條件|x(i)x叫|<=10e6(1)(2)Jacobi迭代法:流程圖結(jié)束程序clear;clc;A=8,-l,l;2/10,01;l/l/-5;b=l;4;3;e=le-6;x0=0;0;0'n=length(A);x=zeros(n/l);k=0;r=max(abs(b);whiler>efori=l:nd=A(i,i);ifabs(d)<ewarning矩陣A輸入有誤);return;e

2、ndsum=0;forj=l:nifj"=isum=sum+A(iJ)*xO(j);endendxl(i)=(b(i)-sum)/A(iJ);endk=k+l;r=max(abs(xl-xO);xO=xl;fprintf('第%(1次迭代:,k)fprintf('n與上次計算結(jié)果的距離%fn,,r)disp(xl);ifk>100warning不收斂,);endendx=xO;程序結(jié)果(1)第1次迭代;與上次計算結(jié)果的距離;0.6000000.12500.4000-0.6000第2次迭代;與上次計算結(jié)果的距離;0.1250000.25000.4350-0.495

3、0第3次迭代;與上次計算結(jié)果的距離;0.0355000.24120.3995-0.4630第4次迭代;與上次計茸結(jié)果的距盅:0.0088500.23280.3981-0.4718第5次迭代;與上次計算結(jié)果的距離:0.0025:20.23370.4006-0.4738第6次迭代:與上次計算結(jié)果的距離:0.0006990.23430.4006-0.4731第二次迭代:與上次計算結(jié)果的距離;0,00Q1840.23420.4005-0.4730第8次迭代:與上次計算結(jié)果的距離:0.0000540.23420.4005-0.4731第9次迭代:與上次計苴結(jié)果的距離:0.0000130.23420.40

4、05-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é)果的距離:2.060000-4.46004.25002.2800第3次迭代:與上次計算結(jié)果的距離:1.505000-4.55602.74502.4670第4次迭代:與上次計算結(jié)果的距離;0.564600-3.99142.62752.0347第5次迭代;與上次計算結(jié)果的距離:0.357

5、300-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é)果的距離:0.000341-3.99993.00021.9999第16次迭代:與上次計苴結(jié)果的距離:0.000155-4.00003.00012.0000第1;次送代;與上次計算結(jié)果的距離:o.000106-4.00003.00002.0000第12次迭代:與上

6、次計算結(jié)果的距離:0.000041-4.00003.00032.0000第1次送代:與上次計算結(jié)果的距澳:0.068570-4.03033.02372.0219第;8次迭代:與上次計算結(jié)果的距禹;0.042216-4.01392.98152.0132翦9次迭代:與上次計算結(jié)果的距離:0.018637-3.99522.99001.9972第10次迭代:與上次計算結(jié)果的距離:0.0L2637-3.99543.00261.9960第11次迭代:與上次計算結(jié)果的胡禹:0.004819-4.00023.00311.9999第12次迭代;與上次計箕結(jié)果的距離:0.0031217.00123.00002.0

7、010第19次迭代:與上次計算結(jié)果的胡離:000027-4.00003.00002.0000旃加次送代;與上次計算結(jié)果的距離:J.000009-4.00003.00002.0000第21次送代:與上次計算結(jié)果的距離:1000006-4.00003.00002.0000第22次送代:與上次計算結(jié)果的距離:1000003-4.00003.00002.0000第23次送代:與上次計算結(jié)果的覽離:3.OOOOOL-4.00003.00002.0000篷乂次迭代:與上次計算結(jié)果的葩離:3.OOOOOL-4.00003.00002.0000Gauss-Seidel迭代法:SJ?clear;clc;%A=8

8、,-14;2,10,01;l,lz-5;%b=l;4;3;A=5,2,l;l,4,22-3,10;b=-12;20;3;m=size(A);ifm=m(2)erro4矩陣A不是方陣力endn=length(b);%初始化N=0%迭代次數(shù)L=zeros(n);%分解A=D+L+U,D是對角陣,l.是下三角陣,U是上三角陣U=zeros(n);D=zeros(n);G=zeros(n)6G=-inv(D+L)*Ud=zeros(n/l)56d=inv(D+L)w,bx=zeros(n4);fori=l:n%初始化L和Uforj=l:nifi<jL(iJ)=A(iJ);endifi>jU(iJ)=A(iJ);endendendfori=l:n%初始化DD(iJ)=A(iJ);endG=-inv(D+L)*U;%初始化Gd=(D+L)b;%初始化d%迭代開始xl=x;

溫馨提示

  • 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)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論