这是国外一本经典教材附带的程序

				% 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')