% 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% (paul.ogrady@may.ie) function [sortedSigEst,SNR] = sortSignalEstimates(M,sigOriginal,sigEst)%May need to resize vectorsif (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;