何秀慧 ，袁 琳

（1.浙江師范大學(xué) 行知學(xué)院，浙江 金華321004；2.天津光電通信技術(shù)有限公司，天津300211）

繼Turbo碼和LDPC碼之后，噴泉碼成為了近年來糾錯(cuò)碼領(lǐng)域新的研究熱點(diǎn)[1]。噴泉碼由于其無比率特性非常適用于時(shí)變信道。目前主要的兩類噴泉碼是LT碼[2]和 Raptor碼[3]。Raptor碼是 LT碼的一種擴(kuò)展，它由預(yù)編碼和LT碼兩部分組成。Shokrollahi等人設(shè)計(jì)的系統(tǒng)Raptor碼已被第三代蜂窩網(wǎng)絡(luò)多媒體廣播/多點(diǎn)傳送服務(wù)標(biāo)準(zhǔn)所采用。Raptor碼最初是為在刪除信道中應(yīng)用而提出的，現(xiàn)已證明在加性高斯白噪聲AWGN(Additive White Gaussian Noise)信道和衰落信道中Raptor碼同樣擁有接近信道容量的性能[4-6]。

[4]和參考文獻(xiàn)[5]中對 Raptor碼在 AWGN信道中采用二進(jìn)制相移鍵控BPSK(Binary Phase Shift Keying)調(diào)制的設(shè)計(jì)問題進(jìn)行了深入研究，給出了Raptor碼在AWGN信道中的置信傳播 BP（Belief Propagation）軟判決譯碼算法，其原理為先進(jìn)行LT碼譯碼再進(jìn)行預(yù)編碼譯碼。最小頻移鍵控MSK（Minimum Shift Keying）由于其屬于恒包絡(luò)調(diào)制和良好的頻譜特性已被廣泛應(yīng)用。本文研究了采用MSK調(diào)制的無線系統(tǒng)中Raptor碼譯碼算法的設(shè)計(jì)問題，給出了一種類似于Turbo碼譯碼的算法。仿真性能表明所提出的Raptor碼譯碼算法優(yōu)于參考文獻(xiàn)[4]中的譯碼算法。

1 Raptor碼及其BP譯碼算法

Raptor碼是一種級(jí)聯(lián)碼，其編碼過程由預(yù)編碼過程和LT碼編碼過程組成。預(yù)編碼通常選用高碼率的線性分組碼，本文中采用了LDPC碼作為預(yù)編碼。預(yù)編碼過程將輸入的信息比特通過傳統(tǒng)的編碼轉(zhuǎn)換為中間比特，中間比特再進(jìn)行編碼比特度分布函數(shù)為Ω(x)的LT碼編碼，從而得到Raptor碼的編碼比特。

相比于LT碼，Raptor碼中由于預(yù)編碼有一定的糾錯(cuò)能力，放寬了對編碼比特度分布函數(shù)Ω(x)的設(shè)計(jì)要求，從而整體降低了編譯碼的復(fù)雜度。圖1給出了Raptor碼的因子圖表示。

假定長度為 k的信息字 S=[s1，…，sk]采用 Raptor碼編碼器產(chǎn)生 n 個(gè)編碼比特 C=[c1，…，cn]，其中 ci∈{0，1}。γ=n/k稱為Raptor碼的譯碼開銷。編碼比特序列采用MSK調(diào)制產(chǎn)生調(diào)制符號(hào)序列 X=[x1，…，xn]，調(diào)制符號(hào)序列經(jīng)過無線信道傳輸，無線信道為均值為 0、方差為σ2的AWGN信道。接收信號(hào)序列為Y=[y1，…，yn]，其中，yi=xi+vi，vi(i=1，…，n)為高斯白噪聲的樣本值。 圖 2給出了本文構(gòu)建的Raptor碼編碼的無線系統(tǒng)圖示。

參考文獻(xiàn)[4]給出了Raptor碼的軟判決BP譯碼算法，其由LT碼的BP譯碼算法和LDPC預(yù)編碼的BP譯碼算法兩部分構(gòu)成。在接收端首先MSK解調(diào)器采用BCJR譯碼算法[7]獲得編碼比特的對數(shù)似然比LLR(Log-Likelihood)L(ci)，接著 LT碼的 BP譯碼器開始工作，用表示第l次迭代時(shí)從第n個(gè)中間比特向第m個(gè)編碼比特傳遞的LLR值消息，表示第l次迭代時(shí)從第m個(gè)編碼比特向第n個(gè)中間比特傳遞的LLR值消息，Ψ (n)表示中間比特或者編碼比特n的鄰點(diǎn)集合，則LT碼BP算法處理步驟如下：

在p次迭代譯碼后，第i個(gè)中間比特的LLR值消息為：

在 q次迭代譯碼后，第 i個(gè)信息比特(i=1，…，k)的LLR值消息為：

2 改進(jìn)的BP譯碼算法

MSK調(diào)制可看成是一種編碼調(diào)制方案，其由二進(jìn)制連續(xù)相位編碼器和一個(gè)無記憶映射器構(gòu)成[8]。在Raptor碼編碼的無線系統(tǒng)中，可以在MSK解調(diào)器、LT碼譯碼器以及預(yù)編碼譯碼器之間進(jìn)行LLR值消息交換，從而達(dá)到進(jìn)一步改善譯碼性能的目的。

具體譯碼過程描述如下，在第l次譯碼迭代時(shí)，首先MSK解調(diào)器采用BCJR譯碼算法獲得編碼比特的LLR值消息：

接著運(yùn)行LT碼BP譯碼處理過程：

LT碼譯碼器提供給MSK解調(diào)器的外部LLR值消息為：

其次運(yùn)行LDPC碼的BP譯碼器：)

則LDPC碼譯碼器提供給LT碼譯碼器的外部LLR值消息為：

在q次迭代譯碼后，第i個(gè)信息比特的LLR值為：

3 仿真實(shí)驗(yàn)

在Matlab軟件中構(gòu)造Raptor碼編碼的無線系統(tǒng)仿真模型。AWGN信道方差取值為1，Raptor碼采用碼率為0.9的規(guī)則LDPC碼作為預(yù)編碼，其變量節(jié)點(diǎn)的度數(shù)為3，校驗(yàn)節(jié)點(diǎn)的度數(shù)為30。長度為9 000的信息字用PEG算法構(gòu)造[9]的LDPC碼進(jìn)行預(yù)編碼產(chǎn)生中間比特，LT碼的編碼比特節(jié)點(diǎn)度分布函數(shù)如下所示[4]：

Raptor碼譯碼器分別采用了參考文獻(xiàn)[4]中BP譯碼算法和本文提出的改進(jìn)BP譯碼算法。參考文獻(xiàn)[4]BP譯碼算法中LT碼譯碼器和LDPC碼譯碼器各執(zhí)行80次迭代譯碼。改進(jìn)BP譯碼算法也進(jìn)行了80次迭代譯碼。圖3給出了在不同譯碼開銷γ下兩種譯碼算法的誤碼率，由圖可知，本文所提算法能獲得更好的譯碼性能。

本文研究了在AWGN信道中傳輸Raptor碼的譯碼器設(shè)計(jì)問題，針對MSK調(diào)制和Raptor碼的特點(diǎn)提出了一種改進(jìn)的Raptor碼譯碼算法。仿真結(jié)果證明了算法的優(yōu)越性，為Raptor碼在無線通信中的應(yīng)用提供了新的依據(jù)。

參考文獻(xiàn)

[1]穆建君,焦曉,曹訓(xùn)志.數(shù)字噴泉碼及其應(yīng)用的研究與展望[J].電子學(xué)報(bào),2009,37(7):1571-1577.

[2]LUBY M.LT codes[C].Proc 43ed Ann IEEE Symp on Foundations of Computer Science,Vancouver,BC,Canada,2002:271-282.

[3]SHOKROLLAHI A.Raptor codes[J].IEEE Trans Inform Theory,2006,52(6):2551-2567.

[4]ETESAMI O,SHOKROLLAHI A.Raptor codes on binary memoryless symmetric channels[J].IEEE Trans Inform Theory,2006.52(5):2033-2051.

[5]CHENG Z,CASTURA J,MAO Y.On the design of raptor codes for binary-input Gaussian cahnnels[J].IEEE Trans Commun,2009,57(11):3269-3277.

[6]袁磊,安建平,李祥明.噴泉碼在無線中繼網(wǎng)絡(luò)中的應(yīng)用[J].信息通信技術(shù),2009,12(6):61-65.

[7]BAHL L,COCKE J,JELINEK F,et al.Optimal decoding of linear codes for minimizing symbol error rate[J].IEEE Trans Inform Theory，1974,20(2):284-287.

[8]RIMOLDI B.A decomposition approach to CPM[J].IEEE Trans Inform Theory，1988,34(2):260-270.

[9]HU X,ELEFTHERIOU E,ARNOLD D.Regular and irregular progressive edge-growth tanner graphs[J].IEEE Trans Inform Theory，2005,51(1):386-398.