信号分析与处理 教学课件 作者 杨西侠 柯晶 5 2FIR.pptVIP

5.3 FIR数字滤波器设计 5.3.1 线性相位FIR数字滤波器的特点 5.3.2 窗函数法设计FIR数字滤波器 5.3.3 频率抽样法设计FIR数字滤波器1. 设计原理 窗函数 主瓣宽度 (2?/N) 旁瓣电平 （dB） 矩形窗 2 ?13 0.9 ?21 三角窗 4 ?25 2.1 ?25 汉宁窗 4 ?31 3.1 ?44 汉明窗 4 ?41 3.3 ?53 布莱克曼窗 6 ?57 5.5 ?74 凯泽窗 ?57 5 ?80 ( ?=7.865) 加窗后滤波器性能指标 过渡带宽度 阻带最小衰减 ? ? (2?/N) （dB） 3. 窗函数法的设计步骤 （1）给定要求的Hd(e jω)，计算FIR digital filter的单位冲激响应hd(n)。 若Hd(e jω)复杂，不便于直接求阶积分，可以对Hd(ejω)进行M点频域抽样，求出上述积分的数值近似值hdp(n)，正好是M点的IDFT，即： 根据频域抽样理论，由频率抽样所求得的hdp(n)是原序列h(n)的周期延拓，则 = IDFT[ Hd (k) ] 若M N，在窗口内，就有hdp(n)有效的逼近hd(n)。 （2）选窗口函数及窗口长度 根据允许的过渡带宽度及阻带衰减要求选择窗函数形状及filter的长度N。（N可通过尝试） （3）按所得窗函数求得 h(n)= hd(n)·w(n) （4）计算 检验各项指标。 窗函数法设计简单实用，但缺点是过渡带及边界频率不易控制，通常需要反复计算。 【例5-7】用窗函数法设计一线性相位数字低通滤波器。给定技术指标： 通带允许起伏 ?1dB 0 ? ? ? 0.3? 阻带衰减 ? 50dB 0.5? ? ? ? ? 解：用窗函数法设计，边界频率不易准确控制，近似取理想低通滤波器的截止频率为 ① 计算hd(n) 以?为中心的偶对称无限长序列，为得到要求的线性相位，取 ? = (N?1)/2 ② 根据阻带衰减要求选窗函数及窗口长度 查表，选hamming窗 ?ω = 3.3 · 2π/N N = 33 ? = 16 ③ 截断求h(n) n= 0, 1, 2, …32 ④ 求 H(e jω) boxcar triang hanning blackman N=22 N=55 频率抽样法是以频率抽样理论为基础的。 据序列傅里叶变换与DFT的关系，一个任意长的序列，对它的频率特性进行N等分间隔抽样，利用IDFT，可以得到一个N点的有限长序列。这个有限长序列是原序列以N为周期的周期序列的主值序列，是原序列的近似，因而它的频率特性也将逼近原序列所对应的频率特性。 一个理想频响Hd(e jω) 单位抽样响应hd(n) N分度间隔抽样 N个H(k) IDFT N点 h(n) H(e jω) ≈ 对Hd(e jω)抽样所得H(k) k = 0, 1, 2, …N-1 内插公式 Hd(e jω) 2．线性相位的约束