張志禹，孫戈，徐高晨，許耀斌

　?。ㄎ靼怖砉ご髮W 自動化與信息工程學院陜西 西安 710048）

       摘要：傳統拋物線方程方法求解起伏地形路徑下的電波傳播無法反映反射和散射效應，對刃峰地形、類似高斯型地形等起伏地形，采用分步傅里葉變換方法求解雙向拋物線方程方法并考慮前向和后向傳播的疊加總場，將計算結果與傳統拋物線方程方法進行比較，預測結果更準確。

　　關鍵詞：電波傳播；雙向拋物線方程；分步傅里葉變換方法

　　中圖分類號：TN011文獻標識碼：ADOI： 10.19358/j.issn.1674-7720.2016.24.025

　　引用格式：張志禹，孫戈，徐高晨，等. 基于寬角雙向拋物線方程的高頻電波傳播預測方法［J］.微型機與應用，2016,35（24）：87-90.

0引言

　　拋物線方程（Parabolic Equation, PE）模型是在研究電波在光滑地表繞射傳播特性時提出的［1］，其理論依據是在一定條件下對二階橢圓型的波動方程進行因式分解，取其傳播方向上一階導數的拋物型方程。目前PE模型求解方法主要分為三類：DMT法［2］、FD法［3］、有限元法［4］。由于障礙物的遮擋是影響電波傳播的重要因素，對于前向傳播的PE模型，并沒有考慮電波傳播遇到障礙物時的反射效應，會給預測帶來一定的偏差。文獻［5］利用雙向PE模型反演對障礙物進行定位，取得了很好的效果。本文利用雙向PE模型分別預測單個起伏地形和多個起伏地形環(huán)境下的高頻電波傳播特性，將數值計算結果與單向PE模型比較，驗證了雙向PE模型具有更準確的預測結果。

1寬角雙向拋物線方程模型

　　1.1雙向拋物線方程模型

　　雙向拋物線方程模型在傳統PE模型基礎上考慮后向傳播，將后向傳播場按距離步進疊加到前向傳播場，從而獲得整個觀測區(qū)域的電波傳播場分布。當電波傳播遇到起伏地形時，使用地形平移模型。地形平移法是對地形屏蔽法的改進，該方法不需要近似處理地形，只需要計算每一步迭代時的地形高度差，再根據地形高度差確定垂直方向平移的計算區(qū)域場值點數，平移方向根據地形斜率判斷。

　　在直角坐標系(x,y,z)下，x表示距離，z表示高度。假設ψ與y無關，則從Maxwell方程組推導的波動方程滿足下式：

　　

　　式（2）中第一式表示前向傳播即為傳統PE模型，第二式表示后向傳播。由于Q計算比較復雜，FeitFleck近似法可將偽微分算子Q展開為：

　　式中p=ksinθ，表示z的Fourier變換域，θ為電波到水平方向的角度。

　　對偽微分算子Q用FeitFleck近似法，考慮前向和后向傳播，則有：

　　

　　式中uf、ub分別表示前向、后向傳播的空間域，Uf、Ub分別表示前向、后向傳播的P域場值，ub反映傳播電波遇到障礙物時的反射效應。

　　1.2單雙刃峰原理

　　以單刃峰環(huán)境為例，單刃峰位于發(fā)射點水平距離xobs處，高為hobs，2WPE算法示意圖如圖1所示，電波從初始場x=0處沿x軸向右傳播，在傳播中遇到障礙物時發(fā)生反射效應，因此電波達到刃峰之前只考慮前向傳播，到達刃峰時分解為前向傳播和后向傳播兩個分量，兩者分別沿x軸的正方向和負方向。通過邊界條件得到刃峰處后向傳播的初始場如式（9）：

　　圖1單刃峰2WPE示意圖

　　整個計算區(qū)域的總場可以表示為：

　　由式(10)看出：電波傳播在單刃峰路徑上時，電波在刃峰處被分解，具有電磁波的繞射和反射效應。當前向傳播到達計算區(qū)域最大處時，后向傳播到達天線發(fā)射處，整個計算才結束。

　　當障礙物為多刃峰時，由于反射作用，在兩個刃峰之間會出現無窮多項前向傳播和后向傳播，因此計算時，要進行反復迭代，圖2為雙刃峰2WPE算法示意圖。

2數值算例

　　下面具體分析單個起伏地形和多個起伏地形環(huán)境下，2WPE模型預測電波傳播特性，并將其結果與單向PE模型比較。在以下算例中假設頻率都為1 000 MHz，電磁波水平極化，天線方向圖為高斯型，波束寬度為3°，水平地面和起伏地形均為理想導體。

　　2.1類似高斯型地形對電波傳播的影響

　　地形函數為：

　　其中x0為包絡半高程到天線發(fā)射點的距離，w為包絡半寬度，h為高斯包絡的垂直高度。

　　取x0=24 km，w=1.25 km，h=500 m。

　　圖3為傳統PE模型和2WPE模型的傳播因子分布剖面圖，傳播因子單位按分貝（dB）計算。從圖中可以看出，傳統的PE模型只能計算前向傳播，而2WPE模型能同時計算前向和后向傳播。觀測點位于障礙物左側區(qū)域時，傳統PE模型結果偏差較大；觀測點在障礙物右側區(qū)域時，兩種模型計算的結果一致。

　　下面具體分析兩種模型的傳播特性，圖4給出傳播因子隨水平距離的變化曲線，兩個圖中觀察點的垂直高度分別為600 m和400 m。綜合分析圖3(a)和圖3(b)，高斯地形左側0~20 km區(qū)域內，傳統PE模型由于本身算法限制而無法計算后向傳播，因此傳播因子非常小，而2WPE模型疊加了反射波，預測結果更加精確。兩個圖中的障礙物左側區(qū)域，2WPE模型得到的傳播因子曲線出現波動，這是前向傳播與后向傳播疊加的結果；在x=25 km后，兩種模型都只有前向傳播，因此曲線基本重合，并且障礙物的遮擋使得障礙物右側的傳播因子急速減小，在很長一段距離內難以恢復。圖4(b) 的觀測高度小于地形高度，因此傳播因子曲線出現分段現象，這與實際情況相吻合。

　　圖5是傳播因子與垂直高度的關系曲線圖。從圖5(a)中看出，在0~100 m范圍內，遇到障礙物后的后向傳播能量集中在100 m高度以上區(qū)域；高度300 m以上時，2WPE模型的傳播因子值大于傳統PE模型計算的值，差值主要是由2WPE模型后向傳播引起的。圖5(b)所示的觀察點x=30 km位于起伏地形右側，不存在后向傳播，兩種模型的計算方法相同，因此曲線吻合。從圖中可以看出，傳播因子在高度為600 m左右處，向上和向下的傳播都在減小。通過設置多個激勵源來擴大傳播區(qū)域可解決此問題。

　　2.2多個起伏地形對電波傳播的影響

　　算例：發(fā)射天線高度為HT=50 m，地表為PEC，水平最大距離Xmax=100 km，最大高度Zmax=500 m，三個單刃峰的水平位置分別為30 km、50 km、70 km，對應的高度分別為100 m、200 m、300 m。

　　圖6為傳播因子在計算區(qū)域的剖面分布圖，傳播路徑上有多個單刃峰。由圖可以看出，其現象與單個刃峰情形類似。比較圖6(a)和圖6(b)發(fā)現，兩圖中左邊兩個刃峰處的傳播因子不同。圖6(a)能夠清晰看到三個刃峰右側附近區(qū)域電波被遮擋，而圖6(b)中左側兩個刃峰處無此現象，這是因為2WPE模型在兩個刃峰之間出現反復的后向傳播，兩個圖中距發(fā)射點70 km處的刃峰右側無后向傳播，因此都出現電磁盲區(qū)的現象。

　　圖7為傳播因子與水平距離和豎直高度的關系曲線圖，圖7(a)觀測點在垂直高度HR=100 m處，從圖中很容易看出傳播因子的數值在刃峰處迅速減小。在水平方向距原點0~70 km范圍內，傳統PE模型只計算前向傳播，而2WPE模型同時考慮了后向傳播，因此2WPE模型得到的傳播因子曲線更加準確；水平方向距原點70 km以后兩種模型都只有前向傳播，因此曲線吻合。此外，由于能量主要集中在前向傳播，后向傳播的影響相對較小，兩種模型預測的結果基本相吻合，在計算精度要求不高的情況下，傳統的單向PE模型能夠更快預測傳播特性而提高效率。圖7(b)為傳播因子與高度的關系曲線圖，其中，觀察點在離原點水平距離40 km處。從圖中看出，在300 m高度以下，傳統PE模型的傳播因子數值偏??；在300 m高度以上，沒有障礙物遮擋，2WPE模型無后向傳播，兩種模型的曲線吻合。

3結論

　　雙向拋物線方程模型克服了傳統拋物線方程模型不能計算后向傳播的缺點，本文用SSFT算法求解寬角雙向拋物線方程模型的推導和實現過程。采用地形平移法處理地形邊界，并且用該模型分析了單個類似高斯型起伏地形、多個刃峰地形等復雜地形對電波傳播的影響，將計算的數值結果與單向拋物線方程模型的結果進行比較，得到更準確的電波傳播預測結果。

參考文獻

　　［1］ LEONTOVICH M A, FOCK V A. Solution of propagation of electromagnetic waves along the earth’s surface by the method of parabolic equations［J］. Journal of Physics, USSR, 1946, 10(1): 1323.

　?。?］ 吳凡. 大氣波導中的拋物型方程法研究［D］. 武漢：武漢理工大學，2008.

　?。?］ 袁曉君，林為干. 電波傳播的拋物方程近似分析法［J］. 電子科技大學學報，1993,22(1)：4348.

　?。?］ 郭建炎，王劍瑩，龍云亮. 基于拋物方程法的粗糙海面電波傳播分析［J］. 通信學報，2009,30(6)：4752.

　　［5］ 王昆, 龍云亮, 劉震宇. 基于雙向拋物方程逆算法的障礙物定位技術研究［J］.電波科學學報，2015, 30(6):11081115.