蝶形、同址和變址計(jì)算方法及公式

蝶形運(yùn)算的特點(diǎn)是，首先每一個(gè)蝶形運(yùn)算都需要兩個(gè)輸入數(shù)據(jù)，計(jì)算結(jié)果也是兩個(gè)數(shù)據(jù)，與其它結(jié)點(diǎn)的數(shù)據(jù)無關(guān)，其它蝶形運(yùn)算也與這兩結(jié)點(diǎn)的數(shù)據(jù)無關(guān)、因此，一個(gè)蝶形 運(yùn)算一旦計(jì)算完畢，原輸入數(shù)據(jù)便失效了。這就意味著輸出數(shù)據(jù)可以立即使用原輸 人數(shù)據(jù)結(jié)點(diǎn)所占用的內(nèi)存。原來的數(shù)據(jù)也就消失了。輸出、輸人數(shù)據(jù)利用同一內(nèi)存單 元的這種蝶形計(jì)算稱為同位(址)計(jì)算。
可以看出， 每一級的蝶形的輸入與輸出在運(yùn)算前后可以存儲(chǔ)在同一地址(原來位置上)的存儲(chǔ)單元中，這種同址運(yùn)算的優(yōu)點(diǎn)是可以節(jié)省存儲(chǔ)單元，從而降低對計(jì)算機(jī)存儲(chǔ)量的要求或降低硬件實(shí)現(xiàn)的成本。

3．變址計(jì)算
從圖3. 19所示的流程圖看出，同址計(jì)算要求輸入x(n)是“混序”排列的。所謂輸入為“混 序”，并不是說輸入是雜亂無章的，實(shí)際上它是有規(guī)律的。如果輸入x(n)的序號用二進(jìn)制碼來 表示，就可以發(fā)現(xiàn)輸入的順序恰好是正序輸入的“碼位倒置”，表3. 3列出了這種規(guī)律。

在實(shí)際運(yùn)算中，按碼位倒置順序輸入數(shù)據(jù)x(n)，特別當(dāng)N較大時(shí)，是很不方便的。因此，數(shù)據(jù)總是按自然順序輸入存儲(chǔ)，然后通過“變址”運(yùn)算將自然順序轉(zhuǎn)換成碼位倒置順序存儲(chǔ)。實(shí)現(xiàn)這種轉(zhuǎn)換的程序可用圖3. 21來說明。

圖中用n表示自然順序的標(biāo)號，用l表示碼位倒置的標(biāo)號。當(dāng)l=n時(shí)，x(n)和x(l)不必互相調(diào)換。當(dāng)l≠n時(shí)， 必須將x(l)和x(n)互相調(diào)換，但只能調(diào)換一次，為此必須規(guī)定每當(dāng)l>n時(shí)，要將x(l)和x(n)相互調(diào)換，即把原來存放x(n)的存儲(chǔ)單元中的數(shù)據(jù)調(diào)入存儲(chǔ)x(l)的存儲(chǔ)單元中，而把原來存儲(chǔ)x(l)的存儲(chǔ)單元中的數(shù)據(jù)調(diào)入到存儲(chǔ)x(n)的存儲(chǔ)單元中。
這樣，按自然序輸入的數(shù)據(jù)x(n)經(jīng)過變址計(jì)算后變成了碼位倒置的排列順序，便可進(jìn)入第一級的蝶形運(yùn)算。
最后介紹一下時(shí)間抽選FFT算法的另外一些形式的流程圖。對于任何流程圖，只要保持 各節(jié)點(diǎn)所連支路及其傳輸系數(shù)不變，則不論節(jié)點(diǎn)位置怎樣排列，所得到的流程圖總是等效的，因而都能得到DFT的正確結(jié)果，只是數(shù)據(jù)的提取和存儲(chǔ)次序不同而已。
把圖3. 19中與x(4)水平相鄰的所有節(jié)點(diǎn)和與x(1)水平相鄰的所有節(jié)點(diǎn)交換，把與x(6)水平相鄰的所有節(jié)點(diǎn)和與 x(3)水平相鄰的所有節(jié)點(diǎn)交換，而與x(2)、x(5)和x(7)水平相鄰各節(jié)點(diǎn)位置不變，就可以從圖3. 19得到圖3.22。圖3.22與圖3.19的區(qū)別只是節(jié)點(diǎn)的排列不同，而支路傳輸比，即WN的 各次冪保持不變。顯然圖3.22所示流程圖的輸入是正序(自然順序)排列的，輸出是碼位倒置 排列的，所以輸出要進(jìn)行變址計(jì)算。圖3. 22所示的流程圖相當(dāng)于最初由庫利和圖基給出的時(shí) 間抽選算法。

另一種形式的流程圖是將節(jié)點(diǎn)排列成輸入 和輸出兩者都是正序排列，但這類流程圖不能進(jìn)行同址計(jì)算，因而需要兩列 長度為N的復(fù)數(shù)存儲(chǔ)器。