实验五FIR数字滤波器的设计.docVIP

试验六 FIR数字滤波器设计 一、试验目标 1．熟悉FIR滤波器设计基础方法 2．掌握用窗函数设计FIR数字滤波器原理和方法。 二、试验内容 1．FIR数字滤波器设计方法 FIR滤波器设计问题在于寻求一系统函数，使其频率响应迫近滤波器要求理想频率响应，其对应单位脉冲响应为。 （1）用窗函数设计FIR滤波器基础原理 设计思想：从时域从发，设计迫近理想。设理想滤波器单位脉冲响应为。以低通线性相位FIR数字滤波器为例。 (6-1) 通常是无限长，且是非因果，不能直接作为FIR滤波器单位脉冲响应。要想得到一个因果有限长滤波器h(n)，最直接方法是截断，即截取为有限长因果序列，并用适宜窗函数进行加权作为FIR滤波器单位脉冲响应。根据线性相位滤波器要求，h(n)必需是偶对称。对称中心必需等于滤波器延时常数，即 (6-2) 用矩形窗设计FIR低通滤波器，所设计滤波器幅度函数在通带和阻带全部展现出振荡现象，且最大波纹大约为幅度9%，这个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，通常采取其它类型窗函数。 （2） 经典窗函数 ① 矩形窗(Rectangle Window) （6-3） ② 三角形窗(Bartlett Window) （6-4） ③ 汉宁(Hanning)窗，又称升余弦窗 （6-5） ④ 汉明(Hamming)窗，又称改善升余弦窗 （6-6） ⑤ 布莱克曼(Blankman)窗，又称二阶升余弦窗 （6-7） ⑥ 凯泽(Kaiser)窗 （6-8） 其中：β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间交换关系，通常说来，β越大，过渡带越宽，阻带越小衰减也越大。I0(·)是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为，β确实定可采取下述经验公式： （6-9） 若滤波器通带和阻带波纹相等即δp=δs时，滤波器节数可经过下式确定： （6-10） 式中： （3）利用窗函数设计FIR滤波器具体步骤以下： 1、确定数字滤波器性能要求，临界频率，滤波器单位脉冲响应长度N。 2、依据性能要求，合理选择单位脉冲响应h(n)奇偶对称性，从而确定理想频率响应幅频特征和相频特征。 3、求理想单位脉冲响应，在实际计算中，可对采样，并对其求IDFT，用替换。 4、选择合适窗函数w（n），依据求所需设计FIR滤波器单位脉冲响应。 5、求，分析其幅频特征，若不满足要求，可合适改变窗函数形式或长度N，反复上述设计过程，以得到满意结果。 2．FIR数字滤波器Matlab实现 MATLAB提供相关函数，函数调用格式： b=fir1(n,wn,’ftype’,window) 其中， n——FIR滤波器阶数，对于高通、带阻滤波器n取偶数。 wn——为滤波器截止频率（归一化频率）； ‘ftype’——为滤波器类型；如’high’为高通，’stop’为带阻等； window——窗函数（列向量、其长度为n+1），缺省时，自动取Hamming窗。 MATLAB提供了多个窗函数： wd=boxcar(N)――返回N点矩形窗函数 wd=triang(N)――返回N点三角窗函数 wd=hanning(N)――返回N点汉宁窗函数 wd=hamming(N)――返回N点汉明窗函数 wd=Blackman(N)――返回N点布莱克曼函数 wd=kaiser(N,beta)――返回给定beta值时N点凯泽窗函数 【实例6-1】依据以下技术指示，设计一个数字FIR低通滤波器。 　　　wp=0.2π,ws=0.3π,Rp=0.25dB,Rs=50dB 因为衰减为50dB，可选择窗口有汉明窗和布莱克曼窗。而汉明窗有较小过分带，所以含有较小阶数，所以选择汉明窗。 解：源程序以下： wp=0.2*pi;ws=0.3*pi;deltaw=ws-wp; N0=ceil(6.6*pi/deltaw);% 查表依据汉明窗设计计算所需滤波器h(n)长度， %ceil(x)取大于等于x最小整数 N=N0+mod(N0+1,2);%为实现第一类偶对称滤波器，