基于Matlab的FIR低通滤波器设计[开题报告]2011-11-07

电子信息学院

本科毕业设计(论文)开题报告

论文题目基于Matlab的FIR低通滤波器设计

学生姓名专业班级性别

一、课题研究意义及现状

所谓数字信号处理技术是指运用数值计算的方法对信号进行分析、变换、综合、估值与识别等的处理。随着电子技术和电子计算机技术的发展，数字信号处理技术受到了越来越广泛的关注，数字信号处理的理论和技术也在不断丰富和完善，新的理论和新技术层出不穷。可以说，数字信号处理是发展最快、应用最广泛、成效最显著的新科学之一，目前已广泛地应用在语音、雷达、声纳、地震、图像、通信、控制、生物医学、遥感遥测、地质勘探、航空航天、故障检测、自动化仪表等领域。

数字滤波器是数字信号处理中的重要组成部分之一，在数字信号处理中，数字滤波占有极其重要的地位。滤波器可以分为IIR滤波器和FIR滤波器两大类，IIR滤波器并不能得到严格的线性相位特性，因此在许多实际应用中为了得到线性相位特征，还必须另外增加相位校正网络，使滤波器设计变得复杂，成本也高，而FIR滤波器在保证幅度特性满足技术要求的同时，还可以很容易的做到有严格的线性相位特性。它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。

目前对数字滤波器的设计有多种方法，其中著名的Matlab软件包，功能强大、使用方便。传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。利用Matlab信号处理工具箱可以快速有效的设计由软件组成的常规数字滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。

本课题在基于Matlab的FIR低通滤波器的设计方法、编程以及实际应用上有着积极的意义。

二、课题研究的主要内容和预期目标

主要研究基于Matlab的FIR低通滤波器。具体是采用窗函数设计法、频率采样法、等波纹逼近法进行设计，并用Matlab软件编写程序进行仿真。

1．研究FIR滤波器的定义、分类、应用以及设计方法。

2．了解FIR滤波器窗函数的设计方法及原理。

3．了解FIR滤波器频率采样的设计方法及原理。

4．了解FIR滤波器等波纹逼近的设计方法及原理。

5．确定滤波器的技术指标，并进行Matlab仿真。

结合上述内容，设计编程，完成利用Matlab软件设计FIR低通滤波器并进行仿真优化。通过设计编程，加深对知识的理解和运用。

三、课题研究的方法及措施

本课题主要研究的是基于Matlab的FIR低通滤波器设计，因此在研究过程中主要是采用理论分析结合软件仿真来实现。总的研究方法是在课题研究之前，掌握一定的理论基础，了解FIR滤波器是如何进行信号各个波段的滤波的，在理论中寻找适合本课题的设计方案；然后，了解Matlab软件应用，使用Matlab软件设计数字低通滤波器并进行仿真优化，实现数字低通滤波器的滤波。具体措施如下：

1. 首先阅读广泛资料，对FIR滤波器的定义、设计方法做一下介绍，分析FIR低通滤

波器的研究意义及作用。

2．Matlab作为工作平台和开发工具，熟悉Matlab软件程序设计的基本原理和实验环境，实现仿真。

3．针对FIR低通滤波器的滤波基本原理进行研究，查阅相关资料书籍。并利用窗函数的设计方法进行滤波设计，并运用Matlab进行模拟仿真。

4．查漏补缺，进行论文的改进。

四、课题研究进度计划

2011/2012（第一学期）：

明确任务，收集资料，确定系统总体方案设计，完成开题报告、文献综述、外文翻译以及开题答辩。

2011/2012（第二学期）：

第1周至第3周：完成理论基础的研究，重点学习基于Matlab的FIR低通滤波器的设计方法，特别是窗函数法的运用。

第3周至第4周：掌握Matlab程序开发方法、调试工具的使用及优化级别的选择等，实现数字滤波器的设计与仿真。编写基于Matlab的FIR低通滤波算法并进行

调试仿真和反复论证。

第5周至第6周：完善设计作品，论文修改。进行方案的补充和修改；撰写毕业论文，修改完善毕业论文。

第7周至第8周：做好论文答辩的PPT资料，准备答辩，并提交所有电子文档材料。

五、参考文献

[1] 丁玉美, 高西全.数字信号处理[M]. 第三版. 西安: 西安电子科技大学出版社, 2008, 6.

[2] 闫胜利. FIR滤波器及设计原理[J]. 长春工程学院学报(自然科学版), 2003, 6, 4(1): 21-24.

[3] 姚齐国. 基于MATLAB的数字滤波器的设计[J]. 江西理工大学学报, 2006, 2, 27(1):

50-52.

[4] 杨守卫. FIR数字滤波器应用分析探讨[J]. 河北省工程咨询院学报, 2011, 7, 29(15):

47-49.

[5] 朱敏. MA TLAB数字信号处理工具箱的开发和应用[J]. 信息与电脑, 2010, 2, 26(8):

154-155.

[6] 姚海燕. FIR数字滤波器设计窗函数法与频率抽样法比较[J].安阳工学院学报, 2007, 6,

12(6): 51-53.

[7] 刘波. MA TLAB信号处理[M]. 北京: 电子工业出版社, 2006, 7.

[8] 邬国扬, 张萱, 吕昂等. 模拟电子技术[M]. 第二版. 西安：西安电子科技大学出版社,

2007, 3.

[9] 郑佳春. 数字滤波器的Matlab设计与DSP的实现[J]. 单片机与嵌入式系统应用, 2003

（7）: 28-28．

[10] 赵惠昌, 王丹琦. 一种基于MATLAB及FPGA的FIR低通滤波器的设计与实现[J]. 电子

设计工程, 2011, 7, 19(14): 131-134.

[11] 张小虹, 黄忠虎, 邱正伦等. 数字信号处理[M].北京：机械工业出版社, 2008, 9.

[12] 姚剑清. 简明数字信号处理[M]. 北京: 人民邮电出版社, 2009, 1.

[13] J.G. Proakis , D.G. Manolakis. Digital Signal Processing-Principles, Algorithms and

Applicat- ions[M]. New Delhi: Prentice-Hall, 2000, 4.

[14]Uwe Meyer-Baese. Digital Signal Processing with Field Programmable Gate Arrays

[J]．Signal And Communication Technology, 2007, 10, 1: 165-165.

指导教师意见

指导教师签名：

年月日

fir低通滤波器matlab编程滤波前后图形

Matlab实现振动信号低通滤波附件txt中的数字是一个实测振动信号，采样频率为5000Hz，试设计一个长度为M=32的FIR低通滤波器，截止频率为600Hz，用此滤波器对此信号进行滤波。要求：（1）计算数字截止频率；（2）给出滤波器系数；（3）绘出原信号波形；（4）绘出滤波后的信号波形；解答过程：第一部分：数字截止频率的计算数字截止频率等于截止频率除以采样频率的一半，即 n=600/5000/2=0.24第二部分：滤波器系数的确定在matlab中输入如下程序，即可得到滤波器系数： n=32 Wn=0.24 b=fir1(n,Wn) 得到的滤波器系数b为 Columns 1 through 9 -0.0008-0.0018-0.0024-0.00140.00210.00750.01100.0077-0.0054Columns 10 through 18

-0.0242-0.0374-0.02990.00870.07560.15370.21660.24070.2166Columns 19 through 27 0.15370.07560.0087-0.0299-0.0374-0.0242-0.00540.00770.0110Columns 28 through 33 0.00750.0021-0.0014-0.0024-0.0018-0.0008 第三部分：原信号波形将附件4中的dat文件利用识别软件读取其中的数据，共1024个点，存在TXT文档中，取名bv.txt，并复制到matlab的work文件夹。在matlab中编写如下程序： x0=load('zhendong.txt');%找到信号数据地址并加载数据。 t=0:1/5000:1023/5000;%将数据的1024个点对应时间加载 figure(1); plot(t,x0); xlabel('t/s'); ylabel('幅值'); 运行之后就得到如下波形，即振动信号的原始波形图： 1.5 1 0.5 幅值

fir低通滤波器设计(完整版)

电子科技大学信息与软件工程学院学院标准实验报告（实验）课程名称数字信号处理电子科技大学教务处制表

电子科技大学实验报告学生姓名：学号：指导教师：实验地点：实验时间：14-18 一、实验室名称：计算机学院机房二、实验项目名称：fir 低通滤波器的设计三、实验学时：四、实验原理： 1. FIR 滤波器 FIR 滤波器是指在有限范围内系统的单位脉冲响应h[k]仅有非零值的滤波器。M 阶FIR 滤波器的系统函数H(z)为 ()[]M k k H z h k z -==∑ 其中H(z)是k z -的M 阶多项式，在有限的z 平面内H(z)有M 个零点，在z 平面原点z=0有M 个极点. FIR 滤波器的频率响应 ()j H e Ω 为 0 ()[]M j jk k H e h k e Ω -Ω ==∑ 它的另外一种表示方法为 () ()()j j j H e H e e φΩΩΩ=

其中 () j H e Ω和()φΩ分别为系统的幅度响应和相位响应。若系统的相位响应()φΩ满足下面的条件 ()φαΩ=-Ω 即系统的群延迟是一个与Ω没有关系的常数α，称为系统H(z)具有严格线性相位。由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR 滤波器设计中一般使用广义线性相位。如果一个离散系统的频率响应 ()j H e Ω 可以表示为 ()()()j j H e A e αβΩ-Ω+=Ω 其中α和β是与Ω无关联的常数，()A Ω是可正可负的实函数，则称系统是广义线性相位的。如果M 阶FIR 滤波器的单位脉冲响应h[k]是实数，则可以证明系统是线性相位的充要条件为 [][]h k h M k =±- 当h[k]满足h[k]=h[M-k],称h[k]偶对称。当h[k]满足h[k]=-h[M-k],称h[k]奇对称。按阶数h[k]又可分为M 奇数和M 偶数，所以线性相位的FIR 滤波器可以有四种类型。 2. 窗函数法设计FIR 滤波器窗函数设计法又称为傅里叶级数法。这种方法首先给出()j d H e Ω， ()j d H e Ω 表示要逼近的理想滤波器的频率响应，则由IDTFT 可得出滤波器的单位脉冲响应为 1 []()2j jk d d h k H e e d π π π ΩΩ-= Ω ? 由于是理想滤波器，故 []d h k 是无限长序列。但是我们所要设计的FIR 滤波器，其h[k]是有限长的。为了能用FIR 滤波器近似理想滤波器，需将理想滤波器的无线长单位脉冲响应 []d h k 分别从左右进行截断。当截断后的单位脉冲响应 []d h k 不是因果系统的时候，可将其右移从而获得因果的FIR 滤波器。

基于MATLAB的FIR数字低通滤波器设计.doc

基于MATLAB的FIR数字低通滤波器设计 ..基于MATLAB的FIR数字低通滤波器设计作者：周龙刚（陕西理工学院物理与电信工程学院通信工程专业2011级4班,陕西汉中723003）指导老师：井敏英[摘要]FIR数字滤波器是数字信号处理的一个重要组成部分，由于FIR数字滤波器具有严格的线性相位，因此在信息的采集和处理过程中得到了广泛的应用。本文介绍了FIR数字滤波器的概念和线性相位的条件，分析了窗函数法、频率采样法和等波纹逼近法设计FIR滤波器的思路和流程。在分析三种设计方法原理的基础上，借助Matlab仿真软件工具箱中的fir1实现窗函数法中的哈明窗设计FIR低通滤波器。[关键词] FIR数字滤波器;线性相位窗函数;法哈明窗;MatlabDesigning FIR low-pass digital filter based on VHDL Zhoulonggang(Grade11,Class4,Major of Communication Engineering，School of Physics and Telecommunication Engineering , Shaanxi University of Technology,Hanzhong 723003,Shaanxi)Tutor:JingYingMinAbstract:FIR digital filter is an important part of digital signal processing,the FIR digitalfilter with linear phase, so it has been widely applied in the collection and processing of information in the course of. This paper introduces the concept of FIR digital filter with linear phase conditions, analysis of the window functionmethod and frequency sampling method and the ripple

fir低通滤波器设计报告

滤波器设计原理本文将介绍数字滤波器的设计基础及用窗函数法设计FIR 滤波器的方法，运用MATLAB 语言实现了低通滤波器的设计以及用CCS软件进行滤波效果的观察。读取语音文件，并加入一定的随机噪声，最后使用窗函数滤波法进行语音滤波，将加噪后的语音文件转换为.dat文件使其能和ccs软件链接，输出个阶段的时域和频域波形。根据数字滤波器冲激响应函数的时域特性。可将数字滤波器分为两种，即无限长冲激响应( IIR) 滤波器和有限长冲激响应(FIR) 滤波器。IIR 滤波器的特征是具有无限持续时间的冲激响应；FIR 滤波器冲激响应只能延续一定时间。其中FIR 滤波器很容易实现严格的线性相位，使信号经过处理后不产生相位失真，舍入误差小，稳定等优点。能够设计具有优良特性的多带通滤波器、微分器和希尔伯特变换器，所以在数字系统、多媒体系统中获得极其广泛的应用。FIR数字滤波器的设计方法有多种，如窗函数设计法、最优化设计和频率取样法等等。而随着MATLAB软件尤其是MATLAB 的信号处理工具箱和Simulink 仿真工具的不断完善，不仅数字滤波器的计算机辅助设计有了可能而且还可以使设计达到最优化。 FIR滤波器的窗函数法的设计采用汉明窗设计低通FIR滤波器使用b=fir1(n,Wn)可得到低通滤波器。其中，0Wn1，Wn=1相当于0.5。其语法格式为 b=fir1(n,Wn); 采用：b=fir1(25, 0.25); 得到归一化系数：

或者在命令行输入fdatool进入滤波器的图形设置界面，如下图所示得到系数（并没有归一化） const int BL = 26; const int16_T B[26] = { -26, 33, 126, 207, 138, -212, -757, -1096, -652, 950, 3513, 6212, 7948, 7948, 6212, 3513, 950, -652, -1096, -757, -212, 138, 207, 126, 33, -26 }; FIR滤波器的设计（Matlab）技术指标为：采用25阶低通滤波器，汉明窗(Hamming Window)函数，截止频率为1000Hz，采样频率为8000Hz，增益40db。下面的程序功能是：读取语音文件，并加入一定的随机噪声，最后使用窗函数滤波法进行语音滤波，将加噪后的语音文件转换为.dat文件使其能和ccs软件链接，输出个阶段的时域和频域波形。

FIR数字滤波器设计与使用

实验报告课程名称：数字信号处理指导老师：刘英成绩：_________________实验名称： FIR数字滤波器设计与使用同组学生姓名：__________ 一、实验目的和要求设计和应用FIR低通滤波器。掌握FIR数字滤波器的窗函数设计法，了解设计参数（窗型、窗长）的影响。二、实验内容和步骤编写MATLAB程序，完成以下工作。 2-1 设计两个FIR低通滤波器，截止频率 C =0.5。（1）用矩形窗，窗长N=41。得出第一个滤波器的单位抽样响应序列h 1(n)。记下h 1 (n) 的各个抽样值，显示h 1 (n)的图形（用stem(.)）。求出该滤波器的频率响应（的N 个抽样）H 1(k)，显示|H 1 (k)|的图形（用plot(.)）。（2）用汉明窗，窗长N=41。得出第二个滤波器的单位抽样响应序列h 2(n)。记下h 2 (n) 的各个抽样值，显示h 2(n)的图形。求出滤波器的频率响应H 2 (k)，显示|H 2 (k)|的图形。（3）由图形，比较h 1(n)与h 2 (n)的差异，|H 1 (k)|与|H 2 (k)|的差异。 2-2 产生长度为200点、均值为零的随机信号序列x(n)（用rand(1,200)0.5）。显示x(n)。求出并显示其幅度谱|X(k)|，观察特征。 2-3 滤波（1）将x(n)作为输入，经过第一个滤波器后的输出序列记为y 1(n)，其幅度谱记为|Y 1 (k)|。显示|X(k)|与|Y 1 (k)|，讨论滤波前后信号的频谱特征。（2）将x(n)作为输入，经过第二个滤波器后的输出序列记为y 2(n)，其幅度谱记为|Y 2 (k)|。比较|Y 1(k)|与|Y 2 (k)|的图形，讨论不同的窗函数设计出的滤波器的滤波效果。 2-4 设计第三个FIR低通滤波器，截止频率 C =0.5。用矩形窗，窗长N=127。用它对x(n)进行滤波。显示输出信号y

基于矩形窗设计FIR数字滤波器

郑州航空工业管理学院《电子信息系统仿真》课程设计 14级专业班级题目基于矩形窗设计FIR数字滤波器姓名学号二О一六年十一月二十五日

第一章 FIR 滤波器的设计原理及方法 FIR 滤波器通常采用窗函数方法来设计。窗设计的基本思想是，首先选择一个适当的理想选频滤波器（它总是具有一个非因果，无限持续时间脉冲响应），然后街区（加窗）它的脉冲响应得到线性相位和因果FIR 滤波器。我们用Hd （e^jw ）表示理想的选频滤波器，它在通带上具有单位增益和线性相位，在阻带上具有零响应。一个带宽wc

matlab的fir高通数字滤波器的设计和分析

摘要无限长脉冲数字滤波器的设计方法只考虑了幅度特性，没有考虑相位特性，所设的滤波器一般是某种确定的非线性相位特性。有限脉冲响应（FIR）滤波器在保证了幅度特性满足技术要求的同时，很容易做到有严格的线性相位特性。本课题利用MATLAB软件实现。MATLAB是“矩阵实验室”（MATrix LABoratoy）的缩写，是一种科学计算软件，它使用方便，输入简捷，运算高效，内容丰富，因此利用MATLAB软件，通过一系列较为系统的函数法，根据已知的技术指标，就可以设计出满足要求的滤波器。关键字:MATLAB；窗函数；FIR带阻数字滤波器；线性相位

目录 1.FIR滤波器简介 (3) 1.1 FIR的特点 (3) 2.2线性相位 (3) 2．主要设计内容 (5) 3.窗函数 (6) 3.1常用窗函数 (6) 3.2窗函数的指标 (9) 4应用窗函数法设计 FIR 数字滤波器的步骤 (10) 4.1数字高通滤波器的设计： (10) 总结 (11) 参考文献 (12) 附录 (13)

1.FIR 滤波器简介数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR ）滤波器和有限冲激响应（FIR ）滤波器。 1.1 FIR 的特点 FIR 滤波器的主要优点为：系统总是稳定的，FIR 滤波器的系统函数可以表示为 (2-1) 易知，H (z ) 在 Z 平面上有 N -1个零点，z =0 是 N -1 阶极点，因此FIR 系统总是稳定的（极点都在单位圆内）。FIR 滤波器的优点之二：容易实现线性相位。当 FIR 系统的单位冲激响应满足时，该系统具有线性相位。（N 为奇数）（2-2）（N 为偶数）（2-3） FIR 滤波器的优点之三：允许设置多通带（或多阻带）滤波器。FIR 滤波器的优点之四：FIR 滤波器可以采用 FFT 方法实现其功能，从而大大提高效率。FIR 滤波器的缺点：由于 FIR 系统只有零点，因此这类系统不像FIR 滤波器不像 IIR 滤波器那样容易取得比较好的通带与阻带衰减特性。要取得较好的衰减特性，一般要求 H (z ) 的阶次较高。综合起来看， FIR 滤波器具有IIR 滤波器没有的许多特点，得到了越来越广泛的应用。 FIR 滤波器的设计方法主要有三种：a.窗函数设计法；b.频率抽样发；c.最小平法抽样法；这里我主要讨论在MATLAB 环境下通过调用信号分析与处理工具箱的几类窗函数来设计滤波器并分析与比较其性能 2.2线性相位一个单一频率的正弦信号通过一个系统，假设它通过这个系统的时间需要t ，则这个信号的输出相位落后原来信号wt 的相位。从这边可以看出，一个正弦信号通过一个系统落后的相位等于它的w *t ；反过来说，如果一个频率为w 的正弦信号通过系统后，它的相位落后delta ，则该信号被延迟了delta /w 的时间。在实 11 1) 1(10)()()()(--=-----=-===∑∑N N n n N N N n n z z f z n h z z n h z H )1()(n N h n h --±=2/)1()(--=N ωω?2/)1(2/)(--=N ωπω?

基于matlab的FIR低通高通带通带阻滤波器设计

基于matlab的FIR低通-高通-带通-带阻滤波器设计

————————————————————————————————作者：————————————————————————————————日期：

北京师范大学课程设计报告课程名称： DSP 设计名称：FIR 低通、高通带通和带阻数字滤波器的设计姓名：学号: 班级：指导教师：起止日期：课程设计任务书

学生班级：学生姓名：学号：设计名称： FIR 低通、高通带通和带阻数字滤波器的设计起止日期：指导教师：设计目标： 1、采用Kaiser 窗设计一个低通FIR 滤波器要求：采样频率为8kHz ；通带：0Hz~1kHz ，带内波动小于5%；阻带：1.5kHz ，带内最小衰减：Rs=40dB 。 2、采用hamming 窗设计一个高通FIR 滤波器要求：通带截至频率wp=rad π6.0, 阻带截止频率ws=rad π4.0，通带最大衰减dB p 25.0=α,阻带最小衰减dB s 50=α 3、采用hamming 设计一个带通滤波器低端阻带截止频率 wls = 0.2*pi ；低端通带截止频率 wlp = 0.35*pi ；高端通带截止频率 whp = 0.65*pi ；高端阻带截止频率 whs = 0.8*pi ； 4、采用Hamming 窗设计一个带阻FIR 滤波器要求：通带：0.35pi~0.65pi ，带内最小衰减Rs=50dB ；阻带：0~0.2pi 和0.8pi~pi ，带内最大衰减：Rp=1dB 。

FIR 低通、高通带通和带阻数字滤波器的设计一、设计目的和意义 1、熟练掌握使用窗函数的设计滤波器的方法，学会设计低通、带通、带阻滤波器。 2、通过对滤波器的设计，了解几种窗函数的性能，学会针对不同的指标选择不同的窗函数。二、设计原理一般，设计线性相位FIR 数字滤波器采用窗函数法或频率抽样法，本设计采用窗函数法，分别采用海明窗和凯泽窗设计带通、带阻和低通。如果所希望的滤波器的理想频率响应函数为)(jw d e H ，如理想的低通，由信号系统的知识知道，在时域系统的冲击响应h d (n)将是无限长的，如图2、图3所示。 H d (w) -w c w c 图2 图3 若时域响应是无限长的，则不可能实现，因此需要对其截断，即设计一个FIR 滤波器频率响应∑-=-=1 0)()(N n jwn jw e n h e H 来逼近)(jw d e H ，即用一个窗函数w(n)来截断h d (n)，如式3所示： )()()(n w n h n h d = （式1）。最简单的截断方法是矩形窗，实际操作中，直接取h d (n)的主要数据即可。 )(n h 作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数为： ∑-=-=1 0)()(N n jwn jw e n h e H （式2）令jw e z =，则 ∑-=-=1 0)()(N n n z n h z H （式3），式中，N 为所选窗函数)(n w 的长度。

FIR低通数字滤波器的设计要点

《DSP技术与应用》课程设计报告课题名称：基于DSP Builder的FIR数字滤波器的设计与实现学院：电子信息工程学院班级：11级电信本01班学号：姓名：

题目基于DSP Builder的FIR数字滤波器的设计与实现摘要 FIR数字滤波器是数字信号处理的一个重要组成部分，由于FIR数字滤波器具有严格的线性相位，因此在信息的采集和处理过程中得到了广泛的应用。本文介绍了FIR数字滤波器的概念和线性相位的条件，分析了窗函数法、频率采样法和等波纹逼近法设计FIR滤波器的思路和流程。在分析三种设计方法原理的基础上，借助Matlab仿真软件工具箱中的fir1、fir2和remez子函数分别实现窗函数法、频率采样法和等波纹逼近法设计FIR滤波器。然后检验滤波器的滤波效果，采用一段音频进行加噪声然后用滤波器滤，对比三段音频效果进而对滤波器的滤波效果进行检验。仿真结果表明，在相频特性上，三种方法设计的FIR滤波器在通带内都具有线性相位；在幅频特性上，相比窗函数法和频率采样法，等波纹逼近法设计FIR滤波器的边界频率精确，通带和阻带衰减控制。

Abstract FIR digital filter is an important part of digital signal processing, the FIR digital filter with linear phase, so it has been widely applied in the collection and processing of information in the course of. This paper introduces the concept of FIR digital filter with linear phase conditions, analysis of the window function method and frequency sampling method and the ripple approximation method of FIR filter design ideas and processes. Based on analyzing the principle of three kinds of design methods, by means of fir1, fir2 and Remez function of Matlab simulation software in the Toolbox window function method and frequency sampling method and respectively realize equiripple approximation method to design FIR filter. Then test the filtering effect of the filter, using an audio add noise and then filter, test three audio effects and comparison of filter filtering effect. Simulation results show that the phase frequency characteristic, three design methods of FIR filter with linear phase are in the pass band; the amplitude frequency characteristics, compared with the window function method and frequency sampling method, equiripple approximation method Design of FIR filter with accurate boundary frequency, the passband and stopband attenuation control.

基于matlab的数字滤波器的设计与仿真开题报告.pdf

一、选题背景、目的和意义背景信号与信号处理时信息科学中近几十年来发展最为迅速的学科之一。长期以来?信号处理技术一直用于转换、产生模拟或数字信号?其中最为频繁应用的领域就是信号的滤波。数字滤波是语音、图像处理、模式识别和谱分析等应用中的一个基本处理部件?它可以满足滤波器对幅度和相位特性的严格要求?避免模拟滤波器无法克服的电压漂移、温度漂移和噪声等问题。目的我的毕业设计的课题是《基于matlab的数字滤波器的设计与仿真》?其主要目的是通过此次课程设计进一步学习和巩固数字信号处理及其相关知识?并学会利用所学的知识能?在设计过程中能综合运用所学知识内容?进一步熟悉和掌握matlab的使用方法;对数字滤波器的原理有较深的了解;为即将进入社会参加工作打下坚实的基础;掌握收集资料、消化资料和综合资料的能力等等。意义从事电子通信业而不能熟练操作使用matlab电子线路设计软件?在工作和学习中将是寸步难行的。在数学、电子、金融等行业?使用matlab等计算机软件对产品进行设计、仿真在很早以前就已经成为了一种趋势?这类软件的问世也极大地提高了设计人员在通信、电子等行业的产品设计质量与效率。众所周知?实际过程中信号传输都要经过调制与解调这一过程?由于消息传过来的原始信号即调制信号具有频谱较低的频谱分量?这种信号在许多信道中不宜传输。因而?在通信系统的发送端通常需要有调制过程?反之在接收端则需要有解调过程。二、国内外的研究综述数字滤波在DSP?数字信号处理?中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应?分为IIR(无限脉冲响应?和FIR?有限脉冲响应?滤波器。如果IIR滤波器和FIR滤波器具有相同的性能?那么通常IIR滤波器可以用较低的阶数获得高的选择性?执行速度更快?所有的储存单元更少?所有既经济又高效。

用窗函数法设计FIR数字低通滤波器要点

河北科技大学课程设计报告学生姓名：学号：专业班级：课程名称：学年学期指导教师： 20 年月

课程设计成绩评定表学生姓名学号成绩专业班级起止时间设计题目指导教师评指导教师：语年月日

目录 1. 窗函数设计低通滤波器 1.1设计目的 (1) 1.2设计原理推导与计算 (1) 1.3设计内容与要求 (2) 1.4设计源程序与运行结果 (3) 1.5思考题 (10) 1.6心得体会 (14) 参考文献 (15)

1.窗函数设计低通滤波器 1.1设计目的 1. 熟悉设计线性相位数字滤波器的一般步骤。 2. 掌握用窗函数法设计FIR 数字滤波器的原理和方法。 3. 熟悉各种窗函数的作用以及各种窗函数对滤波器特性的影响。 4. 学会根据指标要求选择合适的窗函数。 1.2设计原理推导与计算如果所希望的滤波器的理想的频率响应函数为() ωj d e H ，则其对应的单位脉冲响应为 ()() ωπ ωωπ π d e e H n h j j d d ?- = 21 （4.1）窗函数设计法的基本原理是设计设计低通FIR 数字滤波器时，一般以理想低通滤波特性为逼近函数() ωj e H ，即 () ?????≤<≤=-π ωωωωωαω c c j j d ，， e e H 0,其中21-=N α ()() ()[]() a n a n d e e d e e H n h c j j j j d d c c --= = = ??- -- πωωπ ωπ ωαωω ωαω π π ω sin 21 21 用有限长单位脉冲响应序列()n h 逼近()n h d 。由于()n h d 往往是无限长序列，而且是非因果的，所以用窗函数()n ω将()n h d 截断，并进行加权处理，得到： ()()()n n h n h d ω= （4.2） ()n h 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数() ωj e H 为 ()()n j N n j e n h e H ωω ∑-==1 （4.3）式中，N 为所选窗函数()n ω的长度。用窗函数法设计的滤波器性能取决于窗函数()n ω的类型及窗口长度N 的取

FIR数字滤波器设计论文

毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名：日期：

学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期：年月日导师签名：日期：年月日

自适应滤波器的设计开题报告

自适应滤波器的设计与应用一、题目来源来源于其他二、研究目的和意义滤波技术在当今信息处理领域中有着极其重要的应用。滤波是从连续的或离散的输入数据中除去噪音和干扰以提取有用信息的过程，相应的装置就称为滤波器。滤波器实际上是一种选频系统，他对某些频率的信号予以很小的衰减，使该部分信号顺利通过。而对其他不需要的频率信号予以很大的衰减，尽可能阻止这些信号通过。滤波器研究的一个目的就是：如何设计和制造最佳的(或最优的)滤波器。在数字信号处理中，数字滤波是语音和图像处理、模式识别、频谱分析等应用中的一个基本处理算法。在许多应用场合，由于无法预先知道信号和噪声的特性或者它们是随时间变化的，仅仅用 FIR 和 IIR两种具有固定滤波系数的滤波器无法实现最优滤波。在这种情况下，必须设计自适应滤波器，以跟踪信号和噪声的变化。自适应滤波器是利用前一时刻已获得的滤波器参数，自动地调节、更新现时刻的滤波器参数，以适应信号和噪声未知的统计特性，从而实现最优滤波。当在未知统计特性的环境下处理观测信号时，利用自适应滤波器可以获得令人满意的效果，其性能远超过通用方法所设计的固定参数滤波器。

三、阅读的主要参考文献及资料名称 1、《数字信号处理》刘益成（第二版）西安电子科技出版社 2、《数字信号处理》张小虹（第二版）机械工业出版社 3、自适应信号处理[M]．西安：西安电子科技大学出版社，2001． 4．邹理和,数字信号处理, 国防工业出版社,1985 5．丁玉美等, 数字信号处理,西安电子科技大学出版社,1999 6．程佩青, 数字信号处理,清华大学出版社,2001 7. The MathWorks Inc, Signal Processing Toolbox For Use with MATLAB, Sept. 2000 8. vinay K.Ingle, John G.Proakis,数字信号处理及MATLAB实现,陈怀琛等译,电子工业出版社,1998.9 9、《MATLAB编程参考手册》 10、中国期刊网的相关文献 11、赫金，自适应滤波器原理第四版，西安工业出版社，2010-5-1 四、国内外现状和发展趋势与主攻方向自适应滤波器的理论与技术是50年代末和60年代初发展起来的。它是现代信号处理技术的重要组成部分，对复杂信号的处理具有独特的功能。自适应滤波器在数字滤波器中试属于随机数字信号处理的范畴。对于随机数字信号的滤波处理，通常有维纳滤波，卡尔曼滤波和自适应滤波，维纳滤波的权系数是固定的，适用于平稳随机信号；卡尔曼滤波器的权系数是可变的，适用于非平稳随机信号中。但是，只有在对信号和噪声的统计特性先验

利用MATLAB窗函数法设计一个可实现的FIR低通滤波器。

一、实验目的 1.掌握在MATLAB中窗函数的使用方法，了解不同窗函数之间的差别。 2.使用窗函数法设计一个可实现的FIR低通滤波器。 3.观察在相同长度下，不同的窗函数设计出来的滤波器有什么差别。 4.观察同一个窗在不同长度下设计出来的滤波器有什么差别。二、实验条件 PC机，MATLAB7.0 三、实验内容 1）通过help查找窗函数在MATLAB中如何实现

通过example了解MATLAB中窗函数的实现，并且利用矩形窗，汉宁窗，哈明窗，布莱克曼窗和凯塞窗来进行接下来的实验。 2）设计物理可实现的低通滤波器设计思路：因为要设计FIR有限脉冲响应滤波器，通常的理想滤波器的单位脉冲响应h是无限长的，所以需要通过窗来截断它，从而变成可实现的低通滤波器。程序如下：clc;clear all; omga_d=pi/5; omga=0:pi/30:pi; for N=3:4:51; w1= window(@blackman,N); w2 = window(@hamming,N); w3= window(@kaiser,N,2.5); w4= window(@hann,N); w5 = window(@rectwin,N); M=floor(N/2); subplot(311);plot(-M:M,[w1,w2,w3,w4,w5]); axis([-M M 0 1]); legend('Blackman','Hamming','kaiser','hann','rectwin'); n=1:M; hd=sin(n*omga_d)./(n*omga_d)*omga_d/pi; hd=[fliplr(hd),1/omga_d,hd]; h_d1=hd.*w1';h_d2=hd.*w2';h_d3=hd.*w3';h_d4=hd.*w4';h_d5=hd.*w5'; m=1:M; H_d1=2*cos(omga'*m)*h_d1(M+2:N)'+h_d1(M+1); H_d2=2*cos(omga'*m)*h_d2(M+2:N)'+h_d2(M+1); H_d3=2*cos(omga'*m)*h_d3(M+2:N)'+h_d3(M+1); H_d4=2*cos(omga'*m)*h_d4(M+2:N)'+h_d4(M+1); H_d5=2*cos(omga'*m)*h_d5(M+2:N)'+h_d5(M+1); subplot(312);plot(omga,[H_d1,H_d2,H_d3,H_d4,H_d5]); legend('Blackman','Hamming','kaiser','hann','rectwin'); subplot(313);plot(abs([fft(h_d1);fft(h_d2);fft(h_d3);fft(h_d4);fft(h_d5)] )'); pause(); end 程序分析：整个对称窗的长度为N，然而为了在MATLAB中看到窗函数在负值时的形状需将N 变为它的一半，即为2M+1个长度。窗长设置为从3开始以4为间隔一直跳动51。则长度相同的不同窗函数在时域[-M,M]的形状如第一个图所示。对窗函数进行傅里叶变换时，将零点跳过去先构造一个一半的理想滤波器的脉冲响应hd，再将零点位置求导得出的数赋值进去。将生成的hd左右颠倒形成了一个理想的滤波器的脉冲响应。将构造的理想滤波器的脉冲响应依次与之前定义的窗函数相乘，相乘出来的为列向量，用转置将其变成行向量，形成的h_d就是非理想的低通滤波器的脉冲响应序列。因为h_d为对称奇数长度序列，它的DTFT可以是二倍的离散余弦变化，而零点的位置则直接带入求出，两者相加则是H_d。则第二个图表示的是五个矩阵向量在频域的变化，而第三个图表示的是五个非理想低

FIR数字滤波器设计及软件实现

实验五：FIR数字滤波器设计及软件实现一、实验目的：（1）掌握用窗函数法设计FIR数字滤波器的原理和方法。（2）掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。（3）掌握FIR滤波器的快速卷积实现原理。（4）学会调用MATLAB函数设计与实现FIR滤波器。二、实验内容及步骤：（1）认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理；（2）调用信号产生函数xtg产生具有加性噪声的信号xt，并自动显示xt及其频谱，如图1所示；图1 具有加性噪声的信号x(t)及其频谱如图（3）请设计低通滤波器，从高频噪声中提取xt中的单频调幅信号，要求信号幅频失真小于0.1dB，将噪声频谱衰减60dB。先观察xt的频谱，确定滤波器指标参数。（4）根据滤波器指标选择合适的窗函数，计算窗函数的长度N，调用MATLAB 函数fir1设计一个FIR低通滤波器。并编写程序，调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。（4）重复（3），滤波器指标不变，但改用等波纹最佳逼近法，调用MATLAB 函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。友情提示： ○1MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本课本；

○ 2采样频率Fs=1000Hz ，采样周期T=1/Fs ； ○ 3根据图10.6.1(b)和实验要求，可选择滤波器指标参数：通带截止频率fp=120Hz ，阻带截至频率fs=150Hz ，换算成数字频率，通带截止频率p 20.24p f ωπ=T =π，通带最大衰为0.1dB ，阻带截至频率s 20.3s f ωπ=T =π，阻带最小衰为60dB 。] ○ 4实验程序框图如图2所示。图2 实验程序框图三、实验程序： 1、信号产生函数xtg 程序清单： %xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz %载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz. function xt=xtg N=1000;Fs=1000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T; fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=Fc/10;

基于Matlab的FIR低通滤波器设计[开题报告]2011-11-07

电子信息学院本科毕业设计(论文)开题报告论文题目基于Matlab的FIR低通滤波器设计学生姓名专业班级性别一、课题研究意义及现状所谓数字信号处理技术是指运用数值计算的方法对信号进行分析、变换、综合、估值与识别等的处理。随着电子技术和电子计算机技术的发展，数字信号处理技术受到了越来越广泛的关注，数字信号处理的理论和技术也在不断丰富和完善，新的理论和新技术层出不穷。可以说，数字信号处理是发展最快、应用最广泛、成效最显著的新科学之一，目前已广泛地应用在语音、雷达、声纳、地震、图像、通信、控制、生物医学、遥感遥测、地质勘探、航空航天、故障检测、自动化仪表等领域。数字滤波器是数字信号处理中的重要组成部分之一，在数字信号处理中，数字滤波占有极其重要的地位。滤波器可以分为IIR滤波器和FIR滤波器两大类，IIR滤波器并不能得到严格的线性相位特性，因此在许多实际应用中为了得到线性相位特征，还必须另外增加相位校正网络，使滤波器设计变得复杂，成本也高，而FIR滤波器在保证幅度特性满足技术要求的同时，还可以很容易的做到有严格的线性相位特性。它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。目前对数字滤波器的设计有多种方法，其中著名的Matlab软件包，功能强大、使用方便。传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。利用Matlab信号处理工具箱可以快速有效的设计由软件组成的常规数字滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。本课题在基于Matlab的FIR低通滤波器的设计方法、编程以及实际应用上有着积极的意义。

matlab设计低通滤波器

个matlab程序怎么编？（设计低通滤波器）通带边缘频率10khz 阻带边缘频率22khz 阻带衰减75db 采样频率50khz 要求设计这个低通滤波器画出脉冲响应的图形还有滤波器的形状具体程序怎么编？谢谢各位大虾的指点！！！最佳答案 1.1 实验目的 1．了解数字信号处理系统的一般构成； 2．掌握奈奎斯特抽样定理。 1.2 实验仪器 1．YBLD智能综合信号源测试仪1台 2．双踪示波器1台 3．MCOM－TG305数字信号处理与现代通信技术实验箱1台 4．PC机（装有MATLAB、MCOM－TG305配套实验软件）1台 1.3 实验原理一个典型的DSP系统除了数字信号处理部分外，还包括A/D和D/A两部分。这是因为自然界的信号，如声音、图像等大多是模拟信号，因此需要将其数字化后进行数字信号处理，模拟信号的数字化即称为A/D转换。数字信号处理后的数据可能需还原为模拟信号，这就需要进行D/A转换。一个仅包括A/D和D/A两部分的简化数字信号处理系统功能如图1所示。 A/D转换包括三个紧密相关的过程，即抽样、量化和编码。A/D转换中需解决的以下几个重要问题：抽样后输出信号中还有没有原始信号的信息？如果有能不能把它取出来？抽样频率应该如何选择？

奈奎斯特抽样定理（即低通信号的均匀抽样定理）告诉我们，一个频带限制在0至fx以内的低通信号x(t)，如果以fs≥2fx的抽样速率进行均匀抽样，则x(t)可以由抽样后的信号xs(t)完全地确定，即xs(t)包含有x(t)的成分，可以通过适当的低通滤波器不失真地恢复出x(t)。最小抽样速率fs=2fx称为奈奎斯特速率。低通译码编码量化抽样输入信号样点输出滤波输出 A/D（模数转换）D/A（数模转换）图1 低通采样定理演示为方便实现，实验中更换了一种表现形式，即抽样频率固定（10KHz），通过改变输入模拟信号的频率来展示低通抽样定理。我们可以通过研究抽样频率和模拟信号最高频率分量的频率之间的关系，来验证低通抽样定理。 1.4 实验内容 1．软件仿真实验：编写并调试MATLAB程序，分析有关参数，记录有关波形。 2．硬件实验：输入不同频率的正弦信号，观察采样时钟波形、输入信号波形、样点输出波形和滤波输出波形。 1.5 MATLAB参考程序和仿真内容 %*******************************************************************% %f—余弦信号的频率