GS 正交化算法

GS.m

% 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

>> A=GS([1 0 1;2 2 0;3 1 1])
 
A =
 
    1.0000   -0.5000    0.3810
    2.0000    1.0000    0.0952
    3.0000   -0.5000   -0.1905