已经有专门的函数了 B=chol(A)上三角矩阵B和正定矩阵A满足关系A=B'B 所以你要的L就是B'———我自己编了一个，你拿去用好了 function L=Cholesky(A)n=length(A);L=zeros(n);for j=1:length(A)L(j,j)=sqrt(A(j,j)-sum(L(j,1:j-1).^2));for i=(j+1):n L(i,j)=(A...