乔列斯基分解

function L=cholesky_decomposition(A)
    %{
        乔列斯基分解
        param A：实对称正定矩阵
        return L：分解结果，其中A=L*L^T
    %}
    %可以检验是否为实对称正定矩阵

    rows=size(A,1);
    cols=size(A,2);
    L=zeros(rows,cols);

    for jj=1:cols
        %首先计算L主对角线上的元素
        temp=0;
        for kk=1:jj-1
            temp=temp+L(jj,kk)^2;
        end

        L(jj,jj)=sqrt(A(jj,jj)-temp);

        %再计算这一列剩下的元素
        temp=0;
        for kk=1:jj-1
            temp=temp+L(jj,kk)*L(jj:end,kk);
        end
        L(jj:end,jj)=(A(jj:end,jj)-temp)/L(jj,jj);

        %TODO:还可以优化
    end
end