肖金光，周新力，張 燁

(海軍航空工程學(xué)院電子信息工程系，山東煙臺(tái)264001)

1 引言

海面上蒸發(fā)波導(dǎo)(Evaporation Duct，ED)發(fā)生概率較高，能夠陷獲頻率和仰角滿足一定條件的電磁波，形成超視距傳播，對(duì)于雷達(dá)探測(cè)、通信、電子對(duì)抗以及低高度武器或平臺(tái)的攻防都具有重要意義。其中利用波導(dǎo)進(jìn)行通信可行性已被美、澳大利亞等開(kāi)展的試驗(yàn)所證實(shí)，對(duì)其性能的研究和系統(tǒng)設(shè)計(jì)也在持續(xù)開(kāi)展[1－2]。射線追蹤法(Ray Tracing，RT)作為一種基于幾何光學(xué)理論的高頻近似方法，雖然對(duì)衍射和散射的忽略造成了場(chǎng)強(qiáng)計(jì)算誤差，但估計(jì)多徑時(shí)間上具有明顯優(yōu)勢(shì)，因而成為一種非常流行的預(yù)測(cè)無(wú)線信道特性的方法[3]。文獻(xiàn)[4]給出了基于幅度和多徑延遲時(shí)間的信道傳輸函數(shù)，但未考慮海面的反射作用。文獻(xiàn)[5]基于RT，計(jì)算了多徑延遲和到達(dá)信號(hào)幅度，幅度計(jì)算中用粗糙因子修正表面反射造成的場(chǎng)強(qiáng)衰減，但未考慮電波的傳播擴(kuò)散，且給出了波導(dǎo)強(qiáng)度而不是大氣折射率剖面，亦即求解多徑時(shí)間時(shí)未考慮射線在非均勻大氣中傳播速度的不均勻。文獻(xiàn)[6]基于容量有限的Markov生滅過(guò)程研究大氣波導(dǎo)中到達(dá)接收站的射線徑數(shù)，利用MUSIC算法對(duì)多徑到達(dá)角進(jìn)行最佳估計(jì)。

綜上所述，基于RT的蒸發(fā)波導(dǎo)通信研究中，場(chǎng)強(qiáng)計(jì)算和多徑特性估計(jì)存在諸多問(wèn)題有待解決。本文對(duì)蒸發(fā)波導(dǎo)通信帶限信道建模信號(hào)的場(chǎng)強(qiáng)計(jì)算、多徑特性、蒸發(fā)波導(dǎo)超視距通信系統(tǒng)設(shè)計(jì)使用策略、帶限沖激響應(yīng)方法進(jìn)行了研究，并進(jìn)行了數(shù)值試驗(yàn)。

2 基于SRT的蒸發(fā)波導(dǎo)信道建模方法

2.1 分段射線追蹤法

RT認(rèn)為電磁波在非常小的射線管內(nèi)傳播，跟蹤每一條射線傳播軌跡，有傳統(tǒng)積分和泰勒級(jí)數(shù)近似兩種模式。鑒于海洋蒸發(fā)波導(dǎo)高度一般不高于40 m，且本文的研究重點(diǎn)是電波的超視距傳播問(wèn)題，因此選用適于低仰角電波傳播計(jì)算的泰勒級(jí)數(shù)近似法。

對(duì)大氣折指數(shù)n引入平坦地球近似，將電波沿地表的傳播等效為水平傳播，則修正大氣折射指數(shù)為

式中，re為地球半徑。修正大氣折射率M為

設(shè) r1、h1、a1和 r2、h2、a2分別為射線初始位置和到達(dá)位置的距離、高度及仰角，ΔM為折射率變化量，g為修正折射率M在分層內(nèi)的線性變化率，RT的求解公式組為

在某一折射率分層內(nèi)部，上述射線的軌跡就會(huì)有三種情況:單調(diào)向上、單調(diào)向下或者在層內(nèi)拐彎(達(dá)到最高或最低點(diǎn))。如果進(jìn)入別的分層則以進(jìn)入點(diǎn)作為終點(diǎn)，而對(duì)于層內(nèi)拐彎射線(根據(jù)初始點(diǎn)和終點(diǎn)的角度符號(hào)相反可以判斷)，從最高點(diǎn)或最低點(diǎn)處拆成兩段分別計(jì)算。這樣就統(tǒng)一成單調(diào)向上和向下兩種情況，這種大氣分層條件下分段射線追蹤(Subsection Ray Tracing，SRT)計(jì)算方法，既簡(jiǎn)化了運(yùn)算，又便于進(jìn)行傳播路徑長(zhǎng)度和時(shí)間的計(jì)算。

單調(diào)向上和層內(nèi)達(dá)到最高點(diǎn)兩種情況如圖1所示，而單調(diào)向下和層內(nèi)達(dá)到最低點(diǎn)的情況具有與圖1對(duì)稱的形式。

圖1 分段射線追蹤Fig.1 Subsection ray tracing

從而，根據(jù)式(3)～(5)可以由射線初始位置、高度、仰角以及折射率分層情況，求得射線通過(guò)的路徑。

2.2 蒸發(fā)波導(dǎo)模型和電波陷獲傳播條件

中性大氣層結(jié)中蒸發(fā)波導(dǎo)的大氣修正折射率剖面為[7]

式中，表面粗糙參數(shù)z0為1.5×10－4m，d為蒸發(fā)波導(dǎo)高度(單位為m)，M(z0)取為330。若取蒸發(fā)波導(dǎo)高度30 m，天線高度25 m，基于SRT的射線軌跡如圖2所示(只繪出了超視距傳播的射線)。

圖2 射線軌跡Fig.2 Tracks of rays

對(duì)于一定強(qiáng)度/高度的波導(dǎo)，頻率大于某一截止頻率fmin(即波長(zhǎng)小于某一截止波長(zhǎng)，單位Hz)的電磁波才能形成陷獲傳播，如圖3所示。截止頻率與波導(dǎo)高度的近似關(guān)系可以表示為[8]

圖3 最低陷獲頻率Fig.3 Minimum trapped frequency

2.3 海洋蒸發(fā)波導(dǎo)中電波傳播計(jì)算

基于前面給出的大氣分層條件下SRT計(jì)算方法，可以方便地計(jì)算射線的路徑長(zhǎng)度。射線弧與水平位移長(zhǎng)度之差ΔS和射線長(zhǎng)度S分別為

電波在某段路徑上的傳播時(shí)間Δt由射線長(zhǎng)度s及傳播速度v決定:

式中，c0為真空波速。電波強(qiáng)度由初始場(chǎng)、粗糙面反射系數(shù)、空間擴(kuò)散決定，到達(dá)強(qiáng)度由接收球法確定。初始場(chǎng)為E(φ)=f(φ)E0，其中f(φ)為角度 φ 射線對(duì)應(yīng)的天線方向圖增益系數(shù)，全向天線時(shí)為1?？罩袀鞑ド渚€的場(chǎng)強(qiáng)為E2=E1eiks/s，達(dá)到地表射線反射前后的場(chǎng)強(qiáng)關(guān)系為E2=ΓE1，其中Γ為菲涅耳反射系數(shù)，是入射余角的函數(shù)(由射線在表面上的入射余角決定)[9]。因此，假設(shè)射線發(fā)生L次反射，則射線經(jīng)過(guò)路徑S后的場(chǎng)強(qiáng)為

二維空間上某點(diǎn)的場(chǎng)強(qiáng)為所有到達(dá)接收球射線的強(qiáng)度和，接收球半徑R取為[3]

式中， 為輻射源出發(fā)相鄰射線的張角，單位為rad。多徑信號(hào)的接收球法，需要清晰地區(qū)分各路徑，然后針對(duì)每一路徑族，使用接收球法計(jì)算該路徑的信號(hào)強(qiáng)度，不妨稱之為分離多徑信號(hào)族接收球法。

2.4 蒸發(fā)波導(dǎo)帶限信道建模方法

由圖1可知，蒸發(fā)波導(dǎo)中，射線經(jīng)不同路徑到達(dá)接收位置，是典型的多徑信道，其復(fù)低通等效沖激響應(yīng)可以描述為[4]

式中，Ei和τi分別為第i條路徑到達(dá)信號(hào)的幅度和到達(dá)時(shí)間，F(xiàn)t、Fr為發(fā)收天線方向圖量化因子，θt為射線初始角，θr為到達(dá)角。若收發(fā)天線均采用通信常用的全向天線時(shí)，F(xiàn)t和Fr均為1，則式(13)可簡(jiǎn)化為

蒸發(fā)波導(dǎo)中的超視距通信不存在強(qiáng)度占絕對(duì)優(yōu)勢(shì)的視距路徑。為了防止多徑時(shí)延對(duì)接收造成不可逆的影響，信號(hào)符號(hào)周期Ts與最大多徑時(shí)延τmax需要滿足

為避免多徑引起的頻率選擇性衰落，系統(tǒng)相干帶寬Δf可取為

τmax決定了可通信道的最大帶寬，相當(dāng)于頻域加窗，等效于時(shí)域乘以sinc函數(shù)，因此復(fù)低通表示的信道帶限沖激響應(yīng)為

3 信道特性仿真分析與建模

從源點(diǎn)±1°夾角內(nèi)等間隔發(fā)射3000條射線，蒸發(fā)波導(dǎo)高度30 m，海面風(fēng)速5 m/s，發(fā)射天線高度10 m、20 m和25 m時(shí)，電波到達(dá)100 km處的高度與時(shí)延分布關(guān)系如圖4～6所示。

圖4 10 m發(fā)射天線的多徑信號(hào)時(shí)延分布Fig.4 Distribution of multipath signal delay with radiating antenna at 10 m

圖5 20 m發(fā)射天線的多徑信號(hào)時(shí)延分布Fig.5 Distribution of multipath signal delay with radiating antenna at 20 m

圖6 25 m發(fā)射天線的多徑信號(hào)時(shí)延分布Fig.6 Distribution of multipath signal delay with radiating antenna at 25 m

由圖7可知，增加射線數(shù)只是增加了到達(dá)接收位置射線的密度。因?yàn)橛稍c(diǎn)出發(fā)射線的初始角是連續(xù)小角度變化的，相鄰射線經(jīng)歷的傳播路徑差異不大，導(dǎo)致到達(dá)信號(hào)的強(qiáng)度、到達(dá)時(shí)間、到達(dá)角、到達(dá)高度連續(xù)變化。增大射線條數(shù)即減小初始角變化量，表現(xiàn)為到達(dá)信號(hào)參量的插值，這正說(shuō)明了式(12)接收球法進(jìn)行場(chǎng)強(qiáng)計(jì)算的原理。雖然不用置疑的是RT忽略的衍射造成了誤差，但對(duì)于蒸發(fā)波導(dǎo)超視距多徑信號(hào)，重要的是各徑信號(hào)的相對(duì)強(qiáng)度。為了適應(yīng)蒸發(fā)波導(dǎo)高度隨氣象條件的變化，保證較高的可通率，必須使得系統(tǒng)能夠滿足最惡劣條件下的通信。

發(fā)射天線高度仍為20 m，在加倍射線數(shù)基礎(chǔ)上，再將接收距離加倍為200 km，電波到達(dá)高度與時(shí)延分布關(guān)系如圖8所示。

由圖4～6可知，波導(dǎo)高度相同，發(fā)射天線高度越接近波導(dǎo)高度，多徑信號(hào)到達(dá)接收距離時(shí)同一接收高度上的徑數(shù)越少，信號(hào)的最大多徑時(shí)延(即圖中接收高度上所有多徑信號(hào)的時(shí)間跨度)越小，反之亦反。天線越低，射線跳數(shù)增加，多徑情況越復(fù)雜，射線族越難分離。因此，分離多徑信號(hào)族接收球法應(yīng)用于天線高度極低情況時(shí)難度較大，此時(shí)可通過(guò)τmax粗略估計(jì)系統(tǒng)性能，好在一般艦船用通信天線高度不會(huì)低至貼近海面。

發(fā)射天線高度20 m，加倍射線數(shù)目，電波到達(dá)100 km處的高度與時(shí)延分布關(guān)系如圖7所示。

圖8 接收距離加倍時(shí)多徑信號(hào)時(shí)延分布Fig.8 Distribution of multipath signal delay with receiving distance doubled

由圖5和圖8可知，徑數(shù)和最大多徑時(shí)延隨距離增加而增大，但距離和最大多徑時(shí)延增量不是直接的倍數(shù)關(guān)系。對(duì)于圖8，根據(jù)式(12)可知接收球半徑約為0.29 m。則16 m和20 m接收高度上，多徑信號(hào)的相對(duì)功率G延遲剖面如圖9和圖10所示。

圖7 射線數(shù)加倍時(shí)多徑信號(hào)時(shí)延分布Fig.7 Distribution of multipath signal delay with ray number doubled

圖9 16 m高度上相對(duì)功率延遲剖面Fig.9 Relative power delay profile at 16 m

圖10 20 m高度上相對(duì)功率延遲剖面Fig.10 Relative power delay profile at 20 m

由圖9和圖10可知，在16 m和20 m高度上，最大多徑時(shí)延為0.885 ns，兩種情況下主次徑相對(duì)強(qiáng)度差可達(dá)26.8 dB，因此適當(dāng)選擇接收高度，可以使接收多徑信號(hào)的相對(duì)強(qiáng)度達(dá)到最大，更加容易克服多徑延遲對(duì)信號(hào)接收造成的影響，也可以用于指導(dǎo)自適應(yīng)設(shè)計(jì)均衡器，實(shí)現(xiàn)更好的接收性能。

對(duì)于高度20 m上進(jìn)行信號(hào)接收，取信號(hào)的符號(hào)周期 Ts為臨界值，即0.885 ns，由式(16)可知系統(tǒng)相干帶寬Δf為1.1 GHz。工程中通常采取冗余設(shè)計(jì)的辦法降低信號(hào)處理的復(fù)雜度，不妨取符號(hào)周期為4、10倍的τmax以作對(duì)比，相應(yīng)的通信帶寬分別為280 MHz和110 MHz，則帶限沖激響應(yīng)如圖11所示。

圖11 帶限沖激響應(yīng)Fig.11 Bandwidth impulse response

因此，根據(jù)蒸發(fā)波導(dǎo)高度和預(yù)期的通信距離，基于SRT和海洋多徑信號(hào)場(chǎng)強(qiáng)計(jì)算方法，仿真可得最大時(shí)延和功率延遲剖面，從而可得系統(tǒng)可用帶寬和相應(yīng)的帶限沖激響應(yīng)，根據(jù)輻射信號(hào)和帶限沖激響應(yīng)，評(píng)估系統(tǒng)的誤碼性能。當(dāng)然，也可以據(jù)此自適應(yīng)設(shè)計(jì)均衡器，進(jìn)一步提高通信系統(tǒng)的誤碼性能。這對(duì)于基于軟件無(wú)線電技術(shù)建立速率自適應(yīng)數(shù)據(jù)通信，實(shí)現(xiàn)蒸發(fā)波導(dǎo)傳輸信道特性的高效利用具有重要指導(dǎo)意義。

4 結(jié)論

本文提出了一種基于RT進(jìn)行蒸發(fā)波導(dǎo)通信中多徑信號(hào)場(chǎng)強(qiáng)計(jì)算方法，給出了大氣分層條件下SRT計(jì)算方法，綜合考慮了電波空間擴(kuò)散、粗糙海面反射和接收球法，計(jì)算電波到達(dá)強(qiáng)度，解決了大氣波導(dǎo)內(nèi)部多徑信號(hào)強(qiáng)度難以計(jì)算的問(wèn)題;研究了蒸發(fā)波導(dǎo)內(nèi)部多徑信號(hào)的特性;根據(jù)電波形成超視距傳播條件和多徑特性，研究了蒸發(fā)波導(dǎo)超視距通信系統(tǒng)設(shè)計(jì)使用策略;給出了基于多徑特性建立帶限沖激響應(yīng)的方法，為蒸發(fā)波導(dǎo)通信的鏈路預(yù)算提供了有效的途徑，通過(guò)控制路徑損耗，維持正余量，達(dá)到期望的系統(tǒng)差錯(cuò)性能，也可以據(jù)此自適應(yīng)設(shè)計(jì)均衡器，進(jìn)一步提高通信系統(tǒng)的誤碼性能。仿真表明，一整套蒸發(fā)波導(dǎo)通信帶限信道建模方法切實(shí)可行。蒸發(fā)波導(dǎo)通信前景誘人，適于波導(dǎo)高發(fā)區(qū)的島嶼、艦、岸之間移動(dòng)通信節(jié)點(diǎn)的通信，數(shù)據(jù)速率和可通率較高，性價(jià)比優(yōu)勢(shì)明顯。下一步將進(jìn)行系統(tǒng)設(shè)計(jì)并開(kāi)展試驗(yàn)驗(yàn)證。

[1]任香凝，李文計(jì).海面蒸發(fā)波導(dǎo)微波超視距通信可行性分析[J].信息傳輸與接入技術(shù)，2008，34(2):22 －24.REN Xiangning，LI Wenji.Feasibility analysis on beyond－h(huán)orizon microwave communication with evaporation duct over the sea[J].Radio Communications Technology，2008，34(2):22 －24.(in Chinese)

[2]Woods G S.High －capacity，long－ range，over ocean microwave link[J].IEEE Journal of Oceanic Engineering，2009，34(3):323 －330.

[3]呂英華.計(jì)算電磁學(xué)的數(shù)值方法[M].北京:清華大學(xué)出版社，2006:310 －315.LYU Yinghua.Numerical methods in computational electromagnetics[M].Beijing:Tsinghua University Press，2006:310 －315.(in Chinese)

[4]陳紹賀，鄭輝，黃興忠.大氣波導(dǎo)中多徑衰落信道參數(shù)的定量計(jì)算[J].電波科學(xué)學(xué)報(bào)，2010，25(4):638 －645.CHEN Shaohe，ZHENG Hui，HUANG Xingzhong.Numerical computation of parameters of multipath fading channel in the atmospheric duct[J].Chinese Journal of Radio Science，2010，25(4):638 －645.(in Chinese)

[5]趙樓，楊坤德，楊益新.海洋蒸發(fā)波導(dǎo)信道的多徑時(shí)延[J].探測(cè)與控制學(xué)報(bào)2010，32(1):39 －44.ZHAO Lou，YANG Kunde，YANG Yixin.Multi－ path effect of Evaporation Duct[J].Journal of Detection ＆Control，2010，32(1):39 －44.(in Chinese)

[6]趙小龍，黃際英.大氣波導(dǎo)中多徑信道的參數(shù)研究[J].西安電子科技大學(xué)學(xué)報(bào)(自然科學(xué)版)，2008，35(2):314－318.ZHAO Xiaolong，HUANG Jiying.Research on multi －path channel parameters in the Marine atmospheric duc[J].Journal of Xidian University(Natural Science Edition)，2008，35(2):314 －318.(in Chinese)

[7]Sirkova I.Brief review on PE method application to propagation channel modeling in sea environment[J].Central European Journal of Engineering，2012，2(1):19 －38.

[8]Katherine L.A smart climatology of evaporation duct height and surface radar propagation in the INDIAN ocean[D].Monterey，CA:Naval Postgraduate School，2007.

[9]肖金光，周新力，劉曉娣.彎折譜估計(jì)在PE模型電波傳播中的應(yīng)用研究[J].電波科學(xué)學(xué)報(bào)，2014，29(3):559－566.XIAO Jinguang，ZHOU Xinli，LIU Xiaodi.Application of cured wave spectral estimation in wave propagation with parabolic equation model[J].Chinese Journal of Radio Science，2014，29(3):559 －566.(in Chinese)