数值分析matlab程序lu分解

1、程序function[L,U,zhibiao]=LU_Decom(A)%求矩阵A的LU分解，其中，%A为要分解的矩阵%L为单位下三角阵%U为单位下三角阵%flag为指标向量，zhishi='failure'表示计算失败，zhibiao='OK'表示计算成功[n,m]=size(A);%要求所分解的矩阵是方阵；否则停止结算，并输出出错信息ifn~=merror('TherowsandcolumnsofmatrixAmustbeequal!');return;endL=eye(n);U=zeros(n);zhibiao='OK';fork=1:nforj=k:nz=0;forq=1:k-1z=z+L(k,q)*U(q,j);endU(k,j)=A(k,j)-z;endifabs(U(k,k))epszhibiao='failure';return;endfori=k+1:nz=0;z=z+L(i,q)*U(q,k);endL(i,k)=(A(i,k)-z)/U(k,k);endend运算结果：A=[123;456;780][L,U,zhibiao]=LU_Decom(A)OKzhibiaoUL130065032112.170100012、库函数A=[123;456;780][L,U,P]=lu(A)运算结果：0100011005.40038571.0008715.05714.0011429.0001PUL