FIR滤波器加窗效应分析课程设计.docVIP

引 言））同用MATIAB来仿真实现。一、摘要 5 二、FIR 数字滤波器的基本原理 6 2.1关于FIR滤波器 6 2.2 FIR滤波器的优点 6 2.3数字滤波器的设计 7 三、FIR数字滤波器设计的基本方法 8 3.1窗函数法 8 3.2、矩形窗 10 3.4海明窗函数 11 3.5布莱克曼窗 11 3.6三角窗函数 12 3.7凯泽窗 13 3.8多尔夫-切比雪夫窗 13 四、MATLAB仿真滤波实现 14 4.1MATLAB软件简介 14 4.1.1、友好的工作平台和编程环境 14 4.1.2简单易用的程序语言 15 4.1.3强大的科学计算机数据处理能力 15 4.1.4出色的图形处理功能 15 4.1.5应用广泛的模块集合工具箱 15 4.1.6实用的程序接口和发布平台 15 4.2 实验结果分析 16 4.3设计主要用到的MATLAB 函数 20 4.3.1 firl 函数 20 4.3.2 freqz 函数 20 4.3.3 plot函数 20 4.3.4 title函数 21 4.3.5 legend函数 21 4.3.6 xlabel（） ylabel（）函数 21 4.3.7 hold on 21 4.3.8 axis 21 五、课程设计心得 22 六、参考文献 23 一、摘要 数字滤波技术是数字信号处理的一个重要组成部分，滤波器的设计是信号处理的核心问题之一。FIR滤波器典型设计方法是借助于窗函数法和频率采样法。本次课程设计为分析FIR数字滤波器的基本原理，在MATLAB环境下利用窗函数设计FIR滤波器，实现FIR低通滤波器的设计仿真，并对常用的几种窗函数进行比较，同时给出在MATLAB环境下，用窗函数法设计FIR滤波器的过程和实例。仿真结果表明，设计的FIR滤波器的各项性能指标均可达到指定要求，且设计方法简单而有效，能得到较好的幅度特性。因此该方法为快速、高效地设计FIR滤波器提供了一个可靠而有效的途径。通过实验结果表明FIR滤波器准确度高、稳定性好，可以有效的滤除干扰信号，设计结果满足性能指标要求。数字滤波器的应用十分广泛，运行MATLAB语言，能很容易地设计出具有严格要求（如线性相位等）的滤波器。 关键字： FIR 滤波器T 错误！未找到引用源。为左下角图。根据卷积定理，即得实际滤波器的频率响应图形为图中右图。 由图可看出，加矩形窗后使实际频率响应偏离理想频率响应，主要影响有三个方面： （1）理想幅频特性陡直边缘处形成过渡带，过渡带宽取决于矩形窗函数频率响应的主瓣宽度； （2）过渡带两侧形成肩峰和波纹，这是矩形窗函数频率响应的旁瓣引起的，旁瓣相对值越大，旁瓣越多，波纹越多； （3）随窗函数宽度M的增大，矩形窗函数频率响应的主瓣宽度减小，但不改变旁瓣的相对值。 为了改善滤波器的性能，需使窗函数谱满足：主瓣尽可能窄，以使设计出来的滤波器有较陡的过渡带；第一副瓣面积相对主瓣面积尽可能小，即能量尽可能集中在主瓣，外泄少，使设计出来的滤波器的肩峰和余振小逼近于理想滤波器。但是这两个条件是相互矛盾的，实际应用中，折衷处理，兼顾各项指标。 上边只考虑了矩形窗，如果我们使窗的主瓣宽度尽可能地窄，旁瓣尽可能地小，可以获得性能更好的滤波器，通过改变窗的形状来达到这个目的。在数字信号处理的发展过程中形成了不同于矩形窗的很多窗函数，这些窗函数在主瓣和旁瓣特性方面各有特点，可满足不同的要求。为此，用窗函数法设计FIR数字滤波器时，要根据给定的滤波器性能指标选择窗口宽度N和窗函数w(n)。下面具体介绍几类窗函数及其特性。 3.2、矩形窗 矩形窗函数的时域形式可以表示为： 它的频域特性为： 3.3汉宁窗函数 汉宁窗函数的时域形式可以表示为： 它的频域特性为： 其中，为矩形窗函数的幅度频率特性函数。 汉宁窗函数的最大旁瓣值比主瓣值低31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍，为8π/N。 3.4海明窗函数 海明窗函数的时域形式可以表示为: 它的频域特性为: 其中，为矩形窗函数的幅度频率特性函数。 海明窗函数的最大旁瓣值比主瓣值低41dB，但它和汉宁窗函数的主瓣宽度是一样大的。 3.5布莱克曼窗 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，增加N可减少过渡带。 频谱的幅度函数为： +0.04 3.6三角窗函数 三角窗是最简单的频谱函数为非负的一种窗函数。三角窗函数的时域形式可以表示为： 当n为奇数时： 当n为偶数时： 它的频域特性为： 三角窗函数的主瓣宽度为8π/N，比矩形窗函数的主瓣宽度增加了一倍，但是它的旁瓣宽度却小得多。 3.7凯泽窗 以上几种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯泽窗则可自由选择主瓣宽度和旁瓣衰减,如图