一种分段卷积快速算法设计与实现.docVIP

一种分段卷积快速算法的设计与实现 　　摘 要：常用的线性卷积方法要求两个输入序列的持续时间相同，但在实际工程中经常会遇到某个输入序列具有较长持续时间的情况，从而无法达到信号“实时”处理的要求。在这种情况下，分段卷积是一种有效的解决方案。设计了一种分段卷积快速算法模块，在FPGA中采用流水线结构进行实时处理。经检验该方法正确且能很好地满足对信号进行实时处理的要求。 　　关键词：分段卷积；实时处理；重叠相加法;FPGA 　　中图分类号：TN911.72 文献标识码：B 文章编号：1004373X(2008)1714602 　　 　　Design and Implementation of a Fast Algorithm of Partition Convolution 　　CHEN Yue 　　(College of Jincheng,Nanjing University of Aeronautics and Astronautics,Nanjing,211156,China) 　　 　　Abstract：The linear convolution usually needs two finite-time input sequences,but in the actual application,one of the input sequences is very long,so linear convolution can′t process the data in a real-time.On this situation,the partition convolution is an effective solution.This paper proposes a fast algorithm of partition convolution module,using pipeline structure in FPGA to real-time operate the data.The practical results indicate that this method is correct and effective. 　　 　　Keywords：partition convolution;real-time processing;overlap-add method;FPGA 　　 　　1 引 言 　　 　　在理论上，信号??x(n)(长度为M)通过数字系统h(n)(长度为N)后得到输出y(n)(长度为M+N－1)，其中y(n)=x(n)*h(n)。在实际工程中，这两个输入序列的长度可能相差很大，如N??M。在这种情况下，常采用快速线性卷积的方法来完成数据的处理，该方法是先将长序列x(n)存储完毕之后再和短序列h(n)做卷积运算［1,2］。但是，这样做存在以下两个问题： 　　(1) 需要对短序列h(n)补充零点使其长度和长序列x(n)一致，这样必然带来计算机存储量和存储空间的增加； 　　(2) 需要等待长序列x(n)全部产生之后才能将其输入系统并处理，这样必然使得整个系统存在较大的延时，不能实现对信号进行“实时”处理的要求。?? 　　为了解决这个问题，本文设计了一种工程应用型的准实时处理快速线性卷积的方法，即采用分段卷积［3,4］的设计思想将长序列??x(n)分解成若干个较小的段，如每段长为L，得x??k(n)，k=1,2,…M/L。然后对每一个输入段x??k(n)分别计算与h(n)做卷积得到的相应的输出y??k(n)，并将y??k(n)按一定的规则首尾相加，即可得到完整的输出y(n)。由于数字系统的单位抽样响应h(n)?б话惚冉隙?(如FIR数字滤波器)，这样做即可以实现数字信号的“实时”处理要求。 　　根据分段方式的不同，分段卷积法分为重叠相加法(Overlap-add Method)和重叠保留法(Overlap-save Method)两种，本文主要以前者为基础完成设计。 　　 　　2 重叠相加法的基本原理 　　 　　假设??x(n)的长度为M，h(n)的长度为N，且N??M，将长序列x(n)按如下方式分成长度为L的连续有限长序列：?? 　　x??k(n)=x(n), kL≤n≤(k+1)L－1 　　0,其他 ?ИЪ醋佣沃?间不存在重叠。 　　因此，根据卷积的分配律得：?И?y(n)=x(n)*h(n)=∑∞k=0x??k(n)*h(n) 　　=∑∞k=0［x??k(n)*h(n)］=∑∞k=0y??k(n) 　　因此，将原始序列按不重叠分段后，??x(n)和h(n)线性卷积等于各子段x??k(n)和h(n)线性卷积y??k(n)之和。由于x??k