xs.m

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

function r=xs(s,y)
%author lucky zhang
%copyright 2005.6.23
R=y*s';
Y=sum((y.^2)');
S=sum((s.^2)');
r=R./sqrt(Y'*S);
r=abs(r);