平方根法求解对称正定线性方程组的Matlab实现

数值分析实验二：用平方根法和改进的平方根法求线性方程组

1.平方根法

function [] = Square_Root(A,b)
tic
%使用CHOLESKY分解的平方根法
%#ok<*NOPRT> 
n=length(b);%矩阵维度
A(1,1)=sqrt(A(1,1));%计算l(1,1)
for i=2:n
    A(i,1)=A(i,1)/A(1,1);%L的第一列元素
end
for k=2:n
    Sigma1=0;
    for j=1:k-1
        Sigma1=Sigma1+A(k,j)*A(k,j);%lkk减去的求和式子
    end
    A(k,k)=sqrt(A(k,k)-Sigma1);%计算lkk
    if k+1<=n %判断lkk是否为最后一个对角元素
        for i=k+1:n
            Sigma2=0;
            for r=1:k-1
                Sigma2=Sigma2+A(i,r)*A(k,r);%lik减去的求和式子
            end
            A(i,k)=(A(i,k)-Sigma2)/A(k,k);%计算lik
        end
    end
end
L=tril(A);
U=triu(A');
y=L\b;x=U\y;
for i=1:n %将解横向展示
    value=x(i,1);
    fprintf('x%d=%.4f; ',i,value)
end
fprintf('\n')%将光标移到