资源描述：
《实验三 滤波器的结构new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、信息科学与工程学院数字信号处理实验报告实验三滤波器的结构一、实验目的1）认真复习数字滤波器的表示方法，信号流图的基本概念和画法；IIR数字滤波器的基本特性和FIR数字滤波器的基本特性。2）掌握IIR数字滤波器的直接I型、II型、级联型、并联型结构的特点和表示方法；FIR数字滤波器的直接型、级联型、频率抽样型以及快速卷积的结构特点和表示方法；线性相位，FIR滤波器的结构和特点；全通滤波器、梳状滤波器的特点和表示方法。二、实验内容实验代码：clc;n=0:100;a=[126-13];b=[12-43-5];h=impz(b,a)信息科学与工程学院数字信号处理实验报告s
2、tem(h,'.','g');title('冲击响应1')figuren=0:100;a=[2-14-3];b=[14-1-46];h1=impz(b,a)stem(h1,'.','g');title('冲击响应2')实验结果：信息科学与工程学院数字信号处理实验报告(1)直接I型直接2型信息科学与工程学院数字信号处理实验报告（2）直接1型直接2型信息科学与工程学院数字信号处理实验报告实验代码：主程序clc;b1=[16122-4-1];a1=[1-3112718];b2=[381272-2];a2=[1624241451];b3=[210233431164];a3=
3、[367887592671];[sos1,G1]=tf2sos(b1,a1);[sos2,G2]=tf2sos(b2,a2);[sos3,G3]=tf2sos(b3,a3);[C1,B1,A1]=tf2par(b1,a1);[C2,B2,A2]=tf2par(b2,a2);[C3,B3,A3]=tf2par(b3,a3);信息科学与工程学院数字信号处理实验报告tf2par函数：function[C,B,A]=tf2par(b,a)%直接型向并联型转换M=length(b);N=length(a);[r1,p1,C]=residuez(b,a);%先求系统的单根p1,
4、对应的留数r1及直接项Cp=cplxpair(p1,0.000000001);%先配对函数cplxpair由p1找共轭复根p，le-9的误差I=cplxcomp(p1,p);%找p1变为p的排序变化r=r1(I);%变为2阶子系统K=floor(N/2);B=zeros(K,2);A=zeros(K,3);ifK*2==N;fori=1:2:N-2pi=p(i:i+1,:);%取出一对极点ri=r(i:i+1,:);%取出一对对应留数[Bi,Ai]=residuez(ri,pi,[]);B(fix((i+1)/2),:)=real(Bi);A(fix((i+1)/2
5、),:)=real(Ai);end;[Bi,Ai]=residuez(r(N-1),p(N-1),[]);%处理单根B(K,:)=[real(Bi)0];A(K,:)=[real(Ai)0];elsefori=1:2:N-1信息科学与工程学院数字信号处理实验报告pi=p(i:i+1,:);ri=r(i:i+1,:);[Bi,Ai]=residuez(ri,pi,[]);B(fix((i+1)/2),:)=real(Bi);A(fix((i+1)/2),:)=real(Ai);end;end;cplxcomp函数：functionI=cplxcomp(p1,p2)I=
6、[];forj=1:length(p2)fori=1:length(p1)if(abs(p1(i)-p2(j))<0.0001)I=[I,i];end;end;end;I=I';实验结果：sos1=1.0000-3.00002.00001.0000-0.2500-0.12501.00000.00009.00001.00001.00000.5000信息科学与工程学院数字信号处理实验报告G1=0.0625sos2=1.0000-0.333301.00000.500001.00002.00002.00001.00000.50000.25001.00001.00001.00
7、001.00000.50000.5000G2=0.1875sos3=1.00003.00002.00001.00000.50000.25001.00001.00002.00001.00001.00000.33331.00001.00000.50001.00000.66670.3333G3=0.0556C1=-18B1=-10.0500-3.950028.1125-13.3625A1=1.00001.00000.50001.0000-0.2500-0.1250C2=-2B2=-0.8438-1.0313信息科学与工程学院数字信号处理实验报告0.68752.3438
