实验四窗函数法设计FIR数字滤波器论述.docxVIP

数字信号处理实验报告  PAGE \* MERGEFORMAT 16 实验四 窗函数法设计FIR数字滤波器 一、实验目的 1、掌握窗函数法设计FIR数字滤波器的原理及具体方法。 2、掌握频率取样法设计FIR数字滤波器的原理和基本方法。 3、学习利用窗函数法和频率取样法设计低通、带通、高通、带阻数字滤波器。 二、实验环境 计算机、MATLAB软件 三、实验基础理论 窗函数设计FIR滤波器 1.基本原理 窗函数设计法的基本思想为，首先选择一个适当的理想的滤波器，然后用窗函数截取它的单位脉冲响应，得到线性相位和因果的FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器，使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。 2.设计步骤 （1）给定理想滤波器的频率响应，在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽为的低通滤波器由下式给??： 其中为采样延迟，其作用是为了得到一个因果系统。 （2）确定这个滤波器的单位脉冲响应 为了得到一个长度为N的因果的线性相位FIR滤波器，我们令 （3）用窗函数截取得到所设计FIR数字滤波器： 3.窗函数的选择 常用的窗函数有矩形（Rectangular）窗，汉宁（Hanning）窗，海明（Hamming）窗、布莱克曼（Blackman）窗、凯瑟（Kaiser）窗等 表4-1 MATLAB中产生窗函数的命令 MATLAB函数窗函数MATLAB函数窗函数Boxcar矩形窗函数Blackman布莱克曼窗Hanning汉宁窗函数Kaiser凯瑟窗函数Hamming海明窗 表4-2 常用窗函数的特性 窗函数窗函数频率特性加窗后滤波器指标旁瓣峰值dB主瓣宽度过渡带宽最小阻带衰减dB矩形窗-134π/N1.8π/N-21汉宁窗-318π/N6.2π/N-44海明窗-418π/N6.6π/N-53布莱克曼窗-5712π/N11π/N-74凯瑟窗是一种广泛在实际中广泛应用的窗函数，它由下式给定： 其中是修正的零阶贝塞尔函数，参数控制最小阻带衰减，这种窗函数对于相同的N可以提供不同的过渡带宽。由于贝塞尔函数比较复杂，这种窗函数的设计方程很难推导，然而幸运的是，有一些经验设计方程可以直接使用。 已知给定的指标 ，滤波器长度N和凯瑟窗参数可以按如下凯瑟窗方程给出 过渡带带宽： 频率取样设计FIR滤波器 1.基本原理 频率取样法从频域出发，把理想的滤波器等间隔采样得到,将作为实际设计滤波器的： 得到以后可以由来确定唯一确定滤波器的单位脉冲响应,可以由求得： 其中为内插函数： 有求得的频率响应将逼近。 如果我们设计的是线性相位FIR滤波器，则的幅度和相位满足线性相位滤波器的约束条件。 我们将表示为如下形式 当为实数，则 由此得到 即为中心偶对称。在利用线性相位条件可知，对于1型和2型线性相位滤波器： 对于3型和4型线性相位滤波器 2.设计步骤 （1）由给定的理想滤波器给出和。 （2）由求得 （3）根据求得或 四、实验内容 1、设计一个数字低通FIR滤波器，其技术指标如下： 分别采用矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。结合实验结果，分别讨论采用上述方法设计的数字滤波器是否都能满足给定指标要求。 矩形窗 程序代码： wp=0.2*pi;wst=0.3*pi;tr_width=wst-wp; N=ceil(1.8*pi/tr_width) n=0:N-1; wc=(wst+wp)/2; alpha=(N-1)/2; hd=(wc/pi)*sinc((wc/pi)*(n-alpha)); w_boxcar=boxcar(N); h=hd.*w_boxcar; subplot(221); stem(n,hd,filled); axis tight;xlabel(n);ylabel(hd(n)); [Hr,w1]=zerophase(h); subplot(222); plot(w1/pi,Hr); axis;xlabel(\omega/\pi);ylabel(H(\omega)); subplot(223); stem(n,h,filled); axis tight;xlabel(n);ylabel(h(n)); [H,w]=freqz(h,1); subplot(224); plot(w/pi,20*log10(abs(H)/max(H))); axis tight;xlabel(\omega/\pi);ylabel(dB); grid on; MATLAB图形： （2）汉宁窗 程序代码： wp=0.2*pi;