gram_schmidt.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function W=gram_schmidt(s1)

%2005.12.11

s1=s1/norm(s1);

W(:,1)=s1;
W(:,2:31)=randn(31,30);

for i=2:31
for j=1:i-1
q=W(:,i)'*W(:,j);
W(:,i)=W(:,i)-q*W(:,j);
end
W(:,i)=W(:,i)/norm(W(:,i));
end
%fprintf('Done ...\n');