实验六FIR滤波器的设计.docxVIP

实验六FIR滤波器的设计 实验六、用窗函数法设计FIR 数字滤波器 一、实验目的： （1）熟悉基本的窗函数，及其特点。 （2）掌握用窗函数法设计FIR 数字滤波器的原理和方法。 （3）熟悉线性相位FIR 数字滤波器特性。 二.实验原理 （一）FIR滤波器的设计 FIR滤波器具有严格的相位特性，这对于语音信号处理和数据传输是很重要的。目前FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不高的时候是比较灵活方便的。它是从时域出发，用一个窗函数截取一个理想的hd(n)得到h(n)，以有限长序列h(n)近似理想的hd(n)：如果从频域出发，用理想的hd(ejw)在单位圆上等角度取样得到H（k），根据h(k)得到H(z)将逼近理想的Hd(z)这就是频率取样法。 （二）窗函数设计法 同其它的的数字滤波器设计的方法一样，用窗函数设计滤波器也是首先要对滤波器提出性能指标。一般是给定一个理想的频率响应hd(ejw)，使所设计的FIR滤波器的频率响应 h(ejw)去逼近所要求的理想的滤波器的响应hd(ejw)窗函数设计的任务在于寻找一个可实 现（有限长单位脉冲响应）的传递函数。 H(e)??h(n)e?jwn jwn?0N?1去逼近hd(ejw)。我们知道，一个理想的频率响应hd(ejw)的傅里叶变换 1hd(n)?2?2?jwjwnH(e)edw d?0所得到的理想的单位脉冲响应hd(n)往往是一个无限长序列，对hd(n)经过适当的加权、截取处理才得到一个所需要的有限长脉冲响应序列。对应不同的加权、截断，就有不同的窗函 数。所要寻找的滤波器脉冲响应就等于理想脉冲响应和窗函数的乘积。 即:h(n)?hd(n)w(n) 由此可见，窗函数的性质就决定了滤波器的品质。例如：窗函数的主瓣宽度决定了滤波器的过渡带宽；窗函数的旁瓣代销决定了滤波器的阻带衰减。以下是几种常见的窗函数： 1、矩形窗 w(n)?RN(n) 2、Hanning窗 w(n)?0.5[1?cos(2?n)]RN(n) N?11 2?n)]RN(n) N?12?n4?n)?0.08cos()]RN(n) 4、Blackman窗 w(n)?[0.42?0.5cos(N?1N?13、Hamming窗 w(n)?[0.54?0.46cos(5、Kaiser窗 w(n)?I0(?1?[(2n(N?1)I0(?))?1]2 其中I0(?)是零阶贝塞尔函数。Kaiser窗可以通过改变?参数，改变其主瓣宽度和旁瓣大小。 在MATLAB 中产生窗函数十分简单： （1）矩形窗（Rectangle Window） 调用格式：w=boxcar(n)，根据长度n 产生一个矩形窗w。 （2）三角窗（Triangular Window） 调用格式：w=triang(n) ，根据长度n 产生一个三角窗w。 （3）汉宁窗（Hanning Window） 调用格式：w=hanning(n) ，根据长度n 产生一个汉宁窗w。 （4）海明窗（Hamming Window） 调用格式：w=hamming(n) ，根据长度n 产生一个海明窗w。 （5）布拉克曼窗（Blackman Window） 调用格式：w=blackman(n) ，根据长度n 产生一个布拉克曼窗w。 （6）恺撒窗（Kaiser Window） 调用格式：w=kaiser(n,beta) ，根据长度n 和影响窗函数旁瓣的β参数产生一个恺撒窗w。 表1 各种窗函数的性能指标 窗函数 窗 过渡带宽 旁瓣峰值（dB） 加窗后过渡带宽 4?/N -13 1.8?/N 加窗后最大衰减（dB） -21 -44 -53 -74 -80 矩形窗 8?/N -31 6.2?/N Hanning窗 8?/N -47 6.6?/N Hammin窗 12?/N -57 11?/N 10?/N 不确定 不确定 BlackmaKaise窗 （三）Matlab提供了基于窗函数法的FIR滤波器的设计函数fir1和fir2 Matlab中，函数fir1（）和fir2（）利用加窗傅里叶级数法设计FIR滤波器。函数fir1（）用来设计传统的LP（低通）、HP（高通）、BP（带通）、BS（带阻）和多频带FIR滤波器；而函数fir2（）用来设计具有任意幅度响应的的FIR滤波器。 1. fi