张军《数字信号处理》chap5 FIR—DF的设计方法-end.pptVIP

表5-2 凯塞窗的性能 虽然凯塞窗看上去没有初等函数的解析表达式。但是，在设计凯塞窗时，对零阶变形贝塞尔函数可采用无穷级数来表达 （5-50） 这个无穷级数，可用有限项级数去近似，项数多少由要求的精度来确定。 因而采用计算机是很容易求解的。 表5-3 六种窗函数基本参数的比较 最小阻带衰减只由窗形状决定，不受N的影响；而过渡带的宽度则随窗宽的增加而减小。 下面将窗函数法的设计步骤归纳如下： （1） 给定希望逼近的频率响应函数Hd(ejΩ) 。 （2） 求单位脉冲响应hd(n)。 如果Hd(ejΩ)很复杂或不能直接计算积分，则必须用求和代替积分，以便在计算机上计算，也就是要计算离散傅里叶反变换IDTFT。 （3） 由过渡带宽及阻带最小衰减的要求，可选定窗形状， 并估计窗口长度N。按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的窗函数。 （4） 计算所设计的FIR滤波器的单位脉冲响应。 h(n)=hd(n)w(n) 0≤n≤N-1 （5）由h(n)求FIR滤波器的系统函数H(z) （6）求FIR滤波器的频率响应，检验是否满足设计要求， 如不满足，需进行修正。 通常整个设计过程可利用计算机编程来实现，可多选择几种窗函数来试探，从而设计出性能良好的FIR 滤波器。 例5-1：设计一个线性相位FIR低通滤波器通带截止频率为 ，阻带截止频率为 ， 阻带衰减 。 窗口法设计的主要优点是简单，使用方便。窗口函数大多有封闭的公式可循，性能、参数都已有表格、资料可供参考， 计算程序简便， 所以很实用。缺点是通带和阻带的截止频率不易控制。 N=33; b1=fir1(32,1/4,boxcar(33)); %用矩形窗设计FIR低通数字滤波器 5.3 用窗函数法设计FIR滤波器 一、设计方法 设计FIR数字滤波器最简单的方法是窗函数法。这种方法一般是先给定所要求的理想滤波器的频率响应Hd(ejΩ)，要求设计一个FIR滤波器频率响应 , 去逼近理想的频率响应Hd(ejΩ)。然而，窗函数法设计FIR数字滤波器是在时域进行的，因此，必须首先由理想频率响应Hd(ejΩ) 的傅里叶反变换推导出对应的单位脉冲响应hd(n) （5-24） 由于许多理想化的系统均用分段恒定的或分段函数表示的频率响应来定义，因此这种系统具有非因果的和无限长的脉冲响应，即hd(n)一定是无限长的序列，且是非因果的。而我们要设计的是FIR滤波器，其h(n)必定是有限长的，所以要用有限长的h(n)来逼近无限长的hd(n)，最简单且最有效的方法是截断hd (n) 0≤n≤N-1 其他 通常，我们可以把h(n)表示为所需单位脉冲响应与一个有限长的窗口函数序列w(n)的乘积，即 h(n)=hd(n)w(n) （5-26） （5-25） 式中，如果采用如式（5-25）的简单截取，则窗函数为矩形窗。 0≤n≤N-1 其他 （5-27） 矩形窗的波形如图5-3（b）所示。  例如，要求设计一个FIR低通数字滤波器，假设理想低通滤波器的频率响应为 | W |≤ W c W c|W |≤π （5-28） 相应的单位脉冲响应hd(n)为 （5-29） 这是一个中心点在a的偶对称、无限长、非因果序列，hd(n)的波形如图5-3（a）所示。为了构造一个长度为N的线性相位滤波器，只有将hd(n)截取一段，并保证截取的一段对(N-1)/2 对称， 故中心点a必须取a=(N-1)/2。设截取的一段用h(n)表示，如式（5-26）所示，h(n)的波形如图5-3（c）所示。 图 5-3 理想低通的单位脉冲响应及矩形窗 五种常用的窗函数 由复卷积定理可知，时域相乘，频域是周期卷积，故h(n)的频率特性为 （5-30） H(ejΩ)能否逼近Hd(ejΩ) 取决于窗函数的频谱特性W(ejΩ) 这里选用矩形窗RN(n)，其频谱特性为 （5-31） 幅频特性和相频特性为 （5-32） 式中： 其中，WR(W)是周期函数，如图5-4（b）所示。主瓣宽度为4π/N，两侧有许多衰减振荡的旁瓣。通常主瓣定义为原点两边第一个过零点之间的区域。 （5-33） 图5-4 矩