program_2_4.m

来自「这是国外一本经典教材附带的程序」· M 代码 · 共 20 行

% Program 2_4
% Signal Smoothing by a moving-average Filter
%
R = 50;
d = rand(R,1)-0.5;
m = 0:1:R-1;
s = 2*m.*(0.9.^m);
x  = s + d ;
plot(m,d,'r-',m,s,'b--',m,x,'g:')
xlabel('Time index n'); ylabel('Amplitude')
legend('r-','d[n]','b--','s[n]','g:','x[n]');
pause
M = input('Number of input samples = ');
b = ones(M,1)/M;
y = filter(b,1,x);
plot(m,s,'r-',m,y,'b--')
legend('r-','s[n]','b--','y[n]');
xlabel ('Time index n');ylabel('Amplitude')