sortsignalestimates.m

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

% SORTSIGNALESTIMATES Compares estimated signals to originals,
% reverses permutation and calculates Signal to Noise Ratio
%
%  [sortedSigEst,SNR] = sortSignalEstimates(M,sigOriginal,sigEst)
%
%  inputs;
%  M : Number of signals;
%  sigOriginal : Original signals;
%  sigEst = M by T matrix of sources estimates
%
%  outputs;
%  sortedSigEst : Reordered sources estimates
%  SNR : Signal to noise ratio for original signals and their estimates
%
% Code available at http://www-bcl.cs.may.ie/~pogrady/
%
% Function provided as part of Soft-LOST matlab functions
%
% Paul D. O'Grady - 14-11-2003
% ([e-mail]) 
function [sortedSigEst,SNR] = sortSignalEstimates(M,sigOriginal,sigEst)



%May need to resize vectors
if (size(sigOriginal,2) ~= size(sigEst))
  
	xLimit = size(sigEst,2);	%fix scales

	%resize source signals to match FFT versions
	for count = 1:size(sigOriginal,1)
		sigOrig(count,:) = sigOriginal(count,1:xLimit);
	end;
else
  
	sigOrig = sigOriginal;
	
end;



for count1 = 1:M
  
	for count2 = 1:M
	  
	  [SNRTemp, pTemp] = snrPolarity(sigEst(count2,:)',sigOrig(count1,:)');
	  snrValues(count2) = SNRTemp;
	  pValues(count2) = pTemp;
	  
	end;
	
	[value, index] = max(snrValues);
	
	sortedSigEstTemp(count1,:) = pValues(index)*sigEst(index,:);
	SNRValueTemp(count1) =  value;
end;

sortedSigEst = sortedSigEstTemp ;
SNR = SNRValueTemp;