grmsch.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 13 行

function [Q,R]=GrmSch(A)
[m,n]=size(A);
for j=1:n
    V(:,j)=A(:,j);
    for i=1:j-1
        R(i,j)=Q(:,i)'*A(:,j);
        V(:,j)=V(:,j)-R(i,j)*Q(:,i);
    end
    R(j,j)=V(:,j)'*V(:,j);
    R(j,j)=sqrt(R(j,j));
    Q(:,j)=V(:,j)/R(j,j);
end