实验四数字滤波器的结构教案解析.ppt

一、实验目的 加深对数字滤波器分类与结构的了解。 掌握数字滤波器的基本结构及其相互间的转换方法。 学习利用MATLAB语言进行数字滤波器各种结构相互间转换。 二、实验原理及方法 数字滤波器的分类 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据系统的单位冲激响应的特性，又可以分为有限长（FIR）和无限长（IIR）冲激响应滤波器。 一个离散LSI系统可以用系统函数表示： 也可用差分方程表示： 数字滤波器的基本结构 IIR数字滤波器的基本结构分为直接I型、直接II型、级联型和并联型。 FIR数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型。 三、实验用函数 1、tf2sos 功能：将系统传递函数模型(tf)转换为系统函数的二次分式(sos)模型。 调用格式： [sos,g]=tf2sos(num,den)：输入参数为传递函数分子、分母的系数向量，求系统函数二次分式模型的系统矩阵sos、增益系数g。 2、sos2tf 功能：将系统函数的二次分式(sos)模型转换为系统传递函数模型(tf)。 调用格式： [num,den]=tf2sos(sos,g)：输入系统函数二次分式模型的系统矩阵sos、增益系数g，求传递函数分子、分母的系数向量。 四、基本结构的转换 IIR数字滤波器的基本结构与实现 直接型与级联型、并联型之间的转换 已知一个系统的传递函数为： 将其由直接型转换为级联型、并联型。信号流图如下图所示。 从直接型到级联型就是将系统传递函数模型转换为二次分式模型；直接型转换为并联型，就是将系统传递函数转换为极点留数模型。 程序运行结果： sos =1.0000 0.1681 0 1.0000 1.7264 0 1.0000 -0.6681 1.4873 1.0000 -0.4764 0.0724 g = 8 r =10.4743 -9.2372 -33.8809i -9.2372 +33.8809i p =-1.7264 0.2382 + 0.1252i 0.2382 - 0.1252i k =16 FIR数字滤波器的基本结构与实现 横截型与级联型之间的转换 已知一个FIR系统的传递函数为 将其从横截型转换为级联型。信号流图如下所示。 将横截型转换为级联型，就是将系统函数转换为二次分式模型。 程序运行结果： sos =1.0000 0.9500 0 1.0000 0 0 1.0000 -0.5000 1.2500 1.0000 0 0 g =2 级联型的表达式 五、实验内容 输入并运行例题程序，熟悉基本指令的使用。 已知一个系统的级联型表达式为： 信号流图如图所示 将其从级联型转换位直接型和并联型结构，并分别绘制直接型和并联型的信号流图。 六、实验报告要求 简述实验目的和实验原理。 列写练习题的代码并绘制程序产生的图形。 总结实验中你的收获和体会。 * * 实验四 数字滤波器的结构 以上两个公式中，当ak至少有一个不为0时，则在有限z平面上存在极点，表达的是一个IIR数字滤波器；当ak全为0时，系统不存在极点，表示的是一个FIR数字滤波器。FIR滤波器可以看成是IIR数字滤波器的ak全为0的一个特例。 MATLAB程序： b=[8,-4,11,2]; a=[1,1.25,-0.75,0.125]; [sos,g]=tf2sos(b,a) [r,p,k]=residuez(b,a) 级联型： 并联型： MATLAB程序： b=[2,0.9,1.55,2.375]; a=[1]; [sos,g]=tf2sos(b,a) *