MATLAB在FIR数字带通滤波器设计中的应用.docVIP

MA TLA B 在FI R 郝 娟徐沛文 (华中师范大学,湖北武汉 [摘 要] 的参数,提出了利用MATLAB 中的FDATool 工具箱和Simulink [关键词] MATLAB ;FDATool ;Simulink ;数字滤波器 —————————————————————————— —作者简介:郝娟,女,山西忻州人,本科,研究方向:信号处理,通信工程。 1.引言 在通信与电子信息当中,在对信号作分析与处理时,常会遇到有用信号叠加无用噪声的问题。这些噪声信号有的是与信号同时产生的,有的是在传输过程中混入的,在接收的信号中,必须消除或减弱噪声干扰,这是信号处理中十分重要的问题。根据有用信号与噪声的不同特性,消除或减弱噪声,提取有用信号的过程就称为滤波。滤波器的种类很多,实现方法也多种多样。随着数字技术的飞速发展,数字滤波理论也得到了长足的进步。因此,对数字滤波系统硬件实现的要求也越来越高,而软件模拟的方法不仅能及时地提供系统运行的信息,还可以随时改变系统结构从而验证全新的系统,所以软件仿真变得非常实用与重要。利用MATLAB 信号处理箱可以快速有效地实现数字滤波器的设计与仿真。本设计重点是对FIR 滤波器分析与设计,并在此基础上对FIR 滤波实现理论进行研究。 2.FIR 滤波器的原理与设计2.1 滤波器概述 数字滤波在数字信号处理中占有重要的地位,是广泛使用的一种基本线性处理模块,它可以实现模拟器件很难达到的准确线性相位关系特性。数字滤波器分为无限冲激响应滤波器(IIR和有限冲激响应滤波器(FIR。由于FIR 滤波器在阶数相同的条件下运算速度比IIR 滤波器快,同时FIR 滤波器具有不含反馈环路、结构简单以及可以实现的严格线性相位等优点,因而在对相位要求比较严格的条件下,常常采用FIR 数字滤波器。 2.2 FIR 数字滤波器基本原理 FIR 数字滤波器的基本结构为一个分节的延时线,把每一节的输出加权累加,可得到滤波器的输出。FIR 滤波器的冲击响应是有限长的。FIR 数字滤波器的差分方程为:y(n=N-1 i =0Σh(ix(n -i。其中, N 是FIR 滤波器的抽头数,h(i是FIR 滤波器的第i 级抽头系数,x(n表示第n 时刻的输入样本,y (n为滤波器的输出样本。其系统函数为:H(z=N-1 i =0Σh(nz -n 。 可见FIR 系统只有零点,因此这类滤波器不像IIR 滤波器那 特性,自己突出的优点,其一是系统总是稳定的,其二是易实现线性相位,其三是允许设计多通带(或多阻带滤波器。 FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。目前常用的设计方法主要有窗函数法、频率取样法及等波纹逼近法。本文应用窗函数法设计FIR 数字滤波器。 2.3 窗函数法设计原理 窗函数法的核心是从给定的频率特性,通过加窗确定有限长单位脉冲响应序列h(n。加窗的作用是通过把理想滤波器的无限长脉冲响应h d (n乘以窗函数w(n来产生一个被截断的脉冲响应,即:h(n=h d (nw(n。工程中比较常用的窗函数有:矩形窗、三角窗、汉宁窗、汉明窗、布莱克曼(Blackman窗等。这几种窗函数参数的比较见表1所示。 表1几种常见窗函数参数对比 通常窗函数的这些条件很难同时满足,当选用主瓣宽度较窄时,虽然得到较陡的过渡带,但通带和阻带的波动明显增加;当选用最小的旁瓣幅度时,虽能得到匀滑的幅度响应和较小的阻带波动,但过渡带加宽。因此,实际选用的窗函数往往是它们的折中。在保证主瓣宽度达到一定要求的条件下,适当牺牲主瓣宽度来换取旁瓣波动的减少。 3.FIR 数字滤波器的MATLAB 设计3.1 FDATool 界面 滤波器设计与分析工具箱FDATool (Filter Design AnalysisTool是MATLAB 信号处理工具箱里专用的滤波器设计分析工具,它可以设计几乎所有基本的常规滤波器,包 窗函数旁瓣峰值主瓣最小阻带衰减/dB 宽度/π/M 衰减/dB 矩形窗13-4-21三角窗-258-25汉宁窗-318-44汉明窗-418-53布莱克曼窗 -57 12 -74 —49— 。 ,在 。 83.3Hz 图1FDATool带通FIR滤波器设计 本例中,首先在filter type中选择bandpass(带通滤波 器;在design method选项中选择FIR window(FIR滤波器 窗函数法,接着在window specifications选项中选取 hamming;指定filter order项中的specify order=95;由于采 用窗函数法设计,只要给出通带下限截止频率fc1和通带上 限截止频率fc2,选取fc1=70Hz,fc2=84Hz。设置完以后点 击design filter即