Mercurial > hg > aimmat

diff aim-mat/modules/usermodule/mellin/MIpack/GetF0median.m @ 0:74dedb26614d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Initial checkin of AIM-MAT version 1.5 (6.4.2011).
author tomwalters
date Fri, 20 May 2011 12:32:31 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/aim-mat/modules/usermodule/mellin/MIpack/GetF0median.m	Fri May 20 12:32:31 2011 +0100
@@ -0,0 +1,23 @@
+%
+%
+%	GetF0median(F0sequence,NumMed)
+%
+function F0median = GetF0median(F0sequence,NumMed)
+
+if nargin < 2, NumMed = []; end;
+if length(NumMed) == 0,NumMed = 5; end; 
+
+[NumCh, LenF0] = size(F0sequence);
+
+Mtx = zeros(NumMed,LenF0+NumMed);
+
+for nm = 1:NumMed
+  nnl = nm-1+(1:LenF0);
+  Mtx(nm,nnl) = F0sequence(:)';
+end;
+
+aa = median(Mtx,1);
+F0median = aa(fix(NumMed/2)+(1:LenF0));
+
+
+