数字信号处理第7章有限脉冲响应数字滤波器的设计要点解析.ppt

在ω=0~2π区间上N个等间隔的采样频点为   将ω=ωk代入（7.3.5）~（7.3.8）式中，并写成k的函数： （7.3.12） k = 0，1，2，…，N – 1 （7.3.9） （7.3.10） （7.3.11） N = 奇数 N = 偶数 N等于奇数时Hg(k)关于N/2点偶对称。 N等于偶数时，Hg(k)关于N/2点奇对称，且Hg(N/2)=0。  　　设用理想低通作为希望逼近的滤波器Hd(ejω)，截止频率为ωc，采样点数为N，Hg(k)和θ(k)用下列公式计算： N=奇数时， （7.3.13） N=偶数时，   上面公式中的kc是通带内最后一个采样点的序号，所以kc的值取不大于［ωcN/(2π)］的最大整数。另外，对于高通和带阻滤波器，这里N只能取奇数。 （7.3.14） 　　3． 逼近误差及其改进措施 　　如果待逼近的滤波器为Hd(ejω)，对应的单位脉冲响应为hd(n)，则由频率域采样定理知道，在频域0~2π范围等间隔采样N点，利用IDFT得到的h(n)应是hd(n)以N为周期的周期延拓的主值区序列，即  如果Hd(ejω)有间断点，那么相应的单位脉冲响应hd(n)应是无限长的。这样，由于时域混叠及截断，使h(n)与hd(n)有偏差。所以，频域的采样点数N愈大，时域混叠愈小，设计出的滤波器频响特性愈逼近Hd(ejω)。 　　上面是从时域分析其设计误差的来源，下面从频域分析。频域采样定理表明，频率域等间隔采样H(k)，经过IDFT得到h(n)，由（3.3.7）和（3.3.8）式，得到H(ejω)=FT［h(n)］的内插表示形式：   式中  上式表明，在采样频点，ωk=2πk/N, k=0，1，2，…，N－1, Φ(ω－2πk/N)=1，因此采样点处　　　　与H(k)相等，逼近误差为0。在采样点之间，　　　　由N项H(k)Φ(ω－2πk/N)之和形成。频域幅度采样序列Hg(k)及其内插波形Hg(ω)如图7.3.1所示。图7.3.1的上图 7.3.1(a) 中, 实线表示希望逼近的理想幅度函数Hdg(ω), 黑点表示幅度采样序列Hg(k); 下图7.3.1(b)中, 实线Hg(ω）与虚线Ηdg(ω)的误差与Hdg(ω)特性的平滑程度有关，Hdg(ω)特性愈平滑的区域，误差愈小；特性曲线间断点处，误差最大。 　　表现形式为间断点变成倾斜下降的过渡带曲线，过渡带宽度近似为2π/N。通带和阻带内产生震荡波纹，且间断点附近振荡幅度最大，使阻带衰减减小，往往不能满足技术要求。 　　当然，增加N可以使过渡带变窄，但是通带最大衰减和阻带最小衰减随N的增大并无明显改善。且N太大，会增加滤波器的阶数，即增加了运算量和成本。N=15和N=75两种情况下的幅度内插波形Hg(ω)如图7.3.2所示，图中的空心圆和实心圆点分别表示N=15和N=75时的频域幅度采样。 　　运行本书程序库程序fig731-2.m，即可绘制出图7.3.1和图7.3.2，并分别输出两种采样点数（N=15和N=75）的通带最大衰减? p和阻带最小衰减? s。N=15时，通带最大衰减? p=0.8340 dB，阻带最小衰减? s=－15.0788 dB； N=75时，通带最大衰减 ? p=1.0880 dB，阻带最小衰减? s=－16.5815 dB。 所以，直接对理想滤波器的频率响应采样的“基本频率采样设计法”不能满足一般工程对阻带衰减的要求。 图7.3.1 频域幅度采样序列Hg(k)及其内插波形Hg(ω) 图7.3.2 N=15和N=75的幅度内插波形Hg(ω) 　　在窗函数设计法中，通过加大过渡带宽度换取阻带衰减的增加。频率采样法同样满足这一规律。提高阻带衰减的具体方法是在频响间断点附近区间内插一个或几个过渡采样点，使不连续点变成缓慢过渡带，这样，虽然加大了过渡带，但阻带中相邻内插函数的旁瓣正负对消，明显增大了阻带衰减。 　　过渡带采样点的个数与阻带最小衰减? s的关系以及使阻带最小衰减? s最大化的每个过渡带采样值求解都要用优化算法解决。其基本思想是将过渡带采样值设为自由量，用一种优化算法（如线性规划算法）改变它们，最终使阻带最小衰减? s最大。 表7.3.1 过渡带采样点的个数m与滤波器阻带最小衰减? s 　　　　　的经验数据 　　4． 频率采样法设计步骤 　　综上所述，可归纳出频率采样法的设计步骤: 　 （1） 根据阻带最小衰减? s选择过渡带采样点的个数m。 　 （2） 确定过渡带宽度Bt，估算频域采样点数（即滤波器长度）N。如果增加m个过渡带采样点，则过渡带宽度近似变成(m+1)2π/N。当N确定时，m越大，过渡带越宽。如果给定过渡带宽