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
No comments:
Post a Comment