運(yùn)用遺傳算法綜合稀疏陣列

1、 運(yùn)用遺傳算法綜合稀疏陣列RandyL.HauptIEEE高級(jí)會(huì)員摘要：大的天線(xiàn)陣列很難被稀疏用來(lái)獲得較低的旁瓣。對(duì)于非周期陣列的綜合問(wèn)題，傳統(tǒng)的統(tǒng)計(jì)學(xué)方法遠(yuǎn)遠(yuǎn)達(dá)不到最優(yōu)配置的要求。傳統(tǒng)的優(yōu)化方法不適合用來(lái)優(yōu)化多參數(shù)或離散參數(shù)的問(wèn)題。本文將介紹如何利用遺傳算法來(lái)優(yōu)化一個(gè)稀疏陣列，并在一個(gè)周期陣列上利用遺傳算法來(lái)決定哪一個(gè)陣元被稀疏從而抑制陣列的最大相對(duì)旁瓣電平。本文將呈現(xiàn)200個(gè)陣元的線(xiàn)陣和200個(gè)陣元的面陣的仿真結(jié)果，要求稀疏陣列的旁瓣電平低于-20dB，對(duì)于線(xiàn)陣，同時(shí)在掃描角和帶寬上進(jìn)行了優(yōu)化。I簡(jiǎn)介周期性陣列通過(guò)有目的的放置相同權(quán)值的陣元來(lái)產(chǎn)生幅度錐削的低旁瓣，在給定旁瓣要求的前提下利用

2、簡(jiǎn)單的分析方法推出陣元的位置是不可能的1。事實(shí)上，大部分周期陣列的分析方法都是試圖將陣列區(qū)域內(nèi)的陣元密度與幅度錐削的低旁瓣的振幅密度相聯(lián)系起來(lái)，同時(shí)保持陣列孔徑不變2。陣元密度在陣列中心達(dá)到最大并逐漸向邊緣稀疏，通常情況下，旁瓣電平在主瓣附近減小，而在遠(yuǎn)離主瓣的位置上增大3（這個(gè)通常是可接受的）。非周期陣列的綜合方法是在給定條件下達(dá)到均方旁瓣電平或是最大相對(duì)旁瓣電平的要求。稀疏一個(gè)陣列意味著從均勻間隔陣列或者是周期陣列中抽去部分陣元從而在給定孔徑的范圍內(nèi)產(chǎn)生所期望的幅值密度。連接在饋電網(wǎng)絡(luò)上的陣元的狀態(tài)是“開(kāi)”，而連接在匹配負(fù)載或虛負(fù)載上的陣元的狀態(tài)就是“關(guān)”。用稀疏陣來(lái)產(chǎn)生低旁瓣比稀布陣要簡(jiǎn)

3、單的多，稀布陣的陣元位置不確定的，有無(wú)限多的取值可能。稀疏陣有2Q種組合，其中Q是陣元數(shù)目，如果陣列是對(duì)稱(chēng)結(jié)構(gòu)，那么陣元位置的組合數(shù)將顯著減少。稀疏也可以看作是振幅錐度的量化，其中每一個(gè)陣元的振幅用一個(gè)比特來(lái)表示。稀疏一個(gè)大的陣列從而產(chǎn)生低旁瓣涉及到檢查相當(dāng)多的陣元位置組合，目的是找到最好的稀疏方式。只有當(dāng)陣列較小時(shí)無(wú)遺漏的檢查所有的組合才是可行的。大部分的優(yōu)化方法（例如單純形法、Powell方法、共軛梯度法等）不適合于稀疏陣列優(yōu)化，他們只能優(yōu)化一些連續(xù)變量而且會(huì)陷入局部最小值6，此外，這些方法是專(zhuān)們用來(lái)處理連續(xù)參數(shù)問(wèn)題的，而稀疏陣列所處理的是離散參數(shù)問(wèn)題。動(dòng)態(tài)編程法可以?xún)?yōu)化大型參數(shù)組（有很多

4、陣元），但是極易陷入局部最小值。模擬退火算法和遺傳算法是非常適合于稀疏陣列的優(yōu)化方法，它們不受所要優(yōu)化的變量個(gè)數(shù)的限制，盡管處理速度比較慢，但是這些算法可以處理大型陣列。它們是全局性的優(yōu)化方法，擁有能夠探測(cè)到當(dāng)前最小值之外的解決方案的隨機(jī)因子，同時(shí)該算法是收斂的。由于模擬退火算法與遺傳算法的全局特性同時(shí)缺乏輔助信息，這使得它們較其他的非全局性方法有較慢的收斂速度。Ruf利用模擬退火算法優(yōu)化了低冗余線(xiàn)性陣列8，這些天線(xiàn)被應(yīng)用到射電天文學(xué)中評(píng)估場(chǎng)景的光芒。天文學(xué)家喜歡設(shè)計(jì)能夠采樣他們所觀(guān)測(cè)到的空間分布的天線(xiàn)陣列。Ruf的方法能優(yōu)化比以往任何時(shí)候都要大的線(xiàn)性陣列。本文的目的是尋找到一個(gè)稀疏陣，使其最

5、大相對(duì)旁瓣電平最?。╮sll）。本文將介紹一種遺傳算法，用來(lái)在數(shù)值上優(yōu)化一個(gè)線(xiàn)陣和面陣9。遺傳算法是模仿基因重組與進(jìn)化，將優(yōu)化參數(shù)編碼為稱(chēng)作基因的二進(jìn)制串，并執(zhí)行復(fù)制、交叉、自然選擇和變異等遺傳操作來(lái)得到最優(yōu)解。對(duì)于陣列，這些算法得到了較先前優(yōu)化嘗試或是統(tǒng)計(jì)嘗試更好的稀疏結(jié)構(gòu)。其他的優(yōu)化方法不適合于大型陣列，而且統(tǒng)計(jì)方法也找不到最優(yōu)解。本文利用遺傳算法對(duì)200個(gè)陣元的線(xiàn)陣和面陣分別進(jìn)行了優(yōu)化，并且得到了最優(yōu)解。線(xiàn)性陣列通過(guò)優(yōu)化，使其在一系列的掃描角和頻率上的最大相對(duì)旁瓣電平低于了-20dB。II遺傳算法遺傳算法的目的是尋找一組參數(shù)使得函數(shù)的輸出結(jié)果最小，遺傳算法不同于大部分的優(yōu)化方法，因?yàn)樗?/p>

6、如下的一些特點(diǎn)10。1）遺傳算法的處理對(duì)象是經(jīng)過(guò)編碼的參數(shù)，而不是參數(shù)本身。2）遺傳算法采用多點(diǎn)搜索來(lái)代替單點(diǎn)搜索。3）遺傳算法不需要輔助信息。4）遺傳算法采用概率的變遷規(guī)則，而不是確定性規(guī)則。C圖1遺傳算法的流程圖圖1是遺傳算法的流程圖，接下來(lái)的段落將詳細(xì)介紹遺傳算法的每一個(gè)步驟，為了便于查詢(xún)，我們用A-F對(duì)每一步都進(jìn)行了標(biāo)記。所有參數(shù)的數(shù)值都經(jīng)過(guò)了二進(jìn)制編碼（步驟A）,每一個(gè)被編碼的參數(shù)按順序放置，從而形成一個(gè)被稱(chēng)作基因的二進(jìn)制長(zhǎng)串，每一個(gè)基因?qū)?yīng)于其量化參數(shù)所計(jì)算的函數(shù)結(jié)果，因此遺傳算法需要檢查有限的，但卻是有很多的參數(shù)組合方式形成的大量基因。假設(shè)一個(gè)大的基因有N串，每一串有B位參數(shù)，那

7、么一共就有2NB種基因。如果參數(shù)是連續(xù)的，那么由于進(jìn)行二進(jìn)制編碼時(shí)會(huì)產(chǎn)生量化誤差，這將會(huì)限制遺傳算法的性能。從積極的方面來(lái)看，遺傳算法理論上更適合于離散參數(shù)的優(yōu)化問(wèn)題。稀疏陣就是一個(gè)離散參數(shù)問(wèn)題，參數(shù)的每一位代表著陣元的狀態(tài)，例如“on”代表1,“off”代表0。舉個(gè)例子來(lái)看，一個(gè)六元的陣列可以用101101來(lái)表示，其中陣元位置2和5被抽去了，如果一個(gè)線(xiàn)陣關(guān)于其中心對(duì)稱(chēng)，那么其上的2N個(gè)陣元就可以利用一個(gè)N位的基因來(lái)表示，那么上面的例子就可以用101這個(gè)基因來(lái)代替。這個(gè)基因的適應(yīng)度就是遠(yuǎn)場(chǎng)方向圖的最大相對(duì)旁瓣電平，本文的適應(yīng)度函數(shù)是點(diǎn)源陣列的相對(duì)遠(yuǎn)場(chǎng)方向圖，目的就是使相對(duì)旁瓣電平的最大值盡可能

8、的減小，其中陣元是保留還是抽去對(duì)結(jié)果有很大的影響。遺傳算法模仿自然界的基因重組和進(jìn)化10，在自然界里，基因是最基本的構(gòu)建塊。遺傳算法開(kāi)始于對(duì)輸出結(jié)果的隨機(jī)采樣，這種方式得到的基因大都有很差的最大相對(duì)旁瓣電平。結(jié)果只有產(chǎn)生好的輸出結(jié)果的基因被保留下來(lái)，而產(chǎn)生差的輸出結(jié)果的基因?qū)⒈粊G棄。舉個(gè)例子來(lái)看，表I列舉了8個(gè)隨機(jī)生成的基因（步驟B）。表I隨機(jī)生成的8個(gè)基因及其對(duì)應(yīng)的相對(duì)峰值旁瓣電平（RSLL）基因編號(hào)二進(jìn)制編碼串最大相對(duì)峰值旁瓣電平（dB）11111011001-9.0621111111001-14.0831111011110-12.1341101111001-9.7151111011101

9、-11.0961110110010-7.3171111110010-12.7681011110101-8.70這些基因代表著20個(gè)間隔為0.5久的各向同性的點(diǎn)源陣列，每一個(gè)基因可以通過(guò)如下的公式得到其所對(duì)應(yīng)的相對(duì)峰值旁瓣電平、NcosGnnnd+6）1）FF(u)-max2丈aelpat(u)；nFFn-1maxCu1OOOCllDanOOCi1lOlOllDlID11311Hi!miinmiiHiiniiuimt111111111x1moiti11I1N11111KWlOdi11UCMOCHW1C0D115.袖邱砒.期陽(yáng)科tfl禪isbnpicBlftfrWitd-D5i劭G51n龜cp更w

10、EdBllJ*aT-g&一atuiGu圖5(a)主波束為側(cè)射陣圖5(b)主波束掃描到1200用遺傳算法進(jìn)行優(yōu)化時(shí)遇到的另一個(gè)問(wèn)題就是零點(diǎn)干擾，譬如下面的例子,一個(gè)陣元間隔是d=0.5九，頻譜帶寬是10%，掃描范圍是300，已經(jīng)知道了在u=0.025與0.075之間存在干擾，現(xiàn)在的問(wèn)題是在已知角度干擾范圍與帶寬限制和陣列的掃描角的條件下來(lái)稀疏一個(gè)陣列從而獲得盡可能低的旁瓣電平。利用遺傳算法來(lái)優(yōu)化一個(gè)稀疏率為83%、帶寬為10%、角度范圍為u=0.025到0.075之間的陣列，所產(chǎn)生的最大相對(duì)峰值旁瓣電平是-26.dB。圖6(a)是當(dāng)=900且d=0.5九時(shí)稀疏陣的遠(yuǎn)場(chǎng)方向圖，陣元位置標(biāo)在了圖片的

11、上方，圖6(b)是同等優(yōu)化下主波束掃描到1200且d=0.55九時(shí)的遠(yuǎn)場(chǎng)方向圖，所有方向圖中的最大相對(duì)峰值旁瓣電平都低于-16dB。11111111111111111(11111111111111111111111111111111111曲(rfipe創(chuàng)冊(cè)BEpqiltommwc皿匕rHM締I仙ulmin111111111111itnnulinuhnmiinn11111111IlWllUlHllDl11001101：001001Q11000011Q01S3%Mind-ifintn?p*s聞?wù)緉-1擊-ff5).n-iI圖6(a)主波束在側(cè)射陣圖6(b)主波束掃描到1200接下來(lái)的例子是在一個(gè)矩

12、形方格內(nèi)稀疏一個(gè)由20 x10的各向同性點(diǎn)源所組成的平面陣，遠(yuǎn)場(chǎng)方向圖由下面的計(jì)算公式給出：FfQ申）=4比MacosEmmnn=1m=1ihdsinGcosecosEn-ihdsinGsinyx其中：（2）M=y方向的陣元數(shù)。N=x方向的陣元數(shù)。dy=y方向的陣元間隔。dx=x方向的陣元間隔。這個(gè)公式假設(shè)陣列位于x-y平面內(nèi)，而且陣列關(guān)于x軸和y軸都對(duì)稱(chēng)，因而需要對(duì)用于線(xiàn)陣的遺傳算法進(jìn)行一下改變。首先，基因可以被描繪在多個(gè)方面。每個(gè)稀疏陣列的結(jié)構(gòu)被加載到一個(gè)單一的一維數(shù)組中，而且遺傳算法執(zhí)行其在線(xiàn)陣中相同的優(yōu)化方式。另一種可行的方案是將稀疏結(jié)構(gòu)加載到一個(gè)二維的矩陣中，A矩陣提供新的有趣可行的

13、交配方案，由于基因是二維的，因此交叉點(diǎn)也就變成了交叉線(xiàn)，交叉線(xiàn)可以按照水平的、垂直的、對(duì)角線(xiàn)來(lái)畫(huà)出，用于配對(duì)的基因需要交換交叉線(xiàn)兩邊的基因信息。第二個(gè)要考慮的是適應(yīng)度該怎么計(jì)算。下面的例子中，適應(yīng)度是通過(guò)將=0。和=90。的面內(nèi)的最大相對(duì)峰值旁瓣電平相加來(lái)求得的。圖7描述了陣元的狀態(tài)，白色表示存在陣元，黑色表示不存在陣元，一共有108個(gè)陣元，因此陣列稀布率是54%,圖8是陣列的遠(yuǎn)場(chǎng)方向圖，圖9是=0。和=90。的疊繪遠(yuǎn)場(chǎng)方向圖。vr圖820 x10面陣的遠(yuǎn)場(chǎng)方向圖圖720 x10面陣的陣元分布圖白色代表存在陣元，黑色代表不存在陣元（各向同性點(diǎn)源）其他的適應(yīng)度計(jì)算準(zhǔn)則也是可行的。如果適應(yīng)度僅是=

14、00面內(nèi)的最大相對(duì)峰值旁瓣電平之和，那么優(yōu)化后的最佳值為-22.17dB。=00面內(nèi)的最大相對(duì)峰值旁瓣電平是-12.93dB，此時(shí)，陣列的稀布率是58%。適應(yīng)度準(zhǔn)則依賴(lài)于陣列的設(shè)計(jì)要求。nili-nQ,maxsII-2007旳002D代OG08EinftFieiaj圖9=Oo和=9Oo時(shí)疊繪的遠(yuǎn)場(chǎng)方向圖。=Oo面內(nèi)最大旁瓣電平是-20.07dB,Q=90o面內(nèi)最大旁瓣電平是-19.76dB。另一種用來(lái)稀疏面陣的方案是假設(shè)被稀疏陣的每一行是相同的，每一列也是相同的。公式（2）中的amn等價(jià)于aman，其中am是y方向的振幅錐削，an是x方向的振幅錐削。在一個(gè)矩形方格內(nèi)有一個(gè)稀布率是81%的40

15、x40的由各向同性點(diǎn)源所組成的被優(yōu)化過(guò)的面陣，面陣的最大相對(duì)峰值旁瓣電平是-17.2dB,第4、7、35行、列都填充有陣元，第38行、列沒(méi)有陣元?？偨Y(jié)本文介紹了如何利用遺傳算法來(lái)優(yōu)化稀疏線(xiàn)陣和面陣從而獲得盡可能低的相對(duì)峰值旁瓣電平，遺傳算法的優(yōu)點(diǎn)在于它可以?xún)?yōu)化大量的離散參數(shù)。以前利用統(tǒng)計(jì)學(xué)的方法來(lái)表示稀疏陣的幅度錐削的辦法沒(méi)有達(dá)到最佳的稀疏結(jié)果。利用遺傳算法的智能搜索從而找出能產(chǎn)生低旁瓣的最佳稀疏方式。本文論證了多種利用遺傳算法來(lái)優(yōu)化陣列方向圖的方法。許多額外的擴(kuò)展是可行的，包括圓形陣列，由有向陣元組成的面陣，掃描陣列等等。遺傳算法對(duì)于諸如自適應(yīng)零點(diǎn)的這類(lèi)實(shí)時(shí)控制模式不是很適用而且收斂很慢。另

16、一方面，遺傳算法很適用于優(yōu)化陣列設(shè)計(jì)。本文中遺傳算法被用來(lái)優(yōu)化稀疏陣列，但是實(shí)際上，遺傳算法還廣泛應(yīng)用于頻散和天線(xiàn)理論中。參考文獻(xiàn)Y.T.Lo,“Aperiodicarrays,”inAntennaHandbook,Theory,Applications,andDesign,Y.T.LoandS.W.Lee,Eds.NewYorkVanNostrand,1988.R.E.Willey,“Spacetaperingoflinearandplanararrays,”IREAntennasPropagat,vol.AP-10,pp.369-377,July1962.J.Galejs,“Minimiz

17、ationofsidelobesinspacetaperedlineararrays,”IEEETrans.AntennasPropagat.,pp.835-836,Sept.1964.B.D.Steinberg,PrinciplesofApertureandArraySystemDesign.NewYorkWiley,1976.Y.T.LoandS.W.Lee,“Astudyofspacetaperedarrays,ZEEETrans.AntennasPropagat.,vol.AP-14,pp.22-30,Jan.1966.W.H.Pressetal.,NumericalRecipes.NewYork:CambridgeUniversityPress,1992.L.Davis,Ed.,GenericAlgorithmsandSimulatedAnnealing.LosAltos,CA:MorganKaufmann,1987.C.S.Ruf,“Numericalannealingoflow-redundancylineararrays,”ZEEETra