采用频率采样法的FIR滤波器.doc

吹管乐滤波去噪 ——基于频率采样法的FIR滤波器 学生姓名：焦阳 指导老师：胡双红 摘 要 本课程设计主要是设计并。本课程设计仿真平台为MATLAB。首先。可知。Finite Implse Response）H（k）H（k）H（k）H（z） （2-1） 这个公式就为FIR滤波器提供了另外一种结构，这种结构由两部分级联组成。 (2-2) 其中级联的第一部分为梳状滤波器，其结构如图所示： （2-3） X(n) 1 r(n) -Z-N 图2-1 梳状滤波器结构图 第二部分由N各谐振器组成的谐振柜。 它是由N个一阶网络并联而成，而这每一个一阶网络都是一个谐振器 （2-4） 其结构如下图所示： H（k） Hk(z) W-A 图2-2 一阶谐振器 频率抽样型结构特点： （1）它的系数H（k）直接就是在滤波器在处得频率响应。因此，控制 得频率响应是很直接得。 （2）结构有两个主要缺点： a．所有的相乘系数及H（k）都是复数，应将它们先化成二阶实数，这样乘起来较复杂，增加乘法次数，存储量。 b．所有谐振器的极点都是在单位圆上，由决定考虑到系数量化的影响，当系数量化时，极点会移动，有些极点就不能被梳状滤波器的零点所抵消。（零点由延时单元决定，不受量化的影响）系统就不稳定了。 （3）将一阶网络合并为二阶网络： a．第k和第N-k个谐振器合并为一个实系数的二阶网络，因为h(n)是实数，他的DFT也是圆周共轭对称的。 （2-5） 因此，可以将第k和第N-k个谐振器合并为一个二阶网络。 （2-6） 其中： b．第k和第N-k个谐振器合并为一个二阶网络的极点在单位圆内，而不是在单位圆上，因而从频率响应的几何解释可知，它相当于一个有限Q的谐振器。其谐振频率为 。 图2-3 二阶网络结构图 除了共轭复根外，还有实根。 当N=偶数时，有一对实根，它们分别是k=0，k=N/2两个点。 和 （2-7） 图2-4 修正频率抽样结构流图（N=偶数） （2-8） 修正频率抽样结构流图（N=奇数） 图2-5 修正频率抽样结构流图（N=奇数） （2-9） 2.3 频率采样法 2.3.1 设计思路： 这种设计方法是从频域进行设计的一种方法，首先给定一个希望逼近的频率响应。 （2-10） 知道H（k）IDFT定义，可以用这N个采样值H（k） （2-11） （2-12） （2-13） 内插公式： （2-15） 四种线性相位的FIR滤波器如下表2-1所示： 表2-1 四种线性相位的FIR滤波器 2.3.2 逼近误差及其改进措施： （1）采样点上滤波器的实际响应是严格地和理想频率响应数值相等的； （2）在采样点之间的频率响应则是由各采样点的加权内插函数的延伸叠加而成的，因而有一定的逼近误差，误差大小取决于理想频率响应曲线形状。 （3）理想频率响应特性变化越平缓，则内插值越接近理想值，逼近误差越小； （4）如果采样点之间的理想频率特性变化越陡，则内插值与理想值的误差就越大，因而在理想频率特性的不连续点附近，