蓋明明 溫 東 郗海龍 楚 涓

(1.海軍潛艇學(xué)院 青島 266199)(2.北海艦隊(duì)參謀部 青島 266000)

超低頻場強(qiáng)勒讓德函數(shù)的導(dǎo)數(shù)表達(dá)式對比分析*

蓋明明1溫 東1郗海龍2楚 涓1

(1.海軍潛艇學(xué)院 青島 266199)(2.北海艦隊(duì)參謀部 青島 266000)

超低頻電波大氣場強(qiáng)計(jì)算中或包含第一類勒讓德函數(shù)對距離的一階導(dǎo)數(shù)或包含其二階導(dǎo)數(shù),這兩個(gè)導(dǎo)數(shù)可以用無窮級數(shù)表示,但無窮級數(shù)計(jì)算受制于計(jì)算機(jī)的位數(shù)限制,因此在超低頻的高端無法用級數(shù)計(jì)算。除此之外,還有幾個(gè)近似表達(dá)公式。但近似公式的適用范圍和精度尚不明確。所以基于Matlab軟件對無窮級數(shù)公式和三個(gè)常見的近似公式進(jìn)行編程計(jì)算,相互比較,以確定各個(gè)公式的適用條件和精度差別。

超低頻; 大氣場強(qiáng); 勒讓德函數(shù)對距離的一階導(dǎo)數(shù); 勒讓德函數(shù)對距離的二階導(dǎo)數(shù)

Class Number E962

1 引言

超低頻發(fā)射天線目前主要是應(yīng)用低架于地面上的兩端接地的水平導(dǎo)線,作為這種天線的物理長度雖然很長(一般100km左右),但與波長相比仍很短,可以把它看作是水平電偶極子[1]。另外,地面上超低頻電磁波在地面-電離層波導(dǎo)內(nèi)發(fā)射和傳播的,因此,超低頻的輻射和傳播問題,從數(shù)學(xué)上可以歸納為在圓球形地面-電離層波導(dǎo)內(nèi)求解地面上水平電偶極子激發(fā)的電磁場問題。

在波導(dǎo)內(nèi),五個(gè)電磁場分量表達(dá)式中或包含第一類勒讓德函數(shù)Pν0(-cosθ)對地心角θ(距離)的一階導(dǎo)數(shù),或包含其二階導(dǎo)數(shù)。文獻(xiàn)[2～4]中給出勒讓德函數(shù)對距離的一階導(dǎo)數(shù)和二階導(dǎo)數(shù)的無窮級數(shù)公式的無窮級數(shù)形式為[2～4]

這兩個(gè)導(dǎo)數(shù)可以用無窮級數(shù)表示,這是比較精確的表示。但遺憾的是,無窮級數(shù)計(jì)算受制于計(jì)算機(jī)的位數(shù)限制,一般的個(gè)人計(jì)算機(jī)浮點(diǎn)運(yùn)算只提供小數(shù)點(diǎn)后15位精度,因此在超低頻的高端無法用級數(shù)計(jì)算[1]。除此之外,還有對極點(diǎn)波表達(dá)式、平面波導(dǎo)表達(dá)式、駐波波導(dǎo)表達(dá)式等幾個(gè)近似表達(dá)公式。但近似公式的適用范圍尚不明確。所以基于Matlab軟件對無窮級數(shù)公式和常見的三個(gè)近似公式進(jìn)行編程計(jì)算,相互比較,以確定各個(gè)公式的適用條件。

2 勒讓德函數(shù)對距離頭二階導(dǎo)數(shù)級數(shù)表達(dá)式與對極點(diǎn)波表達(dá)式對比分析

文獻(xiàn)[2～3]給出了勒讓德函數(shù)頭二階導(dǎo)數(shù)在對極點(diǎn)附近的近似表達(dá)式[2～3]

以白天為例,分別取頻率30Hz、100Hz在距離0km～20000km范圍進(jìn)行仿真,如圖1和圖2所示,其中實(shí)線‘-’代表級數(shù)公式計(jì)算結(jié)果,點(diǎn)線‘-*’代表駐波公式計(jì)算結(jié)果。

圖2 勒讓德函數(shù)對距離二階導(dǎo)數(shù)級數(shù)表達(dá)式與對極點(diǎn)波表達(dá)式對比分析

圖1和圖2上示出級數(shù)公式和以對極點(diǎn)為中心公式計(jì)算的絕對值曲線的比較。其中圖1(a)的頻率是30Hz,兩條曲線清晰地顯示,大約在3000km以后,它們符合極好,3000km以內(nèi),以對極點(diǎn)為中心公式的數(shù)值小于級數(shù)公式。圖1(b)上是100Hz的同樣的曲線比較,圖1(c)是圖1(b)的部分曲線。從圖1(b)、(c)這兩張圖可以清楚地看出,在大約1500km之后,兩條曲線幾乎完全一致,在大約1500km之內(nèi),以對極點(diǎn)為中心公式的值小于級數(shù)公式的值。這說明以對極點(diǎn)為中心公式的絕對值從離源點(diǎn)很近的距離開始直至對極點(diǎn)與級數(shù)公式的絕對值曲線符合十分良好,而且隨頻率的增高符合的近距離更縮短,符合的程度更接近。圖2與圖1只是在絕對值量級上大得多,而曲線趨勢變化與一階導(dǎo)數(shù)一致。

3 勒讓德函數(shù)對距離頭二階導(dǎo)數(shù)級數(shù)表達(dá)式與平面波導(dǎo)表達(dá)式對比分析

文獻(xiàn)[2～3]中給出勒讓德函數(shù)頭兩階導(dǎo)數(shù)平面近似公式如下[2～3]

以白天為例,分別取頻率30Hz、75Hz、100Hz和150Hz,在距離0km～20000km范圍進(jìn)行仿真,如圖3和圖4所示,其中實(shí)線‘-’代表級數(shù)公式計(jì)算結(jié)果,點(diǎn)線‘-*’代表駐波公式計(jì)算結(jié)果。

圖3 勒讓德函數(shù)對距離一階導(dǎo)數(shù)級數(shù)表達(dá)式和平面波導(dǎo)表達(dá)式對比分析

圖4 勒讓德函數(shù)對距離二階導(dǎo)數(shù)級數(shù)表達(dá)式和平面波導(dǎo)表達(dá)式對比分析

圖3(a)、(b)、(c)、(d)分別顯示20000km之內(nèi)30Hz、75Hz、100Hz和150Hz級數(shù)公式和平面公式計(jì)算的絕對值曲線。級數(shù)曲線有起伏,這是合理的,因?yàn)榧墧?shù)中包含直接波和非直接波,兩個(gè)波的合成必然引起起伏。平面公式是一個(gè)行波,所以是平滑下降的斜線,但由于曲率修正因子的影響,所以接近對極點(diǎn)曲線迅速上升。這說明平面公式在對極點(diǎn)附近不能應(yīng)用。在圖3(a)上示出30Hz級數(shù)計(jì)算和平面公式計(jì)算絕對值曲線的比較。在近距離上,兩條曲線的一致性比100Hz的差一點(diǎn),而且符合良好的距離比100Hz的近得多,這是合理的,因?yàn)?0Hz的衰減比100Hz的小得多,而150Hz符合度更好。圖4與圖3只是在絕對值量級上大得多,而曲線趨勢變化與一階導(dǎo)數(shù)一致。這些結(jié)果說明平面公式可以應(yīng)用于較高頻率,這個(gè)距離隨頻率增高而加大。

4 勒讓德函數(shù)對距離頭二階導(dǎo)數(shù)級數(shù)表達(dá)式與駐波波導(dǎo)表達(dá)式對比分析

文獻(xiàn)中給出勒讓德函數(shù)頭兩階導(dǎo)數(shù)的駐波公式,即兩個(gè)往相反方向傳播的行波公式的合成[2～3],即

+je-jk0Sρi]

因?yàn)槎A導(dǎo)數(shù)曲線變化趨勢與一階導(dǎo)數(shù)一致,只是在絕對值量級上比一階導(dǎo)數(shù)大得多。所以選擇白天,以發(fā)射頻率為100Hz和30Hz的一階導(dǎo)數(shù)為例,得到級數(shù)公式和駐波公式計(jì)算絕對值部分曲線的比較結(jié)果,如圖5所示,其中黑色實(shí)線‘-’代表級數(shù)公式計(jì)算結(jié)果,點(diǎn)線‘-*’代表駐波公式計(jì)算結(jié)果。

從圖5中看出,在大約2000km～19000km范圍內(nèi),這兩個(gè)公式的計(jì)算結(jié)果符合良好。在對極點(diǎn)附近,無論是100Hz還是30Hz駐波公式的絕對值都是急劇上升的,事實(shí)上電磁場在對極點(diǎn)不可能急劇上升的。在源點(diǎn)附近,駐波公式的絕對值小于級數(shù)公式的絕對值。

圖5 勒讓德函數(shù)對距離一階導(dǎo)數(shù)級數(shù)表達(dá)式和駐波波導(dǎo)表達(dá)式對比分析

5 結(jié)語

通過對無窮級數(shù)公式和常見的三個(gè)近似公式進(jìn)行編程計(jì)算,仿真分析相互比較,如果所用計(jì)算機(jī)是浮點(diǎn)16位尾數(shù)的,可以到到如下結(jié)論:

1) 一階導(dǎo)數(shù)的級數(shù)公式理論嚴(yán)密、精度高,是超低頻以下頻率電磁場計(jì)算的首選公式。該公式可以從200km～300km至對極點(diǎn)的距離范圍內(nèi)應(yīng)用,適用的頻率范圍大致在150Hz以下。

2) 當(dāng)工作頻率超過150Hz,可以應(yīng)用平面近似公式,但只適用于300km以外的距離,不能用于對極點(diǎn)附近的電磁場,其最遠(yuǎn)適用距離大致達(dá)15000km,而且這個(gè)距離范圍隨頻率升高而增加,適合于較高頻率。

3) 當(dāng)工作頻率超過150Hz,距離大致從1500km至對極點(diǎn)的范圍內(nèi),可用以對極點(diǎn)為中心的公式,而且隨頻率的增高符合的近距離更縮短,符合的程度更接近,可用于計(jì)算對極點(diǎn)附近的電磁場。

4) 當(dāng)工作頻率超過150Hz,距離大致從2000km開始至19000km范圍內(nèi),可以用駐波公式。

[1] 袁翊.超低頻和極低頻電磁波的傳播及噪聲[M].北京:國防工業(yè)出版社,2011. YUAN Yi. Propagation and Noise of SLF and ELF Electromagnetic Waves[M]. Beijing: National Defense Industry Press,2011.

[2] Casey J P. Extremely Low Frequency (ELF) Propagation Formulas for Dipole Souces Radiating in a Spherical Earth-Ionosphere Waveguide[J]. NUWC-NPT Technical Report,2002(5):732-737.

[3] King W P, Margaret O, Wu T T. Lateral Electromagnetic Waves[M]. Berlin: Springer-verlag press,1992.

[4] 《數(shù)學(xué)手冊》編寫組.數(shù)學(xué)手冊[M].北京:人民教育出版社,1979. 《Mathematics Handbook》Editorial committee. Mathematics Handbook[M]. Beijing: People’s Education Publishing Press,1979.

[5] 奚定平.貝塞爾函數(shù)[M].北京:高等教育出版社,1998. XI Dingping. Bessel Funtion[M]. Beijing: Higher Education Press,1998.

[6] Bakewell. Turbulent Wall-Pressure Fluctuations on a Body of Revolution[J]. The Journal of Acoustical Society of America,1968,43:1358-1363.

[7] Burrows M L. Motion-induced noise in electrode-palrextremely low frequency[J]. IEEE Trans. Com.,1974,22(4):540-542.

[8] 劉翠海,溫東,姜波,等.無線電通信系統(tǒng)仿真及軍事應(yīng)用[M].北京:國防工業(yè)出版社,2013. LIU Cuihai, WEN Dong, JIANG Bo, et al. Radio Communication System Simulation and Military Applications[M]. Beijing: National Defense Industry Press,2013.

[9] Bannister P R. Simplified formulas for ELF propagation at shorter distances[J]. Radio Science,1986,21(3):529-537.

[10] Greifinger C, Greifinger P. On the ionospheric parameters which govern high-latitude ELF propagation in the earth- ionosphere waveguide[J]. Radio Science,1979,14(5):845-889.

Comparative Analysis of Ultra-low Frequency Atmospheric Field Legendre Function Derivative Expressions

GAI Mingming1WEN Dong1XI Hailong2CHU Juan1

(1. Navy Submarine Academy, Qingdao 266199)(2. The Navy North Sea Fleet General Staff, Qingdao 266000)

Ultra-low frequency communication atmosphere field strength calculation includes the first kind of Legendre function to the first derivative of the distance or contains the second derivative, the two derivative can be expressed as an infinite series, but infinite series calculation is under the control of the digits limit of the computer, so series computing can’t be used in high-end of ultra-low frequency. In addition, there are several approximation formula. But the approximate formula of applicable range and precision are unclear. So infinite series formulas and the approximate formula of three common are programmed and calculated by soft ware based on Matlab, comparing each other, the applicable condition and precision difference formula is determined.

ultra-low frequency, atmospheric field strength, Legendre function to the first derivative of the distance, Legendre function to the second derivative of the distance

2016年9月11日,

2016年10月23日

蓋明明,女,碩士研究生,工程師,研究方向:潛艇通信。溫東,男,博士,副教授,研究方向:潛艇通信。郗海龍,男,碩士研究生,工程師,研究方向:軍事信息學(xué)。楚涓,女,碩士研究生,講師,研究方向:潛艇通信。

E962

10.3969/j.issn.1672-9722.2017.03.005