《数字信号处理》实验4 数字滤波器的结构实现.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》 实验目的：本实验包含2个实验子项目，通过matlab编程仿真，掌握IIR/FIR数字滤波器的基本结构及其实现方法。1. 掌握IIR数字滤波器的直接型、级联型及并联型结构及实现方法。2. 掌握FIR数字滤波器的直接型、级联型及频率采样型结构及实现方法。实验四 数字滤波器的结构实现 实验四 数字滤波器的结构实现4.1 IIR数字滤波器的结构实现4.2 FIR数字滤波器的结构实现 4.1 IIR数字滤波器的结构实现一、实验原理1. 无限脉冲响应(IIR)数字滤波器的直接型结构 N阶IIR滤波器的系统函数为： (4.1)其中：ai(i=0~N),bi(i=0~M)----常数系数。直接型结构求系统的输出：Matlab函数：y=filter(b,a,x)其中：x—输入信号；y---输出信号；b,a---H(z)分子/分母多项式系数：b=[b0,b1, b2，…,bM]; a=[ 1, a1, a2，…,aN]; 2. IIR滤波器的级联型结构二阶基本级联结构的IIR滤波器系统函数为： (4.2)(1)直接型结构转换为二阶基本级联结构Matlab函数：sos=tf2sos(b,a)其中：sos---级联型结构二次分式的系数；b,a---直接型结构H(z)分子/分母多项式系数。。(2)计算级联型结构的输出y(n)Matlab函数： y= sosfilt(sos,x) 3. IIR滤波器的并联型结构 (4.3)直接型结构转换为并联结构:Matlab函数：[r,p,k]=residuez(b,a)其中： b,a---H(z)分子/分母多项式系数；r---部分分式系数向量；p—极点向量；k---多项式系数向量。二、实验环境1.计算机1台。2. Windows 7 以上操作系统。3. Matlab 7.0 以上版本软件。 三、实验参考与实验内容1. IIR数字滤波器的直接型结构设：某系统的系统函数H(z)=求：系统单位脉冲响应h(n)及零极点分布图。实验参考程序：参见: 实验项目参考源码\ch4\ch4prog1.m 实验内容：(1)运行此程序，画出h(n)及H(z)的零极点分布图。(2)观察零极点分布图，问：系统极点d1,d2,d3=？零点c1,c2,c3=？ 图4.1 IIR滤波器的直接型结构 2. IIR滤波器的级联型结构实现设：某系统的系统函数H(z)=，将它转换为级联型结构。实验参考程序：参见: 实验项目参考源码\ch4\ch4prog2.m 实验内容：(1)在空格中填入注释。(2)运行此程序，在命令窗口中读取sos，则：sos=?(3)根据sos矩阵，写出H(z)级联型结构的表达式。3. IIR滤波器的并联型结构实现设：某系统的系统函数H(z)=，将它转换为并联型结构。实验参考程序：参见: 实验项目参考源码\ch4\ch4prog3.m 实验内容：(1)在空格中填入注释。(2)运行此程序，在命令窗口中读取r,p,k=?(3)根据r,p,k矩阵，写出H(z)并联型结构的表达式。 四、实验要求1. 简述实验目的。2. 预习实验原理。3. 实验结果及分析：包括注明程序注释、画出实验运行结果波形、回答实验中提出的问题、如果有程序设计要求的，列出所编写的程序清单并简要叙述程序调试过程。返回 4.2 FIR数字滤波器的结构实现一、实验原理1. 有限脉冲响应(FIR)数字滤波器的直接型结构N阶FIR滤波器的系统函数为： (4.4)其中： h(n)---FIR滤波器单位脉冲响应；N---FIR滤波器的阶数。FIR滤波器直接型结构求系统的输出：Matlab函数：y=filter(b,a,x)其中：x—输入信号；y---输出信号；b,a---H(z)分子/分母多项式系数：b=[h(0) h(1) h(2) … h(N-1)]; a=[1