function xmf = medfilt1mit(x,m,nit)
% --------------------------------------------------------------------------------------------
%  medfilt1mit.m: median filter 
%   median of the first/last m values are assumed to the left and right of v
%   nit - number of iteration of the filter
%   Version 1.00, Date: 07/04/2010
% --------------------------------------------------------------------------------------------
%   Maurizio Varanini, Clinical Physiology Institute, CNR, Pisa, Italy
%   For any comment or bug report, please send e-mail to: maurizio.varanini@ifc.cnr.it
% --------------------------------------------------------------------------------------------

if(nargin<3), nit=1; end
if(nargin<2), m=3; end
if(size(x,2)>size(x,1)), x=x'; colV=0; else colV=1; end
n=length(x);
m2=floor(m/2);
xi=median(x(1:min(n,m)));
xf=median(x(end-min(n,m)+1:end));
xt=[xi+zeros(m2,1); x; xf+zeros(m2,1)];
for it=1:nit
    xx=xt;
    xt=medfilt1(xx,m);
    if(all(~(xt-xx))), break, end
end
if(colV), xmf=xt(m2+1:end-m2); else xmf=xt(m2+1:end-m2)'; end
end %== function ================================================================
%