数字信号及图像处理实验7 update.pptxVIP

实验七滤波器设计一、实验内容1：已知：IIR系统的传递函数为：求：1.直接型；（提示：题目给出的是级联形式，调用cas2dir函数，P172）2.包括二阶直接II型子系统的级联形式；（提示：将第1小题求到的a和b系数带入dir2cas函数，P170）3．包括二阶直接II型子系统的并联形式；（提示：将第1小题求到的a和b系数带入dir2par函数，P174）（提示：直接型系数a,b的最后一个零要去掉）

相关理论：IIR滤波器基本结构（第六章P167-180页）直接型：系数组：*4．格型梯形形式。（提示：将第1小题求到的a和b系数带入dir2ladr函数，P198）注意：dir2par程序中：B(K,:)=[real(Brow’)0];A(K,:)=[real(Arow’)0];应改成：B(K,:)=[real(Brow)0];A(K,:)=[real(Arow)0];或者B(K,:)=[real(Brow’);0];A(K,:)=[real(Arow’);0]176页:I=[I,i];

并联型：系数组：系数组之间的转换。级联型：系数组：

1.1直接型转级联型matlab算法实现说明：（P170页）1.计算增益系数b0：b0=b(1)/a(1)，同时bn=bn/b(1)、an=an/a(1)；2.用补零方式使bn和an等长，保证每一个双二阶分式的分子和分母均不为零；3.求bn和an多项式的根：broots=roots(bn)；%broots为K*1的列数组；aroots=roots(an)；%aroots为K*1的列数组；4.将broots和aroots按共轭复数对的次序排列：broots=cplxpair(broots);aroots=cplxpair(aroots);5.将broots和aroots按共轭复数对展开成二阶多项式的形式：fori=1:2:2*KBrow=broots(i:1:i+1,:);%提取共轭复数对Brow=real(poly(Brow));%由共轭复数对求二阶多项式的系数,1*3数组B(fix((i+1)/2),:)=Brow;%将每一组二阶多项式系数组合成系数矩阵BArow=aroots(i:1:i+1,:);%提取共轭复数对Arow=real(poly(Arow));%由共轭复数对求二阶多项式的系数,1*3数组A(fix((i+1)/2),:)=Arow;%将每一组二阶多项式系数组合成系数矩阵Aend函数[b0,B,A]=dir2cas(b,a)，直接型转级联型，an、bn转b0、A、B：*注意：P170页[b0,B,A]=dir2cas(b,a)函数前面的注释部分的说明。

1.2级联型转直接型matlab算法实现说明：（系数组之间的转换，P172页）函数[b,a]=cas2dir(b0,A,B)，直接型转级联型，b0、A、B转an、bn：1.A和B系数组逐级展开：Fori=1:1:Kb=conv(b,B(i,:));%由B的每个二阶多项式系数展开成bn系数；a=conv(a,A(i,:));%由A的每个二阶多项式系数展开成an系数；end2.b系数组乘以增益系数b0：b=b*b0;

1.3直接型转并联型matlab算法实现说明：（P174-175页）1.用留数法将an、bn的多项式转化为留数、极点的分式项，并获得多项式余部C；[r1,p1,C]=residuez(b,a);%helpresiduez显示p1和r1构分式项2.将极点p1按照复共轭对及实数的顺序排列：p=cplxpair(p1eps);3.将留数r1按照p1对应的复共轭、实数顺序排列：p(p1,p)；p将两个混乱的复数数组进行比较，P175页；r=r1(I)；4.将上述得到的留数、极点分式项，按照复共轭对或实数对的方式，两两合并成：fori=1:2:N-2Rrow=r(i:1:i+1,:);%