北交大数字信号研讨实验2.doc

《数字信号处理》课程研究性学习报告 姓名 陶威东 学号 同组成员 张欣悦 指导教师 刘留 时间 2015 年6月1日 数字滤波器设计专题研讨 【目的】 (1) 掌握IIR和FIR数字滤波器的设计方法及各自的特点。 (2) 掌握各种窗函数的时频特性及对滤波器设计的影响。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 1．分析矩形窗、汉纳窗、哈明窗、布莱克曼窗、凯泽窗的频域特性，并进行比较。 【题目分析】 不同窗函数的主瓣宽度，主瓣高度以及过渡带长度和阻带波动都有明显差别。下面我们将通过对时域长度相同的窗函数进行分析来比较不同窗函数的频域特性。取时域长度为50 【仿真结果】 【结果分析】 各种窗有何特点? 通带衰减由小到大的排列为：ｂｌａｃｋｍａｎ，ｈａｍｍｉｎｇ，ｋａｉｓｅｒ； 阻带衰减由小到大的排列为：ｋａｉｓｅｒ，ｈａｍｍｉｎｇ，ｂｌａｃｋｍａｎ； 过渡带由小到大的排列为：ｋａｉｓｅｒ，ｈａｍｍｉｎｇ，ｂｌａｃｋｍａｎ； 在设计滤波器的时候，我们希望得到最大的通带衰减，最小的阻带衰减及最小的过渡带。但在实践 中，这三者是不可能同时满足的。所以，我们要根据实际需求，选择合适的窗函数，以达到最佳的 滤波效果。最后我们还可以看出，主瓣宽度与主瓣幅值的乘积是定值，即主瓣宽度的增加必将导致主瓣幅值的降低。最后我们会发现凯泽窗幅度响应曲线和矩形窗的完全吻合。这是因为凯泽窗是是一个可以变幻形状的窗函数，当系统没有给值时，默认的为矩形窗，所以是完全符合的。 【自主学习内容】 几种不同窗函数的设计 【阅读文献】 [1] 陈后金.数字信号处理[M].北京：高等教育出版社.2008.11 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 对于不同窗函数的选择上，我们需要根据需要在主瓣宽度和过渡带长短之间做出选择 【问题探究】 在谱分析中如何选择窗函数，在滤波器设计中如何选择窗函数？ 提高FIR滤波器阻带衰减是以增加过渡带宽度为代价的。实际中，在满足阻带衰减的前提下，尽可能地选择主瓣宽度较小的窗函数。凯泽窗是可调的，它通过改变窗函数的形状来控制窗函数旁瓣的大小，设计中可根据滤波器的衰减指标来确定窗函数的形状。 【仿真程序】 M=50;L=512; fsam=1000;wsam=2*pi*fsam; w=(-wsam/2+(0:L-1)*wsam/L)/fsam; w1=ones(1,M);X1=fftshift(fft(w1,L));plot(w,abs(X1),g); hold on;xlabel(Normalized frequency); legend(square); title(·幅度响应|);grid on; w2=hanning(M); X2=fftshift(fft(w2,L)); subplot(2,2,1),plot(w,abs(X2),r); hold on;xlabel(Normalized frequency); legend(hanning); title(·幅度响应);grid on; w3=hamming(M); X3=fftshift(fft(w3,L)); subplot(2,2,2),plot(w,abs(X3),y); hold on;xlabel(Normalized frequency); legend(hamming); title(幅度响应|);grid on; w4=blackman(M); X4=fftshift(fft(w4,L)); subplot(2,2,3),plot(w,abs(X4),b); hold on;xlabel(Normalized frequency); legend(blackman); title(幅度响应);grid on; w5=kaiser(M); X5=fftshift(fft(w5,L)); subplot(2,2,4),plot(w,abs(X5),--); xlabel(Normalized frequency); legend(kaiser); title(幅度响应);grid on