基于Blackman窗的FIR数字带通滤波器设计25654.docx

《基于Blackman窗的FIR数字带通滤波器设计25654.docx》由会员分享，可在线阅读，更多相关《基于Blackman窗的FIR数字带通滤波器设计25654.docx（16页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、课 程 设 计课程设计名称:基于Blackman窗的FIR数字带通滤波器设计专 业 班 级 ： 学 生 姓 名 ： 学 号 ： 指 导 教 师 ： 课程设计时间： 2010.6.21-2010.6.25 数字信号处理 专业课程设计任务书学生姓名专业班级学号题 目基于Blackman窗的FIR数字带通滤波器设计课题性质工程技术研究课题来源自拟课题指导教师同组姓名主要内容 设计一个具有如下性能指标的FIR数字带通滤波器：=0.45，=0.65，=0.30，=0.80；窗函数为Blackman窗；分析最后设计结果性能。任务要求1. 温习窗函数法设计数字FIR滤波器的原理和设计方法；2. 求出最后所设

2、计出的滤波器冲激响应系数；3. 绘制所设计滤波器的增益响应曲线。参考文献1 胡广书. 数字信号处理理论、算法与实现M. 北京: 清华大学出版社, 1997.2 R. Lyons. Understanding Digital Signal Processing M. 2nd ed. Prentice Hall PTR., 2004.3 A.V.奥本海姆, R.W.谢弗 and J.R.巴克. 离散时间信号处理M. 第二版. 西安交通大学出版社, 2001.4 S. K. Mitra. Digital Signal Processing: A Computer-Based ApproachM. 3

3、rd ed. McGraw-Hill, 2005.审查意见指导教师签字：李相国教研室主任签字： 年 月 日 说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页1 需求分析数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（iir）滤波器和有限长冲激响应（fir）滤波器。iir数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以iir滤波器的设计可以采取在模拟滤波器设计的基

4、础上进一步变换的方法。fir数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。因此设计FIR滤波器的方法之一可以从时域出发，截取有限长的一段冲击响应作为H(z)的系数，冲击响应长度N就是系统函数H(z)的阶数。只要N足够长，截取的方法合理，总能满足频域的要求。一般这种时域设计、频域检验的方法要反复几个回合才能成功。要设计一个线性相位的FIR数字滤波器，首先要求理想频率响应。是w的周期函数，周期为，可以展开成傅氏级数： 其中是与理想频响对应的理想单位抽样响应序列。但不能用来作为设计FI

5、R DF用的h(n)，因为一般都是无限长、非因果的，物理上无法实现。为了设计出频响类似于理想频响的滤波器，可以考虑用h(n)来近似。窗函数的基本思想：先选取一个理想滤波器（它的单位抽样响应是非因果、无限长的），再截取（或加窗）它的单位抽样响应得到线性相位因果FIR滤波器。这种方法的重点是选择一个合适的窗函数和理想滤波器。设x(n)是一个长序列，w(n)是长度为N的窗函数，用w(n)截断x(n)，得到N点序列xn(n)，即xn(n) = x(n) w(n) 在频域上则有 由此可见，窗函数w(n)不仅仅会影响原信号x(n)在时域上的波形，而且也会影响到频域内的形状。MATLAB信号工具箱主要提供了

6、以下几种窗函数，如表1-1所示：表1-1 MATLAB窗函数窗窗 函 数矩形窗Boxcar巴特利特窗Barlett三角窗Triang布莱克曼窗Blackman海明窗Hamming汉宁窗Hanning凯塞窗Kaiser切比雪夫窗Chebwin加矩形窗后的频谱和理想频谱可得到以下结论：加窗使过渡带变宽，过渡带的带宽取决于窗谱的主瓣宽度。矩形窗情况下的过渡带宽是。N越大，过渡带越窄、越陡；过渡带两旁产生肩峰，肩峰的两侧形成起伏振荡。肩峰幅度取决于窗谱主瓣和旁瓣面积之比。矩形窗情况下是8.95，与N无关。工程上习惯用相对衰耗来描述滤波器，相对衰耗定义为 这样两个肩峰点的相对衰耗分别是0.74dB和-2

7、1dB。其中（-0.0895）对应的点的值定义为阻带最小衰耗。以上的分析可见，滤波器的各种重要指标都是由窗函数决定，因此改进滤波器的关键在于改进窗函数。窗函数谱的两个最重要的指标是：主瓣宽度和旁瓣峰值衰耗。旁瓣峰值衰耗定义为： 旁瓣峰值衰耗20lg(第一旁瓣峰值主瓣峰值) 为了改善滤波器的性能，需使窗函数谱满足：主瓣尽可能窄，以使设计出来的滤波器有较陡的过渡带；第一副瓣面积相对主瓣面积尽可能小，即能量尽可能集中在主瓣，外泄少，使设计出来的滤波器的肩峰和余振小。但上面两个条件是相互矛盾的，实际应用中，折衷处理，兼顾各项指标。本次设计要求以窗函数为Blackman窗为基础，设计一个具有如下性能指标

8、的FIR数字带通滤波器：=0.45，=0.65，=0.30，=0.80。要求：1. 温习窗函数法设计数字FIR滤波器的原理和设计方法；2求出最后所设计出的滤波器冲激响应系数； 3绘制所设计滤波器的增益响应曲线。2 概要设计如果所希望的滤波器的理想的频率响应函数为，则其对应的单位脉冲响应为窗函数设计法的基本原理是用有限长单位脉冲响应序列逼近。由于往往是无限长序列，而且是非因果的，所以用窗函数将截断，并进行加权处理，得到：就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为式中，N为所选窗函数的长度。我们知道，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程

9、中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。主程序框图如图4.1所示。其中幅度特性要求用dB表示。开始读入窗口长度N计算hd(n)调用窗函数子程序求w(n)调用子程序（函数）计算H(k)=DFTh(n)调用绘图子程序（函数）绘制H(k)幅度相位曲线结束图1 主程序框图计算h(n)= hd(n) w(n)布莱克曼窗函数布莱克曼窗函数的时域形式可以表示为 它的频域特性为：其中，为矩形窗函数的幅度频率特性函数。布莱克曼窗函数的最大旁瓣值比主瓣值低57dB，但是主瓣宽度是矩形窗函数的主瓣宽度的3倍，为12/N。Blackman函数：生成海明窗调用方式(1) w = b

10、lackman (n)：输入参数n是窗函数的长度；输出参数w是由窗函数的值组成的n阶向量。(2) w = blackman (n,sflag)：参数sflag用来控制窗函数首尾的两个元素值；其取值为symmetric或periodic；默认值为symmetric。3 运行环境软件环境：1.Windows xp2. MATLAB6.53其他软件根据需要选用硬件环境：计算机：Lenovo Y450处理器型号：Intel酷睿2双核 T6600 处理器参数：主频2.00GHZ主板芯片组：IntelGM45 内存：2GB DDRIII硬盘：西数WDC WD3200BEVT-22ZCTO SATA 显存：

11、512MB 4 开发工具和编程语言开发工具:MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。MATLAB一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂的应用程序（M文件）后再一起运行。编程语言:汇编语言汇编语言(AssemblyLanguage)是面向机器的程序设计语言。在汇编语合中，用助记符(Memoni)代替操作码，用地址符号(Symbol)或标号(Label)代替地址码。这样用符号代替机器语言的二进制码，就把机器

12、语言变成了汇编语言。于是汇编语言亦称为符号语言。使用汇编语言编写的程序，机器不能直接识别，要由一种程序将汇编语言翻译成机器语言，这种起翻译作用的程序叫汇编程序，汇编程序是系统软件中语言处理系统软件。汇编程序把汇编语言翻译成机器语言的过程称为汇编。汇编语言的特点: 1.面向机器的低级语言，通常是为特定的计算机或系列计算机专门设计的。 2.保持了机器语言的优点，具有直接和简捷的特点。 3.可有效地访问、控制计算机的各种硬件设备，如磁盘、存储器、CPU、I/O端口等。 4.目标代码简短，占用内存少，执行速度快，是高效的程序设计语言。 5.经常与高级语言配合使用，应用十分广泛。5 详细设计 算法实现的

13、源程序：（1）理想低通滤波器函数function db,mag,pha,grd,w=freqz_m(b,a);N=1000;H,w=freqz(b,a,N,whole);H=(H(1:1:501);w=(w(1:1:501);mag=abs(H);db=20*log10(mag+eps)/max(mag);pha=angle(H);grd=grpdelay(b,a,w);（2）定义滤波器系数向量函数function hd=ideallp(wc,M);al=(M-1)/2;n=0:(M-1);m=n-al+eps;hd=sin(wc*m)./(pi*m);（3）主程序：% 数字滤波器指标ws1=

14、0.30*pi; wp1=0.45*pi;ws2=0.80*pi; wp2=0.65*pi;As=60;tr_width=min(wp1-ws1),(wp2-ws2);M=ceil(11*pi/tr_width)+1;n=0:1:M-1;wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);% 生成blackman窗w_bla=(blackman(M);h=hd.*w_bla% 频域图像的绘制freqz(h,1)figure(2);subplot(2,2,1),stem(n,hd);title(idael impu

15、lse response)axis(0 M-1 -0.4 0.5);xlabel(n);ylabel(hd(n)xa=0.*n;hold onplot(n,xa,k);hold offsubplot(2,2,2),stem(n,w_bla);title(blackman window)axis(0 M-1 0 1.1);xlabel(n);ylabel(w(n)subplot(2,2,3),stem(n,h);title(actual impulse response)axis(0 M-1 -0.4 0.5);xlabel(n);ylabel(h(n)hold onplot(n,xa,k);h

16、old offn=0:600;k=n/2000;f1=2*pi*450;f2=2*pi*600;x=sin(f1*k)+sin(f2*k);y=filter(h,1,x);x1=sin(f1*k);x2=sin(f2*k);figure(2)subplot(2,2,1);plot(x1);title(x1的图形);grid;axis(0 100*pi -5 5);ylabel(x1);subplot(2,2,2);plot(x2);title(x2的图形);grid;axis(0 100*pi -5 5);ylabel(x2);subplot(2,2,3);plot(x);title(x的图形

17、);grid;axis(0 100*pi -5 5);ylabel(x);subplot(2,2,4);plot(y);title(y的图形);grid;axis(0 100*pi -5 5);ylabel(y);6 调试分析设计带通滤波器时首先要计算出过渡带，然后查表得到所设计窗函数所需要的阶数，不同的窗函数所设计的滤波器的形状各有差异，尤其在主瓣宽度、旁瓣的形状以及主瓣与旁瓣的高度差上有比较明显得差别，实际应用中应根据实际情况，折衷处理，兼顾各项指标。实际中遇到的离散时间信号总是有限长的，因此不可避免地要遇到数据截断问题。而在信号处理中，对离散序列的数据截断是通过序列与窗函数相乘来实现的。

18、7 测试结果测试结果有错，初始参数为ws1=0.35*pi;wp1=0.4*pi;wp2=0.5*pi;ws2=0.55*pi改为要求参数ws1=0.30*pi;wp1=0.45*pi;wp2=0.65*pi;ws2=0.80*pi时，运行出错，错误在主程序第12行：“Rp=-min(db(ws1/delta_w+1:1:wp2/delta_w)：”，由于时间不充分的原因，最后没调试出要求参数。运行结果(图):图1 布莱克曼窗函数波形图图2 布莱克曼窗函频谱形图3布莱克曼窗设计带通滤波器的幅频响应图图4 布莱克曼窗设计带通滤波器的相频响应图图5布莱克曼窗设计带通滤波器的冲击响应图参考文献1刘泉

19、 信号与系统 高等教育出版社 20062 陈怀琛MATLAB及在电子信息课程中的应用(第3版) 电子工业出版社20063 程佩青编. 数字信号处理教程. 北京：清华大学出版社，20014 丁玉美，高西全编著. 数字信号处理.西安：西安电子科技大学出版社，20005刘泉 数字信号处理原理与实现 电子工业出版社 20056 黄文梅，熊桂林，杨勇编著. 信号分析与处理-MATLAB语言及应用.长沙：国防科技大学出版社7 张宗橙，张玲华，曹雪虹编著. 数字图像处理与应用。南京：东南大学出版社，1997 心得体会通过对该课程的学习，我对数字通信理论有了更进一步的理解；通过对该数字滤波器的设计，我也了解了

20、数字滤波器的基本结构和基本特性，而且还掌握了基本的撰写论文的形式和思路。作为一个电子信息工程专业的学生，数字信号处理是我们的重要专业课程，是我们将来从事通信事业的基本保障。通过对作为该论文的重要部分MATLAB的运用，大大提高了我们对集计算，编程与绘图于一体的该应用软件的运用能力。MATLAB包含的几十个工具箱，覆盖了通信，自动控制，信号处理，图象处理，财经，化工，生命科学等科学技术领域，汲取了当今世界这些领域的最新研究成果，已经成为从事科学研究和工程设计不可缺少的工具软件。该课程设计将数字信号处理的有关教学内容和MATLAB语言紧密，有机地结合起来，使我们在学习基础理论知识的同时学会了应用MATLAB，在学习应用MATLAB的同时，加深了对基本知识的理解，增强了我们的计算机应用能力，提高了学习效果。总之，无论是从教学知识掌握出发，还是从对MATLAB的应用出发，通过这次学习，我不但掌握了基于Blackman窗的FIR数字带通滤波器设计的基本知识及其实际应用的技巧，还提高了自己的编程和写报告的能力，收获不小,也巩固了所学知识。