窗函数法设计FIR数字低通滤波器-课程设计初稿(共18页).doc
《窗函数法设计FIR数字低通滤波器-课程设计初稿(共18页).doc》由会员分享,可在线阅读,更多相关《窗函数法设计FIR数字低通滤波器-课程设计初稿(共18页).doc(18页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上数字信号处理课程设计报告 题 目 窗函数法设计FIR数字低通滤波器学 院 信息工程学院 专 业 通信工程 班 级 学 号 学生姓名 指导教师 二 0一二年十二月引言数字化是控制系统的重要发展方向,而数字信号处理已在通信、语音、图像、自动控制、雷达、军事、航空航天等领域广泛应用。数字信号处理方法通常涉及变换、滤波、频谱分析、编码解码等处理。数字滤波是重要环节,它能满足滤波器对幅度和相位特性的严格要求,克服模拟滤波器所无法解决的电压和温度漂移以及噪声等问题。而有限冲激响应FIR滤波器在设计任意幅频特性的同时能够保证严格的线性相位特性。利用FPGA可以重复配置高精度的FIR
2、滤波器,使用VHDL硬件描述语言改变滤波器的系数和阶数,并能实现大量的卷积运算算法。结合MATLAB工具软件的辅助设计,使得FIR滤波器具有快速、灵活、适用性强,硬件资源耗费少等特点。FIR滤波器是最常用的组件之一,它完成信号预调、频带选择和滤波等功能。FIR滤波器在截止频率的边沿陡峭性能虽然不及IIR滤波器,但是,考虑到FIR滤波器严格的线性相位特性和不像IIR滤波器存在稳定性的问题,FIR滤波器能够在数字信号处理领域得到广泛的应用。FIR是有限冲激响应(Finite Impulse Response)的简称。由线性系统理论可知,在某种适度条件下,输入到线性系统的一个冲击完全可以表征系统。当
3、我们处理有限的离散数据时,线形系统的响应(包括对冲击的响应)也是有限的。若线性系统仅是一个空间滤波器,则通过简单地观察它对冲击的响应,我们就可以完全确定该滤波器。通过这种方式确定的滤波器称为有限冲击响应(FIR)滤波器。FIR滤波器是在数字信号处理(DSP)中经常使用的两种基本的滤波器之一。 FIR 滤波器具有严格的相位特性,对于信号处理和数据传输是很重要的。目前 FIR滤波器的设计方法主要有三种:窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。用窗函数设FIR 滤波器的基本思路:从时域出发设计 h(n)逼近理想 hd(n)。设理
4、想滤波器的单位响应在时域表达为hd(n),则Hd(n) 一般是无限长的,且是非因果的,不能直接作为FIR 滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器单位抽样响应 h(n)最直接的方法是先将hd(n)往右平移,再进行截断,即截取为有限长因果序列:h(n)=hd(n)w(n),并用合适的窗函数进行加权作为 FIR 滤波器的单位脉冲响应。按照线性相位滤波器的要求,线性相位FIR数字低通滤波器的单位抽样响应h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数,即用矩形窗设计的FIR 低通滤波器,所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象,且最大波纹大约为幅度的9%,这个现象称为
5、吉布斯(Gibbs)效应。为了消除吉布斯效应,一般采用其他类型的窗函数。MATLAB 设计 FIR 滤波器有多种方法和对应的函数。窗函数设计法不仅在数字滤波器的设计中占有重要的地位,同时可以用于功率谱的估计,从根本上讲,使用窗函数的目的就是消除由无限序列的截短而引起的Gibbs现象所 带来的影响。目 录一、课程设计的目的和意义 1、通过课程设计把自己在大学中所学的知识应用到实践当中。 2、深入了解利用Matlab设计FIR数字滤波器的基本方法。3、在课程设计的过程中掌握程序编译及软件设计的基本方法。4、提高自己对于新知识的学习能力及进行实际操作的能力。5、锻炼自己通过网络及各种资料解决实际问题
6、的能力。二、课程设计题目描述和要求关于窗函数法设计FIR数字低通滤波器题目:对模拟信号进行低通滤波处理,要求通带0f2kHZ内衰减小于1dB,阻带3kHZf上衰减大于30dB。希望对模拟信号采样后用线性相位FIR数字滤波器实现上述滤波,采样频率Fs=10kHZ。要求: 用窗函数法设计满足要求的FIR数字低通滤波器,为了降低运算量,滤波器阶数尽量低。三、设计原理1.基本原理设计低通FIR数字滤波器,寻求一系统函数,使其频率响应逼近滤波器要求的理想频率响应,其对应的单位脉冲响应 (1)其中 如果所希望的滤波器的理想的频率响应函数为,则其对应的单位脉冲响应为 (2)窗函数设计法的基本原理是用有限长单
7、位脉冲响应序列逼近。由于往往是无限长序列,而且是非因果的,所以用窗函数将截断,并进行加权处理,得到: (3)就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数为 (4)式中,N为所选窗函数的长度。2.典型的窗函数(1)矩形窗(Rectangle Window) (5) 其频率响应和幅度响应分别为:, (6)(2)三角形窗(Bartlett Window) (7)其频率响应为: (3)汉宁(Hanning)窗,又称升余弦窗 (8)其频率响应和幅度响应分别为:(4)汉明(Hamming)窗,又称改进的升余弦窗 (9)其幅度响应为: (5)布莱克曼(Blankman)窗,又称二阶升余
8、弦窗 (10)其幅度响应 : (6)凯塞(Kaiser)窗 (11) 其中:是一个可选参数,用来选择主瓣宽度和旁瓣衰减之间的交换关系,一般说来,越大,过渡带越宽,阻带越小衰减也越大。I0()是第一类修正零阶贝塞尔函数。 若阻带最小衰减表示为,的确定可采用下述经验公式: 若滤波器通带和阻带波纹相等即p=s时,滤波器节数可通过下式确定: 式中:我们知道,用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见下表1表1 各种窗函数的基本参数窗函数旁瓣峰值幅
9、度/dB过渡带宽阻带最小衰减/dB矩形窗-134/N-12三角形窗-258/N-25汉宁窗-318/N-44哈明窗-418/N-53不莱克曼窗-5712/N-74凯塞窗(=7.865)-5710/N-80这样选定窗函数类型和长度N之后,求出单位脉冲响应,并按照式(4)求出。是否满足要求,要进行演算。一般在尾部加零使长度满足2的整数次幂,以便用FFT计算。如果要观察细节,补零点数增多即可。如果不满足要求,则要重新选择窗函数类型和长度N ,再次验算,直至满足要求。如果要求线性相位特性,则还必须满足 (12)根据上式中的正、负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据所设计的滤波特性
10、正确选择其中一类,例如,要设计线性相位低通特性,可以选择这一类,而不能选择这一类。我们在设计滤波器时,希望窗谱主瓣尽可能窄,以获得较陡的过渡带,同时尽可能减小最大旁瓣的相对幅度,我们发现这是不可能的,因为这本身就是一个矛盾体,所以在设计滤波器时只是根据实际情况来选择合适的窗函数主程序框图如图下图所示。其中幅度特性要求用dB表示。开始读入窗口长度N计算hd(n)调用窗函数子程序求w(n)调用子程序(函数)计算H(k)=DFTh(n)调用绘图子程序(函数)绘制H(k)幅度相位曲线结束主程序框图计算h(n)= hd(n) w(n)四、FIR 滤波器的窗函数设计基本方法1. FIR滤波器简介 FIR数
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 函数 设计 FIR 数字 滤波器 课程设计 初稿 18
限制150内