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

有限长冲激响应滤波器（FIR）的设计 7.1 线性相位滤波器的特点 7.2 窗函数设计法 7.3 频率抽样设计法 7.4 应用MATLAB设计FIR数字滤波器 有限长单位冲激响应数字滤波器的特点: 有限长单位冲激响应（FIR）可以做成具有严格的线性相位，同时又可以具有任意的幅度特性。 FIR滤波器的单位抽样响应是有限长的，因而FIR滤波器一定是稳定的。 只要经过一定的延时，任何非因果有限长序列都能变成因果的有限长序列，总能用因果系统来实现。 FIR滤波器由于单位冲激响应是有限长的，因而可以用快速傅里叶变换（FFT）算法来实现过滤信号，从而可大大提高运算效率。 但是，要取得很好的衰减特性，FIR滤波器H(z)的阶次比IIR滤波器的要高。 7. 1 线性相位FIR滤波器的特点 7.1.1 线性相位条件 1．I型线性相位滤波器(偶对称，N为奇数) 幅度函数和相位函数的特点： 幅度函数对 偶对称，同时对 也呈偶对称； 相位函数为准确的线性相位。 2．Ⅱ型线性相位滤波器(偶对称，N为偶数) 3．Ⅲ型线性相位滤波器(奇对称，N为奇数) 4．Ⅳ型线性相位滤波器(奇对称，N为偶数) 7. 2 窗函数设计法 7. 2.1 设计方法 以一个截止频率为ωC线性相位的理想矩形幅度特性 的低通滤波器为例。设低通特性的群延时为α，即 依照线性相位滤波器的约束，h(n)必须是偶对称 的，对称中心就为长度的一半(N-1)/2. 因此有 二、加窗处理对频率响应的影响 按照复卷积公式，时域相乘，则在频域上是周期性 卷积关系，即 矩形窗的频率特性 WR(ejω)就是频率内插函数，其幅度函数WR(ω)在 ω＝±2π/N之内为一个主瓣，两侧形成许多衰减振 荡的旁瓣， WR(ω)是周期函数。 理想低通滤波器的频率响应 FIR滤波器的频率响应 加窗处理对理想矩形频率产生的几点影响： 1、使理想频率特性在不连续点处形成一个过渡带，过渡带的宽度等于窗频率响应的主瓣宽度 △ω＝4π/N； 2、在截止频率ω＝ωC±2π/N的地方，H(ω)出现最大的肩峰值，在肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的多少； 3、当截取长度增加时，只会减小过渡带的宽度4π/N,而不会改变肩峰的相对值，这种现象称为吉布斯效应。 7.2.2 各种窗函数 希望窗函数满足两个要求： 1、窗谱主瓣尽可能的窄，以获得较陡的过渡带； 2、尽量减少窗谱的最大旁瓣的相对幅度，也就是使能量尽量集中于主瓣，使肩峰和波纹减小，可增大阻带的衰减。 选用其它形状的窗函数，在边沿处(n=0和n=N-1附近) 的变换比矩形窗要缓慢和平缓，这样可以减小陡峭边 沿所引起的旁瓣分量，使阻带衰减增大； 但是窗谱的主瓣宽度就比矩形窗要宽，这会造成滤波 器幅度函数过渡带的加宽。 1．矩形窗 3．海明（Hamming）窗（又称改进的升余弦窗） 7.2.3 窗函数的设计方法 1、给定要求的频率响应函数Hd(ejω)； 2、对Hd(ejω)求离散傅立叶反变换，得到hd(n); 3、由过渡带宽度及阻带最小衰减的要求，选定窗的形状和N的大小； 4、求得所设计的FIR滤波器的单位抽样响应h(n); 5、对h(n)求傅立叶变换，得到H(ejω)，检验是否满足要求，若不满足，则需要重新设计。 例7.2：设计一个线性相位的FIR滤波器，给定抽样频率为 Ωc＝2π×1.5 ×104(rad/sec),通带截止频率为Ωp＝ 2π×1.5 ×103(rad/sec),阻带起始频率为Ωst＝2π×3 ×104(rad/sec),阻带的衰减不小于－50dB. 解：1、求对应的数字频率。 通带截止频率： 阻带起始频率： 阻带衰减相当于 2、设理想线性相位滤波器Hd(ejω)为 求理想低通滤波器的截止频率 3、确定窗形状及N的大小 由阻带衰减50dB查表，可选海明窗，其阻带最小衰 减－53dB符合要求。 要求的数字频域的过渡带宽度△ω＝0.2π 而海明窗过渡带带宽满足△ω＝6.6π/N 因此 N=6.6π/ △ω=6.6π/ 0.2π=33 τ=(N-1)/2=16 4、由海明窗表达式确定FIR滤波器的h(n) 5、由h(n)确定H(ejω)，再检验指标是否符合要求，如果不满足，则改变N或窗函数的形状来重新计算。 7. 3 频率抽样设计法 频率抽样法是从频域出发，把给定的理想频率响应 加以等间隔抽样，即 7. 4 应用MATLAB设计FIR数字滤波器 7.4.1 与本章内容有关的MATLAB文件 图7.6 滤波器单位抽样响应及其幅频响应曲线 令 （7