第六FIR数字滤波器的设计剖析.ppt

第六章 FIR数字滤波器的设计 6-1 FIR数字滤波器的性质 6-2 FIR滤波器的窗函数设计 6-3 FIR滤波器频率采样法设计 6-4 FIR滤波器的等波纹优化设计 2.线性相位FIR滤波器的幅频特性 3.线性相位FIR滤波器的零点特性 4. IIR滤波器与FIR滤波器的比较 6.2 FIR滤波器的窗函数设计 6.2.1 窗函数设计的基本方法 2. 吉布斯（Gibbs）效应 6.2.2 常用的窗函数1. 三角形窗(Bartlett Window) 2. 汉宁(Hanning)窗，又称升余弦窗 3. 汉明(Hamming)窗，又称改进 的升余弦窗 4. 布莱克曼(Blankman)窗，又称二阶升余弦窗 5. 凯泽(Kaiser)窗 6.2.3 几种常用的理想滤波器 1. 理想高通滤波器 2. 理想带通滤波器 3. 理想带阻滤波器 4．理想线性相位线性差分滤波器 5．理想线性相位希尔伯特（Hilbert）  变换器 6.2.4 窗函数法小结与例子 例6-2-1 设计一线性相位FIR数字低通滤波器,截止频率 ，过渡带宽度 ，阻带衰减 dB。 例6-2-2 用凯塞窗函数设计一线性相位FIR数字高通滤波器，截止频率 ，阻带边频 外阻带衰减 不小于60 dB。 例6-2-3 分别用矩形窗和布拉克曼窗设计一个线性相位的希尔伯特变换器，取N=29。 6.3 FIR滤波器频率采样法设计 2、设计过程和公式 综上，可得到设计公式如下： 3、滤波器的频率响应 6.4 FIR数字滤波器的等波纹优化设计 下面介绍切比雪夫等波纹最佳逼近方法的基本思想和方法。 2、利用切比雪夫等波纹最佳逼近法设计  FIR滤波器的原理 3、 求解方法及Remez算法 Remez算法的计算步骤 解：1）确定节数N 因 60dB，可求得 再求得节数N 因对于高通滤波器，N必须为奇数，故取N=27。 2）计算凯塞窗函数 ，将所求 和N值代入下式即可 3）确定延时值 4）计算理想高通的单位脉冲响应 5）求高通数字滤波器的单位脉冲响应： 其幅频响应和相频响应见下图 解：1） 计算对称中心： 2） 按式（6-2-9）计算理想希尔伯特变换器的单位脉冲响应 3）分别取矩形窗和按式（6-2-14a）计算关于 4）计算实际希尔伯特变换器的单位脉冲响应 对称的按布拉克曼窗； 图6-2-10 用矩形窗设计的希尔伯特变换器（N=29） 图6-2-11 用布拉克曼窗设计的希尔伯特变换器（N=29） 从以上的分析和例子可见窗口法设计的最大优点是十分简单实用，其频域特性容易满足要求。 缺点是在大多数情况，设计所得的频响的边界频率往往不能严格控制，同时等波纹最佳设计相比，由窗函数法设计的滤波器长度往往较大。 1、设计思想 按频域采样定理FIR数字滤波器的传输函数H(z)和单位脉冲响应h(n)可由它的N个频域采样值H(k)唯一确定。 （6-3-1） （6-3-2） 假定要设计的FIR数字滤波器的频率响应为 ， 从频域从发,对理想频响在 间进行N点的等间隔 （6-3-3） 采样得到 令 （6-3-4） 按频域采样定理由Hd(k)求H(z)或h(n)，其流图如下： 为使设计的滤波器为线性相位滤波器，滤波器的频率响应 及其频率采样值H(k)必须满足一定的条件. 一般而言，给定的理想滤波器的幅度响应只给出了 范围频率采样的幅度值，因而上式中的k只能取0到N/2 这(N/2+1)点，另外(N/2-1)点(从N/2+1到N-1)则需根据线性相位的奇偶对称条件求得。 下面以偶对称条件h(n)=h(N-n-1)为例说明设计公式。 设h(n)的频响为 （6-3-5） 为幅度响应 （6-3-6） 对 在 等间隔N点采样得H(k) ,k=0,1,…,N-1 对线性相位滤波器，其相位响应为 令 则 按线性相位的偶对称条件: 若N为奇数, 对π为偶对称: 若N为偶数, 对π为奇对称: 因而对频率采样的幅度值 有: N为奇数 N为偶数 按（6-3-6）式对相位进行频率采样得 ，k=0,1,…,(N-1)/2 而 对N为奇数,(N-1)为偶数,