911制作网站,企业管理平台下载,建设公司网站多少钱,网页设计基础ppt相对于模拟滤波器#xff0c;数字滤波器具有高精度、高可靠性、可编程改变滤波特性、便于集成等一系列优点#xff0c;并且理论上可实现近似理想频率特性的滤波性能。经典的数字滤波器主要包括有限脉冲响应(Finite Impulse Response#xff0c;FIR)滤波器和无限脉冲响应(Inf…相对于模拟滤波器数字滤波器具有高精度、高可靠性、可编程改变滤波特性、便于集成等一系列优点并且理论上可实现近似理想频率特性的滤波性能。经典的数字滤波器主要包括有限脉冲响应(Finite Impulse ResponseFIR)滤波器和无限脉冲响应(Infinite Impulse ResponseIIR)滤波器。相对于IIR滤波器来说当FIR滤波器系数保持线性对称结构时能够在满足幅频响应要求的同时获得严格的线性相位特性。在数据通信、语音信号处理、雷达信号处理等领域通常要求信号在传输过程中不能有明显的相位失真因此FIR滤波器获得了更广泛的应用[1]。因其独特的可编程结构及强大的并行运算能力很适合于实现FIR滤波器设计。基于FPGA的FIR滤波器设计近年来一直是学者们研究的热点主要集中在最优滤波器设计方法、提高滤波器运行速度、滤波器方法等方面[2-7]。常规的系数量化方法仅关注对滤波器本身幅频特性的影响[8]较少考虑在滤波器输出位宽确定的情况下如何提高输出幅度进而最大化输出信噪比的问题。在研究分析传统滤波器系数量化方法的基础上提出了基于滤波器增益的系数量化方法给出了MATLAB及FPGA实现后的ModelSim仿真测试数据结果表明这种方法最多可将滤波输出幅度提高近一倍信噪比提高近6 dB。1 FIR滤波器原理及结构FIR滤波器是指单位脉冲响应长度有限的滤波器。FIR滤波器的突出特点是其单位取样响应h(n)是一个N点长的有限长序列0≤n≤N-1。滤波器的输出y(n)可表示为输入序列x(n)与单位取样响应h(n)的线性卷积。从系统函数很容易看出FIR滤波器只在原点上存在极点这使得FIR滤波器具有全局稳定性。FIR滤波器是由一个抽头延迟线加法器和乘法器的集合构成的每一个乘法器的操作系数就是一个FIR系数。因此FIR的这种结构也被人们称为抽头延迟线结构。图1是FIR滤波器FPGA实现的典型结构图。2 系数量化对FIR滤波器性能的影响2.1 FIR滤波器的MATLAB设计方法在采用FPGA编程实现FIR滤波器之前通常采用MATLAB软件设计滤波器系数。MATLAB提供了丰富的FIR滤波器函数分别对应不同的设计方法如窗函数法、频率取样法、基于“最大误差最小”的最优设计法等。为便于对比采用最优设计方法设计FIR低通滤波器(MATLAB设计函数为firpm)其参数为采样频率32 MHz过渡带0.5 MHz~1.5 MHz。图2为MATLAB设计出的滤波器幅频响应图。由图2可知FIR滤波器的长度N越大滤波器过渡带越接近设计参数通带及阻带滤波性能越好当N41时阻带衰减为26 dB当N61时阻带衰减为38.5 dB当N81时阻带衰减为48 dB当系数长度超过一定值后(本实例为61)滤波器过渡带的改善余量已很小但通带及阻带滤波性能会持续增加。不失一般性下文以长度为61的低通滤波器为例进行讨论。2.2 系数的量化方法及其对滤波性能的影响分析根据MATLAB的FIR滤波器设计原理无论滤波器长度多长滤波器通带内增益均为1(0 dB)滤波器系数为带小数的浮点数。由于FPGA只能处理二进制数据MATLAB设计出的滤波器系数需要转换成二进制数据才能完成FPGA电路实现。滤波器系数是带符号的数据FPGA实现时需要用二进制补码格式表示。设滤波器系数的量化位宽为Q则系数C的表示范围为-2Q-1≤C≤(2Q-1-1)又由于FIR滤波器系数的严格对称特性实际范围为(1-2Q-1)≤C≤(2Q-1-1)。传统的滤波器系数量化方法可分解为3个步骤(1)找出滤波器系数中最大的绝对值数M(2)以M为基准对所有滤波系数进行归一化处理即对所有系数除以M(3)对所有滤波器系数乘以(2Q-1-1)并对处理后的系数取整形成整数系数(4)将整数滤波器系数转换成二进制补码数据。图3为不同量化位宽情况下的滤波器幅频响应图。由图3(a)中的归一化幅频图可知滤波器系数量化位宽越小量化后的滤波性能越差(当量化位宽为6 bit时阻带衰减相比量化前减小了约6 dB)反之量化位宽越大量化前后的滤波性能相差越小。当量化位宽超过一定值(本实例为12 bit)后滤波性能几乎不再变化。不失一般性下文以长度为61、量化位宽为12 bit的低通滤波器为例进行讨论。3 基于增益的FIR滤波器系数量化方法3.1 一般量化方法对系统性能的影响根据前面讨论的滤波器量化方法系数量化过程相当于滤波器增益扩大了(2Q-1-1)/M倍。从图3(b)中的未归一化幅频图可知量化位宽越大滤波器通带内的增益越大。当量化位宽为6 bit时通带增益为54 dB当量化位宽为12 bit时增益为90 dB。FPGA设计FIR滤波器一般采用定点二进制数据进行运算。前面讨论系数量化对系统性能的影响时实际有一个前提条件即滤波运算采用全精度运算滤波运算的结果保留全部有效数据位。当滤波器系数量化位宽为Q时由于滤波器增益扩大了(2Q-1-1)/M倍(一般来说M≤1)为保证全精度运算保留全部有效滤波输出位宽相对于滤波器输入信号而言位宽将增加Qabit。其中对于前面讨论的长度N61量化位宽Q12的低通滤波器而言M0.062 1(2Q-1-1)/M32 963则滤波输出位宽需增加Qa16 bit。假设输入数据位宽为10 bit则滤波器全精度输出位宽为26 bit。3.2 基于增益的系数量化方法根据前面的讨论量化滤波器系数后输出有效位宽所能表示信号的信噪比下降是由于滤波输出信号达不到满量程状态。为此优化滤波器系数量化方法采用基于滤波器增益的方法对系数进行量化。结合前述的滤波器实例具体步骤为(1)对所有滤波器系数乘以21665 536并对处理后的系数取整形成整数系数(2)将整数滤波器系数转换成二进制补码数据。经过改进的量化处理后相当于对滤波增益直接增加了65 536倍(常规方法仅增加了32 963倍)滤波输出的数据幅值能够达到输出数据有效位宽所能表示的满量程值实现了输出位宽一定的条件下输出信号信噪比的最大化。经过这样的量化处理系数的量化位宽为对于上文所讨论的低通滤波器量化位宽需要13 bit。图4为常规量化方法和改进量化方法的滤波器系数(冲激响应)图其中Qa16。从图4可以看出在滤波输出位宽相同的条件下改进的量化方法所产生的滤波器系数(冲激响应)幅值大于常规方法相当于滤波器增益更大。同时也要注意到改进的量化方法所需的滤波器系数位宽由12 bit增加到13 bit。也就是说增加滤波器增益是以增加系数量化位宽为代价的。4 FPGA实现后的测试结果分析接下来采用Xilinx的Spartan6系列FPGA器件XC6SLX25-FT256C为目标平台调用ISE14.7工具中的FIR滤波器核实现前文所讨论的低通滤波器(采样频率为32 MHz过渡带为0.5 MHz~1.5 MHz滤波器长度N61)。仿真工具为ModelSim10.1输入信号为2路单频信号的叠加。两路单频信号频率分别为0.22 MHz和2 MHz功率比为1:1。输入位宽为10 bit滤波器输出位宽为26(Qa16)。图5、图6分别为常规量化方法和改进量化方法的滤波器ModelSim仿真波形图。从图5、图6中可以看出两种方法都能够有效滤除带外干扰信号(2 MHz单频信号)改进量化方法的滤波器输出信号幅度值(16641407)明显大于常规量化方法的幅度值(7910390)相当于增加了滤波增益增加了滤波输出信号的信噪比。5 结论本文详细分析了滤波器系数量化方法对滤波性能的影响。在对常规系数量化方法的步骤、原理进行分析的基础上针对滤波输出数据幅值达不到有效位宽满量程的现象提出了基于滤波器增益的系数量化改进方法。MATLAB仿真及FPGA实现后的ModelSim仿真均表明在滤波输出位宽相同的条件下改进的量化方法可以有效提高滤波器增益进而提高滤波输出信号的信噪比。需要注意的是增加滤波器增益是以增加系数量化位宽、增加FPGA实现后的逻辑资源为代价的。参考文献[1] 高耀鸿.基于FPGA的FIR低通滤波器[D].长沙湖南大学2012.[2] 朱效效蔡俊陆伟.基于优化DA算法滤波器的设计及其FPGA实现[J].电子技术应用201541(2)59-6064.[3] 徐彦凯双凯.提高定点精度的FPGA信号处理算法[J].计算机工程与科学201234(9)197-200.[4] 郭广浩刘志哲孟庆龙等.多级并行流水FIR数字滤波器的设计与验证[J].现代电子技术201538(1)69-72.[5] 李凯勇.基于FPGA高效FIR滤波器设计[J].青海大学学报201735(6)56-60.[6] 张娜李春祎.多相并行FIR滤波器的FPGA高速实现方法[J].无线电通信技术201743(4)86-90.[7] 周龙王红玲崔新莹等.基于FPGA及分布式算法滤波器[J].电子设计工程201725(23)139-142.[8] 曹振吉何敏.基于FPGA和Matlab的FIR数字滤波器[J].现代电子技术201538(19)98-101.作者信息:杜 勇(四川工商学院 电子信息工程学院四川 成都611745)
