信号和系统实验八-FIR数字滤波器的设计.docxVIP

PAGE PAGE 4 信号与系统实验报告 课 程 名 称： 信号与系统实验 实验项目名称： FIR数字滤波器的设计 专 业 班 级： 姓 名： 学 号： 完 成 时 间： 年 月 日 实验目的 进一步理解FIR滤波器的线性相位特性，熟悉4种线性相位FIR滤波器的幅频特性、相频特性和零极点分布情况； 掌握用窗函数法设计FIR数字滤波器的原理及方法，了解各种窗函数对滤波器性能的影响； 掌握用频率采样法设计FIR数字滤波器的原理及方法，了解不同的采样点数和不同的过渡带点数对滤波器性能的影响。 实验原理 FIR滤波器的设计问题在于寻求一系统函数，使其频率响应逼近滤波器要求的理想频率响应，其对应的单位脉冲响应为。 （1）用窗函数设计FIR滤波器的基本原理 设计思想：从时域从发，设计逼近理想。设理想滤波器的单位脉冲响应为。以低通线性相位FIR数字滤波器为例。 (6-1) 一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器h(n)，最直接的方法是截断，即截取为有限长因果序列，并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即 (6-2) 用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，这个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。 （2） 典型的窗函数 ① 矩形窗(Rectangle Window) （6-3） ② 三角形窗(Bartlett Window) （6-4） ③ 汉宁(Hanning)窗，又称升余弦窗 （6-5） ④ 汉明(Hamming)窗，又称改进的升余弦窗 （6-6） ⑤ 布莱克曼(Blankman)窗，又称二阶升余弦窗 （6-7） ⑥ 凯泽(Kaiser)窗 （6-8） 其中：β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，β越大，过渡带越宽，阻带越小衰减也越大。I0(·)是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为，β的确定可采用下述经验公式： （6-9） 若滤波器通带和阻带波纹相等即δp=δs时，滤波器节数可通过下式确定： （6-10） 式中： （3）利用窗函数设计FIR滤波器的具体步骤如下： 1、确定数字滤波器的性能要求，临界频率，滤波器单位脉冲响应长度N。 2、根据性能要求，合理选择单位脉冲响应h(n)的奇偶对称性，从而确定理想频率响应的幅频特性和相频特性。 3、求理想单位脉冲响应，在实际计算中，可对采样，并对其求IDFT的，用代替。 4、选择适当的窗函数w（n），根据求所需设计的FIR滤波器单位脉冲响应。 5、求，分析其幅频特性，若不满足要求，可适当改变窗函数形式或长度N，重复上述设计过程，以得到满意的结果。 实验内容 选择合适的窗函数设计一个线性相位FIR低通滤波器，通带截止频率w,?=0.3mad.阻带截止频率w,?=0.?45mrad,阻带最小衰减a,?=?50dB。要求分别绘制理想冲激响应曲线、窗函数曲线、实际冲激响应曲线、FIR低通滤波器的幅频响应曲线。 2.选择合适的窗函数设计一个线性相位FIR带阻滤波器，下通带的截止频率f2?=1kHz,上通带的截止频率f2?=4.5kHz,通带内的最大衰减a,?=0.5dB,阻带的下限截止频率f?=2kHz,阻带的上限截止频率fe?=3.?5kHz,阻带最小衰减a。=65dB,采样频率F.?=?10kHz。要求分别绘制窗函数曲线、实际滤波器的冲激响应曲线、滤波器的幅频响