(贾)第四章IIR滤波器的设计

1、 第四章第四章 无限长单位脉冲响应无限长单位脉冲响应（IIRIIR）滤波器设计）滤波器设计 许多信息处理过程，如信号的过滤，检测、预测等都要用到滤波器，数字滤波器是数字信号处理中使用得最广泛的一种具有频率选择性的离散线性系统， 数字滤波器的功能（本质）是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列。 数字滤波器是具有频率选择性的离散线性系统，数字滤波器的设计步骤大体可分为以下三步：1.1.根据使用要求确定系统需要的各种性能指标根据使用要求确定系统需要的各种性能指标( (指容限、容差指容限、容差) ) 通常在频域给出滤波器的容限、容差，如图所示。 容限容限-指滤波器通带、阻带的

2、最大衰耗和波动范围。 容差容差-指通带、阻带容许的最大 误差范围。 第四章 IIR滤波器的设计方法在通带内，幅度响应以误差 逼近于1，即1111jHec在阻带中，幅度响应以误差 逼近于零， 即22jHep其中: 、 分别为通带和阻带的截止频率 。cp2. 2. 用稳定的因果系统去逼近设计的系统函数用稳定的因果系统去逼近设计的系统函数 系统函数包括：无限长单位脉冲响应的系统函数和有 限长单位脉冲响应的系统函数，可根据容限选择IIR或FIR 滤波器。 在滤波器设计中逼近逼近是一个重要的环节。 Hz逼近逼近就是给定所要求的滤波器性能后，去寻找一个物理可实现的系统函数，使它的频率特性尽可能近似所要求的

3、滤波器特性，也就是指对理想特性进行逼近，最后得到一个因果、稳定并且可实现的传递函数 。 滤波器的设计实际上就是一个数学逼近的问题。滤波器的设计实际上就是一个数学逼近的问题。 zH3.3.决定决定 的实现方法的实现方法 ( (包括选择运算结构、运算制式及字长)。因为数字滤波器是一个物理可实现的LTI系统，所 以它是用一个有限精度的运算去实现系统函数 。选择的结构形式有：级联、并联、直接、横截、频率采样型。 本章和第五章主要解决第二个步骤的内容，也就是寻求滤波器传递函数设计的问题。 Hz Hz DF -DF - 数字滤波器 AF -AF - 模拟滤波器 单位脉冲响应 单位冲激响应 采样序列 DF的

4、传递函数 AF的传递函数 DF的频率响应 AF的频率响应 nh athaTnh zH ashjHeajH说明一下符号表示的含义：数字滤波器四种类型的理想频率特性：数字滤波器四种类型的理想频率特性：4.1 IIR滤波器设计的特点 IIR滤波器的传递函数可以写成N阶的有理函数： 0111111NNNNiiiiiiiiiiiiAHzzzdcazbz 滤波器设计的核心是求传递函数 ，而 的设计就是确定系数 、 ，或者确定零点 、极点，使得滤波器的 满足给定的性能要求。 zH zHia ib icid zH IIRIIR 指单位脉冲响应 为无限长的滤波器，滤波 器的 有无限个离散值。 nh nh 如果在

5、单位圆内设置一对共轭极点 ，频响在 将有一峰值。0jZr e0 极点 离单位圆愈远，频响在 处的峰值比较平缓。 极点 越接近单位圆，频响在 处的峰值就越尖锐。 0001.累试法 滤波器的幅度特性和相位特性完全由其零点和极点的位置所决定。0一.IIR 滤波器的一般设计方法 零、极点累试法是在频域直接设计，可以完成一些零、极点累试法是在频域直接设计，可以完成一些简单的、阶数不高的滤波器的设计。简单的、阶数不高的滤波器的设计。 例如，实际中常常要求滤除叠加在信号上的50HZ工频干 扰，我们可以在 处安排一个零点， 就可以滤除掉50HZ的工频干扰。1jez11250,TTHz 如果通带太窄，可以把极点

6、向原点平移，如果通带太宽，可以将极点向单位圆移动。我们可以通过几次调整单位圆内极点的位置，去逼近符合设计要求的频响。同样，如果在单位圆上 处设置一个零点，那么频响在 处会出现传输零点，传输零点可以实现陷波作用传输零点可以实现陷波作用。1jze1 就是计算机优化设计，计算机优化设计的思想是借助计算机，使得所设计的滤波器的频响尽可能逼近理想的频响，即: jeH逼近2. 最优化设计法最优化设计一般分为两步进行：最优化设计一般分为两步进行：第一步：选择最优化标准或准则第一步：选择最优化标准或准则 选择一种最佳准则，使得 与 的均方误均方误差最小差最小或者或者最大误差最小最大误差最小。jeHjdeHdj

7、He011NNiiiiiiHabzzzDF的传递函数 通过改变 的系数 、 ,分别计算均方误差 ， 经过多次迭代运算，寻找一组系数 、 ，使得均方误差 为最小的这组系数为最优系数，从而完成最优化设计。 zHia ib E E ib第二步：进行迭代运算，确定最优系数第二步：进行迭代运算，确定最优系数 21MdijjEHHee最小根据最小二乘法准则，要求:ia 3.模仿模拟滤波器的设计 因为模拟滤波器的设计目前已经很完善，AF不仅有简单和严格的设计公式，而且它的设计参数也已经表格化，因此，我们可以借助于模拟滤波器设计的成果来设计数字滤波器。 在模拟系统中，利用工作参数综合法设计滤波器时，无论低通、

8、高通、带通和带阻滤波器，均是先设计一个低通原型，然后经过某种频率变换完成所要求设计的滤波器。即:频率变换低通 （LP)高通（HP)低通原型带通（BP)带阻（BS) 利用模拟滤波器设计数字滤波器，首先利用模拟滤波器的现成结果，在S平面设计出符合要求的模拟滤波器的传递函数 ，再通过一定的映射关系，得到数字滤波器的传递函数 。即： sHa zH( )( )aSHH z线性映射S平面 Z平面 ()()jdH eHj 逼近 使得所设计的数字滤波器的频响去逼近理想模拟滤波器的频响。 在进行IIR 数字滤波器的设计时，要逼近模拟原型低通滤波器，模拟低通滤波器通常以幅度平方函数表示模特性以幅度平方函数表示模特

9、性。二. 最常用的几种模拟原型低通滤波器的逼近方法 jSaaaaasHsHjHjHjH21巴特沃思滤波器（butterworth）最平响应滤波器 巴特沃思低通滤波器的幅度平方函数定义为： NCNaaajjjHjHjH2*2111N 为整数，表示滤波器的阶次； 为截止频率。 c当=0 时，当时，有 又称为滤波器的3bB带宽（或半功率点）。 102jHac212CajH当 时，当 时， 又称为滤波器的3bB带宽 0巴特沃思低通滤波器的巴特沃思低通滤波器的特点：特点： 在 处，即靠近零频处，衰减为0 ，所以巴特沃思滤波器通 带内具有最大平坦的振幅特性，故得名为最平坦响应滤波器。 巴特沃思低通滤波器没

10、有有限零点,零点出现在 处，它属于 “全极点型滤波器全极点型滤波器”。0 2切比雪夫滤波器 (chebychev) 切比雪夫滤波器也是一种全极点型滤波器，它的幅度平方函数 22211CaNjHC其中： 表示通带波纹大小，是小于 1的正数， 越大，波纹越大。 为滤波器的截止频率，但 并不是3db带宽）。cc22maxmax10102minmin20 log10 log10 log13HHdBHHBd带宽 为N阶切比雪夫多项式，定义为：NC1111coscosCCCCCNNNCchch 切比雪夫低通滤波器的特点：1) 通带内等起伏，通带外衰减快；2) 由于过渡带较窄，因此相位特性较差。 3考尔滤波

11、器(cauer)2211aHjKj 112222nSSSSSBsK 考尔滤波器的特点：1) 通带内、外都是等起伏。2) 由于过渡带较窄，因此相位特性较差。三三S S平面到平面到Z Z平面的映射变换平面的映射变换 利用模拟滤波器来设计数字滤波器，就是从已知的AF 的 去设计DF的 ，即: zH aHs( )( )aH SH Z通过变换 这种变换归根结底是一个由S平面到Z平面的变换，并且通常是复变函数的映射变换，这种映射变换映射变换应该满足两个基本的要求两个基本的要求： 的频响应该模仿 的频响 即要求 zH sHa 是因果稳定的映射 指 的因果稳定性通过映射后， 仍应保持因果 稳定。 sHa zH

12、( )( )aH SH Z通过变换4.2 4.2 脉冲响应不变法脉冲响应不变法 根据容限设计好一个模拟滤波器后，就可对此模拟系统进行模仿。数字滤波器从什么角度去模仿模拟滤波器呢？第一种方法是脉冲响应不变法。 1. 脉冲响应不变法 模拟 系统 txta thsHaa tythaa dthxthtxaaaaLTI系统特性可完全由它的冲激响应决定 数字系统 nxn nhzH nynh Kknhkxnhnx 脉冲响应不变法让数字滤波器的脉冲响应和模拟滤波器的脉冲响应在采样点上完全一样。即：( )aTnh nh单位脉冲响应不变法的设计思想：使数字滤波器从时域去模仿模拟滤波器。 sHa 1L tha采样

13、变换Z zH2脉冲响应不变法设计的系统的频率响应模拟系统频响数字系统频响 jHa thFTajeH( )jZeZh n对模拟信号 txa理想采样 txa 1L tha()( )aTnhh n naanTttxtx aastXstdtxe拉氏变换 nstattnTdtex nTasnTnex如果( )aTnnxx ( )nnnnaTXznnzxxz 当 时， STZe sXeXaST 是采样信号的拉氏变换与采样序列Z变换之间的映射关系STZe理想采样的频谱 12maaaSjjjjmTTXXsX 采样信号的频谱是采样信号在虚轴上的拉氏变换 是模拟信号频谱 以采 样频率 为周期的周期延拓ajXajX

14、2ST或 jTj TjaZ eeXjX zX e 理想采样信号的频谱是单位圆上的Z变换，它也等于数字信号的频谱。 因此，我们可以得到采样序列 的频谱 与原信号 的频 谱 之间的映射关系为：)(nx txajXajeX12mjajjmTTXXe 上式表明上式表明：当采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它实际上是完成了完成了由由S S平面平面到到Z Z平面平面的的映射映射，我们利用以上关系同样可以得到数字滤波器的频响 与模拟滤波器的频响 之间的关系：jeHjHa1212()mmaajTmHjjmHjTTTTH e 时表达式表明表达式表明：数字滤波器的频响并不是简单的重现模拟滤波器的频

15、响， 为模拟滤波器频响 的周期重复,是以 为周期的周期延拓。jeHjHas需要强调指出，在周期重复模拟滤波器频响的过程中所存在的问题：首先分析一下 和 以及 的关系： jeHajHajH 注意注意：我们原本是让 来模拟 ，但由于 ，所 以 模拟的是采样信号的频谱 ， 的映射关系反映 的是 和 之间的对应关系，而并不是 和 之 间的关系。jeH( )aTnh nhjeHajHSTZejeHjHajHajeHjHa 如果模拟信号不充分带限，这种周期重复就不可避免的存在着频率混叠问题，如下图所示。 只有当信号是带限信号，并且满足或 时，这时数字滤波器的频响 在折叠频率以内才能不失真的重现模拟滤波器的

16、频响 。即ajH2Sm 、TmmjeH11()jaajjTTTHHHe 但任何一个实际的模拟滤波器都不可能完全是带限的，故不可避免的要出现混叠现象。 如果模拟滤波器频响在折叠频率以上衰减的愈快，则混叠失真就愈小。因此，脉冲响应不变法只适用于带限模拟滤波器。对于带阻和高通滤波器，由于它们高端频率不衰减，因此高频分量将完全混叠在低频响应中，从而使整个频响产生失真。如果设计的高通和带阻滤波器要采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，先滤除掉高于折叠频率以上的频带，然后再使用脉冲响应不变法将其转换为数字滤波器。3模拟滤波器的数字化方法脉冲响应不变法要经历以下四个过程： sHa t

17、ha1L采样( )aTnnhh 变换Z zH 011NiNiaiiiiHa ssb s 用部分分式展开1iNiiASS1L 1NiiiSatAtteuh对 进行采样，得到数字滤波器的单位脉冲响应序列 aht11( )( )( )iiNNiiiiSTS TannnTAAnheehu nu n（1 1）（2 2）（3 3） 111NSiTiiAHzze比较第（1）步和第（4）步 S平面的每一个极点 Z平面是极点 ； 和 的系数相同。因此，只要知道极点 就可以直接得到而不必求解第(2)、(3)步。 变换到 H z aH s H zTiSieZ iS（4 4）从模拟域从模拟域 直接映射到数字域直接映射

18、到数字域)(nhZ 变换 1NaiiiAHsSS 111iNiiS TAH zzeiS4对脉冲响应不变法的修正数字域频响与模拟信号采样后的频谱关系是：12amjmHHjTTe 上式中频响频响 与与T T成反比成反比，通常为了减小脉冲响应不变法的混叠失真，一般选择较高的采样频率，采样周期T就会比较小，如果T很小，数字滤波器的增益就会很高，通常希望 的增益能与采样频率无关(与T无关)。因此实际应用中，脉冲响应不变法要稍做一下修改。jeHjeH令 ， 则 ( )aTnTnhh NiTiSizeTAzH1112()mjamjTHHe 这样，数字滤波器的增益就不会随采样频率的变化而变化。下面举例说明脉冲

19、响应不变法的应用：下面举例说明脉冲响应不变法的应用：例：已知AF的系数函数 用脉冲响应不变法求出相应的数字滤波器的系数函数 aHs22sasab zH解： 1122()()asajbsajbsHAF在 有一个零点；在 有一对共轭极点。OSaPSajb 11112211jbTjbTaTaTTTH zeezeez11221cos12cosaTaTaTzbTTbTzzeee TbjTaTjbTaTaeezeezbTezzcos数字域中 有两个零点 和 ，有一对共轭极点 。可见，脉冲响应不变法脉冲响应不变法对零点没有一一对应的关系对零点没有一一对应的关系。0ZcosaTZTbeaTjbTZee H z

20、从幅频特性频谱中可以看出：由于频谱混叠带来明显的失真。很显然， 在 处比 在 处下降的慢，这主要是因为“混叠”造成的。混叠的主要原因是由于模拟信号频域不充分带限。 jeHaHjT5 5脉冲响应不变法的特点脉冲响应不变法的特点 脉冲响应不变法是一种时域设计法，即： 这种时域模仿使得S平面的极点 ( )aTnh nh但是S平面 Z平面的零点并没有这种单值映射的关系。Z平面的极点为 iS变换到iSiTZe sHa ath1L采样( )aTh nnh 变换Z zH脉冲响应不变法设计思想是 1NiiiaASSHs NiTSizeAzHi111频率坐标是线性变换 。因为模拟和数字频率之间的变换关系是线性的

21、，所以可以设计出满足线性相位的滤波器。只要不发生混叠，变换后数字滤波器的频响可以不失真的恢复出模拟滤波器的频响。即：T ajHHjeT或注意注意：数字域的频响特性都是以 为周期的周期函数，其特性对称于折叠频率，因此， 由 是 的镜像，通常数字滤波器的频率特性我们只考虑 范围的特性即可。2 20 脉冲响应不变法最大的缺点是有频响的混叠，所以只能模仿带限的频响。因此，对衰减特性很好的低通或带通，频响混叠小。对于高通和带阻滤波器，由于高端频率不衰减，所以要加保护滤波器滤掉高于折叠频率以上的部分，再用脉冲响应不变法转化成数字滤波器。0 脉冲响应不变法的最大缺点是产生混叠失真。产生的主要原因是由于S平面

22、到Z平面的变换关系 不是单值的，而是多值的映射关系，由于S平面的多值性造成频率的交叠，产生混叠现象。为了克服这一缺点，我们有必要寻求S平面与Z平面的单值的映射关系。 STeZ下面我们就介绍另外一种变换方法双线性变换法双线性变换法。4.3 双线性变换法 为了寻求S平面与Z平面的一一对应关系，我们先把S平面压缩到某一个中介的 平面的横带内（宽度从 到 ），然后把此横带再变换到整个Z平面，这时S平面与Z平面就建立了一一对应的单值的映射关系，可以消除混叠现象。如下图所示：1STT1. 双线性变换法： 为了完成从 的变换过程，我们寻找它们的 映射关系:1SSZ平面平面平面 虚轴 轴平面jS 压缩到横带内

23、轴的平面TTjS,11这个映射可以通过正切变换来实现： 12Ttg 验证： 当 时， 当 时， 当 时， T12tg 001T12tg 的映射关系：平面平面1SS 21Ttg 221222211111TjtheeeeTtgjjTjTjTjTj由 得到 把此关系解析延拓（扩大）到整个S平面令 ， ，则Sj11Sj11111112222211S TS TS TS TS TS TSTeeeStheee1111S TS TeSe得到 的映射关系平面平面1SS 映射到 S平面Z平面将标准变换关系 代入 ，从而得到 到 的双线性变换的单值映射关系1S TZe1111S TS TeSe平面S平面Z11111

24、1ZSZSZS现在来验证上式中S和Z的变换关系是否满足“映射变换”的两个总要求：(1)虚轴关S平面j与Z平面的系将 、 代入双线性映射关系 中， jZeT11111ZZS12121jjTj tgj tgjeSe 通过验证证实：S平面虚轴 Z平面的单位圆上，表明整个虚轴是单值的对应于单位圆的一周。因此双线性变换确实消除了脉冲响应不变法所存在的混叠误差，所以它的逼近是良好的。()j 映射到jZe 验证S平面与Z平面是否是因果和稳定的jS1111sjjSZSj 2222)(1)(1jjZ当 时， ，S左半平面 Z平面的单位圆内 ;当 时， ，S平面虚轴 Z平面的单位圆上 ;当 时， ，S右半平面 Z

25、平面的单位圆外 。01Z 映射到01Z 映射到01Z 映射到 如果模拟滤波器是稳定的，通过双线性变换后，所得 到的数字滤波器也一定是稳定的。 如果给定模拟滤波器的传递函数 ,变量S与Z之间 有简单的代数关系，只要用代数置换就可以得到数字 滤波器的传递函数。 sHa11111111( )( )ZZaaSzzH zHsH通过验证得到以下结论：可看出：可看出：模拟域与数字域的频率变换关系不是线性的模拟域与数字域的频率变换关系不是线性的 而是非线性关系而是非线性关系。数字滤波器的频响：2()2jaatgjtgjHHHe 上式说明，经过 频率变换，DFDF的频响的频响 等于等于AFAF的频响的频响 。2

26、tgjeHjHa模拟域与数字域的频率关系为： 122Ttgtg2双线性变换的频率响应 S平面到 平面虚轴的映射关系为：1S12Ttg 没有频谱混叠没有频谱混叠 双线性变换的最大优点是避免了频响的混叠效应，因为S平面的虚轴单值对应着Z平面单位圆的一周。从它们的频率变换关系中清楚看出，当 时，终止在折叠频率处，故不会有高于折叠频率的分量出现，避免了由频响产生的混叠失真，但双线性变换虽然避免了频响混叠，但却是以牺牲相位特性为代价换取的。 和和 的关系是非线性变换的关系是非线性变换2tg 如果模拟系统是按时延设计的线性相位的模拟滤波器，经过双线性变换后，数字滤波器将得不到这种线性相位特性，所得到的数字

27、系统的相位不再保持线性相位特性。3双线性变换的优缺点 数字频率和模拟频率之间的关系是例：模拟系统中三阶切比雪夫低通滤波器（Chebyshev)从图中可以看出在靠近截止频率处有频率畸变。 模拟滤波器截止频率为 ，对应的数字域频率为 ，而经过双线性变换后数字域的频率是 ， 。可见，经过非线性变换以后，在截止频率处的频率特性发生了畸变。因此，如果要求相位和频率是线性关系，就不能采用双线性变换法。 cccT 12cctgcc 双线性变换适合分段常数幅频响应的滤波器，而实际中大多数滤波器均为分段常数响应滤波器，如典型的低通、高通、带通、带阻滤波器都是具有分段常数响应的滤波器，所谓分段常数响应的滤波器指分

28、段常数响应的滤波器指在通带内要求逼近一个衰减为0的常数，阻带内要求逼近一个衰减为 的常数，所以我们完全可以采用双线性变换。双线性变换所引起的频率畸变只是出现在数字滤波器靠近截只止频率处。这种畸变完全可以通过预畸变加以校正。 首先对模拟滤波器的临界频率事先加以预畸变，让 ，也就是使 通过双线性变换后正好映射到所需的位置上去。2cctgC2cctgC 设计简单设计简单 双线性变换是IIR滤波器设计中使用最普遍、最有成效的一 种设计方法。由于S S与Z Z之间是简单的代数关系，因此传递函数可 以通过代数置换得到。即 11111111aZZaSZZHH zHs12Cctg2cctg 正好映射到所需的位

29、置上。c4.4 原型变换 模拟滤波器已经形成了许多成熟的设计方法，每种滤波器都有一套准确的计算公式和大量的归一化的设计表格和曲线。所以在模拟滤波器的设计中只要掌握原型变换，就可以通过归一化低通原型的参数，去设计各种实际的低通、高通、带通、带阻滤波器。模拟滤波器的这一套方法可以通过上述讨论的两种变换法而应用于数字滤波器的设计。下面我们讨论各种滤波器设计的问题： 一低通滤波器的设计步骤（由模拟低通到数字低通的设计） 确定数字滤波器的性能要求(容限)，确定各临界频 率 的值。c 将数字域 模拟域2cctg 通过变换关系c把DF的性能要求转换为相应的样本的AF的性能要求，得到模拟滤波器的临界频率 。

30、按照模拟滤波器临界频率 的要求设计模拟滤波 器的传递函数 。 通过脉冲响应不变法或者双线性变换，将 sHa sHa转换为 zHccc例例4.14.1要求用脉冲响应不变法及双线性变换设计一个三阶Butterworth 低通滤波器。设采样周期 (即采样频率 )，3dB截 止频率sT2504SZfkHZCkHf1解：第一步解：第一步: : 的设计的设计 Butterworth滤波器： aHs 22111CNNaaajjHsHs Hs 三阶Butterworth滤波器N=3，有2N=6个极点,分布在半径为 的圆上，将S平面分成6等份。为了保证设计的系统是稳定的， 应该取S左半平面的极点，因此，S左半平

31、面的传递函数是:c 2sHa 10NCNKKaSSHs其中其中： 21122KNKCjSe 1,1,0NK当N=3时，可直接确定三阶Butterworth滤波器的传递函数 为： sHa 32233CCCjjSSCSsaHee 322211CCCSSS第二步：变换法设计滤波器第二步：变换法设计滤波器 脉冲响应不变法脉冲响应不变法 将 用部分分式展开： sHa 6633131322jjccjjSScccsaScHee 将将 直接变换为数字滤波器的传递函数：直接变换为数字滤波器的传递函数： aHs6633113131112211TcTcTcTjcjTZcjjZTcZeezeeeH 663311313

32、1112211cjjjjZZZccccceezeeeH 将 代入 中：5 . 02TfCC zH 21112079.01905.015541.0571.12079.01571.1ZZZZzH 以上得到的是数字三阶Butterworth滤波器的传递函数 ，采用一个一阶基本节和一个二阶基本节并联的结构实现比较方便。 zH除此之外，我们还可采用另外一种方法计算 从上面实例中知道，当求模拟滤波器的传递函数时， 只是一个符号而已，它并不进行计算， 乘以T T后 变成 ，最终传递函数只取决于数字域的参数 。因此只需求出一个相对於采样频率归一化的采样频率归一化的“样本样本”传递传递函函数数即可。这种计算方法

33、可以直接以 作为模拟滤波器的临界频率 来计算(让 ），所得到的 是采样频率归一化后的模拟样本传递函数。 zH sHacCccCcCc sHa 21122CCCSSSaHs2311 22CCCccSSS aHs对采样频率归一化的模拟样本传递函数 展开成部分分式形式： sHa 111100dSAdSAdSAsHa 111110110111ZeAZeAZeAzHddd T=1已归一化双线性变换法设计步骤双线性变换法设计步骤 i) 首先确定数字域临界频率 ;C5 . 02TfCCii) 确定预畸变的模拟滤波器的临界频率 由于双线性变换法频率变换的非线性会引起频率畸变，所以要预畸变。 为了使数字滤波器的

34、临界频率为 ， 模拟滤波器的临界频率应为：C5 . 0C142tgtgCC5 . 021cCtg刚好映射到我们所要求的位置上。得到频率归一化的模拟滤波器的样本传递函数 122123SSSsHaiv）最后将双线性变换的关系式代入 ，得到数字滤波器的传递函数 sHa zH 112312111111 133123ZZaaSZZZZZZHH zHs注意：这里所采用的模拟滤波器 并不是数字滤波器所要模仿的截止频率为 的那个实际滤波器， 只是模拟滤波器的样本传递函数，是归一化的低通原型，它是由模拟滤波器的低通原型到数字滤波器变换的一个中间变换函数。 上面介绍了如何设计一个数字低通滤波器。 sHaZCkHf

35、1 sHaiii）将 代入 2311 22CCCSSSHsa1C sHa 当我们需要设计当我们需要设计高通、带通、带阻数字滤波器时，有两种方法：高通、带通、带阻数字滤波器时，有两种方法：第一种第一种：首先设计一个相应的高通、带通、带阻模拟滤波器，然后再通过脉冲响应不变法或双线性变换法转换为数字滤波器，即:的频带变换模拟到模拟变换法将其数字化变法或者双线性通过脉冲响应不模拟低通原型 模拟的高通、带通、带阻 数字的高通、带通、带阻模拟低通原型 数字的高通、带通、带阻频带变换第二种第二种：直接利用模拟滤波的低通原型，通过一定的频率变换关系，一步完成各种数字滤波器的设计，即上述两种方法的根本区别在于第

36、一种方法是在模拟域进行频带变换，第二种方法是在数字域进行频带变换。第一种方法的实现步骤是：首先确定临界频率 模拟滤波器的临界频率 ；完成模拟高通、带通、带阻滤波器的设计；最后通过脉冲响应不变法或双线性变换法将设计好的 转换成 即可。 第一种方法和讨论的低通滤波器的设计完全一样，因此不再讨论。这里着重介绍第二种方法，因第二种方法更简捷方便，故得到普遍采用。由于脉冲响应不变法对于高通、带阻滤波器不能直接采用，因此第二种方法仅讨论采用双线性变换法直接进行频率变换。c sHa zHc转换成相应的二高通变换 在模拟滤波器的高通设计中，低通至高通的变换就是S变量的倒置变换，数字高通滤波器只要把双线性变换中

37、 即可。1SS低通为 高通为111111ZSZSZS111111ZSZSZS 11111111ZZaahSzzHzHsH将 代入上式复变量S中jZe倒置代换211jjjeSjctge 数字域与模拟域的频率映射关系为数字域与模拟域的频率映射关系为 2ctg 由右图中看出： 当 ，Z平面是映射在 上 当 ，Z平面是映射在 上0 映射为1Z 映射为01Z 说明双线性变换中S变量的倒置变换并不会影响变换后的稳定条件，而且 轴仍然映射到Z平面的单位圆上，只是方向颠倒了。通过这样的变换，就可以直接把模拟低通变换为数字高通，如右图所示。j 高通变换的计算步骤和低通变换是一样的，只是在确定模拟原型预畸变的临界

38、频率时，应该采用2KKctg 例例4.34.3已知采样频率 , ZkHf10sT100ZSkHf5 . 2ZkH5设计一个三阶切比雪夫高通数字滤波器，其通过频率 （但不必考虑 以上的频率分量),通带内损耗不大于1dB。解：切比雪夫低通原型的模平方函数为CNaCjH22211通带损耗 时， ，当N=3时，通过查表可得传递函数：dBAP15089. 0 332230.49130.49131.2380.9883CaCCCHsSSS现在我们来求预畸的截止频率 ，首先确定数字域的截止频率 ，cc ，于是有 ，代入上式得：5 . 02TfCC124CCctgctg 329883. 0238. 14913.

39、 04913. 0SSSsHa 11111111ZZhaaLSZHzHsHZ 3213212041. 06043. 03432. 013311321. 0ZZZZZZzH三三. .带通变换带通变换根据带通滤波器变换的要求 S S平面平面 Z Z平面平面 0S 映射到0jZe0（ 是带通中心频率）S 映射到1Z模拟低通的模拟低通的 数字域的数字域的0 映射到0 映射到0, 满足以上带通变换要求的双线性变换为：满足以上带通变换要求的双线性变换为： 002022cos1( )111jjZeZeZZSF zZZZ验证F(z)是否满足带通变换:将零点 代入F(z) ，S0将极点 代入F(z) ，0jZe

40、1Z jS说明双线性变换满足带通变换的要求。 11cos2202ZZZS 2022022cos112cos11ZZZaaBPSZZZHzHsH直接将 代入 中，求得： sHa将变换关系将变换关系 、 代入代入F(z)=SF(z)=S中有：中有： jZeSjjjeeeeeeejjjjjjjjsincoscoscos211cos200202 数字域数字域 与模拟域与模拟域 之间的映射关系为：之间的映射关系为：sincoscos0由右图看出：在 点正好映射在 上 当 时，正好映射在 上， 当 时，正好映射在 上，因而频率变换关系完全满足带通变换的要求。 000低通幅度响应到带通幅度响应的变换如右图所

41、示带通变换的设计步骤带通变换的设计步骤 在设计带通时，通常给出上下边带的截止频率在设计带通时，通常给出上下边带的截止频率 、 和阶数和阶数N N 。首先根据。首先根据已知的已知的 和和 计算出数字带通的中心频率计算出数字带通的中心频率 和模拟低通的截止频率和模拟低通的截止频率 。12120c1011coscossin 2022coscossin011coscossinC12121212022cossinsinsincoscos模拟滤波器中 和 是一对镜像频率，即 ，且 也是模拟滤波器的截止频率 ，因此由上面两个关系式可解出 和1212 1c01ciiii）求归一化的模拟低通原型）求归一化的模拟

42、低通原型iiiiii）把满足带通变换的双线性变换代入模拟低通）把满足带通变换的双线性变换代入模拟低通 中。中。 sHa sHai i）求模拟低通的截止频率）求模拟低通的截止频率 c 2022022cos112cos11aZZZSZZZHHHzs 例例4.44.4：已知采样频率为已知采样频率为 ， ，要求设计一个三阶巴特沃，要求设计一个三阶巴特沃兹带通滤波器，其上下边带的兹带通滤波器，其上下边带的3dB3dB截止频率分别为截止频率分别为 和和 。ZSkHf100sT100ZkHf5 .371ZkHf5 .12275. 0211Tf25. 0222Tf解：首先确定数字频域的上下边带截止频率为：解：

43、首先确定数字频域的上下边带截止频率为：代入 中，21210sinsinsincos075. 0sin25. 0sin75. 025. 0sincos05 . 00011coscoscos0.5cos0.751sinsin0.75C 再带入 中c将求出的模拟低通的截止频率将求出的模拟低通的截止频率 代入三阶巴特沃兹的传递函数中：代入三阶巴特沃兹的传递函数中：1C 122123SSSsHa 2022cos11ZZSZH zH s4642333121ZZZZ三. 带阻变换 将带通滤波器的频率变换关系倒置一下就可得到带阻滤波器，整个设计方法、步骤和带通滤波器是类似的。带阻变换的关系为：1cos2102

44、2ZZZS0coscossin（ 有 ，分母有变化）j4.5 Z平面变换法 前面讨论的原型变换均是由模拟低通原型来设计各种数字滤波器，这种原型变换的设计方法同样也可以直接在数字域中进行。 如果我们已经有一个数字滤波器的低通原型函数如果我们已经有一个数字滤波器的低通原型函数 ，就可以直，就可以直接在数字域中进行变换，设计出所需要的数字滤波器的接在数字域中进行变换，设计出所需要的数字滤波器的 。即：。即： PHz zH数字域低通原型数字域低通原型 通过一定变换关系LPHPBPBS数字滤波器低通高通带通带阻即：即： zHP 通过映射关系 zH 为了区分变换前后两个不同的Z平面，我们规定规定： 把变换

45、前的把变换前的Z Z平面定义为平面定义为 U U平面平面，变换后的变换后的Z Z平面定义为平面定义为 Z Z平面平面。 U平面的单位圆用 表示，Z平面的单位圆用 表示。jeje从 U ZU Z平面的映射关系用变换函数 表达为：11ug zgU U平面平面 到到 Z Z平面的映射关系平面的映射关系实际上是复平面的保角映射复平面的保角映射。现在我们来分析一下函数 的一般特性：1g zI）变换以后的传递函数应该保持稳定性不变 要求U U平面平面的单位圆 必须映射到必须映射到 Z Z平面平面的单位圆 。jejejjjjggeeee 其中： 是 的相位函数，由上式可得到：jg e1jg e说明 在单位圆

46、上的幅度必须恒等于幅度必须恒等于1 1，这种函数称为全通函数全通函数。1g z全通函数全通函数任何频率的一个信号通过一个网络，输出幅度不变，只是 相位变化，即只对相位进行校正，这种网络称为全通网络， 表示全通网络的网络函数称为全通函数。)(1111Niiiizg zz 全通函数的表达形式全通函数的表达形式： 全通函数 的基本特性： N N 全通函数的阶数全通函数的阶数1zg 是极点，是极点， 可以是实数，也可以是共轭复数。且可以是实数，也可以是共轭复数。且 ， 是零点，所有的零点都是极点的共轭的倒数。是零点，所有的零点都是极点的共轭的倒数。ii1i1i可以证明，Z平面当 由0 0变化到 时，其

47、相位函数 由0 0变到 ， 因此可以因此可以由相位的变化来确定阶数由相位的变化来确定阶数N N。 N 只要选择合适的N和 就可以得到各类需要的变换。下面我们具体讨论各种原型变换：i 一一. . 低通低通 低通低通 jPeHjeH从U平面 Z平面到 及 均是低通滤波器的频响，从低通到低通的变换中，它们的频率都是由 ，区别只是截止频率不同区别只是截止频率不同。如图所示： jPHejHe0 当 ： ：00此时满足变换关系 100jjee1)()(00jjegeg1jjee1)()(jjegeg 当 由0 0变到 时， 相应的也由0 0变到 。根据全通函数全通函数 相相位的变化量为位的变化量为 ，得到满足低通变换的全通函数的阶数全通函数的阶数N=1N=1 。 N所对应的一阶全通函数一阶全通函数应该为： 1111zzzg其中 是实数且1将 及 代入可得到变换所反映的频率变换关系频率变换关系：jez jeu 1jjjeee由上式可以解出 ：221sin21cosarc tg低通原型的截止频率 ，变换后的低通截止频率为 代入上式中可以解出 ： cc22sinsinCCCC变换后得到的新的数字低通滤波器的系统函数是变换后得到的新的数字低通