llh.m

我自己在学习空时编码一书时编写的第一章的仿真程序

function y=llh(A)
n=length(A);
l(1,1)=sqrt(A(1,1));%calcul l_11
l(2:n,1)=A(2:n,1)/l(1,1);%calcul_n1
for i=2:n
    sum=0;
    for j=1:(i-1)
        sum=sum+l(i,j)*conj(l(i,j));
    end
    l(i,i)=sqrt(A(i,i)-sum);
    for j=(i+1):n
        sum=0;
        for k=1:(i-1)
            sum=sum+l(j,k)*conj(l(i,k));
        end
        l(j,i)=(A(j,i)-sum)/l(i,i);
    end
end
y=l;