% G-S正交化
% A = GS(V)
% V 為n*n向量
% A 為n*n正交向量
% Example:
%   A=GS([1 0 1;2 2 0;3 1 1])
% Author：yujunjie
% Author's email: ustb03-07@yahoo.com.cn
function A = GS(V)
n=length(V);
for i=1:n
    k=0;  
    for j=1:i-1
        k=k+(A(:,j)'*V(:,i)/(A(:,j)'*A(:,j)))*A(:,j);
    end
    A(:,i)=V(:,i)-k;
    end