賈建偉

(江蘇聯(lián)合職業(yè)技術(shù)學(xué)院揚(yáng)州分院，江蘇 揚(yáng)州225002)

0 引 言

隨著海洋開發(fā)和海洋科技的發(fā)展，水聲信號處理技術(shù)變得非常活躍，這種技術(shù)在海洋環(huán)境監(jiān)測、海洋資源勘探、軍事演習(xí)等方面有舉足輕重的作用。由于海洋環(huán)境千變?nèi)f化，水聲信號處理往往受到風(fēng)浪流的的影響，因此研究水聲新的定位，以及對信號進(jìn)行處理能夠更有效地提取有用信息進(jìn)行數(shù)據(jù)處理分析。

本文通過研究海上水聲信號的產(chǎn)生和海上水聲信號定位系統(tǒng)的構(gòu)成，設(shè)計(jì)水聲信號處理系統(tǒng)流程，指出陸地上的主控中心與水中的基站如何進(jìn)行信號數(shù)據(jù)傳遞，以及環(huán)境噪聲、傳輸損耗、多徑傳輸、多普勒效應(yīng)等物理特性對信號偏移產(chǎn)生何種影響。最后通過實(shí)驗(yàn)進(jìn)行驗(yàn)證，為了解決當(dāng)子載波在頻域內(nèi)不正交時(shí)的水聲信道傳輸問題，在接收端加入濾波器通過SINR 增益值來說明如何消除多普勒偏移的影響。

1 海上水聲信號的產(chǎn)生

本文利用圖形化編程語言LabView 海上水聲信號的仿真。現(xiàn)實(shí)中，海上水聲信號為頻帶寬度較小的信號，表達(dá)式為［1］:

式中:a(t)為載波幅度;φ(t)為載波相位;f0為信號載頻;T 為工作周期。

本文進(jìn)行聲吶信號仿真的方式為采用正旋單頻信號對FM 進(jìn)行調(diào)頻以及對DSB 進(jìn)行單邊幅度調(diào)制。整個(gè)過程如下［2］:

令接收到的調(diào)制信號為m(t)，載波信號為p(t)=Acos(2πfct)，利用調(diào)制信號去控制載波信號的頻率得到的FM 信號為:

于是得:

設(shè)a(t)=m(t)對載波信號進(jìn)行第2 次調(diào)制，然后用DSB 去調(diào)制sFM(t)，因此獲得模擬水聲信號s(t)。

2 水聲信號的定位以及處理

1)海上水聲信號定位系統(tǒng)的構(gòu)成

水聲信號定位的方式為通過無線電進(jìn)行聲定位。有陸面上的控制中心和海中的基站組成。PC 機(jī)、數(shù)據(jù)采集器、無線收發(fā)器組成了陸上的控制中心，傳感器、無線收發(fā)器、水中控制器主城了海中的基站。在進(jìn)行信號定位時(shí)，陸上的控制中心將各個(gè)控制信號通過無線電裝置傳送到水中的基站，由此可以控制基站中裝置的初始值、校驗(yàn)值以及開啟運(yùn)作模式。在基站接收到工作信號后將進(jìn)行數(shù)據(jù)采集，主要是將傳感器采集的信號進(jìn)行隔離、放大、濾波處理，以及將電壓變換到將模數(shù)轉(zhuǎn)換器的范圍內(nèi)，然后發(fā)送給數(shù)字轉(zhuǎn)換器進(jìn)行抽樣、量化、編碼，然后在接收到主控中心的信號后交給陸上處理，然后在PC機(jī)進(jìn)行分析處理，然后利用達(dá)到時(shí)間差定位方式進(jìn)行定位，實(shí)現(xiàn)方式如圖1 所示。

圖1 水聲信號定位實(shí)現(xiàn)方式構(gòu)成圖Fig.1 Structure diagram of acoustic signal orientation implementation

2)水聲信號處理系統(tǒng)流程

主控中心開始運(yùn)行時(shí)給海中的基站發(fā)送指令，讓其初始化，然后將返回的信息顯示在主控的PC機(jī)上;系統(tǒng)各個(gè)模塊處于待接收數(shù)據(jù)的狀態(tài)，并且對傳感器、數(shù)據(jù)采集卡等裝置進(jìn)行參數(shù)配置，若有觸發(fā)信號則硬件開始工作，在接收到指令后軟件開始工作，進(jìn)行數(shù)據(jù)的采集、存儲、分析;通過無線電將基站的信息反饋到主控中心，并且存儲在數(shù)據(jù)庫中以備對數(shù)據(jù)進(jìn)行分析整理，結(jié)果在PC 機(jī)上顯示出來。在整個(gè)系統(tǒng)中設(shè)置了安全監(jiān)測實(shí)時(shí)診斷系統(tǒng)是否處于正常運(yùn)轉(zhuǎn)狀態(tài)，當(dāng)發(fā)現(xiàn)出現(xiàn)故障，主控中心第一時(shí)間會發(fā)送指令停止運(yùn)行。

圖2 水聲信號處理流程圖Fig.2 Flow chart of the underwater acoustic signal processing

海域的溫度、鹽分、潮流都會影響傳感器對信號數(shù)據(jù)的采集，因此水聲信號處理對于物理特性的仿真非常重要。

3 水聲信道的物理特性

海上水聲信道是一個(gè)時(shí)頻變化，外界環(huán)境影響大的多重干擾信道，其物理特性主要包括環(huán)境噪聲、傳輸損耗、多徑傳輸、多普勒效應(yīng)引發(fā)的信號偏移。

1)環(huán)境噪聲

海上水聲信號往往受到洋流潮汐船舶活動等噪聲影響。通過頻率與分貝之間的轉(zhuǎn)換可知:海上噪聲為1 ～20 kHz 對應(yīng)的分貝為40 ～70 dB，頻率與噪聲成反比，頻率越小，噪聲越大，若噪聲的頻率在1 kHz 以下，噪聲譜在70 dB 以上，由此可得，傳輸信號時(shí)載波的頻率與海上噪聲有關(guān)，一般將基站的載波頻率設(shè)置為1 kHz 以上。

盡管海上噪聲的大小決定了信噪比和發(fā)射功率，但在做物理特性分析時(shí)，一般將1 ～20 kHz 頻率內(nèi)的噪聲近似為高斯噪聲，如圖3 所示。

圖3 高斯噪聲示意圖Fig.3 Schematic diagram of the gauss noise

2)傳輸損耗

聲波在水中傳播時(shí)由于吸收、散射等問題發(fā)生衰落，當(dāng)聲波頻率大于1 kHz 時(shí)，聲波衰落的原因是海水對聲波的吸收，衰減大小與聲波頻率的平方有關(guān)，海水對聲波吸收系數(shù)經(jīng)驗(yàn)值為［3］:

式中:f 為聲波頻率;α 為每公里的分貝數(shù)。

一般情況下，海上通信設(shè)置的載波頻率小于50 kHz，根據(jù)文獻(xiàn)［4］可知水中通信距離和傳輸速率之間的關(guān)系為:

式中:第1 個(gè)R 為通信距離;第2 個(gè)R 為傳輸速率。

聲信號在擴(kuò)散的過程中會按照一定的規(guī)則衰減，當(dāng)傳播距離較遠(yuǎn)時(shí)，距離越大，聲波會按照R－3/2的規(guī)律衰減［5］。由此可得聲信號的損失為:

式中:R 為傳播距離;α 為聲波的損失系數(shù)(dB/km)。

3)多路徑傳輸

由于海底和海面之間的N 次反射或者折射使得聲源信號在抵達(dá)接收器的時(shí)候會有不同的路徑。同時(shí)會影響信號的幅度變小、頻率變低和解碼之間的干擾增強(qiáng)，從而嚴(yán)重影響了主控中心接收到的信號質(zhì)量。

4)多普勒效應(yīng)產(chǎn)生頻率的偏移

信號發(fā)送端和接收端之間產(chǎn)生的相對位移，引起當(dāng)二者靠近時(shí)頻率增大，二者遠(yuǎn)離時(shí)頻率減小。若發(fā)射器端的運(yùn)動速度為vs，接收端的運(yùn)動速度為vr，可知發(fā)射頻率為fs的信號到接收端后頻率變?yōu)椋?］:

式中:vr和vs為從發(fā)射端指向接收端為正方向的向量;fs為發(fā)射頻率;fr為接受頻率;c 為水中的聲速，一般情況下c = 1 500 m/s。

當(dāng)fs= 15 000 Hz，vs= 3 m/s，vr= －3 m/s，接收頻率fr= 15 060 Hz 時(shí)，多普勒頻率偏移為:

信號發(fā)送端和接收端之間沒有產(chǎn)生相對位移時(shí)，由于風(fēng)浪流等因素的影響，信號會因?yàn)檫@些因素而被調(diào)制，同樣會使接收端的頻率發(fā)生偏移，根據(jù)參考文獻(xiàn)［7］可知頻移為:

式中:fw為波浪的運(yùn)動頻率;hw為波浪的高度均方值;w 為風(fēng)速大小;c 為水中的聲速，一般情況下c =1 500 m/s;θ0為聲波到達(dá)接收端時(shí)的入射角。

假設(shè)海上的風(fēng)速v = 20 kn，則信道的多普勒頻率為28.23 ～43.75 Hz。

由上述可知:海上風(fēng)浪流等形成的噪聲對信號的傳播接受造成很大的影響，使得傳輸損耗增大，接收信號的信噪比也減小，低頻率的控制波使得傳輸速率低，并且多普勒效應(yīng)產(chǎn)生的頻率偏移越大，信道的衰減就越快，接收端接收到的信號波形越不能逼真的反應(yīng)發(fā)送端的數(shù)據(jù)。

4 物理特性仿真

在接收端利用加入濾波器來消除多普勒偏移的影響。

載波的頻率f = 8 kHz，風(fēng)速k = 10 kn，聲速c= 1 500 m/s，水的深度50 m，接收端與發(fā)送端的距離200 m，最大延時(shí)τmax= 0.073 3 s，子載波數(shù)N =1 024，子信道帶寬B = 13.28 Hz，時(shí)間間隔Ts=0.148 6 s。從而得到信號發(fā)送端和接收端之間產(chǎn)生的相對位移時(shí)，經(jīng)補(bǔ)償后剩余的多普勒頻率偏移為Fd= 1 Hz 和Fd= 2 Hz，此時(shí)加入消除濾波器后得到的SINR 增益比較如圖4 所示。

圖4 加入消除濾波器后得到的SINR 增益比較圖Fig.4 Comparison of SINR gain

信號發(fā)送端和接收端之間沒有產(chǎn)生的相對位移時(shí)，由于風(fēng)浪流等因素影響SINR 增益。第i 條路徑經(jīng)補(bǔ)償后余下的多普勒頻率偏移fl= Fdcosθl，式中Fd= 1 Hz 和Fd= 2 Hz，各個(gè)聲波到達(dá)接收端時(shí)的入射角θl分別為0.072 3 π，0.548 6π，0.753 1 π，0.570 8 π，0.032 0 π，0.947 8 π。加入消除濾波器后得到的SINR 增益比較如圖5 所示。

圖5 無相對于運(yùn)動時(shí)加入消除濾波器后得到的SINR 增益比較圖Fig.5 Comparison of SINR gain when no relative movement

從圖4 ～圖5 可知，加入消除濾波器后多普勒頻率偏移的SINR 增益增大，即圖中淺色曲線。并且偏移量越大，加入濾波器后的效果越顯著，從而有效地解決了當(dāng)子載波在頻域內(nèi)不正交時(shí)的水聲信道傳輸問題。

5 結(jié) 語

本文通過研究海上水聲信號的產(chǎn)生和海上水聲信號定位系統(tǒng)的構(gòu)成，設(shè)計(jì)水聲信號處理系統(tǒng)流程，指出陸地上的主控中心與水中的基站如何進(jìn)行信號數(shù)據(jù)傳遞，以及環(huán)境噪聲、傳輸損耗、多徑傳輸、多普勒效應(yīng)等物理特性對信號偏移產(chǎn)生何種影響。最后通過實(shí)驗(yàn)驗(yàn)證加入消除濾波器后可以有效地控制多普勒效應(yīng)產(chǎn)生的頻率偏移。

［1］田坦，劉國枝，孫大軍．聲吶技術(shù)［M］．哈爾濱:哈爾濱工業(yè)大學(xué)出版社，2006．TIAN Tan，LIU Guo-zhi，SUN Da-jun．Sonar technology［M］．Harbin:Harbin Institute of Technology Press，2006．

［2］楊政，楊春武，王明洲，等．基于切片譜的主動水聲信號處理方法［J］．艦船科學(xué)技術(shù)，2012，34(10):77 －80．YANG Zheng，YANG Chun-wu，WANG Ming-zhou，et al．Active underwater acoustical signal processing based on slice spectrum［J］．Ship Science and Technology，2012，34(10):77 －80．

［3］Thorp．Deep ocean sound attenuation in the sub-and lowkilocycle-per-second region ［J］．Acoustic Society of American，1965(38):648．

［4］李啟虎．水聲學(xué)研究進(jìn)展［J］．聲學(xué)學(xué)報(bào)，2001，26(4):295 －301．LI Qi-hu．Advances of research work in underwater acoustics［J］．ACTA Acustica，2001，26(4):295 －301．

［5］列·不列霍夫斯基赫，揚(yáng)．雷桑諾夫．海洋聲學(xué)基礎(chǔ)．朱柏賢，譯．北京:海洋出版社，1985．BREKHOVSKIKH L，LYSANOV Y．Ocean acoustics foundation［M］．ZHU Bo-xian，Translate．Beijing:Ocean Press，1985．

［6］許?。曊Z音通信研究［D］．廈門:廈門大學(xué)，2001．XU Jun．Underwater acoustic voice communication［D］．Xiamen:Amoy University，2001．

［7］STOJANOVIC M．Recent advances in high-speed underwater acoustic communicati-ons［C］//IEEE Jour．of Oceanic Engineering，1996，21(2):125 －136．