计算方法_矩阵LU分解法.doc

《计算方法_矩阵LU分解法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、clearall;%A=LU矩阵三角分解法n=input('输入方矩阵的维数:');fori=1:nforj=1:nA(i,j)=input('依次输入矩阵元素：');endend%输入一个n阶方形矩阵forj=1:nL(j,j)=1;%Doolittle分解，L对角元素全为1endforj=1:nU(1,j)=A(1,j);end%U的第一行fori=2:nL(i,1)=A(i,1)/U(1,1);end%L的第一列fork=2:nforj=k:nsum1=0;form=1:k-1sum1=sum1+L

2、(k,m)*U(m,j);end%求和U(k,j)=A(k,j)-sum1;endfori=k+1:nsum2=0;form=1:k-1sum2=sum2+L(i,m)*U(m,k);end%求和L(i,k)=(A(i,k)-sum2)/U(k,k);endendL%输出下三角矩阵LU%输出上三角矩阵U运行结果：（示例）输入方矩阵的维数:4依次输入矩阵元素：1依次输入矩阵元素：1依次输入矩阵元素：2依次输入矩阵元素：3依次输入矩阵元素：0依次输入矩阵元素：2依次输入矩阵元素：1依次输入矩阵元素：2依次输入

3、矩阵元素：1依次输入矩阵元素：-1依次输入矩阵元素：2依次输入矩阵元素：2依次输入矩阵元素：2依次输入矩阵元素：2依次输入矩阵元素：5依次输入矩阵元素：9A=LU分解后则可以求解Ax=b线性方程组，相关计算参考计算方法，这里不再详细介绍。