摘要：在對FFT（快速傅立葉變換）算法進(jìn)行研究的基礎(chǔ)上，描述了用FPGA實(shí)現(xiàn)FFT的方法，并對其中的整體結(jié)構(gòu)、蝶形單元及性能等進(jìn)行了分析。

關(guān)鍵詞：FPGA FFT

傅立葉變換是數(shù)字信號處理中的基本操作，廣泛應(yīng)用于表述及分析離散時(shí)域信號領(lǐng)域。但由于其運(yùn)算量與變換點(diǎn)數(shù)Ｎ的平方成正比關(guān)系，因此，在Ｎ較大時(shí)，直接應(yīng)用ＤＦＴ算法進(jìn)行譜變換是不切合實(shí)際的。然而，快速傅立葉變換技術(shù)的出現(xiàn)使情況發(fā)生了根本性的變化。本文主要描述了采用ＦＰＧＡ來實(shí)現(xiàn)２ｋ／４ｋ／８ｋ點(diǎn)ＦＦＴ的設(shè)計(jì)方法。
１　整體結(jié)構(gòu)

一般情況下，Ｎ點(diǎn)的傅立葉變換對為：

其中，ＷＮ＝ｅｘｐ(－２ ｐｉ／Ｎ)。Ｘ(ｋ)和ｘ(ｎ)都為復(fù)數(shù)。與之相對的快速傅立葉變換有很多種,如ＤＩＴ(時(shí)域抽取法)、ＤＩＦ（頻域抽取法）、Ｃｏｏｌｅｙ－Ｔｕｋｅｙ和Ｗｉｎｏｇｒａｄ等。對于２ｎ傅立葉變換，Ｃｏｏｌｅｙ－Ｔｕｋｅｙ算法可導(dǎo)出ＤＩＴ和ＤＩＦ算法。本文運(yùn)用的基本思想是Ｃｏｏｌｅｙ－Ｔｕｋｅｙ算法，即將高點(diǎn)數(shù)的傅立葉變換通過多重低點(diǎn)數(shù)傅立葉變換來實(shí)現(xiàn)。雖然ＤＩＴ與ＤＩＦ有差別，但由于它們在本質(zhì)上都是一種基于標(biāo)號分解的算法，故在運(yùn)算量和算法復(fù)雜性等方面完全一樣，而沒有性能上的優(yōu)劣之分，所以可以根據(jù)需要任取其中一種，本文主要以ＤＩＴ方法為對象來討論。

Ｎ＝８１９２點(diǎn)ＤＦＴ的運(yùn)算表達(dá)式為：

式中，ｍ＝(４ｎ１＋ｎ２)(２０４８ｋ１＋ｋ２)(ｎ＝４ｎ１＋ｎ２，ｋ＝２０４８ｋ１＋ｋ２)其中ｎ１和ｋ２可?。?１,．．．,２０４７,ｋ１和ｎ２可取０,１,２,３。

由式（３）可知，８ｋ傅立葉變換可由４×２ｋ的傅立葉變換構(gòu)成。同理，４ｋ傅立葉變換可由２×２ｋ的傅立葉變換構(gòu)成。而２ｋ傅立葉變換可由１２８×１６的傅立葉變換構(gòu)成。１２８的傅立葉變換可進(jìn)一步由１６×８的傅立葉變換構(gòu)成，歸根結(jié)底，整個傅立葉變換可由基２、基４的傅立葉變換構(gòu)成。２ｋ的ＦＦＴ可以通過５個基４和１個基２變換來實(shí)現(xiàn)；４ｋ的ＦＦＴ變換可通過６個基４變換來實(shí)現(xiàn)；８ｋ的ＦＦＴ可以通過６個基４和１個基２變換來實(shí)現(xiàn)。也就是說：ＦＦＴ的基本結(jié)構(gòu)可由基２／４模塊、復(fù)數(shù)乘法器、存儲單元和存儲器控制模塊構(gòu)成，其整體結(jié)構(gòu)如所示。

中，ＲＡＭ用來存儲輸入數(shù)據(jù)、運(yùn)算過程中的中間結(jié)果以及運(yùn)算完成后的數(shù)據(jù)，ＲＯＭ用來存儲旋轉(zhuǎn)因子表。蝶形運(yùn)算單元即為基２／４模塊，控制模塊可用于產(chǎn)生控制時(shí)序及地址信號，以控制中間運(yùn)算過程及輸出結(jié)果。
２　蝶形運(yùn)算器的實(shí)現(xiàn)

基４和基２的信號流如所示。圖中，若Ａ＝ｒ０＋ｊ＊ｉ０，Ｂ＝ｒ１＋ｊ＊ｉ１，Ｃ＝ｒ２＋ｊ＊ｉ２，Ｄ＝ｒ３＋ｊ＊ｉ３是要進(jìn)行變換的信號，Ｗｋ０＝ｃ０＋ｊ＊ｓ０＝１，Ｗｋ１＝ｃ１＋ｊ＊ｓ１，Ｗｋ２＝ｃ２＋ｊ＊ｓ２，Ｗｋ３＝ｃ３＋ｊ＊ｓ３為旋轉(zhuǎn)因子，將其分別代入中的基４蝶形運(yùn)算單元，則有：

Ａ′＝[ｒ０＋(ｒ１×ｃ１－ｉ１×ｓ１)＋(ｒ２×ｃ２－ｉ２×ｓ２)＋(ｒ３×ｃ３－ｉ３×ｓ３)]＋ｊ[ｉ０＋(ｉ１×ｃ１＋ｒ１×ｓ１)＋(ｉ２×ｃ２＋ｒ２×ｓ２)＋(ｉ３×ｃ３＋ｒ３×ｓ３)] ?。ǎ矗?/P>

Ｂ′＝[ｒ０＋(ｉ１×ｃ１＋ｒ１×ｓ１)－(ｒ２×ｃ２－ｉ２×ｓ２)－(ｉ３×ｃ３＋ｒ３×ｓ３)]＋ｊ[ｉ０－(ｒ１×ｃ１－ｉ１×ｓ１)－(ｉ２×ｃ２＋ｒ２×ｓ２)＋(ｒ３×ｃ３－ｉ３×ｓ３)] 　(５）

Ｃ′＝[ｒ０－(ｒ１×ｃ１－ｉ１×ｓ１)＋(ｒ２×ｃ２－ｉ２×ｓ２)－(ｒ３×ｃ３－ｉ３×ｓ３)]＋ｊ[ｉ０－(ｉ１×ｃ１＋ｒ１×ｓ１)＋(ｉ２×ｃ２＋ｒ２×ｓ２)－(ｉ３×ｃ３＋ｒ３×ｓ３)] （６）

Ｄ′＝[ｒ０－(ｉ１×ｃ１＋ｒ１×ｓ１)－(ｒ２×ｃ２－ｉ２×ｓ２)＋(ｉ３×ｃ３＋ｒ３×ｓ３)]＋ｊ[ｉ０＋(ｒ１×ｃ１－ｉ１×ｓ１)－(ｉ２×ｃ２＋ｒ２×ｓ２)－(ｒ３×ｃ３－ｉ３×ｓ３)] （７）

而在基２蝶形中，Ｗｋ０和Ｗｋ２的值均為１，這樣，將Ａ，Ｂ，Ｃ和Ｄ的表達(dá)式代入中的基２運(yùn)算的四個等式中，則有：

Ａ′＝ｒ０＋(ｒ１×ｃ１－ｉ１×ｓ１)＋ｊ[ｉ０＋(ｉ１×ｃ１＋ｒ１×ｓ１)] （８）

Ｂ′＝ｒ０－ (ｒ１×ｃ１－ｉ１×ｓ１)＋ｊ[ｉ０－(ｉ１×ｃ１＋ｒ１×ｓ１)] ?。ǎ梗?/P>

Ｃ′＝ｒ２＋(ｒ３×ｃ３－ｉ３×ｓ３)＋ｊ[ｉ０＋(ｉ３×ｃ３＋ｒ３×ｓ３)] （１０）

Ｄ′＝ｒ２－(ｒ３×ｃ３－ｉ３×ｓ３)＋ｊ[ｉ０－(ｉ３×ｃ３＋ｒ３×ｓ３)] （１１）

在上述式（４）～（１１）中有很多類同項(xiàng)，如ｉ１×ｃ１＋ｒ１×ｓ１和ｒ１×ｃ１－ｉ１×ｓ１等，它們僅僅是加減號的不同，其結(jié)構(gòu)和運(yùn)算均類似，這就為簡化電路提供了可能。同時(shí)，在蝶形運(yùn)算中，復(fù)數(shù)乘法可以由實(shí)數(shù)乘法以一定的格式來表示，這也為設(shè)計(jì)復(fù)數(shù)乘法器提供了一種實(shí)現(xiàn)的途徑。

以基４為例，在其運(yùn)算單元中，實(shí)際上只需做三個復(fù)數(shù)乘法運(yùn)算，即只須計(jì)算ＢＷｋ１、ＣＷｋ２和ＤＷｋ３的值即可，這樣在一個基４蝶形單元里面，多只需要３個復(fù)數(shù)乘法器就可以了。在實(shí)際過程中，在不提高時(shí)鐘頻率下，只要將時(shí)序控制好便可利用流水線（Ｐｉｐｅｌｉｎｅ）技術(shù)并只用一個復(fù)數(shù)乘法器就可完成這三個復(fù)數(shù)乘法，大大節(jié)省了硬件資源。
基2和基4蝶形算法的信號流圖
３?。疲疲缘牡刂?/B>

ＦＦＴ變換后輸出的結(jié)果通常為一特定的倒序,因此，幾級變換后對地址的控制必須準(zhǔn)確無誤。

倒序的規(guī)律是和分解的方式密切相關(guān)的，以基８為例，其基本倒序規(guī)則如下：

基８可以用２×２×２三級基２變換來表示，則其輸入順序則可用二進(jìn)制序列（ｎ１ ｎ２ ｎ３）來表示，變換結(jié)束后，其順序?qū)⒆優(yōu)椋ǎ睿?ｎ２ ｎ１），如：Ｘ０１１→ ｘ１１０，即輸入順序?yàn)椋?，輸出時(shí)順序變?yōu)椋丁?/P>

更進(jìn)一步，對于基１６的變換，可由２×２×２×２，４×４，４×２×２等形式來構(gòu)成，相對于不同的分解形式，往往會有不同的倒序方式。以４×４為例，其輸入順序可以用二進(jìn)制序列（ｎ１ ｎ２ ｎ３ ｎ４）來表示變換結(jié)束后，其順序可變?yōu)椋ǎǎ睿?ｎ４）（ｎ１ ｎ２）），如： Ｘ０１１１→ ｘ１１０１。即輸入順序?yàn)椋罚敵鰰r(shí)順序變?yōu)椋保场?/P>

在２ｋ／４ｋ／８ｋ的傅立葉變換中，由于要經(jīng)過多次的基４和基２運(yùn)算，因此，從每次運(yùn)算完成后到進(jìn)入下運(yùn)算前，應(yīng)對運(yùn)算的結(jié)果進(jìn)行倒序，以保證運(yùn)算的正確性。
４　旋轉(zhuǎn)因子

Ｎ點(diǎn)傅立葉變換的旋轉(zhuǎn)因子有著明顯的周期性和對稱性。其周期性表現(xiàn)為：

FFT之所以可使運(yùn)算效率得到提高，就是利用