1、无限长单位脉冲响应无限长单位脉冲响应（IIRIIR）滤波器设计滤波器设计 许多信息处理过程，如信号的过滤，检测、预测等都要用到滤波器，数字滤波器是数字信号处理中使用得最广泛的一种线性系统，是数字信号处理的重要基础。数字滤波器的功能（本质）是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列。实现方法主要有两种：数字信号处理硬件和计算机软件。数字滤波器是具有频率选择性的离散线性系统，数字滤波器的设计步骤大体可分为以下三步：1.1.根据使用要求确定系统需要的各种性能指标（指容限、容差）根据使用要求确定系统需要的各种性能指标（指容限、容差）通常在频域给出滤波器的容限、容差，如图所示。容限
2、容限-指滤波器通带、阻带的最大衰耗和波动范围。容差容差-指通带、阻带容许的最大误差范围。IIR滤波器的设计方法在通带内，幅度响应以误差 逼近于1，即1111jHec在阻带中，幅度响应以误差 而逼近于零，即22jeHp其中，、分别为通带和阻带截止频率。cp2.2.用稳定的因果系统去逼近设计的系统函数用稳定的因果系统去逼近设计的系统函数 系统函数包括：无限长单位脉冲响应的系统函数和有限长单位脉 冲响应的系统函数，可根据容限选择IIR或FIR滤波器。逼近逼近在滤波器设计中逼近是一个重要的环节。逼近就是给定所要求的滤波器性能后，去寻找一个物理可实现的系统函数，使它的频率特性尽可能近似所要求的滤波器特性
3、，也就是指对理想特性进行逼近，最后得到一个因果、稳定并且可实现的传递函数 ，滤波器的设计实际上就是一个数学逼近的问题。zH zH3.3.决定决定 的实现方法的实现方法（包括选择运算结构、运算制式及字长）因为数字滤波器是一种物理可实现的线性时不变系统，所以它是用 一个有限精度的运算去实现系统函数 。选择的结构形式有：级联、并联、直接、横截、频率采样型等结构。本章和第五章主要解决第二个步骤的内容，也就是寻求滤波器传递 函数设计的问题。zH zH说明一下符号表示的含义：DF-DF-数字滤波器 AF-AF-模拟滤波器 单位脉冲响应 单位冲激响应 采样序列 DF的传递函数 AF的传递函数 DF的频率响应
4、 AF的频率响应 nh thanTha zH shajeHjHa数字滤波器四种类型的理想频率特性：数字滤波器四种类型的理想频率特性：NiNiiiinybinxany01)()()(数字滤波器的数学描述：1)差分方程NMzdzcAzbzaZHNiiMiiNiiiMiii一般111110)1()1(1)(2)系统函数简单的数字滤波器1)一阶FIR低通滤波器)1(5.0)(11zzHLP2)一阶FIR高通滤波器)1(21)(11zzHHP3)一阶IIR低通滤波器112112)1()(zzzHLP14)一阶IIR高通滤波器,112)1()(112zzzHHP15)二阶IIR带通滤波器212)1(112
5、1)(zzzzHBP6)二阶IIR带阻滤波器2121BS)1(12121)(zzzzzH4.1 IIR4.1 IIR滤波器设计的特点滤波器设计的特点 nhIIR滤波器的传递函数可以写成N阶的有理函数：0111111NNiiNNiiiiiiiiiiAczHzdzazb z 滤波器设计的核心是求传递函数 ，而 的设计就是确定系数 、，或者确定零点 、极点 ，使得滤波器的 满足给定的性能要求。zH zHia ib icid zH IIR 指单位脉冲响应 为无限长的滤波器，滤波器的 有无限个离散值。一.IIR滤波器的一般设计方法：nh 如果在单位圆内设置一对共轭极点 ，频响在 将有一峰值。0jerZ0
6、 极点 离单位圆愈远，频响在 处的峰值比较平缓。极点 越接近单位圆，频响在 处的峰值就越尖锐。如果通带太窄，可以把极点向原点平移，如果通带太宽，可以将极点向单位圆移动。我们可以通过几次调整单位圆内极点的位置，去逼近符合设计要求的频响。000 同样，如果在单位圆上 处设置一个零点，那么频响在 处会出现传输零点，传输零点可以实现陷波作用。1jez11.累试法 滤波器的幅度特性和相位特性完全由其零点和极点位置所决定。0 零、极点累试法是在频域直接设计，可以完成一些简单的、阶数不高的滤波器的设计。就是计算机优化设计，计算机优化设计的思想是借助计算机，使得所设计的滤波器的频响尽可能逼近理想的频响，即 j
7、eH 逼近jdeH2.最优化设计法最优化设计一般分为两步进行：第一步：选择最优化标准或准则 例如，实际中常常要求滤除叠加在信号上的50HZ工频干扰，我们 可以在 处安排一个零点，就可以1jez,50211HzTT滤除掉 50HZ的工频干扰。011NiNiiiiiHzazbzDF的传递函数 通过改变 的系数 、，分别计算均方误差 ，经过多次迭代运算，寻找一组系数 、，使得均方误差 为最小的这组系数为最优系数，从而完成最优化设计。zHia ib E E ib第二步：进行迭代运算，确定最优系数 21MijwdjwiieHeHE最小选择一种最佳准则，使得 与 的均方误差最小或者最大误差最小。根据最小二
8、乘法准则，要求jeHjdeHia 3.模仿模拟滤波器的设计 因为模拟滤波器的设计目前已经很完善，AF不仅有简单和严格的设计公式，而且它的设计参数也已经表格化，因此，我们可以借助于模拟滤波器设计的成果来设计数字滤波器。在模拟系统中，利用工作参数综合法设计滤波器时，无论低通、高通、带通和带阻滤波器，均是先设计一个低通原型，然后经过某种频率变换完成所要求设计的滤波器。即:利用模拟滤波器设计数字滤波器，首先利用模拟滤波器的现成结果，在S平面设计出符合要求的模拟滤波器的传递函数 ，再通过一定的映射关系，得到数字滤波器的传递函数 。sHa zH 在进行IIR 数字滤波器的设计时，要逼近模拟原型低通滤波器，
9、模拟低通滤波器通常仅考虑幅频特性，习惯上以幅度平方函数来表示模特性。二.最常用的几种模拟原型低通滤波器的逼近方法 jSaaaaasHsHjHjHjH21巴特沃思滤波器（butterworth）最平响应滤波器 巴特沃思低通滤波器的幅度平方函数定义为：NCNaaajjjHjHjH2*2111N 为整数，表示滤波器的阶次；为截止频率。c当=0 时，当时，有 又称为滤波器的3bB带宽（或半功率点）。102jHac212CajH当 时，当 时，又称为滤波器的3bB带宽 0巴特沃思低通滤波器的特点：在 处，即靠近零频处，衰减为0，所以巴特沃思滤波器通带内具 有最大平坦的振幅特性，故得名为最平坦响应滤波器。
10、巴特沃思低通滤波器没有有限零点,零点出现在 处，它属于“全极 点型滤波器”。0 2切比雪夫滤波器(chebychev)切比雪夫滤波器也是一种全极点型滤波器，它的幅度平方函数 CNaCjH22211其中：表示通带波纹大小，是小于 1的正数，越大，波纹越大。为滤波器的截止频率，但 并不是3db带宽）。ccdBHHHHdB2102min2maxminmax101log10log10log203带宽 为N阶切比雪夫多项式，定义为：NC11coscos11CCCCCchNchNNC 切比雪夫低通滤波器的特点：1)通带内等起伏，通带外衰减快；2)由于过渡带较窄，因此相位特性较差。3考尔滤波器(cauer)
11、2211aHjKj 112222nSSSSSBsK 考尔滤波器的特点：1)通带内、外都是等起伏。2)由于过渡带较窄，因此相位特性较差。三S平面到Z平面的映射变换 利用模拟滤波器来设计数字滤波器，就是从已知的模拟滤波器的传递函数 去设计数字滤波器的传递函数 ，即 zH sHa()()aH SH Z通过变换 这种变换归根结底是一个由S平面到Z平面的变换，并且通常是复变函数的映射变换，这种映射变换应该满足两个基本的要求：的频响应该模仿 的频响 即要求 zH sHa 是因果稳定的映射 指 的因果稳定性通过映射后，仍应保持因果 稳定。sHa zH()()aH SH Z通过变换4.2 4.2 脉冲响应不变
12、法脉冲响应不变法 根据容限设计好一个模拟滤波器后，就可对此模拟系统进行模仿。数字滤波器从什么角度去模仿模拟滤波器呢？第一种方法是脉冲响应不变法。1.脉冲响应不变法 模拟 系统 txta thsHaa tythaa dthxthtxaaaaLTI系统特性可以完全由它的冲激响应决定 数字系统 nxn nhzH nynh Kknhkxnhnx 脉冲响应不变法让数字滤波器的脉冲响应和模拟滤波器的脉冲响应在采样点上完全一样。即：nThnha单位脉冲响应不变法的设计思想：使数字滤波器从时域去模仿模拟滤波器。sHa 1L tha采样 nhnTha 变换Z zH2脉冲响应不变法设计的系统的频率响应模拟系统频响
13、数字系统频响 jHa thFTajeHjeZnhZ)(对模拟信号 txa理想采样 txa 1L tha naanTttxtx dtetxsXstaa拉氏变换 ntsadtenTttxnsnTaenTx如果nTxnxa)(nnnanZnTxZnxzX当 时，STeZ sXeXaST 是采样信号的拉氏变换与采样序列Z变换之间的映射关系STeZ 理想采样的频谱 majSaamTjjXTsXjX21采样信号的频谱是采样信号在虚轴上的拉氏变换 是模拟信号频谱 以采 样频率 为周期的周期延拓jXajXaTS2或 jeeZaeXzXjXTjTj理想采样信号的频谱是单位圆上的Z变换，它也等于数字信号的频谱。因
14、此，我们可以得到采样序列 的频谱 与原信号 的频 谱 之间的映射关系为：)(nx txajXajeXmajmTjjXTeX21 上式表明：当采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它实际上是完成了由S平面到Z平面的映射，我们利用以上关系同样可以得到数字滤波器的频响 与模拟滤波器的频响 之间的关系：jeHjHa1212()mmjaaTmH eHjjmHjTTTT 时表达式表明：数字滤波器的频响并不是简单的重现模拟滤波器的频响，为模拟滤波器频响 的周期重复,是以 为周期的周期延拓。jeHjHas需要强调指出，在周期重复模拟滤波器频响的过程中所存在的问题：首先分析一下 和 以及 的关系：j
15、eHjHajHa 注意：我们原本是让 来模拟 ，但实际中由于 ，所以 模拟的是采样信号的频谱 ，的映射关系反映 的是 和 周期延拓后与 之间的对应关系，而并不是 和 之间的关系。jeHjHa nThnhajeHjHaSTeZ jeHjHajHajeHjHa 如果模拟信号不充分带限，这种周期重复就不可避免的存在着频率混叠问题，如下图所示。只有当信号是带限信号，并且满足 、或时，这时数字滤波器的频响 在折叠频率以内才能不失真的重现模拟滤波器的频响 。即jHamS2TmmjeHjHTTjHTeHaaj11 但是任何一个实际的模拟滤波器都不可能完全是带限的，所以不可避免的要出现混叠现象。如果模拟滤波器
16、频响在折叠频率以上衰减的愈快，则混叠失真就愈小。因此，脉冲响应不变法只适用于带限模拟滤波器。对于带阻和高通滤波器，由于它们高端频率不衰减，因此高频分量将完全混叠在低频响应中，从而使整个频响产生失真。如果设计的高通和带阻滤波器要采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，先滤除掉高于折叠频率以上的频带，然后再使用脉冲响应不变法将其转换为数字滤波器。3模拟滤波器的数字化方法脉冲响应不变法要经历以下四个过程：sHa tha1L采样 nhnTha 变换Z zH NiiiNiiiasbsasH101用部分分式展开NiiiSSA11L NitSiatueAthi1对 进行采样，使得数字滤
17、波器的单位脉冲响应序列 tha NinTSiNiTnSianueAnueAnThnhii11（1）（2）（3）4对脉冲响应不变法的修正数字域频响与模拟信号采样后的频谱关系是：majTmjHTeH21上式中频响 与T成反比，通常为了减小脉冲响应不变法的混叠失真，一般选择较高的采样频率，于是采样周期T就会比较小，如果T很小，数字滤波器的增益就会很高，而我们一般希望 的增益能与采样频率无关（与T无关）。因此，在实际应用中，脉冲响应不变法要稍做一下修改。jeHjeH令 ，则 nThTnha NiTiSizeTAzH111majTmjHeH2)(这样，数字滤波器的增益就不会随采样频率的变化而变化。NiT
18、iSzeAizH111比较第（1）步和第（4）步S平面的每一个极点 Z平面是极点 ；与的系数是相同的。因此，只要知道极点 就可以直接得到 ，而不必求解第(2)、(3)步。也就是说我们可以：变换到 zH sHa zHTiSieZ iSiS（4）从模拟域 直接映射到数字域)(nh 变换Z NiiiaSSAsH1 NiTSizeAzHi111下面我们举例说明脉冲响应不变法的应用：例：已知AF的系数函数 用脉冲响应不变法求出相应的数字滤波器的系数函数 sHa22basas zH解：jbasjbassHa2121模拟滤波器在 处有一个零点；在 处有一对共轭极点。aSOjbaSP 11121121zeeT
19、zeeTzHjbTaTjbTaT2211cos21cos1zezbTeTbTezTaTaTa TbjTaTjbTaTaeezeezbTezzcos数字域中 有两个零点 和 ，有一对共轭极点 。可见，脉冲响应不变法对零点没有一一对应的关系。0ZbTeZTacosTbjTaeeZ zH从幅频特性频谱图中可以看出：由于频谱混叠带来明显的失真。很显然，在 处比 在 处下降的慢，这主要是因为“混叠”造成的。混叠的主要原因是由于模拟信号频域的不充分带限。jeHjHaT5脉冲响应不变法的特点 脉冲响应不变法是一种时域设计法，数字滤波器是从时域进行模仿。即 。nThnha但是S平面 Z平面的零点并没有这种单值
20、映射的关系这种时域模仿使得S平面的极点 Z平面的极点为 iS 变换到TiSieZ sHa tha1L采样 nhnTha 变换Z zH脉冲响应不变法设计思想是但我们可以直接从 得到 NiiiaSSAsH1 NiTSizeAzHi111频率坐标是线性变换 即 。因为模拟和数字频率之间的变换关系是线性的，所以可以设 计出满足线性相位的滤波器。只要不发生混叠，变换后数字滤波器的频响 可以不失真的恢复出模拟滤波器的频响。TjHeHajT或注意：数字域的频响特性都是以 为周期的周期函数，其特性对称于折叠频率，因此，由 是 的镜像，通常数字滤波器的频率特性我们只考虑 范围的特性即可。2200 脉冲响应不变法
21、最大的缺点是有频响的混叠，所以只能模仿带限的频响。因此，对衰减特性很好的低通或带通，频响混叠小。对于高通和带阻滤波 器，由于高端频率不衰减，所以要加保护滤波器滤掉高于折叠频率以上的 部分，再用脉冲响应不变法转化成数字滤波器。脉冲响应不变法的最大缺点是产生混叠失真。产生的主要原因是由于S平面到Z平面的变换关系 不是单值的，而是多值的映射关系，由于S平面的多值性造成频率的交叠，产生混叠现象。为了克服这一缺点，我们有必要寻求S平面与Z平面的单值的映射关系。下面我们就介绍另外一种变换方法 双线性变换法。STeZ 4.3 双线性变换法 为了寻求S平面与Z平面的一一对应关系，我们先把S平面压缩到某一个中介
22、的 平面的横带内（宽度从 到 ），然后把此横带再变换到整个Z平面，这时S平面与Z平面就建立了一一对应的单值的映射关系，可以消除混叠现象。如下图所示：1STT1.双线性变换法 为了完成从 这样一个变换过程，我们寻找它们的 映射关系:平面平面平面ZSS1 虚轴 轴平面jS 压缩到横带内轴的平面TTjS,11这个映射可以通过正切变换来实现：21Ttg验证：当 时，当 时，当 时，T12tg001T12tg 的映射关系：平面平面1SS 21Ttg 221222211111TjtheeeeTtgjjTjTjTjTj由 得到 把此关系解析延拓（扩大）到整个S平面令 ，则Sj11SjTSTSTSTSTSTS
23、eeeeeeTSthS11111111222221TSTSeeS1111得到 的映射关系平面平面1SS 映射到 平面S平面Z将标准变换关系 代入 ，从而得到 到 的双线性变换的单值映射关系TSeZ1TSTSeeS1111平面S平面ZSSZZZS111111现在来验证上式中S和Z的变换关系是否满足“映射变换”的两个总要求：j与Z平面的关系S平面虚轴)1(将 、代入双线性映射关系 中，jeZ T11111ZZSjTtgjtgjeeSjj22111 通过验证证实：S平面虚轴()Z平面的单位圆 上，表明整个虚轴是单值的对应于单位圆的一周。因此双线性变换确实消除了脉冲响应不变法所存在的混叠误差，所以它的
24、逼近是良好的。j 映射到jeZ 验证S平面与Z平面是否是因果和稳定的jSjjSSZjs11112222)(1)(1jjZ当 时，S左半平面 Z平面的单位圆内;当 时，S平面虚轴 Z平面的单位圆上;当 时，S右半平面 Z平面的单位圆外。01Z 映射到01Z 映射到01Z 映射到通过验证得到以下结论通过验证得到以下结论：如果模拟滤波器是稳定的，通过双线性变换后，所得到的数字滤波器也 一定是稳定的。如果给定模拟滤波器的传递函数 ，变量S与Z之间有简单的代数关 系，只要用代数置换就可以得到数字滤波器的传递函数。sHa1111)()(1111zzHsHzHaSaZZ说明 在本教材种，引入变换常数C P2
25、47页），主要是为了使模拟滤波器的频率特性与数字滤波器的频率特性在不同频率处有对应关系，通常C=2/T。因此，模拟滤波器和数字滤波器特定频率之间的关系有：2双线性变换的频率响应 S平面到 平面虚轴的映射关系为：1S 21TtgC模拟域与数字域的频率关系为：221tgCTtgC 可以看出：模拟域与数字域的频率变换关系不是线性的，而是非线性关系。221ctgCTtgCc双线性变换的映射关系为：SCSCZZZCS1111数字滤波器的频响)2(2tgCjHjHeHatgCaj 没有频谱混叠 双线性变换的最大优点是避免了频响的混叠效应，因为S平面的虚轴单值 对应着Z平面单位圆的一周。从它们的频率变换关系
26、 中清楚的看出，当 时，终止在折叠频率 处，故不会有高于折叠频率的分量出现，避免了由频 响产生的混叠失真，但双线性变换虽然避免了频响混 叠，但却是以牺牲相位特性为代价换取的。和 的关系是非线性变换2tgC 如果模拟系统是按时延设计的线性相位的模拟滤波器，经过双线性变换后，数字滤波器将得不到这种线性相位特性，所得到的数字系统的相位不再保持线性相位特性。3 3双线性变换的优缺点双线性变换的优缺点 数字频率和模拟频率之间的关系是上式说明，经过 的频率变换后，DF的频响 等于AF的频响 。2tgCjeHjHa例：模拟系统中三阶切比雪夫低通滤波器（Chebyshev)从图中可以看出在靠近截止频率处有频率
27、畸变。模拟滤波器截止频率为 ，对应的数字域频率为 ，而经过双线性变换后数字域的频率是 ，。可见，经过非线性变换以后，在截止频率处的频率特性发生了畸变。因此，如果要求相位和频率是线性关系，就不能采用双线性变换法。cTcccctg12cc 双线性变换适合分段常数幅频响应的滤波器，而实际中大多数滤波器均为分段常数响应滤波器，如典型的低通、高通、带通、带阻滤波器都是具有分段常数响应的滤波器，所谓分段常数响应的滤波器指在通带内要求逼近一个衰减为0的常数，阻带内要求逼近一个衰减为 的常数，所以我们完全可以采用双线性变换。双线性变换所引起的频率畸变只是出现在数字滤波器靠近截只止频率处。这种畸变完全可以通过预
28、畸变加以校正。首先对模拟滤波器的临界频率事先加以预畸变，让 ，也就是使 通过双线性变换后正好映射到所需的位置上去。2cctgC2cctgC 设计简单 双线性变换是IIR滤波器设计中使用最普遍、最有成效的一 种设计方法。由于S与Z之间是简单的代数关系，因此传递函数可 以通过代数置换得到。即 11111111ZZHsHzHaSaZZTcc5.021cCtg2cctgcctg12ccTcc 正好映射到所需的位置上。c4.4 原型变换 模拟滤波器已经形成了许多成熟的设计方法，每种滤波器都有一套准确的计算公式和大量的归一化的设计表格和曲线。所以在模拟滤波器的设计中只要掌握原型变换，就可以通过归一化低通原
29、型的参数，去设计各种实际的低通、高通、带通、带阻滤波器。模拟滤波器的这一套方法可以通过上述讨论的两种变换法而应用于数字滤波器的设计。下面我们讨论各种滤波器设计的问题：一低通滤波器的设计步骤（由模拟低通到数字低通的设计）确定数字滤波器的性能要求（容限），确定各临界频率 的值。c 将数字域 映射到模拟域c2cctgC通过变换关系c把DF的性能要求转换为相应的样本的AF的性能要求，得到模拟滤波器的临界频率c 按照临界频率 的要求设计模拟滤波器的传递函数 通过脉冲响应不变法或者双线性变换，将 c sHa sHa 转换为 zH例例4.14.1 要求用脉冲响应不变法及双线性变换设计一个三阶Butterwo
30、rth低通 滤波器。设采样周期 (即采样频率 )，3dB截止频率sT250ZSkHf4ZCkHf1解：第一步第一步:的设计 Butterworth滤波器的传递函数为：sHa NCNaaajjsHsHsH22111三阶Butterworth滤波器N=3，有2N=6个极点，分布在半径为 的圆上，将S平面分成6等份。为了保证设计的系统是稳定的，的极点应该取S左半平面的极点，因此，S左半平面的传递函数 是c 2sHa sHa 10NKKNCaSSsH其中：21212NKjCKeS1,1,0NK当N=3时，可直接确定三阶Butterworth滤波器的传递函数 为：sHa 32323jeCSjeCSCSC
31、saH322211CCCSSS第二步：第二步：变换法设计滤波器 脉冲响应不变法 将 用部分分式展开为：sHa 2316323163jcSjecjcSjeccScsaH 将 直接变换为数字滤波器的传递函数：sHa 1231163123116311 ZjTcejeTcZjTcejeTcZTceTczH 1231163123116311ZjcejecZjcejecZceczH将上式中共轭复根合并得：22121123cos21623cos26cos231 ZeeZeZZeCCCCCCCCzH 其中数字滤波器的截止频率为 可以看出 只与临界频率 与采样频率 的相对值 有关，而与 、的绝对大小无关，所以只
32、要 一定，所设计的数字滤波器则具有同一个传递函数，这个结论适用于所有的数字滤波器。SCCCCffTfT22 zHCfSfSCffSCffSfCf 例例4.24.2 1kHz 4kHz 则在数字域中的 1MHz 4MHz 设计是相同的。CfSf41SCffCfSf 将 代入 中：5.02TfCC zH 21112079.01905.015541.0571.12079.01571.1ZZZZzH 以上得到的是数字三阶Butterworth滤波器的传递函数 ，采用一个一阶基本节和一个二阶基本节并联的结构实现比较方便。zH除此之外，我们还可以采用另外一种方法计算 从上面的实例中知道，当求模拟滤波器的传
33、递函数 时，只是一个符号而已，它并不进行计算，乘以T后变成 ，最终传递函数只取决于数字域的参数 。因此只需求出一个相对於采样频率归一化的“样本”传递函数即可。这种计算方法可以直接以 作为模拟滤波器的临界频率 来计算(让 ），所得到的 是采样频率归一化后的模拟样本传递函数。zH sHaccCCCcCc sHa CCCSSSsaH22211 322211CCCSSSccsaH得到的 是对采样频率归一化的模拟样本传递函数，再把它展开成部分分式形式：sHa 111100dSAdSAdSAsHa 111110110111ZeAZeAZeAzHddd T=1已归一化双线性变换法设计步骤 i)首先确定数字域
34、临界频率C5.02TfCCii)确定预畸的模拟滤波器的临界频率 由于双线性变换法频率变换的非线性会引起频率畸变，所以要预畸变。为了使数字滤波器的临界频率为 ，模拟滤波器的临界频率应为：C5.0C142tgtgCC5.021cCtg刚好映射到我们所要求的位置上。得到频率归一化的模拟滤波器的样本传递函数 122123SSSsHaiv）最后将双线性变换的关系式代入 ，得到数字滤波器的传递函数 sHa zH 232111333121111111ZZZZZZHsHzHaSaZZ注意：这里所采用的模拟滤波器 并不是数字滤波器所要模仿的截止频率为 的那个实际滤波器，只是模拟滤波器的样本传递函数，是归一化的低
35、通原型，它是由模拟滤波器的低通原型到数字滤波器变换的一个中间变换函数。上面介绍了如何设计一个数字低通滤波器。sHaZCkHf1 sHaiii）将 代入 322211CCCSSSsaH1C sHa 当我们需要设计高通、带通、带阻数字滤波器时，有两种方法：第一种：首先设计一个相应的高通、带通、带阻模拟滤波器，然后再通过脉冲响应不变法或双线性变换法转换为数字滤波器，即的频带变换模拟到模拟变换法将其数字化变法或者双线性通过脉冲响应不模拟低通原型 模拟的高通、带通、带阻 数字的高通、带通、带阻模拟低通原型 数字的高通、带通、带阻频带变换第二种：直接利用模拟滤波的低通原型，通过一定的频率变换关系，一步完成
36、各种数字滤波器的设计，即上述两种方法的根本区别在于第一种方法是在模拟域进行频带变换，第二种方法是在数字域进行频带变换。第一种方法的实现步骤是：首先确定临界频率 模拟滤波器的临界频率 ；完成模拟高通、带通、带阻滤波器的设计；最后通过脉冲响应不变法或者双线性变换法将设计好的 转换成 即可。由于第一种方法和我们讨论的低通滤波器的设计完全一样，因此不再讨论。这里着重介绍第二种方法，因第二种方法更简捷方便，故得到普遍采用。由于 脉冲响应不变法对于高通、带阻滤波器都不能直接采用，或者只能加了保护滤 波器以后使用，因此第二种方法仅讨论采用双线性变换法直接进行频率变换的 情况。c sHa zHc转换成相应的
37、二高通变换 在模拟滤波器的高通设计中，低通至高通的变换就是S变量的倒置变换，只要把双线性变换中 即可。SS1低通为 高通为SSZZZCS111111SSZZZCS111111 11111111zzCHsHzHaLCSaLhZZ将 代入上式复变量S中jeZ 倒置代换jjctgCeeCSjj211数字域与模拟域的频率映射关系为 2ctgC 由右图中看出：当 ，Z平面是映射在 上当 ，Z平面是映射在 上0 映射为1Z 映射为01Z 说明双线性变换中S变量的倒置变换并不会影响变换后的稳定条件，而且 轴仍然映射到Z平面的单位圆上，只是方向颠倒了。通过这样的变换，就 可以直接把模拟低通变换为数字高通，如右
38、图所示。j 高通变换的计算步骤和低通变换是一样的，只是在确定模拟原型预畸变的临界频率时，应该采用2KKctgC 例例4.34.3已知采样频率 ,ZkHf10sT100ZSkHf5.2ZkH5设计一个三阶切比雪夫高通数字滤波器，其通过频率 （但不必考虑 以上的频率分量),通带内损耗不大于1dB。解：切比雪夫低通原型的模平方函数为CNaCjH22211通带损耗 时，当N=3时，通过查表可得传递函数为dBAP15089.0 322339883.0238.14913.04913.0SSSsHCCCCa现在我们来求预畸的截止频率 ，首先确定数字域的截止频率 ，cc ，于是有 ，代入上式得：5.02TfC
39、C142cgctgCC 329883.0238.14913.04913.0SSSsHa 11111111ZZHsHzHLaSLahZZ 3213212041.06043.03432.013311321.0ZZZZZZzH三.带通变换根据带通滤波器变换的要求 S平面 Z平面 0S 映射到0jeZ0（是带通中心频率）S 映射到1Z模拟低通的 数字域的0 映射到0 映射到,0 满足以上带通变换要求的双线性变换为：11cos211)(20200ZZZZZeZeZzFSjj我们可以验证F(z)是否满足带通变换，F(z)有两个零点和两个极点：将零点 代入F(z)，S0 将极点 代入F(z)，0jeZ1Zj
40、S说明双线性变换满足带通变换的要求。11cos2202ZZZS 11cos220211cos2202ZZZHsHzHaLPZZZSaLPBP直接将 代入模拟低通的传递函数 中，求得：sHa将变换关系 、代入F(z)=S中有：jeZ jSjjeeeeeeejjjjjjjjsincoscoscos211cos200202 数字域 与 模拟域 之间 的映射关系为：sincoscos0由右图看出：在 点正好映射在 上 当 时，正好映射在 上，当 时，正好映射在 上，因而频率变换关系完全满足带通变换的要求。000低通幅度响应到带通幅度响应的变换如右图所示带通变换的设计步骤 在设计带通时，通常给出上下边带
41、的截止频率 、和阶数N。首先根据已知的 和 计算出数字带通的中心频率 和模拟低通的截止频率 。12120c1101sincoscos2202sincoscos110sincoscosC2cos2cossinsinsincos212121210模拟滤波器中 和 是一对镜像频率，即 ，而且 也是模拟滤波 器的截止频率 ，因此由上面这两个关系式可以解出 和12211c01cii）求归一化的模拟低通原型iii）把满足带通变换的双线性变换代入模拟低通 中。sHa sHai）求模拟低通的截止频率 c 11cos220211cos2202ZZZHsHzHZZZS 例例4.44.4：已知采样频率为 ，要求设计
42、一个三阶巴特沃兹带通滤波器，其上下边带的3dB截止频率分别为 和 。ZSkHf100sT100ZkHf5.371ZkHf5.12275.0211Tf25.0222Tf解：首先确定数字频域的上下边带截止频率为：代入 中，21210sinsinsincos075.0sin25.0sin75.025.0sincos05.00175.0sin75.0cos5.0cossincoscos110C再带入 中c将求出的模拟低通的截止频率 将代入三阶巴特沃兹的传递函数中：1C 122123SSSsHa 11cos220211cos2202ZZZHsHzHZZZS4642333121ZZZZ三.带阻变换 将带通
43、滤波器的频率变换关系倒置一下就可以得到带阻滤波器，整个设计方法、步骤和带通滤波器是类似的。因此，带阻变换的关系为：1cos21022ZZZS0coscossin（有 ，分母有变化）j4.5 Z平面变换法 前面讨论的原型变换均是由模拟低通原型来设计各种数字滤波器，这种原型变换的设计方法同样也可以直接在数字域中进行。如果我们已经有一个数字滤波器的低通原型函数 ，就可以直接在数字域中进行变换，设计出所需要的数字滤波器的 。也就是：zHP zH数字域低通原型通过一定变换关系滤波器带阻带通高通低通数字BSBPHPLP即：zHP通过映射关系 zH 为了区分变换前后两个不同的Z平面，我们规定：把变换前的Z平
44、面定义为U平面，变换后的Z平面定义为Z平面。U平面的单位圆用 表示，Z平面的单位圆用 表示。jeje从U平面到Z平面的映射关系用变换函数 表达为：11zgug 因为实际的传递函数都是以负幂形式出现，所以 和 都是用负幂次方表示。zuU平面 到 Z平面的映射关系实际上是复平面的保角映射。现在我们来分析一下函数 的一般特性：1zgI）变换以后的传递函数应该保持稳定性不变首先，我们希望变换以后的传递函数应该保持稳定性不变，因此要求U平面的单位圆 必须映射到 Z平面的单位圆 。jejejjjjeegege其中：是 的相位函数，由上式可得到：jeg1jeg说明函数 在单位圆上的幅度必须恒等于1，这种函数
45、称为全通函数。1zg全通函数任何频率的一个信号通过一个网络，输出幅度不变，只是相 位变化，即只对相位进行校正，这种网络称为全通网络，表 示全通网络的网络函数称为全通函数。Niiiizzzg1111全通函数的表达形式为：全通函数 的基本特性：N 全通函数的阶数1zg 是极点，可以是实数，也可以是共轭复数。且 ，因为一 个稳定系统的极点必须在单位圆内。是零点，所有的零点都是极点的共轭的倒数。ii1ii1可以证明，Z平面当 由0变化到 时，其相位函数 由0变到 ，因此 可以由相位的变化来确定阶数N。N 只要选择合适的N和 就可以得到各类需要的变换。下面我们具体讨论各种原型变换：i 一.低通 低通 j
46、PeHjeH从U平面 Z平面到 及 均是低通滤波器的系统函数，从低通到低通的变换中，它们的频率都是由 ，区别只是截止频率不同已。如图所示：jPeHjeH0 当 ：00此时满足变换关系 100jjee1)()(00jjegeg1jjee1)()(jjegeg 当 由0变到 时，相应的也由0变到 。根据全通函数 相位的变化量为 ，得到满足低通变换的全通函数的阶数N等于1。N所对应的一阶全通函数应该为：1111zzzg其中 是实数且1将 及 代入上式，可以得到这个变换所反映的频率变换关系：jez jeujjjeee1由上式可以解出 ：cos12sin122arctg如果低通原型的截止频率 ，变换后的
47、低通截止频率为 ，则代入上式中可以解出 ：cc2sin2sinCCCC变换后得到的新的数字低通滤波器的系统函数是 111111zzzguHzHzguP下面我们讨论 、和 的关系，由频率关系式cos12sin1221ctg 、和 的关系如图所示：时，此变换代表的是频率压缩 时，此变换代表的是频率扩展 由上可知U平面到Z平面的频率变换关系是非线性的。0cc0cc 二.低通 高通 如将 变换为 ，相当于将单位圆上的频响旋转了一个 的角度，进行了一个旋转变换，这时 平面上的低通就变换为相应的高通了。zzz所以只要将 中的 以 替代，就完成了低通到高通的原型变换。1111zzzg1z1 z低通到高通的原
48、型变换为：11111zzzgu在低通到高通的变换中，当低通由通带走向阻带，即 由 时，高通是由通带走向阻带，是由020或 根据全通函数相位 的变化量为 ，则全通函数的阶数N必须为1。所以有 N2cos2cosCCCC三.低通 带通低通到带通的变换应该满足以下两点要求：若带通的中心频率为 ，则对应着低通原型的通带中心，即 点。当带通滤波器由通带走向阻带时，频率 由 ，对应的低通原型的频率 是由 ;同样，当带通滤波器由阻带走向通带时，频率 由 时，对应的低通原型的频率 是由 。000000）（20相当于可以看出，可以看出，低通到带通的变换，由于 变换后在 之间 或 之间形成两个通带，也就是变换函数
49、需要 在单位圆上转两次，因此，全通函数的阶数N应等于2。20由于带通滤波器的变换应满足下面条件：1)1()(11)1()(111111111gzgugzguzz因此，低通到带阻的变换全通函数应取“”号。【例例4 45 5】已知一阶巴特沃思数字低通的传递函数 ，截止频率 。要求设计一阶巴特沃思数字带通滤波器，其上、下截止频率分别 为：。1131)(zzzHP3.0c4.0,6.012解解：首先确定参数，参数的确定请参看P109页表4102)4.06.0(cos2)4.06.0(cos2)(cos2)(cos1212569.1325.0510.0)23.0(2)4.06.0(2212 tgctgt
50、gctgkC 11211111211121211111zkkzkkkkzkkzzzzzzg二阶全通函数为：221.3663.1779.0779.01221.0221.031221.0221.01)()(222222)(11zzzzzzuHzHzguL 系统函数在 处有一对共轭极点，处有两个零点，如图所示。当 时带通中心频率处的频响为：779.0jZ1Z5.0jeZ 1221.3663.1779.0779.0)(5.0jeH11211111212121zkkzkkkkzkkzzg)(1zg把以上参数代入 得：012kk221.01569.11569.111kk四.低通 带阻 低通到带阻的变换，由
