ldl_solve.m

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

function x=LDL_Solve(A,b)
%利用对对称正定矩阵A进行LDL'分解法求解线性方程组Ax=b
n=length(b);
[l,d]=LDL_Factor(A);

y(1)=b(1);
for i=2:n
    z=0;
    for k=1:i-1
        z=z+l(i,k)*y(k);
    end
    y(i)=b(i)-z;
end

x(n)=y(n)/d(n,n);
for i=n-1:-1:1
    z=0;
    for k=i+1:n
        z=z+l(k,i)*x(k);
    end
    x(i)=y(i)/d(i,i)-z;
end