基于多DSPs的FFT并行算法.pdfVIP

V01．3l 第31卷专辑 南昌大学学报(理科版) Suppl 2007年7月 Jo哪al“Nanch鲫gunive璐ity(NaturalScience) Jul．2007 文章编号：1006—04“(2007)专辑一0189一03 基于多DSPs的FFT并行算法 笪丽琴，周朝霞 (南昌大学自动化系，江西南昌330031) 处理和其它很多领域中有着越来越重要的应用。当需要作大点数FF-II时，先将数据在几片DSP8上并行地进行前 几级蝶形运算，然后将结果通过双口存储器(DRAM)传给另一片DsPs进行后几级蝶形运算，以实现大点数的Ff’1． 运算，具有较高的加速比和总效率，适合实时处理系统。 关键词：快速傅立叶变换；多DsP8；并行算法 中图分类号：哪01．6 文献标识码：A 傅立叶变换是数字信号处理中对信号进行分析 时经常采用的一种方法，是一种将时域信号转变为 频域信号的变换形式。在频域分析中，频谱分析是 信号分析的重要内容。它反映了系统性能的好坏。 而(七) ?噼 Cool眵和Tukey(1965)所研究出的计算离散傅氏变图1蝶形运算单元信号流图 换的快速傅立叶变换，将计算量的复杂度从N2级降 由此可以看出，快速傅立叶变换(F丌)是离散 到Nlo＆N级，从而使得FFT在数字信号处理、石油 傅立叶变换(D胛)的一种快速算法，它是利用蝶形 勘探、地震预报、医学断层诊断、编码理论、量子物理 因子WN的周期性与对称性，将N点的DFT分解为 及概率论等领域都得到了广泛的应用。 两个N／2点的D盯，这样DFl总计算量变为原来的 1 FI叮算法 一半，并且可以继续分解下去，这样将长序列的D丌 分解成若干个短序列的D盯计算，来提高D兀的效 按时问抽取的基一2FFT(DIrr)算法是肿最常 率，然后再按一定规则将其合并，从而得到整个的 用的一种表现形式。 DFTr。最小的分解点数称为基数。对于N=2M点的 连续信号x(t)的傅立叶变换其离散形式为： DFT都可以用基2的F丌来实现，这样其计算量可 x(k∞。)：∑1x(nT)e-jk州：∑1x(nT)e—j和 x(k‰)=蛋x(nT)e1。哪。=蛋x(nT)e叫种 n2U noU 以减为(N／2)log：N个加乖运算，在一定程度上减少 (1) 了D兀’的计算量。 令w。=e一臀，并省略符号∞。和T，则上式可写 2 多DSPs实现的F兀并行算法 成： N一1 以N=16点肿为例，如图2所示，可将输人 k=o，l，2，…，N—I(2) x(k)=暑x(n)w守 noU 数据看成一个矩阵A，A表示刚开始数据按循环分 由蝶形因子w。的周期性：w紫=w警N+8’= 配给P(以2台为例)台处理器，每一行表示一台