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

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

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

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

4、(B)zhica=RBRAifzhica0disp(請(qǐng)注意：因?yàn)镽A~=RB，所以此方程組無(wú)解.)returnendifRA==RBifRA==ndisp(請(qǐng)注意：因?yàn)镽A=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)三三實(shí)驗(yàn)過(guò)程：實(shí)驗(yàn)過(guò)程：1（1）編寫高斯消元法的）編寫高斯消元法的MATLAB文件

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