ldl_factor.m

来自「matlab 多个计算方程租的算法汇集 主要讲的是一般方法 还有正定矩阵的解法」· M 代码 · 共 23 行

function [L,D]=LDL_Factor(A)
n=length(A);
L=eye(n);
D=zeros(n);
d=zeros(1,n);
T=zeros(n);
for k=1:n
    d(k)=A(k,k);
    for j=1:k-1
        d(k)=d(k)-L(k,j)*T(k,j);
    end
    if abs(d(k))<1e-10
        return;
    end
    for i=k+1:n
    T(i,k)=A(i,k);
    for j=1:k-1
        T(i,k)=T(i,k)-T(i,j)*L(k,j);
    end
    L(i,k)=T(i,k)/d(k);
end
end
D=diag(d);