LTE系統(tǒng)中FFT的研究與DSP實現(xiàn)
加入技術交流群
使用C語言編寫主函數，匯編語言編寫FFT算法的實現(xiàn)函數。程序中假設輸入數據最大長度為1 024，由于DSP C6455可以直接存取處理32 bit，所以在內存中定義了長度為8 192 bit作為存放輸出序列的內存空間。為了提高運算精確度，輸入數的實部和虛部分別占用一個字，在程序中進行復數相乘操作是采用匯編指令MPYHI。內存定義了長度為2 048 bit的Tempsequence作為存放倒序序列，并且建立了2張旋轉因子查找表，分別為Wr和Wi。
外循環(huán)中，在每次內循環(huán)之前從輸入比特序列中取出32 bit放入一個寄存器，作為一個內循環(huán)的輸入，內循環(huán)結束后，取下一個32 bit輸入比特更新這個寄存器。
內循環(huán)中，計算蝶形過程采用查表的方式。對于每一級，計算出需要的旋轉因子個數以及相同旋轉因子相距的間隔。計算蝶形過程時，首先提取出X(k)，根據相同旋轉因子間隔找到X(k+B)完成蝶形計算。考慮到旋轉因子的對稱性，在內存中存放旋轉因子時只存放一半，剩余的數據根據對稱性進行處理。圖2給出了FFT算法實現(xiàn)計算流程圖。本文引用地址：http://www.ex-cimer.com/article/149952.htm
按時間抽取法的FFT輸入序列是倒序，輸出序列是自然順序；按頻率抽取法的FFT輸入序列是自然順序，輸出序列是倒序的。不管采用哪種方法進行FFT計算，都需要倒序處理。倒序是整個FFT計算的重要部分，進行匯編程序時，按自然順序將輸入數據存入到存儲單元內，通過變址運算，將自然順序的序列按時間抽取法要求進行倒位。
重新排序之前，存儲單元Y中依次存放輸入數據，I表示當前輸入數據比特的順序數的十進制數值，I的取值從0到N-I；J表示當前倒序數的十進制數值。輸入序列的第一個和最后一個數的位置不需要倒序處理，完成倒序的外循環(huán)的次數為N-2。為了保證調換數據的正確性，需要檢測一下是否IJ，只有當IJ，才將Y(I)與Y(J)的內容互換。形成倒序數J以后，就可以實現(xiàn)變址功能，按照自然順序存放在存儲單元的數據重新按照倒序排列。圖3給出了實現(xiàn)倒序的匯編流程圖。
3 性能分析與總結
在DSP軟件實現(xiàn)中，通過指令并行，盡量優(yōu)化程序循環(huán)體，減少或消除程序中的’NOP’指令[6]。通過程序仿真運行，得到統(tǒng)計結果如表1所示。
從表中可以看出，當運用TMS320C64×DSP芯片實現(xiàn)時，由于處理器的超高主頻一般為1 GHz，一個指令周期耗時為1 ns，其運算速率非?？欤耆梢詽M足實時性信號處理。因此，采用旋轉因子查表法的實現(xiàn)方案不僅簡化了程序實現(xiàn)方法，還減少了模塊程序代碼編寫，節(jié)約了系統(tǒng)存儲空間。
本文提出了一種簡單有效的FFT算法實現(xiàn)方案，詳細介紹了算法在DSP的實現(xiàn)方法，并在TMS320C64x芯片上加以實現(xiàn)。程序運行結果表明，該算法能夠滿足TD-LTE系統(tǒng)的需求，具有可行性和高效性。該方案已應用于LTE-TDD無線綜合測試儀表的開發(fā)中。
參考文獻
[1] 丁玉美.數字信號處理[M].西安：西安電子科技大學出版社，2002.
[2] 何方白，張德民.數字信號處理[M].北京：高等教育出版社，2009.
[3] 3GPP TS 36.211 v9.0.0.Evolved universal terrestrial radio access(E-UTRA) physical channels and modulation (Release 9)[S].2009-12.
[4] SAIDI A.Decimation-in-time-frequency FFT algorithm[M]. Manuscript, To be published.1993.
[5] Texas Instruments Incorporated.TMS320C64x/C64x+DSP CPU and instruction set referenceguide[EB/OL].Http://www.ti.com.cn,2008.
[6] Texas Instruments Incorporated.TMS320C6000系列DSP編程工具與指南[M].田黎育，何佩琨，朱夢宇，譯.北京：清華大學出版社，2006.
評論