《數(shù)值分析》實驗報告2_第1頁
《數(shù)值分析》實驗報告2_第2頁
《數(shù)值分析》實驗報告2_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡介

1、數(shù)值分析實驗報告一、問題的提出求解線性方程組的迭代法，即是用某種極限過程去逐步逼近線性方程組的精 確解的過程，迭代法是解大型稀疏矩陣方程組的重要方法。二、實驗名稱運(yùn)用MATLAB編程實現(xiàn)雅可比(Jacobi)迭代和高斯-賽德爾(Gauss-Seidel) 迭代。三、實驗?zāi)康?、熟悉了解雅可比(Jacobi)迭代和高斯-賽德爾(Gauss-Seidel)迭代的 算法。2、學(xué)習(xí)MATLAB軟件的功能。四、基本原理五、實驗環(huán)境操作環(huán)境： Windows10實驗平臺：Matlab7.1軟件六、試驗設(shè)計1、jacobi迭代法算例：課本p54頁例1程序清單Jacobi迭代法的MATLAB函數(shù)文件Jacob

2、i.m如下：function y,n=jacobi(A,b,x0,eps)if nargin=3eps=1.0e-6;elseif nargin=epsx0=y;y=B*x0+f;n=n+1;end實驗結(jié)果及分析： A=8,-3,2;4,11,-1;6,3,12; b=20,33,36; x,n=jacobi(A,b,0,0,0,1.0e-6)x =3.00002.00001.0000 n =162、Gauss-seidel 迭代法算例：課本p54頁例1程序清單：Gauss-Serdel迭代法的MATLAB函數(shù)文件gauseidel.m如下:function y,n=gauseidel(A,b,x0,eps)if nargin=3eps=1.0e-6;elseif nargin=epsx0=y;y=G*x0+f;n=n+1;end (3)實驗結(jié)果及分析： A=8,-3,2;4,11,-1;6,3,12; b=20,33,36; x,n=gauseidel(A,b,0,0,0,1.0e-6)x =3.00002.00001.0000n =9七、結(jié)果說明：高斯-賽德爾迭代法比雅可比迭代法收斂得快一些(達(dá)到相同精度所需迭代 次數(shù)較少)。八、實驗感想通過本次的實驗了解了雅可比(Jacobi)迭代和高斯-賽德爾(Gauss-

溫馨提示

  • 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

提交評論