数字滤波器研究课程论文

研 究 生 课 程 论 文

(2010-2011学年第一学期)

数字滤波器的研究

研究生：

中文题目：数字滤波器的研究

研究生姓名

摘要：无限冲激响应（IIR）数字滤波器和有限冲激响应（FIR）数字滤波器广泛应用于数字信号处理系统中。FIR 滤波器具有很好的线性相位特性，因此线性相位特性在工程实际中如数据通信、图像处理等应用领域具有比较重要的意义。IIR 数字滤波器具结构简单和实现容易的特点, 这对于相位要求不是特别严格的场合有很重要的意义。利用MATLAB 实现IIR 低通数字滤波器的设计,并对被干扰的信号进行滤波,达到了预期结果。

关键词: IIR；FIR；线性相位；数字滤波器；MATLAB

1 引言

数字滤波器是现代数字信号处理的不可或缺的一部分，数字滤波器按其响应形式可分为无限脉冲响应（IIR）滤波器和有限脉冲响应（FIR）滤波器和两种。当系统无严格相位要求时，IIR 数字滤波器与FIR 数字滤波器相比，可用较低的阶数获得较高的选择性[1]。对于相同的设计指标，FIR 滤波器所要求的阶数比IIR滤波器高5～10 倍，而且信号的延迟也较大[2]。MATLAB 作为一种矩阵运算为基础的交互式程序语言，是进行科学研究常用且必不可少的工具。MATLAB着重针对科学计算、工程计算和绘图的需求。它用解释方式工作，键入程序立即得出结果，人机交互性能好，为科学人员所乐于接受。MATLAB 提供了丰富的函数，其中fir1 函数实现了加窗线性相位FIR 滤波器设计的经典方法，fir1 主要用于常用的标准通带滤波器设计，包括低通、带通、高通和带阻数字滤波器[3]。一般来讲，数字滤波器的设计要经过三步：确定设计指标、模拟逼近和数字转换。通常在设计滤波器之前，先根据具体的应用确定一些技术指标，

然后就可以根据数学知识和滤波器的基本原理提出一个滤波器的模型来逼近给定的指标，逼近的结果通常是得到以差分方程或脉冲响应描述的滤波器，最后可以根据这个描述用硬件或软件实现，至此完成一个滤波器设计的全过程[3]。

2.数字滤波器的设计

2.1模拟椭圆滤波器对加噪语音信号的滤波：

设计模拟滤波器的一般方法是：首先根据技术指标确定滤波器的传输函数H(s）,然后综合电路网络实现该传递函数。而设计滤波器传输函数的关键是找到逼近函数，在各种滤波器响应中，椭圆函数响应最优越[4-6]。

首先，打开的是一个摩托车引擎信号的波形文件，它是一个频率范围大致分布在5000-6000Hz的一个低频语音信号。对于该语音信号，施加高斯白噪声，设计一个1dB截止频率6000Hz 最小阻带衰减为100dB的椭圆低通滤波器。其原始信号波形、加噪声后信号波形与滤波后信号波形分别图1和图2：

图1.加噪声后语音信号波形与频谱

图2.原始语音信号与滤波后语音信号波形与频谱对比

通过matlab仿真可以发现，原始语音信号基本可以从高斯白噪声中恢复，但是滤波后的语音信号比原始语音信号有明显的波形延迟，这是模拟椭圆低通滤波器不足的地方。

2.2IIR滤波器的结构与设计

这里采用的是直接II型结构实现IIR 滤波器，如式（1）所示：

（1）

此结构便于准确地实现滤波器的零点、极点，也便于调整滤波器的频率响应性能[7]；另一个优点是所需的存储单元较少，在硬件实现时甚至还可以用一个二阶节进行时分复用，进一步降低了对现场可编程逻辑阵列硬件资源的要求[8]。IIR滤波器的系数在做强噪音下的语音增强，需要通过共振峰对语音进行

端点识别。前两个共振峰对区别不同语音有非常重要的作用[9]；1 kHz以下基本为第一共振峰的范围,对语音感知、语意识别作用比较重要的第二共振峰基本在1 kHz之外[9]。因此对带噪语音信号加入预处理环节,即进行数字滤波(通带下限1 kHz,Butterworth滤波器高通,阻带衰减3 dB，语音信号采样频率为8000 Hz)。利用MATLAB提供的FDATool便可直接得到直接II型各个子系由于硬件当中对小数进行运算耗费资源比较大，需要在计算精度与速度之间优化选择[10]。

下面用双线性变换设计一个0.4dB截止频率为10KHz且在30KHz处有最小阻带衰减为50dB的数字巴特沃兹低通滤波器，其抽样率为100KHz；

设计的具体步骤：

1：先设计一个符合以上指标的模拟巴特沃兹低通滤波器，作为与要设计的数字滤波器的参照滤波器。先用手算计算出滤波器所需的阶数是Ns=7, 与[Ns,wn]=buttord(wp,ws,Rp,Rs,'s');得出的结果一样。

2：双线性变换公式：

)

2

tan(

2ω

T

=

Ω（2）

其中，Ω是模拟的频率，ω是数字的归一化频率，T是积分步长。

用公式（2）预畸所求数字滤波器的的数字频率指标（数字滤波器通带截止频率0.2π，阻带截止频率0.6π），得到一个等价的模拟低通滤波器的频率指标。3：得到原型模拟低通滤波器的指标后作频率响应，从而得模拟的传输函数H(s)，然后对传输函数进行双线性变换，从而得到数字滤波器的传输函数H(z).

结果分析：

图3.模拟原型与数字滤波器幅度相应比较

图4. 模拟原型与数字滤波器相位相应比较

由图1可以看出模拟原型低通滤波器的通带截止频率为10KHz对应着数字低通滤波器的截止频率0.2π rad与双线性变换结果一致。由图2可以看出，相位响应就有很大差别了，这是由双线性变换的频率畸变所致。

2.3 FIR滤波器的结构与设计

对于IIR滤波器必须明确要求推导出来的传输函数是稳定的，而另一方面，对于FIR滤波器的设计，由于滤波器的传输函数是以1 /z的多项式表示的，所以FIR数字滤波器都是稳定的[11]。

利用傅里叶加窗级数法设计线性相位低通FIR滤波器（没有使用到fir1函数）设计指标：通带截止频率4rad/s，阻带截止频率6rad/s，最大通带衰减为0.2dB,最小阻带衰减为42dB，抽样率为18rad/s.

图5.由截短得到的理想FIR低通滤波器的冲击响应