林梅英，許肖梅，陳友淦，張 蘭

(1. 廈門大學(xué)海洋與地球?qū)W院，福建廈門 361102；2. 廈門大學(xué)水聲通信與海洋信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，福建廈門 361005)

0 引 言

水聲信道是一種極其復(fù)雜的隨機(jī)時(shí)-空-頻變信道，表現(xiàn)為強(qiáng)多途、大起伏、窄帶寬及高噪聲干擾背景。高傳輸速率、低誤碼率水聲通信技術(shù)是目前水聲通信(Underwater Acoustic Communication,UWA)研究的重點(diǎn)[1]。要達(dá)到低誤碼率的水聲通信性能，必須采用信道編碼技術(shù)，并要求盡可能采用較短碼長(zhǎng)的碼。

卷積碼、RS碼、低密度奇偶校驗(yàn)(Low-Density Parity-Check, LDPC)碼、Turbo 碼等編碼方案，已在水聲通信中廣泛使用。LDPC碼在這幾種碼中性能最接近香農(nóng)限，且譯碼算法復(fù)雜度較低[2]。而水聲信道要求盡可能用較短碼長(zhǎng)、易于實(shí)時(shí)處理的編碼技術(shù)來(lái)提高通信系統(tǒng)的可靠性，準(zhǔn)循環(huán)低密度奇偶校驗(yàn)(Quasi-Cyclic Low Density Parity Check,QC-LDPC)碼在短碼時(shí)糾錯(cuò)能力比 LDPC碼強(qiáng)，復(fù)雜性更低，是目前水聲信道編碼領(lǐng)域研究的熱點(diǎn)[3]。文獻(xiàn)[3]研究了QC-LDPC編解碼技術(shù)，碼率為1/2，構(gòu)建的LDPC-OFDM水聲通信系統(tǒng)，在某港淺海域取得了不錯(cuò)的效果。

但水聲信道快速時(shí)-頻-幅變化的特征，要求信道編碼的碼率甚至幀長(zhǎng)能夠自適應(yīng)地根據(jù)信道環(huán)境變化實(shí)時(shí)做出相應(yīng)編碼調(diào)整，因而研究碼率兼容編碼就顯得尤為重要。目前碼率兼容碼在水聲信道中的應(yīng)用還較少，文獻(xiàn)[4]在自動(dòng)重傳請(qǐng)求(Automatic Repeat-reQuest, ARQ)淺海水聲系統(tǒng)中采用信息位不變，擴(kuò)展校驗(yàn)位的編碼方法構(gòu)造碼率兼容LDPC(Rate-Compatible LDPC, RC-LDPC)碼。由于QC-LDPC碼的特定結(jié)構(gòu)，本文采用另一種簡(jiǎn)單的擴(kuò)展方法：固定校驗(yàn)位，改變信息位的位數(shù)[5]，以構(gòu)建碼率為1/2、2/3… 5/6等一系列的碼率兼容QC-LDPC碼(RC QC-LDPC)，并在聲速為常數(shù)的均勻介質(zhì)水聲信道、負(fù)聲速梯度水聲信道及5徑淺海水聲(Shallow Water Acoustic, SWA)信道中研究其性能特點(diǎn)，為水聲通信系統(tǒng)參數(shù)的選擇提供參考。

1 系統(tǒng)的仿真模型

水聲自適應(yīng)信道編碼系統(tǒng)模型如圖1所示。仿真中，發(fā)射端發(fā)射信息序列，經(jīng)由某一碼率的 RC QC-LDPC編碼器及BPSK調(diào)制后，進(jìn)入水聲信道模型中傳輸。接收方首先進(jìn)行BPSK解調(diào)，然后進(jìn)行該碼率RC QC-LDPC譯碼，最后對(duì)比信源和信宿計(jì)算誤碼率，由誤碼率的高低來(lái)調(diào)整碼率選擇。

圖1 水聲自適應(yīng)信道編碼通信原理圖Fig.1 Block diagram of the adaptive coding system for the frequency hopping communication in underwater acoustic channel

2 水聲信道模型

實(shí)際水聲信道的傳輸函數(shù)是時(shí)-頻-空變的，不同海域的信道狀況會(huì)有很大不同。但在相干時(shí)間范圍內(nèi)，可認(rèn)為其轉(zhuǎn)移函數(shù)不隨時(shí)間變化。本文將針對(duì)聲速為常數(shù)的等溫層、負(fù)聲速梯度信道模型和典型 5 徑淺海水聲信道模型進(jìn)行RC QC-LDPC碼性能的仿真研究。

Np徑淺海水聲信道一般可以平均幅值αp和時(shí)延tp為特征值，平均幅值αp可由式(1)給出[6]：

而第p途徑的延時(shí)tp=lp/c，c=1500 m/s是水下聲速。所以，Np徑SWA信道的沖擊響應(yīng)為

本文根據(jù)該信道建模方法結(jié)合海域?qū)崪y(cè)參數(shù)進(jìn)行實(shí)際水聲信道的建模工作。在假設(shè)通信帶寬為5 kHz，采樣頻率為10 kHz的條件下，艾宇慧[7]等給出了以下兩種典型海洋水聲信道的傳輸函數(shù)，即：

聲速為常數(shù)的均勻介質(zhì)信道(Invariable Sound Velocity Gradient, ISVG)傳輸函數(shù)為

負(fù)聲速梯度信道(Negative Sound Velocity Gradient, NSVG)傳輸函數(shù)為

而典型5徑淺海水聲信道則采用文獻(xiàn)[6]中典型的準(zhǔn)靜止衰落信道模型進(jìn)行建模?？紤]一個(gè)距離為3 km、水深為75 m的信道模型，假設(shè)發(fā)射端和接收端接近海底，仿真建立5徑淺海水聲信道模型。表1給出了信道每徑相對(duì)時(shí)延和幅度的計(jì)算結(jié)果。調(diào)制方式為BPSK，符號(hào)周期T=2.5 ms。

表1 淺海水聲信道仿真模型Table 1 The profile of shallow water acoustic channel for simulations

3 RC QC-LDPC碼

首先構(gòu)造出一個(gè)低碼率的QC-LDPC碼的校驗(yàn)矩陣(H矩陣)。為了便于編碼，該矩陣的校驗(yàn)部分采用雙對(duì)角結(jié)構(gòu)，信息部分采用稀疏矩陣，通過(guò)增加部分信息位實(shí)現(xiàn)碼率由低到高的變化。利用該方法構(gòu)造的多碼率兼容的QC-LDPC碼，編碼時(shí)只需要增加部分信息位的計(jì)算即可，編碼簡(jiǎn)單且易于多碼率復(fù)用。而譯碼時(shí)以碼率最低的一個(gè)碼字所占的資源為主，只是需要改變部分參數(shù)就可以實(shí)現(xiàn)多碼率兼容的 QC-LDPC碼的譯碼[5]。該方法擴(kuò)展校驗(yàn)矩陣的結(jié)構(gòu)如圖2所示。

圖2 擴(kuò)展信息位方法示意圖Fig.2 Schematic diagram of extended information bits

對(duì)于該類型碼率兼容QC-LDPC碼，本文中設(shè)計(jì)的碼字參數(shù)如表2所示。

表2 信息位擴(kuò)展RC QC-LDPC碼Table 2 RC QC-LDPC codes with information bit expansion

4 仿真分析

仿真實(shí)驗(yàn)中設(shè)計(jì)一組RC QC-LDPC碼：母碼碼率為 1/2，由擴(kuò)展信息位生成碼率分別為 1/2、2/3、3/4、4/5、5/6等5種碼率的RC QC-LDPC碼系列，具體碼長(zhǎng)及信息位長(zhǎng)度見(jiàn)表2；譯碼采用BP方法，迭代次數(shù)為20次；仿真采用Matlab平臺(tái)，仿真每次發(fā)送 20幀數(shù)據(jù)。為研究水聲通信系統(tǒng)中RC QC-LDPC碼在不同信道特征下不同碼率性能，分別在聲速為常數(shù)的ISVG信道、NSVG信道、典型 5 徑淺海水聲信道中仿真研究了RC QC-LDPC碼的性能。假定水聲通信的性能指標(biāo)為誤碼率(BER)低于10-4。仿真結(jié)果如圖3、4、5所示。

由圖3、4可以看出，RC QC-LDPC碼能顯著提高水聲通信系統(tǒng)性能，且碼率越低性能越好，所需要的信噪比也越低。由圖3、4也可發(fā)現(xiàn)，RC QC-LDPC碼在ISVG中的性能比在NSVG差一些，其原因是文獻(xiàn)[7]作者在進(jìn)行水聲信道仿真時(shí)，接收水聽(tīng)器的布放位置靠近海底，從理論上可知，在負(fù)聲速梯度下，聲線往聲速變小的方向即海底方向走，偏下方，導(dǎo)致通信性能比較好。

圖3 RC QC-LDPC編碼在ISVG水聲信道中的BER曲線Fig.3 The BER performance of RC QC-LDPC over ISVG channel

圖4 RC QC-LDPC編碼在NSVG水聲信道中的BER曲線Fig.4 The BER performance of RC QC-LDPC over NSVG channel

圖5 RC QC-LDPC編碼在典型5徑淺海水聲信道中的BER曲線Fig.5 The BER performance of RC QC-LDPC over SWA channel

5徑淺海信道的仿真見(jiàn)圖5。同時(shí)，由于淺海典型5徑信道傳播時(shí)延大，多途干擾強(qiáng)，必須在信道譯碼前加入均衡，才能保證可靠的水聲通信性能(可由圖5中*形曲線看出)。本文自適應(yīng)判決反饋均衡器采用RLS算法，其收斂速度快、適于跟蹤快速變化信道[8]。

RLS算法迭代公式如下[8]：

其中：d(n)為M階自適應(yīng)濾波器的期望輸出；x(n)為n時(shí)刻的輸人向量；ω(n) 為此時(shí)刻濾波器權(quán)系數(shù)向量；e(n)為此時(shí)刻的誤差；K(n)是增益矢量;R(n)是誤差相關(guān)矩陣。仿真中 RLS的參數(shù)設(shè)置為λ= 0 .9999，初始值R(n) = 0 .05×eye(M,M)(eye(n)是產(chǎn)生一個(gè)nxn大小的單位矩陣)，M是均衡器階數(shù)，設(shè)置訓(xùn)練長(zhǎng)度為512 bit，抽頭系數(shù)的個(gè)數(shù)為45。

由圖5可知，所設(shè)計(jì)的RC QC-LDPC碼通過(guò)均衡消除多途帶來(lái)的碼間干擾后，在淺海水聲信道中也是可行的。同樣，在信噪比較高的水聲環(huán)境中可選擇碼率較高的QC-LDPC碼以避免編碼造成冗余數(shù)據(jù)。水聲信道時(shí)延越長(zhǎng)、徑數(shù)越多，滿足水聲通信指標(biāo)所需要的編碼碼率就越低。本文所設(shè)計(jì)的RC QC-LDPC碼在水聲信道中具有很好的適用性，該編碼方案是可行有效的。

為進(jìn)一步明了 RC QC-LDPC碼在淺海水聲通信中的參數(shù)選擇問(wèn)題，本文建立了不同信噪比下滿足通信性能指標(biāo)的編碼碼率的查找表(Look-Up Table, LUT)，以供實(shí)際通信中編碼碼率的選擇。假設(shè)系統(tǒng)性能的要求為誤碼率低于10-4，表3給出了在典型 5徑水聲信道中由不同接收端信噪比(Eb/N0)表征的不同水聲信道的查找表。當(dāng)滿足性能指標(biāo)時(shí)，標(biāo)記為“1”；反之，標(biāo)記為“0”。在使用該表時(shí)，如在典型5徑淺海水聲信道中，當(dāng)估計(jì)到接收端的Eb/N0=9 dB，為了達(dá)到檢測(cè)門限的要求，可選擇的碼率為 1/2和 2/3，最高碼率為1/2(其他碼率不能滿足性能指標(biāo))。因此相對(duì)其它碼率的備選方案，在本次通信傳輸中，系統(tǒng)將集中選擇1/2碼率和2/3碼率這兩種碼，直到獲得更新的Eb/N0估計(jì)結(jié)果。以此類推其他Eb/N0下的情況。在其他水聲信道中也是通過(guò)這樣的方法確定碼率。需要注意的是，在淺海水聲信道，RC QC-LDPC碼的性能是指有經(jīng)過(guò)RLS信道均衡之后的結(jié)果，選擇不同的均衡算法，結(jié)果會(huì)有差異。

表3 淺海水聲信道中RC QC-LDPC 碼性能查找表Table 3 The performance LUT of RC QC-LDPC codes in shallow water acoustic channel

5 結(jié) 論

為保證水聲數(shù)據(jù)可靠傳輸，提高幅-時(shí)-頻-空域動(dòng)態(tài)變化范圍大的水聲信道的信道帶寬、信噪比等綜合利用效率，本文構(gòu)建RC QC-LDPC碼的水聲通信系統(tǒng)。結(jié)合水聲信道特點(diǎn)，研究了 RC QC-LDPC碼在三種不同水聲信道(ISVG、NSVG、SWA)中不同碼率的性能。

仿真結(jié)果表明，RC QC-LDPC碼能有效降低水聲通信的誤碼率，提高水聲通信系統(tǒng)的性能；且碼率越低性能越好，所需要的信噪比也越低。所設(shè)計(jì)的RC QC-LDPC碼結(jié)構(gòu)簡(jiǎn)單、便于實(shí)現(xiàn)，編譯碼靈活，提高了水聲信道的信道利用率，是水聲自適應(yīng)通信中一種不錯(cuò)的選擇。

參考文獻(xiàn)

[1] Stojanovic M. Retrofocusing techniques for high rate acoustic communications[J]. Acoustical Society of America, 2005, 117(3):1173-1185.

[2] Chung S Y, Forney G D, Richardson T J, et al. On the design of low-density parity-check codes within 0.0045 dB of the shannon limit[J]. IEEE Communications Letters, 2001, 5(2): 58-60.

[3] 白櫟旸. LDPC碼及其在淺海OFDM水聲通信系統(tǒng)中的應(yīng)用研究和DSP實(shí)現(xiàn)[D]. 廈門: 廈門大學(xué), 2011.BAI Liyang. The research and DSP implementation of LDPC codes in OFDM underwater acoustic communication system[D].Xiamen: Xiamen University, 2011.

[4] CHEN Y, XU X, ZHANG L, et al. Design and Application of dynamic coding in shallow water acoustic communications[C]//OCEANS, 2012-Yeosu. IEEE, 2012: 1-6.

[5] 李宗臣. 碼率兼容 QC-LDPC 碼的設(shè)計(jì)和實(shí)現(xiàn)[D]. 西安: 西安電子科技大學(xué), 2011.LI Zongcheng. Design and implementation of rate-compatable QC-LDPC codes[D]. Xi’an: Xidian University, 2011.

[6] Vajapeyam M, Vedantam S, Mitra U, et al. Distributed space–time cooperative schemes for underwater acoustic communications[J].IEEE Journal of Oceanic Engineering, 2008, 33(4): 489-501.

[7] 艾宇慧, 高靜. 水聲信道相關(guān)均衡器仿真研究[J]. 聲學(xué)學(xué)報(bào), 1999,24(6): 589-597.AI Yuhui, GAO Jing. A simulation study on correlation equalization of underwater acoustic channels[J]. Acta Acustica, 1999,24(6): 589-597.

[8] 陳海蘭, 胡曉毅, 許茹, 等. LMS算法和RLS算法在水聲信道通信系統(tǒng)中應(yīng)用的比較[J]. 無(wú)線通信, 2007, 242(3): 18-19.CHEN Hailan, HU Xiaoyi, XU Ru, et al. Comparison of LMS algorithm and RLS algorithm using in the underwater acoustic communication system[J]. Modern Electronics Technique, 2007,242(3), 18-19.