高斯消元法matlab實現(xiàn)

1、《數(shù)值分析》實驗報告《數(shù)值分析》實驗報告一、實驗目的與要求一、實驗目的與要求1掌握高斯消去法的基本思路和迭代步驟；2培養(yǎng)編程與上機調(diào)試能力。二、實驗內(nèi)容二、實驗內(nèi)容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=