存在障礙物時(shí)電波傳播拋物線方程分析及其驗(yàn)證?

魏喬菲 尹成友 范啟蒙

(電子工程學(xué)院,脈沖功率激光國家重點(diǎn)實(shí)驗(yàn)室,合肥 230037)

存在障礙物時(shí)電波傳播拋物線方程分析及其驗(yàn)證?

魏喬菲 尹成友?范啟蒙

(電子工程學(xué)院,脈沖功率激光國家重點(diǎn)實(shí)驗(yàn)室,合肥 230037)

(2017年1月14日收到;2017年3月20日收到修改稿)

雙向拋物線方法主要用于起伏地形下電波傳播問題的計(jì)算,該算法本身無法處理地面存在障礙物,尤其是真實(shí)環(huán)境下障礙物與地面為不同媒質(zhì)的情況.因此本文提出一種用于存在障礙物時(shí)電波傳播計(jì)算的拋物線方程新算法.該方法采用區(qū)域分解,對不同障礙物區(qū)域的場值進(jìn)行分區(qū)計(jì)算,并對計(jì)算結(jié)果進(jìn)行相位修正,從而實(shí)現(xiàn)該情況下空間中場值的計(jì)算.在此基礎(chǔ)上,使用矩量法來精確驗(yàn)證拋物線方法的計(jì)算精度.通過實(shí)例分析,證明了存在障礙物時(shí)新算法的精確性,為之后求解真實(shí)環(huán)境下的電波傳播問題提供了參考.

拋物線方程,電波傳播,區(qū)域分解,矩量法

1 引 言

近年來,雙向拋物線方法[1?3]已經(jīng)逐步成為對流層電波傳播問題研究中最常用的方法.該算法可以計(jì)算不規(guī)則地形特征和不同電磁參數(shù)的地表結(jié)構(gòu)對電波傳播的影響,在不規(guī)則地形的處理方面,拋物線方法[4?7]的計(jì)算精度受限于起伏地形的傾角[8?10],傾角超過一定范圍后,雙向拋物線方法的計(jì)算存在較大誤差.這種下邊界條件的處理方式只是將障礙物當(dāng)起伏地形處理,沒有考慮障礙物內(nèi)存在的場值對空間中場值的影響.因此,當(dāng)?shù)孛娲嬖谡系K物時(shí),拋物線方法計(jì)算存在困難.

在拋物線方法的計(jì)算精度驗(yàn)證方面,通常在不規(guī)則地形下使用AREPS軟件[11]或高頻近似法[12?14]來驗(yàn)證拋物線方法計(jì)算的正確性,高頻近似法主要有菲涅耳積分法及幾何光學(xué)結(jié)合一致性幾何繞射理論法兩種方式.這兩種方法都是基于射線追蹤法來計(jì)算空間中場強(qiáng)的分布,因?yàn)檫@兩種高頻近似法在計(jì)算過程中都忽略被障礙物反射的場值,所以主要用于前向拋物線方法計(jì)算的驗(yàn)證.與此同時(shí),障礙物情況越復(fù)雜,空間中繞射場的場值計(jì)算越困難.除此之外,當(dāng)波源傳播角度較大時(shí),為提高近場區(qū)域的驗(yàn)證精度,Omaki等[15]提出用時(shí)域有限差分方法[16]驗(yàn)證拋物線方法計(jì)算結(jié)果的精確性.

結(jié)合上述情況,本文在原有起伏地形下雙向拋物線方法計(jì)算方式的基礎(chǔ)上,提出了一種新的拋物線方程算法,該算法可以較為精確地計(jì)算地面存在障礙物時(shí)空間中的場值大小.當(dāng)?shù)孛娲嬖谡系K物時(shí),根據(jù)區(qū)域分解原理[17]將障礙物處的場值計(jì)算分為兩部分,并分別增加一段重疊區(qū)域,以滿足障礙物上邊界處場值計(jì)算的連續(xù)性.針對不同區(qū)域的場值計(jì)算,使用不同的拋物線計(jì)算方法.除此之外,在場值的迭代計(jì)算中對總場計(jì)算時(shí)的相位差異和障礙物傳播中造成的相位差異兩方面都進(jìn)行了相位修正.在新算法的驗(yàn)證方面,本文用矩量法[18]來驗(yàn)證近距離下雙向拋物線方法求解電波傳播問題的計(jì)算精度.最后通過一系列實(shí)例分析驗(yàn)證了新算法計(jì)算的精確性.

2 障礙物下拋物線方程的處理方法

設(shè)電磁場的時(shí)諧因子為exp(jωt),本文主要計(jì)算討論的是地面存在障礙物時(shí)雙向拋物線方法的具體求解,假設(shè)地面為理想導(dǎo)體,障礙物近似如圖1所示.

圖1 障礙物(矩形)模型Fig.1.The model of obstacle(rectangle).

圖1中ht,hr,h分別代表源的發(fā)射高度、觀察點(diǎn)的接收高度和障礙物的高度.在此基礎(chǔ)上,本節(jié)主要就雙向拋物線方法計(jì)算中的波源設(shè)置、基本原理以及對計(jì)算結(jié)果的相位修正三個方面進(jìn)行分析說明.

2.1 初始場設(shè)置的等效源模型

常用的拋物線方法的波源設(shè)置是利用遠(yuǎn)場與近場之間的關(guān)系反推初始場表達(dá)式,這種初始場u(0,z)的設(shè)置方法存在波束較窄或天線架高不能滿足u(0,z)=0(z<0)時(shí),初始場會延伸到地面以下的問題,并且這種方法沒有建立起初始場與天線電流分布之間的直接聯(lián)系.結(jié)合上述情況,本文使用等效源模型來逼近初始場設(shè)置[19,20].

已知拋物線方程由初始場推導(dǎo)出的自由空間遠(yuǎn)區(qū)場u(x,z)的解析形式如下[21]:

其中H(2)1表示一階第二類漢克爾函數(shù),ρ為源點(diǎn)到場點(diǎn)的距離,k為自由空間的波數(shù).現(xiàn)假設(shè)有一個在y方向無限長的有具體分布數(shù)值橫向磁流帶,電流方向?yàn)閦方向,

則該電流產(chǎn)生的電場為

由(3)和(1)式,可以看出兩者具有相似的形式.根據(jù)對偶定理,同理可推導(dǎo)出等效電流產(chǎn)生的磁場表達(dá)式,由此推出初始場與等效源的關(guān)系如下:

在(4)式的基礎(chǔ)上,可以通過等效源直接逼近拋物線方法的初始場設(shè)置,這種初始場設(shè)置建立了其與天線電流分布之間的直接聯(lián)系,不僅解決了之前方法存在的問題,而且為后面利用矩量法分析驗(yàn)證拋物線方程的計(jì)算精度提供了統(tǒng)一的入射源環(huán)境.

2.2 雙向拋物線方法的基本原理

在球面周向均勻起伏的地形下,通常使用雙向拋物線方法來求電波傳播問題.已知根據(jù)平坦地球模型,可得電磁場滿足二維標(biāo)量波動方程:

其中U表示電磁場量Ey,Hy;m代表修正后的大氣折射率.為解調(diào)(5)式中的相位快變部分,在波動方程的基礎(chǔ)上引入軸向衰減因子如下[1?3]:

其中UF,UB代表前向場和后向場.將(6)式分別代入(5)式,并使用Feit-Fleck寬角近似處理,得到雙向波動方程為

當(dāng)?shù)孛鏋槔硐雽?dǎo)體時(shí),通過分步傅里葉方法求得雙向拋物線方法的步進(jìn)公式如下:

其中F和F?1分別代表傅里葉變換和逆變換.可以看出在忽略傳播方向的前提下,推導(dǎo)出的步進(jìn)公式是完全一致的.為進(jìn)一步分析雙向拋物方程在不規(guī)則地形下的具體實(shí)現(xiàn)方法,以圖1中的障礙物模型為例,文獻(xiàn)[1]中雙向拋物線方法的計(jì)算方式如圖2所示.

圖2 單矩形障礙物下空間場值的處理 (a)前向傳播場;(b)后向傳播場Fig.2.The processing of spatial fi eld under single rectangular obstacle:(a)The forward fi eld;(b)the backward fi eld.

以圖2中對單矩形障礙物的處理方式,可得空間中場值設(shè)定如下:

由(9)式可以看出,當(dāng)?shù)?jì)算到障礙物時(shí),障礙物之下的場值被強(qiáng)行置零,只計(jì)算上空間的場值.文獻(xiàn)[1]中的這種計(jì)算方式將障礙物視為起伏地形,沒有考慮障礙物內(nèi)透射場的存在,也沒有考慮不同極化下障礙物邊界處場值的處理.

2.3 雙向拋物線方法的相位補(bǔ)償

在(9)式中場值分類的基礎(chǔ)上,雙向拋物線方法求解空間中總場如(6)式所示,是將前向和后向波統(tǒng)一到一個坐標(biāo)原點(diǎn)得出的,在疊加計(jì)算時(shí)忽略了前向場與后向場相位不一致的問題.從圖2(b)中對后向傳播場的處理中可以看出后向波源是在x1處產(chǎn)生的,根據(jù)障礙物反射面x1處的邊界條件可得

其中Γ為自由空間向障礙物垂直入射下的反射系數(shù).在迭代求解空間中任意一點(diǎn)x處的總場時(shí),可以看出前向波是從原點(diǎn)出發(fā),后向波是從x1出發(fā),此時(shí)x點(diǎn)的場強(qiáng)應(yīng)該為

對比(11)式與(6)式得,在求解總場時(shí),需要先進(jìn)行相位補(bǔ)償將前向波與后向波的相位統(tǒng)一后疊加.將(11)式代入(6)式得

通過(12)式可以修正雙向拋物線方法總場計(jì)算時(shí)的相位差異,但是在此之前的計(jì)算處理都是將障礙物處理為起伏地形,沒有從根本上處理存在障礙物時(shí)空間中的場值該如何計(jì)算的問題.因此本文提出一種新的拋物線方程計(jì)算方法,將障礙物處場值的計(jì)算進(jìn)行單獨(dú)處理.

3 存在障礙物時(shí)拋物線方程的處理方法

當(dāng)?shù)孛娲嬖谡系K物時(shí),經(jīng)典的拋物線方法是將障礙物與地面視為整體,將起伏地形之下的場值全部置零,并根據(jù)地面媒質(zhì)參數(shù)的不同,使用不同的分步傅里葉變換法.這種計(jì)算方法存在以下兩方面的問題:第一方面是因?yàn)樵摲椒ǖ挠?jì)算精度受限于障礙物的傾角大小;第二方面是沒有考慮障礙物內(nèi)部場值對空間中場值分布的影響.本文提出了一種區(qū)域分解算法.

3.1 區(qū)域分解算法

以圖2中對單矩形障礙物的處理為例,經(jīng)典拋物線方法的計(jì)算區(qū)域是從地面一直到吸收層.本文在此基礎(chǔ)上,對障礙物區(qū)域場值的求解進(jìn)行了分區(qū)處理,當(dāng)步進(jìn)求解到障礙物處,場值的計(jì)算分為兩部分,第一部分是障礙物之上的場值計(jì)算,第二部分是障礙物之下的場值計(jì)算.根據(jù)區(qū)域的不同,采用不同的離散傅里葉變換方法,并對障礙物上邊界處的場值計(jì)算加以修正.根據(jù)區(qū)域分解原理,對場值計(jì)算處理如下.

由圖3可以看出,障礙物區(qū)域的場值計(jì)算主要分為三部分.

1)障礙物上半空間的場值計(jì)算.由圖3可以看出,計(jì)算障礙物上半空間的場值時(shí),該區(qū)域中的下邊界從障礙物的上邊界B點(diǎn)往障礙物之下延伸了一段,這種延伸方式是為了避免對區(qū)域突然截?cái)鄷r(shí)造成的計(jì)算誤差.對重疊區(qū)延伸長度的選擇與計(jì)算精度有關(guān),通常為了計(jì)算方便,直接取障礙物高度的一半.在此基礎(chǔ)上,此時(shí)上空間的下邊界不再是理想導(dǎo)體而是均勻介質(zhì),所以上半空間的迭代計(jì)算相當(dāng)于阻抗邊界下的混合傅里葉求解.對下邊界的處理等效為Leontovich阻抗邊界條件[21]

根據(jù)(13)式,通過離散混合傅里葉變換法就可以求解障礙物之上包括重疊區(qū)域的場值.

圖3 障礙物處場值計(jì)算的區(qū)域分解Fig.3.Domain decomposition of fi eld calculation in obstacle zone.

2)障礙物內(nèi)部空間的場值計(jì)算.此時(shí)介質(zhì)障礙物內(nèi)場值計(jì)算的上邊界也發(fā)生了改變,為了提高計(jì)算精度,直接將障礙物內(nèi)部的場值計(jì)算擴(kuò)充到整個空間中,因?yàn)榈孛嫒允菍?dǎo)體,所以下邊界的處理使用導(dǎo)體邊界條件,傅里葉變換采用正弦變換計(jì)算,最后出障礙物恢復(fù)空間場值時(shí)只提取障礙物高度h以下區(qū)域的場值.

3)求解完障礙物上下空間的場值后,如圖3所示,在求解障礙物u(x±?x,z)時(shí),需要對介質(zhì)障礙物上邊界C處的場值進(jìn)行修正.該點(diǎn)處場值計(jì)算如下:

(14)式是將連接的邊界點(diǎn)處的場值使用兩種區(qū)域算法的平均值代替,這樣確保邊界切向場連續(xù).

綜上,在整個空間的迭代計(jì)算中遇到障礙物處就將計(jì)算區(qū)域按圖3所示分成兩部分計(jì)算,其余平坦區(qū)域的場值計(jì)算方法不變.

3.2 介質(zhì)中拋物線方程的處理及相位補(bǔ)償

在圖3的基礎(chǔ)上,障礙物之下的場值計(jì)算等效于介質(zhì)體內(nèi)的電波傳播問題的求解.因?yàn)榈孛嫒詾槔硐雽?dǎo)體,所以在障礙物內(nèi)的迭代計(jì)算與自由空間的迭代計(jì)算相比區(qū)別在于空間中的折射率發(fā)生了改變,等效為(6)式的解調(diào)因子的改變:

在(16)式步進(jìn)求解的基礎(chǔ)上,還需要考慮進(jìn)出障礙物時(shí)電磁波的透射作用,場值修正為

當(dāng)障礙物為理想導(dǎo)體時(shí),因?yàn)檎系K物內(nèi)不存在場值,所以只需要考慮上半空間的場值求解.又因?yàn)榇藭r(shí)上半空間的下邊界為導(dǎo)體,所以根據(jù)極化模型的不同,混合傅里葉變換法也簡化為正弦或余弦傅里葉變換法,并在恢復(fù)場值時(shí),直接將下半空間的場值賦零.

4 雙向拋物線方程矩量法驗(yàn)證分析

為了驗(yàn)證所提算法的精度,本文提出在近距離下用矩量法精確求解空間中的場強(qiáng)分布,從而驗(yàn)證存在障礙物的情況下新算法計(jì)算結(jié)果的精確性.因?yàn)檎系K物為理想導(dǎo)體,也等效為一種均勻介質(zhì),所以下面主要考慮矩形障礙物為均勻介質(zhì)的情況,并分別從兩種極化角度出發(fā)分析討論矩量法的處理過程.

如圖1所示的障礙物,將地面影響等效到格林函數(shù)G中,只對障礙物采用矩量法分析進(jìn)行處理.由于障礙物為均勻介質(zhì),所以可以用Poggio-Miller-Chang-Harrington-Wu(PMCHW)方程和體等效兩種方法來求解空間中的場強(qiáng)分布,本文選擇PMCHW方程組來聯(lián)立求解上述問題:

其中Hi,Ei,J,M分別代表入射磁場、入射電場以及障礙物表面的等效電流與等效磁流;z1,z2分別指自由空間和均勻介質(zhì)障礙物內(nèi)的波阻抗.L(X),K(X)都是簡化表達(dá)的算子,定義如下:

聯(lián)立(19)式中的算子,通過(18)式中等效電磁流與入射場之間的關(guān)系可以實(shí)現(xiàn)任意極化下空間中的場值求解.因?yàn)榇嬖趯?dǎo)體地面,考慮地面鏡像作用,對于電磁流源,(19)式中的二維格林函數(shù)為

其中ρi,ρ′i代表源點(diǎn)及其鏡像點(diǎn)的位置矢量,G0表示零階第二類漢克爾函數(shù),鏡像格林函數(shù)的正負(fù)取決于等效電磁流的鏡像原理.因?yàn)樵谇蠼膺^程中存在二階導(dǎo)數(shù)的計(jì)算,所以對于縱向電磁流采用點(diǎn)匹配脈沖基求解,橫向電磁流基函數(shù)和檢驗(yàn)源均采用三角基函數(shù)矩量法求解.

在TE波(垂直極化)照射下,障礙物在空間中產(chǎn)生的散射磁場只有y分量.因此使用二維等效電流來逼近拋物線方法的初始場,在垂直極化下,障礙物表面會產(chǎn)生環(huán)向的等效電流,流向同圖1,障礙物表面的等效磁流為y向,進(jìn)而求得空間中總場如下:

在TM波(水平極化)照射下,空間中只有y向電場,因此使用二維的等效磁流來等效初始場,在水平極化下,障礙物表面的等效電流為y向,等效磁流變成環(huán)向,積分方程中的各項(xiàng)分量的求解可以通過對偶原理從(21)式直接求得,惟一的區(qū)別在于考慮鏡像作用時(shí),鏡像電流的等效與鏡像磁流的等效相反.同理求得空間中總電場如下:

5 計(jì)算實(shí)例和結(jié)果分析

本文分析的問題基于圖1所示的矩形障礙物模型,用矩量法計(jì)算近距離下電波傳播的具體場值,從而驗(yàn)證在存在障礙物情況下新算法計(jì)算結(jié)果的精確性,并在障礙物為理想導(dǎo)體的情況下將新算法與文獻(xiàn)[2]的PEtool拋物線方法進(jìn)行對比分析.在下述實(shí)例分析中,矩量法的剖分間隔為0.5 m,新算法的步進(jìn)間隔是1 m,并且所有實(shí)例都是在i5處理器、32位操作系統(tǒng)、4G內(nèi)存的計(jì)算機(jī)設(shè)備環(huán)境下運(yùn)行仿真的.

實(shí)例一:發(fā)射源工作頻率30 MHz,高度為200 m,其發(fā)射波束仰角0?,波束寬度20?,水平極化,輻射源等效輻射功率為W,Z0為自由空間波阻抗.障礙物為雙矩形理想導(dǎo)體,矩形模型形狀一致,第一個矩形模型位于距離發(fā)射源750 m處,第二個矩形模型位于距離發(fā)射源900 m處,障礙物模型寬50 m,高度100 m.計(jì)算得到在接收高度為80 m處場強(qiáng)隨水平距離變化的場強(qiáng)振幅分布如圖4(a)所示,圖4(b)為圖4(a)的局部展開結(jié)果,其中MoM代表矩量法計(jì)算結(jié)果,PE2代表本文提出的存在障礙物情況下基于雙向拋物線方程的新算法的計(jì)算結(jié)果,PEtool代表文獻(xiàn)[2]算法求得的空間場值.

分析圖4(a)可以看出:障礙物后新算法的計(jì)算結(jié)果與矩量法相比幅度基本一致,障礙物之前的場值振蕩較為劇烈,造成這種現(xiàn)象的原因是地面反射波與直達(dá)波疊加相互作用產(chǎn)生的;障礙物之間的場值振蕩是因?yàn)檎系K物之間存在反射場的多次疊加.障礙物之后場值的抬高是因?yàn)殡S距離增加,導(dǎo)體障礙物的繞射場從深陰影區(qū)逐步過渡到直達(dá)波能夠照射到的亮區(qū).由圖4(b)可以看出PEtool算法的計(jì)算結(jié)果與新算法和矩量法相比起伏的相位出現(xiàn)了偏移,幅度也不完全一致.從而說明邊界處理和相位修正對雙向拋物線方法的計(jì)算精度有較大的改善,驗(yàn)證了新算法在障礙物為理想導(dǎo)體情況下的精確性.

圖4 (網(wǎng)刊彩色)水平極化雙矩形模型場值對比 (a)整體場值;(b)局部場值Fig.4.(color online)The contrasts of fi elds for horizontal polarization under double rectangular models:(a)The whole fi elds;(b)the local fi elds.

實(shí)例二:驗(yàn)證水平極化下圖1中單矩形介質(zhì)障礙物模型下新方法計(jì)算結(jié)果的精確性,介質(zhì)的參數(shù)為εr=4,σ=10?5s/m,將發(fā)射源的波束寬度更改為30?,單矩形介質(zhì)障礙物距離波源750 m,障礙物寬度變?yōu)?0 m,其余參數(shù)與實(shí)例一相同,計(jì)算的接收高度在80 m處場強(qiáng)隨水平距離變化的場強(qiáng)振幅分布如圖5所示,圖5(b)為圖5(a)的局部展開結(jié)果.

分析圖5可以看出,當(dāng)障礙物為均勻介質(zhì)時(shí),新算法的計(jì)算結(jié)果與矩量法相比幅度基本一致,尤其是障礙物之后800 m以上,在波形變換凹口處也完全一致,因此可以得出當(dāng)存在障礙物時(shí),使用分區(qū)原理以及介質(zhì)內(nèi)的相位分布修正方法可以有效地計(jì)算此時(shí)空間中的場值大小,也驗(yàn)證了新算法在障礙物為均勻介質(zhì)情況下的精確性.

在該實(shí)例中,矩量法的計(jì)算時(shí)間為335 s,新算法的計(jì)算時(shí)間為5 s,在計(jì)算精度幾乎一樣的情況下,新算法的計(jì)算時(shí)間比矩量法快了將近67倍,從而驗(yàn)證了新算法計(jì)算的有效性.

圖5 (網(wǎng)刊彩色)單矩形介質(zhì)障礙物場值對比 (a)整體場值;(b)局部場值Fig.5.(color online)The contrasts of fi elds under single rectangular dielectric obstacle:(a)The whole if elds;(b)the local fi elds.

實(shí)例三:在實(shí)例一參數(shù)的基礎(chǔ)上將發(fā)射源的波束寬度改為30?,雙矩形寬度改為20 m,其余參數(shù)不變,雙矩形均勻介質(zhì)障礙物的介電常數(shù)同實(shí)例二,計(jì)算水平極化下接收高度為80 m處場強(qiáng)隨水平距離的變化如圖6所示,圖6(b)為圖6(a)的局部展開結(jié)果.

從圖6可以看出,在雙矩形均勻介質(zhì)障礙物模型下,新算法與矩量法的計(jì)算結(jié)果相比,場值變換趨勢基本一致.由圖6(b)中雙矩形障礙物之間的場值比較可以看出,新算法與矩量法的計(jì)算結(jié)果相比誤差非常小,從而體現(xiàn)了新算法計(jì)算的精確性,此時(shí)矩量法所需的時(shí)間為1097 s,而新算法的計(jì)算時(shí)間僅為65 s.當(dāng)?shù)匦巫儞Q更為復(fù)雜時(shí),矩量法的計(jì)算難度會迅速增大甚至無法計(jì)算,新算法的優(yōu)勢性會更進(jìn)一步體現(xiàn).在此基礎(chǔ)上,對于雙向情況下新算法對障礙物內(nèi)的計(jì)算誤差有待進(jìn)一步研究.

圖6 (網(wǎng)刊彩色)雙矩形介質(zhì)障礙物場值的對比 (a)整體場值;(b)局部場值Fig.6.(color online)The contrasts of fi elds under double rectangular dielectric obstacles:(a)The whole if elds;(b)the local fi elds.

6 結(jié) 論

經(jīng)典的拋物線方法只能用于求解起伏地形下的電波傳播問題的計(jì)算,當(dāng)?shù)孛娲嬖谡系K物時(shí),尤其是兩種媒質(zhì)介質(zhì)參數(shù)不同的情況下,雙向拋物線方法無法計(jì)算.本文基于雙向拋物線方法的基本原理,提出了一種新的區(qū)域分解算法,對存在障礙物時(shí)空間中的場值計(jì)算進(jìn)行了區(qū)域分解處理,并對均勻介質(zhì)內(nèi)雙向拋物線方法的求解給出了計(jì)算公式.通過上述方法,使得基于拋物線方程的新算法可以處理地面存在障礙物的情況.在此基礎(chǔ)上,對求解總場時(shí)的前向場與后向場的相位以及出介質(zhì)障礙物邊緣時(shí)上下空間的場值疊加時(shí)的相位都分別進(jìn)行了修正,修正方法有效改善了新算法計(jì)算的空間中場值的相位匹配性.最后,結(jié)合矩量法在近距離下對算法在存在障礙物環(huán)境下的電磁傳播計(jì)算結(jié)果進(jìn)行了驗(yàn)證分析,通過實(shí)例分析,驗(yàn)證了基分區(qū)原理的雙向拋物線方程新算法計(jì)算出場值的精確性以及優(yōu)越性.

[1]Ozlem O 2009 IEEE Trans.Antenn.Propag.57 2706

[2]Ozlem O,Gokhan A,Mustafa K,Levent S 2011 Comput.Phys.Commun.182 2638

[3]Wang K,Long Y L 2012 IEEE Trans.Antenn.Propag.60 4467

[4]Zhang P,Bai L,Wu Z S,Guo L X 2016 IEEE Trans.Antenn.Propag.Mag.58 31

[5]Wang D D,Xi X L,Pu Y R,Liu J F,Zhou L L 2016 IEEE Trans.Antenn.Wireless Propag.Lett.15 734

[6]Yuan X J,Lin W G 1993 Chin.Phys.Lett.10 57

[7]Omaki N,Yun Z Q,Iskander M F 2012 2012 IEEE International Conference on Wireless Information Technology and Systems(ICWITS)Maui,USA,November 11–16,2012 p1

[8]Kuttler J R 1999 IEEE Trans.Antenn.Propag.47 1131

[9]Donohue D J,Kuttler J R 2000 IEEE Trans.Antenn.Propag.48 260

[10]Beilis A,Tappert F D 1979 J.Acoust.Soc.Am.66 811

[11]Wang Y J,Guo L X,Li Q L 2016 11th International Symposium on Antennas,Propagation and EM Theory(ISAPE)Guilin,China,October 18–21,2016 p404

[12]Ozgun O,Sevgi L 2012 Aces J.27 376

[13]Gokhan A,Levent S 2013 IEEE Trans.Antenn.Propag.Mag.55 244

[14]Zhu J,Yin C Y,Wei Q F 2016 J.Microwaves 32 32(in Chinese)[祝杰,尹成友,魏喬菲 2016微波學(xué)報(bào) 32 32]

[15]Omak N,Yun Z Q,Iskander M F 2012 Antennas and Propagation Society International Symposium(APSURSI)Chicago,USA,July 8–14,2012,p1

[16]Lu J,Zhou H C 2016 Chin.Phys.B 25 90203

[17]Pvel V,Pvel P 2007 IEEE Antenn.Wireless Propag.Lett.6 152

[18]Sheng X Q 2004 Computational Electromagnetic Theory(Beijing:Science Press)pp49–53(in Chinese)[盛新慶2004計(jì)算電磁學(xué)要論(北京:科學(xué)出版社)第49—53頁]

[19]Zhu J,Yin C Y 2016 J.Microwaves 32 26(in Chinese)[祝杰,尹成友 2016微波學(xué)報(bào)32 26]

[20]Yin C Y,Zhu J,Wei Q F 2016 37th Progress in Electromagnetics Shanghai,China,August 8–11,2016 p1655

[21]Levy M 2000 Parabolic Equation Methods for Electromagnetic Wave Propagation(London:IEE Press)pp149,287–291

PACS:41.20.JbDOI:10.7498/aps.66.124102

Research and veri fi cation for parabolic equation method of radio wave propagation in obstacle environment?

Wei Qiao-FeiYin Cheng-You?Fan Qi-Meng

(National Key Laboratory of Pulsed Power Laser Technology,Electronic Engineering Institute of PLA,Hefei 230037,China)

14 January 2017;revised manuscript

20 March 2017)

In recent years,the two-way parabolic equation method(2WPE)has been widely utilized for studying the tropospheric ground-wave propagation under the irregular terrain.This algorithm can deal with the in fl uences of the irregular terrain characteristic and the di ff erent electromagnetic parameters of the surface structure on wave propagation.However,there are still some defects in 2WPE method.Firstly,the method considers the irregular terrain and obstacles as a whole,so it cannot deal with the situation where the medium parameters of obstacles and the ground are di ff erent.Secondly,its calculation precision is limited with the inclination of the undulating terrain:if there are obstacles the upper bound of the inclination is easily broken through.Therefore,in this paper,a novel two-way parabolic equation method is proposed for analyzing the radio wave propagation in obstacle environment.According to the principle of domain decomposition,the obstacle zones are divided into two domains in the new algorithm,and the two subdomains are calculated,respectively.Meanwhile,in order to avoid the calculation error caused by the abrupt truncation of the obstacle zone,the fi eld at the upper boundary of obstacles is modi fi ed to ensure the continuity of tangential fi eld.To further improve the accuracy of the new algorithm,according to the historical transmission paths,we exactly retrieve the phases of each forward and backward wave,especially when stepping in and out of the obstacles.Furthermore,the method of moment(MoM)is used to verify the calculation accuracy of the new algorithm in obstacle environment.Although the accuracy of the MoM is very high,it also requires a great deal of calculation resources:it can only be employed to compute the fi elds in short distance.To overcome the difficulty,we use the image principle in the obstacle environment and do not subdivide the ground into segments;therefore the veri fi cation accuracy can be improved.On this basis,to unify the source setting of the new algorithm and the MoM,the equivalent source model is used to set the initial fi eld.Finally,through numerical experiments,the simulation results of both methods agree very well,so the e ff ectiveness of the boundary correction and the phase correction which are presented in this paper are both veri fi ed.The accuracy and superiority of the new algorithm in obstacle environment are also demonstrated.To sum up,the novel two-way parabolic equation method can be used to accurately calculate the fi eld of the space in the obstacle environment,and lays the foundation for the fi eld calculation of radio wave propagation in real environment.

parabolic equation,radio wave propagation,domain decomposition,method of moment

10.7498/aps.66.124102

?總裝備部預(yù)研基金(批準(zhǔn)號:51333020201)資助的課題.

?通信作者.E-mail:cyouyin@sina.com

?2017中國物理學(xué)會Chinese Physical Society

http://wulixb.iphy.ac.cn

*Project supported by the General Equipment Department Pre-Research Foundation,China(Grant No.51333020201).

?Corresponding author.E-mail:cyouyin@sina.com