clear all                  %滑動平均濾波器
n=0:1:20;                 
sn=2.*(n.*0.9.^n);%signal
subplot(311)
stem(n,sn)
axis([-5 25 -0.5 10])
dn=randn(1,21);%噪聲
xn=sn+dn;
subplot(312)
stem(n,xn)
axis([-5 25 -0.5 10])
m=5;
a=1;
b=ones(1,m);
tn=filter(b,a,xn);
yn=1/m.*tn;
subplot(313)
stem(n,yn)
axis([-5 25 -0.5 10])