《基于matlab的数值分析课程设计》.docVIP

数值分析试验报告 矩阵的LU分解 题目：求4阶矩阵的LU分解 方法：杜里特尔分解法 程序： function f=LU_decom(A) [m,n]=size(A) L=eye(n); U=zeros(n); flag=ok; for i=1:n U(1,i)=A(1,i); end for r=2:n L(r,1)=A(r,1)/U(1,1); end for i=2:n for j=i:n z=0; for r=1:i-1 z=z+L(i,r)*U(r,j); end U(i,j)=A(i,j)-z; end if abs(U(i,i))eps flag=failure return; end for k=i+1:n m=0; for q=1:i-1 m=m+L(k,q)*U(q,i); end L(k,i)=(A(k,i)-m)/U(i,i); end end L U end 结果 LU_decom(A) m = 4 n = 4 L