求matlab LU 分解程序要写出详尽的程序及注释

网友回答

matlab有多种LU分解程序下面算一种：
function [L,U]=myLU(A)
%实现对矩阵A的LU分解,L为下三角矩阵
A [n,n]=size(A);
L=zeros(n,n);
U=zeros(n,n);
for i=1:n
L(i,i)=1;
end for k=1:n
for j=k:n
U(k,j)=A(k,j)-sum(L(k,1:k-1).*U(1:k-1,j)');
end for i=k+1:n
L(i,k)=(A(i,k)-sum(L(i,1:k-1).*U(1:k-1,k)'))/U(k,k);
end end 用法,在控制台输入
A=[1 2 3 -4;-3 -4 -12 13;2 10 0 -3;4 14 9 -13];