有限冲激响应滤波器及设计FIR.ppt

* * 1.加窗对滤波器频响产生的影响 （1）出现过渡带，宽度等于WR(ω)的主瓣宽度（对于 矩形窗 ）； （2）通带和阻带内产生波动，其振荡幅度取决于旁瓣 的相对幅度。 2.对窗函数的要求 主瓣宽度小，以获得较陡的过渡带； 与主瓣的幅度相比，旁瓣应尽可能小，把能量尽量集中在主瓣中，以减小通带和阻带中的波纹幅度。 吉布斯（Gibbs）效应 由于对hd(n) 截短，导致了所设计滤波器幅频特性的起 伏波动，波动的幅度强弱完全取决于窗函数的类型，而与窗 的宽度N无关。这种现象称为吉布斯（Gibbs）效应。 N=7 N=21 N=51 N=101 * 三、常用窗函数 （1）矩形窗 主瓣宽度为 * （2）升余弦窗（汉宁窗Hanning Window） 主瓣宽度为 * （3）改进的升余弦窗（海明Hamming窗） 主瓣宽度为 * （4）二阶升余弦窗（布莱克曼Blackman窗） 主瓣宽度为 * Hamming Blackman 图6-10 常用窗函数的时域波形 图6-11 常用窗函数的频谱 Hamming Blackman Rectangle Hanning 图6-12 理想低通加窗后的幅度响应（N=51） Hamming Blackman Rectangle Hanning * （5）凯塞窗（Kaiser Window） 这是一种适应性较强的窗，其窗函数的表示式为 式中，I0(x)是第一类变形零阶贝塞尔函数， 是一个可 自由选择的参数。 凯塞窗可以在主瓣宽度和旁瓣衰减之间自由选择。 * 表6-2 凯塞窗的性能 过渡带 通带波纹/dB 阻带最小衰减/dB 2.120 3.00π/N ±0.27 -30 3.384 4.46π/N ±0.0864 -40 4.538 5.86π/N ±0.0274 -50 5.658 7.24π/N ±0.00868 -60 6.764 8.64π/N ±0.00275 -70 7.865 10.0π/N ±0.000868 -80 8.960 11.4π/N ±0.000275 -90 10.056 12.8π/N ±0.000087 -100 * 表6-3 几种窗函数基本参数的比较 窗函数 旁瓣峰值幅度/dB 过渡带宽 阻带最小衰减/dB 矩形窗 -13 4π/N -21 汉宁窗 -31 8π/N -44 海明窗 -41 8π/N -53 布莱克曼窗 -57 12π/N -74 凯塞窗 -57 10π/N -80 要求：熟悉各种窗函数的技术指标和加窗后对滤波特性 的影响，能根据设计指标正确选择窗函数类型及其长度N。 四、窗函数法设计线性相位FIR滤波器的一般步骤 为了准确控制滤波器通带边缘，常需进行多次设计。 若Hd(ejw)不能用简单函数表示，则可用求和运算代替积分运算。 * 例 题 例6-4 用矩形窗设计一个线性相位带通滤波器 （1）设计N为奇数时的h(n)。  （2）设计N为偶数时的h(n)。  （3）若改用海明窗设计，求以上两种形式的h(n)表达式。 * 解： h(n)=hd(n)RＮ(n) 注意：N取奇、偶，虽然两个表达式形式完全一样，但 在 0≤n≤N-1 上的取值上完全不同。 * 例 6-5 根据下列技术指标，设计一个FIR低通滤波器。 通带截止频率wp=0.2p，通带允许波动Ap=0.25dB；  阻带截止频率ws=0.3p，阻带衰减As=50dB。 解：查表6-3可知，海明窗和布拉克曼窗均可提供大于 50dB的衰减。但海明窗具有较小的过渡带从而具有较小的长 度N。 根据题意，所要设计的滤波器的过渡带为 由表6-3可知，利用海明窗设计的滤波器的过渡带宽 Δw=8p/N，所以低通滤波器单位脉冲响应的长度为 * 3 dB通带截止频率为 由式（6-29）可知，理想低通滤波器的单位脉冲响应为 海明窗为 则所设计的滤波器的单位脉冲响应为 * 例6-6 用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器 设N=11，wc=0.2prad。 解：Matlab程序如下 %汉宁窗 N=11; n=[0:1:N-1];Wc=0.2*pi; hd=ideal_lp(Wc,N);w_han=(hanning(N)); h=hd.*w_han; [db,mag,pha,grd,w]=freqz_m(h,1);plot(w/pi,db); %矩形窗 w_han=(bo