MATLAB Cholesky 分解算法实现及代码优化

B = [2,1,2;4,3,1;6,1,5]; % 原始矩阵 n = size(B,1); % 矩阵维度 L = zeros(n); % 初始化下三角矩阵

for j = 1:n for i = j:n if i == j L(i,j) = sqrt(B(i,j) - sum(L(i,1:j-1).^2)); else L(i,j) = (B(i,j) - L(i,1:j-1)*L(j,1:j-1)') / L(j,j); end end end

L % 输出下三角矩阵