snr.m

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

function r=snr(s,y)
%author lucky zhang
%copyright 2005.6.29
%信号噪声比函数
ri=s*y'/(s*s');
r=ri^2*s*s'/((ri*s-y)*(ri*s-y)');
r=10*log10(r);