cholesky.m

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

function L=Cholesky(A)
%对对称正定矩阵A进行Cholesky分解
n=length(A);
L=zeros(n);
for k=1:n
    delta=A(k,k);
    for j=1:k-1
        delta=delta-L(k,j)^2;
    end
    if delta<1e-10
        return;
    end
    L(k,k)=sqrt(delta);
    for i=k+1:n
        L(i,k)=A(i,k);
        for j=1:k-1
            L(i,k)=L(i,k)-L(i,j)*L(k,j);
        end
        L(i,k)=L(i,k)/L(k,k);
    end
end