iir数字滤波器的设计.pptx

4.1.1 数字滤波器的分类数字滤波器的分类数字滤波器主要有以下三种分类方法：（1）根据滤波器的性能来分，可以分为经典数字滤波器和现代数字滤波器。经典数字滤波器，即一般滤波器，特点是有用信号和干扰信号各占不同的频带，通过一个合适的滤波器选择出有用的频率信号滤除干扰信号，因此也称为选频滤波器。现代数字滤波器的特点是针对信号和干扰的频带相互重叠时，按照随机信号内部的一些统计分布规律，从干扰中提取有用信号。例如：维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。（2）根据数字滤波器的功能来分，可以分为：低通、高通、带通、带阻、多带数字滤波器。（3）根据实现的网络结构或者从单位脉冲响应来分，可以分成：IIR、FIR滤波器。1.数字滤波器的数学描述数字滤波器的数学描述1)差分方程2)系统函数4.1.2数字滤波器的技术指标式中：称为幅频特性，称为相频特性。下面以低通IIR数字滤波器为例，介绍数字滤波器的技术指标。假设数字滤波器的传输函数用下式表示：通带内波动阻带内最小衰减问：其他类型滤波器的频率特性及技术指标？4.1.3 数字滤波器设计方法概述数字滤波器设计方法概述1）按照实际需要确定滤波器的类型和性能指标要求；2）用一个因果稳定系统的H(z)或h(n)去逼近这个性能要求；3）用一个有限精度的运算去实现这个系统函数。包括选择运算结构：如级联型、并联型、卷积型、频率采样型以及快速卷积（FFT）型等。数字滤波器的设计思想都是采用逼近方法。设计的基本步骤是：IIR数字滤波器设计方法可以分为两类：直接设计法和间接设计法。（1）直接设计法直接设计法是直接在频域或者时域进行逼近设计，有零极点累试法、频域逼近法、时域逼近法，通常必须借助计算机进行优化设计，因此也称为最优化设计方法。具体步骤分两步：1)确定一种最优准则，如：最小均方误差准则，即使设计出的实际频率响应的幅度特性与所要求的理想频率响应的均方误差最小。2)在此最佳准则下，通过不断地迭代运算求滤波器的系数，直到满足要求为止。此外还有其他多种误差最小准则。由于计算机技术的发展，最优化设计方法的使用也逐渐增多；但是，由于涉及到最优化设计算法，而且目前已经有商业设计程序可以使用，因此本章不介绍。（2）间接设计法间接设计法是借助于模拟滤波器的设计方法进行的，其设计步骤是：先设计模拟滤波器得到传输函数，然后将按某种方法转换成数字滤波器的系统函数。由于模拟滤波器的设计方法已经发展得很成熟，而且模拟滤波器有简单而严格的设计公式，设计起来方便、准确；因此可将这些理论推广应用到数字域，作为设计数字滤波器的工具。本章后续内容将介绍这种设计方法。模拟滤波器根据幅度特性可以分为低通、高通、带通和带阻滤波器，它们的理想特性如图所示。4.2 4.2 模拟滤波器模拟滤波器的的设计设计模拟滤波器的设计就是根据一组设计规范设计模拟系统函数Ha(s)，使其逼近某个理想滤波器特性。通常都是先设计低通滤波器，再通过频率变换将低通变换为希望类型的滤波器；因此，下面先介绍模拟低通滤波器的技术指标和逼近方法，然后再介绍模拟滤波器的频率变换。4.2.1 模拟低通滤波器的技术指标及逼近方法模拟低通滤波器的技术指标及逼近方法模拟低通滤波器的设计指标用4个参数描述。为通带截止频率；为幅度下降到0.707，也就是下降了3dB时的频率，因此通常被称为3dB截止频率；为阻带截止频率。对于一般滤波器的单位脉冲响应为实数的因果系统有：定义振幅平方函数又S=j，2=-S2A(2)=A(-S2)|S=j问题：如何由A(-S2)|Ha(S)？对于给定的A(-S2)，先在S复平面上标出A(-S2)的极点和零点，由于A(-S2)的极点和零点总是“成对出现”，且对称于S平面的实轴和虚轴；故选用选用A(-S2)在在S左半平面的极点作为左半平面的极点作为Ha(s)的极点，零点的极点，零点可选用任一半可选用任一半。可见只要给定幅度平方函数模型，模拟滤波器就可以设计出来。4.2.2巴特沃兹（巴特沃兹（Butterworth）滤波器逼近）滤波器逼近特点：具有通带内最大平坦的振幅特性，且随f，幅频特性单调。巴特沃兹滤波器(Butterworth滤波器)幅度平方函数：巴特沃兹滤波器幅度平方函数图1中，N增加，通带和阻带的近似性越好，过渡带越陡。通带内，分母/c1，(/c)2N1,增加，A(2)快速减小。=c,幅度衰减相当于3db衰减点。可以证明：理想幅度响应及其逼近响应的差值，与=0处的导数等于0的阶数有关，应尽量使导数等于0的阶次高些。N增加，通带内频率响应变得更为平坦，阻带内的衰减更大；过渡带内频响更趋于斜率为过渡带内频响更趋于斜率为-6NdB/倍频程的渐倍频程的渐近线近线。振幅平方函数的极点：可见，Butterworth滤波器的振幅平方函数有2N个个极点极点，它们均匀对称地分布在|S|=c的圆周上。例：为N=3阶BF振幅平方函数的极点分布，如图。3阶A(-S2)的极点分布考虑到系统的稳定性，可知DF的系统函数是由S平平面面左左半半部部分分的的极极点点（SP0，SP1，SP2）组成的，它们分别为：称为归一化的极点，它只与N有关；因此，当N确定时，归一化的极点可以事先算出并建立一个参数表。可以写成以下分母多项式形式当N确定时，通过查表就可以得到滤波器的归一化原型，例如：当N=3时，查表得：设计的巴特沃兹模拟低通滤波器传递函数为：问题：如何确定滤波器的问题：如何确定滤波器的N、c？根据滤波器的技术指标来确定的。查表得H(p)，则：【例4.2.1】已知模拟低通滤波器通带截止频率fp=0.5Hz，通带最大衰减=2dB，阻带截止频率fr=1.2Hz，阻带最小衰减=30dB，试采用巴特沃斯模型设计该模拟低通滤波器。解：利用上述的方法设计，设计步骤为：(1)确定模拟低通滤波器的阶数N。（2）查表得到5阶归一化低通原型滤波器：b0=1.0000，b1=3.2361，b2=5.2361，b3=5.2361，b4=3.2361取N=5。（3）求3dB截止频率（4）去归一化，得到所设计的低通滤波器：MATLAB提供三个函数用于ButterWorth模拟滤波器的设计。（1）阶数及截止频率的求取函数buttordbuttord函数的语句格式为：n,Wn=buttord(Wp,Wr,Rp,Rr,s)其中：Wp为通带截止频率、Wr为阻带截止频率，单位是rad/s；Rp为通带最大衰减、Rr为阻带最小衰减，单位是dB；选项s表示设计模拟滤波器，省略此参数为设计数字滤波器。当设计带通、带阻滤波器时，Wp，Wr为2元数组。函数的返回值Wn为截止频率，n为满足技术指标的滤波器的最小阶数。（2）设计函数butter有两种语句格式，当设计低通或者带通（带通的Wn为2元数组）滤波器时，语句格式为：b,a=butter(n,Wn,s)当设计其他类型滤波器时，采用ftype说明滤波器的类型，语句格式为：b,a=butter(n,Wn,ftype,s)式中：ftype可以为high,low,stop，分别表示设计的滤波器为高通、低通、带阻滤波器。函数的返回参数b,a为模拟滤波器的传递函数模型：（3）归一化模拟低通原型设计函数z,p,k=buttap(n)，n为阶数，z为零点，p为极点，k为放大系数。【例4.2.2】例4.2.1的MATLAB设计代码为：方法（1）：利用模拟低通原型设计，程序如下：fp=0.5*2*pi;ap=2;fr=1.2*2*pi;as=30;n,Wn=buttord(fp,fr,ap,as,s)z,p,k=buttap(n);%低通原型设计b0=k*real(poly(z);%poly函数是把多项式根转换为多项式系数，real是复数取实部运算a0=real(poly(p);b,a=lp2lp(b0,a0,Wn)；h,f=freqs(b,a);mag=abs(h);mag=20*log10(mag);phase=angle(h);subplot(2,1,1);plot(f/(2*pi),mag);title(N=5ButterworthLowpassFilter);axis(01.5-350);xlabel(f(Hz);ylabel(幅度(dB);grid;subplot(2,1,2);plot(f/(2*pi),phase);axis(01.5-44);xlabel(f(Hz);ylabel(相位(ard);方法（2）：直接设计模拟滤波器，修改程序如下：fp=0.5*2*pi;ap=2;fr=1.2*2*pi;as=30;n,Wn=buttord(fp,fr,ap,as,s)b,a=butter(n,Wn,s)h,f=freqs(b,a);程序运行结果：（1）滤波器的阶数N=5，与上述计算结果相同。（2）截止频率Wn=3.7792rad/sec=0.6015Hz，与计算结果相同。（3）程序运行可得：b=770.9440a=1.000012.229974.7850282.6305660.1398770.9440则所设计的滤波器的传递函数为：（4）滤波器的特性分析结果可以看出在通带截止频率0.5Hz处的衰减小于2dB，阻带截止频率1.2Hz处阻带最小衰减大于30dB，满足指标要求。相位特性曲线在通带内相位特性接近线性。切比雪夫滤波器具有波动性。它有两种类型：（1）切比雪夫型，在通带中具有等波动响应；（2）切比雪夫型，在阻带中具有等波动响应。切比雪夫滤波器在通带范围内是等幅起伏的，所以同同样样的的通通带带衰衰减减，其其阶阶数数较较巴巴特特沃沃兹兹滤滤波波器器要要小小。可根据需要对通带内允许的衰减量（波动范围）提出要求，如：要求波动范围小于1dB等。4.2.3切比雪夫滤波器切比雪夫滤波器切比雪夫型滤波器的振幅平方函数为VN（x）为N阶切比雪夫多项式，定义为式中：表示|Ha(j)|波动范围的参数；当 =0时，c，随/c，0(迅速趋于零)N为偶数，N为奇数，通带内变化范围切比雪夫滤波器的振幅平方特性切比雪夫型的平方幅度响应为：切比雪夫型滤波器通带幅频特性是单调的，而阻带是等波动的。这种滤波器在s平面上既有极点，又有零点。切比雪夫模拟滤波器的设计计算比较复杂，因此只介绍用Matlab设计。MATLAB提供多个函数用于切比雪夫模拟滤波器的设计，下面以型为例进行介绍；对于切比雪夫型滤波器设计，只要把切比雪夫型设计函数中的“1”改为“2”，即为切比雪夫型滤波器设计函数。（1）阶数求取函数cheb1ord，语句格式为：n,Wn=cheb1ord(Wp,Wr,Rp,Rr,s)其中：Wp、Wr为角频率，单位是rad/s。当所设计的滤波器为带通、带阻滤波器时，Wp，Wr为2元数组；返回值Wn为截止频率，n为满足指标滤波器的最小阶数。（2）设计函数cheby1，语句格式有两种：设计低通或者带通（带通的Wn为2元数组）滤波器时格式：b,a=cheby1(n,Rp,Wn,s)设计的是其他类型的滤波器时格式为：b,a=cheby1(n,Rp,Wn,ftype,s)式中：ftype可以为high,low,stop，分别表示设计的滤波器为高通、低通、带阻滤波器。（3）归一化模拟低通原型设计函数z,p,k=cheb1ap(n,Rp)，n为阶数，Rp为通带最大衰减，z为零点，p为极点，k为放大系数。【例4.2.3】例4.2.1的切比雪夫型设计。方法（1）：先设计低通原型，再转换成所要设计的模拟低通滤波器。程序代码：fp=2*pi*0.5;Ap=2;fr=2*pi*1.2;Ar=30;N,Wn=cheb1ord(fp,fr,Ap,Ar,s)z,p,k=cheb1ap(N,Ap);%低通原型A,B,C,D=zp2ss(z,p,k);%转换为状态空间表达式At,Bt,Ct,Dt=lp2lp(A,B,C,D,Wn);%低通原型转换b,a=ss2tf(At,Bt,Ct,Dt);%转换为传输函数表达式h,f=freqs(b,a);mag=abs(h);mag=20*log10(mag);phase=angle(h);subplot(2,1,1);plot(f/(2*pi),mag);title(Cheby1LowpassFilter);axis(01.5-350);xlabel(f(Hz);ylabel(幅度(dB);grid;subplot(2,1,2);plot(f/(2*pi),phase);axis(01.5-44);xlabel(f(Hz);ylabel(相位(ard);程序运行结果：（1）滤波器的阶数N=3（小于巴特沃斯的阶数）。（2）截止频率fn=3.1416ard/s=0.5Hz。（3）传递函数：b=00010.1356a=12.317910.088610.1356得到所设计的滤波器的传递函数为：（4）幅频特性曲线：通带内等波纹波动并且在截止频率0.5Hz处的衰减小于2dB，阻带截止频率1.2Hz处阻带最小衰减大于30dB，满足指标要求；相位特性曲线：在通带内相位特性接近线性。方法（2）：直接采用Matlab函数设计，代码设计部分改为：fp=0.5*2*pi;ap=2;fr=1.2*2*pi;as=30n,Wn=cheb1ord(fp,fr,ap,as,s)b,a=cheby1(n,fp,Wn,s)h,f=freqs(b,a);请读者在上述程序，加入绘制幅频特性、相频特性曲线的代码，调试、运行程序，并与前面结果比较，看是否相同。椭圆滤波器也称考尔滤波器（CauerFilter），其振幅平方函数为：式中：RN（，L）为N阶雅可比椭圆函数L是表示波纹性质的参量4.2.4椭圆滤波器（EllipticFilter）可见，在归一化通带内（-11），在（0，1）间振荡，而超过L后，在间振荡。这一特点使滤波器同时在通使滤波器同时在通带和阻带具有任意衰减量带和阻带具有任意衰减量。的特性曲线椭圆滤波器的振幅平方函数图中和A的定义同切比雪夫滤波器特点：幅值响应在通带和阻带内都是等波纹的幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和波纹要求，椭圆滤波器能获得较其它滤波器更窄的过渡带宽更窄的过渡带宽；就这点而言，椭圆滤波器是最优的。但是其设计计算最复杂，因此通常采用MATLAB设计，设计设计函数函数有：有：n,Wn=ellipord(Wp,Wr,Rp,Rr,s)z,p,k=ellipap(n,Rp,Rr)b,a=ellip(n,Rp,Rr,Wn,s)b,a=ellip(n,Rp,Rr,Wn,ftype,s)实际应用中：在给定相同指标下，三种模型设计都能满足指标要求，但椭圆滤波器阶次最低，切比雪夫次之，巴特沃兹最高，参数的灵敏度则恰恰相反。设计时需根据要求，合理选用。4.2.5 模拟高通、带通、带阻滤波器设计模拟高通、带通、带阻滤波器设计模拟高通、带通、带阻滤波器设计是通过设计模拟低通滤波器再通过频率变换实现的。假设：模拟低通滤波器的传输函数用表示，归一化频率用表示，令：，称为归一化拉斯复变量，为归一化低通滤波器；所要设计类型滤波器的传输函数用表示，归一化频率用表示，令：，称为归一化拉斯复变量，称为归一化传输函数。4.2.5.1低通到高通的频率变换低通到高通的频率变换公式：4.2.5.2.低通到带通的频率变换低通到带通的频率变换带通滤波器归一化边界频率用以下式子计算：，带通滤波器频率转换为归一化低通滤波器频率的关系式：或 归一化低通传输函数到带通传输函数的转换公式4.2.5.3 低通到带阻的频率变换低通到带阻的频率变换 或 4.2.5.4 由模拟低通原型设计模拟高通、带通、带阻由模拟低通原型设计模拟高通、带通、带阻滤波器的滤波器的Matlab函数函数Matlab提供的函数如下：1.低通原型转换为低通：NUMT,DENT=LP2LP(NUM,DEN,Wo)，Wo为截止频率(rad/sec)；2.低通原型转换为高通：NUMT,DENT=LP2HP(NUM,DEN,Wo)，Wo为截止频率(rad/sec)；3.低通原型转换为带通：NUMT,DENT=LP2BP(NUM,DEN,Wo,Bw)，Wo为中心频率(rad/sec)，Bw为通带带宽(rad/sec)；4.低通原型转换为带阻：NUMT,DENT=LP2BS(NUM,DEN,Wo,Bw)，Wo为中心频率(rad/sec)，Bw为阻带带宽(rad/sec)；4.2.5.5 模拟高通、带通、带阻滤波器直接设计的模拟高通、带通、带阻滤波器直接设计的Matlab函数函数Matlab提供了直接设计模拟滤波器函数：（1）b,a=butter(N,Wn,ftype,s)ftype选项确定滤波器类型；可选high、low、stop缺省时，设计的是低通、带通滤波器。当设计带通、带阻滤波器时，Wn为2元数组，结果为2*N阶。（2）切比雪夫型滤波器：b,a=cheby1(N,Rp,Wn,ftype,s)，Rp为通带最大波动（dB）。（3）切比雪夫型滤波器：b,a=cheby2(N,Rr,Wn,ftype,s)，Rr为阻带最小波动（dB）。（4）椭圆滤波器：b,a=ellip(N,Rp，Rr，Wn,ftype,s)，Rp为通带最大波动（dB），Rr为阻带最小波动（dB）。例4.2.43阶巴特沃斯型的归一化低通滤波器的传递函数为试设计下列模拟滤波器：（1）通带为10Hz的低通滤波器；（2）通带下边频为10Hz的高通滤波器；（3）中心频率为10Hz，带宽为2Hz的带通滤波器；（4）中心频率为10Hz，带宽为2Hz的带阻滤波器；又设采样周期为0.01秒，求相应的数字滤波器。解：Matlab程序如下：b=1;a=1221;Wn=10*2*pi;B=2*2*pi;b1,a1=lp2lp(b,a,Wn)；b2,a2=lp2hp(b,a,Wn)b3,a3=lp2bp(b,a,Wn,B)；b4,a4=lp2bs(b,a,Wn,B)运行结果：b1=2.4805e+005a1=1125.667895.72.4805e+005b2=11.2095e-014-4.2767e-013-2.7202e-020a2=1125.667895.72.4805e+005b3=1984.4-5.8081e-012-7.7026e-009-2.9869e-013a3=125.133121592.0042e+0054.8003e+0073.9171e+0086.1529e+010b4=13.5545e-013118442.8096e-0094.6756e+0075.5449e-0066.1529e+010a4=125.133121592.0042e+0054.8003e+0073.9171e+0086.1529e+0104.3 利用模拟滤波器设计利用模拟滤波器设计IIR数字滤波器数字滤波器利用模拟滤波器设计IIR数字滤波器，Ha(s)H（z）S平面Z平面这种映射变换应遵循两个基本原则：1）H（z）的频响要能模仿Ha(s)的频响，即S平面的虚轴应映射到Z平面的单位圆上。2）Ha(S)的因果稳定性映射成H（z）后保持不变，即S平面的左半平面ReS0应映射到Z平面的单位圆以内|Z|1。脉冲响应不变法是从滤波器的脉冲响应出发，使:数字滤波器的单位脉冲响应序列数字滤波器的单位脉冲响应序列h(n)h(n)正好等于模拟正好等于模拟滤波器的冲激响应滤波器的冲激响应h ha a(t)(t)的采样值的采样值，即:h(n)=ha(nT)T为采样周期 假设：Ha(s)及H（z）分别表示 ha(t)的拉氏变换及 h(n)的 Z变换，即:Ha(s)=Lha(t)H(z)=Zh(n)计算H(z)：4.3.1脉冲响应不变法脉冲响应不变法4.3.1.1变换关系变换关系假设模拟滤波器的传递函数若只只有有单单阶阶极极点点，且且分分母母的的阶数高于分子阶数阶数高于分子阶数 NM，则可表达为部分分式形式；其拉氏反变换为：对ha(t)采样得到数字滤波器的单位脉冲响应序列再对h(n)取Z变换，得到数字滤波器的传递函数：比较看到：（1）S平面上的极点S=Si变换到Z平面上是极点（2）Ha(s)与H(Z)中部分分式所对应的系数不变（3）这种Ha(s)到H(Z)的对应变换关系，只有将Ha(s)表达表达为部分分式形式才成立为部分分式形式才成立。4.3.1.2稳定性分析稳定性分析如果模拟滤波器是稳定的，则所有极点Si都在S左半平面，即Resi0,那么变换后H(z)的极点也都在单位圆以内，即：因此，数字滤波器保持稳定。4.3.1.3频率映射关系频率映射关系根据理想采样序列拉氏变换与模拟信号拉氏变换的关系。s平面与z平面的映射关系以上表明，采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它所完成的S平面到Z平面的变换，正是以前所讨论的拉氏变换到Z变换的标准变换关系，即：（1）首先对）首先对Ha(s)作周期延拓，作周期延拓，（2）然后再经过映射关系映射到）然后再经过映射关系映射到 Z 平面上。平面上。S平面上每一条横带，都将重叠地映射到Z平面的整个平面上:每一横带的左半部分映射到每一横带的左半部分映射到Z平面单位圆以内，平面单位圆以内，每一横带的右半部分映射到每一横带的右半部分映射到Z平面单位圆以外，平面单位圆以外，轴上每一段轴上每一段都对应于绕单位圆一周。都对应于绕单位圆一周。S平面Z平面4.3.1.4频率响应分析频率响应分析数字滤波器的频响并不是简单的重现模拟滤波器的频响，而是模拟滤波器频响的周期延拓：如果模拟滤波器的频响带限于折叠频率折叠频率S/2以内，这时数字滤波器的频响才能不失真地重现模拟滤波器的频响（存在于折叠频率S/2以内）但任何一个实际的模拟滤波器，其频响都不可能是真正带限的，因此不可避免地存在频谱的交叠，即混淆。脉冲响应不变法中的频响混淆【例4.3.1】将一个具有如下传递函数的模拟滤波器变换为数字滤波器。模拟滤波器的频率响应为:数字滤波器的频率响应为:4.3.1.5脉冲响应不变法的特点脉冲响应不变法的特点（1）频率变换是线性的，即：与是线性关系。因此如果模拟滤波器的频响带限于折叠频率以内的话，通过变换后数字滤波器的频响可不失真地反映原响应与频率的关系。（2）如果模拟滤波器是稳定的，映射后得到的数字滤波器也是稳定的。（3）缺点：有频谱周期延拓效应，会产生频率混叠。因此只能用于带限的频响特性，如衰减特性很好的低通或带通滤波器。至于高通和带阻滤波器，不适用。因此，脉冲响应不变法只有在一定要满足频率线性关系或严格要求保持瞬态响应时才采用。脉冲响应不变法的主要缺点是频谱交叠产生的混淆，这是从S平面到Z平面的标准变换z=esT的多值对应关系导致的，为了克服这一缺点，设想变换分为两步变换分为两步：第一步：将整个S平面压缩到S1平面的一条横带里；第二步：通过标准变换关系将此横带变换到整个Z平面上去。由此建立S平面与Z平面一一对应的单值关系，消除多值性，也就消除了混淆现象。4.3.2 双线性变换法双线性变换法4.3.2.1变换关系变换关系s平面s1平面z平面双线性变换法的映射关系为了将S平面的j轴压缩到S1平面j1轴上的/T到/T一段上，可通过以下的正切变换正切变换实现：取C=2/T，经过这样的频率变换，当由时,1由-/T经过变化到/T，即S平面的整个j轴被压缩到S1平面的。将这一关系解析扩展至整个S平面，则得到S平面到S1平面的映射关系：表明：（1）从s平面到z平面之间的映射关系都为一一对应的线性变换关系，因此该变换称为双线性变换。（2）当s左半平面映射在单位圆内，右半平面映射到单位圆外；因此，稳定的模拟滤波器通过双线性变换后，所得到的数字滤波器也是稳定的。S平面的虚轴(整个j)对应于Z平面单位圆的一周，S平面的=0处对应于Z平面的=0处，对应的数字滤波器的频率响应终止于折迭频率处，所以双线性变换不不存在混迭效应存在混迭效应。考虑z=ej ,4.3.2.2特点特点（1）s平面与z平面是单值的一一对应关系，即整个j轴单值的对应于单位圆一周，不存在频率混叠。（2）模拟域频率（）与数字域频率（）之间的变换成非线性关系：这将导致如下问题：a.数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变。例如：一个模拟微分器，它的幅度与频率是直线关系，但通过双线性变换后，就不可能得到数字微分器b.线性相位模拟滤波器经双线性变换后，得到的数字滤波器为非线性相位非线性相位。c.要求模拟滤波器的幅频响应必须是分段恒定分段恒定的，故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器。d.设计时需要预畸变即将模拟滤波器的临界频率事事先先加加以以畸畸变变，然后通过双线性变换后正好映射到所需要的频率上。双线性变换时频率的预畸变（3）双线性变换比脉冲响应法的设计计算更直接和简单从模拟滤波器的传递函数可直接通过代数置换得到数字滤波器的系统函数。变换关系式为【例4.3.2】设有一模拟滤波器试用双线性变换法将它转变为数字系统函数H(z)，设采样周期T=2。解：T=2时变换公式变为：4.3.3 用用Matlab实现模拟滤波器变换到实现模拟滤波器变换到IIR数字滤波器数字滤波器（1）脉冲响应不变法的Matlab变换函数语句格式为：BZ,AZ=impinvar(B,A,Fs)（2）双线性变换法的Matlab变换函数语句格式为：BZ,AZ=bilinear(B,A,Fs)其中：B为模拟滤波器系统函数的分子多项式系数矢量，A为模拟滤波器系统函数的分母多项式系数矢量，BZ为数字滤波器系统函数的分子多项式系数矢量，AZ为数字滤波器系统函数的分母多项式系数矢量，Fs为采样频率（单位为Hz）。【例4.3.3】用Matlab实现例4.3.1具有如下传递函数的模拟滤波器变换为数字滤波器，设采样频率1Hz。解：Matlab代码为：b=2;a=143;Fs=1;bz,az=impinvar(b,a,Fs);程序运行结果：bz=00.3181az=1.0000-0.41770.0183程序运行结果与例4.3.1的计算结果一致。【例4.3.4】用Matlab实现例4.3.2。解：Matlab代码：b=1;a=111;Fs=0.5;bz,az=bilinear(b,a,Fs);程序运行结果：bz=0.33330.66670.3333az=1.00000.00000.3333与例4.3.2的计算结果一致。4.4 从模拟低通原型到各种数字滤波器从模拟低通原型到各种数字滤波器的频率变换的频率变换原型变换映射变换原型变换模拟原型模拟低通、高通带通、带阻数字低通、高通带通、带阻也可把前两步合并成一步，直接从模拟低通归一化原型通过一定的频率变换关系，完成各类数字滤波器的设计。下面举例讨论设计各种数字滤波器的变换方法通过模拟低通原型设计数字低通滤波器的步骤：1）确定数字滤波器的性能指标临界频率k、通带阻带衰减。2）由变换关系将k映射到模拟域，得出模拟低通滤波器的临界频率值k。3）根据k设计模拟滤波器的Ha(s)4）把Ha(s)变换成H(z)（数字滤波器传递函数）变换方法有两种：脉冲响应不变法和双线性变换法，下面举例说明。4.4.1 低通变换低通变换【例4.4.1】设采样周期，已知数字低通滤波器通带截止频率，通带最大衰减2dB，阻带截止频率，阻带最小衰减20dB，试采用巴特沃斯模型设计该数字低通滤波器。解法一：脉冲响应不变法1）确定数字低通滤波器的技术指标。，通带最大衰减=2dB；，阻带最小衰减=20dB。2）由于脉冲响应不变法的频率变换关系是，因此有，所以相应的模拟低通滤波器的指标为：，通带最大衰减=2dB；，阻带最小衰减=20dB。因此，取N=3。查表得到模拟低通原型滤波器：3）设计巴特沃兹模拟低通滤波器去归一化得模拟低通滤波器：4）脉冲响应不变法变换可见，H（z）与采样周期T有关，T越小，H（z）的相对增益越大，这是不希望的。改进：H（z）乘以因子T，使H（z）只与有关。解法二：双线性变换法1）首先确定滤波器数字域指标，结果与脉冲响应不变法的步骤1）相同。2）确定模拟低通滤波器的技术指标。根据频率的非线性关系，进行频率预畸变。3)设计模拟低通滤波器因此，取N=3。3dB截止频率：4）进行双线性变换1、从模拟低通原型到数字低通滤波器的程序：fs=4000;Ts=1/fs;fp=500;ap=2;fr=1200;ar=20;%以下用脉冲响应不变法设计wp=2*pi*fp/fs;wr=2*pi*fr/fs;%wp、wr为数字频率Wp=wp/Ts;Wr=wr/Ts;%Wp、Wr为模拟频率N1,Wn=buttord(Wp,Wr,ap,ar,s)z,p,k=buttap(N1);%模拟低通原型设计b0=k*real(poly(z);a0=real(poly(p);B,A=lp2lp(b0,a0,Wn);%模拟低通原型变换到模拟低通滤波器num1,den1=impinvar(B,A,fs);%模拟滤波器经脉冲响应不变法变换到数字滤波器h1,w=freqz(num1,den1);解法三：用MATLAB完成设计%以下用双线性变换法设计Wp=2/Ts*tan(wp/2);Wr=2/Ts*tan(wr/2);%频率预畸变N2,Wn=buttord(Wp,Wr,ap,ar,s)；z,p,k=buttap(N2);b0=k*real(poly(z);a0=real(poly(p);B,A=lp2lp(b0,a0,Wn)num2,den2=bilinear(B,A,fs);%双线性变换h2,w=freqz(num2,den2);f=w/pi*fs/2;figure;plot(f,abs(h1),-.,f,abs(h2),-);text(600,0.66,leftarrow脉冲响应不变法);text(880,0.46,leftarrow双线性变换法);grid;xlabel(频率/Hz)；ylabel(幅度)2、从模拟滤波器到数字滤波器的设计程序：fs=4000;Ts=1/fs;fp=500;ap=2;fr=1200;ar=20;%以下用脉冲响应不变法设计wp=2*pi*fp/fs;wr=2*pi*fr/fs;%wp、wr为数字频率Wp=wp/Ts;Wr=wr/Ts;%Wp、Wr为模拟频率N1,fn=buttord(Wp,Wr,ap,ar,s)B,A=butter(N1,fn,s)；num1,den1=impinvar(B,A,fs);h1,w=freqz(num1,den1);%以下用双线性变换法设计Wp=2/Ts*tan(wp/2);Wr=2/Ts*tan(wr/2);N2,fn=buttord(Wp,Wr,ap,ar,s)B,A=butter(N2,fn,s)num2,den2=bilinear(B,A,fs);h2,w=freqz(num2,den2);f=w/pi*fs/2;figure;plot(f,abs(h1),-.,f,abs(h2),-);text(600,0.66,leftarrow脉冲响应不变法);text(880,0.46,leftarrow双线性变换法);grid;xlabel(频率/Hz)；ylabel(幅度)设计结果满足指标要求。设计高通数字滤波器时，只能采用双线性变换法。在模拟高通滤波器的设计中，低通至高通的变换是s变量的倒置，这一关系应用于双线性变换，则从模拟低通到数字高通的变换，只要将双线性变换式中的s代之以1/s，就可得到模拟低通到数字高通滤波器的变换关系式即：4.4.2 高通变换高通变换由于倒数关系不改变模拟滤波器的稳定性，因此，也不会影响双线变换后的稳定条件，而且轴仍映射在单位圆上，只是方向颠倒了。即如图1.01.00高通变换频率关系这一曲线的形状与双线性变换时的频率非线性关系曲线相对应，只是将坐标倒置，因而通过这一变换后可直接将模拟低通变为数字高通，如图。模拟低通原型到数字高通变换应当明确：所谓高通DF，并不是高到由于数字频域存在折叠频率，对于实数响应的数字滤波器因此有效的数字域仅是有效的数字域仅是 高通变换的计算步骤和低通变换一样。但在确定模拟原型预畸的临界频率时，应采用是镜象部分高通也仅指这一段的高端高通也仅指这一段的高端例例【例4.4.2】设计一数字高通滤波器，它的通带为0.61kHz(fs=2kHz)，通带允许波动0.5dB，阻带衰减在400Hz的频带内至少为20dB。解法一：利用切比雪夫滤波器模型的双线性变换设计1）确定数字高通滤波器的技术指标2）预畸变模拟边界频率确定模拟低通原型指标3）设计切比雪夫模拟低通滤波器查表去归一化可得模拟低通滤波器。用MATLAB编程：fs=2000;Ts=1/fs;fp=600;ap=0.5;fr=400;ar=20;wp=2*pi*fp*Ts;wr=2*pi*fr*Ts;%高通的数字频率Wp=2/Ts*cot(wp/2)%预畸变后的模拟低通频率Wr=2/Ts*cot(wr/2)N,Wn=cheb1ord(Wp,Wr,ap,ar,s);%阶数、截止频率z,p,k=cheb1ap(N,ap);%模拟低通原型设计b0=k*real(poly(z);a0=real(poly(p);B,A=lp2lp(b0,a0,Wn);%低通原型转换为模拟低通最小的滤波器阶数N=4。最后4条命令可用1条模拟滤波器设计命令代替：B,A=cheby1(N,ap,Wn,low,s)程序运行结果得：B=0，0，0，0，2.552591869406931013A=1，3.4798103，1.45003510107，2.516980323581010，2.703842446526551013即：模拟低通滤波器为：4）变换为数字滤波器解法二：直接使用MATLAB设计，代码如下：%数字高通滤波器fs=2000;Ts=1/fs;fp=600;ap=0.5;fr=400;ar=20;wp=2*pi*fp*Ts;wr=2*pi*fr*Ts;Wp=2/Ts*tan(wp/2);Wr=2/Ts*tan(wr/2);%预畸变N,wn=cheb1ord(Wp,Wr,ap,ar,s);B,A=cheby1(N,ap,wn,high,s);%模拟高通滤波器num,den=bilinear(B,A,fs);%变换到数字高通h,w=freqz(num,den);%频率响应f=w/(2*pi)*fs;plot(f,20*log10(abs(h);axis(0,fs/2,-80,10);grid;xlabel(频率/Hz)；ylabel(幅度/dB)；程序运行结果Num=0.03044487003703，-0.12177948014812，0.18266922022218，-0.12177948014812，0.03044487003703den=1，1.38342127667873，1.47205853973165，0.80124615697508，0.22859036591035图 带通原型变换4.4.3 带通变换带通变换如图,如果数字频域上带带通通的的中中心心频频率率为，则带通变换的目的是将:模拟低通将S的原点映射到满足这一要求的变换为：带通的频率变换关系图中点正好映射在上，而映射在，两端，因此满足带通变换的要求。带通变换的频率关系稳定性证明：由于上式完全是实数，所以是映射在S平面轴上。由此证明了，S左半平面映射在单位圆内，而右半平面左半平面映射在单位圆内，而右半平面映射在单位圆外映射在单位圆外，这种变换关系是稳定的变换关系，可用它来完成带通的变换。设由于带通滤波器的边界频率包括：上、下通带的截止频率；上、下阻带的截止频率；【例4.4.3】设采样频率fs=400kHz，试设计一个巴特沃兹带通数字滤波器，其技术指标满足：3dB上下通带截止频率分别110kHz、90kHz，上阻带截止频率为120kHz，阻带最小衰减10dB。解法一：1）确定带通数字滤波器的技术指标。2）求模拟低通的通带、阻带边界频率：3）设计模拟低通滤波器频率增加了约1.05倍，衰减增加了（10-3）dB，故选用二阶巴特沃兹滤波器可满足指标（查表）4）将模拟低通变换为数字带通滤波器解法二：MATLAB设计，代码如下：fs=400;f1=90;f2=110;fr=120;W1=2*fs*tan(2*pi*f1/(2*fs);%预畸变W2=2*f