数值分析matlab程序lu分解2页.doc

1、程序function L,U,zhibiao = LU_Decom( A )%求矩阵A的LU分解，其中，%A为要分解的矩阵%L为单位下三角阵%U为单位下三角阵%flag为指标向量，zhishi=failure表示计算失败，zhibiao=OK表示计算成功n,m=size(A);%要求所分解的矩阵是方阵；否则停止结算，并输出出错信息if n=m error(The rows and columns of matrix A must be equal!); return;endL=eye(n);U=zeros(n);zhibiao=OK;for k=1:n for j=k:n z=0; for q=1:k-1 z=z+L(k,q)*U(q,j); end U(k,j)=A(k,j)-z; end if abs(U(k,k)eps