Gram_

施密特正交化（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