数字信号处理DSP复习第7章.ppt

表7.2.2 6种窗函数的基本参数 MATLAB信号处理工具箱提供了14种窗函数的产生函数，下面列出上述5种窗函数的产生函数及其调用格式： wn=boxcar(N) 　 　％列向量wn中返回长度为N的矩形窗函数w(n) wn=bartlett(N) 　％列向量wn中返回长度为N的三角窗函数w(n) wn=hanning(N) 　％列向量wn中返回长度为N的汉宁窗函数w(n) wn=hamming(N) 　％列向量wn中返回长度为N的哈明窗函数w(n) wn=blackman(N) 　％列向量wn中返回长度为N的布莱克曼窗函数w(n) 7.2.3 用窗函数法设计FIR滤波器的步骤 　　用窗函数法设计FIR滤波器的步骤如下： 　　(1) 根据对过渡带及阻带衰减的指标要求，选择窗函数的类型，并估计窗口长度N。得到窗函数表达式w(n)。 原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。然后根据过渡带宽度估计窗口长度N。待求滤波器的过渡带宽度Bt近似等于窗函数主瓣宽度，且近似与窗口长度N成反比，N≈A/Bt，A取决于窗口类型，例如，矩形窗的A=4π，哈明窗的A=8π等，参数A的近似和精确取值参考表7.2.2。 　　(2) 构造希望逼近的频率响应函数Hd(ejω)，即   对所谓的“标准窗函数法”，就是选择Hd(ejω)为线性相位理想滤波器（理想低通、理想高通、理想带通、理想带阻）。以低通滤波器为例，Hdg(ω)应满足：  　　　　　　　　　　　　　　　　 第7章 有限脉冲响应数字滤波器的设计 第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 (第五章已讲） 7.2 利用窗函数法设计FIR滤波器 7.2 利用窗函数法设计FIR滤波器 7.2.1 窗函数法设计原理 　　设希望逼近的滤波器频率响应函数为Hd(ejω)，其单位脉冲响应是hd(n)。 如果能够由已知的Hd(ejω)求出hd(n)，经过Z变换可得到滤波器的系统函数。但通常以理想滤波器作为Hd(ejω)，其幅度特性逐段恒定，在边界频率处有不连续点，因而hd(n)是无限时宽的，且是非因果序列。 例如，线性相位理想低通滤波器的频率响应函数Hd(ejω)为  　　　　　　　　　　　　　　　　　　　　　  其单位脉冲响应hd(n)为  　　　　　　　　　　　　　　　　　　　　　 线性相位理想高通滤波器的频率响应函数Hd(ejω)为  　　　　　　　　　　　　　　　　　　　　　 其单位脉冲响应hd(n)为  　　　　　　　　　　　　　　　　　　　　　 为了构造一个长度为N的第一类线性相位FIR滤波器，只有将hd(n)截取一段，并保证截取的一段关于 n=(N-1)/2偶对称。设截取的一段用h(n)表示，即  式中, RN(n)是一个矩形序列，长度为N，波形如图7.2.1（b）所示。 由该图可知，当?=（N－1）/2时，截取的一段h(n) 关于n=（N－1）/2偶对称，保证所设计的滤波器具有线性相位。 　　 图7.2.1 窗函数设计法的时域波形（矩形窗，N=30） 这样用一个有限长的序列h(n)去代替hd(n)，肯定会引起误差，表现在频域就是通常所说的吉布斯（Gibbs）效应。该效应引起过渡带加宽以及通带和阻带内的波动，尤其使阻带的衰减小，从而满足不了技术上的要求，如图7.2.2所示。 吉布斯效应是由于将hd(n)直接截断引起的，因此，也称为截断效应。 实际设计的滤波器的单位脉冲响应为h(n)，长度为N，其系统函数为H(z)， 图7.2.2 吉普斯效应 　　RN(n)（矩形序列）就是起对无限长序列的截断作用，可以形象地把RN(n)看做一个窗口，h(n)则是从窗口看到的一段hd(n)序列，所以称h(n)=hd(n)RN(n)为用矩形窗对hd(n)进行加窗处理。 下面分析用矩形窗截断的影响和改进的措施。为了叙述方便，用w(n)表示窗函数，用下标表示窗函数类型，矩形窗记为wR(n)。用N表示窗函数长度。 　　根据傅里叶变换的时域卷积定理，得到（7.2.3）式的傅里叶变换： 式中，Hd(ejω)和WR(ejω)分别是hd(n)和RN(n)的傅里叶变换，即 　　　　　　　　　　　　  WRg(ω)称为矩形窗的幅度函数，如图所示，将图中［－2π/N, 2π/N］区间上的一段波形称为WRg(ω)的主瓣，其余较小的波动称为旁瓣。 将Hd(ej