基于FPGA的FFT算法優(yōu)化及其在磁共振譜儀中的應(yīng)用

2 FFT算法實(shí)現(xiàn)

FFT算法在FPGA上實(shí)現(xiàn)的過程中，信號的值、系統(tǒng)的系數(shù)和運(yùn)算中的結(jié)果都存儲在有限字長的存儲單元中，從而導(dǎo)致了設(shè)計(jì)時(shí)的無限精度轉(zhuǎn)變成實(shí)現(xiàn)時(shí)的有限精度，必將產(chǎn)生相對于原設(shè)計(jì)系統(tǒng)的誤差，嚴(yán)重時(shí)會將由于雙通道不平衡產(chǎn)生的鏡像峰湮滅，從而使整個設(shè)計(jì)失去意義。

為了實(shí)現(xiàn)FFT實(shí)時(shí)運(yùn)算，基于FPGA的FFT信號處理模塊是關(guān)鍵，并且要求此模塊能在頻率至少為210 MHz的系統(tǒng)時(shí)鐘下穩(wěn)定工作。同時(shí)又因?yàn)椴⑿械腇FT設(shè)計(jì)需要占用大量的資源，資源使用率過大會制約布局布線后的時(shí)序收斂。為了平衡資源與速度間的矛盾，整個設(shè)計(jì)通過64點(diǎn)FFT并行模塊的復(fù)用來實(shí)現(xiàn)。

2．1 有限字長效應(yīng)及其優(yōu)化措施

在FFT算法中，采用蝶形計(jì)算，如圖3所示。

對于基2時(shí)間抽選FFT算法，蝶形公式如下：

式中：N為FFT點(diǎn)數(shù)；Nm和P為兩個同迭代次數(shù)m有關(guān)的量；Xm-1為ADC輸出信號經(jīng)過m-2級蝶形運(yùn)算得到的計(jì)算結(jié)果；旋轉(zhuǎn)因子由式(7)給出，每個蝶形包含兩個復(fù)數(shù)乘法。

由于存儲單元有限，必須對計(jì)算結(jié)果截取，進(jìn)行定點(diǎn)化，但是同時(shí)又會引入舍入誤差en，這種現(xiàn)象即為有限字長效應(yīng)。考慮到每個復(fù)數(shù)乘法相當(dāng)于4個實(shí)數(shù)乘法，因此有限字長時(shí)復(fù)數(shù)乘法的實(shí)際乘積可以表示為：

式中e1，e2，e3，e4分別為其對應(yīng)的實(shí)數(shù)乘法的舍入誤差。因此可進(jìn)一步建立如圖4所示的蝶形計(jì)算舍入誤差模型。