淺海中時諧水平電偶極子定位研究*

李 濤 王向軍 嵇 斗

(海軍工程大學電氣工程學院 武漢 430033)

?

淺海中時諧水平電偶極子定位研究*

李 濤 王向軍 嵇 斗

(海軍工程大學電氣工程學院 武漢 430033)

淺海中時諧水平電偶極子定位技術是利用艦船軸頻電場信號對艦船進行定位的理論基礎。文章提出了一種基于遺傳算法的時諧水平電偶極子的位置尋優方法。只需在海水中布置一個電場傳感器,測得相應位置的X、Y、Z方向電場強度模值。在合適的坐標系下,結合淺海中時諧水平電偶極子的電場傳播公式和遺傳優化算法建立位置計算模型,就能精確地確定該偶極子在海水中的位置。給出了仿真計算,結果表明該定位方法具有較高的定位精度和可行性。

時諧水平電偶極子; 淺海; 遺傳算法; 定位

Class Number TM154.3

1 引言

由于腐蝕、船體材料構成的不同、陰極防腐措施及運動艦船切割地磁場等因素,無論艦船是運動的還是靜止的,其周圍海水中都存在著不可避免的艦船電場,艦船電場通常分為兩個部分:靜電場和交變電場[1～2]。水下電場信號可被用在電場兵器、海礦產探測、水下考古、水下搜救等多個領域。隨著技術的發展,可將基于電場信號的非聲探測技術作為傳統水聲探測技術的補充[3～5]。基于艦船靜電場信號對水中目標進行定位和識別的技術已有一些研究成果[5],但由于靜電場信號包含特征較少,易受外界干擾等自身特性的限制[6～7],使用該信號對艦船進行定位與識別受到一定的限制,相比而言,艦船的軸頻電磁場屬于極低頻(ELF)/甚低頻(ULF)電磁場,它具有衰減速度慢、傳播距離遠等特點,是一種可以被遠程探測的艦船物理場[8～10]。因此,軸頻電磁場在水中兵器這一領域中的應用價值就在于它可以被用于遠程攻擊武器的引信中。

把海水看做線性、均勻、各向同性的無限大半空間,艦船在海水中運動時產生的軸頻電磁場可以用低頻時諧電偶極子進行建模[1,8,12]。所以,利用時諧電偶極子產生的電磁場對其進行定位是研究利用艦船在海水中產生的軸頻電場信號對艦船進行跟蹤定位的基礎。目前,已有人提出了一種基于傳感器陣列對電偶極子進行定位的方法,但該方法必須在海水中布放大量的傳感器,且每個傳感器之間有精確的位置排列關系,實際操作難度較大[11]。本文針對時諧水平電偶極子在無限大海水中產生的極低頻電場,在前人推導求得的時諧水平電偶極子的解析表達式的基礎上,利用漢克爾變換式的FFT算法在Matlab上進行編程計算,再利用遺傳優化算法進行建模,在測量出某點X、Y、Z分量電場強度模值的基礎上,成功解算出該時諧水平電偶極子在海水空間中的位置坐標,此方法僅需在海水中布放一個高精度電場傳感器,實際操作簡便且具有較高的精度。

2 時諧水平電偶極子的電場定位模型

如圖1所示,取兩直角坐標系S與S′,其中S坐標系是以時諧水平電偶極子所處位置為坐標原點,且x軸、y軸平行于海平面,z軸垂直于海平面,方向向下。S′坐標系x軸、y軸、z軸與S坐標系相互平行且同向,x軸、y軸位于海平面,海平面上z=0。設海水深度為d,S′以坐標下,設z<0為0區,0d為2區,0區、1區和2區分別為空氣、海水和海底。

在深度為h的海水中放置一極低頻時諧水平電偶極子P,該電偶極子在S坐標系中坐標為(0,0,0),在S′坐標系中該電偶極子的z坐標為h。在深度為l的海水中放置一電場傳感器A,該電場傳感器在S坐標系下坐標未知,將A點與S′坐標系的z軸重合,則可得電場傳感器A在S′坐標系下的坐標為(0,0,l)。

在電場傳感器A處測得該時諧水平電偶極子的電場強度X、Y、Z分量模值分別為E1x、E1y、E1z,將該三分量電場強度模值作為已知量,利用文獻[12～13]淺海中時諧水平電偶極子的電場解析式解算出A點在S坐標系下的坐標為(x1,y1,z1)。利用坐標變換,便可確定出該電偶極子P在S′坐標系下的坐標為(-x1,-y1,l-z1),即實現了電偶極子P的定位。

圖1 電偶極子定位的坐標關系

3 淺海中時諧水平電偶極子的電場解析式

根據文獻[12～13],淺海中時諧水平電偶極子有如下的電場解析式:

(1)

(2)

(3)

4 程序設計

本文仿真計算分為兩步,首先利用合適的積分方法編程求解式(1)～式(3),計算出空間某一點的電場強度模值,然后借助遺傳算法尋優出式中的初始變量,即該點位置坐標。

4.1 漢克爾變換式求解廣義索末菲積分

式(1)～式(3)中大部分參數容易求出,對于式中的廣義索末菲積分求解,文獻[14～16]中提出了一種基于漢克爾變換式的FFT算法用于求解廣義索末菲積分,它具有精度高、速度快等其他數值積分方法無法比擬的優點。因此,本文采用了這種方法來對電場解析式進行數值計算。

n階漢克爾變換的定義式為

(4)

式中Jn(ρξ)為第一類n階貝賽爾函數,積分求解時需用數值濾波計算方法。

Guptasama和Singh線性數值濾波計算公式為

(5)

ξi=(1/ρ)×10[a+(i-1)s],i=1,2,…,n

(6)

數值計算時,計算精度取決于積分區間的長度n、抽樣點的位置ξi和加權系數Wi。通常n越大計算精度越高。本文采用的是Guptasama和Singh給出的61點漢克爾J0變換線性濾波器和47點漢克爾J1變換線性濾波器系數,據Guptasama和Singh測評,該線性濾波器具有很高的計算精度[16]。該濾波系數表參照文獻[16]。

4.2 遺傳優化算法求解源點位置坐標

遺傳算法(簡稱GA)是以自然選擇和遺傳理論為基礎的高效的全局尋優搜索算法,為求解非線性、多模型、多目標等復雜系統優化問題通用框架,Matlab的GA工具箱中有大量的遺傳算法程序可供調用[17～19]。

為確保較好的定位精度,預設誤差范圍d<1×10-13,經多次反復試驗后,確定選擇率為0.7,變異率為0.1,多次迭代,即可輸出場點坐標最優解(x1,y1,z1)。此外程序中還可以跟蹤算法性能,便于算法優化。

最后,利用坐標變換即可求解出該時諧水平電偶極子在S′坐標系下的坐標(-x1,-y1,l-z1)。程序流程圖如圖2所示。

圖2 程序設計流程圖

5 實例計算

假定在淺海中有一時諧水平電偶極子,大小為5C·m,其角頻率為3rad/s,方向沿x軸正方向,深度為30m,在S坐標系下該電偶極子坐標為(0,0,0)。取固定場點A的深度l=50m,A在S坐標系中的坐標為(5000,8000,20)。取空氣的三個電磁參數分別為:σ0=0,ε0=(1/36π)×10-9F/m,μ0=4π×10-7H/m,海水的為:σ1=4s/m,ε1=80ε0,μ1=μ0。

根據式(1)～式(3),利用漢克爾變換式的FFT算法在計算機上求得S坐標系下(5000,8000,20)點電場強度模值為:E1x=4.208×10-10,E1y=1.1099×10-10,E1z=4.02×10-15。圖3為時諧水平電偶極子在固定另外兩個方向坐標的情況下電場模值E1x、E1y、E1z在沿x、y、z方向變化曲線。

圖3 時諧水平電偶極子電場模值E1x,E1y,E1z在沿x,y,z方向變化曲線

場點坐標仿真尋優時,將E1x、E1y、E1z作為已知參考量,由于E1z數值較小,與E1x、E1y差了四個數量級,所以在計算中為了減小z1的誤差,將E1Z放大了104倍使之與E1x,E1y數值接近。設定尋優范圍100

圖4 迭代過程中誤差d和x,y,z坐標的變化過程

6 結語

本文在前人的基礎上,利用漢克爾變換式的FFT算法在計算機上求得淺海中時諧水平電偶極子空間任意位置A的電場模值,并以此作為仿真計算的已知量,利用遺傳優化算法迭代尋找出該點在S坐標系下的坐標,最后利用坐標變換求解出了電偶極子源在S′坐標系下的坐標。運用此方法來實現介質中的電偶極子的定位,使用傳感器較少,操作簡便并具有較高的定位精度,且此方法對垂直時諧偶極子同樣適用,為利用軸頻電場對艦船進行跟蹤定位提供了新的思路。

[1] 龔沈光,盧新城.艦船電場特性初步分析[J].海軍工程大學學報,2008,2(20):1-5.

[2] 嵇斗,王向軍,楊振.一種新的補償陽極優化方法[J].兵工學報,2010,12:1562-1566.

[3] Xiang Jun Wang, Shu Zhang, Yi Liu, et al. The Electrochemical principle of measuring underwater elelctric field[J]. Applied Mechanics and Materials,2012,232:863-868.

[4] Xiang Jun Wang, Yi Liu, Dou Ji. The Research of Electric Field Based on Ring Current Source Model[J]. Applied Mechanics and Materials,2012,622:661-664.

[5] 張華,王向軍,單潮龍,等.基于目標靜電場的水中兵器制導方法研究[J].電子學報,2013,41(3):470-474.

[6] 劉文寶,王向軍,嵇斗.基于電偶極子模型的艦船靜電場深度換算[J].空軍雷達學院學報,2010,6:435-438.

[7] Dou Ji, Xiang Jun Wang, Wen Bao Liu. The Propagation of Horizontal Static Dipole Electromagnetic Field in Shallow Sea[C]//Electromagnetic Field Problems and Applications(ICEF),2012,43.

[8] 汪小娜,肖昌漢,王向軍,等.運動潛艇的感應電場分布[J].海軍工程大學學報,2012,24(4):7-10.

[9] 汪小娜,肖昌漢,王向軍,等.沿任意方向運動的磁偶極子感應電場特性[J].艦船科學技術,2012,12:19-23.

[10] 朱武兵,嵇斗,王向軍,等.淺海中影響運動艦船軸頻電磁場的因素[J].船電技術,2013,11:41-44.

[11] Wang Xiaobei, Xia Li, Wang Xiangjun. Design and Implementation of High-precision underwater elelctric field measurement system[C]//Proceedings of 4thInternational Conference on Modeling Identification an Control. Wuhan,2012:2102-2105.

[12] 毛偉,張寧,林春生.在三層介質中運動的時諧水平電偶極子產生的電磁場[J].電子學報,2009,9:2077-2081.

[13] 雷銀照.時諧電磁場解析方法[M].北京:科學出版社,2000:73-136.

[14] 胡俊,聶在平.索末菲爾德積分新方法_快速漢克爾變換[J].電子學報,1998,26(3):126-128.

[15] 李俊杰,嚴家斌.漢克爾變換精度分析及在電磁法探測中的應用[J].物探化探計算技術,2013,3:288-292.

[16] 阮百堯.均勻水平大地上頻率域垂直磁偶源電磁場數值濾波解法[J].桂林工學院學報,2005,25(1):14-18.

[17] 雷英杰,張善文.MTALAB遺傳算法工具箱及應用[M].西安:西安電子科技大學出版社,2005:117-145.

[18] 張民,王向軍,嵇斗.進化規劃中的變異與收斂[J].海軍工程大學學報,2007,19(1):48-53.

[19] 張民,王向軍,嵇斗.一種新的進化規劃算法[J].海軍工程大學學報,2008,20(3):40-43.

Time-harmonic Horiziontal Dipole in Coastal

LI Tao WANG Xiangjun JI Dou

(College of Electronic Engineering, Naval University of Engineering, Wuhan 430033)

The positioning technology of a time-harmonic horiziontal dipole in coastal is the theoretical basis of ship’s localization based on the ship’s shaft-rate electric field signal. An approach based on genetic algorithm is proposed to search the location of the dipole. Only an electric field sensor in seawater is needed to set up measure the modulus of electric field intensity inX/Y/Zdirections at the corresponding positions. A model for position calculation combine the propagation formula of a time-harmonic horiziontal dipole in coastal with genetic optimal algorithm in an appropriate coordinate system is set up, then the position of the dipole can be determined accurately. The simulation results show the proposed dipole localization method is accurate and feasible.

time-harmonic horiziontal dipole, coastal, genetic algorithm, localization

2014年11月2日,

2014年12月9日

中國博士后科學基金項目(編號:20080441290);中國博士后科學基金特別項目(編號:200902670);國家自然科學基金青年項目(編號:51107145)資助。 作者簡介:李濤,男,碩士研究生,研究方向:電場環境與防護技術。王向軍,男,博士,教授,研究方向:電場環境與防護技術。嵇斗,男,碩士,副教授,研究方向:電場環境與防護技術。

TM154.3

10.3969/j.issn1672-9730.2015.05.039