mgs.m

来自「流形学习程序」· M 代码 · 共 20 行

function [Q, R] = mgs(A);
% --- mgs function: Modified Gram-Schmidt
% Used by HLLE function.

[m, n] = size(A);   % Assume m>=n.
V = A;
R = zeros(n,n);
for i=1:n
    R(i,i) = norm(V(:,i));
    V(:,i) = V(:,i)/R(i,i);
    if (i < n)
        for j = i+1:n
            R(i,j) = V(:,i)' * V(:,j);
            V(:,j) = V(:,j) - R(i,j) * V(:,i);
        end;
     end;
 end;
 Q = V;