median1.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 30 行

function y = median1(s,M);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%	median1.m 
%
%	jmw 
%
%	8/16/93
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% median filter of variable length M (M odd)

L= (M-1)/2;

N = length(s);
tp = zeros(1,M);	% temporary vector
y = zeros(1,N);		% output

for i = (M-L):(N-L),
	for j = 0:L,
		tp(M-j-L) = s(i-j);
	end;
	for j = 1:L,
		tp(M+j-L) = s(i+j);
	end;
	y(i) = median(tp);
end;