用TMS320C54X實(shí)現Vertibi譯碼器

關(guān)鍵詞：Vertibi譯碼 TMS320C54X DSP 度量值 回溯 卷積編碼

引言

隨著(zhù)差錯控制編碼理論的完善和數字電路技術(shù)的發(fā)展，信道編碼已經(jīng)成功地應用于各種通信系統中。其基本做法是：在發(fā)送端將被傳輸的信息序更上附加一些監督碼元，這些多余的碼元與信息序列之間以某種確定的規則相互約束；接收端按照既定的規則檢驗信息碼元與監督碼元之間的關(guān)系，一旦傳輸過(guò)程中發(fā)生錯誤，則信息碼元與監督碼元之間的關(guān)系將受到破壞，從而發(fā)現錯誤，甚至糾正錯誤。按照信息碼元和監督碼元之間的約速方式，可以分為組碼和卷積碼。在GSM和IS-95中主要采用了卷積碼，在每三代移動(dòng)通信中，話(huà)音信號也主要采了卷積編碼。卷積碼的譯碼方法有許多種，其中最重要的就是Vertibi算法。為此，本文討論用DSP實(shí)現Vertibi譯碼器的方法。

1 卷積編碼器簡(jiǎn)介

卷積碼是一種對付突發(fā)錯誤的有效編碼方法，通常記作（n,k,N）。它將k個(gè)信息編碼為n個(gè)比特，編碼效率為Rc=k/n。N為約束長(cháng)度。與分組碼不同，卷積碼中編碼后的n個(gè)碼元不但與當前段的k個(gè)信息有關(guān)，而且與前面N-1段的信息有關(guān)，編碼過(guò)程中相互關(guān)聯(lián)的碼元為Nn個(gè)。其糾錯能力隨著(zhù)N的增加而增大，而差錯率隨著(zhù)N的增加而指數下降。卷積編碼器的結構如圖1所示。

由圖1可知，卷積編碼器包括兩部分：一個(gè)由N段組成的輸入移位寄存器，每段有k段，共有Nk位移位寄存器；n個(gè)模2和相加，其輸入分別對應于n個(gè)基于生成多項式的線(xiàn)性代數方程組。每輸入k個(gè)比特，編碼器輸出n個(gè)比特。在編碼器復雜度相同的情況下，卷積碼的性能優(yōu)于分組碼。

2 Vertibi譯碼的基本原理

卷積碼的譯碼方法主要包括Vertibi算法、Fano算法和堆棧算法等等，其中最重要的和最常用的就是Vertibi算法。Vertibi算法是一種關(guān)于解卷積的最大似然譯碼法。它不是在網(wǎng)格上依次比較所有的可能路徑，而是接收一段，計算一段，保留最有可能的路徑，從而達到整個(gè)碼序列是一個(gè)最大似然序列。

Vertibi算法可以算法描述如下：把在時(shí)刻i、狀態(tài)Sj所對應的網(wǎng)格圖節點(diǎn)記作（Sj，i），給每個(gè)網(wǎng)格圖節點(diǎn)賦值V（Sj，i）。節點(diǎn)值按照如下步驟計算：

①設V(S0,0)=0，i=1。

②在時(shí)刻i,對于進(jìn)入每個(gè)節點(diǎn)的所有路徑計算其不完全路徑的長(cháng)度。

③令V(Sj，i)為在i時(shí)刻，到達與狀態(tài)Sj相對應的節點(diǎn)（Sj,i）的最小不完全路徑長(cháng)度，通過(guò)在前一節點(diǎn)隨機選擇一條路徑就可產(chǎn)生新的結果，非存留支路將從網(wǎng)格圖中刪除。以這種方式，可以從（S0,0）處產(chǎn)生一組最小路徑。

④用L表示輸入編碼段的數目。其中，每段為k比特，m為編碼器中最大寄存器的長(cháng)度。如果iL+m，則令i=i+1，返回第②步。

一旦計算出所有節點(diǎn)值，則從i=L+m時(shí)刻，狀態(tài)S0開(kāi)始，沿網(wǎng)格圖中的幸存路徑反向追尋即可。這樣被定義支路與解碼輸出將是一一對應的。關(guān)于不完全路徑長(cháng)度，硬判斷決解碼采用漢明距離，軟判決解碼采用的是歐幾里德距離。軟判決的特性比硬判決要好2～3dB。

（n,k,N）卷積編碼器共有2 kn個(gè)狀態(tài)，因此Vertibi譯碼器必須具有同樣的2 kn個(gè)狀態(tài)，并且在譯碼過(guò)程中要存儲各狀態(tài)的幸存路徑和長(cháng)度。Vertibi譯碼器的復雜程度隨2 kn指數增加。一般要求n10。

3 Vertibi譯碼器的DSP實(shí)現

譯碼過(guò)程就是根據接收到的數據符號，按最大似然譯碼準則找出編碼器在網(wǎng)絡(luò )圖上所走過(guò)的路徑。Vertibi譯碼的處理過(guò)程如圖2所示。

3.1 度量值更新

度量值的更新包括以下4個(gè)步驟：

①計算每條可能輸入路徑的度量值；

②為每條支路計算總的距離；

③選擇保存最小度量值；

④保存幸存路徑。

Vertibi譯碼時(shí)，每收到一個(gè)符號就進(jìn)行狀態(tài)轉移，需要計算前一個(gè)狀態(tài)到各個(gè)新?tīng)顟B(tài)的分支度量值。我們用DSP設計的譯碼器采用軟判決輸入，度量值用歐氏距離表示。當編碼速率為1/C時(shí)，歐氏距離為：

其中SDn表示接收序列，Gn(j)為網(wǎng)格上每個(gè)路徑狀態(tài)的期望輸入值，j是路徑指示值，C為編碼速率的倒數。將（1）式展開(kāi)得：

對于所有的路徑來(lái)說(shuō)，都是一樣的，2只是個(gè)常數，在進(jìn)行各路徑度量值比較時(shí)，可以不考慮。這樣可簡(jiǎn)化為：

省去膠面的負號，在度量值的比較時(shí)取最大值。對于編碼速率為1/2的卷積碼，分支度量為：

T=SD0G0(j)+SD1G1(j) (4)

當編碼速率為1/3時(shí)，分支度量為：

T=SD0G0（j）+SD1G1(j)+SD2G2(j) （5）

Gn(j)用雙極性表示，即0用+1表示，1用-1表示。這樣分支度量值的計算可以進(jìn)一步簡(jiǎn)化為接收數據的加和減。

下面給出編碼速率為1/3時(shí)，DSP實(shí)現具體程序。

LD *AR1+，16，A ；A=SD(2*i)

ADD*AR1+,16,B,B ;B=SD（2*i）+SD（2*i+1）

ADD*AR1-，16，B，B ;B=SD(2*i)+SD(2*i+)+SD(2*i+2)

STH B,*AR2+ ;temp(0)=SD(2*i)+SD(2*i+2)

SUB*AR1+,16,A,B ;B=SD(2*i)-SD(2*i+)

ADD *AR1-,16,B,B ;B=SD(2*i)-SD(2*i+1)+SD(2*i+2)

STHB,*AR2+ ;temp(1)=SD(2*i)-SD(2*i+1)+SD（2*i+2）

SUB *AR1+，16，A，B ；B=SD(2*i)-SD(2*i+1)

SUB *AR1-,16,B,B ;B=SD(2*i)-SD(2*i+1)-SD(2*i+)

STH B,*AR2+ ;temp(2)=SD(2*i)-SD(2*i+1)-SD(2*i+2)

ADD *AR1+,16,A,B ;B=SD(2*i)+SD(2*i+1)

SUB *AR1+,16,B,B ;B=SD(2*i)+SD(2*i+1)-SD(2*i+2)

STH B,*AR2 ;temp(3)=SD(2*i)+SD(2*i+1)-SD(2*i+2)

加比選單元是Vertibi譯碼器的核心單元。它的主要功能是取出當前狀態(tài)的量度值，分別與其兩個(gè)后續支路的量度相加并比較，選擇羅小的一個(gè)作為后續狀態(tài)的量度，并保存幸存支路。圖3給出了該算法的示意圖。

C54X片內的比較、選擇和存儲單元（CSSU）就是專(zhuān)門(mén)為Viterbi算法設計的加法/比較/選擇（ACS）運算的硬件單元。圖3所示的運算包括加法、比較和選擇三部分操作。其加法運算由DSP的ALU完成。只要將狀態(tài)寄存器ST1中的C16位置成1，ALU就被配置成雙16位工作方式，這樣，就可以在一個(gè)機器周期內執行兩次加法運算。其結果（Old_Met1+D1和Old_Met2+D2）都是16位數，分別存放在累加器的高16位和低16位中。然后，利用CMPS指令對累加器的高16位和低16位進(jìn)行比較，并選擇出較大的一個(gè)數放到指令所指定的存儲單元中。在CMPS指令執行的過(guò)程中，狀態(tài)轉移寄存器TRN自動(dòng)記錄比較的結果，這一點(diǎn)非常有用。實(shí)現一個(gè)蝶式運算的程序如下：

LD *AR2，T ；T=本地距離

DADST *AR5，A ；A=Old_Met(2*j)+T//Old_Met(2*j+1)-T

CMPS A,*AR4+ ;New_Met(j)=(Max(Old_Met(2*j)+T,Old_Met(2*j+1)-T)

；TRN=RTN1

;若（Old_Met(2*j)+T=Old_Met(2*j+1)-T則

TRN[0]=1

CMPS B，*AR3+ ;New_Met(j+2 K-2)=(Max(Old_Met(2*j)-T,Old_Met(2*j+1)+T)

TRN=TRN1

；若（Old_Met(2*j)-T=Old_Met(2*j+1)+T）則TRN[0]=1

3.2 回溯

當接收完1幀數據后，添加尾比特，強迫網(wǎng)格圖的最后一個(gè)狀態(tài)（0狀態(tài)）開(kāi)始，反向追蹤最大似然路徑，完成原始數據的譯碼?；厮莸挠嬎氵^(guò)程如下：

①計算當前狀態(tài)在轉移數據塊中的數據；

②利用BITT指令從T寄存器中讀取當前狀態(tài)相應的轉稱(chēng)比特；

③用上面讀取的轉移比特更新?tīng)顟B(tài)。

溯的主要功能是，從每個(gè)符號間隔的轉移數據中提取正確的位。為了完成該功能，需要保存每個(gè)狀態(tài)相應的轉移比特。該比特表示選定碟形網(wǎng)絡(luò )的是上面分支還是下面分支，所有這些轉移比特構成一塊轉移字。表1給出轉移字塊中狀態(tài)數與相應轉移比特之間的關(guān)系。

表1 一個(gè)符號間隔內數據轉移的狀態(tài)順序

實(shí)現回溯的程序如下：

；以下程序中，累加器A存放狀態(tài)值，累加器B存放臨時(shí)的數據

；K為約束長(cháng)度，MASK=2 K-5-1，ONE=1

RSBX OVM ；關(guān)閉溢出模式

STM *NTRAN_END，AR2；轉移表的結束地址

STM#NBWORDS-1，AR1；要計算的輸出字的序號

MVMM AR1，AR4 ；復制AR1的數值到AR4中

STM #OUTPUT+NBWORDS-1，AR3

；輸出比特的地址指針

LD #0，A 初始0狀態(tài)存入累加器A中

STM#15，BRC ；do i=0,NBWORDS-1

BACK：RPTB TBEND-1 ；do j=0,15

；在轉移字中計算比特位置

SFTL A，-（K-2），B ；B=A>>(K-2)

AND ONE，B ；B=B1=msb of State

ADD A,1,B ;B=B+A1=2*State+msb of State

STLM B,T ;T=B(bit position)

；修正轉移字

SFTL A，-3，B ；B=A/8=State/8

AND MASK，B ；B=BMASK=(K-5)lsb'sof State/8

STLM AR0 ;AR0=轉稱(chēng)字索引

MAR *+AR2（-2K-5） ；修正寄存器值使其復位

MAR *AR2+0 ；加偏移量修正轉移字

BITT*AR2-0 ；測試轉移字中的比特位

ROLTC A

TBEND：STL A，*AR3-

BANZD：BACK，*AR1-

STM #15，BRC ；指向輸出緩沖區的首地址

LD*AR3，A ；將第一個(gè)字裝入累加器A中

RVS：SFTA A，-1，A

STM #15，BRC

RPTB RVS2-1

ROL B

SFTA A，-1，A

RVS2：BANZD RVS，*AR4- ；判斷所有的字是否都計算完

STL B，*AR3+ ；保存剛計算完的字

LD *AR3，A ；裝入下一個(gè)字