FFT(快速傅立葉變換)在無線通信、語音識別、圖像處理和頻譜分析等領域有著廣泛應用。在FFT運算中,核心操作是蝶形運算,而蝶形運算的主要操作是向量旋轉,實現向量旋轉可用複數乘法運算來實現,但複數乘耗費了FFT運算中大量的乘法器資源。CORDIC算法只需簡單的移位與加減運算就能實現向量旋轉,具有使用資源少、硬體規模小等優勢。因此在FFT蝶形運算中用其代替傳統FFT運算中的複數乘法器,可以獲得更好的性能。但傳統CORDIC算法中每次CORDIC迭代方向需由剩餘角度的計算來確定,影響了工作速度。為此,本文根據定點FFT復乘中旋轉因子的旋轉方向可預先確定的特點,對CORDIC算法做了一些改進,在節省資源的同時保證了工作速度。
1 CORDIC算法原理
假設直角坐標系中有一向量A(Xa,Ya),逆時針旋轉?茲角度後得到另一個向量B(Xb,Yb),這個過程可用如下矩陣表示:
針對這一特點,可在CORDIC算法上做一點改進,把旋轉因子所對應的CORDIC旋轉係數預先存在ROM中(人工計算旋轉係數比較麻煩,可用MATLAB編一段程序來計算,並把旋轉係數存為.mif文件以便ROM初始化),而不是把旋轉因子角度預先存在ROM中。這樣,在進行CORDIC運算時,直接從ROM中取出旋轉係數,從而減少計算Zi來確定下一步旋轉方向的步驟,減少CORDIC模塊設計的複雜性,提高了運算速度,並且旋轉係數不比旋轉因子角度佔用的ROM資源多。另外由於旋轉因子需要進行0°、-90°或+90°三種預旋轉,所以預旋轉還要分配兩位二進位數,這樣存儲旋轉係數的ROM就為18位的ROM。
改進的CORDIC算法結構如圖1所示,所有旋轉因子所對應的CORDIC旋轉係數都存儲在ROM中,通過地址產生器的控制實現序列與相應的旋轉因子的復乘運算。與傳統CORDIC算法相比去掉了預旋轉角與已旋轉角之差的計算來確定下一次旋轉方向的結構,不但增加了係數寄存器模塊,而且總體上結構更為簡單。此CORDIC算法還採用流水線結構提高了運算的速度,從當前VLSI的發展趨勢上來看,晶片內的門資源相對富裕,對流水線CORDIC的實現規模約束很小。此外,流水線CORDIC不存在迭代式CORDIC的反饋迴路,使得單元結構更加規則,有利於VLSI實現。