数字信号处理实验——用窗函数设计FIR滤波器.docVIP

实验四 用窗函数设计FIR滤波器 实验目的 熟悉FIR滤波器设计的基本方法。 掌握用窗函数设计FIR数字滤波器的原理及方法，熟悉相应的计算机高级语言编程。 熟悉线性相位FIR滤波器的幅频特性和相位特性。 了解各种不同窗函数对滤波器性能的响应。 实验原理和方法 窗函数法设计的任务在于寻找一个可实现有限长单位脉冲响应的传递函数H(ejw)=h(n)e-jwn 去逼近hd(n)=1/2Hd(ejw)ejwndw 即h(n)=hd(n)w（n） (一)几种常用的窗函数 1、矩形窗 w(n)=RN(n) 2、Hanning窗 w(n)=0.5[1-cos(2πn／N-1)]RN(n) 3、Hamming窗 w(n)=[0.54-0.46cos(2πn／N-1)]R N(n) 4、Blackman窗 w(n)=[0.42-0.5 cos(2πn／N-1)+0.08 cos(4πn／N-1)] RN(n) 5、Kaiser窗 w(n)=I0(β（1-[(2n／(N-1))-1]2）?)／I0(β） （二）窗函数法设计线性相位FIR滤波器的步骤 1、确定数字滤波器的性能要求。确定各临界频率{wk}和滤波器单位脉冲响应长度N。 2、根据性能要求和N值，合理地选择单位脉冲响应h(n)有奇偶对称性，从而确定理想频率响应hd(ejw)的幅频特性和相位特性。 3、用傅里叶反变换公式求得理想单位脉冲响应hd(n)。 4、选择适当的窗函数W（n），求得所设计的FIR滤波器单位脉冲响应。 5、用傅里叶变换求得其频率响应H (ejw)，分析它的幅频特性，若不满足要求，可适当改变窗函数形式或长度N，重复上述过程，直至得到满意的结果。 三、实验内容和步骤 1、分别用矩形窗、Hanning窗、Hamming窗、Blackman窗、Kaiser窗（β=8.5）设计一个长度N=8的线性相位FIR滤波器。 Window=boxcar(8); b=fir1(7,0.4,Window) b = Columns 1 through 7 -0.0883 0.0000 0.2061 0.3822 0.3822 0.2061 0.0000 Column 8 -0.0883 Window=hanning(8); b=fir1(7,0.4,Window) b = Columns 1 through 7 -0.0100 0.0000 0.1501 0.3599 0.3599 0.1501 0.0000 Column 8 -0.0100 Window=hamming(8); b=fir1(7,0.4,Window) b = Columns 1 through 6 -0.0072 0.0000 0.1351 0.3721 0.3721 0.1351 Columns 7 through 8 0.0000 -0.0072 Window=blackman(8); b=fir1(7,0.4,Window) b = Columns 1 through 6 0.0000 0.0000 0.1060 0.3940 0.3940 0.1060 Columns 7 through 8 0.0000 0.0000 Window=kaiser(8,8.5); b=fir1(7,0.4,Window) b = Columns 1 through 6 -0.0001 0.0000 0.1068 0.3933 0.3933 0.1068 Columns 7 through 8 0.0000 -0.0001 2、用Hanning窗设计一个线性相位带通滤波器 Window=hanning(16); b=fir1(15,[0.3 0.5],Window);freqz(b,1) Window=hanning(46); b=fir1(45,[0.3 0.5],Window);freqz(b,1) 分析实验结果：由图可看出当N由15变为45时，主瓣宽度和过渡带宽度都变窄了。 3、改用矩形窗和Blackman窗设计带通滤波器 Window=boxcar(16); b=fir1(15,[0.3 0.5],Window);freqz(b,1) Window=boxcar(46); b=fir1(45,[0.3 0.5],Window);freqz(b,1) Wind