数字的信号处理实验作业.doc
word实验6 数字滤波器的网络结构一、实验目的：1、加深对数字滤波器分类与结构的了解。2、明确数字滤波器的根本结构与其相互间的转换方法。3、掌握用MATLAB语言进展数字滤波器结构间相互转换的子函数与程序编写方法。二、实验原理：1、数字滤波器的分类离散LSI系统对信号的响应过程实际上就是对信号进展滤波的过程。因此，离散LSI系统又称为数字滤波器。数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以与全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器FIR和无限长单位脉冲响应滤波器IIR。一个离散LSI系统可以用系统函数来表示：也可以用差分方程来表示：以上两个公式中，当ak至少有一个不为0时，如此在有限Z平面上存在极点，表达的是以一个IIR数字滤波器；当ak全都为0时，系统不存在极点，表达的是一个FIR数字滤波器。FIR数字滤波器可以看成是IIR数字滤波器的ak全都为0时的一个特例。IIR数字滤波器的根本结构分为直接型、直接型、直接型、级联型和并联型。FIR数字滤波器的根本结构分为横截型又称直接型或卷积型、级联型、线性相位型与频率采样型等。本实验对线性相位型与频率采样型不做讨论，见实验10、12。另外，滤波器的一种新型结构格型结构也逐步投入应用，有全零点FIR系统格型结构、全极点IIR系统格型结构以与全零极点IIR系统格型结构。2、IIR数字滤波器的根本结构与实现1直接型与级联型、并联型的转换例6-1一个系统的传递函数为将其从直接型其信号流图如图6-1所示转换为级联型和并联型。图6-1分析：从直接型转换为级联型，就是将系统的传递函数tf模型转换为二次分式sos模型；从直接型转换为并联型，就是将系统函数的传递函数tf模型转换为极点留数rpk模型。程序清单如下：b=8,-4,11,-2;a=1,-1.25,0.75,-0.125;sos,g=tf2sos(b,a)r,p,k=residuez(b,a)运行结果如下：sos = 1.0000 -0.1900 0 1.0000 -0.2500 0g = 8r = 8.0000 p = 0.2500 k =16由sos和g的数据，可以列写出级联型的表达式：信号流图如图6-2所示：图6-2由r、p、k的数据，可以列写出并联型的表达式：上式中出现了复系数，可采用二阶分割将共轭极点组成分母上的实系数二阶环节。这里使用自定义函数dir2par可以实现滤波器结构从直接型向并联型的转换，且用实系数二阶环节表示。在使用该函数时，调用了另一个自定义函数cplxp以进展复共轭对的正确排序，保证系统二阶环节的分子、分母一定是实数。dir2par函数和cplxp函数定义如下：function I=cplxp(p1,p2) %按共轭条件排列极点留数对 %比拟两个包含同样标量元素但可能具有不同下标的复数对 %本语句必须用在p2=cplxpair(p1)语句之后，以重新排序对应的留数向量I=;for j=1:length(p2) for i=1:length(p1) if(abs(p1(i)-p2(j)<0.0001) I=I,i; end endendI=I'function C,B,A=dir2par(num,den) %直接型到并联型的转换M=length(num);N=length(den);r1,p1,C=residuez(num,den);%先求系统的单根p1对应的留数r1与直接项Cp=cplxpair(p1,10000000*eps);I=cplxp(p1,p);r=r1(I);K=floor(N/2);B=zeros(K,2);A=zeros(K,3);if K*2=N; for i=1:2:N-2; Brow=r(i:1:i+1,:); Arow=p(i:1:i+1,:); Brow,Arow=residuez(Brow,Arow,); B(fix(i+1)/2),:)=real(Brow); A(fix(i+1)/2),:)=real(Arow); endBrow,Arow=residuez(r(N-1),p(N-1),);B(K,:)=real(Brow),0;A(K,:)=real(Arow),0;else for i=1:2:N-1; Brow=r(i:1:i+1,:); Arow=p(i:1:i+1,:); Brow,Arow=residuez(Brow,Arow,); B(fix(i+1)/2),:)=real(Brow); A(fix(i+1)/2),:)=real(Arow); end end将例6-1从直接型转换为并联型的程序改写如下：b=8,-4,11,-2;a=1,-1.25,0.75,-0.125;C,B,A=dir2par(b,a)运行结果如下：C =16B = 8.0000 0A =1.0000 -0.2500 0由A,B,C的数据可以直接写出并联型的表达式：信号流图如图6-3所示：图6-3例6-2一个系统的级联型系数公式为将其从级联型信号流图如图6-4所示转换为直接型和并联型结构。图6-4分析：从级联型转换为直接型，就是将二次分式sos模型转换为系统传递函数tf模型；再使用dirpar.m和cplxp.m函数将直接型转换为并联型。程序清单如下：sos=1 0.9 0 1 -0.25 0 1 -3 2 1 1 0.5; g=0.5; b,a=sos2tf(sos,g) C,B,A=dir2par(b,a)程序运行结果如下： 3.7154 0 1.0000 -0.2500 0由b,a的数据可以直接写出直接型的表达式：信号流图如图6-5所示：图6-5由A,B,C的数据可以写出并联型的表达式：信号流图如图6-6所示：图6-62直接型转换为全零极点IIR系统的格型结构例6-3将例6-1给定的系统传递函数从直接型转换为格型。程序清单如下：b=8,-4,11,-2;a=1,-1.25,0.75,-0.125;K,C=tf2latc(b,a)b,a=latc2tf(K,C)程序运行结果如下：K =C =由K、C参数可以画出格型结构图，如同图6-7所示：图6-73直接型转换为全极点IIR系统的格型结构例6-4将一个全极点IIR系统的传递函数从直接型转换为格型结构。程序清单如下：b=1;a=1,-1.25,0.75,-0.125;K=tf2latc(b,a)b,a=latc2tf(K)程序运行结果如下：a = 1格型结构信号流图如图6-8所示：图6-83、FIR数字滤波器的根本结构与实现1横截型与级联型之间的转换例6-5一个FIR系统的传递函数为将其从横截型信号流图如图6-9所示转换为级联型。图6-9分析：从横截型转换为级联型就是将系统传递函数tf模型转换为二次分式sos模型。程序清单如下：b=2,0.9,1.55,2.375;a=1;sos,g=tf2sos(b,a)b,a=sos2tf(sos,g)程序运行结果如下：sos = 1.0000 0.9500 0 1.0000 0 0 1.0000 -0.5000 1.2500 1.0000 0 0g = 2a = 1 0 0 0由sos和g的数据可以写出级联型的表达式：信号流图如图6-10所示：图6-102横截型转换为全零点FIR系统的格型结构例6-6一个FIR系统的传递函数为将其从横截型转换为全零点FIR系统的格型结构程序清单如下：b=1,2.7917,2,1.375,0.3333;a=1;K=tf2latc(b,a)b,a=latc2tf(K)程序运行结果如下：K =a = 1由K参数画出的格型结构图如图6-11所示：图6-11三、实验内容：1、阅读并输入实验原理中介绍的例题程序，观察输出的数据，结合根本原理理解每一条语句的含义。2、一个IIR系统的传递函数为将其从直接型转换为级联型、并联型和格型结构，并画出各种结构的信号流图。答：直接型：Z-1Z-1Z-11X(n)y(n)1、级联型b=0.1,-0.4,0.4,-0.1;a=1,0.3,0.55,0.2;sos,g=tf2sos(b,a)r,p,k=residuez(b,a)sos = 1.0000 -2.6180 0 1.0000 0.3519 0g =r = 1.1786 p = -0.3519 k =表达式：-1/(1+0.3519 Z-1-1-2)/1-0.0519 Z-1+-2y(n)Z-1Z-1Z-12、并联型：b=0.1,-0.4,0.4,-0.1;a=1,0.3,0.55,0.2;C,B,A=dir2par(b,a)C =B = 1.1786 0A =1.0000 0.3519 0表达式：-1-1-2-1)3、格型：b=0.1,-0.4,0.4,-0.1;a=1,0.3,0.55,0.2;K,C=tf2latc(b,a)b,a=latc2tf(K,C)K =C =b =a =X(n) y(n)3、一个FIR系统的传递函数为将其从横截型转换为级联型和格型结构，并画出各种结构的信号流图。1、级联型：b=0.2,0.88,0.212,0.212,0.885;a=1;sos,g=tf2sos(b,a)b,a=sos2tf(sos,g)sos = 1.0000 5.2595 4.6234 1.0000 0 0 1.0000 -0.8595 0.9571 1.0000 0 0g =b =a = 1 0 0 0表达式：-1-2-1-2Z-1Z-1Z-1Y(n)Xn2、格型：K =b =a = 1YnX(n)四、实验预习：1、认真阅读实验原理局部，明确实验目的，读懂例题程序，了解实验方法。2、根据实验内容预先编写实验程序。3、预习思考题：什么是数字滤波器？数字滤波器是如何分类的？离散LSI系统对信号的响应过程实际上就是对信号进展滤波的过程。因此，离散LSI系统又称为数字滤波器。数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以与全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器FIR和无限长单位脉冲响应滤波器IIR归纳各类数字滤波器的根本结构。IIR数字滤波器的根本结构分为直接型、直接型、直接型、级联型和并联型。FIR数字滤波器的根本结构分为横截型又称直接型或卷积型、级联型、线性相位型与频率采样型等五、实验报告：1、列写调试通过的实验程序，并给出程序运行结果。2、给出预习思考题答案。15 / 15
