《数字信号处理》课程设计报告-DFT在信号频谱分析中的应用窗函数法设计FIR数字低通滤波器.doc

《数字信号处理》课程设计报告 任课教师： 指导教师： 学生学号： 6 学生姓名： 所学专业：电子信息工程 2011年06　月22　日 目录 设计题目 设计目的 设计原理 实现方法（包括MATLAB算法原理等） 设计内容及要求（含有设计源程序） 设计结果及改进建议（画出所有设计曲线，并加以说明） 回答思考题 设计体会 参考文献 设计题目 设计一、DFT在信号频谱分析中的应用 设计八 、窗函数法设计FIR数字低通滤波器 二、设计目的 （一）设计一、DFT在信号频谱分析中应用的原理 1. 熟悉DFT的性质。 2. 加深理解信号频谱的概念及性质。 3. 了解高密度谱与高分辨率频谱的区别。 （二）设计八 、窗函数法设计FIR数字低通滤波器 1. 熟悉设计线性相位数字滤波器的一般步骤。 2. 掌握用窗函数法设计FIR数字滤波器的原理和方法。 3. 熟悉各种窗函数的作用以及各种窗函数对滤波特性的影响。 4. 学会根据指标要求选取合适的窗函数。设计低通FIR滤波器的指标： 通带最大波动 ， 阻带最小衰减 ， 三、设计原理 （一）设计一、DFT在信号频谱分析中应用的原理 所谓信号的频谱分析就是计算信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于直接用计算机进行计算，使其应用受到限制，而DFT是一种时域和频域均离散化的变换，适合数值运算，成为分析离散信号和系统的有力工具。 工程实际中，经常遇到的连续信号Xa(t),其频谱函数Xa(jW)也是连续函数。数字计算机难于处理，因而我们采用DFT来对连续时间信号的傅里叶变换进行逼近，进而分析连续时间信号的频谱。 离散傅里叶变换（DFT）定义：设有限长序列x (n) 长为N（0nN-1），其离散傅里叶变换是一个长为N的频率有限长序列（0kN-1），其正变换为 0kN-1 （） 离散傅里叶变换的实质是：把有限长序列当做周期序列的主值序列进行DFS变换，x(n)、X(k)的长度均为N，都是N个独立值，因此二者具有的信息量是相等的。已知x(n)可以唯一确定X(k),已知X(k)可以唯一确定x(n)。 虽然离散傅里叶变换是两个有限长序列之间的变化，但它们是利用DFS关系推导出来的，因而隐含着周期性。 （二）设计八 、窗函数法设计FIR数字低通滤波器的原理 1、设计步骤 （1）给定所求的频率响应函数Hd(ejw); （2）求hd(n)=IDTFT[Hd(ejw)]; （3）由过渡带宽及阻带最小衰减的要求，选定窗w(n)的形状及N的大小，一般N要通过几次试探而最后确定； （4）求得所设计的FIR滤波器的单位抽样响应h(n)=hd(n)w(n),n=0,1,…，N-1； （5）求H（ejw）=DTFT[h(n)],检验是否满足设计要求，如不满足，则需重新设计。 2、理想滤波器的频率响应函数其对应单位脉冲响应为：。窗函数设计法的基本原理是用有限长单位脉冲响应h（n）逼近 ，用窗函数w（n）将它截断，并进行加权处理。 3、h(n)就是实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数为：。如果要求线性相位特性，则h(n)还必须满足：，根据上式中的正负号和长度N的奇偶性又将线性相位滤波器氛围四类。要根据所设计的滤波特性选择其中一类。例如要设计线性相位低通特性，可以选择h(n)=h(N-1-n)这一类，而不能选择h(n)=-h(N-1-n)这一类。 4、四种窗口函数 （1）矩形窗 （2）汉宁窗（升余弦窗） =0.5 三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍。 汉明窗（改进的升余弦窗） ??? 它是对汉宁窗的改进，在主瓣宽度（对应第一零点的宽度）相同的情况下，旁瓣进一步减小，可使99.96%的能量集中在主瓣内。 （4）布莱克曼窗（三阶升余弦窗） 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，增加N可减少过渡带。 四、实现方法（包括MATLAB算法原理等） （一）设计一、DFT在信号频谱分析中应用的实现方法 1. 用MATLAB语言编写计算序列x(n)的N点DFT的m函数文件dft.m。并与MATLAB中的内部函数文件fft.m作比较。 2. 对离散确定信号 作如下谱分析： 截取使成为有限长序列N()，(长度N自己选)写程序计； 算出的N点DFT ,画出时域序列图xn～n和相应的幅频图。 (2) 将 (1)中补零加长至M点，长度M自己选,（为了比较补零长短的影响，M可以取两次值，一次取较小的整数，一次取较大的整数），编写程序计算的M点DFT, 画出时域序列图和两次补零后相应的DFT幅频图。 (3)利用补零DFT计算 (1)中N点有限长序列频谱并画出相应的幅频图