定點DSP中高精度除法的實現方法

2 高精度算法

由于經典算法在結果上存在精度不高的問題，所以應考慮精度盡可能高的除法算法。

和經典算法一致，本算法也采用移位相減實現除法。但不同的是，在進行移位相減之前算法對除數與被除數進行了一些處理，并且在移位相減時不再分兩種情況。下面具體介紹算法的流程。

首先給出被除數y和除數x。由于定點DSP中所有的數據都以整數的形式存在，所以在這里就把y和x看成是整數，不考慮它們本身的Q值。在圖1所示的流程圖中，A和B是C5416的2個累加器。算法第一步是將x和y相乘，結果存放在B中，以便在最后從B中提取商的符號(商的符號與B的符號相同)。第二步分別對x和y取絕對值以進行無符號數相除。第三步，分別判斷x和y在最高有效位之前共有多少個O，分別記為n和m，即在x和y的分別有(16一n)和(16一m)個有效位。然后將x和y分別左移(n一1)位和(m一2)位，即x=x·2n-1，y=y·2m-2。第四步，將前面所得x和y再進行移位相減，把y存放在A的高位，即AH=y，利用SUBC指令進行1 5次移位相減。在所得的結果當中，A的低位即為商的絕對值，B的符號即為商的符號，兩者結合得到quot值。由于在移位相減之前分別對除數和被除數有一個左移放大的操作，所以在結果中就必須考慮商的Q值。

假設給定的除數和被除數是x0和y0，則進行移位相減的兩數分別為x=x0·2n-1，y=y0·2m-2。根據計算的過程可以得到商的Q值為(14+m-n)，所以最后的結果應該包括quot值和Q值(14+m—n)兩部分。即商=quot·2-Q。

3 實驗結果及討論

前一節中詳細介紹了高精度算法的具體流程，可以發現其與經典算法相比，不同之處在于進行移位相減之前把除數和被除數進行了不同位數的左移，保證進行移位相減的除數的數據位最高位是l，即14位為1，而被除數的第13位為1。正是這個左移使結果的精度得到了提高。

在所引用的經典算法中，當1被除數l

利用高精度算法在TMS320C5416中計算上面的4組數據，所得結果如表2所列。

上面的實驗說明，該算法使計算結果的精度有了很大的提高。取表中第三組數據進行詳細分析，y=479、x=240，利用經典算法結果為整數1。利用改進后的算法，得到的結果有兩部分：quot值=3FDDH，Q值=13。由這兩部分可以算得精確的商值1.995 727 539 062 5，而479/240=1.995 83，計算結果與實際值之間的誤差為O.000105 794 270 83。當Q值=13時，定點數據表示形式允許的誤差為2-13=O.000 122 070 312 5，可見計算結果誤差在允許的誤差范圍之內，也就是說計算結果滿足所推導出的精度的范圍。對比前后兩種算法，可以發現這里的計算結果精度由原來的1提高到了2-13。

通過分析可知，該算法使計算結果在精度上有了非常大的提高，但是在效率上則有所降低。如表3所列：

表3中的時鐘周期指的是DSP的機器時鐘周期。由表3可以發現高精度算法對于不同的數據所需的周期數是不一樣的，而且算法中計算被除數和除數的高位無效位數占用了大部分的周期數。

經典算法在效率上比高精度算法高，主要是因為在高精度算法中高位無效位數的計算效率過低。用C5416計算|被除數|和|除數|的高位無效位數，效率低是由于在C5416中沒有現成的指令來探測數據中哪一位為最高有效位。表3給出的時鐘周期數是在筆者的程序中所得到的結果，由于時間匆忙和優化方面能力的局限性，所以并不能保證在程序優化方面做得很好，所以在C5416中對該部分程序再進行更好的優化以減少執行的時間是很可能的。

在TI公司C6000系列的定點DSP中，有專門的指令LMBD來快速得到累加器中從符號位開始無效0位的個數，只需幾個周期就能完成C5416里高位無效位數的計算，所以在C6000中完成整個算法的效率與經典算法相比肯定能夠提高。并且在C6000系列中對不同的數據進行計算所需的時鐘周期數是一樣的，因此在C6O0O系列定點DSP中該算法有著更好的實用性。

結 語

本文介紹了高精度除法在定點DSP中的實現，通過與經典除法算法的對比，顯示了其高精度的優點。定點TMS320C5416實驗表明該算法使商的精度得到了很大的提高，但是在計算的效率上有待進一步提高。