实验三 滤波器的结构new.doc

实验三 滤波器的结构new.doc信息科学与工程学院数字信号处理实验报告实验三滤波器的结构一、实验目的 1)认真复****数字滤波器的表示方法,信号流图的基本概念和画法; IIR 数字滤波器的基本特性和 FIR 数字滤波器的基本特性。 2)掌握 IIR 数字滤波器的直接 I型、 II型、级联型、并联型结构的特点和表示方法; FIR 数字滤波器的直接型、级联型、频率抽样型以及快速卷积的结构特点和表示方法;线性相位, FIR 滤波器的结构和特点;全通滤波器、梳状滤波器的特点和表示方法。二、实验内容实验代码: clc; n=0:100; a=[1 26 -1 3]; b=[1 2 -43 -5]; h=impz(b,a) 信息科学与工程学院数字信号处理实验报告 stem(h ,'.','g'); title(' 冲击响应 1') figure n=0:100; a=[2 -14 -3]; b=[1 4 -1 -4 6]; h1=impz(b,a) stem(h1,'.','g'); title(' 冲击响应 2')实验结果: 信息科学与工程学院数字信号处理实验报告(1)直接 I型直接 2型信息科学与工程学院数字信号处理实验报告(2)直接 1型直接 2型信息科学与工程学院数字信号处理实验报告实验代码: 主程序 clc; b1=[16 122 -4 -1]; a1=[1 -3 11 27 18]; b2=[3 8 1272 -2]; a2=[16 24 24 145 1]; b3=[2 10 23 34 31 16 4]; a3=[36 78 87 59 267 1]; [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, 对应的留数 r1及直接项 C p=cplxpair(p1,);% 先配对函数 cplxpair 由 p1找共轭复根 p, le-9 的误差 I=p(p1,p);% 找 p1变为 p的排序变化 r=r1(I); %变为 2阶子系统 K=floor(N/2);B=zeros(K,2);A=zeros(K,3); if K*2==N; for i=1:2:N-2 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; [Bi,Ai]=residuez(r(N-1),p(N-1),[]);% 处理单根 B(K,:)=[real(Bi) 0];A(K,:)=[real(Ai) 0]; else for i=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; p 函数: function I=p(p1,p2) I=[]; for j=1:length(p2) for i=1:length(p1) if (abs(p1(i)-p2(j))<) I=[I,i]; end; end; end; I=I'; 实验结果: sos1 = - - - 信息科学与工程学院数字信号处理实验报告 G1 = sos2 = - 0 0 G2 = sos3 = 1