华南农业大学期末考试试卷(A卷)2011学年第1学期考试科目:数字信号处理考试类型:(闭卷)考试考试时间:120分钟学号姓名年级专业题号一二三四………总分得分评阅人一、判断题(本大题共10小题,每小题2分,共20分)1.任何序列都可以用单位脉冲序列的移位加权和来表示。(√)2.序列3()2cos(7)4xnn是周期序列。(√)3.序列y(n)=2x(n)-3,不是移不变系统。(×)4.一个域的离散就必然导致另一个域的连续。(×)5.因果序列的收敛域不包含。(×)6.FIR滤波器设计可利用模拟滤波器设计的结果。(×)7.最小相位延时系统的零点和极点都在单位圆之内。(√)8.若滤波器通带内群延时响应特性是一个常数,则为线性相位系统。(√)9.窗函数的选择原则是在保证阻带衰减的情况下选择主瓣窄的窗函数。(√)10.全通系统的特点是零极点以单位圆镜像对称。(√)二、填空题(本大题共5小题,每小题4分,共20分)1.序列()()xnun的Z变换为1111zzz,)3(nx的Z变换是23111zzzz得分得分2.线性移不变系统是稳定系统的充要条件是Pnhn)(。3.实现FIR线性相位滤波器的条件是h(n)=(1)hNn4.序列x(n)和h(n),长度分别为N和M(NM),二者线性卷积的长度为N+M-1循环卷积与线性卷积的关系是。5.实现一个数字滤波器需要的三种基本运算单元是加法器,乘法器,延迟器。三、简答题(本大题共6小题,每小题5分,共30分)1.简述Z变换、傅里叶变换、离散傅里叶变换三者之间的关系。答:傅氏变换是拉氏变换在虚轴S=jΩ的特例,因而映射到Z平面上为单位圆,即序列在单位圆上的Z变换为序列的傅氏变换;x(n)的N点离散傅里叶变换是x(n)的z变换在单位圆上的N点等间隔抽样;也是是傅里叶变换在区间[0,2π]上的N点等间隔抽样。2.简述按频率抽选法和按时间抽选法两种FFT算法的异同?(包括输入输出顺序、基本碟形、计算量、节点间距离、rWN因子确定等异同)答:输入输出顺序基本碟形计算量节点间距离rWN(蝶形运算两节点的第一个节点为k值,表示成L位二进制数。)DIT输入倒序,输出自然顺序先复乘后加减相同:12m左移L–m位,把右边空出的位置补零,结果为r的二进制数。DIF输入自然顺序,输出倒序先加减后复乘2/2LmmN左移m-1位,把右边空出的位置补零,结果为r的二进制数。3.设某FIR数字滤波器的系统函数为:请画出此滤波器的线性相位结构。得分1.5CM1234113535Hzzzzz()()NclNynyn点圆周卷积是线性卷积以为周期的周期延拓序列的主值序列。解:对系统函数求z反变换,得得,该滤波器单位脉冲响应h(n)偶对称,N为奇数。得线性相位结构如图:4.简述全通系统的特点以及应用。答:全通系统的特点是:零点与极点以单位圆为镜像对称。其应用主要有三方面:1)任一因果稳定系统H(z)都可以表示成全通系统Hap(z)和最小相位系统Hmin(z)的级联2)级联一个全通系统可以使非稳定滤波器变成一个稳定滤波器3)作为相位均衡器,校正系统的非线性相位,而不改变系统的幅度特性5.简述IIR以及FIR数字滤波器的区别。(提示:滤波器特点、设计方法及应用的区别)答:1)从性能上来说,IIR滤波器传输函数的极点可位于单位圆内的任何地方,因此可用较低的阶数获得高的选择性,所用的存贮单元少,所以经济而效率高。但是这个高效率是以相位的非线性为代价的。选择性越好,则相位非线性越严重。相反,FIR滤波器却可以得到严格的线性相位,然而由于FIR滤波器传输函数的极点固定在原点,所以只能用较高的阶数达到高的选择性;对于同样的滤波器设计指标,FIR滤波器所要求的阶数可以比IIR滤波器高5~10倍,结果,成本较高,信号延时也较大;如果按相同的选择性和相同的线性要求来说,则IIR滤波器就必须加全通网络进行相位较正,同样要大增加滤波器的阶数和复杂性。2)从机构上看,IIR滤波器必须采用递归结构,极点位置必须在单位圆内,否则系统将不稳定。相133112345555hnnnnnn反,FIR滤波器主要采用非递归结构,不论在理论上还是在实际的有限精度运算中都不存在稳定性问题,运算误差也较小。此外,FIR滤波器可以采用快速付里叶变换算法,在相同阶数的条件下,运算速度可以快得多。3)从设计工具上看,IIR滤波器可以借助于模拟滤波器的成果,因此一般都有有效的封闭形式的设计公式可供准确计算,计算工作量比较小,对计算工具的要求不高。FIR滤波器设计则一般没有封闭形式的设计公式。窗口法虽然仅仅对窗口函数可以给出计算公式,但计算通带阻带衰减等仍无显式表达式。一般,FIR滤波器的设计只有计算程序可循,因此对计算工具要求较高。另外,IIR滤波器虽然设计简单,但主要是用于设计具有片段常数特性的滤波器,如低通、高通、带通及带阻等,往往脱离不了模拟滤波器的格局。而FIR滤波器则要灵活得多,尤其它能易于适应某些特殊的应用,如构成微分器或积分器,或用于Butterworth、Chebyshev等逼近不可能达到预定指标的情况,例如,由于某些原因要求三角形振幅响应或一些更复杂的幅频响应,因而有更大的适应性和更广阔的天地。6.简述FIR窗函数设计法中,窗函数的选择依据。并分析矩形窗、三角形窗(Bartlett窗)、汉宁窗、海明窗以及布拉克曼窗的特点。答:窗函数的选择依据是:1)希望窗谱的主瓣尽量地窄,以获得较陡的过渡带;2)尽量减少窗谱最大旁瓣的相对幅度,也就是使频域的能量能主要集中在主瓣内。这样使肩峰和波纹减少,就可增大阻带的衰减。相同的阶数下,矩形窗的窗谱主瓣宽度最窄,为N/4,旁瓣幅度最大;三角形窗的窗谱主瓣宽度是矩形窗的两倍,为N/8,旁瓣幅度较矩形窗小;汉宁窗的窗谱主瓣宽度是矩形窗的两倍,为N/8,旁瓣幅度较三角形窗小;海明窗是改进的汉宁窗,其窗谱主瓣宽度也是矩形窗的两倍,为N/8,旁瓣幅度较汉宁窗小;布拉克曼窗的窗谱主瓣宽度是矩形窗的三倍,为N/12,旁瓣幅度最小。因此。用矩形窗设计FIR滤波器时,过渡带最窄,而阻带衰减最小,布莱克曼窗过渡带最宽,但阻带衰减加大。四、分析计算题(本大题共3小题,每小题10分,共30分)1.设一因果系统的传递函数为11210.5()10.70.12zHzzz,(1)系统是否稳定?为什么?(2)写出差分方程。(3)画出系统的极零点分布图。得分1.5CM解:(1)由H(z)的表达式可知,系统有两个极点为:z1=0.3,z2=0.4,所以极点都在单位圆内,因此系统稳定。(2)11210.5()10.70.12zHzzz=112110.310.4zz由于该系统是一个因果,所以,()(2(0.3)(0.4))()nnhnun(3)由11210.5()10.70.12zHzzz,可知b0=1;b1=--0.5;a0=1;a1=0.7;a2=-0.12则差分方程可表示为:()()0.7(1)0.12(2)0.5(1)ynxnxnxnyn(4)零极点分布如图所示:2.下图表示一个5点序列()xn(1)试画出线性卷积()()xnxn(2)试画出5点的圆周卷积()xn⑤()xn;解:()()xnxn3.要求用双线性变换法从三阶巴特沃思模拟滤波器导出一低通数字滤波器,已知3dB截止频率为5Hz,系统抽样频率为1kHz,设T=2s。表巴特沃思滤波器分母多项式12212211NNNNNsasasasas的系数Na1a2a3a4a512345611.414222.61313.23603.863723.41425.23607.46412.61315.23609.14163.23607.46413.8637解:由题意可得3dB截止频率为5Hz,则1025c。去归一化,得322333322322020010001000100050511)10()10(210211)()(ssssssssssHsHcana将()aHs变换成Butterworth数字滤波器:1111112322331122020010001000)()(zzTszzTsassssHzH31121111233)11()11(20)11(20010001000zzzzzz⑤()xn()xnL点圆周卷积是线性卷积以L为周期的周期延拓序列的主值区间