羅亞松許江湖胡洪寧賀靜波陳占偉②（海軍工程大學(xué)電子工程學(xué)院 武漢 430033）

②（海軍青島水警區(qū)后勤部艦艇教練室 青島 266102）

正交頻分復(fù)用傳輸速率最大化自適應(yīng)水聲通信算法研究

羅亞松*①許江湖①胡洪寧①賀靜波①陳占偉②①
①（海軍工程大學(xué)電子工程學(xué)院 武漢 430033）

②（海軍青島水警區(qū)后勤部艦艇教練室 青島 266102）

該文在對正交頻分復(fù)用（OFDM）水聲通信速率最大化問題建?；A(chǔ)上，對注水算法的運(yùn)算量進(jìn)行了合理簡化，實(shí)現(xiàn)了對子載波初始能量的快速分配。針對固定判決門限離散比特分配算法存在的碼元能量嚴(yán)重浪費(fèi)和速率損失問題，提出了基于貪婪的比特分配方案，以較小的計(jì)算量代價(jià)實(shí)現(xiàn)了OFDM通信速率的顯著提高。研究表明，在壓縮信道編碼的不完全信道估計(jì)下，新算法仍具有良好的通信速率及誤碼率性能，適用于時(shí)變水聲信道。

水聲通信；正交頻分復(fù)用（OFDM）；自適應(yīng)；注水算法

1 引言

正交頻分復(fù)用（OFDM）水聲通信技術(shù)在水下系統(tǒng)中得到了廣泛重視與應(yīng)用［1］，基于信道估計(jì)的OFDM 自適應(yīng)調(diào)制技術(shù)已成為惡劣水聲信道環(huán)境中實(shí)現(xiàn)有效數(shù)據(jù)傳輸?shù)闹攸c(diǎn)研究方向［2］。隨著人類水下探索活動(dòng)日益增多，對水聲通信速率的需求也不斷增大，如何充分利用水聲信道資源，在有限的水聲信道帶寬內(nèi)實(shí)現(xiàn) OFDM水聲通信系統(tǒng)的速率最大化，是OFDM水聲通信領(lǐng)域值得關(guān)注的一項(xiàng)重要問題。

本文在對 OFDM水聲通信碼元速率最大化問題的建模分析基礎(chǔ)上，針對注水算法［3］運(yùn)算量較大的問題，提出了簡化計(jì)算量的方法；針對文獻(xiàn)［4］固定判決門限離散比特分配算法所帶來的碼元能量浪費(fèi)及速率損失問題，提出了進(jìn)一步提高碼元速率的二次貪婪比特分配新算法，并對新算法在碼元速率、計(jì)算量、不完全水聲信道估計(jì)下的誤碼性能進(jìn)行了綜合分析與仿真驗(yàn)證。

2 OFDM碼元速率最大化問題建模

OFDM是一種在多個(gè)正交子載波上完成并行數(shù)據(jù)傳輸?shù)耐ㄐ偶夹g(shù)，每個(gè)子載波上傳遞的信息以該子載波上的頻域值 X（k）形式出現(xiàn)，利用IFFT可以快速實(shí)現(xiàn)OFDM信源信息的調(diào)制［5］，OFDM時(shí)域通信信號為

其中，V是以頻率 fs抽樣的IFFT運(yùn)算點(diǎn)個(gè)數(shù)，k是子載波序號，X （k）是第k個(gè)子載波上傳遞的信息值，第k個(gè)子載波的對應(yīng)頻率為

其中， H（k）是水聲信道在子載波k上的信道響應(yīng)，N（k）是子載波k的信道噪聲。根據(jù)香農(nóng)定理，第k個(gè)子載波能獲得的最大通信傳輸速率由式（4）決定，SNRk是子載波k上的等效信噪比，pk是子載波k上的發(fā)射能量，是噪聲 N（k）的能量，為方便說明，將記做 ek，ek可通過現(xiàn)有成熟的信道估計(jì)算法及噪聲估計(jì)算法得到［6-8］， ek越大，說明該子載波上的通信環(huán)境越好，反之，越惡劣。

在水聲信道中，環(huán)境噪聲由地殼震動(dòng)噪聲、航船噪聲、風(fēng)浪噪聲及熱噪聲等多種因素構(gòu)成，子載波k上的噪聲大小可用式（5）衡量［9］，其中，fkL和 fkH分別是子載波k所占帶寬的最低和最高頻率。可見，水聲環(huán)境噪聲具有一定的頻率相關(guān)特性；同時(shí)，由于水聲信道多徑效應(yīng)非常嚴(yán)重，在頻域上即體現(xiàn)為顯著的頻率選擇性衰落［10］，這些都導(dǎo)致了不同子載波信道上通信環(huán)境的巨大差異，即 ek的劇烈起伏。

為各個(gè)子載波平均分配能量 pk的方法并沒有充分考慮這種載波信道差異，會(huì)造成速率-能耗比的低效，為了實(shí)現(xiàn)OFDM碼元通信速率的最大化，需要根據(jù)當(dāng)前水聲信道的估計(jì)狀態(tài)，將有限的碼元發(fā)射能量P（現(xiàn)在多數(shù)水下應(yīng)用系統(tǒng)都是能量受限系統(tǒng)）合理分配給所使用的N個(gè)子載波，碼元速率最大化問題即演變?yōu)楣剑?）所示的條件約束下的函數(shù)極值問題， k0是OFDM占用的初始子載波序號。

3 基于固定門限比特分配的速率最大化算法

為實(shí)現(xiàn)式（6）的最優(yōu)求解，由拉格朗日乘數(shù)法，可定義目標(biāo)函數(shù)為

對式（7）關(guān)于 pk求導(dǎo)，可得

為了使式（6）的目標(biāo)函數(shù)最大化， pk應(yīng)滿足：

其中，W是常數(shù)，根據(jù)式（6）的碼元能量約束條件，對式（9）關(guān)于各個(gè)子載波k求和，可得

由式（9），式（10），利用注水算法即可實(shí)現(xiàn)各子載波分配能量 pk的求解，若 pk≤ 0，說明子載波k的信道條件過于惡劣，不適合傳遞數(shù)據(jù)，這時(shí)實(shí)際使用的子載波數(shù)N應(yīng)減 1，重新利用式（9），式（10）迭代計(jì)算，直到所有 pk＞ 0為止，即可完成子載波能量的優(yōu)化分配，最終使用的子載波數(shù)為N'。

上述注水方法計(jì)算得到的各子載波能量 pk可實(shí)現(xiàn)OFDM碼元速率的最大化，但是，由此計(jì)算得到的每個(gè)子載波最大承載比特?cái)?shù)并非整數(shù)。因此，在實(shí)際應(yīng)用中，需要在注水方法計(jì)算的子載波能量基礎(chǔ)上，對各個(gè)子載波傳遞的比特?cái)?shù)進(jìn)行離散化分配。文獻(xiàn)［4］給出了一種基于固定判決門限的離散比特分配方法，其思路是，假設(shè)各子載波均采用矩形MQAM調(diào)制方式傳遞信息［11］（模數(shù)M可變），根據(jù)水聲通信應(yīng)用所要求的誤碼率指標(biāo) Pe*，由式（11）反推不同模數(shù)M下，滿足 Pe*要求所需的最小等效信噪比SNRM。由于OFDM的通信誤碼率是由最差子載波上的誤碼率決定［12］，因此，為了達(dá)到誤碼率小于的指標(biāo)要求，每個(gè)子載波k上的誤碼率都不應(yīng)超過。根據(jù)上述注水算法得到的分配能量，可計(jì)算子載波k上的信噪比SNRk將SNRk與不同模數(shù)M的SNRM（固定判決門限）進(jìn)行比較，為實(shí)現(xiàn)最大的碼元速率，選擇滿足的最大模數(shù)M，作為子載波k的調(diào)制模式，對應(yīng)子載波k上的比特?cái)?shù)為完成比特離散分配。

4 改進(jìn)OFDM碼元速率優(yōu)化算法

注水算法可根據(jù)當(dāng)前水聲信道狀況 ek的估計(jì)結(jié)果，實(shí)現(xiàn)碼元速率最大化的子載波能量分配，但是，一旦出現(xiàn)某個(gè)子載波的分配能量 pk≤ 0，就需要啟動(dòng)新一輪迭代，并且每次迭代都要重新計(jì)算各子載波上的發(fā)射能量。當(dāng)水聲信道頻域起伏較大，碼元信噪比較低時(shí)，注水算法的計(jì)算復(fù)雜度為 O（N2），這對于具有較大載波數(shù)N的 OFDM水聲通信系統(tǒng)來說，將會(huì)造成巨大的運(yùn)算壓力，不但影響水聲通信實(shí)時(shí)性，降低帶寬利用率，同時(shí)，由于水聲信道具有時(shí)、空變特性，過長的計(jì)算時(shí)間也會(huì)影響能量分配結(jié)果的時(shí)效性。

研究發(fā)現(xiàn)，碼元速率最大化注水算法傾向于為信道條件良好（ek較大）的子載波分配能量，而對于信道惡劣的子載波，分配較少、甚至不分配能量。經(jīng)過第 1輪注水算法迭代計(jì)算出的各子載波能量pk，雖然存在不合情理的負(fù)值，但是分配結(jié)果已經(jīng)能夠較準(zhǔn)確反映注水算法對各個(gè)子載波的分配傾向。這時(shí)可以利用第1輪迭代保留的N+個(gè)子載波正值能量 p+（k _i是保留使用的子載波序號），由式

k_i（12），對碼元能量P進(jìn)行 2次比例分配，可見，此時(shí)既濾除了那些信道條件惡劣、不適宜傳輸?shù)淖虞d波（第1輪分配能量小于0），又在碼元能量約束條件下維持了精確注水算法的分配傾向，同時(shí)將計(jì)算復(fù)雜度降為 O（N ）。

上述簡化注水算法可以用較小的計(jì)算量完成各子載波能量的快速分配，之后，仍需要根據(jù)誤碼率指標(biāo)為每個(gè)子載波分配離散比特?cái)?shù)。上述固定門限比特分配方法只是確保了每個(gè)子載波 k _i上的信噪比 SNRk_i不小于所要求信噪比門限SNRMk_i，而該子載波上超過門限的剩余能量部分既不能用于提高傳輸速率（因?yàn)?SNRk_i＜SNRMk_i*2），也不能提高整個(gè)碼元的誤碼率性能（由其他子載波上的最低誤碼率決定［12］），能量被白白浪費(fèi)。為進(jìn)一步提高碼元速率，本文采用貪婪算法對子載波發(fā)射能量及比特?cái)?shù)進(jìn)行2次修正。算法步驟是：

步驟 1 由式（12），利用簡化注水算法計(jì)算分配給子載波 k _i的正值能量

步驟2 在誤碼率 Pe*要求下，利用固定門限比特分配算法為子載波 k_i選擇MQAM的調(diào)制模數(shù)按式（13），由 Mk_i的判決門限 SNRMk_i，計(jì)算子載波 k _i達(dá)到需要的最低發(fā)射功率，將作為分配給子載波 k _i的初始功率。

步驟 4 由式（15）計(jì)算每個(gè)子載波當(dāng)前的比特?cái)?shù)加 1（調(diào)制模數(shù) Mk_i加倍）后，需要付出的額外能量（仍滿足誤碼率要求）。

步驟 5 為了充分利用有限的碼元能量獲得最大的碼元速率，采用貪婪算法選擇增加1 bit信息所需額外能量 Δpk_i最小的子載波作為下一步的能量和比特分配對象。

相比于固定門限比特分配算法，貪婪比特分配算法能夠?qū)⒏鱾€(gè)子載波的剩余能量統(tǒng)籌利用，以最小的能量代價(jià)換取碼元通信速率的最大化；算法每次迭代即可獲得1 bit的碼元速率提升，消耗的計(jì)算量僅為3次簡單的數(shù)據(jù)更新運(yùn)算（步驟6）和1次比較運(yùn)算（步驟5，可采用排序算法），運(yùn)算效能較高。

5 仿真驗(yàn)證

OFDM 水聲通信信號帶寬設(shè)為 10～20 kHz，fsam=100 kHz， V=1024，通信使用的子載波總數(shù)N= 133。收、發(fā)端深度分別位于30 m和35 m，水深100 m，通信距離1 km，利用BELLHOP算法［13］計(jì)算水聲信道響應(yīng)h。

圖1給出了不同信噪比下，碼元速率最大化精確注水算法和本文簡化注水算法所得到的平均信道容量，可見，注水算法能夠根據(jù)當(dāng)前水聲信道特征更為合理地分配子載波能量，信道容量較子載波能量平均分配模式有顯著提高；同時(shí)，當(dāng)信噪比較低時(shí)，精確注水算法較簡化注水算法有更大的容量優(yōu)勢，但隨著信噪比的增加，當(dāng)信噪比高于0 dB時(shí)（點(diǎn)對點(diǎn)水聲通信通常滿足），本文簡化注水算法能以較小的運(yùn)算量獲得與精確注水算法相當(dāng)?shù)男诺廊萘浚m合于在實(shí)時(shí)性要求高，信道變化較快的水聲通信中使用。

在簡化注水算法子載波能量分配基礎(chǔ)上，圖 2給出了當(dāng)誤碼率指標(biāo)時(shí)，采用固定門限比特分配算法及本文貪婪比特分配算法，每個(gè)OFDM碼元能夠傳遞的比特?cái)?shù)，兩種算法各自的歸一化剩余碼元能量如圖3所示。綜合圖2，圖3可以發(fā)現(xiàn)，固定門限比特分配算法存在有較大的剩余能量，尤其是在信噪比較低、而誤碼率要求較高的惡劣條件時(shí)，碼元能夠發(fā)送的比特?cái)?shù)非常有限（如圖2，圖 3的“○”），原因在于注水算法是在連續(xù)速率流假設(shè)下進(jìn)行最優(yōu)能量分配，當(dāng)利用固定門限算法對子載波比特?cái)?shù)進(jìn)行離散分配時(shí)，注水算法在很多惡劣子載波上的等效信噪比小于1 bit（模2）傳輸?shù)呐袥Q門限 SNR2，因此無法傳遞信息，從而造成大量的能量浪費(fèi)。而本文貪婪比特分配算法能夠?qū)⑹Ｓ嗄芰拷y(tǒng)籌、充分利用，最大限度提高碼元的通信速率，即使在惡劣應(yīng)用條件，也會(huì)將剩余能量集中在某些信道條件相對較好的子載波上使用。圖 2，圖3仿真表明，本文貪婪比特分配算法能夠最大限度地利用碼元能量，在任何信噪比條件下，碼元剩余能量均較小，碼元傳輸比特?cái)?shù)較固定門限比特分配算法有顯著提高，實(shí)現(xiàn)了碼元傳輸速率的優(yōu)化提升。

上述仿真都是在通信發(fā)射端獲得接收端位置良好的水聲信道質(zhì)量估計(jì)前提下得出的結(jié)論，然而在實(shí)際水聲通信應(yīng)用中，接收端處的信道響應(yīng)并不恒定，如果接收端每次都將信道質(zhì)量估計(jì)結(jié)果返回給發(fā)送端，將會(huì)占用極大的通信帶寬。實(shí)際上，水聲信道具有典型的慢時(shí)、空變特性，在較長一段時(shí)間內(nèi)，接收端水聲信道響應(yīng)可以近似認(rèn)為是以當(dāng)前信道估計(jì)結(jié)果為均值的指數(shù)分布規(guī)律變化［14］。基于此，接收端一般采用信道量化聚類壓縮技術(shù)［15］，用較少比特?cái)?shù) m（代表 2m個(gè)聚類中心）向發(fā)射端反饋信道估計(jì)結(jié)果。假設(shè)接收端采用k均值量化聚類壓縮技術(shù)［15］，發(fā)送端根據(jù)收到的m bit反饋信息，以聚類中心作為此時(shí)的信道估計(jì)（≠），利用本文算法處理后的剩余能量如圖4所示（誤碼率= 10-4），可見，本文算法在認(rèn)定的信道質(zhì)量'下都能充分利用碼元能量，剩余能量值均較小，碼元傳遞比特?cái)?shù)與圖2的“*”曲線接近，實(shí)現(xiàn)了傳輸速率提升。但是由于利用分配的子載波能量和比特?cái)?shù)會(huì)造成誤碼率指標(biāo)下降。圖 5給出了時(shí)變水聲信道環(huán)境下，采用不同量化位數(shù)m，本文算法獲得的實(shí)際誤碼率，可見，本文算法只需要較少的量化位數(shù)m，即可達(dá)到與期望值 Pe*=10-4接近的誤碼率，同時(shí)，較小的m會(huì)大幅節(jié)省用于信道估計(jì)反饋的帶寬，變向提高整個(gè)OFDM水聲通信系統(tǒng)的吞吐率，適于時(shí)變水聲信道應(yīng)用。

6 結(jié)束語

圖1 不同能量分配算法的信道容量

圖2 兩種算法碼元傳遞比特?cái)?shù)比較

圖3 兩種算法剩余能量比較圖

圖4 不同比特信道壓縮聚類后，新算法剩余能量曲線

圖5 不同比特信道壓縮聚類后，新算法誤碼率性能曲線

本文研究了一種實(shí)現(xiàn) OFDM 速率最大化的自適應(yīng)水聲通信算法，通過對注水算法的合理簡化，大幅減少了載波能量初始分配所需的運(yùn)算量，采用基于貪婪的子載波離散比特分配方法對固定門限方法的剩余能量進(jìn)行統(tǒng)籌、合理分配，最大限度地提高碼元傳遞比特?cái)?shù)。仿真證明，新算法能夠帶來通信速率的顯著提升，碼元能量得到了充分利用，同時(shí)，在有信道估計(jì)損失情況下，算法仍具有良好的誤碼性能，適合于在慢時(shí)、空變的水聲信道中應(yīng)用。

［1］ Kumar M and Tiwari S. Performance evaluation of conventional and wavelet based OFDM system［J］. Journal of Electronics and Communication， 2013， 67（4）: 348-354.

［2］ Abdallaa A， Ferreiraa R， and Shahparia A. Improved nonlinear tolerance in ultra-dense WDM OFDM system［J］. Optics Communications， 2014， 325（1）: 88-93.

［3］ 趙越， 方旭明， 黃博. 基于負(fù)載均衡的OFDMA雙跳中繼網(wǎng)絡(luò)資源分配策略［J］. 西南交通大學(xué)學(xué)報(bào)， 2013， 48（1）: 94-100. Zhao Yue， Fang Xu-ming， and Huang Bo. Resource allocation scheme based on load balancing for OFDMA two-hop relay networks ［J］. Journal of Southwest Jiaotong University， 2013，48（1）: 94-100.

［4］ Huang Xiao-peng and Lawrence V B. Capacity criterionbased bit and power loading for shallow water acoustic OFDM system with limited feedback［C］. 2011 IEEE 73rd Vehicular Technology Conference， Budapest， 2011: 372-377.

［5］ Xu X K， Wang Z H， Zhou S L， et al.. Parameterizing both path amplitude and delay variations of underwater acoustic channels for block decoding of orthogonal frequency division multiplexing［J］. Acoustical Society of America， 2012， 131（6） : 4672-4679.

［6］ 羅亞松， 胡洪寧， 劉忠， 等. 基于OFDM信號的自適應(yīng)水聲信道估計(jì)算法［J］. 電子科技大學(xué)學(xué)報(bào)， 2014， 43（5）: 678-684. Luo Ya-song， Hu Hong-ning， Liu Zhong， et al.. Self-adjusting underwater acoustic channel estimation algorithm based on OFDM signals［J］. Journal of University of Electronic Science and Technology of China， 2014， 43（5）: 678-684.

［7］ Hai H N， Soh W S， and Motani M. A Bidirectional-Concurrent MAC protocol with packet bursting for underwater acoustic networks［J］. IEEE Journal of Oceanic Engineering，2013， 38（3）: 547-565.

［8］ Liao W H and Huang C C. SF-MAC: a spatially fair MAC protocol for underwater acoustic sensor networks［J］. IEEE Sensors Journal， 2012， 12（6）: 1686-1694.

［9］ Lei Yang and Long Zhou. Adaptive bit loading algorithm for OFDM underwater acoustic communication system［C］. 2011 International Conference on Electronics and Optoelectronics （CEOE 2011）， Dalian， China， 2011: 350-352.

［10］ 王巍， 喬鋼， 王玥， 等. 多輸入多輸出正交頻分復(fù)用淺海水聲通信打孔判決反饋信道估計(jì)算法［J］. 兵工學(xué)報(bào)， 2013， 34（9）: 1116-1124. Wang Wei， Qiao Gang， Wang Yue， et al.. Decision feedback estimation of multiple input/multiple output orthogonal frequency division multiplexing channel based on punching technique via UWA shallow sea［J］. Acta Armamentarii， 2013，34（9）: 1116-1124.

［11］ Xiao Dong， Mo Fu-yuan， and Chen Geng. Study of low bit rate speech codec algorithm in underwater acoustic communication［J］. Chinese Journal of Acoustics， 2013， 32（4）: 411-423.

［12］ 羅亞松， 許江湖 劉忠， 等. 基于判決反饋均衡的水聲信道估計(jì)與通信算法［J］. 上海交通大學(xué)學(xué)報(bào)， 2014， 48（5）: 685-692. Luo Ya-song， Xu Jiang-hu， Liu Zhong， et al.. Self-adjusting underwater acoustic channel estimation and communication algorithms based on decision feedback equalizer［J］. Journal of Shanghai Jiaotong University， 2014， 48（5）: 685-692.

［13］ 張旭， 程琛， 劉艷. 西北太平洋副熱帶模態(tài)水形成區(qū)聲傳播特性分析［J］. 海洋學(xué)報(bào)， 2014， 36（9）: 94-101. Zhang Xu， Cheng Chen， and Liu Yan. Acoustic propagation effect caused by subtropical mode water of northwestern pacific［J］. Acta Oceanologica Sinica， 2014， 36（9）: 94-101.

［14］ Huang Xiao-peng. Capacity criterion-based power loading for underwater acoustic OFDM system with limited feedback［C］. Wireless Communications， Networking and Information Security， Beijing， 2010: 54-58.

［15］ 謝娟英， 王艷娥. 最小方差優(yōu)化初始聚類中心的 K-means算法［J］. 計(jì)算機(jī)工程， 2014， 40（8）: 205-211. Xie Juan-ying and Wang Yan-e. K-means algorithm based on minimum deviation initialized clustering centers［J］. Computer Engineering， 2014， 40（8）: 205-211.

羅亞松： 男，1982年生，講師，主要研究方向?yàn)樗伦越M織網(wǎng)絡(luò).

許江湖： 男，1975年生，講師，主要研究方向?yàn)樗暪こ?

胡洪寧： 男，1980年生，講師，主要研究方向?yàn)橄到y(tǒng)工程.

賀靜波： 男，1978年生，講師，主要研究方向?yàn)樗曤娮訉?

陳占偉： 男，1982年生，工程師，主要研究方向?yàn)闈撏е笓].

Research on Self-adjusting OFDM Underwater Acoustic Communication Algorithm for Transmission Rate Maximization

Luo Ya-song①Xu Jiang-hu①Hu Hong-ning①He Jing-bo①Chen Zhan-wei②
①（Electronics Engineering College， Naval Universit of Engineering， Wuhan 430033， China）②（Coach Division of Logistic Department， Navy Qingdao Marine Area， Qingdao 266102， China）

Based on the modeling of transmission rate maximization problem for Orthogonal Frequency Division Multiplexing （OFDM） underwater acoustic communication， the computing amount for water-filling algorithm is reduced reasonably， which realizes rapid assignments for sub-carriers’ initial energy. To deal with the serious waste of symbol energy existing in the bit assigning algorithm based on fixed threshold， the new greedy algorithm is brought forward， which improves the OFDM transmission rate remarkably with less cost of additional calculation. Results show that the new algorithm behaves favorable properties of transmission rate and Bit Error Rate （BER）even with incomplete channel estimation information suitable for time-variant underwater acoustic channel.

Underwater acoustic communication； Orthogonal Frequency Division Multiplexing （OFDM）； Selfadjusting； Water-filling algorithm

The National Natural Science Foundation of China （61101205）

TN929.3

A

1009-5896（2015）12-2872-05

10.11999/JEIT150440

2015-04-15；改回日期：2015-08-03；網(wǎng)絡(luò)出版：2015-10-13

*通信作者：羅亞松 yours_baggio@sina.com

國家自然科學(xué)基金（61101205）