基于凸優(yōu)化的低復(fù)雜度平面陣列綜合方法

陳金立，鐘濟(jì)陽，曹華松

（南京信息工程大學(xué) 電子與信息工程學(xué)院，江蘇 南京　210044）

?

基于凸優(yōu)化的低復(fù)雜度平面陣列綜合方法

陳金立，鐘濟(jì)陽，曹華松

（南京信息工程大學(xué) 電子與信息工程學(xué)院，江蘇 南京210044）

摘要：針對運用凸優(yōu)化方法綜合平面陣列耗時較長的問題，提出一種基于凸優(yōu)化的低復(fù)雜度平面陣列綜合方法。該方法將基于加權(quán)?1范數(shù)最小化的平面陣列綜合模型中二維波束旁瓣約束條件轉(zhuǎn)換成橫向陣列和縱向陣列的一維波束旁瓣約束問題，從而降低了凸優(yōu)化的約束條件數(shù)，有效加快了算法的收斂速度。仿真結(jié)果表明，該方法雖然在平面陣列綜合的稀疏程度方面要稍劣于傳統(tǒng)凸優(yōu)化方法，但是它能快速獲得具有大范圍低波束旁瓣電平以及更短孔徑的平面陣列。

關(guān)鍵詞：平面陣列；陣列綜合；凸優(yōu)化；低復(fù)雜度

0　引言

稀疏陣列是一種從規(guī)則的柵格中抽去天線單元或接匹配負(fù)載而形成的陣列，其具有孔徑較大、波束較窄以及空間分辨率較高等性能，已被廣泛應(yīng)用于雷達(dá)、電子通信、以及衛(wèi)星廣播電視等領(lǐng)域[1]。然而，稀疏陣列由于其陣元的稀疏布置導(dǎo)致其波束旁瓣電平抬高。稀疏陣列綜合的目的就是通過聯(lián)合優(yōu)化陣列的陣元位置和激勵方式，使得天線陣列能夠以最少的陣元數(shù)目滿足期望的波束輻射特性[2]。

隨著計算機(jī)技術(shù)的飛速發(fā)展，高效的陣列綜合方法已成為研究熱點。目前應(yīng)用于平面稀疏陣列綜合的算法主要有遺傳算法[3]、模擬退火算法[4]、粒子群算法[5]以及蟻群算法[6]等，這些算法從本質(zhì)上來說都是基于隨機(jī)性的自然算法，在綜合過程中需要進(jìn)行大規(guī)模的搜索，導(dǎo)致此類算法的收斂速度慢，并容易陷入局部最優(yōu)解等問題。稀疏陣列在空間上離散、稀疏的分布特點與最近發(fā)展的信號重構(gòu)理論中稀疏信號的特性相類似，可將天線陣列綜合問題看作是空間稀疏信號重構(gòu)的問題，因此，稀疏信號重構(gòu)技術(shù)為陣列綜合問題的解決提供了一個新的途徑。文獻(xiàn)[7]提出了一種基于?1范數(shù)最小化的稀疏線性陣列綜合方法，該方法將稀疏信號重構(gòu)方法如凸優(yōu)化方法應(yīng)用于稀疏陣列綜合問題的求解中。文獻(xiàn)[8]建立了基于迭代加權(quán)?1范數(shù)最小化的陣列綜合模型，利用迭代凸優(yōu)化方法實現(xiàn)線性陣列綜合，即利用上一次迭代獲得的陣列激勵向量求解當(dāng)前迭代的激勵向量，從而能獲得稀疏程度更高的稀疏陣列。

大型平面陣列可以同時測量目標(biāo)的俯仰角和方位角，從而能夠?qū)崿F(xiàn)目標(biāo)的三維空間定位，其應(yīng)用范圍要大于遠(yuǎn)線性陣列，因此研究平面陣列的綜合方法更具有實際意義。文獻(xiàn)[9]利用加權(quán)矩陣的共軛對稱性將波束約束條件轉(zhuǎn)變?yōu)橥购瘮?shù)，從而可利用凸優(yōu)化方法解決平面陣列的綜合問題。文獻(xiàn)[10]利用凸優(yōu)化方法求解一系列基于加權(quán)?1范數(shù)最小化的平面陣列波束綜合問題，在滿足輻射特性的條件下能獲得盡可能稀疏的平面陣列。由于平面陣列通常包含較多的陣元，而且在陣列綜合模型中其二維波束旁瓣的抑制范圍較大，從而導(dǎo)致需要相當(dāng)長的時間來綜合出滿足要求的平面陣列。為此，本文提出了一種基于凸優(yōu)化的低復(fù)雜度平面稀疏陣列綜合方法。該方法首先建立基于加權(quán)?1范數(shù)最小化的平面陣列綜合問題，然后利用Kronecker積的性質(zhì)將整個二維觀測角度上的波形約束條件轉(zhuǎn)換為兩個一維觀測角度上的波形約束問題，從而使得平面陣列綜合的復(fù)雜度與兩個線性陣列綜合問題相當(dāng)。因此，本文方法在滿足期望峰值波束旁瓣電平的條件下能夠快速綜合出孔徑更短的平面陣列，由于該方法綜合后的平面陣列具有大范圍低波束旁瓣電平，使得其在平面陣列的綜合稀疏程度方面會稍劣于傳統(tǒng)凸優(yōu)化方法。

1　平面稀疏陣列模型

考慮如圖1所示的平面陣列，圖中黑點表示天線單元。在陣列區(qū)域內(nèi)，沿x軸及其平行方向各有M個柵格點，沿y軸及其平行方向各有N個柵格點，這些柵格點可布置天線單元，dx和dy分別表示x軸方向和y軸方向的相鄰陣元最小間距。θ和φ分別表示信號的俯仰角和方位角，定義域分別為

圖1　平面稀疏陣列結(jié)構(gòu)

平面陣列的導(dǎo)向矢量A和激勵矢量w分別表示如下：

式中：[?]T表示矢量轉(zhuǎn)置；amn=e-jk(mdxμ+ndyv)是陣列單元(m,n)的方向圖 (m=0,1,2,…,M-1;n=0,1,2,…,N-1)；，λ為波長；μ和 v是方向參數(shù)，分別表示為μ=sin θ cos φ和 v=sin θ sin φ，其取值范圍分別為μ∈[-1,1]和v∈[-1,1]；wmn是第(m,n)個位置上天線單元的激勵值。假設(shè)陣列中的單元都是無方向性的，則平面陣列的方向圖可表示為：

由式（3）可知，影響平面陣列波束方向圖綜合性能的主要因素有陣元的個數(shù)、陣元的位置及其激勵。稀疏陣列綜合的本質(zhì)就是在給定的約束條件下求取上述各參數(shù)的過程。在給定的平面陣列規(guī)模和旁瓣電平約束條件下，以尋找出陣元數(shù)目最少的平面稀疏陣列為目標(biāo)對陣元位置和陣元激勵進(jìn)行綜合，則采用加權(quán)?1范數(shù)優(yōu)化的陣列模型表示如下[10]：

式中：[?]H表示矢量的共軛轉(zhuǎn)置；Q是由一組加權(quán)系數(shù)qk組成的1?范數(shù)加權(quán)矩陣；μ0和 v0為期望波束主瓣所對應(yīng)的方向參數(shù)；Ω為平面陣列的波束旁瓣區(qū)域；psl是給定的峰值旁瓣電平。

2　平面陣列的低復(fù)雜度綜合方法

由式（4）可知，平面陣列綜合問題的求解所需時間與平面陣列的旁瓣區(qū)域Ω有關(guān)，與線性陣列相比，由于平面陣列的觀測角度為二維角度，因此其在旁瓣區(qū)域Ω內(nèi)的采樣數(shù)會呈平方式增長，從而導(dǎo)致陣列綜合問題中約束條件較多，則利用凸優(yōu)化方法求解式（4）會耗時很長。因此為了提高平面陣列綜合的速度，必須簡化式（4）中的波束旁瓣約束條件。

由式（1）可知，平面陣列的導(dǎo)向矢量A可改寫為：

式中：wx和wy分別為平面陣列的橫向陣列和縱向陣列的激勵向量。將式（5）和式（6）代入式（3），根據(jù)Kro?necker積的性質(zhì)可得：

式中：Fx(μ)=(μ)×wx為橫向陣列的波束方向圖；Fy(v)=(v)×wy為縱向陣列的波束方向圖。假設(shè)Fx(μ) 和Fy(v)均為歸一化波束方向圖，即Fx(μ)1，F(xiàn)y(v)1，則根據(jù)式（7），式（4b）可改寫為：

旁瓣電平約束條件式（4c）可等效為：

式中：Ωx和Ωy分別為方向參數(shù) μ和v在波束旁瓣區(qū)域內(nèi)的取值范圍。由于Fx(μ)1，F(xiàn)y(v)1，滿足式（9）的充分條件為：

由式（8）~式（10）可知，基于加權(quán)?1范數(shù)最小化的平面陣列綜合問題，即式（4）可改寫為如下形式：

式中Qx和Qy分別為橫向陣列和縱向陣列的1?范數(shù)加權(quán)矩陣，分別由式（12）和式（13）決定。

由式（11）可知，本文利用Kronecker積的性質(zhì)將整個二維觀測角度上的旁瓣約束條件即式（4）分解成兩個一維觀測角度上的旁瓣約束問題，大大降低了平面陣列綜合時的約束條件數(shù)，即平面陣列綜合問題等效于橫向陣列和縱向陣列的綜合問題，從而使得平面陣列綜合的復(fù)雜度與兩個線性陣列綜合問題相當(dāng)。當(dāng)?shù)螖?shù)達(dá)到或者滿足終止條件（14）和（15）時平面陣列綜合完成：

式中ξx和ξy為誤差值。令：

3　仿真實驗

本文對基于凸優(yōu)化的平面陣列綜合問題中的約束條件進(jìn)行了簡化，提出一種低復(fù)雜度的平面陣列綜合方法。下面設(shè)計了運用傳統(tǒng)凸優(yōu)化方法和本文方法進(jìn)行平面陣列綜合的對比實驗。以下仿真實驗均使用配置為Intel?Core?i5?4570處理器、主頻為3.2 GHz、內(nèi)存為4 GB的計算機(jī)。

仿真參數(shù)設(shè)置如下：平面陣列x軸方向的陣列孔徑為4.5λ，y軸方向的陣列孔徑為8λ，x軸方向陣元間距dx=，y軸方向陣元間距dy=，則初始化平面陣列的橫向陣元數(shù)M=19，縱向陣元數(shù)N=33。要求陣列綜合后波束方向圖峰值旁瓣電平小于-20 dB，設(shè)置方向參數(shù) μ和 v在其取值范圍內(nèi)的采樣數(shù)分別為Lμ=41和Lv=41，誤差最小值 ξx=ξy=10-3。

圖2（a）和圖2（b）分別是利用文獻(xiàn)[10]中的迭代凸優(yōu)化方法和本文方法進(jìn)行平面陣列綜合獲得的陣元位置分布圖。圖3（a）和圖3（b）分別為文獻(xiàn)[10]方法和本文方法綜合后的平面陣列波束方向圖。

圖2　利用文獻(xiàn)[10]方法和本文方法綜合后的陣元位置分布圖

由圖2和圖3可知，兩種方法綜合后的平面陣列的波束旁瓣電平均控制在-20 dB以下；文獻(xiàn)[10]方法需要3 168 s的綜合時間獲得陣元數(shù)為41的平面稀疏陣列，這是由于二維觀測角度所對應(yīng)的旁瓣約束條件較多，影響了平面陣列綜合的收斂速度；而本文方法僅需3 s即可快速完成平面陣列的綜合，綜合后的平面稀疏陣列陣元數(shù)為48，由于本文方法綜合后的平面陣列具有大范圍低波束旁瓣電平以及更短的陣列孔徑，這使得其在平面陣列的綜合稀疏程度方面要稍劣于文獻(xiàn)[10]中的方法。平面陣列具有大范圍低波束旁瓣電平，可以降低外界干擾信號的影響，而更短的陣列孔徑可降低平面陣列對工作環(huán)境的限制，使得布陣更具有靈活性。

圖3　文獻(xiàn)[10]方法和本文方法綜合后的平面陣列波束方向圖

4　結(jié)語

由于平面陣列包含陣元數(shù)較多且其二維波束旁瓣的抑制范圍較大致使凸優(yōu)化約束條件數(shù)劇增，從而導(dǎo)致平面陣列綜合的耗時很長。本文提出了一種基于凸優(yōu)化的低復(fù)雜度平面陣列綜合方法，利用Kronecker積的性質(zhì)將平面陣列綜合中的二維波束旁瓣約束條件轉(zhuǎn)化為一維波束旁瓣的約束條件，大大減少了凸優(yōu)化的約束條件數(shù)，從而加快了平面陣列綜合的速度。本文提出的平面陣列綜合方法在收斂速度、陣列孔徑長度以及低波束旁瓣范圍等方面均優(yōu)于傳統(tǒng)的凸優(yōu)化方法。

參考文獻(xiàn)

[1]吳海洲，王鵬毅，郭肅麗.全空域相控陣測控系統(tǒng)波束形成分析[J].無線電工程，2011，41（11）：13?15．

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

[3]彭祥龍.用遺傳算法優(yōu)化任意稀布率的平面陣列[J].電訊技術(shù)，2007，47（3）：153?158.

[4]廖先華，楊建紅，張立軍，等.基于模擬退火算法的平面稀疏陣優(yōu)化[J].現(xiàn)代雷達(dá)，2012，34（10）：57?59.

[5]DELIGKARIS K V，ZAHARIS Z D.Thinned planar array de?sign using Boolean PSO with velocity mutation[J].IEEE Tran?sactions on Magnetics，2009，45（3）：1490?1493.

[6]RAJO?LGLESIAS E，QUEVEDO?TERUEL O.Linear array syn?thesis using an ant?colony?optimization?based algorithm[J]. IEEE Transactions on Antennas and Propagation，2007，49 （2）：70?79.

[7]CEN L，SER W，CEN W，et al.Linear sparse array synthesis via convex optimization[C]//Proceedings of 2010 IEEE Interna?tional SymposiumonCircuitsandSystems.Paris：IEEE，2010：4233?4236.

[8]PRISCO G，D′URSO M.Maximally sparse arrays via sequen?tial convex optimization[J].IEEE Antennas and Wireless Pro?pagation Letters，2012，11：192?195.

[9]NAI S E，SER W，YU Z L，et al.Beampattern synthesis for linear and planar arrays with antenna selection by convex opti?mization[J].IEEE Transactions on Antennas and Propagation，2010，58（12）：3923?3930.

[10]FUCHS B.Synthesis of sparse arrays with focused or shaped beampattern via sequential convex optimizations[J].IEEE Transactions on Antennas and Propagation，2012，60（7）：3499?3503.

中圖分類號：TN957?34

文獻(xiàn)標(biāo)識碼：A

文章編號：1004?373X（2016）03?0009?04

doi：10.16652/j.issn.1004?373x.2016.03.003

收稿日期：2015?10?16

基金項目：國家自然科學(xué)基金（61302188，61372066）；江蘇省自然科學(xué)基金（BK20131005）

作者簡介：陳金立（1982—），男，浙江寧波人，講師，博士。主要研究方向為陣列綜合和陣列信號處理。

Low complexity planar array synthesis method based on convex optimization

CHEN Jinli，ZHONG Jiyang，CAO Huasong
（College of Electronic and Information Engineering，Nanjing University of Information Science and Technology，Nanjing 210044，China）

Abstract：Since the convex optimization method used to synthesize the planar array has long time consuming，a low com?plexity planar array synthesis method based on convex optimization is proposed.The two?dimensional beam sidelobe constraint condition in the planar array synthesis model based on weighted ?1norm minimization is transformed into the one?dimensional beam sidelobe constraint condition for horizontal array and longitudinal array by the method，which can reduce the number of convex optimization constraint conditions，and effectively promote the convergence speed of the algorithm.The simulation results show that the proposed method performs slightly worse than the traditional convex optimization method in the aspect of sparsity degree of synthesis planar array，but it can fast obtain the planar array with shorter aperture and broader low?beam sidelobe level.

Keywords：planar array；array synthesis；convex optimization；low complexity