数值分析matlab程序lu分解

《数值分析matlab程序lu分解》由会员分享，可在线阅读，更多相关《数值分析matlab程序lu分解（2页珍藏版）》请在装配图网上搜索。

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 zhibiao=failure; return; end for i=k+1:n z=0; z=z+L(i,q)*U(q,k); end L(i,k)=(A(i,k)-z)/U(k,k);endend运算结果：A=1 2 3;4 5 6;7 8 0 L,U,zhibiao = LU_Decom( A )2、库函数A=1 2 3;4 5 6;7 8 0L,U,P=lu(A)运算结果：