基于matlab的fir滤波器设计与仿真 毕业设计论文

郑州轻工业学院课程设计说明书题目：基于MATLAB的FIR滤波器设计院（系）：电气信息工程学院专业班级：电子信息工程10-01郑州轻工业学院课程设计任务书题目基于MATLAB的FIR数字滤波器的设计专业、班级电子信息工程10级1班学号541001030136 姓名。

主要内容、基本要求、主要参考资料等：主要内容：采用窗函数法，分别利用五种窗函数（矩形窗、汉明窗、汉宁窗、三角形窗、布拉克曼窗），进行FIR数字滤波器的设计。

基本要求：根据给定的各类滤波器的技术指标，分别设计实现数字高通滤波器、数字带通滤波器和数字带阻滤波器，并据此进行分析总结：1、在相同的技术指标要求下，用不同的窗函数实现有何异同。

2、在相同的技术指标要求下，用FIR和IIR数字滤波器实现有何异同。

主要参考资料：1、《数字信号处理教程（第三版）》，程佩青著，清华大学出版社，2007。

2、《数字信号处理教程——MATLAB释义与实现（第2版）》，陈怀琛著，电子工业出版社，2008。

完成期限： 2013.6.24~2013.6.28指导教师签名：课程负责人签名：目录摘要 (I)1概述 (1)1.1课题背景及意义 (1)1.2 数字滤波器概述 (1)1.2.1 数字滤波器简介 (1)1.2.2 数字滤波器分类 (2)1.2.3 FIR滤波器概述 (2)2 FIR数字滤波器设计方法研究 (4)2.1 FIR滤波器的设计原理 (4)2.2 几种常见的窗函数 (5)2.2.1 矩形窗 (5)2.2.2 汉宁窗 (6)2.2.3 海明窗 (6)2.2.4 三角形窗函数 (6)2.2.5 布拉克曼窗 (7)2.3 五种窗函数基本参数比较 (8)2.4 窗函数设计FIR过程 (9)3 FIR滤波器设计 (10)3.1 要用到的两个主要MATLAB函数 (10)3.1.1 fir1函数介绍 (10)3.1.2 freqz函数介绍 (10)3.2 高通FIR滤波器设计 (10)3.2.1 基于矩形窗高通FIR滤波器 (11)3.2.2 基于海明窗高通FIR滤波器 (12)3.2.3 基于汉宁窗高通FIR滤波器 (13)3.2.4 基于三角形窗高通FIR滤波器 (14)3.2.5 基于布拉克曼窗高通FIR滤波器 (16)3.3 带通FIR滤波器设计 (17)3.3.1 基于矩形窗带通FIR滤波器 (17)3.3.2 基于海明窗带通FIR滤波器 (18)3.3.3 基于汉宁窗带通FIR滤波器 (20)3.3.4 基于三角形窗带通FIR滤波器 (21)3.3.5基于布拉克曼窗带通FIR滤波器 (22)3.4 带阻FIR滤波器设计 (23)3.4.1 基于矩形窗带阻FIR滤波器 (24)3.4.2 基于海明窗带阻FIR滤波器 (25)3.4.3 基于汉宁窗带阻FIR滤波器 (26)3.4.4 基于三角形窗带阻FIR滤波器 (27)3.4.5 基于布拉克曼窗带阻FIR滤波器 (29)4 结论 (31)5 总结 (32)致谢 (33)参考文献 (34)基于MATLAB 的FIR数字滤波器设计摘要滤波器的设计是信号处理的核心问题之一。

本设计根据FIR滤波器的原理，简要介绍了利用各种窗函数设计FIR数字滤波器的方法，并且分析了不同窗函数形式对滤波器性能的影响。

采用Matlab软件对FIR数字滤波器进行仿真设计，简化了设计中繁琐的计算。

本文采用窗函数法通过调用Matlab函数设计FIR数字滤波器。

分析表明所用的矩形窗，海明窗，汉宁窗，三角窗，布拉克曼窗，五种窗函数中布拉克曼窗性能最好。

最后运用matlab进行了仿真。

关键词：matlab 数字信号处理 FIR滤波器1概述1.1课题背景及意义随着数字化时代的到来, 数字信号处理已经成为当今一门极其重要的技术, 并且在通信、语音、图像、自动控制等众多领域得到了广泛的应用。

而滤波器在实际的信号处理中起到了重要的作用，它是去除信号中噪声的基本手段，因此在数字信号的处理中滤波是一个基本问题。

FIR滤波器是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限的，其在通信、图像处理、模式识别等领域都有着广泛的应用。

几乎在所有的工程技术领域中都会涉及到信号的处理问题，其信号表现形式有电、磁、机械以及热、光、声等。

信号处理的目的一般是对信号进行分析、变换、综合、估值与识别等。

如何在较强的噪声背景下提取出真正的信号或信号的特征，并将其应用于工程实际是信号处理的首要任务。

数字信号处理中一个非常重要且应用普遍的技术就是数字滤波。

数字滤波器有FIR数字滤波器和IIR数字滤波器，IIR数字滤波器的设计方法是利用模拟滤波器成熟的理论及设计图表进行设计的，因而保留了一些典型模拟滤波器优良的幅度特性，但设计中只考虑了幅度特性，没考虑相位特性，所设计的滤波器一般是某种确定的非线性相位特性。

为了得到线性相位特性，对IIR 滤波器必须另外加相位校正网络，使滤波器设计变得复杂，成本也高，又难以得到严格的线性相位特性。

而FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性，同时为了使FIR数字滤波器的设计更优化，因而研究FIR 数字滤波器的优化设计具有重要的理论意义。

MATLAB 是一种面向科学和工程计算的语言, 它集数值分析、矩阵运算、信号处理和图形显示于一体,功能强大、使用方便。

相比传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，利用matlab信号处理工具箱（signal processing toolbox）可以快速有效的设计由软件组成的常规数字滤波器的设计方法，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量。

1.2 数字滤波器概述1.2.1 数字滤波器简介数字滤波在DSP中占有重要地位。

数字滤波器按实现的网络结构或者从单位脉冲响应，分为IIR（无限脉冲响应）和FIR（有限脉冲响应）滤波器。

如果IRR滤波器和FIR滤波器具有相同的性能，那么通常IIR滤波器可以用较低的阶数获得高的选择性，执行速度更快，所有的存储单元更少，所以既经济又高效。

数字滤波器精确度高，使用灵活，可靠性高，具有模拟设备没有的许多优点，已广泛地应用与各个科学技术领域，例如数字电视，语音，通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。

随着信息时代数字时代的到来，数字滤波技术已经成为一门及其重要的科学和技术领域。

以往的滤波器大多采用模拟电路技术，但是模拟电路技术存在很多难以解决的问题，而采用数字则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前的发展方向。

数字滤波在通信。

图像编码，语音编码，雷达等许多领域有着十分广泛的应用。

目前，数字信号滤波器的设计在图像处理，数据压缩等方面的应用取得了令人瞩目的进展和成就。

近年来迅速发展。

1.2.2 数字滤波器分类按功能分：低通、高通、带通、带阻、全通滤波器。

按实现的网络结构或单位抽样响应分：无限脉冲响应滤波器（IIR滤波器）、有限脉冲响应滤波器（FIR滤波器），其中 FIR数字滤波器网络不存在反馈支路，其单位冲激响应为有限长。

IIR数字滤波器网络存在反馈支路，即信号流图中存在环路，其单位冲激响应为无限长。

另外，按线性系统它可以被分为线性与非线性，按因果性分因果与非因果等。

其中，线性时不变的数字滤波器是最基本的类型；而由于数字系统可以对延时器加以利用，因此可以引入一定程度的非因果性，获得比传统的因果滤波器更灵活强大的特性；相对于IIR滤波器，FIR滤波器有着易于实现和系统绝对稳定的优势，因此得到广泛的应用。

1.2.3 FIR滤波器概述根据冲激响应的时域特性, 数字滤波器可分为无限长冲激响应滤波器(IIR) 和有限长冲激响应滤波器(FIR) 。

相比于IIR滤波器，FIR的突出优点是: （1）具有精确的线性相位（2）总是稳定的（3）硬件容易实现（4）滤波器的过渡过程具有有限区间但其在满足同样阻带衰减的情况下需要的阶数较高。

滤波器的阶数越高, 占用的运算时间越多, 因此在满足指标要求的情况下应尽量减少滤波器的阶数。

FIR 滤波器的基本结构可以理解为一个分节的延时线, 把每一节的输出加权累加, 可得到滤波器的输出。

FIR 滤波器的冲激响应h (n) 是有限长的, 数学上M 阶FIR滤波器可以表示为:10()()()N i y n h i x n i -==-∑ （1-1）其系统函数为:（1-2）由于FIR 滤波器的单位冲激响应是一个有限长序列，系统函数是的次多项式，它在平面有个零点，同时在原点有阶重极点。

因此，永远稳定。

2 FIR 数字滤波器设计方法研究2.1 FIR 滤波器的设计原理FIR 滤波器的系统输入输出差分方程为：y[n]=∑=-1-0k )()(N k n x k h (2-1)所以FIR 滤波器的系统函数为：(2-2)由于FIR 滤波器的单位脉冲响应是一个有限长序列，是的次多项式，它在Z 平面上有个零点，同时在原点有阶重极点。

因此，永远稳定。

FIR 滤波器设计的任务是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。

由于FIR 滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。

FIR 滤波器目前常用的设计方法有窗函数法和频率采样法，窗函数法是从时域进行设计，而频率采样法是从频域进行设计。

窗函数法由于简单、物理意义清晰，因而得到了较为广泛的应用。

窗函数法设计的基本思想是：首先根据技术指标要求，选取合适的阶数 N 和窗函数的类型，使其幅频特性逼近理想滤波器幅频特性。

其次，因为理想滤波器的是无限长的，所以需要对进行截断，数学上称这种方法为窗函数法。

(2-3)(2-4)(2-5) 其中是希望逼近的滤波器的单位脉冲响应，是希望逼近的滤波器的幅频特性，是所设计的滤波器的单位脉冲响应，是一个矩阵序列，长度为N 。

在截断的过程中要保证FIR 滤波器的稳定性与线性相位的特性。

由于截断效应使所设计滤波器的幅频特性存在误差，称之为吉布斯效应。

最后，验证所设计的滤波器是否满足要求，若不满足，重新设计。

常用的窗函数有矩形窗、三角窗、汉宁窗和凯塞窗等。

窗函数的选择一般要满足以下三个条件：(1)具有较低的旁瓣幅度，尤其是第一旁瓣幅度。

(2)旁瓣幅度下降速度要快，以利于增加阻带衰减。

(3)主瓣宽度要窄，以获得较陡的过渡带。

2.2 几种常见的窗函数2.2.1 矩形窗矩形窗函数的时域形式可以表示为：1,01()()0,N n N w n R n ≤≤-⎧==⎨⎩其他 (2-6)它的频域特性为：()1j j 2sin 2e esin 2N N W ωωωω-⎛⎫- ⎪⎝⎭⎛⎫ ⎪⎝⎭=⎛⎫ ⎪⎝⎭ (2-7)矩形窗使用最多，习惯上不加窗就是使信号通过了矩形窗。