数字信号处理第六章.ppt

1、第6章 IIR数字滤波器设计,数字滤波器的设计过程一般可归纳为以下四个步骤： (1)按实际要求，确定滤波器的性能指标。 (2)寻找一个稳定因果系统去逼近这个性能指标，即求出数字滤波器的系统函数或单位冲激响应，系统有IIR和FIR两种形式。 (3)采用适当的结构和合适的字长去实现此数字滤波器系统。 (4)验证所设计的系统是否满足给定的性能指标，不满足时对第(2)、(3)步进行修改。,滤波器的性能要求，往往以频率响应幅度特性的容许误差来表征。以低通滤波器为例，频率响应有通带、过渡带及阻带三个范围。在通带内，幅度响应以误差1 逼近于1，即,在阻带中，幅度响应以误差小于2 逼近于零，即,低通滤波器特性

2、,目前，IIR数字滤波器设计主要有以下两种方法： 利用模拟滤波器理论设计数字滤波器，根据模拟滤 波器理论设计出满足要求的模拟滤波器的Ha(s)，然 后再根据Ha(s)求得相应的数字滤波器的H(z)，使H(z) 的频率响应逼近Ha(s)的频率响应，同时H(z)也必须 保持Ha(s)的因果性和稳定性。 利用优化技术设计数字滤波器，在某种最优化准则 意义上逼近所需要的频率响应。这种方法需要进行大 量的迭代运算，必须借助计算机完成。,6.1 模拟低通滤波器的设计方法 模拟滤波器的设计，就是用模拟系统的系统函数去逼近所要求的理想特性。标准的模拟低通滤波器的设计公式有巴特沃思、切比雪夫、贝塞尔、考尔(又称

3、椭圆函数)滤波器等，它们都是根据幅度平方函数来确定的。,6.1.1 幅度平方函数 理想低通滤波器，其模拟理想低通滤波器的幅度响应特性可用幅度平方函数表示,式中，Ha(s)为所设计的模拟滤波器的系统函数，它是s的有理函数；Ha(j)是其稳态响应，|Ha(j)|为滤波器的稳态振幅特性。,(2)虚轴上的零点一定是二阶的。,(1)Ha(j)和Ha(-j)的零极点是镜像对称的。,为保证所设计的滤波器稳定，其极点必须落在S平面左半平面，所以落在S平面左半平面的极点属于 ，落在S右半平面的极点属于 。 零点分布与滤波器相位有关，如果要求滤波器具有最小相位，则应选取S平面左半平面的零点。,6.1.2 巴特沃思

4、低通滤波器设计 巴特沃思低通滤通器的幅度平方函数为,式中，N为正整数，称为滤波器的阶数。c为通带截止频率，当c时,所以c称为3dB衰减点,其幅度平方函数特性如图，不管N为多少，都 通过: 3dB衰减点,当=0时， 当=c时， 当c时，N 越大， 所以N值越大，滤波器通带和阻带的近似性越好。,这种函数具有以下特点：通带内具有最大平坦幅度特性，在正频率范围内随频率升高而单调下降；阶次越高，特性越接近矩形；没有零点。,为了根据幅度平方函数来构造模拟滤波器的系统函数，从每一对极点中选出一个位于S平面左半平面的极点，以保证模拟滤波器的稳定性和因果性。因此，可以直接写出模拟滤波器的系统函数为,其极点为,式

5、中，A0为归一化常数，一般 ，sk为s平面左半平面的极点。,例 求三阶巴特沃斯低通滤波器的系统函数，c=2rad/s。 解 按幅度平方函数求解。因为,则,所有的极点为,选出左半平面的三个极点：,6.1.3 切比雪夫低通滤波器设计 切比雪夫滤波器的幅度特性就在一个频带中(通带或阻带)具有等波纹特性；一种是在通带中是等波纹的，在阻带中是单调的，称为切比雪夫型；一种是在通带内是单调的，在阻带内是等波纹的，称为切比雪夫型。,切比雪夫型滤波器的幅度平方函数为,式中，为小于1的正数，表示通带波动程度，值越大波动也越大；N为正整数，表示滤波器的阶次；/c可以看作以截止频率作为基准频率的归一化频率，TN(x)

6、为切比雪夫多项式：,6.1.4 椭圆滤波器 椭圆滤波器的特点是在通带和阻带的范围内都具有等波动的幅度特性。由于这种滤波器由雅可比椭圆函数来决定，故称之为椭圆滤波器。它的幅度平方函数表示为,JN(x)为N阶雅可比椭圆函数,与巴特沃斯和切比雪夫滤波器相比，由于椭圆滤波器在通带和阻带中都具有等波动特性，所以在相同技术指标要求条件下，它的过渡带最陡，或者说它的阶次可以最低。但是，这种滤波器在参数变化时对特性的影响(灵敏度)也最大。,6.2 冲激响应不变变换法 6.2.1 变换原理 冲激响应不变变换法的基本准则是使数字滤波器的单位冲激响应h(n)等于模拟滤波器的冲激响应ha(t)的抽样值，即,模拟滤波器

7、的系统函数通常是有理函数形式，并且分母的阶次高于分子的阶次。如果Ha(s)仅含有单极点，则,其拉氏反变换：,可见H(z)可由Ha(s)通过下列对应关系得到：, Ha(s)与H(z)的各部分分式的系数相同； 它们的极点根据 实现映射； Ha(s)与H(z)的零点之间没有对应关系。,6.2.2 混叠失真 稳定的模拟滤波器，其全部极点均位于S平面的左半平面内，当其映射到Z平面时，全部极点在单位圆内部，即,只要模拟滤波器是稳定的，用冲激响应不变变换法设计的数字滤波器一定是稳定的。 按冲激响应不变变换法从S平面到Z平面的映射不是单值关系。,S平面上每一条宽为2/T的横带都重叠地映射到 全部z平面上，每一

8、条横带的左半部分映射到z平面 单位圆内，横带右半部分映射到z平面单位圆外。,抽样信号z变换与模拟信号拉氏变换之间关系为： 数字滤波器与模拟滤波器频率响应之间的关系为： 如果模拟滤波器的频率响应限带于抽样频率之内，即 则数字滤波器的频率响应将无失真地重现模拟滤波器 的频率响应,但是任何一个实际的模拟滤波器，其频率响应 不可能真正是限带的，因此会不可避免地产生频谱 混叠。这样，数字滤波器的频率响应就不同于原模 拟滤波器的频率响应而有一定失真。可以说混叠是 冲激响应不变变换法的最大缺点，只有当模拟滤波 器在s/2以上的频率衰减很大时，这个失真才很 小，这时采用冲激响应不变变换法设计的数字滤波 器才能

9、满足精度要求。,在实际应用中冲激响应不变变换法要适当进行修改，当抽样时间间隔很小时，数字滤波器具有很高增益，为了使数字滤波器的增益不随抽样频率而变化，可作如下修改,例6.3 已知模拟滤波器系统函数，用冲激响应不变变换法求数字滤波器系统函数。,6.3 双线性变换法 冲激响应不变变换法有可能产生频谱混叠，原因是s平面与z平面是非单值映射。有必要寻求s平面与z平面的单值映射关系，克服混叠现象。 6.3.1 变换原理 设描述模拟滤波器的微分方程为,方程两边同时取拉氏变换，得到模拟滤波器的系统函数,用梯形法近似逼近积分，得差分方程,式中,对上式两边取z变换，整理得,比较可知s平面与z平面的双线性变换映射

10、关系,6.3.2 逼近情况讨论,s平面左半平面映射到z平面单位圆内；s平面右半平面映射到z平面单位圆外；s平面虚轴映射到z平面单位圆上。,由以上讨论分析可得如下结论： (1)如果模拟滤波器是稳定的，用双线性变换法设计的数字滤波器也一定是稳定的。 (2)双线性变换法中，s平面与z平面是单值映射关系，因此消除了冲激响应不变变换法所存在的混叠问题，但是存在着严重的非线性。 因此，双线性变换法克服了冲激响应不变变换法 所存在的混叠问题，是以引入频率失真为代价的， 当这种失真是可以允许的或者能够得到补偿时，才 能采用双线性变换法设计数字滤波器。,“预畸变”方法是补偿双线性变换中频率非线性关系的有效方法。

11、 首先，进行预畸变。即根据数字滤波器的截止频率求出模拟滤波器的截止频率。,然后，设计低通模拟滤波器，通带截止频率为c、 阻带截止频率为s,把所设计的模拟滤波器截止频率映射成数字频率,这说明，对数字截止频率 进行预畸变， 得到 ；再根据 设计模拟滤波器 ； 经双线性方法映射成的数字滤波器的截止频 率，刚好是 ，即变换前、后一样。,6.3.3 冲激响应不变变换法与双线性变换法的比较 冲激响应不变变换法会造成频响混叠，不宜用来设计高通、带阻滤波器，适用于基本上是限带的滤波器，如低通和带通滤波器。 双线性变换法克服了频率响应的混叠现象，但频率变换关系产生了非线性。对具有分段常数幅频特性的选频滤波器来说

12、，频率非线性失真问题可以用预畸变方法解决。并且双线性变换法同冲激响应不变变换法相比，它具有计算简单和易于实现的特点。因此，实际工作中广发采用双线性变换法来设计IIR数字滤波器。,例6.5 已知某滤波器性能指标为抽样频率 ，通带01.3kHz，阻带2.64kHz，将其转换成数字域频率。 解：通带截止频率,阻带起始频率,总结利用模拟滤波器设计IIR数字低通滤波器的步骤：,1.确定数字低通滤波器的技术指标；,2.转换成模拟低通滤波器的技术指标；,3.按照技术指标设计模拟低通滤波器；,4.转换成低通数字滤波器；,T的 选择,T的 选择,对于脉冲响应不变法：数字滤波器的周期为 ，因此一定满足,对于双线性

13、变换法：不存在混叠现象，T也可以任选,解：（一）用冲激响应不变法设计数字低通滤波器 （1）数字低通滤波器的技术指标为： p = 0.2 rad , p = -1dB s = 0.3 rad , s = -15dB （2）模拟低通滤波器的技术指标为：(T=1s) p = 0.2 rad/s , p = -1dB s = 0.3 rad/s , s = -15dB,例：设计数字低通滤波器，要求在通带内频率低于0.2rad 时，容许幅度误差小于1dB；在频率0.3到之间的阻带误差大 于15dB。指定模拟滤波器采用Butterworth模拟低通滤波器。 试分别用冲激响应不变法和双线性变换法设计滤波器。

14、,（3）设计Butterworth模拟低通滤波器。先计算阶数N及3dB截止频率c。,解上述方程组得N=5.3046，取N=6.将N=6代入上式，求得 。,（4）由N6查表，得到归一化传输函数为：,去归一化（将p = s/c代入）：,（5）由冲激响应不变法将Ha(s)转换成H(z)，得,（二）用双线性变换法设计数字低通滤波器 （1）数字低通滤波器的技术指标为： p = 0.2 rad , p = -1dB s = 0.3 rad , s = -15dB （2）模拟低通滤波器的技术指标为：(T=1s) = (2/T)tan(p/2) , p = 2tan0.1 = 0.65 rad/s , p =

15、 -1dB s = 2tan0.15 = 1.019 rad/s,s = -15dB,（3）设计Butterworth模拟低通滤波器。先计算阶数N及3dB截止频率c。,解上述方程组得N=5.3046，取N=6.将N=6代入上式，求得 。,（4）由N6查表，得到归一化传输函数为：,去归一化（将p = s/c代入）：,（5）由双线性变换将Ha(s)转换成H(z)。得,6.4 频带变换,6.4.1 模拟频带变换 设模拟低通滤波器,低通-低通,低通-高通,低通-带通,低通-带阻,6.4.2 数字频带变换 设原型数字滤波器截止频率为,低通-低通,低通-高通,低通-带通,低通-带阻,6.4.3 频带变换原

16、理 数字频带变换是在数字域进行频带变换,映射关系为,要求该变换必须满足下列条件： (1) 必须是 的有理函数。 (2)要求z平面单位圆内部必须映射到Z平面单位圆内部。 (3)变换前后两个函数要满足一定的频带变换要求。,例6.6 已知某低通滤波器截止频率为50Hz,抽样频率为500Hz，用双线性变换法设计的数字滤波器系统函数如下,求变换出截止频率为200Hz，抽样频率仍为500Hz的数字高通滤波器。,解：由数字低通变高通的公式,例6.7 用模拟频带变换法，由二阶巴特沃思函数设计截止频率为200Hz,抽样频率为500Hz的数字高通滤波器。,无论是经模拟频带变换，还是经数字频带变换，设计同一性能数字滤波器最终结果是相同的。,6.5 IIR数字滤波器的计算机优化设计 以上数字滤波器设计是根据某些近似准则来获得物理可实现的线性滤波系统。但是预期的频率响应与实际的频率响应总存在一定误差。 近代滤波器设计法可使设计的滤波器频率响尽可能逼近于所需频率响应。 常用的方法有最小均方