利用窗函数法设计线性相位FIR滤波器.doc

《利用窗函数法设计线性相位FIR滤波器.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西安工程大学课程设计说明书题目：利用窗函数法设计线性相位FIR滤波器学院（系）：计算机科学学院年级专业：电子信息科学与技术01班学号：******学生姓名：***日期：20131231-201413论文摘要数字信号处理学科的一项重大进展是关于数字滤波器设计方法的研究。而FIR数字滤波器可以方便地实现线性相位且其群时延不随频率变化的，在任何幅度特性下都能稳定的工作，因此在数字信号处理中占有非常重要的地位。在现代电子系统中，FIR数字滤波器以其良好的线性特性被广泛使用。FIR数字滤波器传统的设计方法有窗函数法、频率抽样法和等波纹逼近法。用窗

2、函数设计FIR数字滤波器就是用有限长的脉冲相应逼近序列，其基本设计思想为：首先选定一个理想的选频滤波器，然后截取它的脉冲响应得到线性相位。本文是以窗函数设计方法为基础的。首先介绍了满足线性相位条件的FIR数字滤波器的基础知识，简述了FIR数字滤波器的设计方法与流程。本课题主要应用MATLAB软件利用窗函数法设计满足线性相位的FIR数字滤波器，并对所设计的滤波器进行仿真，得到滤波器的单位脉冲响应、频率响应特性。关键字：数字滤波器;　C语言。窗口函数设计FIR滤波器论文㈠、窗函数法的基本思想窗函数设计的基本思想是要选取某一种合适的理想频率选择性滤波器，

3、然后将它的脉冲响应截断以得到一个线性相位和因果的FIR滤波器。因此这种方法的重点在于选择某种合适的窗函数和一种理想滤波器。对于给定的滤波器技术指标，选择滤波器长度和具有最窄主瓣宽度和尽可能小的旁瓣衰减的某个窗函数。任何数字滤波器的频率响应都是ω的周期函数，它的傅立叶级数展开式为：(1-1)式中，其中的为滤波器的归一化的截止频率。傅立叶系数ha(n)实际上就是理想数字滤波器的冲激响应。获得有限冲激响应数字滤波器的一种可能方法就是把无穷级数截取为有限项级数来近似，而吉布斯(Gibbs)现象使得直接截取法不甚令人满意。窗函数法就是用被称为窗函数的有限加权

4、系列{ω(n)}来修正式(1-2)的傅立叶级数，以求得要求的有限冲激响应序列h(n)，即有：(1-2)是有限长序列，当n>N-1及n<0时，。㈡、窗函数的选择工程中比较常用的窗函数有：矩形窗函数、三角形((Bartlett)窗函数、汉宁(Harming)窗函数、汉明(Hamming)窗函数、布莱克曼(Blackman)窗函数和凯塞(Kaiser)窗函数。这几种窗函数的比较见下表所示：窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/dB矩形窗-134π/N-12三角窗-258π/N-25汉宁窗-318π/N-44海明窗-418π/N-53布莱克曼窗-5

5、712π/N-74凯塞窗-5710π/N-80　六种窗函数的基本参数几种常用窗函数对比，窗函数的选择原则是：①具有较低的旁瓣幅度，尤其是第一旁瓣幅度；②旁瓣幅度下降速度要大，以利增加阻带衰减；③主瓣的宽度要窄，以获得较陡的过渡带。通常上述三点很难同时满足。当选用主瓣宽度较窄时，虽然得到较陡的过渡带，但通带和阻带的波动明显增加：当选用最小的旁瓣幅度时，虽能得到匀滑的幅度响应和较小的阻带波动，但过渡带加宽。因此，实际选用的窗函数往往是它们的折衷。在保证主瓣宽度达到一定要求的条件下，适当牺牲主瓣宽度来换取旁瓣波动的减少。就几种窗函数的比较在本课题中我们采

6、用了汉明窗。1、矩形窗2、汉宁窗3、三角窗4、海明窗5、布拉克曼窗6、凯塞窗其中是I。（β）第一类零阶修正贝塞耳函数。口是控制窗函数形状的参数，β越大，w（n）窗越窄，频谱的旁瓣越小，但主瓣也相应加宽。口的典型值为4≤p≤9。P=0时，凯塞窗变为矩形窗；β=5．44时，凯塞窗与海明窗接近；p=8.5时，凯塞窗与布拉克曼窗接近。㈣、FIR滤波器的设计方法1．选择窗函数的类型和长度（1）.根据阻带最小衰减选择窗函数的类型原则是：在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。（2）.根据过渡带的宽度选择窗函数的长度2．按性能指标要求，构造希望频

7、率响应函数(1-12)近似为过渡带中心频率，幅度函数衰减一半（-6dB)(1-13)（1）.确定期望滤波器的单位脉冲响应(1-14)（2）.加窗得到设计结果(1-15)㈤、利用窗函数设计C语言实现：1子函数语句voldfirwin(n，band，fln，fhn，wn，h)2．形参说明n—整型变量。滤波器的阶数。band—整型变量。滤波器的类型。取值为1、2、3和4，分别对应低通、高通、带通和带阻滤波器。fln——双精度实型变量。fhn——双精度实型变量。对于低通和高通滤渡器，fln：通带边界频率对于带通和带阻滤渡器，fln：通带下边界频率，fhn：

8、通带上边界频率wn一整型变量。窗函数的类型；取值l到7，分别对应矩形窗，图基窗，三角窗，汉宁窗，海明窗，布拉克曼窗和凯塞窗