引言

在需要硬件實現對數運算的場合，其精度和速度是必須考慮的問題。目前硬件實現對數變換的方法主要有查表法、泰勒公式展開法和線性近似法。查表法所需要的存儲單元隨著精度的增加或輸入值范圍的增大而成指數增加；泰勒公式展開法需要乘法器，面積大不易實現；線性近似法的精度有限，且需要誤差校正電路，實現較難。

本文利用CORD IC算法在FPGA上實現了高速自然對數變換器。CORD IC算法即坐標旋轉數字計算方法最初由J.D.Volder于1959年提出，其基本思想是用一系列與運算基數相關的角度的不斷偏擺從而逼近所需旋轉的角度。1971年J.S.Walter提出統一的CORD IC算法，把圓周旋轉、直線旋轉和雙曲旋轉統一到同一個CORD IC迭代方程里，為同一硬件實現多功能運算提供了前提。由于它將許多復雜的算術運算化成簡單的加法和移位操作，在不影響運算速度和精度的情況下，極大的降低了硬件設計的復雜性，節約了硬件資源。利用CORD IC算法可以直接實現乘法、除法、正余弦函數，反正切函數、雙曲函數等，對輸入進行適當的初始化可以實現正切、雙曲正切、對數和指數等函數。

在數字信號處理領域用的較多的是DSP和FPGA。DSP的優勢源于多數信號處理算法的乘-累加運算（MAC）都是非常密集的。FPGA通過多極流水線架構也能夠用來實現MAC單元，并且FPGA技術可以通過一個芯片上的多級MAC單元來提供更多的帶寬，速度可以比數字信號處理芯片快，并且功耗較低。CORD IC算法完全由移位和加法操作完成，因此利用FPGA可以實現更高的運算速度。本文采用流水線結構在FPGA上實現基于CORD IC的對數變換，可以達到80MHz的處理速度。

CORD IC算法實現自然對數運算

CORD IC算法最初是用于計算三角函數的，后來由于其算法的簡單、硬件易于實現等多種優勢，而被廣泛的用于多種初等函數的運算中（包括三角函數、乘除法運算、指數運算、對數運算等）。本文主要利用CORD IC算法的雙曲旋轉法實現自然對數運算。

在雙曲坐標系下，CORD IC算法的迭代方程為：

由于 ，所以迭代序列必須從n=1開始，為保證迭代序列收斂，因此迭代序列n的取值從第4項開始每隔3n+1項必須重復一次，即n=1，2，3，4，4，5，…，40，40，…。

在向量模式下，經n次迭代后的輸出方程為：

因為

所以令x=t+1，y=t-1

則

所以對于t，如果我們要求lnt，只要做如下初始化：

X=t+1，y=t-1，z=0

則輸出z=0.5ln（t），只需要在CORD IC之后做一次左移即可。

如（2）式所示，為保證迭代序列的收斂，|tanh-1（y0/x0）|≤1.1182，因此|y/x|max≈0.8069，n→∞，反雙曲正切的定義域為（-1，1），可見函數的輸入范圍受到了限制。解決的方法是增加n為負數的迭代，改進的算法公式為：

當n≤0時

當n》0時

收斂的范圍變成

|tanh-1（y/x）|≤θmax

其中

當M=5時，θmax=12.4264，函數tanh-1的范圍是［-12.4264，12.4264］。也就是說此時y/x可以接近于［-1，1］，幾乎覆蓋tanh-1的整個定義域。因此硬件實現過程中可以從-5開始迭代。

對數運算的FPGA實現

CORD IC算法完全由移位和相加完成，很容易在硬件上實現。由于FPGA具有并行處理能力，利用FPGA實現對數變換，速度可以比數字信號處理芯片快，以滿足某些高速處理的要求。本文采用的FPGA芯片是Altera公司的cyclone系列芯片EP1C6Q240C8。該芯片內部共有邏輯單元5980個，支持近12萬門的設計，內部嵌有約12Kbyte的RAM，包含2個生成時鐘的鎖相環，最大用戶I/O數185個，滿足設計要求。

實現方法

· 預處理單元

欲利用CORD IC方法求自然對數，必須對輸入進行初始化，經迭代運算后才能得到對數運算值。該對數變換器的輸入為16位數，在預處理單元中將輸入分別加減一，并將位寬擴大為40位，最高位作為符號位覆值給第一次迭代的x0和y0，如圖1所示，圖中s代表符號位。擴大位寬可以提高輸出精度。

圖1初始迭代值x0，y0

· CORD IC單元

CORD IC單元是實現對數變換器的核心。本文利用流水線結構實現CORD IC算法，其結構如圖2所示。在設計中，采用由28級CORD IC運算單元組成的流水線結構，為擴大輸入范圍，從n=-5開始迭代，移位序列為［7，6，5，4，3，2，1，2…28］。前6級根據（3）式進行迭代，后22級根據（4）式進行迭代。經過28級流水線運算后，y變為0，z左移一位就是要求的對數值。每一級電路結構主要包括2個移位寄存器和3個加減法器，這些移位寄存器各自有不同的固定的移位次數，加減法選擇由該級中y的最高位（符號位）決定。θn為第n次迭代的旋轉角度，并作為常數直接連到了累加器上，不需要存儲空間和讀取時間。

圖2 CORDIC流水線結構

· 后處理單元

由CORD IC得到的z=1/2ln（t），因此將結果左移一位，并截取高16位作為最終的輸出。其中最高位為符號位，最大輸入值65535的對數值為11.0903，對應的輸出為7FFF，其余輸出均除以對應的值即得到相應的對數值。

實驗結果

在Quartus II 5.1軟件環境下使用Verilog HDL語言完成了上述各算法，并在cyclone系列芯片EP1C6Q240C8上實現。圖3 為對數運算時序仿真圖。表1為對數運算結果與理論值的比較。

圖3對數運算時序仿真圖

由表1可看出，該對數運算器的輸出誤差為10-4數量級。由于采用流水線結構，能夠在執行進程的同時輸入數據，從而極大的提高了程序的運行效率。該設計需要30個時鐘周期獲取第一個計算結果，而只需要一個時鐘周期來獲取隨后的計算結果。利用Quartus Ⅱ5.1軟件進行時序分析，該運算器的最高頻率可達到80MHz。該運算器適用于高速大數據量的數據處理。

表1對數運算結果與理論值的比較

結語

利用對數變換可以將乘除法變換為加減法實現，有利于乘除法在硬件中的實現。由于CORD IC算法完全由移位和相加運算完成，降低了復雜性，易于硬件的實現。利用CORD IC算法在FPGA上設計了一種自然對數變換器。實現過程中采用流水線結構，提高了系統的運行效率。實驗結果表明該對數運算器的輸出誤差為10-4數量級，最高頻率可達到80MHz。該運算器適用于高速大數據量的數據處理。

責任編輯：gt