DSP课程设计说明书 基于MATLAB的语音信号分析及滤波文档格式.docx

1、按照滤波器的功能数字滤波器可以分为低通滤波器（LPF）,高通滤波器（HPF）,带通滤波器（BPF）和带阻滤波器（BSF）。按照单位取样响应 h（n）的时域特性经典数字滤波器可分为无限冲激响应（IIR，InfiniteImpulseResponse）滤波器和有限冲激响应（FIR，FiniteImpulseResponse）滤波器。如果单位取样响应是时宽无限的 h（n）， no n 则称之为IIR 滤波器；而如果单位取样响应是时宽有限的 h（n）， n1 n n2 ，则称之为FIR 滤波器。这两类滤波器无论是在性能上还是在设计方法上都有着很大的区别。FIR 滤波器可以对给定的频率特性直接进行设计，

2、而 IIR 滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法来进行设计。2.2 IIR 滤波器的设计原理利用模拟滤波器来设计 IIR 数字滤波器是最常见的方法，其设计过程是按照技术要求设计模拟滤波器，得到其系统函数Ha（s） ，然后经过一定的变换，得到满足要求的数字滤波器的H（z） 。这种设计方法的优点是模拟滤波器的设计理论很成熟，有很多现成的公式、表格、图形和曲线可以直接使用，从而使 IIR 数字滤波器的设计变得很简单。常用的 IIR 滤波器设计方法主要有以下几种：1. 以模拟滤波器函数为基础的变换法：即先设计一满足指定条件的模拟滤波器 H（s），再将该模拟滤波器转化为数字滤波器

3、 H（z）。2. 直接设计法：在 z 平面内，根据零、极点对系统特性的影响，调整零极点位置得 H（z）。3. 最优化设计法（计算机辅助设计）：在某种最小化误差准则下，建立差分方程系数ak 、bk 对理想特性的逼近方程，使用迭代方法解方程组得到最佳逼近系统。由于此方法计算量大，需要借助于计算机进行设计。在此主要介绍由模拟滤波器转换为数字滤波器的设计方法。2.2.1 IIR 模拟低通滤波器的设计为了从模拟滤波器设计 IIR 数字滤波器，必须先设计一个满足技术指标的模拟原型滤波器。设计“模拟原型”滤波器有多种方法，如模拟低通逼近有巴特沃斯（Butter worth）型，切比雪夫（Chebyshev）

4、型或椭圆（Elliptic）型。低通滤波器是最基本的，至于高通、带通、带阻等滤波器可以用频率变换的方法由低通滤波器变换得到。模拟滤波器设计就是将一组规定的设计要求，转化为相应的模拟系统函数Ha （s） ，使其逼近某个理想滤波器的特性。这里主要介绍巴特沃斯（Butterworth）低通滤波器的设计原理。巴特沃斯低通滤波器的幅度平方函数定义为：| H （ jl） |2 =11+ C 2l2 N（2-1）其中 C 为一常数参数，N 为滤波器阶数，l为归一化低通截止频率,l= W / W p 。式中 N 为整数，是滤波器的阶次。巴特沃斯低通滤波器在通带内具有最大平坦的振幅特性，这就是说，N 阶低通滤波

5、器在W = 0 处幅度平方函数的前 2N-1 阶导数等于零，在阻带内的逼近是单调变化的。巴特沃斯低通滤波器的振幅特性如图 2-1 所示。H （ jW） 210.707N=2N=4N=10W pW图 2-1 巴特沃斯低通滤波器的振幅特性滤波器的特性完全由其阶数 N 决定。当 N 增加时，滤波器的特性曲线变得更陡峭，这时虽然由（2-1）式决定了在W = W p 处的幅度函数总是衰减 3dB，但是它们将在通带的更大范围内接近于 1，在阻带内更迅速的接近于零，因而振幅特性更接近于理想的矩形频率特性。滤波器的振幅特性对参数 N 的依赖关系如图 2-1 所示。设归一化巴特沃斯低通滤波器的归一化频率为l，归

6、一化传递函数为 H （ p） ，其中 p =jl，则由（2-1）式得：H （ jl） 2l= pj=11+ C 2 （-1）N p2 N由于Ha （s）Ha（-s） = A2 （W）W=- js=1+ （s jWc）2 N（2-2）所以巴特沃斯滤波器属于全极点滤波器。1. 常用设计巴特沃斯低通滤波器指标lp ：通带截止频率；ap ：通带衰减，单位：dB；ls ：阻带起始频率；as ：阻带衰减，单位：dB。H （ jW）apasW pWsW图 2-2 巴特沃斯低通滤波器指标说明：（1） 衰减在这里以分贝（dB）为单位；即a= 10 lg= 10 lg 1+ C 2l2 N （2） 当a= 3dB

7、 时W p = WC 为通常意义上的截止频率。（3） 在滤波器设计中常选用归一化的频率l= W / WC ，即l = W p = 1 ,l = WspWsWpp2. 巴特沃斯低通滤波器设计实质根据设计指标要求lp ，ap ，ls ，as 确定归一化巴特沃斯低通滤波器幅度平方函数中的待定系数 C 及滤波器的阶数 N；然后再根据幅度平方函数确定巴特沃斯低通滤波器的传递函数 H（s）。（1） 将实际频率W 归一化得l = W p= 1 ，l = Ws ，再根据已知的apa| H （ jl） |2 =1sp ，s ，幅度平方函数（2） 求 C 和 N1+ C 2l2 N确定 C 和 N。由a（l） =

8、 10 lg1= 10 lg （1+ C 2l2 N ）并带入 l ，a ，l，a 得p10 lg （1+ C 2l2 N ）= apss10 lg （1+ C 2l2 N ）=assC 2l2 Na= 1010 -110ps因为lp = 1 ，所以C 2l2 N =as 10 -1C 2 = 1010 -1asas 由l2 N = 1010 -1 = 1010 -1 两边取对数得：sC 2ap 1010 -1N = lg algls其中a =as ap a这样可以求出 C 和 N。注意：当a = 3dB 时， C 2 = 1010 -1 = 100.3 -1 = 1，即 C=1，此时巴特沃斯

9、滤波器只剩下一个参数 N。（3） 确定巴特沃斯滤波器的传递函数 H（p）。H （ p）H （- p） = G（ jl） 2l= p =1+（ p ）2 N=11+ （-1）N p2 N由1+ （-1）N p2 N = 0 ，解得极点为：2 Nkp = e j 2 k +N -1p,k = 1, 2,L, 2N将 p 左半平面的极点赋予 H （ p） 即H （ p） =1其中 pk= e j 2 k +N -1p,（p - p1 ）（p - p2 ）L（p - pN ）k = 1, 2,L, N为了便于设计，工程上已将当lp = 1 时，各阶巴特沃斯低通滤波器系统函数设计成表格供查阅，该表如表

10、2-1 所示。在表 2-1 中的函数被称为归一化巴特沃斯原型低通滤波器系统函数。阶次归一化系统函数s + 12s 2 +2s + 13s3 + 2s 2 + 2s + 14s4 + 2.6s3 + 3.4s2 + 2.6s +15s5 + 3.2361s4 + 5.2361s3 + 5.2361s2 + 3.2361s +1表 2-1 归一化巴特沃斯模拟低通滤波器系统函数表（4） 去掉归一化影响上面设计中采用归一化的频率即lp = 1 ，而实际中截止频率为W p ，所以要进行如下的变量代换：p = jl= j W = sW pW p即H （s） = H （ p）p= s 综上，归纳出设计巴特沃斯低通滤波器的方法如下：（1） 计算归一化频率l = W p = 1 ，l = Ws 。psa p （2） 根据设计要求按照C 2 = 1010 -1和 N =lg a其中a =计算巴特沃斯滤波器的参数 C 和阶次 N；注意当ap = 3dB 时 C=1。（3） 利用 N 查表获得归一化巴特沃斯低通原型滤波器的系统函数 H （ p） ；（4） 令 H （ p） 中的 p = s 得到截止频率为W的巴特沃斯低通滤波器的系统函数。2.2.2 IIR 数字低通滤波器的设计（双线性变换法）从模拟滤波器设计 IIR 数字滤波器就是要由模拟滤波器的系统函数 H