满秩方阵的LDU分解

%本函数将一个满秩方阵按LDU方式分解
%本函数需要调用 Doolittle分解

%请确保该文件存在工作目录下
function [L,D,U]=LDU(A)
    [L U]=Doolittle(A);
    b=size(A);
    n=b(1);
    D=zeros(n,n);
   
    for i=1:n
        D(i,i)=U(i,i);
    end
    U=inv(D)*U;
end

 

或者把DooLittle函数整合到这个函数中来。