fir数字滤波器分布式算法的原理及fpga实现

《fir数字滤波器分布式算法的原理及fpga实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、FIR数字滤波器分布式算法的原理及FPGA实现摘要：在利用FPGA实现数字信号处理方面，分布式算法发挥着关键作用，与传统的乘积-积结构相比，具有并行处理的高效性特点。详细研究了基于FPGA、采用分布式算法实现FIR数字滤波器的原理和方法，并通过XilinxISE在Modelsim下进行了仿真。   关键词：分布式算法DALUTFPGAFIR数字滤波器正在迅速地代替传统的由R、L、C元件和运算放大器组成的模块滤波器并且日益成为DSP的一种主要处理环节。FPGA也在逐渐取代ASIC和PDSP，用作前端数字信号处理的运算（如：FIR滤波、CORDIC算法或FFT）。乘累加运算是实现大多数D

2、SP算法的重要途径，而分布式算法则能够大大提高乘累加运算的效能。1传统的乘累加结构（MACMultiplyingAccumulator）FIR数字滤波器基本理论FIR滤波器被称为有限长脉冲响应滤波器，与IIR数字滤波器相对应，它的单位脉冲响应h(n)只有有限个数据点。输入信号经过线性时不变系系统输出的过程是一个输入信号与单位脉冲响应进行线性卷积的过程，即：式中，x(n)是输入信号，y(n)是卷积输出，h(n)是系统的单位脉冲响应。可以看出，每次采样y(n)需要进行L次乘法和L-1次加法操作实现乘累加之和，其中L是滤波器单位脉冲响应h(n)的长度。可以发现，当L很大时，每计算一个点，则

3、需要很长的延迟时间。2乘累加运算的位宽分配DSP算法最主要的就是进行乘累加运算。假设采样信号的位宽用N来表示，则N位与N位的乘累结果需要2N位的寄存器来保存；如果两个操作数都是有符号数，则乘积只有2N-1个有效位，因为产生了两个符号位。为了使累加器的结果不产生溢出(计算结果超出了范围，产生了一个进位)，需要对累加器进行冗余设计，也就是说要在累加器2N的位宽上多设计出K位，累加器的长度M计算方式如下（L为滤波器的长度）：对于无符号数：M=2N+K=2N+log2L对于有符号数：M=2N＋K=2N+log2L-13乘累加运算的分布式算法原理分析得益于XilinxFPGA查找表结构的潜能，

4、分布式算法在滤波器设计方面显示出了很高的效率，自20世纪90年代初以来越来越受到人们的重要。分布式算法是基于查找表的一种计算方法，在利用FPGA实现数字信号处理方面发挥着重要的作用，可以大大提高信号的处理效率。它主要应用于数字滤波、频率转换等数字信号处理的乘累加运算。分布式算法推导如下：设Ak是已知常数（如滤波器系数Taps、FFT中的正弦/余弦基本函数等），Xk(n)是变量，可以看作是n时刻的第k个采样输入数据，y(n)代表n时刻的系统响应。那么它们的内积为：其中，Xk(n)变量可以写成下面的格式：式中，B为数据格式的字长,Xkb是变量的二进制位，只有“0”和“1”两种状态。将（2

5、）式代入（1）式得：（这种方式把MAC运算转化为根据输入信号x寻址ROM的加法运算，需要2^N的ROM空间,N为系数的个数，和2^n相乘只需左右移位，系数Ak如果是负数，就使用补码表示，不对，减法怎么做？）4FPGA实现过程中查找表的构造方法根据以上论述，括号中的每一乘积项代表着输入变量的某一位与常量的二进制“与”操作，加号代表着算术和操作，指数因子对括号中的值加权。如果事先构造一个查找表，该表存储着括号中所有可能的组合值，就可以通过所有输入变量相对应位的组合向量（XNb,X(N-1)b,...x1b）对该表进行寻址，该查找表称为DALUT。DALUT的构造规则如表1所示。5采用分布

6、式算法实现FIR数字滤波器为了说明问题，以一个三个系数的FIR数字滤波器为例设计分布式算法，字宽也设置为三位。设FIR数字滤波器系数为：h(0)=5,h(1)=2,h(2)=3。在进行FPGA设计时，该表以组件Component形式构建，设置为ROM结构，提供输入寻址端口table_in[2..0],输出端口table_out[3..0]。FPGA算法的结构图如图2所示。   算法实现中的几个关键问题为：（1）采用状态机实现分布式算法的状态转移状态机的实现如图3所示，设置三个状态s0、s1、s2。状态s0完成数据的装入，数据寄存器需要成对出现，一个完成数据的延迟，另一个完成数据的移位

7、，并将状态转移到s1；状态s1完成查找表功能、数据移位和分布式算法的乘累加运算，数据移位一个数据宽带后将状态转移到s2;状态s2完成数据的输出，并将状态转移到s0。利用状态机可以条理清楚地简化计算过程，在算法实现时发挥着关键的作用。   （2）系统时钟与数据输入时钟的关系根据上述的状态转移关系，可以得出：每输入一个数据，在下一次数据输入之前，需要在状态s1停留一个数据宽带（三位）的时钟时间，在s2停留一个时钟的数据输出时间。也就是说，系统时钟频率应是数据输