一DSP定點(diǎn)算數(shù)運(yùn)算1數(shù)的定標(biāo)在定點(diǎn)DSP芯片中,采用

1、一 定點(diǎn)算數(shù)運(yùn)算 1 數(shù)的定標(biāo) 在定點(diǎn)DSP芯片中，采用定點(diǎn)數(shù)進(jìn)行數(shù)值運(yùn)算，其操作數(shù)一般采用整型數(shù)來(lái)表示。一個(gè)整型數(shù)的最大表示范圍取決于DSP芯片所給定的字長(zhǎng)，一般為16位或24位。顯然，字長(zhǎng)越長(zhǎng)，所能表示的數(shù)的范圍越大，精度也越高。如無(wú)特別說(shuō)明，本書均以16位字長(zhǎng)為例。 DSP芯片的數(shù)以2的補(bǔ)碼形式表示。每個(gè)16位數(shù)用一個(gè)符號(hào)位來(lái)表示數(shù)的正負(fù)，0表示數(shù)值為正，l則表示數(shù)值為負(fù)。其余15位表示數(shù)值的大小。因此， 二進(jìn)制數(shù)0010000000000011b=8195 二進(jìn)制數(shù)1111111111111100b= -4 對(duì)DSP芯片而言，參與數(shù)值運(yùn)算的數(shù)就是16位的整型數(shù)。但在許多情況下，數(shù)學(xué)運(yùn)算

2、過(guò)程中的數(shù)不一定都是整數(shù)。那么，DSP芯片是如何處理小數(shù)的呢？應(yīng)該說(shuō)，DSP芯片本身無(wú)能為力。那么是不是說(shuō)DSP芯片就不能處理各種小數(shù)呢？當(dāng)然不是。這其中的關(guān)鍵就是由程序員來(lái)確定一個(gè)數(shù)的小數(shù)點(diǎn)處于16位中的哪一位。這就是數(shù)的定標(biāo)。 通過(guò)設(shè)定小數(shù)點(diǎn)在16位數(shù)中的不同位置，就可以表示不同大小和不同精度的小數(shù)了。數(shù)的定標(biāo)有Q表示法和S表示法兩種。表1.1列出了一個(gè)16位數(shù)的16種Q表示、S表示及它們所能表示的十進(jìn)制數(shù)值范圍。 從表1.1可以看出，同樣一個(gè)16位數(shù)，若小數(shù)點(diǎn)設(shè)定的位置不同，它所表示的數(shù)也就不同。例如， 16進(jìn)制數(shù)2000H=8192，用Q0表示 16進(jìn)制數(shù)2000H=0.25，用Q15

3、表示 但對(duì)于DSP芯片來(lái)說(shuō)，處理方法是完全相同的。 從表1.1還可以看出，不同的Q所表示的數(shù)不僅范圍不同，而且精度也不相同。Q越大，數(shù)值范圍越小，但精度越高；相反，Q越小，數(shù)值范圍越大，但精度就越低。例如，Q0 的數(shù)值范圍是一32768到+32767，其精度為1，而Q15的數(shù)值范圍為-1到0.9999695，精度為1/32768=0.00003051。因此，對(duì)定點(diǎn)數(shù)而言，數(shù)值范圍與精度是一對(duì)矛盾，一個(gè)變量要想能夠表示比較大的數(shù)值范圍，必須以犧牲精度為代價(jià)；而想精度提高，則數(shù)的表示范圍就相應(yīng)地減小。在實(shí)際的定點(diǎn)算法中，為了達(dá)到最佳的性能，必須充分考慮到這一點(diǎn)。 浮點(diǎn)數(shù)與定點(diǎn)數(shù)的轉(zhuǎn)換關(guān)系可表示為：

4、 浮點(diǎn)數(shù)(x)轉(zhuǎn)換為定點(diǎn)數(shù)(xq)：xq=(int)x* 2Q 定點(diǎn)數(shù)(xq)轉(zhuǎn)換為浮點(diǎn)數(shù)(x)：x=(float)xq*2-Q 例如，浮點(diǎn)數(shù)x=0.5，定標(biāo)Q=15，則定點(diǎn)數(shù)xq=L0.5*32768J=16384，式中LJ表示下取整。反之，一個(gè)用Q=15表示的定點(diǎn)數(shù)16384，其浮點(diǎn)數(shù)為163幼*2-15=16384/32768=0.5。浮點(diǎn)數(shù)轉(zhuǎn)換為定點(diǎn)數(shù)時(shí)，為了降低截尾誤差，在取整前可以先加上0.5。 表1.1 Q表示、S表示及數(shù)值范圍 Q表示 S表示 十進(jìn)制數(shù)表示范圍 Q15 S0.15 -1x0.9999695 Q14 S1.14 -2x1.9999390 Q13 S2.13 -4x

5、3.9998779 Q12 S3.12 -8x7.9997559 Q11 S4.11 -16x15.9995117 Q10 S5.10 -32x31.9990234 Q9 S6.9 -64x63.9980469 Q8 S7.8 -128x127.9960938 Q7 S8.7 -256x255.9921875 Q6 S9.6 -512x511.9804375 Q5 S10.5 -1024x1023.96875 Q4 S11.4 -2048x2047.9375 Q3 S12.3 -4096x4095.875 Q2 S13.2 -8192x8191.75 Q1 S14.1 -16384x16383

6、.5 Q0 S15.0 -32768x32767 2 高級(jí)語(yǔ)言：從浮點(diǎn)到定點(diǎn) 我們?cè)诰帉慏SP模擬算法時(shí)，為了方便，一般都是采用高級(jí)語(yǔ)言(如C語(yǔ)言)來(lái)編寫模擬程序。程序中所用的變量一般既有整型數(shù)，又有浮點(diǎn)數(shù)。如例1.1程序中的變量i是整型數(shù)，而pi是浮點(diǎn)數(shù)，hamwindow則是浮點(diǎn)數(shù)組。 例1.1 256點(diǎn)漢明窗計(jì)算 int i；+ float pi=3.14l59； float hamwindow256； for(i=0；i(Qx-Qz)，若Qx=Qz z=temp(Qz-Qx)，若Qx32767，因此 Qx=1，Qy=0，Qz=0，則定點(diǎn)加法為： x=30000；y=20000； tem

7、p=200001=35000； 因?yàn)閦的Q值為0，所以定點(diǎn)值z(mì)=35000就是浮點(diǎn)值，這里z是一個(gè)長(zhǎng)整型數(shù)。當(dāng)加法或加法的結(jié)果超過(guò)16位表示范圍時(shí)，如果程序員事先能夠了解到這種情況，并且需要保持運(yùn)算精度時(shí)，則必須保持32位結(jié)果。如果程序中是按照16位數(shù)進(jìn)行運(yùn)算的，則超過(guò)16位實(shí)際上就是出現(xiàn)了溢出。如果不采取適當(dāng)?shù)拇胧?，則數(shù)據(jù)溢出會(huì)導(dǎo)致運(yùn)算精度的嚴(yán)重惡化。一般的定點(diǎn)DSP芯片都沒(méi)有溢出保護(hù)功能，當(dāng)溢出保護(hù)功能有效時(shí)，一旦出現(xiàn)溢出，則累加器ACC的結(jié)果為最大的飽和值(上溢為7FFFH，下溢為8001H)，從而達(dá)到防止溢出引起精度嚴(yán)重惡化的目的。 2.2乘法運(yùn)算的C語(yǔ)言定點(diǎn)模擬 設(shè)浮點(diǎn)乘法運(yùn)算的表達(dá)

8、式為： float x，y，z； z=xy； 假設(shè)經(jīng)過(guò)統(tǒng)計(jì)后x的定標(biāo)值為Qx，y的定標(biāo)值為Qy，乘積z的定標(biāo)值為Qz，則 z=xy zq*2-Qx=xq*yq*2-(Qx+Qy) zq=(xqyq)2Qz-(Qx+Qy) 所以定點(diǎn)表示的乘法為： int x，y，z； long temp； temp=(long)x； z=(temp*y)(Qx+Qy-Qz)； 例1.5定點(diǎn)乘法。 設(shè)x=18.4，y=36.8，則浮點(diǎn)運(yùn)算值為=18.4*36.8=677.12； 根據(jù)上節(jié)，得Qx=10，Qy=9，Qz=5，所以 x=18841；y=18841； temp=18841L； z=(18841L*188

9、41)(10+9-5)=354983281L14=21666； 因?yàn)閦的定標(biāo)值為5，故定點(diǎn)z=21666，即為浮點(diǎn)的z=21666/32=677.08。 2.3除法運(yùn)算的C語(yǔ)言定點(diǎn)摸擬 設(shè)浮點(diǎn)除法運(yùn)算的表達(dá)式為： float x，y，z； z=x/y； 假設(shè)經(jīng)過(guò)統(tǒng)計(jì)后被除數(shù)x的定標(biāo)值為Qx，除數(shù)y的定標(biāo)值為Qy，商z的定標(biāo)值為Qz，則 z=x/y zq*2-Qz=(xq*2-Qx)/(yq*2-Qy) zq=(xq*2(Qz-Qx+Qy)/yq 所以定點(diǎn)表示的除法為： int x，y，z； long temp； temp=(long)x； z=(temp(Qz-Qx+Qy)/y； 例1.6定點(diǎn)

10、除法。 設(shè)x=18.4，y=36.8，浮點(diǎn)運(yùn)算值為z=x/y=18.4/36.8=0.5； 根據(jù)上節(jié)，得Qx=10，Qy=9，Qz=15；所以有 z=18841，y=18841； temp=(long)18841； z=(18841L(15-10+9)/18841=3O8690944L/18841=16384； 因?yàn)樯蘻的定標(biāo)值為15，所以定點(diǎn)z=16384，即為浮點(diǎn)z=16384/215=0.5。 2.4程序變量的Q值確定 在前面幾節(jié)介紹的例子中，由于x，y，z的值都是已知的，因此從浮點(diǎn)變?yōu)槎c(diǎn)時(shí)Q值很好確定。在實(shí)際的DSP應(yīng)用中，程序中參與運(yùn)算的都是變量，那么如何確定浮點(diǎn)程序中變量的Q值呢

11、？從前面的分析可以知道，確定變量的Q值實(shí)際上就是確定變量的動(dòng)態(tài)范圍，動(dòng)態(tài)范圍確定了，則Q值也就確定了。 設(shè)變量的絕對(duì)值的最大值為|max|，注意|max|必須小于或等于32767。取一個(gè)整數(shù)n，使?jié)M足 2n-1|max|2n 則有 2-Q=2-15*2n=2-(15-n) Q=15-n 例如，某變量的值在-1至+1之間，即|max|1，因此n=0，Q=15-n=15。 既然確定了變量的|max|就可以確定其Q值，那么變量的|max|又是如何確定的呢？一般來(lái)說(shuō)，確定變量的|max|有兩種方法。一種是理論分析法，另一種是統(tǒng)計(jì)分析法。 1. 理論分析法 有些變量的動(dòng)態(tài)范圍通過(guò)理論分析是可以確定的。例

12、如： (1)三角函數(shù)。y=sin(x)或y=cos(x)，由三角函數(shù)知識(shí)可知，|y|=1。 (2)漢明窗。y(n)=0.54一0.46cosnn/(N-1)，0=n=N-1。因?yàn)?1=cos2n/(N-1)=1，所以0.08=y(n)=1.0。 (3)FIR卷積。y(n)=h(k)x(n-k)，設(shè)|h(k)|=1.0，且x(n)是模擬信號(hào)12位量化值，即有|x(n)|=211，則|y(n)|=211。 (4)理論已經(jīng)證明，在自相關(guān)線性預(yù)測(cè)編碼(LPC)的程序設(shè)計(jì)中，反射系數(shù)ki滿足下列不等式：|ki|15； for(i=0；i(n-1)；i)x1n-i-2=xinlength-i-1； 主程序

13、與浮點(diǎn)的完全一樣。“ 3 DSP定點(diǎn)算術(shù)運(yùn)算 定點(diǎn)DSP芯片的數(shù)值表示基于2的補(bǔ)碼表示形式。每個(gè)16位數(shù)用l個(gè)符號(hào)位、i個(gè)整數(shù)位和15-i個(gè)小數(shù)位來(lái)表示。因此： 00000010.10100000 表示的值為： 212-12-3=2.625 這個(gè)數(shù)可用Q8格式(8個(gè)小數(shù)位)來(lái)表示，其表示的數(shù)值范圍為-128至l27.996，一個(gè)Q8定點(diǎn)數(shù)的小數(shù)精度為1/256=0.004。 雖然特殊情況(如動(dòng)態(tài)范圍和精度要求)必須使用混合表示法。但是，更通常的是全部以Q15格式表示的小數(shù)或以Q0格式表示的整數(shù)來(lái)工作。這一點(diǎn)對(duì)于主要是乘法和累加的信號(hào)處理算法特別現(xiàn)實(shí)，小數(shù)乘以小數(shù)得小數(shù)，整數(shù)乘以整數(shù)得整數(shù)。當(dāng)然

14、，乘積累加時(shí)可能會(huì)出現(xiàn)溢出現(xiàn)象，在這種情況下，程序員應(yīng)當(dāng)了解數(shù)學(xué)里面的物理過(guò)程以注意可能的溢出情況。下面我們來(lái)討論乘法、加法和除法的DSP定點(diǎn)運(yùn)算，匯編程序以TMS320C25為例。 3.1定點(diǎn)乘法 兩個(gè)定點(diǎn)數(shù)相乘時(shí)可以分為下列三種情況： 1. 小數(shù)乘小數(shù) 例1.9 Q15*Q15=Q30 0.5*0.5=0.25 0.100000000000000；Q15 * 0.100000000000000；Q15 - 00.010000000000000000000000000000=0.25；Q30 兩個(gè)Q15的小數(shù)相乘后得到一個(gè)Q30的小數(shù)，即有兩個(gè)符號(hào)位。一般情況下相乘后得到的滿精度數(shù)不必全部保

15、留，而只需保留16位單精度數(shù)。由于相乘后得到的高16位不滿15位的小數(shù)據(jù)度，為了達(dá)到15位精度，可將乘積左移一位，下面是上述乘法的TMS320C25程序： LT OP1；OP1=4000H(0.5/Q15) MPY OP2；oP2=4000H(0.5/Ql5) PAC SACH ANS，1；ANS=2000H(0.25/Q15)二 整數(shù)乘整數(shù) 例1.10 Q0*Q0=Q0 17*(-5)=-85 0000000000010001=l7 *1111111111111011=-5-11111111111111111111111110101011=-85 3. 混合表示法 許多情況下，運(yùn)算過(guò)程中為了

16、既滿足數(shù)值的動(dòng)態(tài)范圍又保證一定的精度，就必須采用Q0與Q15之間的表示法。比如，數(shù)值1.2345，顯然Q15無(wú)法表示，而若用Q0表示，則最接近的數(shù)是1，精度無(wú)法保證。因此，數(shù)1.2345最佳的表示法是Q14。 例1.11 1.5*0.75= 1.125 01.10000000000000=1.5；Q14 *00.11000000000000=0.75；Q14 - 0001.0010000000000000000000000000=1.125 Q28 Q14的最大值不大于2，因此，兩個(gè)Q14數(shù)相乘得到的乘積不大于4。 一般地，若一個(gè)數(shù)的整數(shù)位為i位，小數(shù)位為j位，另一個(gè)數(shù)的整數(shù)位為m位，小數(shù)位為

17、n位，則這兩個(gè)數(shù)的乘積為(i+m)位整數(shù)位和(j+n)位小數(shù)位。這個(gè)乘積的最高16位可能的精度為(im)整數(shù)位和(15- i- m)小數(shù)位。 但是，若事先了解數(shù)的動(dòng)態(tài)范圍，就可以增加數(shù)的精度。例如，程序員了解到上述乘積不會(huì)大于1.8，就可以用Q14數(shù)表示乘積，而不是理論上的最佳情況Q13。例3.11的TMS320C25程序如下： LT OP1；OP1 = 6000H(1.5/Ql4) MPY OP2；OP2 = 3000H(0.75/Q14) PAC SACH ANS，1；ANS=2400H(1.125/Q13) 上述方法，為了精度均對(duì)乘的結(jié)果舍位，結(jié)果所產(chǎn)生的誤差相當(dāng)于減去一個(gè)LSB(最低位

18、)。采用下面簡(jiǎn)單的舍人方法，可使誤差減少二分之一。 LT OP1 MPY OP2 PAC ADD ONE，14(上舍入) SACH ANS，1 上述程序說(shuō)明，不管ANS為正或負(fù)，所產(chǎn)生的誤差是l/2 LSB，其中存儲(chǔ)單元ONE的值為1。 3.2定點(diǎn)加法 乘的過(guò)程中，程序員可不考慮溢出而只需調(diào)整運(yùn)算中的小數(shù)點(diǎn)。而加法則是一個(gè)更加復(fù)雜的過(guò)程。首先，加法運(yùn)算必須用相同的Q點(diǎn)表示，其次，程序員或者允許其結(jié)果有足夠的高位以適應(yīng)位的增長(zhǎng)，或者必須準(zhǔn)備解決溢出問(wèn)題。如果操作數(shù)僅為16位長(zhǎng)，其結(jié)果可用雙精度數(shù)表示。下面舉例說(shuō)明16位數(shù)相加的兩種途徑。 1.保留32位結(jié)果 LAC OP1；(Q15) ADD O

19、P2；(Ql5) SACH ANSHI ；(高16位結(jié)果) SACL ANSLO ：(低16位結(jié)果) 2.調(diào)整小數(shù)點(diǎn)保留16位結(jié)果 LAC OP1，15；(Q14數(shù)用ACCH表示) ADD OP2，15；(Q14數(shù)用ACCH表示) SACH ANS；(Q14) 加法運(yùn)算最可能出現(xiàn)的問(wèn)題是運(yùn)算結(jié)果溢出。TMS320提供了檢查溢出的專用指令BV，此外，使用溢出保護(hù)功能可使累加結(jié)果溢出時(shí)累加器飽和為最大的整數(shù)或負(fù)數(shù)。當(dāng)然，即使如此，運(yùn)算精度還是大大降低。因此，最好的方法是完全理解基本的物理過(guò)程并注意選擇數(shù)的表達(dá)方式。 3.3定點(diǎn)除法 在通用DSP芯片中，一般不提供單周期的除法指令，為此必須采用除法子

20、程序來(lái)實(shí)現(xiàn)。二進(jìn)制除法是乘法的逆運(yùn)算。乘法包括一系列的移位和加法，而除法可分解為一系列的減法和移位。下面我們來(lái)說(shuō)明除法的實(shí)現(xiàn)過(guò)程。 設(shè)累加器為8位，且除法運(yùn)算為10除以3。除的過(guò)程包括與被除法有關(guān)的除數(shù)逐步移位，在每一步進(jìn)行減法運(yùn)算，如果能減則將位插入商中。 (1)除數(shù)的最低有效位對(duì)齊被除數(shù)的最高有效位。 0000l0l0 - 00011000 - 11110010 (2)由于減法結(jié)果為負(fù)，放棄減法結(jié)果，將被除數(shù)左移一位，再減。 00010100 - 00011000 - 11111000 (3)結(jié)果仍為負(fù)，放棄減法結(jié)果，被除數(shù)左移一位，再減。 00101000 - 00011000 - 00

21、010000 (4)結(jié)果為正，將減法結(jié)果左移一位后加1，作最后一次減。 00100001 - 00011000 - 00001001 (5)結(jié)果為正，將結(jié)果左移一位加1 得最后結(jié)果。高4位代表余數(shù)，低4位表示商。 00010011 即，商為0011= 三 余數(shù)為0001= 1。 TMS320沒(méi)有專門的除法指令，但使用條件減指令SUBC可以完成有效靈活的除法功能。使用這一指令的唯一限制是兩個(gè)操作數(shù)必須為正。程序員必須事先了解其可能的運(yùn)算數(shù)的特性，如其商是否可以用小數(shù)表示及商的精度是否可被計(jì)算出來(lái)。這里每一種考慮可影響如何使用SUBC指令的問(wèn)題。下面我們給出兩種不同情況下的TMS320C25除法程

22、序。 (1)分子小于分母 DIV_A： LT NUMERA MPY DENOM PAC SACH TEMSGN；取商的符號(hào) LAC DENOM ABS SACL DENOM；使分母為正 ZALH NUMERA； 分子為正 ABS RPTK 14 SUBC DENOM；除循環(huán)15次 SACL QUOT LAC TEMSGN BGEZ A1；若符號(hào)為正，則完成 ZAC SUB QUOT SACL QUOT；若為負(fù)，則商為負(fù) A1： RET 這個(gè)程序中，分子在NUMERA中，分母在DENOM中，商存在QUOT中，TEMSGN為暫存單元。 (2)規(guī)定商的精度 DIV_B: LT NUMERA MPY

23、DENOM PAC SACH TEMSGN；取商的符號(hào) LAC DENOM ABS SACL DENOM; 使分母為正 LACK 15 ADD FRAC SACL FRAC；計(jì)算循環(huán)計(jì)數(shù)器 LAC NUMERA ABS ; 使分子為正 RPT FRAC SUBC DENOM; 除循環(huán)16FRAC次 SACL QUOT LAC TEMSGN BGEZ B1;若符號(hào)為正，則完成 ZAC SUB QUOT SACL QUOT；若為負(fù)，則商為負(fù) B1： RET 與DIV_A相同，這個(gè)程序中，分子在NUMERA中，分母在DENOM中，商存在QUOT中，TEMSGN為暫存單元。FRAC中規(guī)定商的精度，如商

24、的精度為Q13，則調(diào)用程序前FRAC單元中的值應(yīng)為13。 4 非線性運(yùn)算的定點(diǎn)快速實(shí)現(xiàn) 在數(shù)值運(yùn)算中，除基本的加減乘除運(yùn)算外，還有其它許多非線性運(yùn)算，如，對(duì)數(shù)運(yùn)算，開方運(yùn)算，指數(shù)運(yùn)算，三角函數(shù)運(yùn)算等，實(shí)現(xiàn)這些非線性運(yùn)算的方法一般有：(1)調(diào)用DSP編譯系統(tǒng)的庫(kù)函數(shù)；(2)查表法；(3)混合法。下面我們分別介紹這三種方法。 1.調(diào)用DSP編譯系統(tǒng)的庫(kù)函數(shù) TMS320C2X/C5X的C編譯器提供了比較豐富的運(yùn)行支持庫(kù)函數(shù)。在這些庫(kù)函數(shù)中，包含了諸如對(duì)數(shù)、開方、三角函數(shù)、指數(shù)等常用的非線性函數(shù)。在C程序中(也可在匯編程序中)只要采用與庫(kù)函數(shù)相同的變量定義，就可以直接調(diào)用。例如，在庫(kù)函數(shù)中，定義了以

25、10為底的常用對(duì)數(shù)log10()： #includemath.h double，log10(double x)； 在C程序中按如下方式調(diào)用： float x，y; X=10.0; y=log10(x)； 從上例可以看出，庫(kù)函數(shù)中的常用對(duì)數(shù)log10()要求的輸入值為浮點(diǎn)數(shù)，返回值也為浮點(diǎn)數(shù)，運(yùn)算的精度完全可以保證。直接調(diào)用庫(kù)函數(shù)非常方便，但由于運(yùn)算量大，很難在實(shí)時(shí)DSP中得到應(yīng)用。 2.查表法 在實(shí)時(shí)DSP應(yīng)用中實(shí)現(xiàn)非線性運(yùn)算，一般都采取適當(dāng)降低運(yùn)算精度來(lái)提高程序的運(yùn)算速度。查表法是快速實(shí)現(xiàn)非線性運(yùn)算最常用的方法。采用這種方法必須根據(jù)自變量的范圍和精度要求制作一張表格。顯然輸人的范圍越大，精度

26、要求越高，則所需的表格就越大，即存儲(chǔ)量也越大。查表法求值所需的計(jì)算就是根據(jù)輸入值確定表的地址，根據(jù)地址就可得到相應(yīng)的值，因而運(yùn)算量較小。查表法比較適合于非線性函數(shù)是周期函數(shù)或已知非線性函數(shù)輸入值范圍這兩種情況、例1.12和例1. 13分別說(shuō)明這兩種情況。 例1.12 已知正弦函數(shù)y=cos(x)，制作一個(gè)512點(diǎn)表格，并說(shuō)明查表方法。由于正弦函數(shù)是周期函數(shù)，函數(shù)值在-1至+1之間，用查表法比較合適。由于Q15的表示范圍為1-至32767/32768之間，原則上講-1至1的范圍必須用Q14表示。但一般從方便和總體精度考慮，類似情況仍用Q15表示，此時(shí)+1用32767來(lái)表示。 (1)產(chǎn)生5l2點(diǎn)值

27、的C語(yǔ)言程序如下所示。 #define N 512 #define pi 3.14l59 int sin_tab5l2； void main() int i; for(i=0；iN；i+)sin_tabi=(int)(32767*sin(2*pi*i/N)； (2)查表 查表實(shí)際上就是根據(jù)輸人值確定表的地址。設(shè)輸入x在02之間，則x對(duì)應(yīng)于512點(diǎn)表的地址為：index=(int)(512*x/2)，則y=sin(x)=sin_tabindex如果x用Q12定點(diǎn)數(shù)表示，將512/2用Q8表示為20861，則計(jì)算正弦表的地址的公式為。 index=(x*20861L)20； 例1.12用查表法求以

28、2為底的對(duì)數(shù)，已知自變量值范圍為0.5-1，要求將自變量范圍均勻劃分為10等分。試制作這個(gè)表格并說(shuō)明查表方法。 (1)作表： y=log2(x)，由于x在0.5到1之間，因此y在-1到0之間，x和y均可用Q15表示。由于對(duì)x均勻劃分為10段，因此，10段對(duì)應(yīng)于輸入x的范圍如表3.2所示。若每一段的對(duì)數(shù)值都取第一點(diǎn)的對(duì)數(shù)值，則表中第一段的對(duì)數(shù)值為y0(Q15)=(int)(log(O.5)*32768)，第二段的對(duì)數(shù)值為y1(Q15)=(int)(log2(0.55)*32768)，依次類推，如表3.2所示。 (2)查表： 查表時(shí)，先根據(jù)輸人值計(jì)算表的地址，計(jì)算方法為： index=(x-163

29、84)*20)15； 式中， index就是查表用的地址。例如，已知輸人x=26869，則index=6，因此，y= -10549。 表1.2 logtab0 10點(diǎn)對(duì)數(shù)表 地址 輸入值 對(duì)數(shù)值(Q15) 0 0.50-0.55 -32768 1 0.55-0.60 -28262 2 0.60-0.65 -24149 3 0.65-0.70 -20365 4 0.70-0.75 -16862 5 0.75-0.80 -13600 6 0.80-0.85 -10549 7 0.85-0.90 -7683 8 0.90-0.95 -4981 9 0.95-1.00 -2425 3.混合法 (1)提

30、高查表法的精度 上述方法查表所得結(jié)果的精度隨表的大小而變化，表越大，則精度越高，但存儲(chǔ)量也越大。當(dāng)系統(tǒng)的存儲(chǔ)量有限而精度要求也較高時(shí)，查表法就不太適合。那么能否在適當(dāng)增加運(yùn)算量的情況下提高非線性運(yùn)算的精度呢？下面介紹一種查表結(jié)合少量運(yùn)算來(lái)計(jì)算非線性函數(shù)的混合法，這種方法適用于在輸入變量的范圍內(nèi)函數(shù)呈單調(diào)變化的情形?；旌戏ㄊ窃诓楸淼幕A(chǔ)上來(lái)用計(jì)算的方法以提高當(dāng)輸入值處于表格兩點(diǎn)之間時(shí)的精度。提高精度的一個(gè)簡(jiǎn)便方法是采用折線近似法，如圖1.1所示。 圖1.1提高精度的折線近似法” 仍以求以2為底的對(duì)數(shù)為例(例1.12)。設(shè)輸入值為x，則精確的對(duì)數(shù)值為y，在表格值的兩點(diǎn)之間作一直線，用y作為y的近似

31、值，則有： y=y0y 其中y0由查表求得?，F(xiàn)在只需在查表求得y0的基礎(chǔ)上增加y既可。y的計(jì)算方法如下： y=(x/x0)y=x(y0/x0) 其中y0/x0對(duì)每一段來(lái)說(shuō)是一個(gè)恒定值，可作一個(gè)表格直接查得。此外計(jì)算此時(shí)需用到每段橫坐標(biāo)的起始值，這個(gè)值也可作一個(gè)表格。這佯共有三個(gè)大小均為10的表格，分別為存儲(chǔ)每段起點(diǎn)對(duì)數(shù)值的表logtab0、存儲(chǔ)每段y0/x0值的表logtab1和存儲(chǔ)每段輸入起始值x0的表logtab2，表logtab1和表logtab2可用下列兩個(gè)數(shù)組表示。 int logtab110=22529，20567，18920，17517，16308， 15255，14330，13

32、511，12780，12124；/*y0/x0：Q13*/ int logtab210=16384，18022，19660，21299，22938， 24576，26214，27853，29491，31130；/*x0：Q15*/ 綜上所述，采用混合法計(jì)算對(duì)數(shù)值的方法可歸納為： (1)根據(jù)輸人值，計(jì)算查表地址：index=(x-16384)*20)15； (2)查表得y0=logtab0index； (3)計(jì)算x=x-logtab2index； (4)計(jì)算y=(x*logtab1index)13； (5)計(jì)算得結(jié)果y=y0y。 例1.13已知x=0.54，求log2(x)。 0.54的精確對(duì)數(shù)

33、值為y=log2(0.54)=-0.889。 混合法求對(duì)數(shù)值的過(guò)程為： (1)定標(biāo)Q15，定標(biāo)值x=0.54*32768=17694； (2)表地址index=(x-16384)*20)15=0； (3)查表得y0=logtab00=-32768； (4)計(jì)算x=x-logtab20=17694-16384=1310； (5)計(jì)算y=(xlogtab1013=(13l0*22529L)13=3602 (6)計(jì)算結(jié)果y=y0y=-327683602=-29166。 結(jié)果y為Q15定標(biāo)，析算成浮點(diǎn)數(shù)為-29166/32768=-0.89，可見精度較高。 (2)擴(kuò)大自變量范圍 如上所述，查表法比較適

34、用于周期函數(shù)或自變量的動(dòng)態(tài)范圍不是太大的情形。對(duì)于像對(duì)數(shù)這樣的非線性函數(shù)，輸入值和函數(shù)值的變化范圍都很大。如果輸入值的變化范圍很大，則作表就比較困難。那么能否比較好地解決這個(gè)問(wèn)題，即不便表格太大，又能得到比較高的精度呢？下面我們來(lái)討論一種切實(shí)可行的方法。 設(shè)x是一個(gè)大于0.5的數(shù)，則x可以表示為下列形式： x=m*2e 式中，0.5=m 16; （2）無(wú)符號(hào)結(jié)果： unsign m1,m2,result; result=(unsigned long)m1*(unsigned long)m2) 16; TMS320C2X/C5X的C編譯器將浮點(diǎn)數(shù)表示為IEEE單精度格式。單精度數(shù)和雙精度數(shù)都表示

35、位32位，兩者沒(méi)有任何區(qū)別。在TMS320C2X/C5X的浮點(diǎn)庫(kù)中提供了一組浮點(diǎn)數(shù)學(xué)庫(kù)函數(shù)，如加法、減法。乘法、除法、比較、整數(shù)和浮點(diǎn)數(shù)轉(zhuǎn)換、標(biāo)準(zhǔn)的錯(cuò)誤處理等。這些函數(shù)也不遵循標(biāo)準(zhǔn)的C調(diào)用規(guī)則。調(diào)用這些函數(shù)時(shí)，C編譯器先將參數(shù)壓人運(yùn)行堆棧，然后調(diào)用浮點(diǎn)庫(kù)函數(shù)。函數(shù)執(zhí)行時(shí)，首先將參數(shù)從堆棧中彈出，然后執(zhí)行函數(shù)運(yùn)算，最后將運(yùn)算結(jié)果壓人堆棧返回。 有些浮點(diǎn)庫(kù)函數(shù)需要整型或長(zhǎng)整型參數(shù)或返回整型或長(zhǎng)整型結(jié)果，對(duì)這些函數(shù)，用累加器的低16位傳遞或返回整型數(shù)，而用累加器的所有32位傳遞或返回長(zhǎng)整型數(shù)。 六、TMS320C2X/C5X C語(yǔ)言程序開發(fā)舉例 本節(jié)我們以TMS320C2X為例，說(shuō)明定點(diǎn)DSP芯片C程序的開發(fā)過(guò)程。軟件開發(fā)過(guò)程與浮點(diǎn)DSP芯片的開發(fā)過(guò)程相類似，主要分以下幾個(gè)步驟： 1. 用編輯器（如EDIT、PE2等）編輯一個(gè)或多個(gè)C程序，如example1.c，example2.c。 2. 用一步編譯程序dspcl.exe對(duì)C程序編譯匯編形成目標(biāo)文件，如example1.obj，example2.obj： dspcl_v25_g_mn_o2 example1.c dspc_