基于FEKO的非合作目標動態極化散射特性實時仿真*

吳盛源, 張小寬, 袁俊超, 林存坤

(1 空軍工程大學防空反導學院, 西安 710051; 2 95972部隊, 甘肅酒泉 735000)

基于FEKO的非合作目標動態極化散射特性實時仿真*

吳盛源1, 張小寬1, 袁俊超1, 林存坤2

(1 空軍工程大學防空反導學院, 西安 710051; 2 95972部隊, 甘肅酒泉 735000)

針對非合作目標動態測量難以開展的問題,研究了基于FEKO的雷達目標動態極化散射特性實時仿真方法。首先給出了隨時間變化的方位角和俯仰角的求解公式,然后介紹了利用FEKO實時仿真計算目標動態極化散射矩陣的方法及其步驟。仿真結果表明,機動目標的動態極化散射特性敏感于雷達視線角,雷達視線角變化越劇烈,插值法獲得的極化數據越偏離實時仿真數據。驗證了該實時仿真方法的有效性和準確性。

動態;實時;極化散射特性;FEKO;雷達視線角

0 引言

雷達目標的極化散射矩陣表征了全部的目標散射特性信息,是極化信息處理的基礎[1-2]。準確獲取目標的極化散射矩陣是研究極化檢測、極化抗干擾、極化目標識別的前提[3-8]。對于運動目標,不能用單一的極化散射矩陣表示,需要根據觀測視線的變化用時變的極化散射矩陣序列表示[9-10]。

獲取非合作目標動態極化散射特性的主要方法是:首先通過微波暗室獲得間隔1°的全空域數據,再根據目標機動航跡獲得隨時間變化的雷達視線角,最后利用插值的方法獲得動態極化散射矩陣。暗室縮比模型測量方法成本高,耗費時間長,采用插值的方法將產生誤差,無法獲得連續的動態極化散射特性?；谲浖姆抡娣椒ㄓ捎谄淇芍貜托院?、成本低等優點,越來越受到重視。因此,文中基于商業電磁仿真軟件FEKO,研究了雷達目標動態極化散射特性的軟件仿真方法。

1 雷達視線角

目標運動過程中,在雷達坐標系中坐標的變化將引起目標坐標系中雷達視線角的變化,而目標的極化散射特性取決于雷達視線角。因此,必須根據目標飛行航跡和目標姿態的變化,進行坐標變換,求出目標坐標系中隨時間變化的方位角和俯仰角,雷達坐標系和目標坐標系的定義以及具體的轉換過程可參考文獻[11-13]。

(1)

(2)

2 仿真方法

2.1 幾何建模

FEKO本身可以利用CADFEKO等自身軟件模塊來建立幾何模型;對于復雜目標建模,FEKO還可以通過與多個專業三維設計軟件(如catia,ansys等)接口,導入它們建立的三維模型,也可以直接導入劃分好網格的文件。

2.2 網格剖分

一般來說,網格單元邊長取1/8～1/10波長,而對于采用PO算法則可以適當取大一些邊長值。對于已經劃分網格的模型,可以在FEKO中重新定義網格,來滿足計算要求。網格劃分的大小決定了計算能否成功以及計算時間的長短,所以應盡量以節省時間而又能成功計算為準則。

2.3 算法選擇

FEKO[14-15]是一個以矩量法(MoM)為基礎的商業電磁仿真軟件,并且集成了多層快速多級子算法(MLFMM)、幾何繞射理論(GTD)、一致性繞射理論(UTD)和物理光學法(PO)等多種算法。提取目標極化散射矩陣需要獲取準確的幅度和相位信息,數值方法能求解非常復雜的散射體,同時具有很高的求解精度,因此應當采用精確數值計算方法。

隨著對矩陣方程快速求解問題的研究,研究者開發出了基于MoM的快速多級子算法(MLFMM),矩量法需要計算各個基函數間的相互作用,而MLFMM首先將基函數歸類分組,然后計算不同基函數組間的相互作用,如圖1所示。

MoM將N個基函數看成是相互孤立的,造成了N2規模的內存需求和N3規模的CPU消耗時間。而針對同一個求解問題,MLFMM算法只需要N·log(N)規模的內存需求和N·[log(N)]2規模的CPU消耗時間。

2.4 動態參數設置

提取目標的靜態極化散射矩陣,在CADFEKO中完成目標幾何建模、頻率設置、網格剖分、姿態角設置、單雙站設置和算法選擇后,首先設置水平極化發射,仿真生成的.out文件中EPHI即為散射矩陣中的Shh元素,ETHETA即為散射矩陣中的Svh元素。然后設置垂直極化發射,仿真生成的.out文件中EPHI即為散射矩陣中的元素Shv,ETHETA即為散射矩陣中的元素Svv。

提取運動目標的動態極化散射矩陣,首先在MATLAB中根據坐標轉換公式,以一個脈沖重復周期為間隔進行采點,獲取隨目標實時變化的雷達視線角,并保存在一個.dat文件中,.dat文件中有兩列數,第一列表示方位角,第二列表示俯仰角。EDITFEKO中的for卡和next卡提供了強大的循環調用功能,可以實時調用隨時間變化的雷達視線角進行計算,求得目標動態極化散射矩陣。在CADFEKO中完成目標幾何建模、頻率設置、網格剖分和算法選擇的基礎上,然后在EDITFEKO中設置參數的具體實現步驟如下:

1)在“** Sources”下面,添加for卡,for卡中設置循環變量名為“i”,以及循環的起始值、終止值和循環步長。

2)添加#卡,定義方位角變量為“phi”,讀取.dat文件,設置第i行、第一列的數為方位角;定義俯仰角變量為“theta”,讀取.dat文件,設置第i行、第二列數為俯仰角。

3)將原有的A0卡刪除,添加新的A0卡,在Number ofθangles后輸入:1,在Number ofφangles后輸入:1,在Initialθvalue后輸入:#theta,在Initialφvalue后輸入:#phi,設置幅度為1,相位為0,設置所需的極化方式。

4)添加FF卡,即設置遠場求解,并根據需要設置單雙站仿真。雙站角的設置方法與入射角的設置方法類似。

5)點擊RUN FEKO,運行FEKO進行仿真計算。

對于同時全極化測量體制雷達,首先設置水平極化發射,獲得動態極化散射矩陣第一列元素,然后設置垂直極化發射,獲得動態極化散射矩陣第二列元素。

3 動態極化特性仿真分析

目標朝向雷達站飛行并做偏航機動,其機動航跡為拋物線,且機頭方向為拋物線的切線方向,起點坐標為(80 km,0 km,5 km),終點坐標為(5 km,6 km,5 km)。圖2～圖4給出了其機動曲線、雷達視線角變化,目標動態多極化RCS如圖5所示,圖6對比了文中提出的實時仿真方法與插值法獲得的極化比。

由圖5可知,利用文中提出的仿真方法,能仿真獲得任意時刻的目標極化散射特性。說明對于高速機動目標,不能用單一的極化散射矩陣描述極化特性,必須研究其實時動態極化散射特性;由圖6可知,插值法獲得的極化比在實時仿真獲得的極化比上下波動,在前140 s,雷達視線角變化緩慢,極化比變化緩慢,兩種方法的差值較小,最大差值為4。在140 s之后,雷達視線角變化劇烈,極化比劇烈變化,兩種方法的差值較大,最大差值達到25。說明當飛機做超機動時,采用插值法研究其動態極化散射特性將帶來較大的誤差,采用文中提出的實時仿真方法,能夠保證獲得的極化散射特性的準確性。

4 結論

文中研究了基于FEKO軟件的目標動態極化散射特性提取方法,并利用該方法仿真獲得偏航機動下的目標實時動態極化散射特性。仿真結果表明,文中提出的軟件實時仿真方法可以提取準確的目標動態極化散射矩陣,是研究非合作目標動態極化散射特性的有效方法。

[1] 戴幻堯, 劉勇, 周波, 等. 極化雷達測量體制研究進展 [J]. 現代雷達, 2014, 36(1): 1-5.

[2] 李永禎, 李棉全, 程旭, 等. 雷達極化測量體制研究綜述 [J]. 系統工程與電子技術, 2013, 35(9): 1873-1876.

[3] TOUZI R, HURLEY J, VACHON P W. Optimization of the degree of polarization for enhanced ship detection using polarimetric RADARSAT-2 [J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(10): 5403-5424.

[4] XIAO J J, NEHORAI A. Joint transmitter and receiver polarization optimization for scattering estimation in clutter [J]. IEEE Transactions on Signal Process, 2009, 57(10): 4142-4147.

[5] 張強. 有源誘餌的雙極化識別 [J]. 電訊技術, 2015, 55(2): 168-174.

[6] 常宇亮. 瞬態極化雷達測量、檢測與抗干技術研究 [D]. 長沙: 國防科學技術大學, 2010.

[7] ZHANG Qinyu, CAO Bin, WANG Jian, et al. Polarization filtering technique based on oblique projections [J]. Science China Information Sciences, 2010, 53(5): 1056-1066.

[8] 劉勇. 動態目標極化特性測量與極化雷達抗干擾新方法研究 [D]. 長沙: 國防科學技術大學, 2011.

[9] LIU Yong, LI Yongzhen, WANG Xuesong. Instantaneous measurement for radar target polarization scattering matrix [J]. Journal of Systems Engineering and Electronics, 2010, 21(6): 968-974.

[10] 王雪松. 寬帶極化信息處理的研究 [D]. 長沙: 國防科學技術大學, 1999.

[11] 黃培康, 殷紅成, 許小劍. 雷達目標特性 [M]. 北京: 電子工業出版社, 2005: 130-132.

[12] 周超, 張小寬, 張敬偉, 等. 典型隱身飛機動態RCS時間序列研究 [J]. 空軍工程大學學報(自然科學版), 2013, 14(3): 15-18.

[13] 周超, 張小寬, 吳國成. 基于坐標轉換目標動態RCS時間序列研究 [J]. 火力與指揮控制, 2014, 39(3): 56-59.

[14] 李亭, 韓紅斌, 劉宇. 基于FEKO的雷達目標電磁散射特性研究 [J]. 現代電子技術, 2015, 38(21): 39-41.

[15] 郭鵬, 白亮, 武夢潔, 等. 基于FEKO的雷達散射截面實時計算 [J]. 航空科學技術, 2013(6): 72-75.

Real-timeSimulationoftheDynamicPolarizationScatteringCharacteristicsofNoncooperativeTargetBasedonFEKO

WU Shengyuan1, ZHANG Xiaokuan1, YUAN Junchao1, LIN Cunkun2

(1 Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China; 2 No.95972 Unit, Gansu Jiuquan 735000, China)

The real time simulation method of dynamic polarization scattering characteristics of radar target based on FEKO was studied to solve the problem that the dynamic measurement of noncooperative targets was difficult to carry out. Firstly, the the formula of azimuth and elevation angle that varied with time variation were given. Secondly, the method and the procedure of the real time simulation and calculation of dynamic polarization scattering matrix of the targets using FEKO were introduced. The simulation results showed that the dynamic polarization scattering characteristics of maneuvering targets were more sensitive than radar aspect angle, and the more severe the change of radar aspect angle, the bigger the deviation of the polarization data obtained by interpolation from the real time simulation data was. The validity and accuracy of the real time simulation method were verified.

dynamic; real time; polarization scattering characteristics; FEKO; radar aspect angle

TN974

A

2016-08-29

重點實驗室基金(STES201401-2)資助

吳盛源(1991-),男,福建漳州人,碩士研究生,研究方向:雷達目標特性及其應用。