一、引　　言
　　天線的最大相對旁瓣電平是評價天線性能的一個重要參數，在給定天線形狀與陣元個數的前提下，如何通過恰當地選擇各陣元的間距、饋電流幅值及相位來最大限度地降低旁瓣電平是陣列天線綜合中的一類重要課題.對于形狀復雜的大陣列天線，傳統的解析法(如道爾夫——切比雪夫綜合法等)難以計算，采用數值分析方法較為適宜.由于天線最優化問題中的目標函數或約束條件大多呈多參數、非線性、不可微甚至不連續，因此基于梯度尋優技術的傳統數值優化方法無法有效地求得工程滿意解；近年來一種模擬自然進化的遺傳算法開始應用于計算電磁學領域［1，2］，該算法只要求待解問題是可計算的，并無可微性等其它限制，同時，由于該算法采用了優化的隨機搜索技術，能以較大的概率和較快的速率求得全局最優解.本文運用遺傳算法對一個具有1024個陣元的矩形平面陣列的陣元間距及饋電流幅值進行了優化，使該方陣的最大相對旁瓣電平由均勻方陣的-13.27dB降至-34.56dB.結果表明，遺傳算法對解決天線系統中大量復雜的最優化問題具有廣闊的應用前景.

二、陣列天線旁瓣電平的優化問題
　　1.陣列天線的方向圖函數
　　考慮由2Nx行2Ny列陣元構成的矩形平面陣列.各陣元的相位相同，間距dxi、dyj及歸一化電流振幅Ixm、Iyn可以不同，但關于x軸和y軸對稱，如圖1.設陣元的方向圖函數為cosθ，則此面陣的方向圖函數為［4］：

　　(1)

這里，k=2π/λ，λ為波長.第m行n列陣元的電流振幅由下式計算：

I_mn=I_xmI_ynI₀(2)

式中，I0為電流振幅基數.
 

圖1　矩形平面陣列結構示意圖

　　2.用遺傳算法優化陣列天線旁瓣電平
　　(1)編碼方案與計算流程
　　應用遺傳算法時首先應對解參數進行編碼，陣元的歸一化電流及間距的編碼方案由下面兩式給出：

　(3)
　(4)

式中，t為調節參數，在本文的實例中取為t=12.5;b為二進制碼向量，其元素取值為0或1；運用遺傳算法的計算流程圖如圖2所示.

圖2　遺傳算法計算流程圖

　　(2)自適應變異算子
　　為避免迭代后期出現封閉競爭而導致的進化滯緩及未成熟收斂現象，本文提出了根據每代最佳個體進化的狀況動態調整變異比例的自適應變異算子.設Pm為變異概率，其在每代中的取值由下式確定：

P_m=P_m0(R^S+1)　　(5)

式中Pm0為初定的變異比率；R為最佳個體連續未進化的代數；S為變異參數，本文的計算實例中取為S=2.
　　對比研究表明，采用自適應變異算子能較好的提高進化速率及優化結果.

三、計算實例
　　選擇一個32行×32列的矩形平面陣列作為研究實例，對其最大相對旁瓣電平進行了優化.選取dxi=dyi、Ixm=Iym,根據方向圖相乘原理和本題的對稱性，只須計算一行中位于y軸一側的16個陣元的優化分布.該陣列的嘗試解的規模達到了(32×32)16=1.4615×1048，但遺傳算法運用其特有的選擇優化策略，僅通過對其中的28200個嘗試解進行評估后，便給出了最大相對旁瓣電平的優化值：-34.56dB，此結果較之均勻面陣的-13.27dB［4］降低了21.29dB，優化的效果是明顯的.圖3給出了經優化后的非均勻面陣在X-Z主面內的方向圖.與之對應的陣元間距及歸一化電流見表1.

圖3　已優化的非均勻面陣X-Z主面方向圖

表1　經優化后的陣元電流振幅及間距
 

　　圖4描繪了本例中經優化后的平面陣列在三維空間中的立體方向圖.從圖中可以清晰地看到，除了在X軸(Φ=0°)和Y軸(Φ=90°)處各有一列副瓣低于-34.56dB的波峰之外，在其它區域內沒有電平高于-35dB的副瓣波峰；而且，在副瓣被大幅度降低的情況下，主瓣并未明顯展寬.計算得知，具有相同陣元數的均勻面陣的主瓣寬度為3.23°，而經優化后的面陣的主瓣寬度為4.15°，僅比優化前增加了0.92°，這表明該面陣仍具有良好的方向性.

圖4　優化后的平面陣列立體方向圖

四、結　　論
　　用遺傳算法求解大型陣列天線的優化問題，具有算法結構簡單、搜索效率高、優化結果好的特點.文中提出的自適應變異算子能夠避免封閉競爭導致的進化滯緩及未成熟收斂現象，對豐富遺傳算法的應用作了有益的探討.