实验九-FIR滤波器设计.pptVIP

二、窗函数设计法 1、设计方法 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 取矩形窗： 加窗处理后对频率响应的影响： 时域乘积相当于频域卷积 理想滤波器的频率响应： 加窗函数的影响： 不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。 加窗函数的影响： 2、各种窗函数 窗函数的要求： 3、窗函数法的设计步骤 例：用窗函数法设计线性相位FIR低通数字滤波器，通带截止频率Wp=0.5π，阻带截止频率Ws=0.6π，阻带衰减不小于40dB，通带衰减不大于3dB. 根据阻带衰减，采用汉宁窗设计滤波器。MATLAB程序： Wp=0.5*pi; Ws=0.6*pi; wdel=Ws-Wp;%过度带宽度 N=ceil(8*pi/wdel);%计算滤波器长度 Wn=(Wp+Ws)/2;%截止频率 window=hanning(N+1); b=fir1(N,Wn/pi,window); freqz(b,1,512) figure；n=0:N;stem(n,b) * * w(n)：窗函数序列 要选择合适的形状和长度 其理想单位抽样响应： 中心点为 α 的偶对称无限长非因果序列 则FIR滤波器的单位抽样响应： 按第一类线性相位条件，得 而矩形窗的频率响应： 矩形窗的幅度函数 N p 2 N p 2 - N p 4 N p 4 - N p p - 主瓣 旁瓣 其幅度函数： 则FIR滤波器的频率响应： 在 处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少 改变N只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例。其相对比例由窗函数形状决定，称为Gibbs效应 窗谱主瓣尽可能窄，以获得较陡的过渡带 尽量减少窗谱最大旁瓣的相对幅度 以减小肩峰和波纹 但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 -21 -25 -44 -53 -74 -80 0.9 2.1 3.1 3.3 5.5 5 2 4 4 4 6 -13 -25 -31 -41 -57 -57 矩形窗 三角形窗 汉宁窗 海明窗 布拉克曼窗 凯泽窗 阻带最小衰减 /dB 过渡带宽 主瓣宽度 旁瓣峰值 /dB 加窗后滤波器性能指标 窗谱性能指标 窗函数 阻带最小衰减只由窗形状决定 过渡带宽则与窗形状和窗宽N都有关 给定理想的频率响应函数 及技术指标 求出理想的单位抽样响应 根据阻带衰减选择窗函数 计算频率响应 ，验算指标是否满足要求 根据过渡带宽度确定N值 求所设计的FIR滤波器的单位抽样响应 MATLAB工具箱提供了采用窗函数法设计FIR滤波器的函数fir1 其调用格式为： b=fir1(n,Wn) b=fir1(n,Wn,’ftype’) b=fir1(n,Wn,window) b=fir1(n,Wn,’ftype’,window) 说明：n 表示滤波器的阶数， Wn 为滤波器截止频率，取值范围0～1; ‘ftype’ 参数代表所设计的滤波器类型：‘high’——高通， ‘stop’——带阻滤波器，默认时为低通和带通滤波器； ‘window’ 为窗函数，是长度为n+1的列向量，默认时窗函数自动取Hamming窗。 MATLAB工具箱提供了采用窗函数法设计FIR滤波器的函数fir1 ‘window’ 为窗函数类型： Bartlett : 巴特利特 窗 Blackman ：布莱克曼窗 Boxcar ：矩形窗 Chebwin ：契比雪夫窗 Hamming ：海明窗 Hanning ：汉宁窗 Kaiser ：凯泽窗 triang ：三角窗