高斯消去法求矩阵的逆.doc

《高斯消去法求矩阵的逆.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高斯消去法求矩阵的逆。a=input('请输入线性方程组的系数矩阵');b=input('请输入线性方程组的右端列向量');[Row,Col]=size(a);ifRow~=Col'线性方程组的系数矩阵非方阵,程序终止'returnendifdet(a)==0'线性方程组的系数矩阵行列式为0'returnend%Gauss消元法fork=1:n-1ifa(k,k)==0'对角线元素a(%1d,%1d)为零，程序终止'returnendfori=k+1:nl(i,k)=a(i,k)/a(k,k);forj=k+1:na(i,

2、j)=a(i,j)-l(i,j)*a(k,j);endb(i)=b(i)-l(i,k)*b(k);endend%回代过程ifa(n,n)==0'对角线元素a（%1d,%1d)为零,程序终止'returnendx(n)=b(n)/a(n,n);'线性方程组解为'x(n)fori=n-1:-1:1x(i)=b(i);forj=k+1:nx(i)=x(i)-a(i,j)*x(j);endx(i)=x(i)/a(i,i);x(i)end