胡耀輝 韓樹平 劉建波 張 奇 張延厚

①(海軍潛艇學(xué)院 青島 266199)

②(中國人民解放軍92020部隊 青島 266071)

③(中國人民解放軍91323部隊 濟(jì)源 454650)

1 引言

直接序列擴(kuò)頻(Direct Sequence Spread Spectrum, DSSS)通信技術(shù)通過擴(kuò)展帶寬降低信號譜密度，兼具低檢測概率和抗干擾特性，廣泛應(yīng)用于穩(wěn)健水聲通信[1,2]。然而，在獲得較高擴(kuò)頻增益的同時，其對水聲信道多普勒效應(yīng)產(chǎn)生的相位抖動較為敏感且通信速率較低，面對AUV等水下移動平臺復(fù)雜運動狀態(tài)下的可靠通信時存在一定挑戰(zhàn)。提升水下擴(kuò)頻通信穩(wěn)健性和頻譜利用率已成為近年來水聲擴(kuò)頻通信領(lǐng)域的兩個重要研究方向。

在抗載波相位跳變實現(xiàn)水下穩(wěn)健擴(kuò)頻通信方面，傳統(tǒng)做法是采用鎖相環(huán)技術(shù)在完成捕獲鎖定后動態(tài)跟蹤載波頻率相位的變化[3]，然而其在低信噪比下性能受到限制。近年來，穩(wěn)健擴(kuò)頻通信研究多集中于時域差分?jǐn)U頻方向，且主要針對海洋復(fù)雜動力環(huán)境和弱時變多普勒運動擴(kuò)頻，針對水下AUV復(fù)雜運動狀態(tài)下的時變多普勒擴(kuò)頻通信研究還較少。文獻(xiàn)[4]使用大量數(shù)據(jù)揭示了海洋動力環(huán)境對水聲通信的影響；殷敬偉等人[5]采用單矢量差分能量檢測器解決了矢量水聽器轉(zhuǎn)動漂移引起的相位跳變，采用能量法獲得了較穩(wěn)健矢量處理增益；杜鵬宇等人[6]提出一種勻速背景下的擴(kuò)頻通信實時多普勒估計方法，并使用一種適用于移動擴(kuò)頻通信的改進(jìn)差分能量檢測器，實現(xiàn)單擴(kuò)頻碼寬度內(nèi)相位穩(wěn)定下的穩(wěn)定解碼，然而并未考慮復(fù)雜相對運動多普勒時變情況下的穩(wěn)健通信，且采用重復(fù)擴(kuò)頻碼進(jìn)行差分調(diào)制，雖然可實現(xiàn)低信噪比水聲穩(wěn)健通信，但其時域循環(huán)特征明顯，頻譜利用率偏低。

在提升擴(kuò)頻通信頻譜利用率方面，多種擴(kuò)頻調(diào)制技術(shù)先后被提出，主要包括多進(jìn)制擴(kuò)頻通信[7,8]、并行組合擴(kuò)頻通信[9]及近年來新興的碼索引調(diào)制通信技術(shù)[10,11]等。其中多進(jìn)制擴(kuò)頻通信將多個信息比特映射為M個擴(kuò)頻通信碼組中的1個進(jìn)行通信，又稱M元擴(kuò)頻通信，它為之后的擴(kuò)頻調(diào)制技術(shù)發(fā)展奠定了基礎(chǔ)；并行組合擴(kuò)頻通信是在多進(jìn)制擴(kuò)頻通信基礎(chǔ)上發(fā)展而來的，在擴(kuò)頻碼組中選取多個進(jìn)行疊加，可獲得更高頻帶利用率；新興的碼索引調(diào)制技術(shù)是將傳輸信息比特分為映射組和基帶調(diào)制組，并分別映射為擴(kuò)頻碼索引和MPSK調(diào)制極性符號，調(diào)制符號的實部和虛部再分別選擇各自激活的擴(kuò)頻碼進(jìn)行擴(kuò)頻。在相同頻譜效率時，碼索引調(diào)制相對于多進(jìn)制擴(kuò)頻明顯減少擴(kuò)頻碼的使用個數(shù)，同時較并行組合擴(kuò)頻降低了信噪比要求。另外，碼索引調(diào)制相對于傳統(tǒng)DSSS在每個傳輸時隙增加了映射比特，因而提升了系統(tǒng)的頻譜利用率和能量效率[12]。

由于多普勒效應(yīng)廣泛存在于海上水聲通信過程，針對削弱擴(kuò)頻增益、帶來劇烈相位跳變的大多普勒情況，一般需要開展多普勒估計與補償。目前水聲通信多普勒估計方法主要有：采用對多普勒敏感的訓(xùn)練序列開展的時頻2維搜索法、采用對多普勒效應(yīng)容忍度高的調(diào)頻信號開展的塊估計法、采用與編碼技術(shù)結(jié)合的迭代估計算法等。其中，時頻2維搜索在追求模糊度函數(shù)最大化的過程計算量大，不易應(yīng)用于實時通信過程；編碼迭代技術(shù)要求信噪比較高，不適用低信噪比的水聲通信過程；而基于調(diào)頻信號的塊估計法，在多普勒時變較弱的水聲多普勒粗估計中應(yīng)用較多，難以適用于復(fù)雜運動狀態(tài)下的水聲通信。目前，隨著AUV裝備技術(shù)的發(fā)展其機動性能得到不斷提升[13,14]，高機動性AUV移動速度最高可達(dá)8～15 kn、加速度可達(dá)0.5～1 m/s2、加速度變化率則可達(dá)0.05～0.15 m/s3，其復(fù)雜運動背景下的交互穩(wěn)健通信需要實現(xiàn)多普勒估計補償?shù)耐瑫r提高通信系統(tǒng)對多普勒的魯棒性。

本文提出一種基于2維組合擴(kuò)頻碼組的幀結(jié)構(gòu)，傳輸信息比特分為映射調(diào)制組和極性調(diào)制組，其中映射調(diào)制組根據(jù)擴(kuò)頻碼組中重復(fù)碼、與重復(fù)碼并行的排列碼、2維組合碼各自設(shè)置的映射索引表激活對應(yīng)擴(kuò)頻碼，極性調(diào)制組則對應(yīng)調(diào)制2維并行組合擴(kuò)頻碼中的其中1維擴(kuò)頻碼。幀結(jié)構(gòu)中重復(fù)碼用于碼組間的多普勒動態(tài)估計，由于其并非采用固定的擴(kuò)頻序列，降低了時域循環(huán)特征。該結(jié)構(gòu)擴(kuò)頻碼組可用于估計多普勒時變的收發(fā)相對運動過程，在接收端多普勒處理過程采用頻率壓縮-能量接收器(Frequency Compression-Energy Detector, FCED)提高多普勒處理容限，應(yīng)對多普勒估計誤差。最后，提出一種組合差分解碼算法，實現(xiàn)2維組合擴(kuò)頻間的極性差分解碼。經(jīng)理論及仿真試驗分析，該水聲通信方法具備一定的抗復(fù)雜時變多普勒能力。

下文中N為擴(kuò)頻碼索引庫的擴(kuò)頻碼數(shù)量，n為幀結(jié)構(gòu)中擴(kuò)頻碼組的碼數(shù)量，fc為 載頻，fs為采樣率，B為通信帶寬，L為擴(kuò)頻碼長度，β為多普勒因子，c為聲速(仿真計算中取1500 m/s)，a為相對運動加速度。

2 組合差分?jǐn)U頻水聲通信

2.1 碼元幀結(jié)構(gòu)

如圖1所示，擴(kuò)頻碼組前兩個并行擴(kuò)頻碼中疊加有相同重復(fù)擴(kuò)頻碼PNs，碼組間以時間順序分為前端重復(fù)碼和后端重復(fù)碼，用于時域相關(guān)并開展多普勒粗估計。需要注意的是，與PNs疊加的兩個擴(kuò)頻碼是存在順序關(guān)系的排列，其映射對應(yīng)的信息比特數(shù)不同于后面的疊加組合，需單獨設(shè)立索引表進(jìn)行映射，這里標(biāo)注為PNp。后面8對并行組合擴(kuò)頻碼均標(biāo)注為PNc。激活擴(kuò)頻碼根據(jù)圖2的信息比特-擴(kuò)頻碼索引關(guān)系進(jìn)行選擇、極化和疊加。

圖1 2維組合差分?jǐn)U頻碼元幀結(jié)構(gòu)及多普勒塊估計示意

圖2 發(fā)射端擴(kuò)頻碼組索引映射關(guān)系

該方法每組并行擴(kuò)頻碼調(diào)制信息比特關(guān)系見表1。以本文通信參數(shù)背景為例，同樣擴(kuò)頻碼組長度下該方法較傳統(tǒng)差分?jǐn)U頻傳輸比特信息量(10 bit)提高8.7倍。

表1 每組碼元調(diào)制信息比特數(shù)

2.2 時變多普勒估計方法

對圖1中擴(kuò)頻碼組中重復(fù)碼進(jìn)行時域相關(guān)處理，得到碼組兩端時刻的收發(fā)端徑向相對速度V1,V2，為下一步的M元能量接收處理及頻率壓縮提供基礎(chǔ)。下面從理論上分析其可行性。這里假設(shè)單個擴(kuò)頻碼時間寬度內(nèi)多普勒因子近似視為不變量。設(shè)相鄰重復(fù)碼時域波形為PNs(t)，相鄰擴(kuò)頻碼多普勒因子為β1和β2，h1和h2分別表示相鄰擴(kuò)頻碼時間寬度內(nèi)的信道沖擊響應(yīng)函數(shù)，則其對應(yīng)接收信號經(jīng)下變頻至基帶處理后，s1和s2可分別表示為

由于相鄰擴(kuò)頻碼多普勒因子變化量為相對小量，式(5)中( 1+β1)/(1+β2)近似取1。由擴(kuò)頻碼自相關(guān)性質(zhì)可知，當(dāng)式(6)成立時，R12產(chǎn)生相關(guān)峰值。這里也可以理解為，當(dāng)擴(kuò)頻碼寬度內(nèi)相對速度變化量aL/B小于擴(kuò)頻增益的相對速度容限時，相鄰重復(fù)碼可產(chǎn)生相關(guān)峰。以本文提出的特例情況分析，擴(kuò)頻增益在能量法計算的相對速度容限取0.4 m/s(分析見下節(jié))，加速度小于4.7 m/s2均可滿足，實際海上運動平臺目前均滿足此要求。這里也可看出，帶寬越寬、碼長越短越易滿足此條件。

這里幀長度不宜過長，否則當(dāng)幀后端重復(fù)碼最大可能縮放量超過1個擴(kuò)頻碼長時，在重復(fù)碼與相鄰組合碼相同的情況下無法正確估計多普勒因子。這里假設(shè)通過重復(fù)碼相關(guān)獲得的圖1中碼組兩端相對速度分別為V1, V2。

為了考察以上算法的效果，這里結(jié)合本文特例參數(shù)對重復(fù)碼相關(guān)求相對速度結(jié)果進(jìn)行仿真模擬，0 dB信噪比背景下得出不同相對速度和加速度背景下相對速度估計絕對誤差如圖4所示，重復(fù)碼相關(guān)增益變化如圖5所示。

由圖4可以看出，重復(fù)碼估計誤差隨相對速度增加有波動增加趨勢，這主要是數(shù)字采樣誤差帶來的；同時，加速度超過1.5 m/s2時，因重復(fù)碼前后多普勒因子變化較大，相對速度絕對估計誤差猛然增大。由圖5可看出，解擴(kuò)增益隨相對速度和加速度的增加呈減弱趨勢，且在加速度1.5 m/s2處存在相對低值區(qū)?？偟膩砜矗嚓P(guān)增益變化較小，圖中增益極小值較極大值僅下降3 %，對解擴(kuò)影響不大。

圖3 重復(fù)碼相關(guān)時域窗截取過程示意

圖4 重復(fù)碼相對速度估計誤差

圖5 重復(fù)碼相關(guān)增益隨相對速度及加速度變化情況

圖6、圖7為將采樣率提高1倍后的相對速度估計誤差圖與相關(guān)增益圖?？梢钥闯?，采樣率加倍后相對速度估計誤差隨相對速度增加波動趨勢熨平趨緩，較原采樣率仿真結(jié)果相對速度誤差絕對值整體有所減小。然而，在實際工程應(yīng)用中，受硬件性能的限制，無法無限制地增加采樣率，且低信噪比下相關(guān)峰值定位精度難以保證高準(zhǔn)確度，這樣重復(fù)碼估計相對速度必然帶入一定誤差。

圖6 采樣率加倍后重復(fù)碼相對速度估計誤差

圖7 采樣率加倍后重復(fù)碼相關(guān)增益

為了適應(yīng)低信噪比環(huán)境，對于每幀信號寬度內(nèi)收發(fā)相關(guān)運動視為恒加速度運動，對每幀6組重復(fù)碼估計出的相對速度進(jìn)行最小二乘擬合處理，具體分析見第3節(jié)。對碼組內(nèi)多普勒估計時，認(rèn)為單個擴(kuò)頻碼寬度內(nèi)是恒多普勒勻速運動。多普勒估計后，需進(jìn)行擴(kuò)頻碼段截取并分別進(jìn)行下變頻處理，按照估計相對速度截取碼位的同時，考慮到估計誤差及其隨時域累積問題，這里按照2 m/s的極限誤差并采用累積計算，對下變頻及對應(yīng)解擴(kuò)處理的時域截取做擴(kuò)窗處理，本文特例擴(kuò)窗長度取0.07 L。

2.3 頻率壓縮-能量接收器

為了應(yīng)對上節(jié)中復(fù)雜相對運動引起的多普勒估計誤差，本文接收端采用頻率壓縮-能量接收器(FCED)，增加系統(tǒng)的多普勒容限，提高對多普勒估計誤差的魯棒性。

2.3.1 能量檢測接收器

擴(kuò)頻碼在階數(shù)較低時對時域壓縮擴(kuò)展有一定的容限，在 Δv較小時不明顯影響其擴(kuò)頻增益[6]。分別使用碼索引庫N個擴(kuò)頻碼對式(11)進(jìn)行解擴(kuò)處理，得到相關(guān)峰值最大值對應(yīng)的擴(kuò)頻碼，進(jìn)而完成索引解碼。這里采用虛擬時反聚焦各聲線能量，進(jìn)一步提高解擴(kuò)增益。

由式(10)對相位跳變的近似處理可以看出，對于能量檢測接收器，載波頻率越低、估計相對運動速度誤差越小、擴(kuò)頻碼寬度越小，越有利于減小單個擴(kuò)頻碼寬度內(nèi)的載波跳變，保持?jǐn)U頻增益。這里以本文特例參數(shù)進(jìn)行仿真分析，擴(kuò)頻信號采用能量檢測接收器后擴(kuò)頻增益對多普勒頻偏的容限情況如圖8所示?？梢钥闯?，隨著相對速度估計誤差的增大，擴(kuò)頻碼內(nèi)載波跳變不再穩(wěn)定，使式(10)的近似處理難以成立，造成相關(guān)增益也明顯降低。然而，能量接收器較傳統(tǒng)下變頻處理較明顯地提升了多普勒估計誤差容限。設(shè)以標(biāo)準(zhǔn)化相關(guān)系數(shù)0.78為閾值，則擴(kuò)頻索引解碼多普勒容限較傳統(tǒng)下變頻處理獲得3 dB左右增益。

另外，圖8僅反映解擴(kuò)增益絕對值的影響，其極性受載波相位跳變影響也十分明顯。圖9為相對速度誤差0.3 m/s時對應(yīng)的傳統(tǒng)直接擴(kuò)頻信號解擴(kuò)效果和能量接收器解擴(kuò)效果(極性均為1)，可以看出，能量法在索引解碼時對于0.3 m/s的多普勒估計誤差有較好的容限；傳統(tǒng)處理結(jié)果含載波跳變項，擴(kuò)頻結(jié)果周期性變化甚至出現(xiàn)倒現(xiàn)象，難以準(zhǔn)確實現(xiàn)解擴(kuò)解碼。

圖8 不同接收處理方法擴(kuò)頻增益隨相對速度誤差變化過程

圖9 不同下變頻處理方法擴(kuò)頻增益效果

2.3.2 頻率壓縮算法

頻率壓縮捕獲(Frequency Compression Search, FCS)[15,16]是近年來興起的一種擴(kuò)頻碼動態(tài)時頻2維壓縮捕獲中的預(yù)處理方法，該方法將本地不同頻點的載波信號進(jìn)行非相干疊加，疊加后的信號作為本地載波與輸入信號進(jìn)行下變頻處理，然后進(jìn)行相干捕獲。本文所提出的頻率壓縮FC用于多普勒估計誤差背景下的快速碼索引庫激活碼搜索，具體操作方法是：

步驟1 在接收端設(shè)置多個NCO生成多個頻點本振信號，頻點設(shè)置以擴(kuò)頻碼寬度內(nèi)多普勒估計補償后載頻fc+fd為中心點，以多普勒容限df為步長間隔取多個頻點，這里fd為采用2.2節(jié)方法得出的多普勒頻偏粗估計值。

步驟2 對多個頻點本地載波進(jìn)行非相干疊加處理，而后對數(shù)字信號進(jìn)行下變頻處理和解擴(kuò)處理，選取擴(kuò)頻增益最大的前兩個碼即解碼索引碼。

步驟3 對數(shù)字接收信號利用各頻點本地載波分別進(jìn)行2次下變頻處理，利用步驟2得到的解碼索引碼各頻點下變頻結(jié)果進(jìn)行解擴(kuò)，通過解擴(kuò)增益最大值找到對應(yīng)的頻點即帶有df/2誤差的載波估計值，該估計值將用于組合差分解碼。

本文仿真試驗中接收端采用5個NCO進(jìn)行接收端處理，其生成載波頻率分別為fc+fd,fc+fd±df及fc+fd±2df。根據(jù)圖8模擬結(jié)果，能量接收器對多普勒容限定為±0.4 m/s，對應(yīng)的df為9 Hz。圖10為經(jīng)頻率壓縮算法后得出的擴(kuò)頻增益隨相對速度估計誤差變化過程?？梢钥闯?，經(jīng)頻率壓縮算法處理，系統(tǒng)對多普勒估計誤差的容限增大到[-2.0 m/s 2.0 m/s]，較單純的能量接收器擴(kuò)展了5倍。該值在近海能較好地覆蓋2.2節(jié)大時變多普勒估計算法的估計誤差和中低海況下一般海洋水文動力環(huán)境的波動影響。

圖10 采用頻率壓縮處理前后能量接收器擴(kuò)頻增益隨相對速度誤差變化過程

2.4 組合差分算法

這樣可以判斷dB的極性值，從而完成組合差分的極性解調(diào)解碼。在實際操作中，由于組合差分前FC-ED已選定了最大化解擴(kuò)輸出的NCO頻點，因此引起載波相位?的變化對應(yīng)的最大頻率偏差為df/2，對應(yīng)本文中相對速度估誤差為0.4 m/s。

圖11為綜合重復(fù)碼相關(guān)多普勒粗估計、FC-ED算法、組合差分算法的接收端數(shù)據(jù)處理流程。接收數(shù)據(jù)經(jīng)FC-ED算法后進(jìn)行映射解碼、經(jīng)組合差分算法后得出極性解碼，最后完成接收端譯碼。

圖11 組合差分?jǐn)U頻通信接收端數(shù)據(jù)處理流程

3 仿真研究

首先考察本文提出的FC-ED算法在靜態(tài)和弱相對運動背景下的處理性能進(jìn)行仿真分析。這里結(jié)合傳統(tǒng)并行組合擴(kuò)頻進(jìn)行性能對比，以下處理未引入多普勒估計。圖12給出了收發(fā)靜態(tài)和以收發(fā)端0.3 m/s相對運動背景下兩種方法的誤碼率情況。這里傳統(tǒng)并行組合差分采用擴(kuò)頻碼庫數(shù)量為9，并行疊加發(fā)送3個擴(kuò)頻碼，對應(yīng)通信速度與本文特例一致，即頻譜利用率相同?？梢钥闯觯⑿薪M合擴(kuò)頻通信在相對運動背景下已難以有效解碼，而本文提出的組合差分接收器受多普勒影響較小。靜態(tài)背景下，本文提出的組合差分接收器在相同信噪比條件下較并行組合擴(kuò)頻模式誤碼率存在優(yōu)勢，這主要是2維擴(kuò)頻降低疊加干擾，配合較大擴(kuò)頻碼索引庫，犧牲一定計算力換取低信噪比性能；比較組合差分接收器在靜態(tài)與相對運動背景下的解碼效果，相對運動背景下雖然采用了能量接收，但擴(kuò)頻增益在多普勒影響下也會受到削弱影響，在極低信噪比下解碼效果不及靜態(tài)通信效果。

圖12 不同多普勒背景下兩種通信方法解碼性能對比

考慮到水下移動平臺在高速運動狀態(tài)下噪聲背景復(fù)雜，為考察本文提出多普勒估計方法在低信噪比、復(fù)雜運動姿態(tài)背景下估計效果，這里對低信噪比下大多普勒時變和勻變加速度相對運動過程進(jìn)行仿真分析。

以-8 dB信噪比進(jìn)行勻加速背景下的重復(fù)碼相對速度估計仿真，得到相對速度誤差與相關(guān)增益變化結(jié)果如圖13、圖14所示。圖中相對速度指重復(fù)碼中前碼初始時刻收發(fā)端相對速度，加速度指一對重復(fù)碼寬度內(nèi)勻加速運動對應(yīng)的相對運動加速度。可以看出，低信噪比條件下相對速度估計誤差波動劇烈，規(guī)律性較弱，但總體絕對誤差范圍在1.2 m/s以下，在M元能量接收器與頻率壓縮算法的容限以內(nèi)；相關(guān)增益變化劇烈，但最小值較最大值僅下降25 %左右，在一定信噪比條件下仍可以實現(xiàn)有效估計。

圖13 -8 dB信噪比下重復(fù)碼相對速度估計誤差

圖14 -8 dB信噪比下重復(fù)碼相關(guān)增益

這里以4 m/s為初始速度，1 m/s2為恒加速度進(jìn)行大多普勒時變相對速度估計，一幀以內(nèi)的相對速度估計效果如圖15所示?？梢钥闯?，低信噪比下相對速度估計誤差擺動更為劇烈，幅度更大；在進(jìn)行最小二乘擬合后，高信噪比和低信噪比條件下的相對速度估計誤差均有所減小，且低信噪比相對速度誤差減小更為明顯，誤差峰值得以明顯削弱。

圖15 恒加速度相對速度估計過程

下面考察多普勒估計方法對加速度勻變化背景下的變加速場景適用性情況，這里取仿真初始速度2 m/s，初始加速度1 m/s2，加速度變化率0.2 m/s3，1幀以內(nèi)的相對速度估計效果如圖16所示?？梢钥闯?，水下運動平臺推動馬力變化帶來的速度非線性時變較弱，低信噪比下雖波動較強，但取最小二乘擬合后，誤差明顯回落。

圖16 變加速度相對速度估計過程

針對組合差分算法極性解碼過程中不同信噪比下算法的表現(xiàn)進(jìn)行仿真分析，這里任取一對已知組合疊加碼，對其中一擴(kuò)頻碼進(jìn)行順序極性調(diào)制，進(jìn)行104bit極性碼通信仿真，獲得仿真誤碼率效果如圖17所示?？梢钥闯觯鄬λ俣日`差在0.4 m/s以下、信噪比高于-13 dB時誤碼率低于10-3，驗證了2.4節(jié)的分析。

圖17 組合差分算法隨相對速度誤差誤碼率仿真

4 試驗分析

以2019年11月青島國家深海試驗基地擴(kuò)頻通信海試數(shù)據(jù)為基礎(chǔ)進(jìn)行試驗分析。試驗水深8.9 m，試驗當(dāng)日海況3級，風(fēng)浪涌浪較大，聲速剖面為微弱負(fù)梯度，聲速從海面至海底由1510.5 m/s減小至1509.5 m/s。海試過程中的TCD幀結(jié)構(gòu)如圖1所示，共發(fā)送38幀數(shù)據(jù)，海試過程收發(fā)船平均距離3.02 km，接收船自由漂移未錨定。海試過程信道估計結(jié)果及通信信號時域過程分別如圖18、圖19所示，可以看出，通信過程信道多徑明顯，背景噪聲復(fù)雜，瞬時噪聲多發(fā)。

圖18 海試信道估計

圖19 信號時域過程

試驗過程接收信號的信噪比由式(14)估算

其中，Pr為通信信號持續(xù)時間內(nèi)接收信號平均功率，Pμ為通信間隙噪聲信號平均功率。38幀信號的接收信噪比如圖20所示。

圖20 海試通信信噪比過程

可以看出，海試過程背景噪聲功率變化強烈，同樣聲源功率背景下信噪比跨度可達(dá)10 dB以上。本次試驗接收信號信噪比主要集中在-4～1 dB，兩幀信號信噪比達(dá)-10 dB左右，整體信噪比較低。

擴(kuò)頻通信試驗數(shù)據(jù)分別采用鎖相環(huán)技術(shù)動態(tài)跟蹤方式和本文提出的能量檢測接收器-頻率壓縮擴(kuò)大多普勒容限方法進(jìn)行解擴(kuò)解碼，采用鎖相環(huán)技術(shù)解碼誤碼率為2.4 %，而采用本文提出的能量檢測接收-頻率壓縮提高多普勒容限技術(shù)解碼誤碼率為0.7 %。由于試驗過程為準(zhǔn)靜態(tài)水聲通信，環(huán)境引起的波動干擾在多普勒容限內(nèi)，未利用重復(fù)碼估計相對速度。這里鎖相環(huán)技術(shù)的誤差主要是低信噪比條件下部分低信噪比幀跟蹤過程失鎖、重復(fù)追蹤過程帶來的，而能量檢測接收-頻率壓縮方法誤碼主要來源于低信噪比條件下的瞬時噪聲影響了個別擴(kuò)頻碼的解擴(kuò)效果。

可以看出，本文提出的頻率壓縮-能量接收器(FC-ED)對近海復(fù)雜海況、收發(fā)端漂移擺動、較低信噪比條件有一定的適應(yīng)能力，然而在復(fù)雜背景噪聲處理方面仍需開展其他有效預(yù)處理工作。

為了考察本文擴(kuò)頻通信方法在近海多途及低信噪比背景下平臺相對運動及大多普勒時變條件下的擴(kuò)頻通信效果，下面通過對接收端信號采樣數(shù)據(jù)進(jìn)行插值，模擬大多普勒時變過程。這里需要指出的是，在實際相對運動過程中，水聲信道是緩慢變化的，通過文獻(xiàn)[6]的移動水聲試驗結(jié)論可知，在相鄰擴(kuò)頻碼的較短時間跨度內(nèi)，信道變化極小、高度相關(guān)，這里靜態(tài)采用圖18信道不會對試驗效果產(chǎn)生顯著影響。

對于勻加速度相對運動過程，這里對海試接收數(shù)據(jù)各幀分別進(jìn)行添加白噪聲處理(對已小于目標(biāo)信噪比的幀不做處理)，對接收數(shù)據(jù)按照不同假想相對速度進(jìn)行插值處理，取初始速度2 m/s，加速度范圍取0～2 m/s2，信噪比取-15～0 dB，進(jìn)行1000次蒙特卡羅仿真，得到平均誤碼率結(jié)果如圖21(a)所示?？梢钥闯?，本文提出的通信方法在信噪比-8 dB以上時誤碼率低于10-2，且受相對加速度變化影響較小。

圖21 組合差分?jǐn)U頻水聲通信方法海試推算解碼效果

針對勻變加速相對運動過程，對海試數(shù)據(jù)處理類似勻加速相對運動過程，這里取初始速度2 m/s，初始加速度0.5 m/s2，加速度變化率取0～0.3 m/s3，信噪比取-15～0 dB，進(jìn)行1000次蒙特卡羅仿真，得到平均誤碼率結(jié)果如圖21(b)所示?？梢钥闯?，在信噪比高于-7 dB時，誤碼率低于10-2，且受加速度變化率影響不明顯。

這里可以看出，本文提出的通信方法在較低信噪比下仍可以實現(xiàn)較高解碼準(zhǔn)確率，能夠適應(yīng)水下平臺復(fù)雜運動背景。另外，由于試驗采用的擴(kuò)頻碼碼長較短，相關(guān)增益有限，且有部分增益需抵消低信噪比下變速運動中相關(guān)增益的下降，該通信方法在極低信噪比下解碼效果迅速下降?？紤]到長擴(kuò)頻碼雖有較高擴(kuò)頻增益卻對多普勒效應(yīng)極為敏感，在水下移動通信系統(tǒng)設(shè)計時應(yīng)綜合考量系統(tǒng)通信噪聲環(huán)境、隱蔽性需求、水下移動平臺產(chǎn)生的相對運動速度范圍、多普勒時變可能范圍等因素，合理設(shè)置擴(kuò)頻碼索引庫內(nèi)擴(kuò)頻碼長度。

5 結(jié)束語

本文針對水下移動平臺時變多普勒水聲擴(kuò)頻通信提出一種2維組合差分?jǐn)U頻通信方法，采用設(shè)計TCD幀結(jié)構(gòu)中重復(fù)碼進(jìn)行瞬時相對速度粗估計；利用接收端頻率壓縮-能量接收器(FC-ED)提高系統(tǒng)對相對速度估計誤差的容限，并應(yīng)對部分環(huán)境波動引起的多普勒現(xiàn)象；提出一種組合差分算法，在能量檢測索引解碼基礎(chǔ)上，加入極性調(diào)制解碼，進(jìn)一步提高通信速率。仿真和海試推算表明，該通信方法在較低信噪比下具備較好抗多普勒時變干擾能力。在后續(xù)研究中，將結(jié)合海洋噪聲預(yù)處理技術(shù)和糾錯碼技術(shù)，進(jìn)一步降低指令傳輸中的誤碼率。