《测试信号分析与处理第2版宋爱国刘文波王爱民电子课件第7章数字滤波器的设计》由会员分享,可在线阅读,更多相关《测试信号分析与处理第2版宋爱国刘文波王爱民电子课件第7章数字滤波器的设计(175页珍藏版)》请在金锄头文库上搜索。
1、,第七章 数字滤波器的设计,概述,7.1,IIR数字滤波器设计,7.2,FIR数字滤波器的设计,7.3,智能仪器中常用的数字滤波算法简介,7.4,,内容提要,本章主要介绍IIR数字滤波器和FIR数字滤波器的常用设计方法,简单介绍智能仪器仪表中常用的数字滤波算法。,,第一节 概述,一个数字滤波器可以用一个N阶差分方程来描述,即,(7-1),或者以它的系统函数H(z)来描述,(7-2),,第一节 概述,数字滤波器的设计大致包括以下4个步骤: 1)根据任务需要,确定数字滤波器应达到的性能指标,如通带截止频率 、阻带截止频率 、通带起伏 等。此外还必须确定采样周期T或采样频率 。 2)确定数字滤波器的
2、系统函数H(z)或h(n),使其频率特性满足技术指标要求。,,第一节 概述,3)用一个有限精度的运算去实现H(z)或h(n),包括选择合理的网络结构、恰当的有效字长,以及有效数字的处理方法等。 4)确定工程实现方法。用实际数字系统(通用计算机软件或专用数字滤波器硬件)实现H(z)或h(n)。,,第二节 IIR数字滤波器设计,设计IIR数字滤波器的方法可以归纳为两类: (一)模拟数字转换法 (二)是直接设计法,,利用模拟滤波器成熟的理论和方法来设计IIR数字低通滤波器的设计过程是:按照技术要求设计一个模拟低通滤波器,得到模拟低通滤波器的传递函数 ,再按一定的转换关系将 转换成数字低通滤波器的系统
3、函数 。这样设计的关键问题就是寻找这种转换关系,将S平面上的 转换成Z平面上的 。为了保证转换后的 稳定且满足技术要求,对转换关系提出两点要求:,,1)因果稳定的模拟滤波器转换成数字滤波器,仍是因果稳定的。模拟滤波器因果稳定要求其传递函数 的极点全部位于S平面的左半平面;数字滤波器因果稳定则要求 的极点全部在单位圆内。因此,转换关系应是S平面的左半平面映射至Z平面的单位圆内部。 2)数字滤波器的频率响应模仿模拟滤波器的频率响应,S平面的虚轴映射为Z平面的单位圆,相应的频率之间成线性关系。 将传递函数 从S平面转换到Z平面的方法有多种,工程上常用的是冲激响应不变法和双线性变换法。,,主要内容,,
4、一、用冲激响应不变法设计IIR数字低通滤波器,冲激响应不变法是使数字滤波器的单位冲激响应h(n)等于模拟滤波器的单位冲激响应的等间隔采样。即,(7-3),式中,T为采样间隔。 设模拟滤波器 只有单阶极点(若有多重极点,则求拉氏反变换会复杂一些),且分母多项式的阶次高于分子多项式的阶次,将 用部分分式表示,(7-4),,一、用冲激响应不变法设计IIR数字低通滤波器,其中, 为 的单阶极点。将 进行拉氏反变换得到,式中, 是单位阶跃函数。对 进行等间隔采样,采样间隔为T,得到,对上式进行Z变换,得到数字滤波器的系统函数H(z),(7-5),,一、用冲激响应不变法设计IIR数字低通滤波器,对比式(7
5、-4)与式(7-5), 在S平面上的极点 ,根据 的关系映射为 在Z平面上的极点 ,系数 不变化。 模拟信号 的傅里叶变换 和其采样信号 的傅里叶变换 之间的关系,(7-6),,一、用冲激响应不变法设计IIR数字低通滤波器,将 代入上式,得,(7-7),因为数字滤波器的系统函数是由采样信号 得来的,因此有,(7-8),上式表明将模拟信号 的拉氏变换在S平面上沿虚轴按照周期 延拓后,再映射到Z平面上,就得H(z)。,,一、用冲激响应不变法设计IIR数字低通滤波器,设,即,(7-9),由第三章的知识可知,S平面的虚轴( )映射成Z平面的单位圆(r1),S平面左半平面( )映射到Z平面的单位圆内(
6、),S平面右半平面( )映射到Z平面的单位圆外( )。这说明如果 因果稳定,转换后得到的H(z)仍是因果稳定的。,,一、用冲激响应不变法设计IIR数字低通滤波器,另外,注意到 是一个周期函数,可写成,,M为任意整数,当模拟频率 从 变化到 时,数字频率 则从 变化到 ,且按照式(7-9), ,即 与 之间成线性关系。 但是,从模拟信号 到采样信号 ,其拉氏变换要按照式(7-7),以 为周期,沿虚轴方向进行周期化。如果原模拟信号 的频带不是限于 之间,则会在 的奇数倍附近产生频率混叠,从而映射到Z平面,在 附近产生频率混叠。冲激响应不变法的频率混叠现象如图7-1所示。,,一、用冲激响应不变法设计
7、IIR数字低通滤波器,这种频率混叠现象会使设计出的数字滤波器在 附近的频率特性,程度不同地偏离模拟滤波器在附近的频率特性,严重时使数字滤波器不满足给定的技术指标。为此,希望设计的滤波器是带限滤波器,如果不是带限的,例如高通滤波器、带阻滤波器,需要在高通带阻滤波器之前加保护滤波器,滤除高于折叠频率 以上的频带,以免产生频率混叠现象。但这样会增加系统的成本和复杂性,因此,高通与带阻滤波器不适合用这种方法设计。,图7-1 冲激响应不变法的频率混叠现象,,一、用冲激响应不变法设计IIR数字低通滤波器,假设 没有频率混叠现象,即满足,0,,按照式(7-8),并将关系式 , 代入,得,说明用冲激响应不变法
8、设计的数字滤波器可以很好地重现原模拟滤波器的频响。上式中, 的幅频特性与采样间隔成反比,这样当T较小时, 就会有太高的增益。为避免这一现象,令 ,那么,,一、用冲激响应不变法设计IIR数字低通滤波器,此时, 综上所述,冲激响应不变法的优点是频率坐标变换是线性的,即 ,如果不考虑频率混叠现象,用这种方法设计的数字滤波器会很好地重现原模拟滤波器的频率特性。另一个优点是数字滤波器的单位冲激响应完全模仿模拟滤波器的单位冲激响应,时域特性逼近好。缺点是会产生频率混叠现象,适合低通、带通滤波器的设计,不适合高通、带阻滤波器的设计。,,一、用冲激响应不变法设计IIR数字低通滤波器,例7-1 已知模拟滤波器的
9、传递函数 为 用冲激响应不变法将 转换成数字滤波器的系统函数H(z)。 解:首先将写成部分分式 极点为 , 那么H(z)的极点为 ,,,一、用冲激响应不变法设计IIR数字低通滤波器,按照式(7-5),经过整理,得到 式中,T是采样间隔,T的选取应按照滤波器最高截止频率的2倍以上选取,若T选取过大,则会使 附近频率混叠现象严重。这里选取T1s和T0.1s两种情况,以便进行比较。 设T1s时用 表示,T0.1s时用 表示,则,,一、用冲激响应不变法设计IIR数字低通滤波器,它们的幅频特性如图7-2所示。图7-2a表示模拟滤波器的幅频特性,图7-2b表示T1s,转换成数字滤波器的幅频特性,图7-2c
10、表示T0.1s,转换成数字滤波器的幅频特性。很明显,T0.1s时,它的幅频特性和模拟滤波器的幅频特性很近似,只是在折叠频率附近有很轻的混叠现象。而对于T1s情况,频率混叠现象很严重。,,一、用冲激响应不变法设计IIR数字低通滤波器,a),b),,一、用冲激响应不变法设计IIR数字低通滤波器,也可以利用MATLAB的函数impinvar实现冲激响应不变法模拟到数字的滤波器转换。,c) 图 7-2 例7-1的幅频特性,,二、用双线性变换法设计IIR数字低通滤波器,冲激响应不变法的主要缺点是会产生频率混叠现象,使数字滤波器的频响偏移模拟滤波器的频响。为了克服这一缺点,可以采用非线性频率压缩方法,将整
11、个S平面频率轴上的频率范围压缩到 之间,再用 转换到Z平面上。设 , ,经过非线性频率压缩后用 , 表示,这里用正切变换实现频率压缩:,(7-10),,二、用双线性变换法设计IIR数字低通滤波器,式中,T是采样间隔,当 从 经过0变换到 时, 则由 经过0变换到 ,实现了S平面上整个虚轴完全压缩到 平面上的虚轴的 之间的转换。这样便有 再通过 转换到Z平面上,得到,(7-11),(7-12),(7-13),,二、用双线性变换法设计IIR数字低通滤波器,式(7-12)或式(7-13)称为双线性变换。从S平面映射到 平面,再从 平面映射到Z平面,其映射过程如图7-3所示。由于从S平面到 平面具有非
12、线性频率压缩的功能,因此不可能产生频率混叠现象,这是双线性变换法比较冲激响应不变法最大的优点。另外,从 平面转换到Z平面仍然采用标准转换关系 , 平面的 之间水平带的左半平面映射到Z平面单位圆内部,虚轴映射成单位圆。这样,当 因果稳定,转换成的H(z)也是因果稳定的。,,二、用双线性变换法设计IIR数字低通滤波器,下面分析模拟频率 和数字频率 之间的关系。 令 , ,代入式(7-12)中,有,图 7-3 双线性变换法的映射关系,(7-14),,二、用双线性变换法设计IIR数字低通滤波器,上式说明,S平面上 与Z平面的 成非线性正切关系,如图7-4所示。在 0附近接近线性关系;当 增加时, 增加
13、得愈来愈快;当 趋近 时, 趋近于 。正是因为这种非线性关系,消除了频率混叠现象。,图7-4 双线性变换法的频率变换关系,,二、用双线性变换法设计IIR数字低通滤波器,与 之间的非线性关系是双线性变换法的缺点,直接影响数字滤波器频响逼真的模仿模拟滤波器的频响,幅频特性和相频特性失真的情况如图7-5所示。这种非线性影响的实质问题是:如果 的刻度是均匀的,则映射到Z平面的 的刻度是不均匀的,而是随 增加愈来愈密。但是,如果模拟滤波器的频响具有片段常数特性,则转换到Z平面数字滤波器仍具有片段常数特性,主要是特性转折点频率值与模拟滤波器特性转折点的频率值成非线性关系。当然,对于不是片段常数的相位特性仍
14、有非线性失真。因此,双线性变换法适合片断常数特性的滤波器的设计。实际应用中,一般设计滤波器通带和阻带均要求是片段常数,因此双线性变换法得到了广泛的应用。在设计时要注意边界频率如通带截止频率、阻带下截止频率等的转换关系要用式(7-14)计算。,,二、用双线性变换法设计IIR数字低通滤波器,图7-5 双线性变换法幅度和相位特性的非线性映射,,二、用双线性变换法设计IIR数字低通滤波器,双线性变换法可由简单的代数公式式(7-12)将 直接转换成H(z),这是该变换法的优点。但当阶数稍高时,将H(z)整理成需要的形式,也不是一件简单的工作。为简化设计,可以将模拟滤波器各系数和经双线性变换法得到的数字滤波器的各系数之间关系,列成表格供设计时使用。 双线性变换也可以利用MATLAB函数bi