数字滤波器研究课程论文.doc

研究生课程论文

(2010-2011学年第一学期)

数字滤波器的研究

研究生：

提交日期：研究生签名：

学号

学院

电子与信息学院

课程编号

S0809074

课程名称

现代数字信号处理

学位类别

硕士

任课教师

教师评语：

成绩评定：分任课教师签名：年月日

PAGE2

PAGE1

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

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

IIR滤波器的结构与设计

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

（1）

此结构便于准确地实现滤波器的零点、极点，也便于调整滤波器的频率响应性能[7]；另一个优点是所需的存储单元较少，在硬件实现时甚至还可以用一个二阶节进行时分复用，进一步降低了对现场可编程逻辑阵列硬件资源的要求[8]。IIR滤波器的系数在做强噪音下的语音增强，需要通过共振峰对语音进行端点识别。前两个共振峰对区别不同语音有非常重要的作用[9]；1kHz以下基本为第一共振峰的范围,对语音感知、语意识别作用比较重要的第二共振峰基本在1kHz之外[9]。因此对带噪语音信号加入预处理环节,即进行数字滤波(通带下限1kHz,Butterworth滤波器高通,阻带衰减3dB，语音信号采样频率为8000Hz)。利用MATLAB提供的FDATool便可直接得到直接II型各个子系由于硬件当中对小数进行运算耗费资源比较大，需要在计算精度与速度之间优化选择[10]。

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

设计的具体步骤：

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

2：双线性变换公式：（2）

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

用公式（2）预畸所求数字滤波器的的数字频率指标（数字滤波器通带截止频率0.2π，阻带截止频率0.6π），得到一个等价的模拟低通滤波器的频率指标。

3：得到原型模拟低通滤波器的指标后作频率响应，从而得模拟的传输函数H(s)，然后对传输函数进行双线性变换，从而得到数字滤波器的传输函数H(z).

结果分析：

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

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

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

2.3FIR滤波器的结构与设计

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

利用傅里叶加窗级数法设计线性相位低通FIR滤波器（没有使用到fir1函数）

设计指标：通带截止频率4rad/s，阻带截止频率6rad/s，最大通带衰减为0.2dB,最小阻带衰减为42dB，抽样率为18rad/s.

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

图6.各种窗函数设计的低通滤波器

结果分析：

Hamming窗函数的长度是30

通带波纹是Rp=0.0394dB

最小阻带衰减时As=52dB

Hanning窗函数的长度是29

通带波纹是Rp=0.0711dB

最小阻带衰减时As=44dB

Blackman窗函数的长度是52

通带波纹是Rp=0.0027dB

最小阻带衰减时As=75dB

均满足滤波器的性能指标。

结论

在电子系统中，具有严格的线性相位特性的FIR数字滤波器被广泛使用，是信号处理的基本构成之一。利用matlab的强大运算功能，基于matlab信号处理工具箱（signalprocessingtoolbox）的数字滤波器设计法可以快速有效的设计由软件组成的常规数字滤波器，设计方便、快捷，极大的减轻了工作量。在设计过程中可以对比滤波器特性，随时更改程序参数，以达到滤波器设计的最优化。利用matlab设计数字滤波器在信号处理软件和微机保护中，有着广泛的应用前景。本文通过调用simulink中的功能模