孟則宇，杜 彪,2，解 磊

(1.中國(guó)電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081；2.射電天文技術(shù)聯(lián)合實(shí)驗(yàn)室，河北 石家莊 050081)

基于CMA-ES優(yōu)化方法的波紋喇叭優(yōu)化技術(shù)

孟則宇1，杜 彪1,2，解 磊1

(1.中國(guó)電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081；2.射電天文技術(shù)聯(lián)合實(shí)驗(yàn)室，河北 石家莊 050081)

介紹了旋轉(zhuǎn)體時(shí)域有限差分法(BOR-FDTD)，導(dǎo)出了電磁場(chǎng)迭代計(jì)算公式。給出了自適應(yīng)協(xié)方差矩陣進(jìn)化策略(CMA-ES)的基本原理和步驟。提出了一種基于BOR-FDTD和CMA-ES的波紋喇叭優(yōu)化設(shè)計(jì)技術(shù)，并將該項(xiàng)技術(shù)用于平方公里陣(Square Kilometre Array，SKA)望遠(yuǎn)鏡天線Band 4(2.8～5.18 GHz)波紋喇叭饋源的優(yōu)化設(shè)計(jì)。計(jì)算結(jié)果表明，該饋源在工作頻帶內(nèi)反射損耗基本在-20 dB以下，天線口徑效率均優(yōu)于86.5%，且口徑效率隨頻率的變化較小。

SKA；反射面天線;波紋喇叭；BOR-FDTD；CMA-ES

0 引言

作為一種快速高效的數(shù)值方法，模式匹配法[1]已廣泛應(yīng)用于波紋喇叭[2]的分析和設(shè)計(jì)。然而大張角波紋喇叭的空間場(chǎng)和內(nèi)場(chǎng)之間存在較強(qiáng)的耦合效應(yīng)，采用此方法分析會(huì)引入較大的誤差。常規(guī)的三維全波分析方法(矩量法、有限元法和時(shí)域有限差分法等)能夠較為精確地分析大張角波紋喇叭的性能，但其計(jì)算量大，分析時(shí)間較長(zhǎng)，難于進(jìn)行優(yōu)化設(shè)計(jì)。針對(duì)以上問(wèn)題，本文采用BOR-FDTD方法[3]對(duì)大張角波紋喇叭的性能進(jìn)行了快速、高效和精確地分析，導(dǎo)出了電場(chǎng)和磁場(chǎng)的迭代計(jì)算公式。

CMA-ES是一種新型高效的全局優(yōu)化算法，其已經(jīng)在微帶貼片天線[11]、便攜式陶瓷諧振器天線[5]和平面超寬帶天線等[6]的優(yōu)化設(shè)計(jì)獲得了應(yīng)用，然而其在波紋喇叭優(yōu)化設(shè)計(jì)問(wèn)題中的應(yīng)用尚無(wú)人探討?；贐OR-FDTD分析方法和CMA-ES優(yōu)化方法，本文研究了波紋喇叭的快速優(yōu)化設(shè)計(jì)技術(shù)，并應(yīng)用于SKA天線[7]的 Band 4 波紋喇叭饋源的優(yōu)化設(shè)計(jì)。仿真結(jié)果表明，該饋源在工作頻帶內(nèi)反射損耗基本小于-20 dB，天線口徑效率均優(yōu)于86.5%，且在工作頻帶內(nèi)天線口徑效率變化較小。

1 算法基本原理

1.1 BOR-FDTD 的基本原理

波紋喇叭在結(jié)構(gòu)上具有旋轉(zhuǎn)對(duì)稱(chēng)性，因此可以采用BOR-FDTD對(duì)其性能進(jìn)行精確和高效的分析[8]。BOR-FDTD通過(guò)把電場(chǎng)和磁場(chǎng)進(jìn)行傅里葉展開(kāi)，利用時(shí)域有限差分法的基本原理求解波動(dòng)方程，將計(jì)算區(qū)域從三維空間轉(zhuǎn)換到一個(gè)二維的旋轉(zhuǎn)對(duì)稱(chēng)面上，可以極大地節(jié)省計(jì)算資源和計(jì)算時(shí)間。

利用待分析結(jié)構(gòu)在φ方向上的旋轉(zhuǎn)對(duì)稱(chēng)性，將電場(chǎng)和磁場(chǎng)進(jìn)行傅里葉展開(kāi)：

(1)

(2)

(3)

(4)

進(jìn)一步將式(3)和式(4)中的微分形式使用差分形式代替，可以得到差分形式的BOR-FDTD迭代公式。ρ方向的電場(chǎng)和磁場(chǎng)的迭代公式為：

(5)

(6)

式中，(i,j)為剖分網(wǎng)格的空間坐標(biāo)，上標(biāo)中的n代表時(shí)間步。同理可得φ和z方向上的電磁場(chǎng)迭代公式。

在波紋喇叭的分析中，端口采用主模TE11模式電場(chǎng)進(jìn)行激勵(lì)[9]，利用PML吸收邊界模擬無(wú)限空間的電磁場(chǎng)傳播，經(jīng)過(guò)迭代計(jì)算得到計(jì)算區(qū)域每個(gè)網(wǎng)格的電磁場(chǎng)分量的數(shù)值，繼而由近遠(yuǎn)場(chǎng)變換得到波紋喇叭的輻射方向圖[10]。

1.2 CMA-ES的基本原理

CMA-ES算法[11]與一般進(jìn)化算法相比其特點(diǎn)為：采用協(xié)方差矩陣自適應(yīng)地調(diào)整種群的分布，然后通過(guò)種群的分布情況利用多變量正態(tài)分布采點(diǎn)來(lái)產(chǎn)生算法種群，并且還引入了進(jìn)化路徑累積的概念。下面對(duì)其基本原理和步驟進(jìn)行簡(jiǎn)要介紹。

1.2.1 參數(shù)初始化

在優(yōu)化開(kāi)始前，首先需要設(shè)置種群的相關(guān)參數(shù)：均值m(0)，全局步長(zhǎng)σ(0)，協(xié)方差矩陣C(0)，子代個(gè)體數(shù)λ，父代個(gè)體數(shù)μ，進(jìn)化路徑Pσ(0)和Pc(0)，最大迭代次數(shù)，以及算法的一些自適應(yīng)參數(shù)。

1.2.2 種群采樣

采樣公式如下：

xk(g+1)=m(g)+σk(g)Nk(0,C(g))=

m(g)+σk(g)B(g)D(g)Nk(0,I)

(k=1,2,3,...,λ),

(7)

1.2.3 種群的評(píng)價(jià)與選擇

算法通過(guò)對(duì)個(gè)體的適應(yīng)度函數(shù)值進(jìn)行評(píng)估，進(jìn)行(μ,λ)截?cái)噙x擇，組成當(dāng)前最優(yōu)子群。

1.2.4 均值移動(dòng)

當(dāng)前的最優(yōu)子群通過(guò)加權(quán)重組得到新的分布均值。

(8)

式中,xi:λ(g+1)為g+1代中排名在第i位的個(gè)體。

1.2.5 協(xié)方差矩陣自適應(yīng)

先對(duì)協(xié)方差矩陣進(jìn)化路徑Pc更新：

m(g))/σ(g),

(9)

協(xié)方差矩陣C的調(diào)整公式如下：

C(g+1)=(1-c1-cμ)C(g)+c1(Pc(g+1)Pc(g+1)T+

(10)

式中，cc為Pc的更新學(xué)習(xí)速率，hσ為Heaviside函數(shù)，用于控制||Pc||的過(guò)大增長(zhǎng)，μeff為方差有效選擇質(zhì)量，且1≤μeff≤μ，c1和cμ分別為C的“秩1”和“秩μ”的更新學(xué)習(xí)速率。

1.2.6 步長(zhǎng)控制

步長(zhǎng)控制也稱(chēng)為累積步長(zhǎng)自適應(yīng)調(diào)整，首先對(duì)步長(zhǎng)進(jìn)化路徑Pσ的更新如下：

(11)

步長(zhǎng)σ的更新公式為：

(12)

式中，cσ為Pσ的更新學(xué)習(xí)率，dσ為接近于1的阻尼系數(shù)，E(||N(0,I)||)為歸一化進(jìn)化路徑在隨機(jī)選擇下的期望長(zhǎng)度。

在CMA-ES算法中以上步驟循環(huán)往復(fù)，引導(dǎo)種群收斂于全局最優(yōu)解。

2 基于BOR-FDTD和CMA-ES的波紋喇叭優(yōu)化設(shè)計(jì)

2.1 優(yōu)化設(shè)計(jì)流程

圖1給出了基于BOR-FDTD和CMA-ES的波紋喇叭優(yōu)化設(shè)計(jì)流程。首先對(duì)BOR-FDTD和CMA-ES進(jìn)行參數(shù)初始化，定義結(jié)構(gòu)變量，確定初始值。使用BOR-FDTD對(duì)建立的模型進(jìn)行分析，得到照射電平、等化誤差、交叉極化電平和反射損耗等性能參數(shù)，進(jìn)而得到目標(biāo)函數(shù)值。再利用CMA-ES對(duì)目標(biāo)函數(shù)進(jìn)行優(yōu)化，直到滿足指標(biāo)要求為止。

圖1 分析優(yōu)化設(shè)計(jì)流程圖

2.2 波紋喇叭的優(yōu)化模型

以SKA天線的band 4饋源為優(yōu)化實(shí)例，其性能要求為：工作頻帶為2.8～ 5.18 GHz，照射角為60°，邊緣照射電平控制在-12～ -17 dB，峰值交叉極化電平小于-20 dB，反射損耗小于-20 dB。

圖2為波紋喇叭結(jié)構(gòu)示意圖，在此模型的設(shè)計(jì)中主要分為3部分結(jié)構(gòu)：① 圓波導(dǎo)開(kāi)口處的臺(tái)階結(jié)構(gòu)，用于減少反射損耗，設(shè)置每個(gè)臺(tái)階長(zhǎng)度為b(i)，高度為v(i)；② 4個(gè)環(huán)加載形式的徑向槽結(jié)構(gòu)，每個(gè)槽設(shè)置槽深為d(i)，槽寬為w(i)；③ 在圓波導(dǎo)最外層的軸向槽結(jié)構(gòu)用于減少漏失損耗，其寬度為h，深度為a。在優(yōu)化中將這些結(jié)構(gòu)的尺寸作為優(yōu)化的參數(shù)，通過(guò)分析和優(yōu)化程序確定結(jié)構(gòu)的尺寸，以最終滿足所預(yù)期的優(yōu)化目標(biāo)。

圖2 波紋喇叭結(jié)構(gòu)

波紋喇叭的設(shè)計(jì)目標(biāo)是通過(guò)對(duì)其結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化設(shè)計(jì)，盡可能地提高天線效率ηf，然而若要精確計(jì)算SKA天線的效率費(fèi)時(shí)費(fèi)力，因此本文采用文獻(xiàn)[12]中給出的饋源效率因子計(jì)算公式來(lái)估算SKA天線的效率：

(13)

式中，ψ0為照射角，co(ψ)和xp(ψ)分別為波紋喇叭的45°面主極化和交叉極化輻射方向圖。

在優(yōu)化天線效率的同時(shí)要保證饋源工作頻帶內(nèi)反射系數(shù)低于RLmax，照射角內(nèi)峰值交叉極化電平(相對(duì)主極化最大值歸一)小于XPLmax。設(shè)M為工作頻帶內(nèi)頻率采樣點(diǎn)數(shù)，XPLm和RLm為第m個(gè)頻點(diǎn)的峰值交叉極化和端口反射系數(shù)，則該問(wèn)題可以抽象為以下目標(biāo)函數(shù)的最小化問(wèn)題：

f=1-ηf，

(14)

同時(shí)，應(yīng)滿足以下約束條件：

(15)

(16)

式中，RLmax和XPLmax均為-20 dB。采用精確罰函數(shù)法將原問(wèn)題轉(zhuǎn)化為無(wú)約束優(yōu)化問(wèn)題，從而該優(yōu)化問(wèn)題的目標(biāo)函數(shù)為：

(17)

式中，M為一個(gè)處理約束的大數(shù)，從而保證不滿足約束的解的目標(biāo)函數(shù)值一定大于滿足約束的解。

3 優(yōu)化設(shè)計(jì)結(jié)果

最終優(yōu)化得到的波紋喇叭開(kāi)口處采用3個(gè)臺(tái)階結(jié)構(gòu)，波紋槽數(shù)量為4個(gè)。如圖2所示，其長(zhǎng)度L相當(dāng)于2.8 GHz時(shí)波長(zhǎng)的1.16倍，開(kāi)口直徑D相當(dāng)于2.64λ，優(yōu)化出的饋源結(jié)構(gòu)緊湊。

圖3給出了優(yōu)化后的波紋喇叭在60°照射角內(nèi)的方向圖，可見(jiàn)喇叭輻射方向圖等化較好，交叉極化電平也能滿足指標(biāo)要求。

圖4為2.8～5.18GHz頻帶內(nèi)的反射損耗，在2.8GHz時(shí)反射損耗為18.3dB，在大部分工作頻帶內(nèi)反射損耗小于-20dB。

將優(yōu)化后的饋源帶入到SKA賦形雙偏置格里高利天線[13]中計(jì)算出的口徑效率如圖5所示?？梢?jiàn)工作頻帶內(nèi)天線口徑效率均大于86.5%，并且在工作頻帶內(nèi)天線效率波動(dòng)較小，滿足設(shè)計(jì)要求。

圖3 60°照射角內(nèi)方向圖

圖4 工作頻帶內(nèi)反射損耗

圖5 工作頻帶內(nèi)天線口徑效率

4 結(jié)束語(yǔ)

基于BOR-FDTD和CMA-ES提出了一種波紋喇叭優(yōu)化設(shè)計(jì)技術(shù)，并成功應(yīng)用于SKABand4 饋源的優(yōu)化設(shè)計(jì)。優(yōu)化結(jié)果表明，在工作頻帶內(nèi)反射損耗基本小于-20dB，天線口徑效率均優(yōu)于86.5%，并且口徑效率曲線波動(dòng)較小。證明了該優(yōu)化設(shè)計(jì)方法的正確性和有效性。本文所提出的優(yōu)化設(shè)計(jì)技術(shù)可廣泛用于其他旋轉(zhuǎn)對(duì)稱(chēng)結(jié)構(gòu)饋源喇叭的優(yōu)化設(shè)計(jì)。

[1]JamesG.DesignofWide-BandCompactCorrugatedHorns[J].IEEETransactionsonAntennas&Propagation,1984,32(10):1134-1138.

[2] 許德森，張文靜，秦順友.小口徑大張角波紋喇叭的優(yōu)良設(shè)計(jì) [J].無(wú)線電通信技術(shù)，2003，29(4)：49-50.

[3]DavidsonDB,ZiolkowskiRW.Body-of-revolutionFinite-differenceTime-domainModelingofSpace-timeFocusingbyaThree-dimensionalLens[J].IEEEAntennasandPropagation,1994,11(4): 1471-1490.

[4]GregoryMD,WernerDH.Multi-BandandWidebandAntennaDesignUsingPortSubstitutionandCMA-ES[C]∥2013IEEEAntennasandPropagationSocietyInternationalSymposium,2013: 598-599.

[5]FangXS,ChowCK,LeungKW,etal.NewSingle-/Dual-ModeDesignFormulasoftheRectangularDielectricResonatorAntennaUsingCovarianceMatrixAdaptationEvolutionaryStrategy[J].IEEEAntennasandWirelessPropagationLetters,2011,10:734-737.

[6]DaherEB,HoorfarA.NovelDesignofaPlanarUWBAntennawithEnhancedGainUsingaMixed-ParameterCMA-ESAlgorithm[C]∥RadioScienceMeeting,2014USNC-URSI:61-65.

[7] 杜 彪.大口徑反射面天線技術(shù)綜述[J].無(wú)線電通信技術(shù),2016,42(1):01-08.

[8]ChenY,MittraR,HarmsP.Finite-differenceTime-domainAlgorithmforSolvingMaxwell’sEquationsinRotationallySymmetricGeometries[J].IEEETransactionsonMicrowaveTheory&Techniques,1996,44(6):832-839.

[9]PratherDW,ShiS.FormulationandApplicationofTheFinite-differenceTime-domainMethodforTheAnalysisofAxiallySymmetricDiffractiveOpticalElements[J].IEEETrans.Microw.TheoryTech.,1999,16:1131-1142.

[10]FarahatN,YuW,MittraR.AFastNear-to-Far-FieldTransformationinBodyofRevolutionFinite-DifferenceTime-DomainMethod[J].IEEETransactionsonAntennas&Propagation,2003,51(9):2534-2540.

[11]GregoryMD,WernerDH.Multi-BandandWidebandAntennaDesignUsingPortSubstitutionandCMA-ES[C]∥IEEEInternationalSymposiumonAntennasandPropagation&Usnc/ursiNationalRadioScienceMeeting,2013:598-599.

[12]KildalPS.FactorizationoftheFeedEffiencyofParaboloidsandCassegrainAntenna[J].IEEETransactionsonAntennasandPropagation,1985,33(8):903-908.

[13]劉勝文，杜 彪，伍 洋.雙鏡天線的主副鏡擴(kuò)展研究與設(shè)計(jì)[J].無(wú)線電通信技術(shù),2015,42(1):61-64.

Design and Optimization of Corrugated Horn Based on CMA-ES

MENG Ze-yu1,DU Biao1,2，XIE Lei1

(1.The 54th Research Institute of CETC,Shijiazhuang Hebei 050081，China；2.Joint Laboratory for Radio Astronomy Technology ,NAOC & CETC54,Shijianzhuang Hebei 050081，China)

Body-of-revolution finite-difference time-domain (BOR-FDTD) algorithm is presented,and the formula for iterative computation of electromagnetic field is obtained.The principle and procedures of covariance matrix adaptation evolution strategy (CMA-ES) is introduced to optimize the structure parameters of the corrugated horn.An optimization method of corrugated horn based on BOR-FDTD and CMA-ES is proposed and is applied to the design of Band 4 feed for the SKA Dish.Optimization results show that the reflection loss is less than -20 dB over the operating frequency range,the aperture efficiency is better than 86.5%,and the efficiency variation with the frequency is small.

SKA;Dish;corrugated horn;BOR-FDTD;CMA-ES

10.3969/j.issn.1003-3114.2017.01.17

孟則宇,杜 彪,解 磊.基于CMA-ES優(yōu)化方法的波紋喇叭優(yōu)化技術(shù)[J].無(wú)線電通信技術(shù)，2017，43(1):68-72.

2016-09-23

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(2013CB837902)；國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(SS2014AA122001)；國(guó)家自然科學(xué)基金項(xiàng)目(11261140641)

孟則宇(1992—)，男，碩士研究生，主要研究方向：電磁場(chǎng)與微波技術(shù)。杜 彪(1962—)，男，博士，研究員，現(xiàn)任中國(guó)電子科技集團(tuán)公司第五十四研究所首席專(zhuān)家、所副總工程師，主要研究方向：射電望遠(yuǎn)鏡天線、衛(wèi)星通信地球站天線、微波天線、饋源系統(tǒng)和陣列天線等。

TN820

A

1003-3114(2017)01-68-5