测试信号分析与处理宋爱国第7章节数字滤波器的设计.ppt

第七章 数字滤波器的设计,,概述,7.1,,IIR数字滤波器设计,7.2,,FIR数字滤波器的设计,7.3,,智能仪器中常用的数字滤波算法简介,7.4,内容提要,本章主要介绍IIR数字滤波器和FIR数字滤波器的常用设计方法，简单介绍智能仪器仪表中常用的数字滤波算法第一节 概述,一个数字滤波器可以用一个N阶差分方程来描述，即,,（7-1）,或者以它的系统函数H(z)来描述,,（7-2）,第一节 概述,数字滤波器的设计大致包括以下4个步骤： 1）根据任务需要，确定数字滤波器应达到的性能指标，如通带截止频率 、阻带截止频率 、通带起伏 等此外还必须确定采样周期T或采样频率 2）确定数字滤波器的系统函数H(z)或h(n)，使其频率特性满足技术指标要求第一节 概述,3）用一个有限精度的运算去实现H(z)或h(n)，包括选择合理的网络结构、恰当的有效字长，以及有效数字的处理方法等 4）确定工程实现方法用实际数字系统（通用计算机软件或专用数字滤波器硬件）实现H(z)或h(n)第二节 IIR数字滤波器设计,设计IIR数字滤波器的方法可以归纳为两类： (一)模拟——数字转换法 (二)是直接设计法,,,,利用模拟滤波器成熟的理论和方法来设计IIR数字低通滤波器的设计过程是：按照技术要求设计一个模拟低通滤波器，得到模拟低通滤波器的传递函数 ，再按一定的转换关系将 转换成数字低通滤波器的系统函数 。

这样设计的关键问题就是寻找这种转换关系，将S平面上的 转换成Z平面上的 为了保证转换后的 稳定且满足技术要求，对转换关系提出两点要求：,,,,,,,,1）因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的模拟滤波器因果稳定要求其传递函数 的极点全部位于S平面的左半平面；数字滤波器因果稳定则要求 的极点全部在单位圆内因此，转换关系应是S平面的左半平面映射至Z平面的单位圆内部 2）数字滤波器的频率响应模仿模拟滤波器的频率响应，S平面的虚轴映射为Z平面的单位圆，相应的频率之间成线性关系 将传递函数 从S平面转换到Z平面的方法有多种，工程上常用的是冲激响应不变法和双线性变换法主要内容,一、用冲激响应不变法设计IIR数字低通滤波器,冲激响应不变法是使数字滤波器的单位冲激响应h(n)等于模拟滤波器的单位冲激响应的等间隔采样即,,（7-3）,式中，T为采样间隔 设模拟滤波器 只有单阶极点（若有多重极点，则求拉氏反变换会复杂一些），且分母多项式的阶次高于分子多项式的阶次，将 用部分分式表示,,,,（7-4）,一、用冲激响应不变法设计IIR数字低通滤波器,其中， 为 的单阶极点。

将 进行拉氏反变换得到,,,,,,式中， 是单位阶跃函数对 进行等间隔采样，采样间隔为T，得到,,,,对上式进行Z变换，得到数字滤波器的系统函数H（z）,,（7-5）,一、用冲激响应不变法设计IIR数字低通滤波器,对比式（7-4）与式（7-5）， 在S平面上的极点 ，根据 的关系映射为 在Z平面上的极点 ，系数 不变化 模拟信号 的傅里叶变换 和其采样信号 的傅里叶变换 之间的关系,,,,,,,,,,,,（7-6）,一、用冲激响应不变法设计IIR数字低通滤波器,将 代入上式，得,,,（7-7）,因为数字滤波器的系统函数是由采样信号 得来的，因此有,,,（7-8）,上式表明将模拟信号 的拉氏变换在S平面上沿虚轴按照周期 延拓后，再映射到Z平面上，就得H（z）一、用冲激响应不变法设计IIR数字低通滤波器,设,,,即,,（7-9）,由第三章的知识可知，S平面的虚轴（ ）映射成Z平面的单位圆（r＝1），S平面左半平面（ ）映射到Z平面的单位圆内（ ），S平面右半平面（ ）映射到Z平面的单位圆外（ ）这说明如果 因果稳定，转换后得到的H（z）仍是因果稳定的。

一、用冲激响应不变法设计IIR数字低通滤波器,另外，注意到 是一个周期函数，可写成,,,，M为任意整数,当模拟频率 从 变化到 时，数字频率 则从 变化 到，且按照式（7-9）， ，即 与 之间成线性关系 但是，从模拟信号 到采样信号 ，其拉氏变换要按照式（7-7），以 为周期，沿虚轴方向进行周期化如果原模拟信号 的频带不是限于 之间，则会在 的奇数倍附近产生频率混叠，从而映射到Z平面，在 附近产生频率混叠冲激响应不变法的频率混叠现象如图7-1所示一、用冲激响应不变法设计IIR数字低通滤波器,这种频率混叠现象会使设计出的数字滤波器在 附近的频率特性，程度不同地偏离模拟滤波器在附近的频率特性，严重时使数字滤波器不满足给定的技术指标为此，希望设计的滤波器是带限滤波器，如果不是带限的，例如高通滤波器、带阻滤波器，需要在高通带阻滤波器之前加保护滤波器，滤除高于折叠频率 以上的频带，以免产生频率混叠现象但这样会增加系统的成本和复杂性，因此，高通与带阻滤波器不适合用这种方法设计图7-1 冲激响应不变法的频率混叠现象,,,一、用冲激响应不变法设计IIR数字低通滤波器,假设 没有频率混叠现象，即满足,,,＝0，,,按照式（7-8），并将关系式 ， 代入，得,,,,说明用冲激响应不变法设计的数字滤波器可以很好地重现原模拟滤波器的频响。

上式中， 的幅频特性与采样间隔成反比，这样当T较小时， 就会有太高的增益为避免这一现象，令 ，那么,,,,一、用冲激响应不变法设计IIR数字低通滤波器,此时， 综上所述，冲激响应不变法的优点是频率坐标变换是线性的，即 ，如果不考虑频率混叠现象，用这种方法设计的数字滤波器会很好地重现原模拟滤波器的频率特性另一个优点是数字滤波器的单位冲激响应完全模仿模拟滤波器的单位冲激响应，时域特性逼近好缺点是会产生频率混叠现象，适合低通、带通滤波器的设计，不适合高通、带阻滤波器的设计一、用冲激响应不变法设计IIR数字低通滤波器,例7-1 已知模拟滤波器的传递函数 为 用冲激响应不变法将 转换成数字滤波器的系统函数H(z) 解：首先将写成部分分式 极点为 ， 那么H（z）的极点为 ，,,,,,,,,,一、用冲激响应不变法设计IIR数字低通滤波器,按照式（7-5），经过整理，得到 式中，T是采样间隔，T的选取应按照滤波器最高截止频率的2倍以上选取，若T选取过大，则会使 附近频率混叠现象严重这里选取T＝1s和T＝0.1s两种情况，以便进行比较 设T＝1s时用 表示，T＝0.1s时用 表示，则,,,,,,一、用冲激响应不变法设计IIR数字低通滤波器,它们的幅频特性如图7-2所示。

图7-2a表示模拟滤波器的幅频特性，图7-2b表示T＝1s，转换成数字滤波器的幅频特性，图7-2c表示T＝0.1s，转换成数字滤波器的幅频特性很明显，T＝0.1s时，它的幅频特性和模拟滤波器的幅频特性很近似，只是在折叠频率附近有很轻的混叠现象而对于T＝1s情况，频率混叠现象很严重一、用冲激响应不变法设计IIR数字低通滤波器,a）,b),一、用冲激响应不变法设计IIR数字低通滤波器,也可以利用MATLAB的函数impinvar实现冲激响应不变法模拟到数字的滤波器转换c) 图 7-2 例7-1的幅频特性,,,,,,二、用双线性变换法设计IIR数字低通滤波器,冲激响应不变法的主要缺点是会产生频率混叠现象，使数字滤波器的频响偏移模拟滤波器的频响为了克服这一缺点，可以采用非线性频率压缩方法，将整个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平面到 平面具有非线性频率压缩的功能，因此不可能产生频率混叠现象，这是双线性变换法比较冲激响应不变法最大的优点另外，从 平面转换到Z平面仍然采用标准转换关系 ， 平面的 之间水平带的左半平面映射到Z平面单位圆内部，虚轴映射成单位圆这样，当 因果稳定，转换成的H（z）也是因果稳定的二、用双线性变换法设计IIR数字低通滤波器,下面分析模拟频率 和数字频率 之间的关系 令 ， ，代入式（7-12）中，有,,图 7-3 双线性变换法的映射关系,,,,,,,（7-14）,二、用双线性变换法设计IIR数字低通滤波器,上式说明，S平面上 与Z平面的 成非线性正切关系，如图7-4所示在 ＝0附近接近线性关系；当 增加时， 增加得愈来愈快；当 趋近 时， 趋近于 正是因为这种非线性关系，消除了频率混叠现象图7-4 双线性变换法的频率变换关系,二、用双线性变换法设计IIR数字低通滤波器,与 之间的非线性关系是双线性变换法的缺点，直接影响数字滤波器频响逼真的模仿模拟滤波器的频响，幅频特性和相频特性失真的情况如图7-5所示。

这种非线性影响的实质问题是：如果 的刻度是均匀的，则映射到Z平面的 的刻度是不均匀的，而是随 增加愈来愈密但是，如果模拟滤波器的频响具有片段常数特性，则转换到Z平面数字滤波器仍具有片段常数特性，主要是特性转折点频率值与模拟滤波器特性转折点的频率值成非线性关系当然，对于不是片段常数的相位特性仍有非线性失真因此，双线性变换法适合片断常数特性的滤波器的设计实际应用中，一般设计滤波器通带和阻带均要求是片段常数，因此双线性变换法得到了广泛的应用在设计时要注意边界频率如通带截止频率、阻带下截止频率等的转换关系要用式（7-14）计算二、用双线性变换法设计IIR数字低通滤波器,,图7-5 双线性变换法幅度和相位特性的非线性映射,二、用双线性变换法设计IIR数字低通滤波器,双线性变换法可由简单的代数公式式（7-12）将 直接转换成H（z），这是该变换法的优点但当阶数稍高时，将H（z）整理成需要的形式，也不是一件简单的工作为简化设计，可以将模拟滤波器各系数和经双线性变换法得到的数字滤波器的各系数之间关系，列成表格供设计时使用 双线性变换也可以利用MATLAB函数bilinear完成bilinear的语法为 [zd,pd,kd]=bilinear(z,p,k,fs),,二、用双线性变换法设计IIR数字低通滤波器,例7-2 试分别用冲激响应不变法和双线性变换法将图7-6所示的RC低通滤波器转换成数字滤波器。

图7-6 RC低通滤波器,解:首先按照图7-6写出该滤波器的传递函数 为 ， 利用冲激响应不变法转换，数字滤波器的系统函数 为,,,,,,二、用双线性变换法设计IIR数字低通滤波器,利用双线性变换法转换，数字滤波器的系统函数 为,,,,，,,和 的网络结构分别如图7-7a、b所示a）,b）,图7-7 例7-2中,和,的网络结构,,,,,二、用双线性变换法设计IIR数字低通滤波器,设＝1000，T＝0.001s和0.002s， 和 的归一化幅频特性如图7-8所示图7-8a是模拟滤波器幅频特性，是一个低通滤波器，但拖了很长的尾巴图7-8b是采用冲激响应不变法转换成的数字滤波器幅频特性，图中 处对应的模拟频率与采样间隔T有关，当T＝0.。