matlab结合矩形窗设计fir滤波器.doc
通信专业课程设计一太原科技大学课程设计(论文)设计论文题目MATLAB结合矩形窗设计FIR滤波器姓名___学号___班级__学院___指导教师________2011年1月17日太原科技大学课程设计(论文)任务书学院(直属系)电子信息工程学院时间2010年12月31日学生姓名指导教师设计(论文)题目MATLAB结合矩形窗设计FIR滤波器主要研究内容1、了解基于矩形窗的FIR滤波器的原理、方法、性能。2、了解MATLAB的基本功能,并设计出不同特性的矩形窗滤波器。3、分析信号滤波前后的时域和频域波形,并能与理论分析结果进行对比。研究方法在MATLAB环境下利用矩形窗函数法进行信号分析与滤波器的设计。主要技术指标或研究目标产生一个连续信号,包含低频,中频,高频分量,对其进行采样,进行频谱分析,使用矩形窗设计不同特性的数字滤波器对信号进行滤波处理。教研室意见教研室主任(专业负责人)签字年月日MATLAB结合矩形窗设计FIR滤波器I目录摘要II第1章绪论1第2章FIR滤波器设计原理221FIR滤波器222窗函数223矩形窗3第3章仿真设计与结果分析531设计流程图532产生原始信号并分析频谱633使用矩形窗设计不同特性的数字滤波器834信号滤波处理1035结果分析14第4章总结17参考文献17附录18MATLAB结合矩形窗设计FIR滤波器IIMATLAB结合矩形窗设计FIR滤波器摘要MATLAB是矩阵实验室MATRIXLABORATORY之意。她不仅具备卓越的数值计算能力,还提供了专业水平的符号计算,文字处理,可视化建模仿真和实时控制等功能。本文主要基于MATLAB平台,实现信号分析与滤波器的设计。并改变滤波器参数或特性低通、高通、带通或带阻,实现不同的滤波要求。本设计产生一个连续信号,包含低频,中频,高频分量,对其进行采样,进行频谱分析,使用矩形窗设计不同特性的数字滤波器对信号进行滤波处理,分析所设计滤波器画出了频率特性曲线,并对信号进行滤波,画出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化,分析的结果表明初步实现了设计目标。关键词MATLAB,信号分析,FIR滤波器,矩形窗MATLAB结合矩形窗设计FIR滤波器1第1章绪论数字信号处理技术及其应用,目前正以惊人的速度向前发展着。随着大规模集成电路的出现和数字部件的成本下降、体积缩小及运算速度提高,数字信号处理的应用日益广泛。目前已制成多种专用数字滤波器,取样率可高达兆赫。高速专用快速傅里叶变换处理机已有商品出售。简单的数字滤波器已制成集成电路片。目前几乎所有的语音带宽压缩系统都倾于全数字化,因为目前它是最实际可行的方法。除了专用数字信号处理硬件有所发展之外,还出现了可编程序的数字信号处理专用计算机,这种计算机的构造特别适于解决数字信号处理问题。它目前应用于实时信号处理以及设计、模拟专用数字硬件【1】。数字信号处理的重要性仍在不断提高,毫无停滞的迹象。其主要研究用数字序列或符号序列表示信号,并用数字计算方法对这些序列进行处理,以便把信号变换成符合某种需要的形式。数字信号处理的主要内容包括频谱分析、数字滤波与信号的识别等。而数字滤波器大略可分为有限冲激响应型和无限冲激响应型两类,可用硬件和软件两种方式实现。在硬件实现方式中,它由加法器、乘法器等单元所组成,这与电阻器、电感器和电容器所构成的模拟滤波器完全不同,而软件上可以借助MATLAB来仿真【2】。本文主要基于MATLAB平台,实现信号分析与滤波器的设计。并改变滤波器参数或特性低通、高通、带通或带阻,实现不同的滤波要求。设计产生一个连续信号,包含低频,中频,高频分量,对其进行采样,进行频谱分析,使用矩形窗设计不同特性的数字滤波器对信号进行滤波处理,分析所设计滤波器画出了频率特性曲线,并对信号进行滤波,画出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化,分析的结果表明初步实现了设计目标。MATLAB结合矩形窗设计FIR滤波器2第2章FIR滤波器设计原理21FIR滤波器FIR滤波器的数学表达式为(21)10NIYNHXNI式中N为FIR滤波器的抽头数;XN为第N时刻的输入样本;HI为FIR滤波器第I级抽头系数。普通的直接型FIR滤波器结构如图21所示。XN1/Z1/Z1/ZH0H1HN1HNYNH2图21卷积码编码器的一般形式在自适应处理、数据通信等领域中往往要求信号在传输过程中不能有明显的相位失真,FIR滤波器可以做到线性相位满足此要求。FIR滤波器实质上是一个分节的延迟线,把每一节的输出加权累加,得到滤波器的输出。对于FIR滤波器的单位脉冲响应HI只要满足以下2个条件之一,则为线性相位滤波器。221,ODSYMETRDUALHNII线性相位的FIR滤波器具有中心对称的特性,其对称中心在N/2处。FIR滤波器的结构主要是非递归结构,没有输出到输入的反馈。并且FIR滤波器很容易获得严格的线性相位特性,避免被处理信号产生相位失真【2】。而线性相位体现在时域中仅仅是HN在时间的延时,这个特点在图像信号处理、数据传输等波形传递系统中是非常重要的。此外,他不会发生阻塞现象,能避免强信号淹没弱信号,因此特别适合信号强弱相差悬殊的情况。其主要的不足之处是,其较好的性能是以较高的阶数为代价换来的。因此,在保证相同性能的前提下,努力降低其阶数是FIR数字滤波器设计的重要因素之一。MATLAB结合矩形窗设计FIR滤波器322窗函数FIR滤波器的设计方法有窗函数法、频率取样法和最优化设计法。其中窗函数法是设计FIR滤波器最简单有效的方法,也是最常用的方法。在本设计中,所用到的滤波器的系数都是借助于窗函数法完成的。窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。利用加窗函数进行截断和平滑,实现一个物理可实现且具有线性相位的FIR滤波器的设计目的。FIR滤波器的窗函数法设计过程为(23)XKIDTFDTFIJDDHEHHHE式中HDEJΩ为逼近的理想滤波器频率响应;HDK为理想滤波器的单位脉冲响应,是无限长序列【3】。为获取实际应用的FIR滤波,需将HDK截断,用有限长的HK近似表示,用窗函数可以避免吉布斯现象HKHDKΩK,最后得到实际FIR滤波的频率响应HEJΩ。设计常用的窗函数有矩形窗、汉宁窗、海明窗、凯撒窗等。矩形窗是一种比较容易实现的窗,本设计选择矩形窗实现在MATLAB下,这些窗函数分别为1.矩形窗WBOXCARN,产生一个N点的矩形窗函数。2.三角窗WTRIANGN,产生一个N点的三角窗函数。3.汉明窗WHAMMINGN,产生一个N点的汉明窗函数。4.汉宁窗WHANNINGN,产生一个N点的汉宁窗函数。4.布莱克曼窗WBLACKMANN,产生一个N点的布莱克曼窗函数。6.凯泽窗WKAISERN,BETA,产生一个N点的凯泽窗数,其中BETA为影响窗函数旁瓣的参数。23矩形窗这是一种最简单的窗函数,从阻带衰减的观点来看也是性能最差的一种,其窗函数为(24)1,01NMWNELS它的频率响应函数是MATLAB结合矩形窗设计FIR滤波器4(212SINS