数字滤波器设计

.04210.01630.00540.00150.00030.00010.00000.0000B=1.0e+007*Columns1through18000000000000000000Columns19through36000000000000000000Columns37through4400000003.1532采用窗函数法(分别用汉宁窗、哈明窗、布莱克曼窗函数)设计满足要求的FIR低通滤波器,求出h(n),并画出幅频响应损耗函数曲线.汉宁窗程序Fs=20000;fp=4000;fs=4500;m=[1100];wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Rp=0.5;As=50;Bt=ws-wp;N0=ceil(6.6*pi/Bt);N=N0+mod(N0+1,2);wc=(wp+ws)/2/pi;hn=fir1(N-1,wc,hanning(N));freqz(hn,1,512);plot(w,20*log(abs(hn)));grid;axis([0,1,-1000,100]);xlabel('频率/kHz');ylabel('幅值/dB');滤波器h(n)的波形及其幅频响应损哈明窗程序Fs=20000;fp=4000;fs=4500;m=[1100];wp=2*pi*fp/Fs;ws=2*pi*fs/Fs;Rp=0.5;Rs=50;Bt=ws-wp;N0=ceil(6.6*pi/Bt);N=N0+mod(N0+1,2);wc=(wp+ws)/2/pi;hn=fir1(N-1,wc,hamming(N));freqz(hn,1,512);plot(w,20*log(abs(hn)));gridon;axis([0,1,-1000,100]);xlabel('频率/kHz');ylabel('幅值/dB');滤波器h(n)的波形及其幅频响应损耗函数曲线