基于MATLAB与FPGAFIR滤波器设计与仿真.docVIP

基于MATLAB与FPGAFIR滤波器设计与仿真 　　【摘要】数字滤波器是数字信号处理领域内的重要组成部分。通过原理分析和结构对比，本文结合MATLAB工具软件介绍了线性相位的FIR数字滤波器的设计方法，并在Xilinx的FPGA器件上完成设计实现。最后，使用MATLAB和ISE软件对数据进行了分析，证实了设计实现的正确性与可行性。 　　【关键词】FPGA；MATLAB；FIR；数字滤波器；ISE 　　数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机械振动信号、遥感遥测信号等。上述这些信号大部分是模拟信号，也有小部分是数字信号。数字滤波器是对数字信号进行滤波处理以得到期望的响应特性的离散时间系统。数字滤波器工作在数字域，它处理的对象是经由采样器件将模拟信号转换而得到的数字信号。它的设计方法有很多，比较常用的有窗函数法、频率抽样法和最佳一直逼近法，但这些设计方法在设计滤波器尤其是高阶滤波器时工作量很大，而利用MATLAB信号处理工具箱进行数字滤波器的设计工作，能够减少工作量，提高设计效率。文中使用MATLAB这一高效工具完成一种FIR滤波器的前期参数设计工作，随后给出了FPGA实现该FIR滤波器的设计方案，最后通过MATLAB和ISE的联合仿真，验证了该FIR滤波器满足设计要求。 　　1.FIR滤波器基本原理 　　根据冲激响应的时域特性，数字滤波器可分为无限长冲激响应滤波器（IIR）和有限长冲激响应滤波器（FIR）。相比于IIR滤波器，FIR的突出优点是：具有精确的线性相位；总是稳定的；硬件容易实现；滤波器的过渡过程具有有限区间。 　　FIR滤波器的基本结构可以理解为一个分节的延时线，把每一节的输出加权累加，可得到滤波器的输出。FIR滤波器的冲激响应h（n）是有限长的，数学上N阶FIR滤波器可以表示为： 　　2.FIR滤波器的结构形式 　　FIR数字滤波器的实现一般有直接型、转置型、线性相位型、级联和频率采样等。 　　（1）直接型FIR数字滤波器 　　直接型FIR数字滤波器也称卷积型或横截型FIR数字滤波器。因为差分方程是信号的卷积形式，所以称为卷积型FIR数字滤波器。因为FIR数字滤波器的输入是一条x（n）延时链的横向结构，所以称为横截型FIR数字滤波器。图1给出了一个N阶的直接型FIR数字滤波器的图解。 　　（2）转置型FIR数字滤波器 　　直接型FIR数字滤波器的一种变换可称之为转置型FIR数字滤波器，可以根据图1来构造：l）互换滤波器的输入和输出；2）颠倒信号流向的方向；3）用差分放大器代替加法器，反之亦可。转置型FIR数字滤波器，这种FIR数字滤波器有两个特点：一是所有的乘法器的被乘数都是相同的，二是无需额外寄存器来缓存采样的数据。 　　（3）线性相位的FIR数字滤波器 　　FIR数字滤波器系统的最主要特性就是它能够构成具有线性相位的数字滤波器。所谓的线性相位特性就是指数字滤波器对不同频率的输入信号所产生的相移与输入信号的频率是直线关系。因此，在FIR数字滤波器通带内的信号通过FIR数字滤波器后，除了由相频特性的斜率所造成的延迟之外，可以不失真的保留通带内的全部信号。线性相位的FIR数字滤波器的对称性还可以降低滤波器中所需的乘法器的个数。如图2所示，可以看到在“对称”的结构中，每一个滤波周期均提供了一个乘法器的预算资源，使得乘法器的数量是图1中（直接型FIR数字滤波器）的一半。 　　（4）级联FIR数字滤波 　　级联FIR数字滤波器在需要控制FIR数字滤波器的系统函数的零点时，可将式（3）分解成二阶的实系数因子的模式： 　　这样就可以得到二阶级联结构。这种结构中，每一节都控制着一对一的零点，因此在需要控制系统的传输零点时可采用。但相应的FIR数字滤波器系数会增加，乘法运算的次数会增加，而且需要乘法器的个数更多，因此还会需要更多的存储器，运算时间也会比直接型FIR数字滤波器增加。 　　（5）频率采样FIR数字滤波器 　　频率采样FIR数字滤波器结构，是用系数将FIR数字滤波器参数化的实现的结构。一个有限长的序列可以通过相同长度的频域采样值而唯一确定。 　　将系统函数在单位圆上N等分后取样，是单位取样响应h（n）的离散傅里叶变换H（k）。H（k）跟系统函数的关系可以用内插公式来表示： 　　频率采样FIR数字滤波器的优点在于其选频性好，适用于窄带滤波；不同的FIR数字滤波器，若其长度相同，则可通过改变系统的系数后，用同一个网络实现，复用性好。频率采样FIR数字滤波器不足之处在于在具体实现时难免存在误差，零、极点可能不能够正好抵消，造成系统的不稳