黃 穎

(南京郵電大學(xué) 電磁場與微波技術(shù)系，江蘇 南京 210003)

傅里葉分布拋物方程法的準三維研究

黃 穎

(南京郵電大學(xué) 電磁場與微波技術(shù)系，江蘇 南京 210003)

文章介紹了二維拋物方程的傅里葉分步法，提出了一種基于二維方法的準三維研究方法。將三維圖形分解為兩個二維圖形。首先通過二維方法解決圓錐的主視圖(三角形)和俯視圖(圓形)的電波傳播問題，將主視圖和俯視圖的計算結(jié)果分別與參考文獻[5]和HFSS中的模型進行對比，從而驗證其正確性。該方法比二維的研究范圍更廣，同時比直接使用三維公式求解簡單。

電波傳播；拋物方程；傅里葉分步法；準三維

0 引言

拋物方程法(Parabolic Equation,PE)是近年來研究電波傳播常用的一種數(shù)值方法。通過對波動方程近似得到關(guān)于傳播方向一階導(dǎo)數(shù)的拋物型方程[1]。

目前，求解拋物方程的方法主要有傅里葉分步法 (Fourier Split-Step, SSFT)和有限差分分解法(Finite Difference, FD)。SSFT算法和FD算法均為步進迭代計算的方法，但各自有不同的特點。SSFT算法對步進的限制非常寬松，允許相對較大的水平步長，因而計算速度很快，且數(shù)值穩(wěn)定性高。雖然FD算法對于復(fù)雜的地表邊界條件的處理比較簡單，但步進會受到電波波長的限制，網(wǎng)格劃分較細，計算速度慢。且FD算法需要大量的矩陣運算，對于計算機硬件也有一定的要求。此外，PE方法是一個定解問題，也就是說，若給出已知的初始場條件和邊界條件，就能夠通過數(shù)值算法求得定解。因此在求解遠距離、大范圍的電波傳播問題時，拋物方程主要采用 SSFT 來進行求解[2]。

目前采用SSFT研究大型不規(guī)則地面的電波傳播問題都是基于二維。雖然SSFT算法計算速度很快，且數(shù)值穩(wěn)定性高，但二維模型無法反應(yīng)橫向媒質(zhì)產(chǎn)生的橫向散射等效應(yīng)，且三維模型不需要專門尋找三維地形對電波傳播產(chǎn)生反射或繞射的區(qū)域，也不需要設(shè)置各種傳播機制的判據(jù)，因此該方法計算量相對較小，求解精度更高[1-7]。近年來對三維拋物方程電波傳播問題的研究越來越受重視[8-9]，但三維模型的分析求解過程比二維復(fù)雜，本文提出一種準三維方法，即通過研究三維圖形的主視圖和俯視圖將三維問題視作兩個二維問題，分別對這兩個二維問題進行研究，最終將其傳播因子數(shù)值在相同高度處進行相加，即準三維。

本文首先介紹二維傅里葉分布拋物方程法的求解公式，然后提出用準三維方法研究圓錐形障礙物的電波傳播問題，從視圖的角度將圓錐問題分解為三角形和圓形。將三角形的傳播因子曲線與文獻[5]相比較，驗證其正確性；并通過在HFSS中建立模型驗證圓形障礙物電波傳播的正確性。該方法比二維的研究范圍更廣，同時比直接使用三維公式求解簡單。

1 理論分析

拋物方程是已知源點x=0處的波，求解與其相距不遠處的下一點波的大小。通過步進迭代法可以得到距離x=0處的下一點x+Δx的波動大小近似為：

(1)

其中，u是波幅度，x和z分別表示橫向和縱向(傳播方向)坐標；k0是自由空間中的波數(shù)；F代表快速傅里葉變換(FFT)，p表示轉(zhuǎn)換變量，p=ksinθ中θ是距離水平面的角度；m是由m=n2-1+2z/ae折射指數(shù)變換而來，其中n是折射指數(shù)，ae是地球半徑[9]。

圖1是利用文獻[5]的參數(shù)，將三維模型轉(zhuǎn)化為二維模型的過程，其中(a)是目標物體圓錐，圓錐的底面是半徑為2.5 km的圓，高是50 m，源的高度為10 m；(b)和(c)分別是圓錐在(x,z)面(主視圖)和(x,y)面(俯視圖)的三角形和圓形。另外，圖1中俯視圖的高度設(shè)置為點源處，因此圖1(c)中的半徑為2 km。

對圖1(b)和(c)進行分析時，其計算區(qū)域不同，采用相同天線，假設(shè)電波在平面地表上傳播，其參數(shù)是：頻率為100 MHz、水平極化、垂直波數(shù)寬度為10°的高斯天線。迭代步長Δx=50 m，物體表面為理想導(dǎo)體表面，地面近似為平地面。

圖1 模型圖

2 三角形模型的分析

圖1中三角形模型的計算區(qū)域是30 km×500 m。計算中使用的天線參數(shù)如上節(jié)所述。對主視圖中三角形障礙物的電波傳播進行分析計算。

主視圖中的發(fā)射天線高度設(shè)為10 m。圖2所示是底邊長為5 km，高度為50 m，起點在坐標原點的三角形障礙物電波傳播的偽彩圖。該圖提供了足夠好的傾斜地形采樣。峰值衍射和陰影引起的場變化尖銳，特點正確明顯。

圖2 三角形障礙物電波傳播的偽彩圖

圖3為傳播距離為5 km處的傳播因子圖。

圖3 傳播距離為5 km處的傳播因子

從圖3的曲線可看出本文中的計算結(jié)果與文獻[5]的結(jié)果非常吻合。

3 圓形模型的分析

俯視圖中發(fā)射天線高度設(shè)為2.5 km，圓形模型的計算區(qū)域是5 km×5 km。其天線的參數(shù)與三角形模型的相同。用傅里葉分布拋物方程法得到在傳播距離為5 km處圓形障礙物的傳播因子。

為了驗證使用拋物方程法(PE)計算圓形障礙物電波傳播的正確性，本文利用了HFSS的仿真和計算。在HFSS中建立立方體輻射邊界，使用偶極子天線作為發(fā)射源，將圓柱視作障礙物，材質(zhì)為理想導(dǎo)體，圓柱的高度與立方體的邊長相等，立方體邊界設(shè)為輻射邊界。觀察與偶極子天線所在面相對的面的求解線上的傳播因子，來驗證拋物方程法的模型。偶極子天線距離圓柱中心的長度為L，圓柱底面半徑為R，為了與PE方法的尺寸相對應(yīng)，設(shè)定比例大小R/L=r/l=0.8，求解線距離圓柱中心的大小為L，通過HFSS仿真來計算其傳播因子。圖4為圓形障礙物使用拋物方程法得到的偽彩圖。

圖4 PE圓形障礙物的偽彩圖

從圖4可以看出，由于圓形障礙物的尺寸很大，電波在傳播過程中幾乎被障礙物阻擋。

圖5 100 MHz時PE和HFSS的傳播因子比較

圖5為采用PE和HFSS得到的傳播因子示意圖，從圖中可以看出PE和HFSS的結(jié)果比較相近，說明使用PE計算圓形障礙物電波傳播的正確性。

將比例大小R/L固定，把天線替換成頻率為300 MHz偶極子天線。圖6為300 MHz時PE方法和HFSS傳播因子。

圖6 300 MHz時PE和HFSS的傳播因子比較

從圖6看出PE和HFSS的曲線數(shù)值和趨勢幾乎吻合。根據(jù)繞射原理，當頻率越高時波長越短，繞射現(xiàn)象越不明顯，所以300 MHz頻率下的衰減比100 MHz時嚴重，因此，圖6的衰減值比圖5大。其比較結(jié)果與理論相符合，所以用傅里葉分布拋物方程法求解的圓形障礙物的傳播因子是合理的。

4 準三維的應(yīng)用

下面將三角形與圓形的結(jié)果相結(jié)合來研究準三維圓錐在頻率為100 MHz時的傳播因子。將圖1中的源每抬高1 m得到不同的俯視圖，分別計算不同高度處的半徑以及不同的圓形障礙物的傳播因子，最終將其在圓心處的傳播因子與三角形等同高度處相加，得到圖7所示的準圓錐在0～50 m處的傳播因子示意圖。

圖7 在0～50 m處的準圓錐的傳播因子示意圖

如圖7所示，隨著圓錐的高度增加，俯視圖中圓半徑逐漸減小，衰減值也隨之減小，高度10 m～45 m之間衰減值趨于平穩(wěn)，當半徑逐漸趨于0時，曲線急劇下降，衰減越少，障礙物對波的阻擋越小。

通過上面的結(jié)果可以得出本文提出的準三維的方法可以預(yù)測圓錐障礙物在電波傳播中的繞射現(xiàn)象，從而驗證了此方案的可行性。

在用拋物方程法計算三角形障礙物時的總步數(shù)是600步，時間為3.921 47 s；計算圓形障礙物時是100步，100 MHz的時間為2.916 79 s，300 MHz的時間為2.383 34 s。仿真代碼是使用MATLAB軟件運行，使用的是2.1 GHz-CPU的個人筆記本。

5 結(jié)論

本文提出了一種基于二維拋物方程法的準三維電波傳播方法，將三維圓錐分解為兩個二維的三角形和圓形，并將三角形和圓形的傳播因子的數(shù)值在相同高度處相加，得到圓錐的準三維傳播因子。通過對比文獻[5]和HFSS的結(jié)果驗證了兩個二維模型的正確性。該方法比二維的研究范圍更廣，同時比直接使用三維公式求解簡單。

[1] 李方, 察豪. 拋物方程法研究不規(guī)則地形對電波傳播的影響[J]. 火力與指揮控制, 2010, 35(8):114-116.

[2] 張青洪, 廖成, 盛楠,等. 拋物方程方法的亞網(wǎng)格模型及其應(yīng)用研究[J]. 電子與信息學(xué)報, 2014, 36(8):2005-2009.

[3] Zhi Chao, GUIZHEN L. Research of the application of the Fourier transform algorithm for parabolic equation in radio wave diffraction[C]. Environmental Electromagnetics, IEEE, 2015:94-96.

[4] CHEN Y, QIN W. A variable range step technique for propagation predictions over large irregular terrain using the Fourier split-step parabolic equation[C]. Asia-Pacific Microwave Conference, IEEE, 2015:1-3.

[5] 胡繪斌, 毛鈞杰, 柴舜連. 大氣折射率剖面在寬角拋物方程中的應(yīng)用[J]. 微波學(xué)報, 2006, 22(1):5-8.

[6] DONOHUE D J, KUTTLER J R. Propagation modeling over terrain using the parabolic wave equation[J]. IEEE Transactions on Antennas & Propagation, 2000, 48(2):260-277.

[7] HVIID J T, ANDERSEN J B, TOFTGARD J, et al. Terrain-based propagation model for rural area-an integral equation approach[J]. Antennas & Propagation IEEE Transactions on, 1995, 43(1):41-46.

[8] 倪軍林. 矢量拋物線方程在電波傳播中的應(yīng)用[D]. 杭州:杭州電子科技大學(xué), 2009.

[9] JANASWAMY R. Path loss predictions in the presence of buildings on flat terrain: a 3-D vector parabolic equation approach[J]. IEEE Transactions on Antennas & Propagation, 2003, 51(8):1716-1728.

Research of quasi three dimensional of Fourier split-step parabolic equation approach

Huang Ying

(Department of Electromagnetic Field and Microwave Technology, Nanjing University of Posts and Telecommunications, Nanjing 210003,China)

The two-dimensional (2-D) Fourier split-step parabolic equation approach is introduced, and the quasi-three-dimensional (quasi-3-D) wave propagation method based on 2-D is proposed in this paper. The 3-D patterning is divided into two 2-D patterning. The results of the main view (triangle) and the top view (circular) of the cone are calculated by 2-D approach, which are compared with the results of the reference [5] and the HFSS to verify their correctness respectively. The quasi-3-D is more extensive than the 2-D approach, and it is simpler than the direct use of the 3-D formula.

wave propagation; parabolic equation; Fourier split-step; quasi-three-dimensional

TN92

A

10.19358/j.issn.1674- 7720.2017.03.005

黃穎.傅里葉分布拋物方程法的準三維研究[J].微型機與應(yīng)用，2017,36(3)：16-18，22.

2016-10-12)

黃穎(1991-)，通信作者，女，碩士研究生，主要研究方向：拋物方程法在電波傳播預(yù)測計算中的應(yīng)用研究。E-mail：1143937804@qq.com。