基于Matlab的FIR带通滤波器设计与实现

摘要：借助Matlab的FDATOOL滤波器设计分析软件，设计了一种FIR数字带通滤波器，并对一段含噪语音信号进行滤波。利用汇编语言编程，在DSP上实现了该滤波器。实验结果表明，该数字带通滤波器精确，稳定性好，易于移植，具有很强的实用性与灵活性。

　　关键词：FIR数字带通滤波器；语音去噪；DSP；Matlab

　　引言

　　随着电子技术的飞速发展，人们正逐步进入数字化时代，数字滤波器越来越受到人们的关注，并且在近代电信设备和各类控制系统中的应用极为广泛，如语音处理、图像处理、通信、电视、雷达、生物医学信号处理等。数字滤波器根据其冲激响应函数的时域特性，可分为有无限长单位冲激响应（IIR）滤波器和有限长单位冲击响应（FIR）滤波器。IIR滤波器虽然可以利用模拟滤波器设计的结果，而模拟滤波器的设计有大量的图表可查，方便简单，但它的相位是非线性的。图像处理以及数据传输都要求信道具有线性相位，同时又可以具有任意的幅度特性。此外，FIR滤波器的单位冲激响应是有限长的，因而滤波器一定是稳定的。随着Matlab软件和信号处理工具箱的不断完善，可以利用Matl ab信号处理工具箱快速有效地实现数字滤波器的设计、分析和仿真。本文首先使用Matlab设计了一个FIR数字带通滤波器，并对一段含噪语音信号进行滤波，然后在TI公司生产的TMS320C5402芯片上具体实现了该滤波器。

　　1 FIR数字带通滤波器设计原理

　　设FIR滤波器单位脉冲响应h（n）长度为N，其系统函数H（z）为：

　　FIR滤波器的设计任务是选择有限长度的h（n），使传输函数满足技术要求。一般是先给定所要求的理想滤波器频率响应，导出hd（n）。因为理想滤波器的冲激响应hd（n）是无限长的非因果序列，而所要设计的是hd（n）是有限长的FIR滤波器，所以要用有限长序列hd（n）来逼近无限长序列h（n）。设：

　　由此可见，窗函数不仅影响原信号在时域内的波形，而且也影响频域内的波形。

　　2 用Matlab辅助设计与仿真

　　设计一个FIR滤波器，关键是要得到正确的h（n）系数。把生成的滤波器系数加载到目标DSP程序中有两种方式：

　　（1）把滤波器系数输人到一个C语言的头文件，在所建工程中添加该头文件；

　　（2）直接把生成的滤波器系数加载到DSP程序的一个变量中。

本文采用第二种方法。选择菜单File→Export打开Export对话框，选择Coefficient File（ASCII），指定输出文件名，生成*fcf文件