基于faool的窗函数法fir滤波器设计教学.docxVIP

基于faool的窗函数法fir滤波器设计教学 1. 基于fpga设计的数字滤波器设计 “数字预处理”是电子信息安全专业最重要的专业课程。许多大学教授都有这样的课程。其主要内容是阐述数字信号处理基本理论及其两个主要应用, 即:谱分析和数字滤波器设计。该课程的特点是:概念抽象, 数学计算量大, 并且涉及大量的公式推导和矩阵计算, 计算复杂。学生在学习这门课程时普遍感到很吃力, 对其中的分析方法与基本理论不能很好地理解与掌握。如何帮助学生理解与掌握课程中的基本概念、基本原理和基本分析方法, 提高综合运用所学知识解决实际问题的能力, 是本课程教学中所要解决的主要问题。 数字滤波器设计是数字信号处理课程核心内容之一。对于滤波器的设计, 比如FIR滤波器, 很多教材介绍窗函数设计方法, 在理想滤波器单位脉冲响应基础上, 通过加窗来得到所要设计滤波器的系数。但从理论推导到滤波器系数计算, 学生所见到的只是一堆公式和数字, 很难直观地理解窗长、窗型等不同因素对FIR滤波器设计的影响度, 不利于学生有效掌握其设计思想和应用。 本文介绍了一种采用FDATool工具辅助窗函数法FIR滤波器设计教学的方法。利用FDATool强大的图形功能, 可以使学生能够更直观的理解窗函数法FIR滤波器设计重点内容, 进而提高学生的学习兴趣, 加强理论和实际结合的能力, 培养学生的综合设计能力与实际工作能力。 2. 基于窗函数法设计菲尔滤波 2.1 数字滤波器的设计原则 窗函数设计FIR数字滤波器的基本思想是根据给定的滤波器技术指标选择窗型和确定滤波器阶数N。其核心是由给定的频率特性通过对理想滤波器的单位脉冲响应加窗来确定有限长单位脉冲响应h (n) 。下面介绍窗函数设计FIR数字低通滤波器的过程: 设:理想低通的传输函数为: Ηd(ejω)={1,|ω|≤ωc0,ωc＜ω≤πHd(ejω)={1,|ω|≤ωc0,ωc＜ω≤π 则其脉冲响应: hd(n)=ΙDFΤ[Ηd(ejω)]=12πωc∫-ωcejωndω=sin(ωcn)πnhd(n)=IDFT[Hd(ejω)]=12π∫?ωcωcejωndω=sin(ωcn)πn 对理想低通的脉冲响应加窗得所设计滤波器的单位脉冲响应: h(n) =hd(n)w(n) , 其中w(n) 是长度为N的窗函数。 在时域h(n) 以 (N-1) /2为对称, 保证FIR数字滤波器为线性相位;对h(n) 右移 (N-1) /2以保证所设计FIR数字滤波器的因果性。 为了省去每次的移位, 可对理想低通在频域作相移, 即: Ηd(ejω)={e-jωα0≤|ω|≤ωc0ωc≤|ω|≤πα=(Ν-1)/2Hd(ejω)={e?jωα0≤|ω|≤ωc0ωc≤|ω|≤πα=(N?1)/2 如前所述, 可得所设计滤波器的单位脉冲响应: h(n)=hd(n)w(n)=sin(ωc(n-a))π(n-a)w(n)h(n)=hd(n)w(n)=sin(ωc(n?a))π(n?a)w(n) 在频域检查H(ejω) =DFT[h(n) ]是否满足技术滤波器设计的指标要求, 若不满足, 需修改设计, 直至满足技术指标要求为止。 2.2 过渡带的设计 由离散时间付立叶变换的性质可知:时域加窗等于频域卷积。所以, 窗函数法所设计FIR滤波器的幅度特性一定是矩形窗函数的幅度特性与理想低通滤波器的幅度特性的卷积。这对理想频率响应产生了以下两点影响: ①使理想频率特性不连续点ω=ωc处, 形成了一个过渡带。过渡带的宽度等于矩形窗的频率响应RN(ω) 的主瓣宽度△ω= 4π/N; ②在截止频率ωc的两边ω=ωc±2π/N处 (即过渡带的两边) ,H(ejω) 的幅度特性出现最大的肩峰值, 肩峰的两侧形成起伏振荡, 其振荡幅度取决于旁瓣的相对幅度, 而振荡的快慢, 则取决于RN(ω) 波动的快慢。这直接影响滤波器的性能。 由上述可知, 对FIR滤波器设计中的窗函数有如下要求: ①窗谱主瓣尽可能的窄, 以获得较陡的过渡带; ②尽量减小窗谱的最大旁瓣的相对幅度, 也就是能量尽量集中在主瓣, 这样使肩峰和波纹减小, 可以增大阻带的衰减。 窗函数法FIR滤波器设计中常用窗函数的窗谱性能指标性如表1。 3. fir滤波器 FDATool是MATLAB数字信号处理工具箱中一种图形化的滤波器设计与分析工具, 使用该工具可以快速设计各种类型的滤波器, 并计算出滤波器系数, 还可以画出滤波器的幅度/相位相应、群/相位延迟和零极点分布图等。在FDATool设计界面下, 可分别设计FIR和IIR滤波器。 在FDATool图形界面下, 根据滤波器的频带参数, 通带截止频率ωc以及滤波器的阶数等, 即可以设计出符合要求的滤波器。 下面给出截止频率的ωc=0.2π线性相位FIR数字低通滤波器的