program_2_5.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

% Program 2_5
% Illustration of Median Filtering
%
N = input('Median Filter Length = ');
R = 50; a = rand(1,R)-0.4;
b = round(a); % Generate impulse noise
m = 0:R-1;
s = 2*m.*(0.9.^m); % Generate signal
x = s + b; % Impulse noise corrupted signal
y = medfilt1(x,3); % Median filtering
subplot(2,1,1)
stem(m,x);axis([0 50 -1 8]);
xlabel('n');ylabel('Amplitude');
title('Impulse Noise Corrupted Signal');
subplot(2,1,2)
stem(m,y);
xlabel('n');ylabel('Amplitude');
title('Output of Median Filter');