黃金鑫,周志權(quán),曹逸飛,趙 揚,2*

(1.哈爾濱工業(yè)大學(威海) 信息科學與工程學院, 威海 264209,中國;2.威海市智能光聲檢測與傳感技術(shù)重點實驗室, 威海 264209,中國)

0 引 言

天空地海潛一體化通信在軍事力量儲備和國民經(jīng)濟生產(chǎn)中的地位日益突出,其中空中平臺和水下目標之間的通信技術(shù)是構(gòu)建天空地海潛一體化通信的重要組成部分。在現(xiàn)有的空海通信方案中,潛行器接收飛機的信息一般分為上浮延展天線和飛機投放中繼浮標兩種方式[1-2],但這兩種方式在實際應(yīng)用中存在暴露目標的可能性,從而降低了潛行器的隱蔽性。

激光致聲空-水跨介質(zhì)通信通過熱膨脹或光擊穿機制在水氣交界面進行光聲轉(zhuǎn)換,可以將激光在空氣信道和聲波在海洋信道的傳輸優(yōu)勢結(jié)合起來,擺脫空-水界面的限制[3-13]。該通信方式具有無需在介質(zhì)中部署任何物理換能器及中繼器的特點,增強了通信目標的隱蔽性,是未來空中平臺與水下目標之間通信的重要方式之一。

自1962年前蘇聯(lián)PROKHOROV和美國WHITE的研究團隊先后發(fā)現(xiàn)了濃縮介質(zhì)在脈沖激光的作用下產(chǎn)生聲波的現(xiàn)象以來,國內(nèi)外相關(guān)機構(gòu)逐步開展光聲轉(zhuǎn)換技術(shù)在水下目標探測、空-水跨介質(zhì)通信等領(lǐng)域的研究工作[14-15]。BLACKMON等人[16-18]從理論上對激光致聲的線性特性和非線性特性進行了分析與研究,并從實驗角度驗證了該技術(shù)在通信領(lǐng)域的應(yīng)用可行性。研究發(fā)現(xiàn):激光致聲的非線性特征受環(huán)境干擾較大,實際情況下閾值難以控制,不利于保持通信過程的可靠性。因此,目前研究多利用激光致聲的線性特性,即根據(jù)激光致聲熱膨脹效應(yīng)來進行空中-水下跨介質(zhì)通信,例如ANTONELLI等人基于該效應(yīng)根據(jù)所需的聲波頻段,調(diào)制長脈沖激光并實現(xiàn)了通信功能[19]。此外,PENG等人[20]基于該效應(yīng)提出了一種控制激光器重復(fù)頻率進而實現(xiàn)空中-水下跨介質(zhì)通信的方法。

本文作者在上述研究的基礎(chǔ)上,依次針對激光致聲空中-水下跨介質(zhì)通信速率的需求,提出了一種多進制多頻移鍵控(n-multifrequency shift keying,n-MFSK)調(diào)制方式的激光致聲空-水跨介質(zhì)通信方法,分別在采用長脈沖法和重復(fù)頻率法的調(diào)制方式基礎(chǔ)上,結(jié)合現(xiàn)有激光器的關(guān)鍵技術(shù)參數(shù),通過仿真計算對通信速率的提升效果進行了評價與分析,為激光致聲在空中-水下跨介質(zhì)通信領(lǐng)域中的應(yīng)用提供了一定的參考依據(jù)。

1 熱膨脹機制

液體中的激光致聲,即在液體中激光能量達到一定強度激發(fā)產(chǎn)生聲波的過程。根據(jù)激光相互作用區(qū)域的能量密度與時空分布,可把液體中光聲效應(yīng)機制主要分為熱膨脹、汽化與光擊穿3種。與汽化機制和光擊穿機制相比,熱膨脹機制物理過程簡單,實驗上容易實現(xiàn)。相對于另外兩種致聲機制,熱膨脹機制致聲原理的研究相對比較成熟,本文中關(guān)于液體中激光致聲空-水跨介質(zhì)研究將針對熱膨脹機制進行。

在熱膨脹機制下,假設(shè)液體為均勻的理想流體,激光激發(fā)水下聲波的轉(zhuǎn)換方程[21]:

(1)

其中:

s=(αc2βE0/cp)·exp(-αz)

(2)

(3)

式中:s為聲壓;c為聲速;x、y、z為空間坐標方向距離;t為時間;α為水的吸收系數(shù);E0為表面激光能量密度;cp為單位質(zhì)量定壓比熱容;H(x,y,z,t)為單位時間內(nèi)液體吸收并進行能量轉(zhuǎn)化的電磁能量密度;β為液體的體積熱膨脹系數(shù)(ρ0為密度,T為溫度,假設(shè)T在過程中不變,β為常數(shù))。由參考文獻[16],將式(1)傅里葉變換取實部,得到其激光垂直入射表達式為:

(4)

其中,

(5)

式中:r為接收點與聲源之間的距離;ω為頻率;A為給定激光偏振情況下激光對液體的透射率;i為入射角;φ為折射角;n1為空氣中的折射率;n2為液體中的折射率;I0為光強;a為光束半徑;I(ω)為激光頻譜;E(ω)為單位時間吸收并轉(zhuǎn)化為熱量的電磁波能量密度頻譜,可表示為:

(6)

式中:k為波數(shù),與激光波長有關(guān);μ為水對激光的吸收系數(shù);θ為觀測角,即r與垂直向下方向的夾角;τ0為垂直特性延遲時間;τ1為水平特性延遲時間。

2 n-MFSK調(diào)制

n-MFSK調(diào)制是在多進制頻移鍵控(n-frequency shift keying,n-FSK)調(diào)制基礎(chǔ)上,進一步提高通信速率并改善誤比特率。本文中針對激光致聲空-水跨介質(zhì)長脈沖法和重復(fù)頻率法,提供n-MFSK調(diào)制方案,整體激光致聲通信原理圖如圖1所示。

圖1 激光致聲通信調(diào)制解調(diào)原理圖Fig.1 Modulation and demodulation principles diagram for laser induced acoustic communication

2.1 長脈沖法n-MFSK調(diào)制

長脈沖法是利用某種調(diào)制方法對長脈沖激光進行調(diào)制的方法,本文中根據(jù)參考文獻[16]來選取熱膨脹效應(yīng)下的光聲轉(zhuǎn)換仿真參數(shù),具體參數(shù)如表1所示。使用Nd∶glass激光器,長脈沖激光波形近似為方形脈沖,在邊緣逐漸下降,該波形可在極端聚焦情況下減少非線性效應(yīng)波形突變的影響。

表1 長脈沖法仿真參數(shù)[16]Table 1 Simulation parameters of long-pulse-duration laser method[16]

由于水聲信道存在多徑效應(yīng),需選取合適的符號間隔以減輕符號間干擾,故本文中選取n-FSK調(diào)制,該調(diào)制系統(tǒng)由多臺激光器構(gòu)成激光陣列,通過控制激光器發(fā)射頻率減少由多徑引起的符號間干擾。以2-FSK為例,將其中一個調(diào)制頻率記為1,另一個調(diào)制頻率記為0,兩個調(diào)制頻率中心頻率都為最佳中心頻率。每臺激光器,符號持續(xù)時間為2 ms,符號間隔為98 ms,通過機械斬波器控制激光頻率變化,例如,傳輸速率為100 bits/s時,需由10臺激光器組成陣列,以激光器發(fā)射頻率并行發(fā)射。該方案可通過增加調(diào)制頻率數(shù)量,以減少激光器數(shù)量并提升傳輸速率,但所需帶寬隨傳輸速率呈非線性增加。理論上可采用增加符號持續(xù)時間,減少所占帶寬,為此需減少激光器發(fā)射功率,以增加脈沖持續(xù)時間,這會導(dǎo)致聲波在水中的聲壓級(sound pressure level,SPL)和傳輸距離減少。此外,隨著調(diào)制頻率數(shù)量增加,會降低對多徑效應(yīng)的抑制,可通過信道編碼減少多徑效應(yīng),但會降低傳輸速率。

為此,本文作者提出了另一種調(diào)制方式,即n-MFSK調(diào)制,在一個符號間隔期間發(fā)射多個頻率,以增加調(diào)制頻率數(shù)量。以2-MFSK為例,將兩臺相同激光器基頻設(shè)為f1和f2,其中,當兩臺激光器均不發(fā)射激光時,記為00;當激光器1發(fā)射頻率為f1、激光器2不發(fā)射時,記為01;當激光器1不發(fā)射、激光器2發(fā)射頻率為f2時,記為10;當激光器1發(fā)射頻率為f1、激光器2發(fā)射頻率為f2時,記為11。當傳輸速率達到100 bits/s時,需要5臺激光器組成陣列。2-FSK和2-MFSK調(diào)制的時域和頻域結(jié)果分別如圖2a和圖2b所示。

圖2 長脈沖法時域和頻域仿真結(jié)果Fig.2 Long-pulse-duration laser time domain and frequency domain simulation results

該方案的優(yōu)點為通過激光器頻率組合,增加每個符號傳輸比特數(shù)來提高傳輸速率,并且在使用開關(guān)鍵控時,可提高帶寬利用率。

根據(jù)表1中的參數(shù),計算不同調(diào)制方法下長脈沖法水下通信距離為[17]:

(7)

式中:Xs為聲源級;Xt為聲傳播損失;Xn為背景噪聲;RSNR為接收器檢測信噪比;α為海水吸收系數(shù)。

表2為長脈沖法n-FSK和n-MFSK調(diào)制方法通信速率、帶寬、聲壓級和水中傳輸距離的比較。其中背景噪聲取3級海況,RSNR=10 dB,α=0.02 dB/km,聲壓級作用范圍以137 dB為基準(依據(jù)表1計算)。隨著調(diào)制頻率數(shù)量增加,n-MFSK調(diào)制可提高通信速率,并改善頻帶利用率,但聲壓級和水中傳輸距離會隨符號持續(xù)時間內(nèi)調(diào)制頻率數(shù)量增加而減小。比較4-FSK和2-MFSK調(diào)制,兩者可達到相同的通信速率,且2-MFSK只需要1 kHz的帶寬,然而在最壞情況下,當一個符號持續(xù)時間傳輸兩個頻率時,激光能量必須除以2,因此聲壓級減小到131 dB。這導(dǎo)致水中傳輸距離從4-FSK的708 m減小到355 m,隨著一個符號時間內(nèi)傳輸頻率數(shù)量增加,水中傳輸距離將大幅度減少。因此,當采用n-MFSK調(diào)制時,需權(quán)衡通信速率與水中傳輸距離的關(guān)系,根據(jù)實際需求選擇合適的調(diào)制方式。

表2 長脈沖法n-FSK和n-MFSK調(diào)制的比較Table 2 Comparison of n-FSK and n-MFSK modulation by long-pulse-duration laser method

2.2 采用重復(fù)頻率法的n-MFSK調(diào)制

重復(fù)頻率法是基于激光致聲的熱膨脹效應(yīng)通過控制激光器的重復(fù)頻率產(chǎn)生調(diào)制聲信號的方法。由于單激光脈沖產(chǎn)生的聲信號頻譜所占帶寬大,頻帶利用率低,其信號能量會在水下發(fā)射嚴重衰減,易產(chǎn)生畸變,不利于傳輸編碼。因此,需采用產(chǎn)生窄帶且中心頻率可控的聲信號,即高重復(fù)率的激光脈沖信號。

當激光垂直射向水面時,激光脈沖重復(fù)率為fr=1/T,T為激光脈沖周期時,激光脈沖的頻譜可表示為[20]:

(8)

式中:I0(ω)為單脈沖的頻譜;G(ω)頻譜由于諧波關(guān)系呈現(xiàn)梳狀結(jié)構(gòu);N為階數(shù);fr為相鄰譜峰的間距,且其周期可控。將式(8)代入式(4)可得[21]:

(9)

式中:E(ω)與μ、a、θ等參數(shù)有關(guān),其中激光器參數(shù)μ和a是固定的,故只需控制fr使s(r,ω)窄帶特性最好即可。選擇激光器發(fā)射重復(fù)率fr時,要使聲信號高次諧波被抑制且中心頻率處的能量強,因此,激光器的重復(fù)頻率需滿足[20]:

(10)

式中:s(fr)為中心頻率的幅度;s(fmax)為頻譜幅度的最大值;max|s(nfr)|為高次諧波頻譜幅度的最大值,由此,可確定激光器發(fā)射重復(fù)率的范圍。根據(jù)表3中的參數(shù),采用調(diào)QNd∶YAG高重頻激光器,并采用靈敏度為-170 dB的光纖水聽器對聲波進行接收,聯(lián)合式(8)～式(10)計算fr的范圍可得:

表3 重復(fù)頻率法仿真參數(shù)[20]Table 3 Simulation parameters of high repetitive rate method[20]

12 kHz≤fr≤31 kHz

(11)

目前現(xiàn)有的激光器難以從實驗的角度進行高速通信研究,由于相鄰激光脈沖產(chǎn)生的聲信號具有很好的重復(fù)性,故可采用單個或多個實驗信號進行重頻組合,獲得所需的高重復(fù)率激光。圖3a中將單脈沖激光致聲信號以fr=21 kHz重頻組合,圖3b是依據(jù)式(9)計算后得到結(jié)果和單脈沖重頻組合信號頻譜進行對比。兩者頻譜特征基本吻合,驗證了該方案的可行性。

圖3 a—單脈沖激光信號時域波形 b—重頻組合信號時域波形 c—I0(ω)G(ω)與單脈沖重頻組合頻域波形Fig.3 a—single pulse laser signal time domain waveform b—re-frequency combination signal time domain waveform c—spectrogram with single laser signal re-frequency combination and I0(ω)G(ω)

通過控制激光器的重復(fù)頻率使其有規(guī)律地變化,從而完成n-FSK和n-MFSK調(diào)制以搭載信息,n-FSK調(diào)制中一個碼元信號以重復(fù)頻率fr的N個激光脈沖激發(fā)聲信號表示,在符號間隔(0,T0)中,碼元信號可表示為:

si(t)=Re[IFFT(s(fr,i))],(i=1,2,…,n)

(12)

故n-FSK調(diào)制信號可表示為:

(i=1,2,…,n)

(13)

式中:IFFT(s(fr,i))表示對s(fr,i)進行的傅里葉逆變換(inverse fast Fourier transform,IFFT);gτ(t-mT0)是脈寬為τ、周期為T0的門函數(shù)。

n-MFSK調(diào)制中每個碼元含有零個、一個或多個調(diào)制頻率,單個碼元信號可表示為:

(i=1,2,…,n)

(14)

式中:ni取值為0或1,取0代表該碼元含有調(diào)制頻率fr,i,取1代表不含調(diào)制頻率fr,i。

n-MFSK調(diào)制信號可表示為:

(i=1,2,…,n)

(15)

通過FSK和MFSK調(diào)制信號搭載信息,實現(xiàn)激光信號的編碼傳輸。圖4a和圖4b分別對應(yīng)4-FSK和2-MFSK調(diào)制時域和頻域波形,符號間隔為50 ms。

圖4 重復(fù)頻率法時域波形Fig.4 High repetitive rate method time domain waveforms

4-FSK調(diào)制頻率fr,i根據(jù)式(12)重復(fù)頻率范圍等間隔取值,分別取為13 kHz、17 kHz、21 kHz和25 kHz,代表4個碼元,每個碼元用二進制依次表示為{00,01,10,11},碼元帶寬約為500 Hz。2-MFSK調(diào)制頻率fr,i分別取為17 kHz和25 kHz,第1個碼元不包含頻率信息為fr,0,第2個碼元含有1個17 kHz頻率信息為fr,1,第3個碼元含有1個25 kHz頻率信息為fr,2,第4個碼元含有兩個17 kHz和25 kHz頻率信息為fr,1,2,每個碼元用二進制依次表示為{00,01,10,11}。如圖3所示,2-MFSK調(diào)制可減少基頻數(shù)量,同時提高通信速率并提高帶寬利用率。

表4 重復(fù)頻率法n-FSK和n-MFSK調(diào)制的比較Table 4 Comparison of high repetitive rate method n-FSK and n-MFSK modulation

表4為重復(fù)頻率法n-FSK和n-MFSK調(diào)制的比較。參數(shù)計算方式與長脈沖法相同,聲壓級作用范圍以139 dB為基準(依據(jù)表2計算)。對比n-FSK和n-MFSK調(diào)制,在相同調(diào)制頻率數(shù)量時,n-MFSK調(diào)制可提高通信速率,并提高頻帶利用率,與此同時代價是水中通信距離的損失。

3 結(jié) 論

對提升激光致聲空-水跨介質(zhì)通信速率的方法開展了仿真研究工作,得到以下研究結(jié)果。

(a)n-FSK調(diào)制可有效減少符號間的多徑干擾,長脈沖法采用激光器陣列的方式實現(xiàn),而重復(fù)頻率法采用激光器頻率切換的方式實現(xiàn),兩種方法都可通過增加調(diào)制頻率數(shù)量提高通信速率,并保持聲壓級和水下通信范圍不變,但所需帶寬會隨傳輸速率呈非線性增加。

(b)n-MFSK調(diào)制可在n-FSK調(diào)制的基礎(chǔ)上進一步提高通信速率,長脈沖法主要通過激光陣列的方式,并采用頻率疊加的方式實現(xiàn),重復(fù)頻率法通過控制激光器頻率變化時域上分配時間段的方式實現(xiàn)。隨著調(diào)制頻率數(shù)量增加,n-MFSK調(diào)制可提高通信速率,并改善頻帶利用率,但聲壓級和水中傳輸距離會隨符號持續(xù)時間內(nèi)調(diào)制頻率數(shù)量的增加而減小。

(c)n-FSK調(diào)制可保證通信速率提升,不會對水中傳輸距離造成損失。相比n-FSK,n-MFSK調(diào)制可提高通信速率并改善頻帶利用率,但水中傳輸距離會受調(diào)制頻率數(shù)量影響。實際使用時,需在通信速率與水中傳輸距離之間進行權(quán)衡,以實際需求作為調(diào)制方式選擇基準。

本文作者提出了一種n-MFSK調(diào)制方式,進一步提升激光致聲空-水跨介質(zhì)通信速率,分別采用長脈沖法和重復(fù)頻率法的調(diào)制方式,結(jié)合現(xiàn)有激光器的技術(shù)指標分析了通信性能,計算了通信速率的提升效果,研究結(jié)果對激光致聲在空-水跨介質(zhì)通信領(lǐng)域中的實際應(yīng)用有一定價值。提供的可選擇方案有助于未來實現(xiàn)空中平臺與水下目標之間的通信,長脈沖激光陣列的方式需搭載于機載平臺,且占據(jù)較大空間,而重復(fù)頻率激光頻率切換的方式,目前受激光器的重復(fù)頻率等性能指標的限制,隨著激光器工藝的提升可作為未來空-水跨介質(zhì)通信的方案之一。

下一步研究將針對實際海試環(huán)境因素影響,特別是環(huán)境噪聲對于通信誤比特率的影響,開展環(huán)境噪聲對于通信質(zhì)量影響的評價工作,包括海洋背景噪聲、潛行器本身的螺旋槳運行的機械噪聲等。