资源描述：
《实验七滤波器设计》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验七滤波器设计一．实验目的：掌握滤波器结构类型转换的Matlab实现方法；掌握用滤波器直接形式、级联形式、并联形式、格型形式以及频率采样形式的Matlab实现。了解和学习组合滤波的结构类型特性；进一步理解和巩固理论知识，提高分析和解决实际问题的能力。二．实验原理：1．滤波器数学模型之间的转换2．滤波器结构类型之间的转换三．实验内容：（一）已知：IIR系统的传递函数为：，求：1.直接I型和直接II型；2.包括二阶直接II型子系统的级联形式；3．包括二阶直接II型子系统的并联形式；4．格型梯形形式。解：（1）先定义cas2dir.m文件clearall;
2、clc;b0=2.45;B=[1,-1.3,01,-3.8,01,0.3,0.15];A=[1,0.83,-0.151,1.26,0.51,-0.625,0];[b,a]=cas2dir(b0,B,A)b=2.4500-11.76008.72201.75661.815400a=1.00001.46500.0896-0.6464-0.21630.04690（2）先定义dir2cas.m文件clearall;clc;b=[2.45,-11.76,8.722,1.7566,1.8154];a=[1,1.465,0.0896,-0.6464,-0.2163,0
3、.0469];[b0,B,A]=dir2cas(b,a)b0=2.4500B=1.00000.30000.15001.0000001.0000-5.10004.9400A=1.00001.26000.50011.00000.982701.0000-0.77770.0954（3）先定义dir2par.m文件与cplxcomp.m文件clearall;clc;b=[2.45,-11.76,8.722,1.7566,1.8154];a=[1,1.465,0.0896,-0.6464,-0.2163,0.0469];[C,B,A]=dir2par(b,a)r=
4、-23.7795+8.7951i-23.7795-8.7951i52.8612-5.85233.0000C=[]B=-47.5589-24.312647.0089-13.82233.00000A=1.00001.26000.50011.00000.8300-0.15001.0000-0.62500（4）先定义dir2ladr.m文件clearall;clc;b=[2.45,-11.76,8.722,1.7566,1.8154];a=[1,1.465,0.0896,-0.6464,-0.2163,0.0469];[K,C]=dir2ladr(b,a)K=
5、00000.0469C=2.4500-11.76008.72201.75661.81540（二）已知线性系统由下列系统函数描述：求：1。线性相位形式；2．频率采样形式。解：（1）clearall;clc;b=[1/32,1/16,1/8,1/4,1/2,1,1/2,1/4,1/8,1/16,1/32];broots=roots(b)B1=real(poly([broots(1),broots(2),broots(5),broots(6)]))B2=real(poly([broots(3),broots(4),broots(7),broots(8)]))
6、broots=1.0403+1.6421i1.0403-1.6421i-0.8748+1.6437i-0.8748-1.6437i-1.83090.2753+0.4346i0.2753-0.4346i-0.5462-0.2523+0.4741i-0.2523-0.4741iB1=1.0000-0.52510.03836.9269-1.9046B2=1.00002.02043.79060.6761-0.5213（2）先定义dir2fs.m文件clearall;clc;h=[1/32,1/16,1/8,1/4,1/2,1,1/2,1/4,1/8,1/16,
7、1/32];[C,B,A]=dir2fs(h)C=3.74311.76581.09200.78060.68092.9375B=00.95950-0.841300.65490-0.415400.1423A=1.0000-1.68251.00001.0000-0.83081.00001.00000.28461.00001.00001.30971.00001.00001.91901.00001.0000-1.00000（三）一个FIR滤波器的格型参数为：求出它的脉冲响应，并求出其直接形式和格型形式结构。解：先定义latcfilt.m文件与impseq.m文件
8、clearall;clc;[x,n]=impseq(0,0,3);formatlongK=[2
