利用MATLAB和窗函数法设计FIR数字滤波器.pdfVIP

2007年第9期 福建电脑 177 黎民红 (南半大学湖南街阳421001) 【摘要J：传统的量字溥盘蓦_陵计方法繁琐且蛄果不直现．奉文利用MAlI仙鼻有强大的科学计算和固形显示遗一优 点，与宙矗敷法谩计理诗相嫱合共同设计nR敷字滤连嚣，不但使设计结果曼加直观，而且提高了溥凌器的谩计精度，^L而更 好地达到顸期挂秉。 I关键词J：窗函敷法；MA肛JAB；nR滤波嚣 l_引官 目莞塑 在数字信号处理中．数字滤波器十分重要并已获得广泛应 用，效字滤渡器与模拟滤渡器比较，具有精度高、稳定、体积小、 Ⅲ十口 重量轻、灵活、不要求阻抗匹配以及实现模拟滤波器无法宴现的 机#g 攀蔓～ 特殊滤波功能等优点．使数字滤波器成为信号处理领域的一十 重要研究课题【l】。现代数字滤波器可阻用软件或硬件两种方式 为了满足工程上需要．可以通过改变窗函数的形状来改善 米实现．软件方式实现的优点是可以通过滤波器参数的改变去 滤渡器的幅度频率特性。宙函数的选择原则是： 调整滤波器的性能。 (1)具有较低的旁瓣幅度，尤其是第一旁瓣的幅度； MAll^B是具有很强的科学计算和图形显示功能的软件系 圆旁瓣的幅度下降的速率要快．有利于增加阻带衰碱； 统．可以对数字滤波器进行精确设计．并且方便地进行FPr频谱 (3)主瓣的宽度要窄．这样可以得到较窄的速度带。 分析与频谱田显示。因而．恢宴现具有编程简单．对设备硬件要求 通常以上几种难以同时满足．实际中选用的窗函数往往是 不高．易于实现．能快速地在普通Pc机上完成信号的频谱分析和 它们的折中。 频谱图显示等特点．能够对所计的数字滤渡器进行快速地检验 2．3宙函散{击的设计步骤 和分析。 用窗函数法设计nR滤波器的步骤如下^ 2．F1R滤波器的设忡 (1)根据过渡带宽及阻带衰碱要求避择窗函数的类型井估 Z1FTR游渡器简介 计窗口长度Nf或阶效M=N一1k宙函数类型可根据最小阻带衰 数字滤波器是一个能够完成特定任务的离散时间系统．它 减AS独立选择．因为窗口长度N对最小阻带衰减^s没有影 可以利用有限精度算法来实现。当采用一个因果稳定的离散线 响。在N和窗函数类型确定后。即可调用MA|nJ^B中的窗函数 性移不变系统的系鲈函数去逼近滤渡器所要求的性能指标时． 求出窗函数。 由于系统函数有无限长单位冲激响应ⅡIR1系统函数和有限长单 ∞根据待隶滤渡器的理想频率响应求出理想单位脉冲响 位冲激响应∞R)系统函数两种，相应地数字滤渡器也就有无限 应h和)。如果给出待求滤波器的频率响应为凰㈣则理想的单 长单位冲激响应0IR)滤渡器和有限长单位冲激响应∞R】滤渡器位脉冲响应可以用下面的傅里叶反变换式球出： 两种剐。 1_F 玎R滤渡器的设计问题实质上是确定能满足所要求的转移 hd(n)2玄上，Hd(e归)e’啊m f31 序列或脉冲响应的常数的问题．设计方法主要有窗函数法、频率 在一般情况下^肛犀不能用封闭公式表示的，需要采用数 采样法和等波纹最佳逼近法等。nR滤波嚣的基本结构可虬理 值方法表示。从删到∞珈采样N点．采用离散傅垦叶反变按 解为一个分节的延时线．把母一节的输出加权累加．可得到滤波 皿兀’即可求出。 嚣的输出。HR滤波器的冲激响应h(n)是有限长的，数学上M阶0)计算滤波器的单位脉冲响应h(n)。它是理想单位脉冲响 F1R滤渡器可以表示为m： 应和窗函数的乘积．即 幽 h(Ⅱ)=h ’ (4)