function [meaMi,meaMa] = meanMiMaSc(v,nel,percmi,percma)
% ---------------------------------------------------------------------------------------------
% [meaMi,meaMa] = meanMiMaSc(v,nel,percmi,percma)
% ---------------------------------------------------------------------------------------------
% meanmaxsc.m: Compute the average value of the minima and of the maxima
% computed on data intervals of the input vector.
% Distribution tails can be excluded using parameters "perci" and "percf"
% "v" = input data vector
% "nel" = data interval length
% "percmi" = % of min values (mins and maxs) to be excluded; (if positive)
% number of min values (if negative)
% "percma" = % of max values to be excluded (if positive)
% number of max (mins and maxs) values (if negative)
%
% Version 1.00, Date: 01/05/2000
% ---------------------------------------------------------------------------------------------
% 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), percmi=5; end
if(nargin<4), percma=percmi; end
if(nel>length(v)), nel=length(v); end
j=1;
mini=zeros(floor(length(v)/nel),1);
maxi=zeros(floor(length(v)/nel),1);
for i=1:nel:length(v)-nel+1
mini(j)=min(v(i:i+nel-1));
maxi(j)=max(v(i:i+nel-1));
j=j+1;
end
if(percmi<0), ii=1-percmi;
else ii=1+floor(length(maxi)*percmi/100); end
if(percma<0), fi=length(maxi)+percma;
else fi=length(maxi)-floor(length(maxi)*percma/100); end
omaxi=sort(maxi);
omini=sort(mini);
meaMi=mean(omini(ii:fi));
meaMa=mean(omaxi(ii:fi));
end %== function ================================================================
%