LLT.m

% 平方根(LLT)分解
% 注释的地方需要修改
clear

% 以下A,b即原方程组化为矩阵形式Ax=b的A,b
A=[1,1,2;1,2,0;2,0,11];     % 系数矩阵A，逗号分隔行，分号分隔列
b=[5; 8; 7];                % 向量b

[n,m]=size(A);        
L=zeros(n,m);         
  
for i=1:n 
    L(i,i)=A(i,i);
    for k=1:i-1
        L(i,i)=L(i,i)-L(i,k)^2;
    end
    L(i,i)=L(i,i)^0.5;

    for j=i+1:n
        L(j,i)=A(j,i);
        for k=1:i-1
            L(j,i)=L(j,i)-L(j,k)*L(i,k);
        end
        L(j,i)=L(j,i)/L(i,i);
    end
end

y=L\b;
x=L'\y