高斯消元法matlab實現(xiàn)_第1頁
已閱讀1頁，還剩4頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、《數(shù)值分析》實驗報告《數(shù)值分析》實驗報告一、實驗目的與要求一、實驗目的與要求1掌握高斯消去法的基本思路和迭代步驟；2培養(yǎng)編程與上機調試能力。二、實驗內容二、實驗內容1編寫用高斯消元法解線性方程組的MATLAB程序，并求解下面的線性方程組，然后用逆矩陣解方程組的方法驗證.（1）（2）1231231230.1012.3043.5551.1831.3473.7124.6232.1372.8351.0725.6433.035xxxxxxxxx?

2、???????????????12312312352828321361xxxxxxxxx??????????????2編寫用列主元高斯消元法解線性方程組的MATLAB程序，并求解下面的線性方程組，然后用逆矩陣解方程組的方法驗證.（1）（2）1231231230.1012.3043.5551.1831.3473.7124.6232.1372.8351.0725.6433.035xxxxxxxxx????????????????123123

3、12352828321361xxxxxxxxx??????????????三MATLAB計算源程序計算源程序1.用高斯消元法解線性方程組用高斯消元法解線性方程組的MATLAB程序程序bAX?輸入的量：系數(shù)矩陣和常系數(shù)向量；Ab輸出的量：系數(shù)矩陣和增廣矩陣的秩RARB方程組中未知量的個數(shù)ABn和有關方程組解及其解的信息.Xfunction[RARBnX]=gaus(Ab)B=[Ab]n=length(b)RA=rank(A)RB=rank

4、(B)zhica=RBRAifzhica0disp(請注意：因為RA~=RB，所以此方程組無解.)returnendifRA==RBifRA==ndisp(請注意：因為RA=RB=n，所以此方程組有唯一解.)X=zeros(n1)C=zeros(1n1)fp=1:n1fk=p1:nm=B(kp)B(pp)B(kp:n1)=B(kp:n1)mB(pp:n1)三三實驗過程：實驗過程：1（1）編寫高斯消元法的）編寫高斯消元法的MATLAB文件

5、如下：文件如下：clearA=[0.1012.3043.5551.3473.7124.6232.8351.0725.643]b=[1.1832.1373.035][RARBnX]=gaus(Ab)運行結果為：運行結果為：請注意：因為RA=RB=n，所以此方程組有唯一解.RA=3RB=3n=3X=0.39820.01380.3351（2）編寫高斯消元法）編寫高斯消元法MATLAB文件如下：文件如下：clearA=[521283136]b=

溫馨提示

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

評論

0/150

提交評論