IIR数字滤波器要点(共15页).doc

《IIR数字滤波器要点(共15页).doc》由会员分享，可在线阅读，更多相关《IIR数字滤波器要点(共15页).doc（15页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上 IIR低通数字滤波器设计 院系：信院电子系 专业：电子信息工程 学号 姓名 专心-专注-专业目录 234摘要数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。数字滤波器技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统中，使用最多最复杂的要算数字

2、滤波器了。数字滤波器是数字滤波器除了具有数字信号处理的固有优点外，还有灵活性强、稳定性好、滤波精度高等优点。第1节 IIR数字滤波器的设计1.1 IIR数字滤波器的主要技术指标在设计IIR数字滤波器之前，首先了解数字滤波器的技术指标，才能将其转化成满足要求的模拟滤波器技术指标，再将模拟滤波器转化成数字滤波器。数字滤波器的主要技术指标为：(1) 特征频率参数 滤波器的频率参数主要有：通带截止频率：为通带与过渡带的边界点，在该点信号增益到规定的下限。阻带截频：为阻带于过渡带的边界点，在该点信号衰减到规定的下限。转折频率：为 信号功率衰减到(约为3)时的频率，但在多数情况下也常以作为通带或阻带截频。

3、当电路没有损耗时，固有频率，就是其谐振频率。(2) 增益和衰减 滤波器在通带内的增益并不是常数。对低通滤波器来说，通带增益一般指时的增益；对于高通而言，通带增益指时的增益；对于带通而言，通带增益则指中心频率处的增益。对带阻滤波器而言，应给出阻带衰减，衰减定义为增益的倒数。通带增益变化指通带内各点增益的最大变化量，如果以为单位，则指增益值的变化量。(3) 阻尼系数与品质因数 阻尼系数是表征滤波器对角频率为信号的阻尼作用，同时是滤波器中表示能量衰减的一项重要指标，也是与传递函数的极点实部大小相关的一项系数。阻尼系数的倒数1/称为品质因数Q，是评价带通和带阻滤波器频率选择特性的一个重要指标，它可以表

4、示为。其中为带通或带阻的带宽，为中心频率，在多数情况下中心频率等于固有频率。1.2 IIR数字滤波器的设计过程IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方法进行。其设计步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。这种方法比较容易一些，因为模拟滤波器的设计方法已经非常成熟，不仅有完整的设计公式，还有完善的图表供查阅；另外一种直接在频率或者时域内进行，由于需要解联立方程，设计时需要计算机做辅助设计。其设计步骤是：先设计过渡模拟滤波器得到系统函数，然后将按某种方法转换成数字滤波器的系统函数。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：(1)确定数字低通滤

5、波器的技术指标：通带边界频率、通带最大衰减、阻带截止频率、阻带最小衰减。(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。(3)按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。(4)用所选的转换方法，将模拟滤波器转换成数字低通滤波器系统函数。IIR数字滤波器的设计流程图如下：数字滤波器技术指标模拟滤波器 技术指标 数字滤波器 模拟滤波器 模拟滤波器设计方法变换变换图1.1 IIR数字滤波器的设计步骤流程图1.3双线性变换法设计IIR数字滤波器一 设计原理脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克服这一缺点，可

6、以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-/T/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，为了将S平面的整个虚轴j压缩到S1平面j1轴上的-/T到/T段上，可以通过以下的正切变换实现 式中,T仍是采样间隔。当1由-/T经过0变化到/T时，由-经过0变化到+，也即映射了整个j轴。将式（2-6）写成 将此关系解析延拓到整个S平面和S1平面，令j=s，j1=s1，

7、则得 再将S1平面通过以下标准变换关系映射到Z平面 从而得到S平面和Z平面的单值映射关系为： 这种变换都是两个线性函数之比，因此称为双线性变换式与式的双线性变换符合映射变换应满足的两点要求。首先,把z=ej，可得 即S平面的虚轴映射到Z平面的单位圆。其次，将s=+j代入式（2-12），得 因此 由此看出，当0时，|z|0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。第2节 IIR数字滤波器的实现2.1切比雪夫（Chebyshev）数字

8、滤波器的设计切比雪夫型滤波器在通带内幅度特性是等波纹的，在阻带内是单调的。切比雪夫滤波器在通带内是单调的，在阻带内是等波纹的。切比雪夫型低通滤波器的幅度平方函数为： 上式结果是一个正实数，式中的是一个小于1的正数，它与通带的波纹有关，越大，波纹越大，为通带的截止频率，时N阶切比雪夫多项式，定义为： 切比雪夫滤波器的传递函数参数确定：，N。当N大于或等于1时，切比雪夫多项式递推公式为：切比雪夫滤波器的阶数N等于通带内等幅波动的次数，即最大值和最小值的总数，可由下式确定， 其中均为与衰减有关的参数，Ap为通带内的最大衰减，As为阻带内的容许衰减，它们的单位为。因为滤波器的幅度平方函数为： 进而化简

9、得到：，又因为，由式(4-7)得到：故滤波器的阶数N为： 若要求阻带截止频率上的衰减越大，或过渡带内幅度特性越陡，则所需的阶数N越高。从式中可知切比雪夫滤波器的幅度平方函数也只有极点没有零点，且只需求出左半平面的极点，所以一旦求出，N，就能确定。利用MATLAB设计切比雪夫滤波器的实例。A、现在以MATLAB的函数cheblap(语法：z,p,k=cheblap(n，rp)，其中n为阶数，rp为通带的幅度差，z,p,k分别代表零点，极点，增益)为工具设计切比雪夫滤波器。设计技术指标为：通带截止频率，通带的最大衰减，阻带截止频率，阻带的最小衰减。程序如下：a=10000;b=40000;AP=3

10、;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); %求零、极点，增益syms rad; %定义基本符号对象hs1=K/(i*rad-P(1)/(i*rad-P(2)/(i*rad-P(3);hs2=10*log10(abs(hs1)2);ezplot(hs2,-15,15); %绘制符号表达式的二维曲线grid on; %画坐标网格线仿真图形如图1.2所示： 图1.2切比雪夫型带通滤波器响应曲线B、实际上利用MATLA

11、B的cheby1 函数设计出切比雪夫I型低通的模拟滤波器，然后采用变换的方法得到数字的高通、低通、带通和带阻滤波器。利用b,a=cheby1(n,Rp,Wn,ftype)方式设计出阶数为n、截止频率为Wn、带通波纹最大衰减为Rp 的数字低通滤波器。其中参数ftype 的形式决定了滤波器的形式，当它为high时得到滤波器为n 阶的、截止频率为Wn 的高通滤波器。若Wn是一个含有两个元素向量w1 w2，则cheby1 函数返回值是阶数为2n 的带通滤波器，滤波器的系统函数有理多项式的系数，通带范围是w1Ww2。利用 MATLAB 中的cheby2 设计各种形式的滤波器的方法和利用cheby1 设计

12、滤波器的方法相同，只是cheby2 设计出的滤波器在带阻是等波纹的，在带通是单调的；而cheby1 则恰好相反。对于采样频率为1000 Hz 的采样信号，设计一个阶数为9 阶、截止频率为300 Hz 的低通chebyshevII 数字滤波器，其中滤波器在带阻的波纹为20 db，其设计程序如下，频率响应如图1.3程序如下：b,a=cheby2(9,20,300/500);freqz(b,a,512,1000);仿真图形如图1.3所示： 图1.3 切比雪夫低通II 数字滤波器对于采样频率为1000 Hz 的采样信号，设计一个阶数为9 阶、截止频率为300 Hz 的低通cheby1 数字滤波器，其中

13、滤波器在带通的波纹为20db，其设计程序如下，频率响应如图1.3程序如下：b,a=cheby1(9,20,300/500);freqz(b,a,512,1000);仿真图形如图1.4所示： 图1.4切比雪夫低通I 数字滤波器从仿真图像可以看出，所设计的滤波器技术指标都基本达到要求，滤波器的设计是成功的。结 论通过对数字滤波器的设计，熟悉了MATLAB的运行环境，初步掌握了MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用；熟悉了滤波器设计的一般原理，对滤波器有了一个感性的认识；学会了数字滤波器设计的一般步骤；加深了对滤波器设计中产生误差的原因以及双线性变换法优缺点的理解和

14、认识。总之，使理论联系了实际，巩固并深化了对课本基本知识的认识和理解，使理论得以升华。致 谢本论文是在赵春雨老师的支持悉心指导下完成的。她严肃的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和激励着我。从课题的选择到项目的最终完成，令老师都始终给予我细心的指导和不懈的支持。在此谨向赵老师致以诚挚的谢意和崇高的敬意。我还要感谢在一起讨论的室友，正是由于你们的帮助和支持，我才能克服一个一个的困难和疑惑，直至本文的顺利完成。在论文即将完成之际，我的心情无法平静，从开始进入课题到论文的顺利完成，各位可敬的师长、同学、朋友给了我太多精神的鼓舞和学术上的帮助，再次表示衷心感谢!参考文献1 倪养华

15、数字信号处理与实现上海：上海交通大学出版社，19982 飞思科技产品研发中心MATLAB7辅助信号处理技术与应用北京：电子工业出版社， 20053 楼顺天基于MATLAB的系统分析与设计一信号处理西安：西安电子科技大学出版社，1998.4 张葛祥MATLAB仿真技术与应用北京：清华大学出版社，20035 胡广书 数字信号处理、理论、算法与实现M 北京：清华大学出版社，19976 陈希林一种LabWindows/CVI与MATLAB混合编程的实现方法J.微计算机信息，20057 刘波MATLAB信号处理北京：电子工业出版社，20068 施阳MATLAB语言工具箱西安：西北工业大学出版社，19999 高西全数字信号处理西安电子科技大学出版社 第三版，2008