機載天線方向圖回歸研究

機載天線方向圖回歸研究

張建華李振亞

(電子工程學(xué)院,安徽 合肥 230037)

摘要利用工程軟件建立無人機及其天線的一體化仿真模型,綜合考慮飛行姿態(tài)、機載天線對地面接收天線的極化匹配和方向性等因素,研究了機載天線對地面目標的干擾方向圖.采用遺傳算法,結(jié)合機載天線仿真的數(shù)值結(jié)果,設(shè)計適應(yīng)度函數(shù),進行方向圖回歸研究,得到了便于實際應(yīng)用的解析公式.

關(guān)鍵詞機載天線;方向圖;天線極化;極化匹配;曲線回歸

中圖分類號TN82

文獻標志碼A

文章編號1005-0388(2015)04-0777-06

AbstractBased on the established simulation model that integrates an unmanned aerial vehicle (UAV) with its airborne antenna, the interference patterns of the airborne antenna for ground target were studied with engineering software, considering the flight attitude, the polarization match of the airborne antenna for ground receiving antennas and the direction pattern of the airborne antenna. Combined with the numerical results of the airborne antenna simulation, a fitness function was proposed for the regression of the interference patterns by genetic algorithm. An analytical formula is obtained which is believed to be suitable for practical application.

收稿日期：2014-09-01

作者簡介

Regression analysis of airborne antenna pattern

ZHANG JianhuaLI Zhenya

(ElectronicEngineeringInstitute,HefeiAnhui230037,China)

Key wordsairborne antenna; radiation pattern; antenna polarization; polarization matching; curvilinear regression

引言

無人機可以作為電子對抗系統(tǒng)的作戰(zhàn)平臺,飛臨敵方縱深數(shù)千米高空實施干擾,具有抵近作戰(zhàn)優(yōu)勢[1-2],與地面設(shè)備相比,減小了路徑傳輸損耗,天線升空可以增大電波傳播視距.機載天線作為無人機系統(tǒng)的重要組成部件之一,其性能直接影響系統(tǒng)的整體效能.對于無人機運用的非天線專業(yè)人員,在進行無人機航路規(guī)劃、干擾效能評估時,迫切需要能夠直接進行解析計算的有關(guān)機載天線方向圖的理論公式,因此有必要進行機載天線方向圖回歸的研究,這涉及到以下幾方面的問題:

聯(lián)系人： 張建華 E-mail: zhangjh12@gmail.com

一是對天線而言機體結(jié)構(gòu)是一個復(fù)雜的電磁散射體,其影響將使天線電參數(shù)發(fā)生改變[3-5],研究方法有:并行時域有限差分法[3]、一致性幾何繞射理論[4]、矩量法結(jié)合物理光學(xué)法[5]、有限元法[6]、時域有限差分法結(jié)合幾何繞射理論[7]等方法,本文采用工程軟件CATIA[8]和FEKO[9]相結(jié)合的方法,建立機載天線的仿真模型,計算機載天線輻射場的全面電參數(shù).二是飛行姿態(tài)的變化造成機載發(fā)射天線與被干擾天線之間極化失配,且失配度處于不斷變化之中,本文依據(jù)文獻[10]的極化匹配理論,對機載天線在各種飛行姿態(tài)下的極化匹配因子進行了推導(dǎo).三是如何進行方向圖回歸計算以及構(gòu)建方向圖回歸適應(yīng)度函數(shù),遺傳算法已被廣泛應(yīng)用于天線相關(guān)參數(shù)的優(yōu)化[11-15],如采用遺傳算法進行天線陣方向圖綜合[11],優(yōu)化陣列分布降低旁瓣電平和主瓣寬度[12],優(yōu)化機載天線位置降低其電磁耦合[13],對天線結(jié)構(gòu)進行優(yōu)化設(shè)計[14-15]等,因此本文提出了基于改進的量子遺傳算法 (Improved Quantum Genetic Algorithm,IQGA)的天線方向圖回歸算法,獲得了機載天線對于地面目標的干擾方向圖回歸解析式.

1機載天線一體化建模

機載天線處在一個復(fù)雜的電磁環(huán)境中,由于機體的電磁散射作用,使得天線的電參數(shù)發(fā)生變化.機載天線電參數(shù)除了與天線自身結(jié)構(gòu)有關(guān)外,還與無人機機體的外形結(jié)構(gòu)和天線的安裝位置有關(guān).一方面,機體本身可以作為機載天線的一部分,天線工作時會激勵起表面感應(yīng)電流分布；另一方面,機載天線在機體不同的安裝位置導(dǎo)致機體對天線繞射場和反射場的產(chǎn)生,這些綜合作用將影響天線的性能.

采用CATIA[8]對無人機及其機載天線進行一體化建模,CATIA具有很強的曲面構(gòu)造功能,便于對無人機進行建模,無人機機頭、機身和垂尾等復(fù)雜曲面采用高級曲線擬合,可以較準確地建立機載平臺模型.采用CATIA對無人機主體曲面擬合建模后,保存成后綴.model文件,然后結(jié)合FEKO[9]軟件進行了機載天線的仿真分析,即采用FEKO5.4導(dǎo)入該.model文件.FEKO軟件基于矩量法,擁有高效的多層快速多極子法,并將矩量法與高頻分析方法相結(jié)合,可以快速、準確地分析各類天線.

采用CATIA+FEKO的方法,建立了機載天線和機身平臺的一體化仿真模型后,就可以計算獲得機載天線的方向圖、增益、輸入阻抗和軸比等天線電參數(shù).

2機載天線的干擾方向圖函數(shù)

機載天線在實際運用中,我們關(guān)注其對不同位置處目標的干擾效能,為此,需要研究目標處可能被接收的干擾功率密度.下面將機載天線與收發(fā)天線的極化匹配因子相結(jié)合進行研究.

2.1干擾方向圖函數(shù)的定義

設(shè)發(fā)射機輸出功率為PT,則輻射至干擾目標處的功率密度為

(1)

式中: GT為機載天線的增益; FT為機載天線的歸一化方向函數(shù); r為機載天線至目標處的距離; Δ為俯視角; φ為方位角.

對于地面目標,收發(fā)天線之間的距離r與無人機飛行高度h和俯視角Δ的關(guān)系為r=h/sinΔ,代入式(1)可得

(2)

目標處可能被接收的干擾功率密度除了與式(2)的功率密度有關(guān)外,還與收發(fā)天線間的極化匹配因子有關(guān),因此將極化匹配因子ρ 乘以式(2),可得目標處可被接收的功率密度為

(3)

式中,α、β 和γ分別為滾轉(zhuǎn)角、俯仰角和偏航角,如圖1所示.式中極化匹配因子ρ將在2.2節(jié)推導(dǎo).

圖1　無人機坐標系統(tǒng)

式(3)中的參數(shù)可以分為兩大類:一類是與機載天線方向性、極化特性以及無人機飛行姿態(tài)變化無關(guān)的參量;另一類則是與其有關(guān)的參量,我們將這一類與方向角、飛行姿態(tài)角有關(guān)的量定義為機載天線的干擾功率密度方向圖函數(shù),簡稱干擾方向圖函數(shù)φ,表達式為

(4)

將干擾方向圖函數(shù)g代入式(3),可得目標處可被接收的功率密度

(5)

干擾方向圖函數(shù)式(4)中極化匹配因子ρ(Δ,φ,α,β,γ)由第2.2節(jié)推導(dǎo),FT(Δ,φ)隨飛行姿態(tài)的變化將在第2.3節(jié)中進行分析.

2.2機載天線極化匹配因子的分析

將天線實際接收的功率與在同方向、同強度且極化完全匹配條件下天線接收的功率之比定義為極化匹配因子.極化匹配因子可用如下表達式來計算[10,16]

(6)

(7)

(8)

(9)

設(shè)地面接收天線的振子軸線方向為aR,文獻[18]導(dǎo)出極化匹配因子為

(10)

(11)

ξ(tan2Δ+cos2φ);

(12)

(13)

式中,

(14)

2.3歸一化方向函數(shù)FT隨飛行姿態(tài)的變化

隨著飛行姿態(tài)的改變,歸一化方向函數(shù)FT是一個變化的函數(shù).設(shè)以天線為參考系的球坐標中極角和方位角為(θ1,φ1),在z軸垂直于地面的球坐標中極角和方位角為(θ,φ),其中極角θ與俯視角Δ關(guān)系為θ=Δ+π/2.當無人機正常飛行時兩種坐標系是一致的,當飛行姿態(tài)改變時,天線發(fā)生轉(zhuǎn)動,同一個函數(shù)值FT(θ1,φ1)指向了FT(θ,φ),利用旋轉(zhuǎn)矩陣T可得角度變換關(guān)系為

(15)

這里需要特別指出的是,式(15)角度關(guān)系表達的是同一個函數(shù)值FT由方向(θ1,φ1)指向了不同的方向(θ,φ),即FT(θ1,φ1)=FT(θ,φ),由前面的建模仿真可獲得FT(θ1,φ1).將式(7)～(9)分別代入式(15),可得角度變換關(guān)系如下:

1) 滾轉(zhuǎn)姿態(tài)的角度關(guān)系

θ1=cos-1(cosαcosθ-sinαsinθsinφ);

(16)

φ1= angle(sinθcosφ,sinαcosθ+cosαsinθsinφ).

(17)

式中,函數(shù)angle表示對任意兩個變量var1、var2的計算式,

angle(var1,var2)=arctan(var2/var1).

(18)

2) 俯仰姿態(tài)的角度關(guān)系

θ1=arccos(sinβsinθcosφ+cosβcosθ);

(19)

φ1= angle(cosβsinθcosφ-sinβcosθ,sinθsinφ).

(20)

3) 偏航姿態(tài)的角度關(guān)系為θ1=θ,φ1=φ-γ.

3機載天線的干擾方向圖回歸研究

干擾方向圖綜合反映了目標處可被接收的干擾功率密度與機載天線的方向性、無人機飛行姿態(tài)、極化匹配特性、目標位置等參量的關(guān)系,對于運用無人機的非天線專業(yè)人員,為了方便他們對無人機載天線的效能進行定量評估,下面將干擾方向圖結(jié)果回歸于便于引用的解析公式.

機載天線干擾方向圖回歸的研究思路是:利用一體化建模及FEKO計算產(chǎn)生的數(shù)據(jù),結(jié)合極化匹配因子,采用IQGA,設(shè)計適應(yīng)度函數(shù),對方向圖中相關(guān)的回歸因子進行優(yōu)化,從而得到回歸公式,圖2為天線方向圖回歸的計算流程.

圖2　方向圖回歸的IQGA算法流程

3.1適應(yīng)度函數(shù)的設(shè)計

考慮到機載天線對地面目標的干擾方向圖既與x有關(guān)又與y有關(guān),是(x,y)的復(fù)雜函數(shù),設(shè)計適應(yīng)度函數(shù)為

(21)

式中: (xp,yq)為地面剖分點的坐標;g(xp,yq)為通過建模仿真獲得的干擾方向圖數(shù)值結(jié)果;Mx、My為研究區(qū)域剖分的總點數(shù);ai,j為待求系數(shù);n為回歸的最高次冪.

項數(shù)n的選取需綜合考慮以下兩個因素:

1) 回歸結(jié)果的精確度;

2) 工程應(yīng)用中的繁簡程度.

一般而言,n值越大精確度越高,但是其回歸公式也將隨項數(shù)越多而變得繁瑣,因此回歸中需兼顧上述兩方面的要求.

3.2干擾方向圖的回歸

為了不涉密,下面以公開出版文獻[19]中的對數(shù)周期天線為例,對于無人機常采用的盤旋飛行,取滾轉(zhuǎn)角為17°,設(shè)飛行高度為3 km,地面天線取垂直極化天線aR=(0,0,1)T,而對于水平極化天線aR=(0,1,0)T,由于機載天線傾斜45°安裝,因此其干擾方向圖與垂直極化天線aR=(0,0,1)T的干擾方向圖相差無幾,此處不再贅述.

利用仿真結(jié)果以及適應(yīng)度函數(shù)式(21),通過編程計算,可得歸一化的干擾方向圖如圖3所示,圖中以無人機當前所在位置為坐標中心點,無人機機身方向為x軸,圖中等值線數(shù)字為方向圖函數(shù)的分貝數(shù),回歸公式為

(22)

(a) FEKO仿真　　　　　　　　　　　(b)回歸公式計算 圖3　滾轉(zhuǎn)姿態(tài)時的機載天線干擾方向圖

i,j01234560-0.4430.106 4.12e-2-1.5e-21.22e-35.69e-5-7.01e-610.276 6.01e-2-3.61e-23.34e-3-1.59e-43.38e-6020.106-2.19e-2 6.18e-3-3.09e-4 3.87e-6003-4.0e-2 2.59e-3-3.84e-4 9.93e-60004 4.73e-3-1.33e-4 8.02e-600005-2.45e-4 2.49e-6000006 4.78e-6000000

由圖3可見,公式(22)在主要作用區(qū)域較好地計算了干擾功率密度方向圖函數(shù),但在邊緣區(qū)域精度較差,因此，公式(22)的適用范圍是

(23)

假如要增大回歸公式的適用范圍,可以采用增加回歸公式冪次的方法,或者采用分割地面區(qū)域使用不同回歸系數(shù)的方法.

圖3表明,對于常采用的盤旋飛行,對于地面垂直極化天線,最大干擾功率密度出現(xiàn)在無人機前方4.2 km附近,偏離航線1.5 km左右,3 dB作用區(qū)域范圍是1.7 km

為了清楚顯示回歸結(jié)果與原結(jié)果的差異,圖3和式(22)都是歸一化的.通過仿真計算,最大的g因子為gmax=0.098,因此,由式(22)乘以0.098即構(gòu)成干擾方向圖函數(shù)

(24)

4結(jié)論

本文提出了一種基于IQGA的機載天線干擾方向圖回歸算法,構(gòu)建了方向圖回歸適應(yīng)度函數(shù),仿真結(jié)果表明,所得到的回歸公式較精確地表達了機載天線的干擾方向圖.

回歸公式具有普適性,針對不同的無人機和天線只是式中系數(shù)ai,j具體的數(shù)值不同.該解析公式為無人機運用人員進行效能評估提供了便利實用的計算公式.

參考文獻

[1]沈文亮, 張卓鴻. 無人機在電子對抗中的應(yīng)用研究[J]. 艦船電子對抗, 2013, 36(6): 14-18.

SHENWenliang,ZHANGZhuohong.ResearchintotheapplicationofUAVtoelectroniccountermeasure[J].ShipboardElectronicCountermeasure, 2013, 36 (6): 14-18. (inChinese)

[2]葉瑞芳, 武坦然, 任翔宇. 外軍電子戰(zhàn)無人機發(fā)展前景[J]. 航天電子對抗, 2013, 29(2):12-15.

YERuifang,WUTanran,RENXiangyu.ThefuturedevelopmentofforeignmilitaryelectronicwarfareUVA[J].AerospaceElectronicWarfare, 2013, 29(2): 12-15. (inChinese)

[3]雷繼兆, 梁昌洪, 丁偉, 等. 機載天線輻射特性的并行FDTD分析[J]. 電波科學(xué)學(xué)報, 2008, 23(6): 1139-1143.

LEIJizhao,LIANGChanghong,DINGWei,etal.AnalysisofradiationcharactersofairborneantennaswithparallelFDTD[J].ChineseJournalofRadioScience, 2008, 23(6): 1139-1143. (inChinese)

[4]高軍, 曹祥玉, 劉濤. 機載通信天線輻射特性分析[J]. 西安電子科技大學(xué)學(xué)報, 2008, 35(4): 754-768.

GAOJun,CAOXiangyu,LIUTao.Analysisoftheradiationcharacteristicoftheairbornecommunicationantenna[J].JournalofXidianUniversity, 2008, 35(4): 754-768. (inChinese)

[5]LIUZL,WANGX,WANGCF.Installedperformancemodelingofcomplexantennaarraymountedonextremelylarge-scaleplatformusingfastMoM-POhybridframework[J].IEEETransonAntennasandPropagation, 2014, 62(7): 3852-3858.

[6]NEVEUN,HONGYK,LEEJJ,etal.Miniaturehexaferriteaxial-modehelicalantennaforunmannedaerialvehicleapplications[J].IEEETransonMagnetics, 2013, 49(7): 4265-4268.

[7]雷繼兆, 梁昌洪, 丁偉, 等. 并行FDTD-UTD方法分析機載相控陣天線[J]. 電波科學(xué)學(xué)報, 2009, 24(1): 60-64.

LEIJizhao,LIANGChanghong,DINGWei,etal.Analysisofairbornephased-arrayantennasusinghybridmethodofparallelFDTDandUTD[J].ChineseJournalofRadioScience, 2009, 24(1): 60-64. (inChinese)

[8]劉宏新, 徐高偉, 孟永超.CATIA三維設(shè)計基礎(chǔ)與應(yīng)用[M]. 北京: 機械工業(yè)出版社, 2014.

[9]范麗思, 崔耀中.FEKO5.4實例教程[M]. 北京: 國防工業(yè)出版社, 2012.

[10]STUTAMANWL,THIELEGA.AntennatheoryandDesign[M]. 2nded.NewYork:JohnWiley&SonsInc, 1998: 397-400.

[11]石力, 陳鑫, 吳瑋琦. 基于自適應(yīng)遺傳算法的方向圖綜合[J]. 電波科學(xué)學(xué)報, 2014, 29(1): 169-177.

SHILi,CHENXin,WUWeiqi.Patternsynthesisofantennabasedonanadaptivegeneticalgorithm[J].ChineseJournalofRadioScience, 2014, 29(1): 169-177. (inChinese)

[12]BHARGAVA,GUPTAN.Multiobjectivegeneticoptimizationofnonuniformlineararraywithlowsidelobesandbeamwidth[J].IEEEAntennasandWirelessPropagationLetters, 2013, 12: 1547-1549.

[13]KOPEREM,WOODWD,SCHNEIDERSW.Aircraftantennacouplingminimizationusinggeneticalgorithmsandapproximations[J].IEEETransonAerospaceandElectronicSystems, 2004, 40(2):742-751.

[14]孫思揚, 呂英華, 張金玲, 等. 基于遺傳算法的超寬帶微帶天線優(yōu)化設(shè)計[J]. 電波科學(xué)學(xué)報, 2011, 26(1): 62-66.

SUNSiyang,LüYinghua,ZHANGJinling,etal.DesignandoptimizationofUWBantennabasedongeneticalgorithm[J].ChineseJournalofRadioScience, 2011, 26(1): 62-66. (inChinese)

[15]DEBA,ROYJS,GUPTAB.Performancecomparisonofdifferentialevolution,particleswarmoptimizationandgeneticalgorithminthedesignofcircularlypolarizedmicrostripantennas[J].IEEETransonAntennasandPropagation, 2014, 62(8): 3920-3928.

[16]康行健. 天線原理與設(shè)計[M]. 北京: 北京理工大學(xué)出版社, 1993.

[17]束咸榮, 何炳發(fā). 論天線收發(fā)方向圖的非互易性[J]. 電波科學(xué)學(xué)報, 2009, 24(2): 293-296.

SHUXianrong,HEBingfa.Onnon-reciprocitybetweenreceivedandtransmittingpatternsofantenna[J].ChineseJournalofRadioScience, 2009, 24(2): 293-296. (inChinese)

[18]張建華, 周曉暉, 蔣鵬. 機載天線極化失配及方向性的分析計算[J]. 電波科學(xué)學(xué)報, 2007, 22(3): 497-501.

ZHANGJianhua,ZHOUXiaohui,JIANGPeng.Analysisofpolarizationmismatchingandradiationpatternofaircraftantenna[J].ChineseJournalofRadioScience, 2007, 22(3): 497-501. (inChinese)

[19]宋錚, 張建華, 黃冶. 天線與電波傳播[M]. 2版. 西安: 西安電子科大出版社, 2011: 125-131.

張建華(1962-),女,江蘇人,教授,研究方向為電磁場理論與天線技術(shù).

李振亞(1989-),男,河南人,碩士研究生,研究方向為天線技術(shù).

王玉文, 房艷燕, 董志偉, 等. 太赫茲波沿大氣層傾斜路徑的傳輸衰減[J]. 電波科學(xué)學(xué)報,2015,30(4):783-788. doi: 10.13443/j.cjors. 2014090804

WANG Yuwen, FANG Yanyan, DONG Zhiwei, et al. Terahertz transmission attenuation along the atmospheric slant path [J]. Chinese Journal of Radio Science,2015,30(4):783-788. (in Chinese). doi: 10.13443/j.cjors.2014090804