基于窗函数法的FIR数字低通滤波器设计docx.docx
- 文档编号:6609888
- 上传时间:2023-01-08
- 格式:DOCX
- 页数:12
- 大小:40.85KB
基于窗函数法的FIR数字低通滤波器设计docx.docx
《基于窗函数法的FIR数字低通滤波器设计docx.docx》由会员分享,可在线阅读,更多相关《基于窗函数法的FIR数字低通滤波器设计docx.docx(12页珍藏版)》请在冰豆网上搜索。
基于窗函数法的FIR数字低通滤波器设计docx
基于窗函数法的FIR数字低通滤波器设计
摘要
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
根据其单位冲激响应函数的时域特性可分为两类:
无限冲激响应(IIR)滤波器和有限冲激响应(FIR)滤波器。
与IIR滤波器相比,FIR的实现是非递归的,总是稳定的;更重要的是,FIR滤波器在满足幅频响应要求的同时,可以获得严格的线性相位特性。
因此,它在高保真的信号处理,如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。
滤波器的设计是信号处理的核心问题之一。
根据FIR滤波器的原理,提出了FIR滤波器的窗函数设计法,给出了在MATLAB环境下,用窗函数法设计FIR滤波器的过程和设计实例。
通过利用不同的窗函数方法设计FIR滤波器,对所设计的滤波器进行分析比较,得出各种方法设计的滤波器的优缺点及其不同的使用场合,从而可以在设计滤波器时能够正确的选择FIR数字滤波器的窗函数的选取及设计方法。
关键词:
FIR滤波器,MATLAB,窗函数
摘要I
1概述1
1.1FIR滤波器简介1
1.2窗函数设计法1
2设计原理3
2.1基本原理3
2.2典型的窗函数4
3几种数字低通滤波器的窗函数设计7
3.1采用矩形窗设计FIR数字低通滤波器7
3.2采用汉明窗设计FIR数字低通滤波器7
3.3采用布莱克曼窗设计FIR数字低通滤波器9
参考文献10
附录12
1概述
1.1FIR滤波器简介
FIR数字滤波器设计最简单的方法是窗函数法,通常也称为傅立叶级数法。
它是在时域进行的,因而必须由理想滤波器的频率响应Hd(eJW)推导出其单位冲激响应纪(〃),在设计一个FIR数字滤波器的单位冲激响应力(〃)去逼近纪(〃)。
根据冲激响应的时域特性,数字滤波器可分为无限长冲激响应(IIR)和有限长冲激响应滤波器(FIR),FIR的突出优点是:
系统总是稳定的、易于实现线性相位、允许设计多通带(或多阻带)滤波器,但与IIR相比,在满足同样阻带衰减的情况下需要的阶数较高,滤波器的阶数越高,占用的运算时间越多,因此在满足指标要求的情况下应尽量减少滤波器的阶数。
1.2窗函数设计法
窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲应响应序列的设计方法,通常在设计滤波器之前,应该先根据具体的工程应用确定滤波器的技术指标,在大多数实际应用中,数字滤波器常常被用来实现选频操作,所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。
用窗函数法设计FIR滤波器的步骤如下:
(1)根据过渡带宽及阻带衰减要求,选择窗函数的类型并估计窗口长度N(或阶数M=N-1),窗函数类型可根据最小阻带衰减As独立选择,因为窗口长度N对最小阻带衰减As没有影响,在确定窗函数类型以后,可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N,设待求滤波器的过渡带宽为△(«),它与窗口长度N近似成反比,窗函数类型确定后,其计算公式也确定了,不过这些公式是近似的,得出的窗口长度还要在计算中逐步修正,原则是在保证阻带衰减满足要求的情况下,尽量选择较小的N,在N和窗函数类型确定后,即可调用MATLAB中的窗函数求出窗函数W(n)。
(2)根据待求滤波器的理想频率响应求出理想单位脉冲响应hd(n),如果给出待求滤波器频率应为Hd(』w),则理想的单位脉冲响应可以用下面的傅里叶反变换式求出:
侦〃)=
(1)
在一般情况下,hd(n)是不能用封闭公式表示的,需要采用数值方法表示;从w=0到w=2兀采样N点,采用离散傅里叶反变换(IDFT)即可求出。
(3)计算滤波器的单位脉冲响应h(n),它是理想单位脉冲响应和窗函数的乘积。
(4)算技术指标是否满足要求,为了计算数字滤波器在频域中的特性,可调用freqz子程序,如果不满足要求,可根据具体情况,调整窗函数类型或长度,直到满足要求为止。
使用窗函数法设计时要满足以下两个条件:
1窗谱主瓣尽可能地窄,以获得较陡的过渡带;
2尽量减少窗谱的最大旁瓣的相对幅度,也就是使能量尽量集中于主瓣,减小峰肩和纹波,进行增加阻带的衰减。
窗函数的选择原则是:
[1]具有较低的旁瓣幅度,尤其是第一旁瓣的幅度;
[2]旁瓣的幅度下降的速率要快,以利于增加阻带的衰减;
[3]主瓣的宽度要窄,这样可以得到比较窄的过渡带。
通常上述的几点难以同时满足。
实际中设计FIR数字滤波器往往要求是线性相位的,因此要求w(〃)满足线性相位的条件,即要求w(n)满足:
w(〃)=w(N—1一〃)
(2)
所以,窗函数不仅有截短的作用,而且能够起到平滑的作用,在很多领域得到了应用。
2设计原理
2.1基本原理
设计低通FIR数字滤波器,寻求一系统函数H(z),使其频率响应H")逼近滤波器要求的理想频率响应H"),其对应的单位脉冲响应知(〃)
ei(0
厂冲加|<(0<
0,mc<|(y|<
(3)
甘r-HN
其中a=
2
")=券Ico
引如广_sin[外(〃-。
)]c>c>LCCU—/\
Infji\n-a)
如果所希望的滤波器的理想的频率响应函数为Hd
eit0),则其对应的单位脉冲响
应为
(4)
窗函数设计法的基本原理是用有限长单位脉冲响应序列/z(")逼近如(〃)。
由于如(〃)往往是无限长序列,而且是非因果的,所以用窗函数必7)将如(〃)截断,并进行加权处理,得到:
(5)
/?
(〃)=hd
/z(")就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数
h")为
(6)
丑")=勇枷/
n=0
式中,N为所选窗函数刃(")的长度。
2.2典型的窗函数
(7)
(1)矩形窗(RectangleWindow)
w(〃)=Rn(〃)
其频率响应和幅度响应分别为:
川(此)=sin(N%2)/咛,吼㈣='出(此/2)
sin(刃/2)
sin(刃/2)
(8)
(2)三角形窗(BartlettWindow)
2n
N-l
(3)
n-T
2-
N—1
2n
其频率响应为:
N-l
(9)
W(/)=Z[Sin(AW4)]2/" Nsin((z>/2) 汉宁(Hanning)窗,又称升余弦窗 (10) =|[1-cos(^^)]Rn(") ZTV—1 其频率响应和幅度响应分别为: W(e勺={0.5吼(®)+0.25[吼-—)+WR^+肴)诉 5。 5+。 .25[住-"住+万二)] (4) 汉明(Hamming)窗,又称改进的升余弦窗 (11) w(〃)=[0.54-0.46cos(B今&(〃) N—1 其幅度响应为: Wo)=0.54Wr(刃)+0.23[Wr一趾丁)+Wr(口+-^-)] (5) 布莱克曼(Blankman)窗,又称二阶升余弦窗 (12) 2n7T4^771 w(〃)=[0.42-0.5cos(^—)+0.08cos(—)]RN(〃) N-1N-1 其幅度响应: 271271 W㈣=0.42W.例+0.25[Wr-—-)+WR^+—)] N—lN-1 4勿4勿 +0.04[吟-—)+吟(刃+—)] (6)凯塞(Kaiser)窗 响=""2曾顼2,0割景_1 IM (13) 其中: B是一个可选参数,用来选择主瓣宽度和旁瓣衰减之间的交换关系,一 般说来,B越大,过渡带越宽,阻带越小衰减也越大。 Io(・)是第一类修正零阶贝 塞尔函数。 若阻带最小衰减表示为&=-20log】。 5,,3的确定可采用下述经验公式: 04<21
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 函数 FIR 数字 滤波器 设计 docx