淮北煤炭师范学院2009届学士学位论文基于MA TLAB的数字滤波器设计学院、专业物理与电子信息学院电子信息科学与技术研究方向基于MATLAB的数字滤波器设计学生姓名耿博学号200513432024指导教师姓名邹锋指导教师职称讲师2009 年4 月18基于MATLAB的数字滤波器设计耿博（淮北煤炭师范学院物理与电子信息学院235000）摘要随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。

数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。

在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。

数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。

对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。

本文根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出了IIR和FIR对语音滤波的效果。

关键词数字滤波器；IIR ；FIR ；MATLABThe Design of Digital Filter based on MATLABGeng BoSchool of Physics and Electronics Information, Huaibei Coal Industry Teachers‟ College, 235000ABSTRACT Along with the information age and the digital world arrival, the digital signal processing has become a now extremely important discipline and the area of technology.The digital signal processing in the correspondence, the multitudinous domains the pronunciation such as the image, the automatic control, the radar, the military, the aerospace, the medical service and the domestic electric appliances and so on have obtained the widespread application.In the digital signal processing application, the digital filter are extremely important and have obtained the widespread application.The digital filter are the digital signal processing important content, the digital filter may divide into IIR and the FIR two main kinds. As for the IIR digital filter design, we need the help of analog prototype filter, and then transform analog filter into digital filter. In the paper we use the design of the pulse response invariable method, the bilinear method and full function design; as for the FIR filter, we can design it directly based on the giving frequency, in the paper it uses the design of the window function.This article according to the IIR filter and the FIR filter characteristic, uses the bilinearity method of transformation under the MATLAB bad boundary to design IIR and to design the FIR numeral filter separately with the window box number, and carries on the analysis to the gathering pronunciation signal, and finally gives IIR and FIR to the pronunciation filter effect.Keywords Digtial Filter;IIR;FIR;MATLAB目次1 引言 (1)2 数字滤波器及MATLAB语言概述 (2)2.1 数字滤波器的定义和分类 (2)2.2 常用滤波器的性能指标 (3)2.3 MATLAB概述 (6)3 IIR滤波器设计 (8)3.1 双线性变换法设计IIR数字滤波器 (8)3.2 脉冲响应不变法 (12)3.3 完全设计函数法 (15)3.4 语音滤波实例 (16)4 FIR滤波器设计 (21)4.1 窗函数法 (21)4.2 FIR滤波器滤波实例 (25)5 总结 (29)参考文献 (30)致谢 (31)1 引言数字滤波在通信、图像编码、语音编码、雷达等许多领域中有着十分广泛的应用。

目前，数字信号滤波器的设计在图像处理、数据压缩等方面的应用取得了令人瞩目的进展和成就。

近年来迅速发展起来。

MATLAB是美国MathWorks公司推出的一套用于工程计算的可视化高性能语言与软件环境。

MATLAB为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。

它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。

MATLAB推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理(signalproeessing)、图像处理(imageproeessing)、小波(wavelet)等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。

数字滤波器与模拟滤波器相比，具有精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配以及能实现模拟滤波器无法进行的特殊滤波等优点]1[。

本文主要介绍有限冲激响应数字滤波器（FIR）和无限冲激响应数字滤波器（IIR）的设计原理、方法、步骤以及在MATLAB中的实现，并以实例形式列出设计程序和仿真结果。

2数字滤波器及MATLAB语言概述2.1数字滤波器的定义和分类数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散时间线性非时变系统，其输入是一组数字量，其输出是经过变换的另一组数字量。

因此，数字滤波器本身既可以是用数字硬件装配成的一台完成给定运算的专用的数字计算机，也可以将所需要的运算编成程序，让通用计算机来执行。

数字滤波器，输入输出均为数字信号，通过一定的运算关系，改变输入信号中所含频率成分的相对比例，或则滤除某些频率成分的器件]2[。

对于数字滤波器而言，若系统函数为H(z),其脉冲响应为h(n),输入时间序列为x(n),则它们在时域内的关系式如下：y(n)=h(n)﹡x(n) （2-1）在Z域内，输入和输出存在如下关系：Y(z)= H(z)X(z) （2-2）式中，X(z)、Y(z)分别为x(n)和y(n)的Z变换。

在频域内，输入和输出则存在如下关系：Y(jω)=H(jω)X(jω) （2-3）式中，H(jω)是数字滤波器的频率特性；X(jω)Y(jω)分别为x(n)和y(n)的频谱，而ω为数字角频率。

数字滤波器可以有很多种分类方法，但总体上可分为两大类。

一类称为经典滤波器，即一般的滤波器，其特点是输入信号中的有用成分和希望滤除的成分占用不同的频带，通过合适的选频滤波器可以实现滤波]3[。

例如，若输入信号中有干扰，信号和干扰的频带互不重叠，则可滤出信号中的干扰得到纯信号。

但是，如果输入信号中信号和干扰的频带相重叠，则干扰就不能被有效的滤出。

另一类称为现代滤波器，如维纳滤波器、卡尔曼滤波器等，其输入信号中有用信号和希望滤除的成分频带重叠。

对于经典滤波器，从频域上也可以分为低通、高通、带通和带阻滤波器。

从时域特性上看，数字滤波器还可以分为有限冲激响应数字滤波器（FIR）和无限冲激响应数字滤波器（IIR）。

对于有限冲激响应数字滤波器（FIR），其输出y(n)只取决于有限个过去和现在的输入，x(n),x(n-1),…,x(n-m),滤波器的输入输出关系可表示为y(n)=∑=-Mrrrnxb)(（2-4）对于无限冲激响应数字滤波器（IIR ），它的输出不仅取决于过去和现在的输入，而且还取决于过去的输出，其差分方程为y(n)+∑=-N k k k n y a 1)(=∑=-Mr r r n x b 0)( (2-5)该差分方程的单位冲激响应是无限延续的。

2.2 常用滤波器的性能指标滤波器性能一般用系统频率特性)(ωj eH 来说明，常用的性能指标主要有以下三个参数：1. 幅度平方函数2*()()*()()()()()j j j j j j z e H e H eH e H e H eH z H z ωωωωωω-==⋅== (2-6)该性能指标主要用来说明系统的幅频特性。

2.相位函数()()Re[()]Im[()]()j j j j j j e H e H e j H e H e e ωωωωωβ=+= (2-7)其中：⎭⎬⎫⎩⎨⎧=)](Re[)](Im[)(ωωωβj j j e H e H arctg e (2-8)该指标主要用来说明系统的相位特性。

3.群延时ωβωτωd ed j )]([)(-= （2-9）定义为相位对角频率导数的负值，说明了滤波器对不同的频率成分的平均延时。

当要求在通带内的群延迟是常数时，滤波器相位响应特性应该是线性的。

实际设计中所能得到的滤波器的频率特性与理想滤波器的频率特性之间存在着一些显著的差别，现以低通滤波器的频率特性为例进行说明。