基于自適應混合差分進化算法的大規模稀布子陣優化設計

倪濤

【摘 要】本文對大規模寬帶稀布線陣采用子陣劃分的方法進行優化設計。將大規模陣列劃分為若干個子陣降低了優化的計算量;通過每個子陣和子陣間的不等間距設計有效地降低了柵瓣。以寬帶類八木印刷天線為例，采用自適應混合差分進化算法同時優化陣列的單元間距和激勵系數，分別對筆形波束、掃描波束以及單脈沖波束做了優化設計，結果表明該方法是可以有效解決大規模稀布子陣優化問題。

【關鍵詞】稀布陣列;線陣;子陣;自適應混合差分進化算法

中圖分類號： TN820 文獻標識碼： A 文章編號： 2095-2457（2019）11-0020-004

【Abstract】In this paper，? large-scale broadband sparse linear array is optimized by the method of sub-array partitioning. Dividing a large-scale array into several sub-arrays reduces the amount of computation for optimization; effectively reducing the grating lobes by unequal spacing between each sub-array and sub-arrays. Taking the broadband printed Quasi-Yagi antenna? as an example， the self-adaptive hybrid differential evolution algorithm is used to optimize both? the element spacing and excitation of the array. The pencil beam， scanning beam and monopulse beam are optimized respectively. The results show that the method can effectively solve the problem of large-scale sparse sub-array optimization.

【Key words】Sparse array; Linear array; Sub-array; Self-Adaptive Hybrid Differential Evolution

0 引言

不等間距天線陣列一直以來收到了學者們的廣泛關注。相比于周期的陣列，不等間距陣列打破了陣因子的周期性，陣列隨機性的優勢被逐步挖掘。首先，對于寬帶天線陣列，采用不等間距陣列可以有效地抑制陣列在高頻掃描工作時的柵瓣。第二，對于等幅同相的周期陣列，其副瓣電平理論值為-13.4dB，采用不等間距陣列可以有效地降低副瓣電平。第三，對于相同口徑的陣列，采用不等間距陣列可以降低單元數，從而提升陣列口徑的效率[1-2]。

近年來，針對不等間距陣列的研究主要關注在以下兩類：（1）稀疏陣：單元為等間距排列，但其激勵系數為隨機分布，部分單元雖然在陣列中間但是不工作;（2）稀布陣：單元間距為隨機分布。稀布陣由于其單元間距可以作為設計變量，因此大大提升了陣列設計的自由度，因此受到了廣泛的關注。對于大規模寬帶陣列，陣元間距優化的計算量比較大，如何快速優化一個陣列獲得期望的輻射方向圖是非常值得研究的問題。對于大規模陣列，采用子陣劃分技術是一種非常有效的手段，不僅可以有效地降低設計問題的維度，也可以降低工程造價的難度。對于該類陣列的優化是一個非線性不可導的數學問題，近年來，進化類算法如粒子群算法、差分進化算法、蟻群算法[3-4]等已經廣泛應用于該類問題。在文獻[5]中，作者提出了一種基于新型的自適應策略，將簡化的三點二次差值算子應用于差分進化算法的參數控制中，提出了一種自適應混合差分進化（Self-Adaptive Hybrid Differential Evolution，? 簡稱SAHDE）算法，該算法已經廣泛應用于很多天線陣列設計問題中，其收斂速度更快、結果更優。

在本文中，采用SAHDE算法對大規模寬帶稀布線陣進行優化設計。為了降低設計的變量數目和減少計算量，將大規模陣列劃分為若干個子陣;為了抑制陣列的柵瓣，每個子陣的間距是不等間距，每個子陣間的間距也為不等間距。采用寬帶類八木印刷天線作為單元天線，其工作帶寬為8GHz-12GHz。為了獲得較低的副瓣電平，單元間距和激勵系數同時作為優化變量。筆形波束、掃描波束以及單脈沖波束在本文中均分別做了優化設計，結果表明該方法是一種有效解決大規模稀布子陣優化問題的方法。

1 稀布子陣原理分析

圖1給出了一維線陣的結構示意圖。對于給定頻率fre，N元稀布陣列的陣列方向圖由下式給出

式中，Ii是第i個單元的激勵電流;Di-1是第i個單元和第一個單元之間的距離且D0=0;fre表示頻率;ffre（θ）是在給定頻率fre下的單元方向圖;kfre表示頻率fre的傳播常數;N是陣列單元的總數。

優化目標設計為在寬頻帶內陣列的副瓣電平，適應度函數定義為：

根據表達式（2），優化變量的數量為2*N-1。基于子陣的思想，N元稀布陣列可以分為幾個稀布陣列，原理圖如圖1（b）所示。因此，表達式（1）可以寫為：

優化目標設計仍為在寬頻帶內陣列的副瓣電平，因此，該問題的適應度函數定義為：

其中M是每個子陣的單元數量，MM是劃分子陣的數量。Ii是每個子陣中第i個單元的激勵電流;Isubi是第i個子陣的激勵電流。Dsubi-1是第i個子陣和第一個子陣之間的距離且Dsub0=0;根據表達式（4），優化變量的數量是2*（M+MM）-2。對于天線陣列，單元的數量N=M*MM。通過表達式（2）和（4）分別優化適應度;優化變量的數量分別為2*M*MM-1和2*（M+MM）-4;對于大規模陣列，2*（M+MM）-2遠小于2*M*MM-1，更容易獲得收斂結果。

因此，對于陣列副瓣電平的優化設計便轉換為對表達式（4）最小值的優化。

2 單元設計

單元采用印刷的類八木天線，詳細參數如圖2所示。貼片印刷在相對介電常數為10.2，厚度為0.635mm的介質基板上。單元的S11仿真結果如圖3（a）所示，在7.7GHz-12GHz時S11小于-10dB。單元帶寬超過40%。XOZ平面內不同頻率點的電場強度方向圖如圖3（b）所示。在頻帶8-12GHz內選取了9個采樣點，采樣間隔為0.5GHz。

3 陣列優化設計

3.1 筆形波束

在本節中對筆形波束進行了優化。單元的數量N=1000，并且該陣列被分成MM個子陣，每個子陣的數量是M，因此N=M*MM。單元間距和激勵是優化變量。該方向圖由式（3）進行計算計算。適應度函數和約束條件如下：

其中HPBW（frei）desired是第i個頻點期望獲得的半功率波束寬度;HPBW（frei）G是優化過程陣列獲得的波束寬度，λ是所有頻點的最大波長。為了保持天線陣列的孔徑恒定，將DM-1和DsubMM-1設為常數。激勵幅度限制在[0.2，1]。

對于波束寬度的約束采用罰函數法，H[.]是Heaviside階躍函數。罰因子：KK=1000。所需的HPBW設定為0.08°。SAHDE算法的控制參數設置為：種群規模為NP=100.優化代數Gmax=100，其余參數和文獻【5】中保持一致。在表1中，給出了稀布陣列的不同劃分方法，組1，M=1000，MM=1表示相等的間距和激勵，組2-9表示八種不同的分割方式，其中間距和激勵同時被優化。第1組的PSLL在9個采樣頻率點處的值分別為-13.31，-13.38，-13.27，-13.42，-5.73，-2.67，-2.69，-3.42，-3.36dB。在較低頻率點，陣列的PSLL接近-13.4dB，并且柵瓣將出現在較高頻率點。根據圖5，組2-9的副瓣電平優于組1，同時HPBW在所有頻率點均小于或等于0.08°。對于不同的劃分方法，優化變量的數量并沒有不同，并且子陣的空間色散在工作頻帶中也是不同的。，因此分工方式取決于實際工程應用。

3.2 掃描波束

在相控陣中，主波束在空間內自由掃描， 在本節中，對主波束掃到特定角度的陣列，優化其陣列參數。該陣列方向圖表示如下：

其中是掃描角度，適應度和約束條件由表達式（5）給出。 所需的半功率波瓣寬度設定為0.1°。 控制參數設置和3.1節中相同。如圖6所示，第一條曲線是主波束掃描角度為0°，副瓣電平小于-24.04 dB。 2-4條曲線分別代表主波束掃到10°，30°和45°，此時副瓣電平分別小于-23.08dB，-21.16dB和-18.77dB，同時波束寬度小于0.1°。

3.3 單脈沖波束

單脈沖雷達系統因其高精度角度測量而被廣泛使用。單脈沖天線要求通過相同的陣列配置生成和差方向圖。在這部分中，稀布線陣的結構示意圖如圖7所示。

陣列在全局坐標坐標中是對稱的，并且每個子陣在局部坐標中也是對稱的。陣列方向圖表示如下：

由于陣列的對稱性，對于和波束，Isubii=Isub-ii，Iii=I-ii時，那么

對于差波束，Isubii=Isub-ii，Iii=I-ii時，

約束條件如下所示：

目標函數如下所示：

其中MM=10，M=25，單元數為1000;w1=0.4，w2=0.6，kk=1000;HPBW 和HPBW 分別設置為0.08和0.07。當和波束獨立優化時，HPBW小于0.08°，PSLL小于-22.64dB。當差波束獨立優化時，HPBW在所有頻率點都小于0.07°，并且PSLL小于-20.39dB。PSLL的結果如圖8所示。對于實際工程問題，一組激勵對于和波束和差波束需要同時獲得較低的副瓣電平。當同時優化和差方向圖時，結果如圖8所示，和方向圖的副瓣電平小于-19.36dB，半功率波瓣寬度也小于0.08°，差方向圖以和方向圖的最大值作歸一化，差波束的半功率波瓣寬度小于0.07°，副瓣電平于-18.75dB。

4 總結

在本文中，采用SAHDE算法對大規模寬帶稀布線陣進行優化設計。通過稀布的方式有效地抑制了陣列的柵瓣，通過采用子陣劃分的方式有效地降低了優化設計的變量數目和計算量。采用工作帶寬為8GHz-12GHz寬帶類八木印刷天線作為單元天線，以單元間距和激勵系數同時作為優化變量，以陣列的副瓣電平為優化目標;分別對筆形波束、掃描波束以及單脈沖波束三種波束做了優化設計，獲得陣列的的優化結果。

【參考文獻】

[1]B.P.Kumar and itle and G.R.Branner，“Design of unequally for performance improvement，”IEEE Trans.Antennas Propag.， Vol.47，No.3，pp.511-523，Mar.1999.

[2]R.L.Haupt，“Unit circle representation of aperiodic arrays，” IEEE Trans.Antennas Propag.，Vol.43，No.10，pp.1152-1155，Oct.1995.

[3]R.L.Haupt，“Thinned arrays using genetic algorithms，”IEEE Trans.Antennas Propag.，Vol.42，No.7，pp.993-999，Jul.1994.

[4]K.S.Cheng，Z.S.He，and C.L.Han，“A modified real GA for the sparse linear array synthesis with multiple constraints，”IEEE Trans.Antennas Propag.，Vol.54，No.7，pp.2169-2173，Jul.2006.

[5]T.Ni，Y.-C.Jiao，et? al.Worst-case? tolerance? synthesis? for low-sidelobe? sparse? linear? arrays? using? a? novel self-adaptive? hybrid? differential? evolution? algorithm[J].Progress? in? Electromagnetics? Research? B.2016，66：91-105.作者簡介：倪濤（1988—），陜西安康人，工程師，主要從事天線陣列研究。