《无限脉冲响应数字滤波器的设计ppt培训课件》由会员分享，可在线阅读，更多相关《无限脉冲响应数字滤波器的设计ppt培训课件（54页珍藏版）》请在金锄头文库上搜索。
1、无限脉冲响应数字滤波器的设计,学习目标 理解数字滤波器的基本概念 掌握冲激响应不变法 掌握双线性变换法 掌握Butterworth低通滤波器特点 理解利用模拟滤波器设计IIR数字滤波器的设计过程,7.1 数字滤波器的基本概念,1. 数字滤波器的分类数字滤波器从实现的网络结构或者从单位脉冲响应分类，可以分成无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。它们的系统函数分别为：,理想低通、高通、带通、带阻滤波器幅度特性,按滤波器的功能分类：低通高通带通带阻,P316317,2数字滤波器的技术要求我们通常用的数字滤波器一般属于选频滤波器。假设数字滤波器的传输函数H(e j)用下式表示：,
2、低通滤波器的技术要求,P348 图7.1,通带内和阻带内允许的衰减一般用dB数表示,通带内允许的最大衰减,如将|H(ej0)|归一化为1，则：,阻带内允许的最小衰减,P349 图7.2,3. 数字滤波器设计方法概述IIR滤波器和FIR滤波器的设计方法是很不相同的。IIR滤波器设计方法有两类：经常用的一类设计方法是借助于模拟滤波器的设计方法进行的。设计步骤是：先设计模拟滤波器得到传输函数Ha(s)，然后将Ha(s)按某种方法转换成数字滤波器的系统函数H(z)。另一类设计方法是直接在时域或频域中进行设计，设计时需要通过计算机作辅助设计。,P351 图7.3,7.5 模拟滤波器的设计,模拟滤波器的理
3、论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。,设计方法: 先设计模拟低通滤波器，再通过频率变换的方法将低通滤波器转换为其它类型的滤波器。,各种理想滤波器的幅频特性,1.模拟低通滤波器的设计指标及逼近方法,模拟低通滤波器的设计指标有Rp和As一般用dB数表示。,对于单调下降的幅度特性，可表示成：,滤波器的技术指标给定后，需要设计一个传输函数Ha(s)，希望其幅度平方函数满足给
4、定的指标Rp和As，一般滤波器的单位冲激响应为实数，因此,要求是稳定的，即极点落在s平面的左半平面。,要求极点落在s平面的右半平面。,2.巴特沃斯低通滤波器的设计方法巴特沃斯低通滤波器的幅度平方函数|Ha(j)|2用下式表示：,图7.4 巴特沃斯幅度特性和N的关系,将幅度平方函数|Ha(j)|2写成s的函数：,此式表明幅度平方函数有2N个极点，极点sk用下式表示：,三阶巴特沃斯滤波器极点分布,为得到稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(-s)。 Ha(s)的表示式为,设N=3，极点有6个，它们分别为,取s平面左半平面的极点s0,s1
5、,s2组成Ha(s)：,由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。这里采用对3dB截止频率c归一化为1，归一化后的Han(s)表示为,查P370 表7.2 7.3 7.4,令,则N由下式表示：,关于3dB截止频率c，如果技术指标中没有给出，可以按照下式求出,总结以上，低通巴特沃斯滤波器的设计步骤如下：(1)根据技术指标p,Rp,st和As，求出滤波器的阶数N。(2)确定滤波器极点。 (3)确定系统函数Ha(s)。,表 巴特沃斯归一化低通滤波器参数,查P370 表7.2 7.3 7.4,例 已知通带截止频率fc=5kHz，通带最大衰减Rp=2dB，阻带截止频率fs=12kHz
6、，阻带最小衰减As=30dB，按照以上技术指标设计巴特沃斯低通滤波器。解 (1) 求阶数N,(2) 求极点,(3) 求系统函数,上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。这里不如直接查表简单，由N=5，直接查表得到：极点：-0.3090j0.9511,-0.8090j0.5878; -1.0000,式中 a0=1.0000,a1=3.2361,a2=5.2361,a3=5.2361,a4=3.2361,对于实际滤波器，其系统函数应将上式中的S换为,7.7 冲激响应不变法（脉冲响应不变法）,为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求：(1)
7、H(z)的频响应该模仿Ha(s)的频响，即：s平面的虚轴 映射到z平面的单位圆 上。(2) Ha(s)的因果稳定性通过映射后H(z)仍应保持，即：s平面的左半平面应该映射到Z平面的单位圆内。,方法：,先设计模拟低通滤波器的传输函数,设模拟滤波器的传输函数为Ha(s),其单位冲激响应是ha(t),将Ha(s)用部分分式表示：,将Ha(s)进行逆拉氏变换得到ha(t)：,对ha(t)进行等间隔采样，采样间隔为T，得到：,对上式进行Z变换，得到数字滤波器的系统函数H(z)：,由上式可以看出：冲激响应不变法将s平面的极点 映射到z平面的映像极点为 。可以证明，用冲激响应不变法将模拟滤波器 转换成数字滤
8、波器H(z)时，整个s平面映射到Z平面的映射关系为：,表明：数字频率与模拟频率之间是线性关系，这是冲激响应不变法的优点之一。=0, r=1 0, r1,图7-15 z=esT,s平面与z平面之间的映射关系,数字滤波器的频率响应和模拟滤波器的频率响应间的关系为：,即，数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。,图 脉冲响应不变法的频率混叠现象,只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，即：,才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不产生混叠失真，即,但是，任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生周期延拓分量的频谱交叠
9、，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。,混叠：,（1）一定存在，且在 的奇数倍处； （2）高通，带阻不适合用。,缺点：,频谱混叠,优点：,脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好；,例 已知模拟滤波器的传输函数Ha(s)为用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。,极点为,H(z)的极点为,解 首先将Ha(s)写成部分分式：,设T=1s时用H1(z)表示，T=0.1s时用H2(z)表示，则,高频抑制不够充分，已发生混叠。,7.9 双线性变换法,、变换原理,脉冲响应
10、不变法的主要缺点是产生频率响应的混叠失真。这是因为从s平面到平面是多值的映射关系所造成的。,为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-/T/T之间，再用z=esT转换到Z平面上。,第一步先将整个s平面压缩映射到s1平面的-/T/T一条横带里；,这样就使s平面与Z平面建立了一一对应的单值映射关系，消除了多值变换性，也就消除了频谱混叠现象 。,模拟频率和数字频率之间的关系,式中T仍是采样间隔，当1从-/T经过0变化到/T时，则由-经过0变化到+，实现了s平面上整个虚轴完全压缩到s1平面上虚轴的/T之间的转换。,为了将s平面的整个虚轴j压缩到s1平面j1轴上的-/
11、T到/T段上，可以通过以下的正切变换实现:,再将s1平面通过以下标准变换关系映射到Z平面：,从而得到s平面和Z平面的单值映射关系为：,将此关系解析延拓到整个s平面和s1平面，令j=s，j1=s1，则得,这种变换是s平面与Z平面之间的单值映射关系，都是两个线性函数之比，因此称为双线性变换。,如果为了使AF和DF的某一频率有对应关系，可引入待定常数c，,（1）令s=j,z=e j，有,双线性变换法的频率变换关系,二、逼近的情况,（2）令s=+j,则：,可见，当0时，|z|1。 即：s的左半平面映射到Z平面的单位圆内，s的右半平面映射到Z平面的单位圆外，s的虚轴映射到Z平面的单位圆上。因此，稳定的模
12、拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。,即：s平面上与Z平面的成非线性的正切关系，从而消除了频谱混叠现象。,三、优缺点分析：,双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为s平面与Z平面是单值的一一对应关系。s平面整个j轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。,s平面上与Z平面的成非线性的正切关系，,在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率=处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。,双
13、线性变换的这个特点是靠频率的严重非线性关系而得到的。由于这种频率之间的非线性变换关系，就产生了新的问题：,首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；,其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变。,对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸变来加以校正
14、。也就是将临界模拟频率事先加以畸变， 然后经变换后正好映射到所需要的数字频率上。,双线性变换法比起脉冲响应不变法来，在设计和运算上也比较直接和简单。由于双线性变换法中，s到z之间的变换是简单的代数关系，所以可以直接将变换关系式代入到模拟系统传递函数，得到数字滤波器的系统函数，即：,频率响应也可用直接代换的方法得到,需要特别强调的是，若模拟滤波器Ha(s)为低通滤波器，应用,变换得到的数字滤波器H(z)也是低通滤波器；,也是高通滤波器；若为带通、带阻滤波器也是如此。,在IIR数字滤波器的设计中，当强调模仿滤波器的瞬态响应时，采用脉冲响应不变法较好；而在其余情况下，大多采用双线性变换法。,例 试分
15、别用脉冲响应不变法和双线性变换法将下图所示的RC低通滤波器转换成数字滤波器。解 首先按照下图写出该滤波器的传输函数Ha(s)为,利用脉冲响应不变法转换，数字滤波器的系统函 数H1(z)为,利用双线性变换法转换，数字滤波器的系统函数H2(z)为,H1(z)和H2(z)的网络结构分别如图所示。,H1(z)和H2(z)的网络结构(a)H1(z); (b)H2(z),数字滤波器H1(z)和H2(z)的幅频特性,利用模拟滤波器设计IIR数字低通滤波器的步骤(1)确定数字低通滤波器的技术指标：通带截止频率c、通带衰减RP、阻带截止频率st、阻带衰减As。(2)将数字低通滤波器的技术指标转换成模拟低通滤波器
16、的技术指标。,采用脉冲响应不变法和双线性变换法，其边界频率的转换关系分别为,(3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。(4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。,例 设计低通数字滤波器，要求在通带内频率低于0.2rad时，容许幅度误差在1dB以内；在频率0.3到之间的阻带衰减大于15dB。指定模拟滤波器采用巴特沃斯低通滤波器。试分别用脉冲响应不变法和双线性变换法设计滤波器。,解 (1) 用脉冲响应不变法设计数字低通滤波器。 数字低通的技术指标为 p=0.2rad,Rp=1dB;st=0.3rad,As=15dB 模拟低通的技术指标为T=1s,p=0.2rad/s,Rp=1dB;st=0.3rad/s,As=15dB,设计巴特沃斯低通滤波器。先计算阶数N。,
