基于FEKO的非合作目標(biāo)動(dòng)態(tài)極化散射特性實(shí)時(shí)仿真*

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

(1 空軍工程大學(xué)防空反導(dǎo)學(xué)院, 西安 710051; 2 95972部隊(duì), 甘肅酒泉 735000)

基于FEKO的非合作目標(biāo)動(dòng)態(tài)極化散射特性實(shí)時(shí)仿真*

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

(1 空軍工程大學(xué)防空反導(dǎo)學(xué)院, 西安 710051; 2 95972部隊(duì), 甘肅酒泉 735000)

針對(duì)非合作目標(biāo)動(dòng)態(tài)測(cè)量難以開展的問題,研究了基于FEKO的雷達(dá)目標(biāo)動(dòng)態(tài)極化散射特性實(shí)時(shí)仿真方法。首先給出了隨時(shí)間變化的方位角和俯仰角的求解公式,然后介紹了利用FEKO實(shí)時(shí)仿真計(jì)算目標(biāo)動(dòng)態(tài)極化散射矩陣的方法及其步驟。仿真結(jié)果表明,機(jī)動(dòng)目標(biāo)的動(dòng)態(tài)極化散射特性敏感于雷達(dá)視線角,雷達(dá)視線角變化越劇烈,插值法獲得的極化數(shù)據(jù)越偏離實(shí)時(shí)仿真數(shù)據(jù)。驗(yàn)證了該實(shí)時(shí)仿真方法的有效性和準(zhǔn)確性。

動(dòng)態(tài);實(shí)時(shí);極化散射特性;FEKO;雷達(dá)視線角

0 引言

雷達(dá)目標(biāo)的極化散射矩陣表征了全部的目標(biāo)散射特性信息,是極化信息處理的基礎(chǔ)[1-2]。準(zhǔn)確獲取目標(biāo)的極化散射矩陣是研究極化檢測(cè)、極化抗干擾、極化目標(biāo)識(shí)別的前提[3-8]。對(duì)于運(yùn)動(dòng)目標(biāo),不能用單一的極化散射矩陣表示,需要根據(jù)觀測(cè)視線的變化用時(shí)變的極化散射矩陣序列表示[9-10]。

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

1 雷達(dá)視線角

目標(biāo)運(yùn)動(dòng)過程中,在雷達(dá)坐標(biāo)系中坐標(biāo)的變化將引起目標(biāo)坐標(biāo)系中雷達(dá)視線角的變化,而目標(biāo)的極化散射特性取決于雷達(dá)視線角。因此,必須根據(jù)目標(biāo)飛行航跡和目標(biāo)姿態(tài)的變化,進(jìn)行坐標(biāo)變換,求出目標(biāo)坐標(biāo)系中隨時(shí)間變化的方位角和俯仰角,雷達(dá)坐標(biāo)系和目標(biāo)坐標(biāo)系的定義以及具體的轉(zhuǎn)換過程可參考文獻(xiàn)[11-13]。

(1)

(2)

2 仿真方法

2.1 幾何建模

FEKO本身可以利用CADFEKO等自身軟件模塊來建立幾何模型;對(duì)于復(fù)雜目標(biāo)建模,FEKO還可以通過與多個(gè)專業(yè)三維設(shè)計(jì)軟件(如catia,ansys等)接口,導(dǎo)入它們建立的三維模型,也可以直接導(dǎo)入劃分好網(wǎng)格的文件。

2.2 網(wǎng)格剖分

一般來說,網(wǎng)格單元邊長取1/8～1/10波長,而對(duì)于采用PO算法則可以適當(dāng)取大一些邊長值。對(duì)于已經(jīng)劃分網(wǎng)格的模型,可以在FEKO中重新定義網(wǎng)格,來滿足計(jì)算要求。網(wǎng)格劃分的大小決定了計(jì)算能否成功以及計(jì)算時(shí)間的長短,所以應(yīng)盡量以節(jié)省時(shí)間而又能成功計(jì)算為準(zhǔn)則。

2.3 算法選擇

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

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

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

2.4 動(dòng)態(tài)參數(shù)設(shè)置

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

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

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

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

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

4)添加FF卡,即設(shè)置遠(yuǎn)場(chǎng)求解,并根據(jù)需要設(shè)置單雙站仿真。雙站角的設(shè)置方法與入射角的設(shè)置方法類似。

5)點(diǎn)擊RUN FEKO,運(yùn)行FEKO進(jìn)行仿真計(jì)算。

對(duì)于同時(shí)全極化測(cè)量體制雷達(dá),首先設(shè)置水平極化發(fā)射,獲得動(dòng)態(tài)極化散射矩陣第一列元素,然后設(shè)置垂直極化發(fā)射,獲得動(dòng)態(tài)極化散射矩陣第二列元素。

3 動(dòng)態(tài)極化特性仿真分析

目標(biāo)朝向雷達(dá)站飛行并做偏航機(jī)動(dòng),其機(jī)動(dòng)航跡為拋物線,且機(jī)頭方向?yàn)閽佄锞€的切線方向,起點(diǎn)坐標(biāo)為(80 km,0 km,5 km),終點(diǎn)坐標(biāo)為(5 km,6 km,5 km)。圖2～圖4給出了其機(jī)動(dòng)曲線、雷達(dá)視線角變化,目標(biāo)動(dòng)態(tài)多極化RCS如圖5所示,圖6對(duì)比了文中提出的實(shí)時(shí)仿真方法與插值法獲得的極化比。

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

4 結(jié)論

文中研究了基于FEKO軟件的目標(biāo)動(dòng)態(tài)極化散射特性提取方法,并利用該方法仿真獲得偏航機(jī)動(dòng)下的目標(biāo)實(shí)時(shí)動(dòng)態(tài)極化散射特性。仿真結(jié)果表明,文中提出的軟件實(shí)時(shí)仿真方法可以提取準(zhǔn)確的目標(biāo)動(dòng)態(tài)極化散射矩陣,是研究非合作目標(biāo)動(dòng)態(tài)極化散射特性的有效方法。

[1] 戴幻堯, 劉勇, 周波, 等. 極化雷達(dá)測(cè)量體制研究進(jìn)展 [J]. 現(xiàn)代雷達(dá), 2014, 36(1): 1-5.

[2] 李永禎, 李棉全, 程旭, 等. 雷達(dá)極化測(cè)量體制研究綜述 [J]. 系統(tǒng)工程與電子技術(shù), 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] 張強(qiáng). 有源誘餌的雙極化識(shí)別 [J]. 電訊技術(shù), 2015, 55(2): 168-174.

[6] 常宇亮. 瞬態(tài)極化雷達(dá)測(cè)量、檢測(cè)與抗干技術(shù)研究 [D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 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òng)態(tài)目標(biāo)極化特性測(cè)量與極化雷達(dá)抗干擾新方法研究 [D]. 長沙: 國防科學(xué)技術(shù)大學(xué), 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]. 長沙: 國防科學(xué)技術(shù)大學(xué), 1999.

[11] 黃培康, 殷紅成, 許小劍. 雷達(dá)目標(biāo)特性 [M]. 北京: 電子工業(yè)出版社, 2005: 130-132.

[12] 周超, 張小寬, 張敬偉, 等. 典型隱身飛機(jī)動(dòng)態(tài)RCS時(shí)間序列研究 [J]. 空軍工程大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 14(3): 15-18.

[13] 周超, 張小寬, 吳國成. 基于坐標(biāo)轉(zhuǎn)換目標(biāo)動(dòng)態(tài)RCS時(shí)間序列研究 [J]. 火力與指揮控制, 2014, 39(3): 56-59.

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

[15] 郭鵬, 白亮, 武夢(mèng)潔, 等. 基于FEKO的雷達(dá)散射截面實(shí)時(shí)計(jì)算 [J]. 航空科學(xué)技術(shù), 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

重點(diǎn)實(shí)驗(yàn)室基金(STES201401-2)資助

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