snrpolarity.m

来自「语音信号处理前分离信号噪声,matlab环境下运行,使信号处理效果更好」· M 代码 · 共 43 行

% SNRPOLARITY Signal to Noise ratio
%
%  [a,p] = snr(est,org)
%
%  est =  the estimate vector
%  org = the original vector
%  a = the SNR in dB
%  p = the est polarity with minimum snr 
%
%  tests for est and inverted est, function returns largest SNR
%
%  20Log10(Orig/Err)
%
% Paul D. O'Grady - 14-11-2003
% ([e-mail])

function [a,p] = snrPolarity(est,org)

orgNorm = org/norm(org);
estNorm = est/norm(est);

SNR(1) = norm(orgNorm)/norm(estNorm-orgNorm);

estNorm = (est*-1/norm(est*-1));
SNR(2) = norm(orgNorm)/norm(estNorm-orgNorm);

[SNRResult polarity] = max(SNR);


%return results
a = 20*log10(SNRResult);
if (polarity == 2)
p = -1;
else
p = 1;
end