施密特正交化(matlab代码)
function basis = Gram_Schimidt( rank_n , basis )
%作者:192152王旭
%功能:利用Granm_Schimidt算法解正交basis
%输入要求:rank_n是基底的阶数,n维空间有n个向量作为基底
%输入为 rank_n=3 , basis =[ 1 0 0 ; 4 2 15 ; 0 0 3 ]
%输出为 basis = [ 1 0 0 ; 0 2 15 ; 0 -90/229 12/229]
%--------------------------------------------------------------------------
%测试需要
%输入为 rank_n=3 , basis =[ 1 0 0 ; 4 2 15 ; 0 0 3 ]
%basis = Granm_Schimidt( rank_n , basis )
%输出为 basis = [ 1 0 0 ; 0 2 15 ; 0 -90/229 12/229]
%--------------------------------------------------------------------------
for i = 2:rank_n
for j = 1:i-1
basis(i,:) = basis(i,:) - dot( basis(i,:) , basis(j,:) ) / dot( basis(j,:) , basis(j,:) ) * ( basis(j,:) )
end
end
basis
end