FIR滤波器窗函数设计要点
数字信号处理实验——用窗函数设计FIR滤波器

实验四 用窗函数设计FIR 滤波器一、 实验目的1、熟悉FIR 滤波器设计的基本方法。
2、掌握用窗函数设计FIR 数字滤波器的原理及方法,熟悉相应的计算机高级语言编程。
3、熟悉线性相位FIR 滤波器的幅频特性和相位特性。
4、了解各种不同窗函数对滤波器性能的响应。
二、 实验原理和方法窗函数法设计的任务在于寻找一个可实现有限长单位脉冲响应的传递函数H(e jw )=∑-=10N n h(n)e -jwn 去逼近h d (n)=1/2π⎰π20H d (e jw )e jwn dw即h(n)=h d (n)w (n ) (一)几种常用的窗函数1、矩形窗 w(n)=R N (n)2、Hanning 窗 w(n)=0.5[1-cos(2πn /N-1)]R N (n)3、Hamming 窗 w(n)=[0.54-0.46cos(2πn /N-1)]R N (n)4、Blackman 窗 w(n)=[0.42-0.5 cos(2πn /N-1)+0.08 cos(4πn /N-1)] R N (n)5、Kaiser 窗 w(n)=I 0(β(1-[(2n /(N-1))-1]2)½)/I 0(β)(二)窗函数法设计线性相位FIR 滤波器的步骤1、确定数字滤波器的性能要求。
确定各临界频率{w k }和滤波器单位脉冲响应长度N 。
2、根据性能要求和N 值,合理地选择单位脉冲响应h(n)有奇偶对称性,从而确定理想频率响应h d (e jw)的幅频特性和相位特性。
3、用傅里叶反变换公式求得理想单位脉冲响应h d (n)。
4、选择适当的窗函数W (n ),求得所设计的FIR 滤波器单位脉冲响应。
5、用傅里叶变换求得其频率响应H (e jw),分析它的幅频特性,若不满足要求,可适当改变窗函数形式或长度N ,重复上述过程,直至得到满意的结果。
三、实验内容和步骤1、分别用矩形窗、Hanning 窗、Hamming 窗、Blackman 窗、Kaiser 窗(β=8.5)设计一个长度N=8的线性相位FIR 滤波器。
实验7 窗函数法设计FIR数字滤波器

实验7窗函数法设计FIR数字滤波器一、实验目的掌握窗函数法设计F1R数字滤波器的原理和具体方法二、实验设备与环境计算机、Mat1ab软件环境三、实验基础理论1>基本原理窗函数设计法的基本思想为,首先选择一个适当的理想的滤波器Hd(,3),然后用窗函数截取它的单位脉冲响应%(九),得到线性相位和因果的FIR滤波器,这种方法的重点是选择一个合适的窗函数和理想滤波器,使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。
2、设计步骤(1)给定理想滤波器的频率响应Hd("3),在通带上具有单位增益和线性相位,在阻带上具有零响应。
一个带宽为g(3c<Tr)的低通滤波器由下式给定h(e j^=(eW∣ω∣≤ωc虱)一1Oωc<∣ω∣<π其中α为采样延迟,其作用是为了得到因果的系统。
(2)确定这个滤波器的单位脉冲响应为了得到一个h(n)长度为N的因果的线性相位FIR滤波器,我们令N-Ia=-2-(3)用窗函数截取hd(τι)得到所设计FIR数字滤波器h(n)h(n)=h d(n)w(n)3、窗函数的选择常用的窗函数有矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗等。
Mat1ab提供了一些函数用于产生窗函数,如下表所示:在设计过程中我们需要根据给定的滤波器技术指标,选择滤波器长度N 和窗函数3(n)°表7.2列出了常用的窗函数的一些特性,可供设计时参考。
其中幻是修正的零阶贝塞尔函数,参数B 控制最小阻带衰减,这种窗函数对于相同的N 可以提供不同的过渡带宽。
由于贝塞尔函数比较更杂,这种窗函数的设计方程很难推导,然而幸运的是,有一些经验设计方程可以直接使用。
已知给定的指标叫Msc,Rp 和4,滤波器长度N 和凯瑟窗参数B 可以按如下凯瑟窗设计方程给出过渡带宽:∆ω=ωst -ωp入一7.95 2.285∆ω_(0.1102(4-8.7) ,P=iθ.5842(4-21)04+0.07886(4-21), 四、实验内容1、设计一个数字低通FIR 滤波器,其技术指标如下ωp =0.2τr,RP=0.25dBωst =0.3τr,A s =50dB分别采用矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。
FIR滤波器设计要点

FIR滤波器设计要点FIR (Finite Impulse Response) 滤波器是一种数字滤波器,其设计要点包括滤波器类型选择、滤波器系数设计、频率响应规格、窗函数和滤波器长度选择等。
以下是对这些要点的详细介绍。
1.滤波器类型选择:在设计FIR滤波器之前,需要确定滤波器的类型。
常见的FIR滤波器类型包括低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
不同类型的滤波器适用于不同的应用场景,因此在选择滤波器类型时需要考虑系统的需求。
2.滤波器系数设计:FIR滤波器的核心是滤波器系数。
滤波器系数决定了滤波器的频率响应和滤波特性。
常用的设计方法包括窗函数法、最小均方误差法和频率抽样法等。
窗函数法是最常用的设计方法,其基本思想是通过选择合适的窗函数来得到滤波器系数。
3.频率响应规格:在设计FIR滤波器时,需要明确所需的频率响应规格,包括通带增益、阻带衰减、过渡带宽等。
这些规格直接影响了滤波器的性能,因此需要根据具体应用场景来确定。
4.窗函数选择:窗函数在FIR滤波器设计中起到了重要的作用。
常用的窗函数包括矩形窗、汉宁窗、汉明窗、布莱克曼窗等。
不同的窗函数具有不同的特性,选择合适的窗函数可以得到优良的滤波器性能。
5.滤波器长度选择:滤波器长度决定了滤波器的频率分辨率和时间分辨率。
滤波器长度越长,频率响应越尖锐,但计算复杂度也越高。
因此,在设计FIR滤波器时需要权衡计算复杂度和性能要求,选择合适的滤波器长度。
6.优化设计:7.实现方式:总之,设计FIR滤波器要点包括滤波器类型选择、滤波器系数设计、频率响应规格、窗函数和滤波器长度选择等。
设计者需要根据具体的应用场景和性能要求来进行合理的设计和优化,以满足系统的需求。
用窗函数设计FIR滤波器

实验六 用窗函数设计FIR 滤波器1.实验目的(1) 熟悉FIR 滤波器设计的方法和原理(2) 掌握用窗函数法设计FIR 滤波器的方法和原理,熟悉滤波器的特性 (3) 了解各种窗函数滤波器特性的影响2.实验原理FIR 滤波器的设计方法主要有三种:窗函数法、频率取样法、切比雪夫等波纹逼近法。
FIR 滤波器的设计是要寻求一系统函数)(z H ,使其频率响应)(ωj eH 逼近滤波器要求的理想频率响应,其对应的单位脉冲响应)(n h d 。
(1)用窗函数设计FIR 滤波器的基本方法在时域用一个窗函数截取理想的)(n h d 得到)(n h ,以有限长序列)(n h 近似逼近理想的)(n h d ;在频域用理想的)(ωj d e H 在单位圆上等角度取样得到h(k),根据h(k)得到H(z)将逼近理想的Hd(z)。
设理想滤波器)(ωj d e H 的单位脉冲响应为)(n h d 。
以低通线性相位FIR 数字滤波器为例。
⎰∑--∞-∞===ππωωωωωπd e e H n h e n h e H jn j d d jn n dj d )(21)()()()(n h d 一般是无限长的、非因果的,不能直接作为FIR 滤波器的单位脉冲响应。
要想得到一个因果的有限长的滤波器h(n),最直接的方法是截断)()()(n w n h n h d =,即截取为有限长因果序列,并用合适的窗函数进行加权作为FIR 滤波器的单位脉冲响应。
按照线性相位滤波器的要求,h(n)必须是偶对称的。
对称中心必须等于滤波器的延时常数,即⎩⎨⎧-==2/)1()()()(N a n w n h n h d 用矩形窗设计的FIR 低通滤波器,所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象,且最大波纹大约为幅度的9%,(现象称为吉布斯(Gibbs )效应)。
(2)典型的窗函数(a )矩形窗(Rectangle Window))()(n R n w N =其频率响应和幅度响应分别为:21)2/sin()2/sin()(--=N j j eN e W ωωωω,)2/sin()2/sin()(ωωωN W R =在matlab 中调用w=boxcar(N)函数,N 为窗函数的长度 (b )三角形窗(Bartlett Window)⎪⎩⎪⎨⎧-≤<----≤≤-=121,122210,12)(N n N N n N n N n n w其频率响应为:212])2/sin()4/sin([2)(--=N j j e N N e W ωωωω在matlab 中调用w=triang(N)函数,N 为窗函数的长度(c )汉宁(Hanning)窗,又称升余弦窗)()]12cos(1[21)(n R N n n w N --=π其频率响应和幅度响应分别为:)]12()12([25.0)(5.0)()()]}12()12([25.0)(5.0{)()21(-++--+==-++--+=---N W N W W W e W eN W N W W e W R R R aj N j R R R j πωπωωωωπωπωωωωω在matlab 中调用w=hanning(N)函数,N 为窗函数的长度 (d )汉明(Hamming)窗,又称改进的升余弦窗)()]12cos(46.054.0[)(n R N n n w N --=π其幅度响应为:)]12()12([23.0)(54.0)(-++--+=N W N W W W R R R πωπωωω 在matlab 中调用w=hamming(N)函数,N 为窗函数的长度(e )布莱克曼(Blankman)窗,又称二阶升余弦窗)()]14cos(08.0)12cos(5.042.0[)(n R N n N n n w N -+--=ππ 其幅度响应为:)]14()14([04.0)]12()12([25.0)(42.0)(-++--+-++--+=N W N W N W N W W W R R R R R πωπωπωπωωω在matlab 中调用w=blackman(N)函数,N 为窗函数的长度(f )凯泽(Kaiser)窗10,)())]1/(21[1()(020-≤≤---=N n I N n I n w ββ其中:β是一个可选参数,用来选择主瓣宽度和旁瓣衰减之间的交换关系,一般说来,β越大,过渡带越宽,阻带越小衰减也越大。
fir滤波器窗函数设计法

fir滤波器窗函数设计法
FIR (Finite Impulse Response)滤波器的窗函数设计法是一种经典的数字滤波器设计方法。
它通过选择一个合适的窗函数来对滤波器的频率响应进行加权,从而实现对信号的滤波。
窗函数设计法的基本步骤如下:
1. 确定滤波器的设计规格:包括截止频率、通带和阻带的幅频响应要求等。
2. 根据设计规格,计算出滤波器的理想频率响应:可以使用理想滤波器的频率响应作为目标。
3. 选择一个合适的窗函数:常用的窗函数有矩形窗、汉宁窗、汉明窗等。
选择窗函数的关键是考虑到主瓣宽度和副瓣衰减的平衡。
4. 将选择的窗函数应用于理想频率响应上,得到加权后的频率响应。
5. 对加权后的频率响应进行反变换,得到滤波器的时域响应。
6. 根据需要,对时域响应进行截断或零增益处理,以满足设计规格。
7. 最后,根据计算得到的滤波器系数,可以通过巴特沃斯频率抽样公式将其转换为巴特沃斯直接型或传输函数型,以便在数字系统中实现滤波。
需要注意的是,FIR滤波器的窗函数设计法是一种近似方法,设计的滤波器无法完全符合理想要求。
设计过程中需要权衡主瓣宽度和副瓣衰减等因素,以及选择合适的截断或零增益处理方式,以获得满
足实际需求的滤波器性能。
实验六用窗函数法设计FIR滤波器分析解析

实验六用窗函数法设计FIR滤波器分析解析一、引言数字滤波器是数字信号处理中的重要组成部分。
滤波器可以用于去除噪声、调整频率响应以及提取感兴趣的信号。
有许多方法可以设计数字滤波器,包括窗函数法、频域法和优化法等。
本实验将重点介绍窗函数法设计FIR滤波器的原理和过程。
二、窗函数法设计FIR滤波器窗函数法是设计FIR滤波器的一种常用方法。
其基本原理是将滤波器的频率响应与理想滤波器的频率响应进行乘积。
理想滤波器的频率响应通常为矩形函数,而窗函数则用于提取有限长度的理想滤波器的频率响应。
窗函数的选择在FIR滤波器的设计中起着重要的作用。
常用的窗函数包括矩形窗、汉宁窗、汉明窗、布莱克曼窗等。
对于每种窗函数,都有不同的特性和性能指标,如主瓣宽度、副瓣抑制比等。
根据不同的应用需求,可以选择合适的窗函数。
窗函数法设计FIR滤波器的具体步骤如下:1.确定滤波器的阶数N。
阶数N决定了滤波器的复杂度,一般情况下,阶数越低,滤波器的简单度越高,但频率响应的近似程度也会降低。
2.确定滤波器的截止频率。
根据应用需求,确定滤波器的截止频率,并选择合适的窗函数。
3.根据窗函数长度和截止频率计算理想滤波器的频率响应。
根据所选窗函数的特性,计算理想滤波器的频率响应。
4.根据理想滤波器的频率响应和窗函数的频率响应,得到所需的FIR滤波器的频率响应。
将理想滤波器的频率响应与窗函数的频率响应进行乘积,即可得到所需滤波器的频率响应。
5.对所得到的频率响应进行逆傅里叶变换,得到时域的滤波器系数。
6.实现滤波器。
利用所得到的滤波器系数,可以通过卷积运算实现滤波器。
三、实验结果与分析本实验以Matlab软件为平台,利用窗函数法设计了一个低通滤波器。
滤波器的阶数为16,截止频率为500Hz,采样频率为1000Hz,选择了汉宁窗。
根据上述步骤,计算得到了所需的滤波器的频率响应和时域的滤波器系数。
利用这些系数,通过卷积运算,实现了滤波器。
为了验证滤波器的性能,将滤波器应用于输入信号,观察输出信号的变化。
归纳窗函数法设计fir滤波器的设计过程

归纳窗函数法设计fir滤波器的设计过程“归纳窗函数法设计FIR滤波器的设计过程”是一个涉及数字信号处理领域的内容。
FIR滤波器是一种常用的数字滤波器,广泛应用于信号处理、通信系统等领域。
设计FIR滤波器的过程中,窗函数是其中一个重要的设计工具。
本文将从基本概念、设计过程和实例分析三个方面来详细介绍归纳窗函数法设计FIR滤波器的设计过程。
一、基本概念归纳窗函数法是一种常用的FIR滤波器设计方法。
其基本思想是通过选取合适的窗函数来对滤波器的频率响应进行近似,从而实现对信号的滤波。
在设计过程中,需要明确以下几个基本概念:1.1 FIR滤波器FIR滤波器(Finite Impulse Response Filter)是一种无穷冲激响应的滤波器,其输出序列只与输入序列的有限个值有关。
FIR滤波器由一组滤波器系数决定,称为滤波器的冲激响应。
常见的FIR滤波器类型有低通滤波器、高通滤波器、带通滤波器等。
1.2 窗函数窗函数(Window Function)是信号处理中常用的一种数学函数,用于将信号在时域或频域上进行截断、平滑或加权。
在归纳窗函数法中,利用窗函数对滤波器的频率响应进行调整。
常见的窗函数有矩形窗、汉宁窗、汉明窗等。
二、设计过程下面将详细介绍归纳窗函数法设计FIR滤波器的设计过程。
2.1 确定滤波器的规格首先,需要明确滤波器的规格,包括滤波器的类型(低通、高通、带通等)、截止频率、通带和阻带的衰减要求等。
2.2 确定滤波器的阶数根据滤波器的规格,可以确定滤波器的阶数。
阶数决定了滤波器的复杂度和滤波器系数的数量,一般通过经验公式进行估算。
2.3 选择窗函数根据滤波器的阶数和设计要求,选择合适的窗函数。
不同的窗函数具有不同的频域特性,适用于不同的应用场合。
一般而言,矩形窗适用于近似理想滤波器,汉宁窗适用于要求较窄过渡带的滤波器,汉明窗适用于要求极窄过渡带的滤波器。
2.4 确定窗函数的参数根据滤波器的阶数和设计要求,确定窗函数的参数。
FIR滤波器窗函数设计.

课题名称:FIR滤波器窗函数设计FIR 滤波器窗函数设计引言:数字滤波器(Digital Filter)是指输入、输出都是离散时间信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。
在许多数字信号处理系统中,如图像信号处理等,有限冲激响应(FIR )滤波器是最常用的组件之一,它完成信号预调、频带选择和滤波等功能。
FIR 滤波器虽然在截止频率的边沿陡峭性能上不及无限冲激响应(IIR )滤波器,但是却具有严格的线性相位特性,稳定性好,能设计成多通带(或多阻带)滤波器组,所以能够在数字信号处理领域得到广泛的应用。
一、 数字滤波器的分类(1) 根据系统响应函数的时间特性分为两类1. FIR (Finite Impulse Response )数字滤波器网络()0,0[][]0,Mn k k b n M y n b x n k h n n =≤≤⎧=-⇔⎨⎩∑ 其他 特点:不存在反馈支路,其单位冲激响应为有限长。
2. IIR (Infinite Impulse Response )数字滤波器网络01[][][]M Nk k k k y n b x n k a y n k ===---∑∑ 特点:存在反馈支路,即信号流图中存在环路,其单位冲激响应为无限长。
(2) FIR 数字滤波器IIR 数字滤波器的区别1. 从性能上来说,IIR 滤波器传递函数包括零点和极点两组可调因素,对极点的惟一限制是在单位圆内。
因此可用较低的阶数获得高的选择性,所用的存储单元少,计算量小,效率高。
但是这个高效率是以相位的非线性为代价的。
选择性越好,则相位非线性越严重。
FIR 滤波器传递函数的极点固定在原点,是不能动的,它只能靠改变零点位置来改变它的性能。
所以要达到高的选择性,必须用较高的阶数;对于同样的滤波器设计指标,FIR 滤波器所要求的阶数可能比IIR 滤波器高5-10倍,但是 FIR 滤波器可以得到严格的线性相位。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
课题名称:FIR滤波器窗函数设计FIR 滤波器窗函数设计引言:数字滤波器(Digital Filter)是指输入、输出都是离散时间信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。
在许多数字信号处理系统中,如图像信号处理等,有限冲激响应(FIR )滤波器是最常用的组件之一,它完成信号预调、频带选择和滤波等功能。
FIR 滤波器虽然在截止频率的边沿陡峭性能上不及无限冲激响应(IIR )滤波器,但是却具有严格的线性相位特性,稳定性好,能设计成多通带(或多阻带)滤波器组,所以能够在数字信号处理领域得到广泛的应用。
一、 数字滤波器的分类(1) 根据系统响应函数的时间特性分为两类1. FIR (Finite Impulse Response )数字滤波器网络()0,0[][]0,M n k k b n M y n b x n k h n n =≤≤⎧=-⇔⎨⎩∑ 其他特点:不存在反馈支路,其单位冲激响应为有限长。
2. IIR (Infinite Impulse Response )数字滤波器网络01[][][]M Nk k k k y n b x n k a y n k ===---∑∑ 特点:存在反馈支路,即信号流图中存在环路,其单位冲激响应为无限长。
(2) FIR 数字滤波器IIR 数字滤波器的区别1. 从性能上来说,IIR 滤波器传递函数包括零点和极点两组可调因素,对极点的惟一限制是在单位圆内。
因此可用较低的阶数获得高的选择性,所用的存储单元少,计算量小,效率高。
但是这个高效率是以相位的非线性为代价的。
选择性越好,则相位非线性越严重。
FIR 滤波器传递函数的极点固定在原点,是不能动的,它只能靠改变零点位置来改变它的性能。
所以要达到高的选择性,必须用较高的阶数;对于同样的滤波器设计指标,FIR 滤波器所要求的阶数可能比IIR 滤波器高5-10倍,但是 FIR 滤波器可以得到严格的线性相位。
2. 从结构上看,IIR 滤波器必须采用递归结构,极点位置必须在单位圆内,否则系统将不稳定。
相反,FIR 滤波器只要采用非递归结构,不论在理论上还是在实际的有限精度运算中都不存在稳定性问题,因此造成的频率特性误差也较小。
此外FIR 滤波器可以采用快速傅里叶变换算法,在相同阶数的条件下,运算速度可以快得多。
3. 从设计工具看,IIR 滤波器可以借助于模拟滤波器的成果,因此一般都有有效的封闭形式的设计公式可供准确计算,计算工作量比较小,对计算工具的要求不高。
从上面的简单比较可以看到IIR 与FIR 滤波器各有所长,所以在实际应用时应该从多方面考虑来加以选择。
从使用要求上来看,在对相位要求不敏感的场合,如语言通信等,选用IIR 较为合适,这样可以充分发挥其经济高效的特点;对于图像信号处理,数据传输等以波形携带信息的系统,则对线性相位要求较高。
如果有条件,采用FIR 滤波器较好。
另外,不论IIR 和FIR ,阶数越高,信号延迟越大;同时在IIR 滤波器中,阶数越高,系数的精度要求越高,否则很容易造成有限字长的误差使极点移到单位园外。
因此在阶数选择上是综合考虑的。
二、 FIR 数字滤波器的设计(一) FIR 滤波器原理FIR 滤波器的系统输入输出差分方程为:10[]()()N k y n h k x n k -==-∑ 所以FIR 滤波器的系统函数为: 1()()()()N nn Y z H z h n z X z --===∑ 由于FIR 滤波器的单位脉冲响应h(n)是一个有限长序列,H (z)是的(N −1)次多项式,它在Z 平面上有(N −1)个零点,同时在原点有(N −1)阶重极点。
因此,H(z)永远稳定。
FIR 滤波器设计的任务是选择有限长度的h(n),使传输函数满足一定的幅度特性和线性相位要求。
由于FIR 滤波器很容易实现严格的线性相位,所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。
(二) FIR 滤波器的窗函数设计原理 窗函数法的设计思想是按照所要求的理想滤波器频率响应,设计一个FIR 滤波器,使之频率响应来逼近。
先由的傅里叶反变换导出理想滤波器的冲激响应序列,即:由于是矩形频率特性,所以是一无限长的序列,且是非因果的,而要计的FIR 滤波器的冲激响应序列是有限长的,所以要用有限长的序列h(n)来逼近无限长的序列,最有效的方法是截断,或者说用一个有限长度的窗口函数w(n)序列来截取,即: 。
按照复卷积公式,在时域中的乘积关系可表示成在频域中的周期性卷积关系,即可得所设计的FIR 滤波器的频率响应:其中,为截断窗函数的频率特性。
由此可见,实际的FIR 数字滤波器的频率响应逼近理想滤波器频率响应的好坏,完全取决于窗函数的频率特性。
如果w(n)具有下列形式:⎩⎨⎧<≤≥<=N n N n n n w 0,1,0,0)( w(n)相当于一个矩形,我们称之为矩形窗。
即我们可采用矩形窗函数w(n)将无限脉冲响应截取一段来近似为。
经过加矩形窗后所得的滤波器实际频率响应能否很好地逼近理想频率响应呢?下图给出了理想滤波器加矩形窗后的情况。
理想低通滤波器的频率响应如图中左上角图,矩形窗的频率响应为左下角图。
根据卷积定理,即得实际滤波器的频率响应图形为图中右图。
由图可看出,加矩形窗后使实际频率响应偏离理想频率响应,主要影响有三个方面:(1)理想幅频特性陡直边缘处形成过渡带,过渡带宽取决于矩形窗函数频率响应的主瓣宽度。
(2)过渡带两侧形成肩峰和波纹,这是矩形窗函数频率响应的旁瓣引起的,旁瓣相对值越大,旁瓣越多,波纹越多。
(3)随窗函数宽度N 的增大,矩形窗函数频率响应的主瓣宽度减小,但不改变旁瓣的相对值。
为了改善滤波器的性能,需使窗函数谱满足:主瓣尽可能窄,以使设计出来的滤波器有较陡的过渡带;第一副瓣面积相对主瓣面积尽可能小,即能量尽可能集中在主瓣,外泄少,使设计出来的滤波器的肩峰和余振小逼近于理想滤波器。
但是这两个条件是相互矛盾的,实际应用中,折衷处理,兼顾各项指标。
(三) MATLAB 信号处理中提供的窗函数上边只考虑了矩形窗,如果我们使窗的主瓣宽度尽可能地窄,旁瓣尽可能地小,可以获得性能更好的滤波器,通过改变窗的形状来达到这个目的。
在数字信号处理的发展过程中形成了不同于矩形窗的很多窗函数,这些窗函数在主瓣和旁瓣特性方面各有特点,可满足不同的要求。
为此,用窗函数法设计FIR 数字滤波器时,要根据给定的滤波器性能指标选择窗口宽度N 和窗函数w(n)。
下面具体介绍几类类窗函数及其特性。
1. 矩形窗矩形窗函数的时域形式可以表示为:1,01()()0,N n N w n R n ≤≤-⎧==⎨⎩其他 它的频域特性为:()1j j 2sin 2e e sin 2N N W ωωωω-⎛⎫- ⎪⎝⎭⎛⎫ ⎪⎝⎭=⎛⎫ ⎪⎝⎭2. 汉宁窗函数汉宁窗函数的时域形式可以表示为:⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛+-=1π2cos 15.0)(n k k w N k ,,2,1 = 它的频域特性为:()()⎪⎭⎫ ⎝⎛--⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-++⎪⎭⎫ ⎝⎛--+=21j e 1π21π225.05.0N R R R N W N W W W ωωωωω其中,)(ωR W 为矩形窗函数的幅度频率特性函数。
汉宁窗函数的最大旁瓣值比主瓣值低31dB ,但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍,为8π/N 。
3. 海明窗函数海明窗函数的时域形式可以表示为:⎪⎭⎫ ⎝⎛--=1π2cos 46.054.0)(N k k w N k ,,2,1 = 它的频域特性为:⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-++⎪⎭⎫ ⎝⎛--+=1π21π223.0)(54.0)(N W N W W W R R R ωωωω 其中,()R W ω为矩形窗函数的幅度频率特性函数。
海明窗函数的最大旁瓣值比主瓣值低41dB ,但它和汉宁窗函数的主瓣宽度是一样大的。
4. 三角窗函数 三角窗是最简单的频谱函数为非负的一种窗函数。
三角窗函数的时域形式可以表示为:当n 为奇数时: ⎪⎩⎪⎨⎧≤≤+++-+≤≤+=n k n n k n n k n k k w 21,1)1(2211,12)( 当n 为偶数时:⎪⎩⎪⎨⎧≤≤+-≤≤-=n k n n k n n k n k k w 2,)1(221,12)( 它的频域特性为:()()221j j 2sin 41sin 12e e ⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛--=⎪⎭⎫ ⎝⎛--ωωωωN N W N R三角窗函数的主瓣宽度为8π/N ,比矩形窗函数的主瓣宽度增加了一倍,但是它的旁瓣宽度却小得多。
(四) 各窗函数的图形及幅频特性1. 用MATLAB 编程绘制各种窗函数的形状,窗函数的长度为21。
程序如下:clf;Nwin=21;n=0:Nwin-1;figure(1)for ii=1:4switch iicase 1w=boxcar(Nwin); %矩形窗stext='矩形窗';case 2w=hanning(Nwin); %汉宁窗stext='汉宁窗';case 3w=hamming(Nwin); %海明窗stext='海明窗';case 4w=triang(Nwin); %三角窗stext='三角窗';endposplot=['2,2,' int2str(ii)];subplot(posplot);stem(n,w); %绘出窗函数hold onplot (n,w,'r'); %绘出包络线xlabel('n');ylabel('w(n)');title(stext);hold off;grid on;end程序运行结果如下图:2.用MATLAB编程,采用512个频率点绘制各窗函数的幅频特性。
程序如下:clf;Nf=512; %窗函数复数频率特性的数据点数Nwin=20; %窗函数数据长度figure(1)for ii=1:4switch iicase 1w=boxcar(Nwin); %矩形窗stext='矩形窗';case 2w=hanning(Nwin); %汉宁窗stext='汉宁窗';case 3w=hamming(Nwin); %海明窗stext='海明窗';case 4w=triang(Nwin); %三角窗stext='三角窗';end[y,f]=freqz(w,1,Nf); %求解窗函数的幅频特性mag=abs(y); %求得窗函数幅频特性posplot=['2,2,' int2str(ii)];subplot(posplot);plot(f/pi,20* log10(mag/max(mag))); %绘制窗函数的幅频特性xlabel('归一化频率');ylabel('振幅/dB');title(stext);grid on;end程序运行结果如下图:由上图各窗函数的幅频特性可以看到各种窗函数都具有明显的主瓣和旁瓣。