mgrmsch.m

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

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