N=51；

clear

N=51;wc=0.4;

h=fir1(50,wc,boxcar(N)); subplot(211); stem(h);

xlabel('n');ylabel('h'); [H,W]=freqz(h,1) subplot(212); plot(W/pi,abs(H)); title('矩形窗振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')

0.60.40.20-0.2010203040n矩形窗振幅特性/dB5060hH(w)1.510.5000.10.20.30.40.50.6相对频率0.70.80.91

N=76

clear

N=76;wc=0.4;

h=fir1(75,wc,boxcar(N));

subplot(211); stem(h);

xlabel('n');ylabel('h'); [H,W]=freqz(h,1) subplot(212); plot(W/pi,abs(H)); title('矩形窗振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')

0.60.40.20-0.2010204050n矩形窗振幅特性/dB30607080hH(w)1.510.5000.10.20.30.40.50.6相对频率0.70.80.91

N=26

clear

N=26;wc=0.4;

h=fir1(25,wc,hanning(N)); subplot(211); stem(h);

xlabel('n');ylabel('h'); [H,W]=freqz(h,1)

subplot(212); plot(W/pi,abs(H)); title('汉宁窗振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')

0.60.40.20-0.205101520n汉宁窗振幅特性/dB2530hH(w)1.510.5000.10.20.30.40.50.6相对频率0.70.80.91

N=51

clear

N=51;wc=0.4;

h=fir1(50,wc,hanning(N)); subplot(211); stem(h);

xlabel('n');ylabel('h'); [H,W]=freqz(h,1) subplot(212); plot(W/pi,abs(H)); title('汉宁窗振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')

0.60.40.20-0.2010203040n汉宁窗振幅特性/dB5060hH(w)1.510.5000.10.20.30.40.50.6相对频率0.70.80.91

N=76

clear

N=76;wc=0.4;

h=fir1(75,wc,hanning(N)); subplot(211); stem(h);

xlabel('n');ylabel('h'); [H,W]=freqz(h,1) subplot(212); plot(W/pi,abs(H)); title('汉宁窗振幅特性/dB'); xlabel('相对频率');ylabel('H(w)')