tomwalters@0: %
tomwalters@0: %
tomwalters@0: %	GetF0median(F0sequence,NumMed)
tomwalters@0: %
tomwalters@0: function F0median = GetF0median(F0sequence,NumMed)
tomwalters@0: 
tomwalters@0: if nargin < 2, NumMed = []; end;
tomwalters@0: if length(NumMed) == 0,NumMed = 5; end; 
tomwalters@0: 
tomwalters@0: [NumCh, LenF0] = size(F0sequence);
tomwalters@0: 
tomwalters@0: Mtx = zeros(NumMed,LenF0+NumMed);
tomwalters@0: 
tomwalters@0: for nm = 1:NumMed
tomwalters@0:   nnl = nm-1+(1:LenF0);
tomwalters@0:   Mtx(nm,nnl) = F0sequence(:)';
tomwalters@0: end;
tomwalters@0: 
tomwalters@0: aa = median(Mtx,1);
tomwalters@0: F0median = aa(fix(NumMed/2)+(1:LenF0));
tomwalters@0: 
tomwalters@0: 
tomwalters@0: