第7章有限脉冲响应数字滤波器的设计C概要.pptVIP

第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用等波纹最佳逼近法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 2）　FIR滤波器具有线性相位的条件是h(n)为实序列，且满足h(n)=-h(N－n－1)，其频响函数应满足的条件是： 　频域采样定理表明，频率域等间隔采样H(k)，经过IDFT得到h(n)，得到H(ejω)=FT［h(n)］的内插表示形式： 在采样点上实际频率响应严格等于理想频率响应 采样点之间实际频响由各采样点加权内插函数延伸而成，因而有误差，误差大小和理想频率响应的曲线形状有关 改进：在理想频率响应的不连续的边缘点加上一些过渡的采样点 第7章 有限脉冲响应数字滤波器的设计 图7.3.1 频域幅度采样序列Hg(k)及其内插波形Hg(ω) Hdg(ω)特性愈平滑的区域，误差愈小；特性曲线间断点处，误差最大。 第7章 有限脉冲响应数字滤波器的设计 　　表现形式为间断点变成倾斜下降的过渡带曲线，过渡带宽度近似为2π/N。通带和阻带内产生震荡波纹，且间断点附近振荡幅度最大，使阻带衰减减小，往往不能满足技术要求。 　　当然，增加N可以使过渡带变窄，但是通带最大衰减和阻带最小衰减随N的增大并无明显改善。且N太大，会增加滤波器的阶数，即增加了运算量和成本。 N=15和N=75两种情况下的幅度内插波形Hg(ω)如图7.3.2所示，图中的空心圆和实心圆点分别表示N=15和N=75时的频域幅度采样。 第7章 有限脉冲响应数字滤波器的设计 图7.3.2 N=15和N=75的幅度内插波形Hg(ω) 第7章 有限脉冲响应数字滤波器的设计 　　在窗函数设计法中，通过加大过渡带宽度换取阻带衰减的增加。频率采样法同样满足这一规律。 提高阻带衰减的具体方法：在频响间断点附近区间内插一个或几个过渡采样点，使不连续点变成缓慢过渡带，这样，虽然加大了过渡带，但阻带中相邻内插函数的旁瓣正负对消，明显增大了阻带衰减。 　　过渡带采样点的个数与阻带最小衰减? s的关系以及使阻带最小衰减? s最大化的每个过渡带采样值求解都要用优化算法解决。其基本思想是将过渡带采样值设为自由量，用一种优化算法（如线性规划算法）改变它们，最终使阻带最小衰减? s最大。 第7章 有限脉冲响应数字滤波器的设计 　将过渡带采样点的个数m与滤波器阻带最小衰减? s的经验数据列于表7.3.1中，我们可以根据给定的阻带最小衰减? s选择过渡带采样点的个数m。 表7.3.1 过渡带采样点的个数m与滤波器阻带最小衰减? s的经验数据 第7章 有限脉冲响应数字滤波器的设计 　　四． 频率采样法设计步骤 　　 　 （1） 根据阻带最小衰减? s选择过渡带采样点的个数m。 　 （2） 确定过渡带宽度Bt，估算频域采样点数（即滤波器长度）N。 如果增加m个过渡带采样点，则过渡带宽度近似变成(m+1)2π/N。当N确定时，m越大，过渡带越宽。如果给定过渡带宽度Bt，则要求(m+1)2π/N≤Bt ，滤波器长度N必须满足如下估算公式: （7.3.15） 第7章 有限脉冲响应数字滤波器的设计 　　（3） 构造一个希望逼近的频率响应函数: 设计标准型片断常数特性的FIR数字滤波器时，一般构造幅度特性函数Hdg(ω)为相应的理想频响特性，且满足表7.1.1要求的对称性。 　 （4） 按照（7.3.1）式进行频域采样: （7.3.16） （7.3.17） 第7章 有限脉冲响应数字滤波器的设计 并加入过渡带采样。过渡带采样值可以设置为经验值，或用累试法确定，也可以采用优化算法估算。 　 （5） 对H(k)进行N点IDFT，得到第一类线性相位FIR数字滤波器的单位脉冲响应： 　　　　　　　　　　　　　　　　　 （7.3.18） 　　（6） 检验设计结果。如果阻带最小衰减未达到指标要求，则要改变过渡带采样值，直到满足指标要求为止。如果滤波器边界频率未达到指标要求，则要微调Hdg(ω)的边界频率。 上述设计过程中的计算相当繁琐，所以通常借助计算机设计。 第7章 有限脉冲响应数字滤波器的设计 　　【例7.3.1】 用频率采样法设计第一类线性相位低通FIR数字滤波器，要求通带截止频率ωp=π/3，阻带最小衰减大于40 dB，过渡带宽度Bt≤π/16。 　　解 查表7.3.1， ? s=40 dB时，过渡带采样点数m=1。 根据m=1和Bt≤π/16估算滤波器长度: N≥(m+1)2π/Bt=64 留一点富余量，取N=65。 构造Hd(ejω)=Hdg(ω)e－jω(N－