用该窗函数设计的fir数字滤波器.pptVIP

同样，根据公式（7.2.1）： 将Hd(ejω)写成Hd(ejω)=Hdg(ω)e－jω?， 理想低通滤波器的幅度特性函数（如图7.2.3（a）所示）为 7.2.3 用窗函数法设计FIR滤波器的步骤 用窗函数法设计FIR滤波器的步骤如下： 根据对过渡带及阻带衰减的指标要求，选择窗函数的 类型，并估计窗口长度N。 先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄(过渡带宽度Bt近似等于窗函数主瓣宽度) 的窗函数。 然后规定根据选用的窗函数的过渡带宽度（精确值或近似值）小于或等于实际的过渡带宽度，求出窗口长度N。 (2) 构造希望逼近的理想频率响应函数Hd(ejω)，即   对所谓的“标准窗函数法”，就是选择Hd(ejω)为线性 相位理想滤波器（理想低通、理想高通、理想带通、 理想带阻）。以低通滤波器为例，Hdg(ω)应满足：  　　　　　　　　　　　　　　　　 （7.2.19） 如果设计指标给定通带边界频率和阻带边界频率 ωp和ωs， 一般取 （7.2.20） (3) 计算理想的单位脉冲响应hd(n)。 如果给出待求滤波器的频响函数为Hd(ejω)，那么单位 脉冲响应用下式求出： 如果Hd(ejω)较复杂，或者不能用封闭公式表示，则不 能用上式求出hd(n)。我们可以对Hd(ejω)从ω=0到2π 采样M点，采样值为 然后对上式进行M点IDFT(IFFT)，得到： 　　　　　　 （7.2.21） k=0 ,1, 2, …, M－1 （7.2.22） 根据频域采样理论，hdM(n)与hd(n)应满足如下关系： 因此，如果M选得较大，可以保证在窗口内hdM(n)有效 逼近hd(n)。 对（7.2.19）式给出的线性相位理想低通滤波器作为 Hd(ejω)，由（7.2.21）式求出单位脉冲响应hd(n)： 为保证线性相位特性， ? =(N－1)/2。 (4) 加窗得到设计最终结果：h(n)=hd(n)w(n)。 【例7.2.1】 用窗函数法设计线性相位高通FIRDF，要求通 带截止频率ωp=π/2 rad，阻带截止频率ωs=π/4 rad，通带 最大衰减? p=1 dB，阻带最小衰减? s=40 dB。 解 ： （1）选择窗函数w(n)，计算窗函数长度N。已知阻带最小 衰减? s=40 dB，由表（7.2.2）可知汉宁窗和哈明窗均满足 要求，我们选择汉宁窗。 本例中窗函数过渡带宽度Bt≤ωp－ωs=π/4， 汉宁窗的精 确过渡带宽度Bt=6.2π/N，所以要求Bt=6.2π/N≤π/4，解 得N≥24.8。对高通滤波器N必须取奇数，取N=25。由式 （7.2.11）, 有 （2）构造理想滤波器的频响特性Hd(ejω):  式中 （3）求出理想滤波器的单位脉冲响应hd(n):  将τ=12代入得  δ(n－12)对应全通滤波器，　　　　　　是截止频率为3π/8的理想低通滤波器的单位脉冲响应，二者之差就是理想高通滤波器的单位脉冲响应。这就是求理想高通滤波器的单位脉冲响应的另一个公式。 （4） 加窗: 图7.2.2 吉普斯效应 下面讨论这种截断效应的产生，以及如何构造窗函数w(n), 用来减少截断效应，设计一个能满足技术要求的FIR线性相位滤波器。 以上就是用窗函数法设计FIR滤波器的思想。 另外，我们知道Hd(ejω)是一个以2π为周期的函数，可 以展为傅里叶级数，即   傅里叶级数的系数为hd(n)，当然就是Hd(ejω)对应的单位 脉冲响应。 设计FIR滤波器就是根据要求找到N个傅里叶级数系数 h(n)，n=1, 2, … , N－1，以N项傅氏级数去近似代替无 限项傅氏级数，这样在一些频率不连续点附近会引起较 大误差，这种误差就是前面说的截断效应，如图7.2.2所 示。 因此，窗函数法也称为傅氏级数法。显然，选取傅氏级 数的项数愈多，引起的误差就愈小，但项数增多即h(n) 长度增加，也使成本和滤波计算量加大，所以应在满足 技术要求的条件下，尽量减小h(n)的长度。　　 在（7.2.3）式中，RN(n) 就是起到对无限长序列的截断作 用，可以形象地把RN(n)看做一个窗口，h(n)则是从窗口 看到的一段hd(n)序列，所以称h(n)=hd(n)RN(n)为用矩形 窗对hd(n)进行加窗处理。 下面分析用矩形窗截断的