diff Code/Descriptors/Matlab/MPEG7/FromWeb/VoiceSauce/func_pickpeaks.m @ 4:92ca03a8fa99 tip

Update to ICASSP 2013 benchmark

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Code/Descriptors/Matlab/MPEG7/FromWeb/VoiceSauce/func_pickpeaks.m	Wed Feb 13 11:02:39 2013 +0000
@@ -0,0 +1,49 @@
+function [ymax,inx] = func_pickpeaks(y, winlen)
+% [ymax,inx] = func_pickpeaks(y, winlen)
+% Input:  y - from wavread
+%         winlen - size of window to use
+% Output: ymax - peak values
+%         inx - peak positions
+% Notes:  
+%
+% Author: Yen-Liang Shue and Markus Iseli, Speech Processing and Auditory Perception Laboratory, UCLA
+% Copyright UCLA SPAPL 2009
+
+% first, get all maxima
+ymin = min(y);
+y = y - ymin;
+dy = [y(1); transpose(diff(y))];   % diff() gives a shift by 1: insert one sample
+inx1 = diff(dy < 0);    % is +1 for maxima, -1 for minima
+inx1(1) = 0; % do not accept maxima at begining
+inx1(end) = 0; % do not accept maxima at end
+inx = inx1 > 0;         % choose maxima only
+ymax= y(inx);
+inx = find(inx);
+%plot(y);
+%hold on;
+%plot(inx,ymax,'or');
+%hold off;
+
+nofmax = length(ymax);
+if nofmax==1
+   return;
+end
+% now filter maxima with window of length winlen
+for cnt = 1 : nofmax
+   arr = inx(1:cnt);
+   cmp = inx(cnt)-winlen;
+   arr2 = arr>cmp;
+   %ymax(arr2)
+   [m, mi] = max(ymax(arr2));
+   ymax(arr2)=-60000;
+   ymax(mi+length(arr2)-sum(arr2))=m;
+   %ymax(arr2)
+end
+temp = find(ymax>0);
+inx  = inx(temp);
+ymax = ymax(temp);
+%plot(y);
+%hold on;
+%plot(inx,ymax,'or');
+%hold off;
+ymax = ymax + ymin;