IIR数字滤波器要点

IIIR低通数字滤波器设计院系：信院电子系专业：电子信息工程学号姓名I目录摘要...........................................................2第1节IIR数字滤波器的设计.....................................31.1IIR数字滤波器的主要技术指标...........................41.2IIR数字滤波器的设计过程...............................41.3双线性变换法设计IIR数字滤波器..........................5第2节IIR数字滤波器的实现.....................................62.1切比雪夫数字滤波器的设计................................6结论..........................................................12致谢..........................................................13参考文献.......................................................142摘要数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。数字滤波器技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统中，使用最多最复杂的要算数字滤波器了。数字滤波器是数字滤波器除了具有数字信号处理的固有优点外，还有灵活性强、稳定性好、滤波精度高等优点。3第1节IIR数字滤波器的设计1.1IIR数字滤波器的主要技术指标在设计IIR数字滤波器之前，首先了解数字滤波器的技术指标，才能将其转化成满足要求的模拟滤波器技术指标，再将模拟滤波器转化成数字滤波器。数字滤波器的主要技术指标为：(1)特征频率参数滤波器的频率参数主要有：①通带截止频率/2ppf：为通带与过渡带的边界点，在该点信号增益到规定的下限。②阻带截频/2rrf：为阻带于过渡带的边界点，在该点信号衰减到规定的下限。③转折频率/2ccf：为信号功率衰减到1/2(约为3dB)时的频率，但在多数情况下也常以cf作为通带或阻带截频。④当电路没有损耗时，固有频率/2f，就是其谐振频率。(2)增益和衰减滤波器在通带内的增益并不是常数。①对低通滤波器来说，通带增益pK一般指0时的增益；对于高通而言，通带增益pK指时的增益；对于带通而言，通带增益pK则指中心频率处的增益。②对带阻滤波器而言，应给出阻带衰减，衰减定义为增益的倒数。③通带增益变化pK指通带内各点增益的最大变化量，如果pK以dB为单位，则指增益dB值的变化量。(3)阻尼系数与品质因数阻尼系数是表征滤波器对角频率为信号的阻尼作用，同时是滤波器中表示能量衰减的一项重要指标，也是与传递函数的极点实部大小相关的一项系数。阻尼系数的倒数1/称为品质因数Q，是评价带通和带阻滤波器频率选择特性的一个重要指标，它可以表示为Q。其中为带通或带阻的3dB带宽，为中心频率，在多数情况下中心频率等于固有频率。1.2IIR数字滤波器的设计过程IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方法4进行。其设计步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。这种方法比较容易一些，因为模拟滤波器的设计方法已经非常成熟，不仅有完整的设计公式，还有完善的图表供查阅；另外一种直接在频率或者时域内进行，由于需要解联立方程，设计时需要计算机做辅助设计。其设计步骤是：先设计过渡模拟滤波器得到系统函数)(sHa，然后将)(sHa按某种方法转换成数字滤波器的系统函数)(zH。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：(1)确定数字低通滤波器的技术指标：通带边界频率p、通带最大衰减p、阻带截止频率s、阻带最小衰减s。(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。(3)按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。(4)用所选的转换方法，将模拟滤波器)(sHa转换成数字低通滤波器系统函数)(zH。IIR数字滤波器的设计流程图如下：图1.1IIR数字滤波器的设计步骤流程图)(ZfS数字滤波器技术指标模拟滤波器技术指标数字滤波器)(ZH模拟滤波器)(SHa模拟滤波器设计方法变换)(g变换51.3双线性变换法设计IIR数字滤波器一设计原理脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T～π/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现2tan21TT式中,T仍是采样间隔。当Ω1由-π/T经过0变化到π/T时，Ω由-∞经过0变化到+∞，也即映射了整个jΩ轴。将式（2-6）写成2/2/2/2/111.2TjTjTjTjeeeeTj将此关系解析延拓到整个S平面和S1平面，令jΩ=s，jΩ1=s1，则得TsTsTsTsTsTsiiiiieeTTsTeeeeTs11.22tanh2.2112/2/2/2/再将S1平面通过以下标准变换关系映射到Z平面Tsez1从而得到S平面和Z平面的单值映射关系为：11112zzTssTsTsTsTz222121这种变换都是两个线性函数之比，因此称为双线性变换6式与式的双线性变换符合映射变换应满足的两点要求。首先,把z=ejω，可得jTjeeTsjj2tan2112即S平面的虚轴映射到Z平面的单位圆。其次，将s=σ+jΩ代入式（2-12），得jTjTz22因此222222TTz由此看出，当σ0时，|z|1；当σ0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。第2节IIR数字滤波器的实现2.1切比雪夫（Chebyshev）数字滤波器的设计切比雪夫Ⅰ型滤波器在通带内幅度特性是等波纹的，在阻带内是单调的。切比雪夫Ⅱ滤波器在通带内是单调的，在阻带内是等波纹的。切比雪夫Ⅰ型低通滤波器的幅度平方函数为：7)(11)(222cNajjCjH上式结果是一个正实数，式中的是一个小于1的正数，它与通带的波纹有关，越大，波纹越大，c为通带的截止频率，()NCx时N阶切比雪夫多项式，定义为：1))(coscos(1))((11)(xxNxxNchchNxC切比雪夫滤波器的传递函数参数确定：，c，N。当N大于或等于1时，切比雪夫多项式递推公式为：)()()(11xCxNCxCNNN切比雪夫滤波器的阶数N等于通带内等幅波动的次数，即最大值和最小值的总数，可由下式确定0.10.5(101)pA，0.10.5(101)sA其中,均为与衰减有关的参数，Ap为通带内的最大衰减，As为阻带内的容许衰减，它们的单位为dB。因为滤波器的幅度平方函数为：222211)(11)(cNajjCjH进而化简得到：()NcjCj，又因为(/)1c，由式(4-7)得到：))/((/j1ccNjjNchchjC）（8故滤波器的阶数N为：)()(11echchN若要求阻带截止频率上的衰减越大，或过渡带内幅度特性越陡，则所需的阶数N越高。从式中可知切比雪夫滤波器的幅度平方函数也只有极点没有零点，且只需求出左半平面的极点，所以一旦求出，c，N，就能确定()aHs。利用MATLAB设计切比雪夫滤波器的实例。A、现在以MATLAB的函数cheblap(语法：［z,p,k］=cheblap(n，rp)，其中n为阶数，rp为通带的幅度差，z,p,k分别代表零点，极点，增益)为工具设计切比雪夫Ⅰ滤波器。设计技术指标为：通带截止频率410/crads，通带的最大衰减3pAdB，阻带截止频率4410/srads，阻带的最小衰减35sAdB。程序如下：a=10000;b=40000;AP=3;AS=35;%技术指标T1=sqrt(10^(0.1*AP)-1);T2=sqrt(10^(0.1*AS)-1);N=ceil(acosh(T2/T1)/acosh(b/a));%求阶数[Z,P,K]=cheb1ap(N,AP);%求零、极点，增益symsrad;%定义基本符号对象hs1=K/(i*rad-P(1))/(i*rad-P(2))/(i*rad-P(3));hs2=10*log10((abs(hs1))^2);ezplot(hs2,[-15,15]);%绘制符号表达式的二维曲线gridon;%画坐标网格线仿真图形如图1.2所示：9-15-10-5051015-80-70-60-50-40-30-20-100rad(10log(20378958691192777514264974195041/(324518553658426726783156020576256abs(radi+2689731717386215/18014398509481984-(1017604581116467i)/1125899906842624)2abs(radi+2689731717386215/18014398509481984+(1017604581116467i)/1125899906842624)2abs(2689731717386217/9007199254740992+radi)2)))/log(10)图1.2切比雪夫Ⅰ型带通滤波器响应曲线B、实际上利用MATLAB的cheby1函数设计出切比雪夫I型低通的模拟滤波器，然后采用变换的方法得到数字的高通、低通、带通和带阻滤波器。利用[b,a]=cheby1(n,Rp,Wn,‘ftype’)方式设计出阶数为n、截止频率为Wn、带通波纹最大衰减为Rp的数字低通滤波器。其中参数ftype的形式决定了滤波器的形式，当它为‘high’时得到滤波器为n阶的、截止频率为Wn的高通滤波器。若Wn是一个含有两个元素向量[w1w2]，则cheby1函数返回值是阶数为2n的带通滤波器，滤波器的系统函数有理多项式的系数，通带范围是w1Ww2。利用MATLAB中的cheby2设计各种形式的滤波器的方法和利用cheby1设计滤波器的方法相同，只是cheby2设计出的滤波器在带阻是等波纹的，在带通是单调的；而cheby1则恰好相反。对于采样频率为1000Hz的采样信号，设计一个阶数为9阶、截止频率为300Hz的低通chebyshevII数字滤波器，其中滤波器在带阻的波纹为20db，其设计程序如下，频率响应如图1.3程序如下：[b,a]=cheby2(9,20,300/500);freqz(