FIR低通滤波器的设计 数字信号处理课程设计 毕业设计.doc

PAGE / NUMPAGES 目录 TOC \o 1-3 \h \z \u 1 技术要求 1 2 基本原理 1 2.1 FIR滤波器的设计 1 窗函数设计法 2 2.3 各种窗函数 4 3 滤波器的实现 4 3.1 Simulink模块实现滤波器功能 4 3.2 窗函数法编程实现滤波器功能 4 3.2.1 窗函数的实现 4 3.2.2 程序流程图 4 3.2.3 滤波器仿真结果 4 4 源程序代码 4 5 结论 4 6 心得体会 4 7 思考题 4 8 参考文献 4 附录 4 附录A Matlab相关函数及参数说明 4 附录B 调试分析 4 FIR低通滤波器的设计 1 技术要求 用窗函数法设计线性相位FIR低通滤波器。要求通带截止频率ωc=π/4,单位脉冲响应hn的长度N=23。绘出hn及其幅频响应特性曲线。 2 基本原理 随着信息时代和数字世界的到来,数字信号处理已成为当今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像,自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中,滤波占有极其重要的地位。数字滤波器是谱分析、雷达信号处理、通信信号处理应用中的基本处理算法,在系统设计中,滤波器的好坏将直接影响系统的性能。 现代数字滤波器可以用软件或设计专用的数字处理硬件两种方式来实现,用软件来实现数字滤波器优点是随着滤波器参数的改变,很容易改变滤波器的性能。 数字滤波器从功能上看,可分为低通、高通、带通、带阻滤波器。数字滤波器根据其单位脉冲响应可分为IIRInfinite Impulse Response无限长冲激响应滤波器和FIRFinite Impulse Response有限长冲激响应滤波器两类。IIR滤波器可以用较少的阶数获得很高的选择特性,但在有限精度的运算中可能出现不稳定现象,而且相位特性不好控制。在许多实际应用中为了保证滤波后的信号不产生相位失真。一般均采用FIR滤波器。 2.1 FIR滤波器的设计 FIR 滤波器能够在保证幅度特性满足技术要求的同时,易做成严格的线性相位特性,且FIR滤波器的单位抽样响应是有限长的,因而滤波器一定是稳定的,而且可以用快速傅里叶变换算法实现,大大提高了运算速率。同时只要经过一定的延时,任何非因果有限长序列都能变成因果的有限长序列,所以系统总能用因果系统来实现。但FIR必须用很长的冲激响应滤波器才能很好地逼近锐截止的滤波器,需要很大的运算量,要取得很好的衰减特性,需要较高的阶次。 有限长单位冲激响应〔FIR〕滤波器有以下特点： 1 系统的单位冲激响应h n在有限个n值处不为零； 2 系统函数Hz在|z|0处收敛,极点全部在z = 0处〔因果系统〕； 3 结构上主要是非递归结构,没有输出到输入的反馈,但有些结构中〔例如频率抽样结构〕也包含有反馈的递归部分。 长度为N的FIR滤波器的单位冲激响应hn的系统函数为： H(z)= 其差分方程为： y FIR滤波器直接型结构如图1所示： Z-1Z-1xn Z-1 Z-1 xn Z-1 Z-1 hN-1h1 hN-1 h1 h0 hN yn yn 图1 FIR滤波器直接型结构 FIR滤波器的设计任务是选择有限长的单位冲激响应,使传输函数满足技术要求。主要设计方法有窗函数法、频率取样法和等波纹逼近法等3种。 窗函数设计法 基于窗函数的FIR数字滤波器的设计方法通常也称之为傅里叶级数法,是用一定宽度窗函数街区无限脉冲响应序列,获得有限长的脉冲响应序列,从而得到FIR滤波器。它是在时域进行的,有理想滤波器的频率响应Hdejw推导出其单位冲激响应hdn,再设计一个FIR数字滤波器的单位冲激响应hn去逼近hdn,表示为： h 由此得到的离散滤波器的系统传递函数Hdz为： H 该hdn为无限长序列,因此Hdz是物理不可实现的。为了是系统变为物理可实现的,且使实际的FIR滤波器尽可能逼近理想滤波器的频率响应,用一个有限长度的窗函数将无限脉冲响应hdn截取一断hn来近似表示hdn,可得：hn=hdnwn,从而有： H 式中N表示窗口长度,这样Hz就是物理可实现的系统。并且从想形象为FIR滤波器的充要条件可知,为了获得线性相位FIR数字滤波器的冲激响应hn,则序列hn应该有τ=(N- 加窗处理对理想矩形频率响应产生了以下几点影响： 〔1〕加窗处理使理想频率特性在不连续点外边沿加宽,形成一个过渡带,过渡带的宽度等于窗的频率响应WRω的主瓣宽度Δω=4π/N。注意,这里所指的过渡带是两个肩峰之间的宽度,与滤波器真正的过渡带不同。 〔2〕在截止频率ωC两边ω=ωC±2π/N的地方〔即过渡带两边〕,Hω出现最大的肩峰值,肩峰的两侧形成起伏振荡,其振荡幅度取决于旁瓣的相对幅度,而振荡的多少,取决于旁瓣的多少。 〔3〕增加截取