chol_solve.m

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

function x=Chol_Solve(A,b)
%利用对称正定矩阵之Cholesky分解求解线性方程组Ax=b
n=length(b);
l=Cholesky(A);
x=ones(1,n);
y=ones(1,n);
for i=1:n
    z=0;
    for k=1:i-1
        z=z+l(i,k)*y(k);
    end
    y(i)=(b(i)-z)/l(i,i);
end
for i=n:-1:1
    z=0;
    for k=i+1:n
        z=z+l(k,i)*x(k);
    end
    x(i)=(y(i)-z)/l(i,i);
end