数字信号处理课设报告课案.docVIP

数字信号处理课程设计 题目 无限冲激响应(IIR)数字滤波器的直接型 专业 电子信息工程 姓名 卢想 学号 20130210410117 时间 2016年 6 月 22日 目 录 一、课程设计目的 1 二、课程设计原理与计算方法 2 三、课程设计内容 3 四、程序与调试 4 4.1激励信号1. 4 4.2.激励信号2 5 4.3.扩展函数impseq 6 五、课程设计总结 7 参考文献 7 一、课程设计目的 1、掌握IIR滤波器的直接II型、级联型和并联型三种结构的基本原理和特点。 2、掌握利用Matlab实现IIR滤波器的三种结构的程序设计方法，并能够进行三者之间的相互转换。 3、掌握滤波器频响特性的绘制方法 二、课程设计原理与计算方法 按照结构划分数字滤波器，有递归式和非递归式两种。递归式数字滤波器的差分方程为 　　　　　　　　　　　　　　　　　　　　 (1) 其中至少有一个.非递归式数字滤波器的差分方程为 (2) 可以看出递归式数字滤波器的响应不仅与激励有关，而且与以前的输出信号有关；而非递归式数字滤波器的响应仅只与激励有关。 按照单位样值响应划分数字滤波器，则有无限冲激响应(IIR)和有限冲激响应(FIR)之分。IIR滤波器是递归式的，差分方程如(1)式所示，FIR滤波器一般是非递归式的，差分方程如(2)式所示。 １、直接II型（也称为正准型结构） 根据(1)式，IIR滤波器的传输函数为 　　　　　　　　　　　　　　　　　　　　　　　　　 (3) 其中已假设(1)式中的，对于其它情况，则可令相应的某些系数为零。 令　　　　　　 则有　　　　　　　　　　　　　　　　　　　　　　　　　　(4) 由此可以得到相应的时域中激励与响应之间的关系为 (5) 其中是与(4)式中的相应的中间函数序列。 x(n) y2(n) b0 y(n) z-1 a1 b1 z-1 a2 b2 z-1 aN-1 bN-1 z-1 aN bN 　　图 6-1 直接II型结构信号流 由(5)式确定的直接II型的信号流图如图6-1所示，其中将中间的两条延时链合并为一条，实际的信号流将按（5）式分成两个延时链独立运行。编程时，用三个数组分别存放系数和，。由图可以看出，沿中间的延时链自上向下传播过程中将逐级向右移位，每一级向左边与的乘积按累加，再和相加，得到的沿中间的延时链又自上向下传播逐级向右移位，每一级向右边与的乘积按累加，其结果就是对应于激励的响应。然后向右移位一个单位时间，输入激励计算下一个时刻的响应。 直接II型结构具有简单直观的典型网络结构形式，在计算机上很容易实现。但是它对系数的量值变化比较敏感，直接确定了系统零、极点的位置，从而影响到系统的性能。尤其当阶数N较高时，系统对系数的字长效应很敏感，产生的误差也较大。 三、课程设计内容 设三阶滤波器的传输函数为 　　 激励信号为 　　 　　 根据所给定的滤波器系统函数和参数值，用直接型实 IIR数字滤波器结构，分别使激励和通过该滤波器，求出相应于不同激励的响应，作出输出响应信号的时间曲线和幅频、相频特性曲线，并判断为何种滤波器(低通、高通、带通)。 四、程序与调试 4.1激励信号1. b=[1?0.5814?0.2114?0]；a=[1?-0.3984?0.2475?-0.04322]；delta=impseq(0,-10,10)；hdir=filter(b,a,delta)；n=[0:40]；x1=exp(-0.5.*n).*cos(n.*(pi./2))；y1=filter(b,a,x1)；subplot(2,2,1)；stem(x1,.);title(x1)subplot(222);plot(y1);title(直接型)subplot(223);plot(angle(y1));title(直接型相频特性)subplot(224);plot(abs(y1));title(直接型幅频特性) 运行图1