星載稀疏相控陣天線的多目標優(yōu)化設(shè)計

劉 恒，趙宏偉，劉 波

（西安空間無線電技術(shù)研究所 陜西 西安　710100）

星載稀疏相控陣天線的多目標優(yōu)化設(shè)計

劉 恒，趙宏偉，劉 波

（西安空間無線電技術(shù)研究所 陜西 西安710100）

針對星載相控陣天線高增益、低副瓣的設(shè)計要求，提出了對星載六邊形相控陣天線進行稀疏布陣的方法。根據(jù)需求建立了星載相控陣天線高增益、低副瓣的兩目標的優(yōu)化模型，使用NASA-II算法對其進行多目標設(shè)計優(yōu)化。為了快速計算六邊形天線陣的方向圖，通過在孔徑中添加虛擬陣元把六邊形陣列轉(zhuǎn)化為可以實現(xiàn)二維傅里葉變換的矩形陣列，大大節(jié)省了算法優(yōu)化時間。結(jié)果表明，星載稀疏相控陣天線多目標優(yōu)化思路和方法在工程應(yīng)用中的有效性和合理性，為設(shè)計者提供了多種可行的選擇方案。

稀疏陣列；相控陣天線；FFT；NASA-II算法

星載稀疏相控陣天線是指從一個陣元周期分布的陣列中去除掉一些天線陣元，或者把這些陣元連接到匹配負載上。既可以減少陣列天線成本和重量，還可以獲得與滿陣排布相當?shù)恼ㄊ?，當陣元均勻激勵時，稀疏陣列天線可以獲得比滿陣布置更低的副瓣電平［1］。在過去幾十年中，稀疏相控陣天線的研究一直是一項重要和有意義的工作。由于在同等陣元間距下，六邊形相控陣天線在掃描時比矩形柵格陣列出現(xiàn)柵瓣的時間更晚；且在等幅激勵時近似圓形的六邊形平面相控陣能夠獲得相對低的副瓣電平，因此六邊形相控陣天線陣被廣泛應(yīng)用在衛(wèi)星通信中［2-3］。在星載稀疏相控陣天線設(shè)計中，希望主瓣增益盡可能的高而副瓣電平盡可能的低，同時最大副瓣電平和主瓣增益分別低于和高于某一給定的約束值。因此，星載稀疏相控陣天線的設(shè)計是典型的多目標優(yōu)化問題。

隨著計算機技術(shù)的發(fā)展，高效的稀疏陣列優(yōu)化算法已成研究熱點，主要有遺傳算法［4-5］、粒子群算法［6-7］、進化差分算法［8-9］、蟻群算法［10］等各種優(yōu)化算法的混合算法［11-12］。各種智能進化優(yōu)化算法。在這些文獻中通常把多個目標函數(shù)通過不同的加權(quán)權(quán)重構(gòu)造成單個目標函數(shù)進行優(yōu)化。在實際設(shè)計中，權(quán)重的分配及單目標函數(shù)構(gòu)造方式帶有很強的主觀性，而決策者往往需要多種可供選擇的方案進行比較，就需要不斷的修改權(quán)重及單目標函數(shù)進行很多次優(yōu)化。在多目標優(yōu)化中，各目標之間無法比較，存在沖突，導(dǎo)致不存在對所有目標函數(shù)都是最優(yōu)解。Goldberg將Pareto知識［13］引入進化算法中來解決多目標優(yōu)化問題的新思路之后，越來越多學(xué)者們提出了基于Pareto支配關(guān)系的進化算法，使之成為求解多目標優(yōu)化問題的一種趨勢。其中以非支配遺傳算法（NSGA-II）［14］最具有代表性。

本文應(yīng)用NSGA-II算法對衛(wèi)星上廣泛應(yīng)用的正三角形排列的六邊形平面相控陣天線進行稀疏布陣優(yōu)化，以達到高增益、低副瓣的輻射特性。首先利用六邊形陣列的幾何對稱陣性確定了六邊形平面相控陣天線中的獨立單元數(shù)，實現(xiàn)對大型六邊形相控陣天線優(yōu)化變量進行降維。并根據(jù)六邊形陣列的排列特征，通過在陣列中添加虛擬單元的方式把六邊形轉(zhuǎn)換為矩形陣列，從而利用FFT快速計算陣列方向圖。然后建立了增益和峰值副瓣電平兩目標優(yōu)化模型的解析表達式。利用NASA-II算法對1027個單元（L=18）的六邊形相控陣天線進行多目標優(yōu)化設(shè)計，以降低陣列天線的峰值副瓣電平、增高主瓣增益為目的，得到了最優(yōu)解集的Pareto前端，取得了滿意的結(jié)果。說明了星載稀疏相控陣天線多目標優(yōu)化思路和方法在工程應(yīng)用中的合理性和有效性，為設(shè)計者提供了多種可行的選擇方案。

1　星載相控陣天線模型

1.1六邊形天線方向圖

星載六邊形相控陣天線的幾何結(jié)構(gòu)如圖1所示，它是由按等邊三角形排列的天線單元組成?？梢钥闯隽呅纹矫嫦嗫仃囂炀€是由若干行橫向等間距的直線陣在縱向交替緊密排布組成，其每行的單元數(shù)量隨著遠離中心行依次遞減。對于一個由L個同心六邊形環(huán)陣組成六邊形天線，其第一行的單元數(shù)為L+1個；隨著行數(shù)的增加，每行的單元數(shù)依次增一個，到L+1行的單元數(shù)增加到最多的2L+1個；再隨著行數(shù)繼續(xù)增加，每行的單元數(shù)則依次減少一個，到2L+1行的單元數(shù)減少到L+1個。下面將利用六邊形平面天線陣的這一結(jié)構(gòu)特性推導(dǎo)它的輻射方向圖的解析表達式。

圖1　六邊形平面陣列天線結(jié)構(gòu)圖

圖2　可見空間區(qū)域示意圖

由于六邊形相控陣天線的每行在y方向是交替排布的，如果通過在每行單元中間添加一個虛擬“單元”，這樣就把六邊形平面陣列變換為矩形柵格的平面陣列。每行沿x軸方向間距dx=0.5d，每列沿y軸方向間距dy=0.5d·tan60°，d為正三角形排列的單元間距。根據(jù)六邊形平面陣列的結(jié)構(gòu)和添加的虛擬單元的位置可知，對于一個由L個同心六邊形環(huán)陣通過在每一行單元中間添加虛擬單元變化為一個2L+1行4L+1列的矩形柵格平面陣。由陣列天線理論可知，在不考慮互耦的情況下，輻射單元為理想點源，六邊形陣列方向圖滿足乘積定理可以表示為：

式中：L為同心六邊形環(huán)陣的數(shù)量，k為波數(shù)，θ，φ分別為球坐標系的下俯仰角和方位角，Amn是第（m，n）個陣元激勵幅度。為了方便，把由M×N個陣元激勵組成的矩陣記著A，對于矩陣A與六邊形平面陣列不匹配的點，如添加的虛擬陣元和位于六邊形陣列孔徑外的點，其激勵幅度Amn=0。對于與陣列匹配的點，在稀疏優(yōu)化過程中，當Amn=1，表示單元在工作；當Amn=0，表示單元被剔除。一個由L個同心六邊形環(huán)陣組成六邊形平面天線陣的單元總數(shù)可由下式確定：

由于大型六邊形陣列天線單元數(shù)成百上千，而智能進化算法在對優(yōu)化變量超過200個的時候，效果變的不太理想，因此需要對六邊形陣列天線獨立單元進行降維。根據(jù)對稱性，總共有6個對稱軸，位于對稱軸上的陣元有6個對稱陣元，而非對稱軸上的陣元有12個對稱陣元。在陣列稀疏優(yōu)化過程中如果某一個陣元被選定剔除，則與其對稱的所有的陣元都將剔除。如圖1所示位于圖中P標記和y軸的30°扇形區(qū)內(nèi)的單元為該天線陣的獨立單元。一個具有L個同心六邊形環(huán)陣的六邊形相控陣天線的獨立單元數(shù)可由下式計算，其中int（·）為取整函數(shù)。

1.2方向圖與FFT的關(guān)系

由于在進化過程中要反復(fù)計算方向圖函數(shù)的采樣點來評估個體的優(yōu)劣，如果能將降低計算方向圖的時間就能夠減少評估個體的時間，從而大大提高計算速度。設(shè)一個輻射陣元為理想點源，沿x軸方向的陣元間距為dx，y軸方向的陣元間距為dy的M×N元平面陣，在不考慮互耦的情況下，陣列方向圖滿足乘積定理可以表示為：

式中：Amn是第（m，n）個陣元激勵幅度，k=2π/λ，λ為自由空間波長，u=sinθcosφ，v=sinθsinφ是方向余弦，θ，φ分別為球坐標系的下俯仰角和方位角。令p=Mkdxu/2π+1，q=Nkdyv/2π+ 1，則式（4）變換為

由式 （4）可以看出陣因子AF與陣元激勵A(yù)之間存在IFFT關(guān)系A(chǔ)F=MN×IFFT（A），這樣根據(jù)具體的問題，將（θ，φ）域的方向圖特性映射到（p，q）域中，就可以在（p，q）進行優(yōu)化。

由u=sinθcosφ，v=sinθsinφ可得u2+v2=sin2θ≤1，在u，v坐標系下，可見區(qū)域與一圓區(qū)域相對應(yīng)。對于不同間距的平面陣列，利用二維FFT計算出方向圖的可見區(qū)域如圖2所示。當dx≤λ/2，dy≤λ/2時，可見區(qū)域為半徑為1的圓A；當dx＞λ/2，dy＜λ/2時，可見區(qū)域為半徑為1的圓與矩形D的交集；當dx＜λ/2，dy＞λ/2時，可見區(qū)域為半徑為1的圓與矩形C的交集；當dx＞λ/2，dy＞λ/2時，可見區(qū)域為矩形C與矩形D的交集。

1.3目標函數(shù)及約束

星載相控陣天線為了得到高增益低副瓣的輻射特性，稀疏布陣綜合優(yōu)化的兩個目標函數(shù)為最大副瓣電平MSLL與主瓣點的增益G。

其中，x為待優(yōu)化的相控陣天線獨立單元的狀態(tài)，η為天線的效率，MSLL和D分別為陣列結(jié)構(gòu)x計算得到的最大副瓣電平和最大方向系數(shù)。u0=sinθ0cosφ0，v0=sinθ0sinφ0是主瓣點（θ0，φ0）的方向余弦。二重積分的可視空間為u2+v2≤1?？紤]到實際需求，通常最大副瓣電平低于SLLmax=-20 dB，所以約束條件g（x）=MSLL-SLLmax＜0。

2　多目標優(yōu)化設(shè)計

對于帶約束的多目標優(yōu)化問題，其數(shù)學(xué)上可以表示為：

其中F（x）為目標函數(shù)，g（x）為約束函數(shù)，n和k分別為目標函數(shù)和約束函數(shù)的個數(shù)。多目標優(yōu)化問題是一個解的集合，這些解在全部目標函數(shù)而言是無法比較優(yōu)劣的，其特點是無法在改進任何目標函數(shù)的同時不削弱至少一個其他目標函數(shù)，這些解的集合稱為Pareto非支配解集合。

2.1Pareto最優(yōu)解排序

在單目標優(yōu)化算法中，兩個解的優(yōu)劣程度可以通過他們的唯一適應(yīng)度值大小確定。在多目標優(yōu)化算法中，所有的目標函數(shù)是同樣的重要的，即在優(yōu)化過程中不能為了改善某些目標函數(shù)而惡化其他任何一個目標函數(shù)。設(shè)x1，x2是所描述的多目標問題的解，當下列條件成立時稱x1限制性支配x2，也稱與x2相比，x1是Pareto占優(yōu)。

1）x1屬于可行解而不屬于x2。

2）x1和x2都不屬于可行解且在約束函數(shù) gi（x1）≤gi（x2）（i=1，2，…，k）。

3）x1和 x2都屬于可行解且在目標函數(shù)Fi（x1）≤Fi（x2）（i=1，2，…，m）。

對于多目標優(yōu)化的非支配解，經(jīng)過非支配排序和擁擠度計算［14］，種群中的每個個體i都有兩個屬性；非支配序 irank和擁擠度idistance。判斷種群中任意兩個個體優(yōu)劣的標準為，如果兩個個體非支配序不同，非支配序小的更優(yōu)；如果兩個個體的非支配序相同，則擁擠度較大的個體更優(yōu)。

2.2NSGA-II算法

將進化算法與Pareto概論結(jié)合，發(fā)展出很多基于Pareto最優(yōu)概論的多目標進化算法，其中具有代表性的是NSGAII。由于NSGA-II在對稀疏相控陣天線進行優(yōu)化的時可以直接進行二進制編碼，不需要進行編碼處理，是稀疏相控陣天線多目標優(yōu)化的最佳工具。本文直接采用的文獻［14］的二進制編碼NSGA-II算法。由于多目標優(yōu)化中非支配解的優(yōu)劣無法判斷，因此進化后的種群規(guī)模會增加，經(jīng)過交叉、變異后的NP個子代與NP個父代混合進行非支配排序，種群規(guī)模在進化后增加為2NP，然后通過非支配排和擁擠度排序后選擇前NP個個體進入下一代。根據(jù)以上策略，NSGA-II算法步驟如下：

1）初始化種群，生成種群規(guī)模為NP個父代。

2）計算種群中每個個體的目標函數(shù)值和約束函數(shù)值。

3）通過遺傳算法的交叉、變異生成群規(guī)模為NP個子代。

4）計算子代的目標函數(shù)及約束函數(shù)。

5）對父代和子代混合的2NP個個體進行非支配排序，選出前NP個性能最佳的個體。

6）從步驟3）重復(fù)直至滿足迭代結(jié)束條件。

3　實例分析

在實際應(yīng)用中，經(jīng)過計算和分析，我們?nèi)》N群規(guī)模NP等于二進制序列長度的二倍。陣元間距d=λ/2，二維IFFT運算采樣點數(shù)K×K=512×512。NSGA-II算法參數(shù)的選擇使用文獻［15］推薦的。單點交叉，交叉概率取0.9，變異概率取0.01，進化代數(shù)取5000次。假設(shè)天線的效率η=1，則增益G等于方向性系數(shù)D，以由1027個單元（L=18）的六邊形相控陣天線為例進行優(yōu)化實例。

圖3　Pareto最優(yōu)解集前端分布

圖3給出了在5 000次兩目標函數(shù)評價后，1 027個單元的組成六邊形平面相控陣天線稀疏布陣的Pareto最優(yōu)解集前端分布?？梢钥闯龇逯蹈卑觌娖较孪逓?（-27.8 dB，33.2 dB），約束的上限為（-20 dB，34.1 dB）。圖4和圖5分別給出了峰值副瓣電平上限和下線的孔徑陣元分布和歸一化遠場方向圖，其可視空間為u2+v2≤1?？梢钥闯?，優(yōu)化得到的六邊形稀疏陣列符合陣列優(yōu)化規(guī)律，即在優(yōu)化陣列中，陣元的稀疏總是發(fā)生在陣列邊緣，而陣列中心的陣元一般不會被稀疏掉。

表給出了從求解得到的Pareto最優(yōu)解中選擇4個最優(yōu)解對應(yīng)的變量值以及對應(yīng)的目標函數(shù)值。其中包括兩個邊界點和兩個中間點。表格中的數(shù)據(jù)表明，峰值副瓣電平PSLL最低為-27.8 dB，方向性系數(shù)達到最小值33.2 dB，而當PSLL最高為-20 dB，方向性系數(shù)達到最大值34.1 dB，這兩個目標函數(shù)不能同時取得最優(yōu)值。說明在星載相控陣天線優(yōu)化設(shè)計中，最大副瓣電平和增益確實是一組存在矛盾的目標，適合用多目標優(yōu)化方法求解這類問題。相對于單目標優(yōu)化算法，采用多目標方法能從多個不同的角度提供不同的最優(yōu)解，為決策者提供更多的指導(dǎo)信息和幫助。

圖4　方向性系數(shù)最好的解

圖5　峰值副瓣電平最好的解

表1　最優(yōu)集合中4個最優(yōu)解對應(yīng)的優(yōu)化結(jié)果

4　結(jié)　論

文中以星載稀疏相控陣天線優(yōu)化布陣為研究對象，采用多目標優(yōu)化算法，將陣列的優(yōu)化布陣轉(zhuǎn)化為多目標優(yōu)化的數(shù)學(xué)問題。為了降低大型六邊形陣列的優(yōu)化變量維數(shù)，通過對陣列中對稱的單元進行等效處理，獨立單元數(shù)降低為陣列單元的十分之一。同時把六邊形陣列轉(zhuǎn)化為可以實現(xiàn)二維傅里葉變換的矩形陣列模型，使計算量成倍減少，也可用于解決其他種類的矩形陣和三角陣平面陣列天線的稀疏布陣優(yōu)化，具有很好工程實用價值，為大型陣列天線的稀疏優(yōu)化提供了有效途徑。星載天線的高增益、低副瓣為優(yōu)化目標，充分體現(xiàn)了星載稀疏相控陣天線多目標優(yōu)化思路和方法的合理性和有效性，為決策者提供了多種可行的選擇方案，未來可考慮更多的目標及約束進行進一步的應(yīng)用研究。

［1］W.P.M.N.Keizer.Synthesis of Thinned planar circular and square Arrays using density tapering［J］.IEEE Transactions on Antennas and Propagation，2014，62（4）：1555-1563.

［2］Dietrich F J，Metzen P，Monte P.The Globalstar Cellular satellite system［J］.IEEE trans on Antennas，1998，46（6）：935-942.

［3］李東風(fēng)，龔中麟.六邊形平面天線陣優(yōu)化稀疏布陣研究［J］.電子學(xué)報，2002，30（3）：376-380.

［4］Haupt R L.Thinned arrays using genetic algorithms［J］.IEEE Transactions on Antennas and Propagation，1994，42（7）：993-999.

［5］王玲玲，方大綱.運用遺傳算法綜合稀疏陣列［J］.電子學(xué)報，2003，31（12A）：2135-2138.

［6］N.Jin，Y.Rahmat-Samii.Advances in particle swarm optimization for antenna design：real-number，binary，singleobjective and multiple objective implementations［J］.IEEE Trans on Antennas and Propagation，2007（55）：556-567.

［7］M.Donelli，A.Martini，A.Massa.A hybrid approach based on PSO and Hadamard difference sets for the synthesis of square thinned arrays［J］.IEEE Trans on Antennas and Propagation，2009，57：2491-2495.

［8］D.G.Kurup，M.Himdi，A.Rydbery，Synthesis of uniform amplitude unequally spaced arrays using the differential evolution algorithm［J］.IEEE Trans on Antennas and Propagation，2003，51：2210-2217.

［9］L.Zhang，Y C Jiao，Z B Weng，F(xiàn) S Zhang.Design of planar thinnedarraysusingaBooleandifferentialevolution algorithm［J］.IET Microwaves，Antennas&Propagation，2010，4（12）：2172-2178.

［10］ó.Quevedo-Teruel，E.Rajo-Iglesias，Ant colony optimization in thinned array synthesis with minimum sidelobe level，IEEE Antennas Wireless Propagation Letters，2006，5：349-352.

［11］Leonardo F.Yepes，David H.Convarrubias，Miguel A. Alonso，et al.Hybrid Sparse Linear Array Synthesis Applied to Phased Antenna Arrays［J］.IEEE Antenna and Wireless Propagation Letters，2014（13）：185-188.

［12］OvidioMarioBucci，StefanoPerna，DanielePinchera. Synthesis of Isophoric Sparse Array Allowing Zoomable Beams and Arbitrary Coverage in Satellite Communications ［J］.IEEETransactionsonAntennasandPropagation，2015，63（4）：1445-1457.

［13］Goldberg D E.Genetic Algorithm in Search，Optimization and Machine Learning.Boston，MA，USA：Addison-Wesley Longman Publishing Co.，Inc.，1989.

［14］K.Deb，A.Pratap，S.Agarwal，T.Meyarivan，A fast and elitist multiobjective genetic algorithm：NSGA-II［J］.IEEE trans on Evol.Comput.，2002，6（2）：182-197.

［15］Pareto Optimization Microwave Filter Design Using Multiobjective Differential Evolution［J］.IEEE Trans.Antenna propagation，2010，58（1）：132-142.

Multi-objective design optimization of satellite thinned array antenna

LIU Heng，ZHAO Hong-wei，LIU Bo
（China Academy of Space Technology（Xi'an），Xi'an 710100，China）

Considering the design purpose of high gain，low side lobes satellite phased array antennas，a thinning process is proposed for hexagonal array.The two objective optimization modes of high gain and low side lobes of satellite antenna are set up according to the requirements.And the NASA-II is applied at this two objective design optimization.Virtual elements are added in the array for translating the hexagon to rectangle which could be convenient to calculate the array factor of hexagon by using the FFT.The simulated results show that the availability and rationality of the multi objective optimization in the satellite thinned phased array antenna，and it could offer several feasible scheme for the designer.

thinned array；phased array antenna；fast fourier transform algorithm；NASA-algorithm

TN959

A

1674－6236（2016）09-0122-04

2015-11-23稿件編號：201511212

國家自然科學(xué)基金（61201089）；國家重點實驗室基金（9140C530101130C53013）

劉 恒（1986—），男，湖南耒陽人，博士。研究方向：陣列天線優(yōu)化與設(shè)計。