mroot_squaring.m

来自「改进平方根法求解对称正定线性方程组」· M 代码 · 共 30 行

% 改进平方根法解对称系数阵线性方程组
function b=MRoot_squaring(A,b,n)
% 求解LR分解的各元素，并以重址方式存储在矩阵A中
for i=2:n
    for j=i:n
        S=A(i,j);
        for k=1:i-1
            S=S-A(k,i)*A(k,j)/A(k,k);
        end
        A(i,j)=S;
    end
end
% 解Ly=b  y也存放在b里
for i=2:n
    S=b(i,1);
    for j=1:i-1
        S=S-A(j,i)*b(j)/A(j,j);
    end
    b(i)=S;
end
%解Rx=y
b(n)=b(n)/A(n,n);
for i=n-1:-1:1
    S=b(i)/A(i,i);
    for j=i+1:n
        S=S-A(i,j)*b(j)/A(i,i);
    end
    b(i)=S;
end
return