基于MATLAB的FIR数字低通滤波器设计综述Word格式.docx
- 文档编号:12946047
- 上传时间:2022-10-01
- 格式:DOCX
- 页数:25
- 大小:304.20KB
基于MATLAB的FIR数字低通滤波器设计综述Word格式.docx
《基于MATLAB的FIR数字低通滤波器设计综述Word格式.docx》由会员分享,可在线阅读,更多相关《基于MATLAB的FIR数字低通滤波器设计综述Word格式.docx(25页珍藏版)》请在冰豆网上搜索。
Tutor:
JingYingMin
Abstract:
FIRdigitalfilterisanimportantpartofdigitalsignalprocessing,theFIRdigitalfilterwithlinearphase,soithasbeenwidelyappliedinthecollectionandprocessingofinformationinthecourseof.ThispaperintroducestheconceptofFIRdigitalfilterwithlinearphaseconditions,analysisofthewindowfunctionmethodandfrequencysamplingmethodandtherippleapproximationmethodofFIRfilterdesignideasandprocesses.Basedonanalyzingtheprincipleofthreekindsofdesignmethod,usingMatlabsimulationsoftwarefir1toolboxindesignandimplementationofFIRHammingwindowwindowfunctionmethodinthelowpassfilter.
Keywords:
FIRdigitalfilter,linearphase,thewindowfunctionmethod,Hammingwindow,Matlab
引言..............................................................................................................................................................1
引言
随着信息科学和计算机技术的不断发展,数字信号处理(DSP,DigitalSignalProcessing)的理论和技术也得到了飞速的发展,并逐渐成为一门重要的学科,它的重要性在日常通信、图像处理、遥感、声纳、生物医学、地震、消费电子、国防军事、医疗方面等显得尤为突出。
在我们面临的信息革命中,数字信号处理几乎涉及了所有的工程技术领域[1]。
数字信号处理是一种将信号以数字形式进行处理的一种理论和技术,它的目的是将真实世界中的一些信号进行分析并滤波,最后得出其中的有用的信号。
数字滤波器是数字信号处理的一种,一般根据单位脉冲响应h(n)分为无限脉冲响应(IIR)和有限脉冲响应(FIR)系统。
IIR数字滤波器的设计方法简单,特别是采用双线性变换法来设计的数字滤波器不存在频域混叠的现象,但是IIR滤波器存在一个较为明显的缺憾,就是它的相位响应一般都是非线性的,而在传输频带内的相位响应如果不是线性的,就会造成有用信号的传输失真,而FIR数字滤波器不仅可以设计成任意的幅度响应,而且可以设计成在通频带内具有良好的线性相位响应[2]。
FIR数字滤波器的单位脉冲响应h(n)有限长,所以FIR数字滤波器是稳定的,不存在稳定性的问题,且可以通过快速傅里叶变换(FFT)的算法来实现信号滤波,大大的提高的运算效率。
因此,FIR数字滤波器日益引起了人们的关注。
FIR数字滤波器的设计方法有很多,比较常用的有窗函数设计法、频率采样设计法、等波纹逼近法等。
本课题通过运用窗函数设计FIR数字低通滤波器,并实现对给定的信号进行滤波。
窗函数设计法是最基本的数字滤波方法,是利用傅里叶反变换(IDTFT)计算给定的频响的理想单位脉冲响应,再加以窗函数进行截断和平滑[3]。
Matlab软件的信号处理工具箱提供了FIR数字滤波器设计的子函数,运用Matlab软件设计可以避免繁杂的数学运算,而且具有丰富的绘图功能,可以方便地查看所设计的数字滤波器的幅度响应和相位响应是否满足设计要求。
因此,本课题在理论分析各种FIR数字滤波器设计方法的基础上,运用Matlab软件进行仿真分析。
1.基本原理
1.1FIR数字滤波器概述
一般来说一个经典的数字滤波器是一个线性时不变系统,其数学模型可以用Z域系统函数
来表示:
(1-1)
其中
均为滤波器参数。
在式(1-1)中,当
值不全为零值时,Z域系统函数
必定含有一个或一个以上的极值点,此时单位脉冲响应为无限长,对于一个稳定的数字滤波器来说,Z域系统函数
必须在单位圆内,因而把含有极值点的Z域系统函数
的数字滤波器称为无限脉冲响应数字滤波器(InfiniteImpulseResponse),即IIR数字滤波器[4]。
而当
值全为零时,Z域的系统函数
只有一个零点,式(1-1)表示的系统函数
可以写成:
(1-2)
公式(1-2)表明,FIR滤波器的系统函数是
的
阶多项式,在有限
平面
上有
个零点,而在
平面原点
处有
阶极点。
公式(1-2)表示的系统,其单位脉冲响应可以表示为:
(1-3)
在式(1-3)中,只有当
,
才有非零值,所以数字滤波器的脉冲响应是有限长的,因此在数字信号处理中把这种数字滤波器称为有限脉冲响应数字滤波器(FiniteImpuleseResponse),即FIR数字滤波器。
FIR数字滤波器最突出的两个优点是:
(1)只要对
附加一定的条件,就很容易获得严格的线性相位。
(2)由于
的极点位于原点
处,所以FIR数字滤波器不存在稳定性问题。
1.2FIR数字滤波器线性相位定义
设FIR数字滤波器脉冲响应的长度为N,则其频率响应可以表示为:
(1-4)
上式通过欧拉恒等式展开可得到
的相位特性
,有两种线性相位特性,通常称为第一类线性相位和第二类线性相位。
第一类线性相位特性:
是一个与
无关的常数
第二类线性相位特性:
是起始相位
严格地说第二种情况时的
是不具有线性相位特性的,但上述两种情况都满足群延迟是一个常数,仍可以视为具有线性相位的,在第二类线性相位中
是常用的一种情况[5]。
1.3FIR数字滤波器线性相位时域条件
对于第一类线性相位,即
,通过一系列的运算整理之后可得到一个三角函数求和公式:
(1-5)
式中正弦函数
为奇对称,当
时,对称中心为
需要满足关于
偶对称,即要求:
(1-6)
对于第二类线性相位,即
时,通过运算得到公式:
(1-7)
函数
为偶对称,当
时,对称中心也为
。
若要使上式成立,则要使
关于
奇对称,即要求:
(1-8)
从上述分析看来,线性相位FIR数字滤波器的时域约束条件是指满足线性相位时对
的约束条件,对于第一类线性相位,冲激响应
满足(1-6)式;
对于第二类线性相位,冲激响应h(n)满足(1-8)式。
2.系统设计
FIR数字滤波器的设计方法主要有窗函数设计法、频率采样设计法以及等波纹逼近设计法三种,其中窗函数设计法是最常用的,其次是频率采样法,但这两种方法在设计中还会存在一些不足之处,所以需要优化的设计方法,而等波纹逼近法很好的弥补了窗函数法和频率采样法的不足[6]。
2.1FIR数字滤波器的窗函数设计方法
2.1.1窗函数法的设计思路
窗函数设计法是FIR数字滤波器里最简单的一种设计法,又叫傅里叶级数法,为了设计简单方便,通常选择所希望逼近的滤波器的频率响应函数
为具有片段常数特性的理想滤波器,寻找一组
,确定其频率响应
,然后用
来逼近
[1]。
窗函数法设计FIR滤波器是在时域中进行的,那么可以通过傅里叶反变换得到得到频率响应
,即:
(2-1)
在实际中,
一般是处于逐段恒定的,在边界频率处有不连续点,因而单位脉冲响应
是无限长的非因果序列,不能直接作为FIR数字滤波器的单位脉冲响应,因此需要对
进行截断,转换为有限长的一段因果序列,也就是用一个有限长度的窗函数序列
来截取
,即
,并将非因果序列转变为一个因果序列。
截取的长度和加权窗函数的形状都直接影响到逼近精度。
以截止频率为
,相位为零的理想低通滤波器为例,其频率特性为:
(2-2)
通过傅里叶反变换得到对应的
为:
(2-3)
此时的
是一个无限长的非因果序列,我们需要对其进行截断,变成一个有限长的因果序列。
可以先把
向右平移
个点,得到
(2-4)
相应的传输函数
(2-5)
然后对
截取从0到
的N个点,N为窗函数的长度,所得的结果
表示为:
(2-6)
表示窗函数,一般用下标来表示窗函数的类型,矩形窗记为
2.1.2常见窗函数介绍
常见的窗函数有矩形窗(RectangleWindow)、三角形窗(BartlerrWindow)、汉宁(Hanning)窗——升余弦窗、哈明(Hamming)窗——改进的升余弦窗、布莱克曼(Blackman)窗、凯塞—贝塞尔窗(Kaiser-BaselWindow)[7]。
矩形窗的窗函数为:
(2-7)
其频谱的幅度函数为
(2-8)
矩形窗的主瓣宽度为
,用矩形窗设计的FIR数字滤波器的过渡带宽度近似为
三角形窗的窗函数为:
(2-9)
(2-10)
三角窗的主瓣宽度为
汉宁窗窗函数为:
(2-11)
汉宁窗的频谱的幅度函数为
(2-12)
汉宁窗的主瓣宽度为
,汉宁窗在其两个端点都为零,实际中这两个端点的数据是不可用的。
哈明窗的窗函数为
(2-13)
其幅度函数为
(2-14)
哈明窗是一种改进的余弦窗,能量更加集中在主瓣,是一种高效的窗函数,主瓣宽度与汉宁窗的相同。
布莱克曼窗窗函数为
(2-15)
(2-16)
该窗函数位移不同,幅度函数也不同,会使旁瓣进一步抵消,主瓣宽度为
凯塞窗是一种最优窗函数,不同于前面五种窗函数,凯塞窗是一种参数
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 MATLAB FIR 数字 滤波器 设计 综述