fir数字滤波器 实验四窗函数法设计FIR数字滤波器.doc

fir数字滤波器 实验四窗函数法设计FIR数字滤波器 实验四 窗函数法设计FIR数字滤波器 一、实验目的 1、掌握窗函数法设计FIR数字滤波器的原理及具体方法。 2、掌握频率取样法设计FIR数字滤波器的原理和基本方法。 3、学习利用窗函数法和频率取样法设计低通、带通、高通、带阻数字滤波器。 二、实验环境… 一、HRM实践教学论坛日程 二、HRM沙盘模拟大赛日程 三、主讲嘉宾介绍 杨河清教授 日本流通经济大学经济学博士，国 家级教学名师、国务院特贴专家、中国 杰出社会科学家。1999年至2012年任首 都经济贸易大学劳动经济学院第一任院 长，现任首… 浅谈混凝土施工缝的留置及施工中的注意事项 作者：马甲子 刘芳 来源：《房地产导刊》2014年第05期 【摘 要】混凝土施工缝问题在建筑施工中是经常遇到的问题，要掌握好混凝土施工缝的设置原则和留置方法，并在施工中保证施工缝的质量。 【关键词】混凝土施工… 实验四 窗函数法设计FIR数字滤波器 一、实验目的 1、掌握窗函数法设计FIR数字滤波器的原理及具体方法。 2、掌握频率取样法设计FIR数字滤波器的原理和基本方法。 3、学习利用窗函数法和频率取样法设计低通、带通、高通、带阻数字滤波器。 二、实验环境 计算机、MATLAB软件 三、实验基础理论 窗函数设计FIR滤波器 1.基本原理 窗函数设计法的基本思想为，首先选择一个适当的理想的滤波器Hd(ej?)，然后 用窗函数截取它的单位脉冲响应hd(n)，得到线性相位和因果的FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器，使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。 2.设计步骤 （1）给定理想滤波器的频率响应Hd(ej?)，在通带上具有单位增益和线性相位， 在阻带上具有零响应。一个带宽为?c(?c??)的低通滤波器由下式给定： 其中?为采样延迟，其作用是为了得到一个因果系统。 （2）确定这个滤波器的单位脉冲响应 Hd(ej?)?e?ja?,|?|??c,Hd(ej?)?0,?c?|?|?? hd(n)? sin(?c(n?a)) ?(n?a) 为了得到一个h(n)长度为N的因果的线性相位FIR滤波器，我们令 a? N?1 2 （3）用窗函数截取hd(n)得到所设计FIR数字滤波器：h(n)?hd(n)RN(n) 3.窗函数的选择 常用的窗函数有矩形（Rectangular）窗，汉宁（Hanning）窗，海明（Hamming）窗、布莱克曼（Blackman）窗、凯瑟（Kaiser）窗等 表4-1 MATLAB中产生窗函数的命令 表4-2 常用窗函数的特性 ?I0? ?(n)?? I0[?] 其中I0?x?是修正的零阶贝塞尔函数，参数?控制最小阻带衰减，这种窗函数对于相同的N可以提供不同的过渡带宽。由于贝塞尔函数比较复杂，这种窗函数的设计方程很难推导，然而幸运的是，有一些经验设计方程可以直接使用。 已知给定的指标?p,?st,Rp和As ，滤波器长度N和凯瑟窗参数?可以按如下凯瑟窗方程给出 过渡带带宽：????st??p N? As?7.95 ?1 2.285?? ??0.1102(As?8.7),As?50???0.4 ??0.5842(As?21)?0.07886(As?21),21?As?50 频率取样设计FIR滤波器 1.基本原理 频率取样法从频域出发，把理想的滤波器Hd(ej?)等间隔采样得到Hd(k),将Hd(k)作为实际设计滤波器的H(k)： H(k)?Hd(k)?H(ej?)| ?? 2? kN k?0,1,?,N?1 j? 得到H(k)以后可以由H(k)来确定唯一确定滤波器的单位脉冲响应h(n),H(e)可以由H(k)求得： h(n)?IDFT[H(k)]H(e)??H(k)?(?? j? k?0N?1 2? k)N 其中?(x)为内插函数： sin(N?/2)?j?N2?1 ?（?)??e Nsin(?/2) 有H(k)求得的频率响应H(ej?)将逼近Hd(ej?)。 如果我们设计的是线性相位FIR滤波器，则H(k)的幅度和相位满足线性相位滤波器的约束条件。 我们将H(k)表示为如下形式 H(k)=H(k)ej?(k)?Hr(k)ej?(k) 当h(n)为实数，则 H(k)?H*(N?k) 由此得到 Hr(k)?Hr(N?k) 即Hr(k)以k?N/2为中心偶对称。在利用线性相位条件可知，对于1型和2型线性相位滤波器： ?N?12?k ?()?2N? ?(k)?? ?(N?1)2?(N?k)?N?2 ?N?1? k?0,?,???2? N?1??k???1,?N?1??2? 对于3型和4型线性相位滤波器 ?