FIR数字滤波器的（三角）窗函数法设计.docVIP

课程设计说明书 NO. PAGE 21 沈 阳 大 学 语音信号的数字滤波 ——FIR数字滤波器的（三角）窗函数法设计 课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；掌握利用数字滤波器对信号进行滤波的方法。并能够对设计结果加以分析。 2.设计步骤及要求 2.1课程设计的要求 巩固和加深对理论课中知识的理解，提高学生对所学理论知识的综合运用能力。通过课程设计，使学生具有下列基本技能：培养学生查阅参考资料、手册的自学能力，学会自己分析、解决问题的方法；要求学生能够熟练地用Matlab语言编写数字信号处理的应用程序；了解各种窗函数对滤波器的影响。 2.2 FIR滤波器的原理和分类 2.2.1 FIR滤波器设计步骤流程图 图1所示 图1 设计步骤流程 2.2.2 FIR 滤波器的基本结构 FIR滤波器的系统函数为： 。FIR滤波器有以下特点： (1)系统的单位冲激响应h(n)在有限个n值处不为零； (2)系统函数H(z)在|z|0处收敛，极点全部在z=0处(稳定系统)； 结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中含有反馈的递归。 2.2.3 FIR滤波器实现的基本结构有： (1)横截型（卷积型、直接型） a.一般FIR滤波器的横截型结构： 给定差分方程为：? 。 b.线性相位FIR滤波器的横截型结构 ① N为奇数时线性相位FIR滤波器实现结构如2图所示： ??② N为偶数时线性相位FIR滤波器实现结构如3图所示： 图2 N为基数时 图3 N为偶数时 （2）级联型 将H(z)分解为若干个实系数一阶或二阶因子相乘： 级联结构如下图4所示： 图4 级联结构 （3）频率取样型 ? 若FIR滤波器的冲激响应为有限长(N点)序列h(n)，则有如5图所示的关系： 图5 关系图 ??? 因此，对h(n)可以利用DFT得到H(k)，然后利用内插公式： （1） 来表示系统函数，这就为FIR滤波器提供了另外一种结构：频率抽样结构， ? 这种结构由两部分级联而成：分析系统函数 （4）快速卷积结构 若FIR滤波器的单位冲激响应h(n)是一个N1点有限长序列，输入x(n)是一个N2点有限 长序列，那么输出y(n)是x(n)与h(n)的线性卷积，它是一个L＝N1+N2-1点的有限长序 列。而圆周卷积可以用DFT和IDFT的方法来计算，得到FIR滤波器的快速卷积结构如图 6所示： 图6 快速卷积结构 2.2.4 线性相位FIR数字滤波器的条件和特点 （1）线性相位条件 对于长度为N的h(n)，传输函数为 （2） （3） 式中，称为幅度特性， 称为相位特性。其中，这里不同于， 为的实函数，可能取负值，而总是正值。线性相位是指是 的线性函数，即 =，为常数 （4） 如果满足下式 =是起始相位 （5） 此时不具有线性相位，但以上两种情况都满足群时延是一个常数，即 （6） 满足第一类线性相位的条件是： （7）