數(shù)值線性代數(shù)第二版徐樹方高立張平文上機習(xí)題第一章實驗報告

1、上機習(xí)題1.先用你所熟悉的的計算機語言將不選主元和列主元Gauss消去法編寫成通用的子程序；然后用你編寫的程序求解84階方程組；最后將你的計算結(jié)果與方程的精確解進(jìn)行比較，并就此談?wù)勀銓auss消去法的看法。Sol：（1）先用matlab將不選主元和列主元Gauss消去法編寫成通用的子程序，得到：PUL不選主元Gauss消去法：得到滿足??)(AGaussLAUL?ULLUA?列主元Gauss消去法：得到滿足??)(AGaussColP

2、UL?PULLUPA?（2）用前代法解，得??PbbLy?y用回代法解，得yUx?x求解程序為（可缺省，缺省時默認(rèn)為單位矩陣）??PULbAGaussx?P（3）計算腳本為ex1_1代碼代碼%算法算法1.1.31.1.3（計算三角分解：（計算三角分解：GaussGauss消去法）消去法）functionfunction[LU]=GaussLA(A)[LU]=GaussLA(A)n=length(A)n=length(A)ffk=1:n1

3、k=1:n1A(k1:nk)=A(k1:nk)A(kk)A(k1:nk)=A(k1:nk)A(kk)A(k1:nk1:n)=A(k1:nk1:n)A(k1:nk)A(kk1:n)A(k1:nk1:n)=A(k1:nk1:n)A(k1:nk)A(kk1:n)endendU=triu(A)U=triu(A)ffi=1:n1i=1:n1temp=P(i:)temp=P(i:)P(i:)=P(u(i):)P(i:)=P(u(i):)P(u(i)

4、:)=tempP(u(i):)=tempendendendend%高斯消去法解線性方程組高斯消去法解線性方程組functionfunctionx=Gauss(AbLUP)x=Gauss(AbLUP)ififnargin5nargin5P=eye(length(A))P=eye(length(A))endendn=length(A)n=length(A)b=Pbb=Pbffj=1:n1j=1:n1b(j)=b(j)L(jj)b(j)=b(