陳軒++黃心漢
收稿日期:2013-06-18
作者簡介:陳 軒(1990—),男,湖北漢川人,碩士研究生,研究方向:智能控制與智能應(yīng)用。
文章編號:1003-6199(2014)02-0089-04
摘 要:區(qū)域防空雷達網(wǎng)是防空作戰(zhàn)空情預(yù)警的發(fā)展趨勢,為了提高區(qū)域雷達網(wǎng)探測能力和抗綜合電子干擾、抗隱身技術(shù)與隱身飛機的威脅,抗低空、超低空突防及抗反輻射導(dǎo)彈(ARM)的能力,研究雷達組網(wǎng)的問題,介紹免疫遺傳算法的基本原理和特點,提出基于免疫遺傳算法的雷達組網(wǎng)方法,通過計算機仿真實驗證明方法的可行性。
關(guān)鍵詞:免疫遺傳算法;雷達組網(wǎng);覆蓋系數(shù)
中圖分類號:TP39文獻標(biāo)識碼:A
Approach to Radar Netting Based on Immune Genetic Algorithm
CHEN Xuan, HUANG Xinhan
(School of Automation, Huazhong University of Science and Technology, Wuhan,Hubei 430074,China)
Abstract:Regional airdefense radar netting is the tendency of air defense warfare in future, in order to improve the detective performance and the performance of anti synthesized electronic interference, anti stealth technique and stealth aircraft, anti low, ultra low altitude penetration and anti anti-radiation missiles(ARM), radar netting is studied. The principle and characteristics of immune genetic algorithm is introduced, and immune genetic algorithm is put forward in radar netting. The simulation experiment by computers indicates valid of this method.
Key words:immune genetic algorithm; radar network; coverage coefficient
1 引 言
隨著新型空襲兵器和航天技術(shù)的不斷發(fā)展,單臺雷達無論在探測能力上,還是電子防御功能上都有較大的局限性,因此雷達組網(wǎng)技術(shù)在現(xiàn)代戰(zhàn)爭和人類對宇宙的探測中起著舉足輕重的作用,該技術(shù)是通過將多功能、多類型、多頻段的雷達進行組網(wǎng),實現(xiàn)實時數(shù)據(jù)傳輸,資源共享,并在此基礎(chǔ)上對數(shù)據(jù)實時處理,這已被證明是一種有效的方法[1]。雷達組網(wǎng)系統(tǒng)的關(guān)鍵問題是如何對多臺雷達進行最佳組網(wǎng)以獲得最優(yōu)的防御能力。目前進行雷達組網(wǎng)的方法很多,常用的有效用函數(shù)法或?qū)<曳?,它們是根?jù)雷達覆蓋防守區(qū)域面積、雷達部署任務(wù)、單臺雷達探測距離、地形、銜接角、遮蔽角、起伏角等因素進行加權(quán)求和得到陣地的效用值,但是這些都不能很好地解決執(zhí)行速度的問題。將免疫遺傳算法應(yīng)用于雷達組網(wǎng),能較快地使布陣接近最優(yōu)解,從而避免了采用窮舉的方法帶來執(zhí)行速度慢的問題,并且克服了遺傳算法未成熟收斂和局部搜索能力差的缺陷。
2 免疫遺傳算法原理和設(shè)計
2.1 免疫遺傳算法原理
生物免疫系統(tǒng)對抗原會自動產(chǎn)生相應(yīng)的抗體來防御,這一過程被稱為免疫應(yīng)答。在此過程中,部分抗體作為記憶細胞保存下來,當(dāng)同類抗原再次侵入時,記憶細胞被激活并產(chǎn)生大量抗體,使再次應(yīng)答比初次應(yīng)答更迅速,體現(xiàn)了免疫系統(tǒng)的記憶功能。同時,抗體與抗體之間也相互促進和抑制,以維持抗體的多樣性及免疫平衡,這種平衡是依濃度機制完成的,即抗體的濃度越高,則越受抑制,反之亦然,體現(xiàn)了免疫系統(tǒng)的自我調(diào)節(jié)功能??贵w的濃度計算是系統(tǒng)保持種群多樣性的基本手段之一[2]。
傳統(tǒng)的遺傳算法是一種具有“生成+檢測”的迭代過程的搜索算法,而免疫遺傳算法IGA (Immune Genetic Algorithm,IGA)則是一種借鑒生物免疫系統(tǒng)的自適應(yīng)識別和排除侵入機體的抗原性異物的功能,將生物免疫系統(tǒng)的學(xué)習(xí)、記憶、多樣性和識別特點引入的遺傳算法。在解決實際問題時,目標(biāo)函數(shù)和約束條件作為抗原輸入,隨后產(chǎn)生初始抗體群,計算抗原和抗體的親和力用來描述可行解與最優(yōu)解的逼近程度,并通過一系列遺傳操作的計算,在保持抗體多樣性的情況下,找出針對該抗原的抗體,即問題的解[3]。免疫遺傳算法與傳統(tǒng)的遺傳算法相比,具有如下顯著特點:
1)具有免疫記憶功能:可加快搜索速度,提高總體搜索能力,確保快速收斂于全局最優(yōu)解;
2)具有抗體的多樣性保持功能:可提高全局搜索能力,避免未成熟收斂;
3)具有自我調(diào)節(jié)功能:可提高局部搜索能力。
2.2 免疫遺傳算法設(shè)計
免疫遺傳算法在標(biāo)準(zhǔn)遺傳算法的基礎(chǔ)上增加了抗體濃度概率計算、抗體的促進與抑制等模塊來提高解的多樣性。該算法因為將免疫系統(tǒng)中抗體多樣性維持機制引入了遺傳算法,使得其性能比標(biāo)準(zhǔn)遺傳算法更進了一步。在解決實際問題時,目標(biāo)函數(shù)和約束條件作為抗原輸入,隨后產(chǎn)生初始抗體群,并通過一系列遺傳操作及抗體親和度的計算,在保持抗體多樣性的情況下,找出針對該抗原的抗體,也就是問題的解[3]。免疫遺傳算法的流程圖如圖1所示,其基本的步驟如下:
計算技術(shù)與自動化2014年6月
第33第2期陳 軒等:基于免疫遺傳算法的雷達組網(wǎng)方法
1)算法初始化。輸入抗原及設(shè)定參數(shù):輸入目標(biāo)函數(shù)以及約束條件,作為抗原的輸入;設(shè)定種群規(guī)模、選擇概率、交叉概率、變異概率等參數(shù)。
2)初始抗體產(chǎn)生。在第一次迭代時,抗體通常在解空間中用隨機的方法產(chǎn)生。
3)親和度及濃度的計算。計算各抗體和抗原的親和度以及各抗體的濃度。
4)終止條件判斷。判斷是否滿足終止條件,如果滿足終止條件,將與抗原親和度最高的抗體加入細胞記憶數(shù)據(jù)庫,然后終止;否則繼續(xù)。
5)選擇、交叉、變異操作。根據(jù)設(shè)置的選擇概率、交叉概率和變異概率選擇抗體,對抗體進行交叉和變異操作。
6)根據(jù)以上的操作更新群體以后轉(zhuǎn)入步驟3)。
圖1 免疫遺傳算法流程圖
3 采用免疫遺傳算法的雷達組網(wǎng)
3.1 雷達組網(wǎng)問題描述
在采用免疫遺傳算法進行雷達組網(wǎng)的方案中,雷達的覆蓋系數(shù)為:
ρ=[(∪Ni=1Si)∩Sarea]/Sarea(1)
式中Si為第i臺雷達的探測范圍,N為雷達總數(shù),Sarea給定的責(zé)任區(qū)域。ρ∈[0,1] 表示N臺雷達所覆蓋的有效責(zé)任區(qū)域占總責(zé)任區(qū)域的比重。maxρ為本次雷達組網(wǎng)的目標(biāo)函數(shù),也就是免疫遺傳算法中的抗原,即
maxρ=max [(∪Ni=1Si)∩Sarea]/Sarea (2)
3.2 抗體編碼
雷達坐標(biāo)是所求問題解的信息,為了縮小抗體空間,提高搜索效率,采用了對雷達組網(wǎng)比較直觀的抗體編碼方式,假設(shè)一共有5臺雷達進行組網(wǎng),則N種雷達布陣的方案可以表示為如下的結(jié)構(gòu):
X11Y11X12Y12X13Y13X14Y14X15Y15第1種雷達組網(wǎng)方案
X21Y21X22Y22X23Y23X24Y24X25Y25
第2種雷達組網(wǎng)方案
……
Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5
第n種雷達組網(wǎng)方案
其中,Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5(n表示雷達組網(wǎng)的任意一種方案)分別是第1臺、第2臺、第3臺、第4臺和第5臺雷達的橫坐標(biāo)和縱坐標(biāo)。雷達布陣的每一種方案對應(yīng)為免疫遺傳算法中抗體種群中的每一個個體。
3.3 親和力計算
親和力是指發(fā)生免疫系統(tǒng)識別的抗體和抗原的結(jié)構(gòu)越互補,結(jié)合越可能發(fā)生,而把彼此結(jié)合的強度稱之為親和力。抗原需要盡可能好的與入侵的抗體相結(jié)合。二者匹配程度越好,結(jié)合就越好,抗原和抗體親和力就越大。對于雷達組網(wǎng)問題,抗原對應(yīng)的是雷達組網(wǎng)的最大覆蓋系數(shù),由于雷達組網(wǎng)的總責(zé)任區(qū)域的面積是一定的,可以定義抗體Ab與抗原Ag的之間的親和力為:
(Ab,Ag)=1Tb-Tg(3)
式中:Tb,Ta分別為抗體Ab,Ag對應(yīng)的雷達網(wǎng)的覆蓋系數(shù),其中Tg也是所求的最大的雷達網(wǎng)覆蓋系數(shù)。
定義抗體Ab1與抗體Ab2之間的親和力為:
App(Ab1,Ab2)=1|Tb1-Tb2|(4)
式中:Tb1、Tb2分別為抗體Ab1與抗體Ab2對應(yīng)的雷達網(wǎng)覆蓋系數(shù)[4]。
3.4 算法選擇、交叉、變異算子
免疫遺傳算法能使抗體保持多樣性并且最終能夠收斂到最優(yōu)解的主要操作,就是因為在算法中有選擇、交叉和變異算子的存在,從而使整個抗體群沿著適應(yīng)度較好的方向搜索。
1)選擇算子
選擇是根據(jù)適者生存原則選擇下一代抗體,在基于免疫遺傳算法的雷達組網(wǎng)中選擇的是下一代的組網(wǎng)方案。采用如下的選擇算子:
PS(xi)=αρ(xi)∑ni=1ρ(xi)+(1-α)1Ne-ciβ(5)
式中:ρ(xi)是以式(1)為適應(yīng)度函數(shù)的雷達組網(wǎng)覆蓋系數(shù);ci是抗體xi的濃度,即群體中相似抗體所占的比重[5],即
ci=和抗體i親和度大于λ的抗體數(shù)N(6)
其中λ為相似常數(shù),一般取為0.9~1;
α和β是常數(shù)調(diào)節(jié)因子; N是該種群內(nèi)抗體的總數(shù)。
2)交叉
交叉是在選中的抗體中,對兩個不同的個體按交叉概率Pc隨機選中相同的位置進行基因交換(一般交叉概率取值為0.15~0.75),從而產(chǎn)生新的抗體,也就是新的雷達組網(wǎng)方案。如果對于5臺組網(wǎng)雷達,隨機選擇的是抗體1和抗體2進行交叉,抗體1和抗體2的編碼如下:
抗體1:X11Y11X12Y12X13Y13X14Y14X15Y15
抗體2:X21Y21X22Y22X23Y23X24Y24X25Y25
按照交叉概率產(chǎn)生的交叉點為4,則交叉以后的抗體1和抗體2的編碼分別是:
抗體1:X11Y11X12Y12X23Y23X24Y24X25Y25
抗體2:X21Y21X22Y22X13Y13X14Y14X15Y15
3)變異
變異是在選中的抗體中,對個體中的某些基因以一定的變異概率對某些抗體的某些位執(zhí)行異向轉(zhuǎn)化(一般變異概率的取值為0.01~0.2)。在變異的時候,對于交叉過程中產(chǎn)生的抗體方案,隨機產(chǎn)生一個介于[-MAX/2, +MAX/2]的數(shù)字rand,其中MAX為橫坐標(biāo)和縱坐標(biāo)可取的最大值。變異以后坐標(biāo)值x與原來值x之間的關(guān)系如下:
x,=x+rand(7)
變異以后如果x′>MAX,則取x′=MAX;如果x,<0,則取x′= 0。
單靠變異不能在求解中得到好處,但是,它能保證算法過程不會產(chǎn)生無法進化的單一群體。因為在所有的個體一樣的時候,選擇和交叉無法產(chǎn)生新的基因,這時只能靠變異產(chǎn)生新的基因,即變異增加了全局優(yōu)化的特質(zhì)。
3.5 基于濃度的種群更新
由式(5)可見,本文中的選擇算子是采用基于抗體的雷達組網(wǎng)覆蓋系數(shù)以及抗體濃度的選擇概率Ps(xi),從而有效地保證了抗體的多樣性,避免了算法的早熟問題,能夠獲得更好的覆蓋系數(shù)。
4 仿真實驗
根據(jù)前面的算法流程,選取目標(biāo)函數(shù)為N臺雷達所覆蓋的有效責(zé)任區(qū)域占總責(zé)任區(qū)域的比重。雷達的臺數(shù)分別為3臺,4臺,5臺,初始抗體的數(shù)量為N=50,即一共有50種方案。利用matlab仿真的結(jié)果如圖2所示。圖中矩形框表示的是組網(wǎng)雷達陣地的有效范圍(長為200公里,寬為100公里),雷達的類型一共有兩種,即探測范圍為50公里的雷達和探測范圍為30公里的雷達。參照雷達組網(wǎng)的實際要求,在仿真的過程中,最多只提供2臺性能優(yōu)越的雷達,效能較差的雷達數(shù)量不限。
(a)3臺雷達的面積覆蓋率
(b)4臺雷達的面積覆蓋率
(c)5臺雷達的面積覆蓋率
圖2 matlab仿真結(jié)果
在圖2中當(dāng)雷達數(shù)量為3臺時,組網(wǎng)雷達的面積覆蓋率為87.84%(圖2(a));當(dāng)為4臺時,組網(wǎng)雷達的面積覆蓋率為91.53%(圖2(b));當(dāng)為5臺時,組網(wǎng)雷達的覆蓋率為94.27%(圖2(c))。將仿真的結(jié)果與傳統(tǒng)遺傳算法的組網(wǎng)覆蓋率對比[6],如表1所示。
表1 雷達組網(wǎng)覆蓋率對比
3臺雷達
4臺雷達
5臺雷達
傳統(tǒng)遺傳算法覆蓋率
81.5%
86.1%
88.2%
免疫遺傳算法覆蓋率
87.84%
91.53%
94.27%
由表1可以看出,基于免疫遺傳算法雷達組網(wǎng)的覆蓋率要高很多,將免疫遺傳算法應(yīng)用于雷達組網(wǎng)是一種有效的方法。
5 結(jié) 論
雷達組網(wǎng)能夠有效地提高雷達系統(tǒng)的整體性能,更好的適應(yīng)高科技電子戰(zhàn)、信息戰(zhàn)。將免疫遺傳算法應(yīng)用于雷達組網(wǎng)系統(tǒng),能夠有效地提高雷達網(wǎng)的面積覆蓋率。免疫遺傳算法與傳統(tǒng)的遺傳算法相比,能夠克服傳統(tǒng)遺傳算法的未成熟收斂,提高全局搜索能力。
參考文獻
[1] 彭獲由. 區(qū)域性雷達組網(wǎng)[J]. 電子科學(xué)技術(shù)評論, 1992(3) :1- 6
[2] LUO WENJIAN, CAO XIANBIN,WANG XUFA. An Immune Genetic Algorithm Based on Immune Regulation[A]. Proceedings of the 2002 Congress on Evolutionary ComputationCEC2002. May 12-17, 2002. Honolulu, USA. Vol.2: 801-806.
[3] 段玉波, 任偉建, 霍鳳財, 等. 一種新的免疫遺傳算法及其應(yīng)用[J], 2005, 20(10):1185-1188.
[4] 謝剛, 武斌, 謝克明. 基于免疫遺傳算法的TSP優(yōu)化問題求解[J]. 太原理工大學(xué)學(xué)報, 2007, 38(3): 199-201.
[5] 劉旭旺, 齊微. 免疫遺傳算法的研究及其在TSP中的應(yīng)用[EB/OL]. 中國科技論文在線,http://www.paper.edu.cn/releasepaper/content/200706-566, 2007-06-28
[6] 莫建文, 韓傳久, 張彤. 遺傳算法在雷達組網(wǎng)最優(yōu)化布陣中算法研究[J]. 信息技術(shù), 2005(3): 66-68.
X11Y11X12Y12X13Y13X14Y14X15Y15第1種雷達組網(wǎng)方案
X21Y21X22Y22X23Y23X24Y24X25Y25
第2種雷達組網(wǎng)方案
……
Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5
第n種雷達組網(wǎng)方案
其中,Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5(n表示雷達組網(wǎng)的任意一種方案)分別是第1臺、第2臺、第3臺、第4臺和第5臺雷達的橫坐標(biāo)和縱坐標(biāo)。雷達布陣的每一種方案對應(yīng)為免疫遺傳算法中抗體種群中的每一個個體。
3.3 親和力計算
親和力是指發(fā)生免疫系統(tǒng)識別的抗體和抗原的結(jié)構(gòu)越互補,結(jié)合越可能發(fā)生,而把彼此結(jié)合的強度稱之為親和力??乖枰M可能好的與入侵的抗體相結(jié)合。二者匹配程度越好,結(jié)合就越好,抗原和抗體親和力就越大。對于雷達組網(wǎng)問題,抗原對應(yīng)的是雷達組網(wǎng)的最大覆蓋系數(shù),由于雷達組網(wǎng)的總責(zé)任區(qū)域的面積是一定的,可以定義抗體Ab與抗原Ag的之間的親和力為:
(Ab,Ag)=1Tb-Tg(3)
式中:Tb,Ta分別為抗體Ab,Ag對應(yīng)的雷達網(wǎng)的覆蓋系數(shù),其中Tg也是所求的最大的雷達網(wǎng)覆蓋系數(shù)。
定義抗體Ab1與抗體Ab2之間的親和力為:
App(Ab1,Ab2)=1|Tb1-Tb2|(4)
式中:Tb1、Tb2分別為抗體Ab1與抗體Ab2對應(yīng)的雷達網(wǎng)覆蓋系數(shù)[4]。
3.4 算法選擇、交叉、變異算子
免疫遺傳算法能使抗體保持多樣性并且最終能夠收斂到最優(yōu)解的主要操作,就是因為在算法中有選擇、交叉和變異算子的存在,從而使整個抗體群沿著適應(yīng)度較好的方向搜索。
1)選擇算子
選擇是根據(jù)適者生存原則選擇下一代抗體,在基于免疫遺傳算法的雷達組網(wǎng)中選擇的是下一代的組網(wǎng)方案。采用如下的選擇算子:
PS(xi)=αρ(xi)∑ni=1ρ(xi)+(1-α)1Ne-ciβ(5)
式中:ρ(xi)是以式(1)為適應(yīng)度函數(shù)的雷達組網(wǎng)覆蓋系數(shù);ci是抗體xi的濃度,即群體中相似抗體所占的比重[5],即
ci=和抗體i親和度大于λ的抗體數(shù)N(6)
其中λ為相似常數(shù),一般取為0.9~1;
α和β是常數(shù)調(diào)節(jié)因子; N是該種群內(nèi)抗體的總數(shù)。
2)交叉
交叉是在選中的抗體中,對兩個不同的個體按交叉概率Pc隨機選中相同的位置進行基因交換(一般交叉概率取值為0.15~0.75),從而產(chǎn)生新的抗體,也就是新的雷達組網(wǎng)方案。如果對于5臺組網(wǎng)雷達,隨機選擇的是抗體1和抗體2進行交叉,抗體1和抗體2的編碼如下:
抗體1:X11Y11X12Y12X13Y13X14Y14X15Y15
抗體2:X21Y21X22Y22X23Y23X24Y24X25Y25
按照交叉概率產(chǎn)生的交叉點為4,則交叉以后的抗體1和抗體2的編碼分別是:
抗體1:X11Y11X12Y12X23Y23X24Y24X25Y25
抗體2:X21Y21X22Y22X13Y13X14Y14X15Y15
3)變異
變異是在選中的抗體中,對個體中的某些基因以一定的變異概率對某些抗體的某些位執(zhí)行異向轉(zhuǎn)化(一般變異概率的取值為0.01~0.2)。在變異的時候,對于交叉過程中產(chǎn)生的抗體方案,隨機產(chǎn)生一個介于[-MAX/2, +MAX/2]的數(shù)字rand,其中MAX為橫坐標(biāo)和縱坐標(biāo)可取的最大值。變異以后坐標(biāo)值x與原來值x之間的關(guān)系如下:
x,=x+rand(7)
變異以后如果x′>MAX,則取x′=MAX;如果x,<0,則取x′= 0。
單靠變異不能在求解中得到好處,但是,它能保證算法過程不會產(chǎn)生無法進化的單一群體。因為在所有的個體一樣的時候,選擇和交叉無法產(chǎn)生新的基因,這時只能靠變異產(chǎn)生新的基因,即變異增加了全局優(yōu)化的特質(zhì)。
3.5 基于濃度的種群更新
由式(5)可見,本文中的選擇算子是采用基于抗體的雷達組網(wǎng)覆蓋系數(shù)以及抗體濃度的選擇概率Ps(xi),從而有效地保證了抗體的多樣性,避免了算法的早熟問題,能夠獲得更好的覆蓋系數(shù)。
4 仿真實驗
根據(jù)前面的算法流程,選取目標(biāo)函數(shù)為N臺雷達所覆蓋的有效責(zé)任區(qū)域占總責(zé)任區(qū)域的比重。雷達的臺數(shù)分別為3臺,4臺,5臺,初始抗體的數(shù)量為N=50,即一共有50種方案。利用matlab仿真的結(jié)果如圖2所示。圖中矩形框表示的是組網(wǎng)雷達陣地的有效范圍(長為200公里,寬為100公里),雷達的類型一共有兩種,即探測范圍為50公里的雷達和探測范圍為30公里的雷達。參照雷達組網(wǎng)的實際要求,在仿真的過程中,最多只提供2臺性能優(yōu)越的雷達,效能較差的雷達數(shù)量不限。
(a)3臺雷達的面積覆蓋率
(b)4臺雷達的面積覆蓋率
(c)5臺雷達的面積覆蓋率
圖2 matlab仿真結(jié)果
在圖2中當(dāng)雷達數(shù)量為3臺時,組網(wǎng)雷達的面積覆蓋率為87.84%(圖2(a));當(dāng)為4臺時,組網(wǎng)雷達的面積覆蓋率為91.53%(圖2(b));當(dāng)為5臺時,組網(wǎng)雷達的覆蓋率為94.27%(圖2(c))。將仿真的結(jié)果與傳統(tǒng)遺傳算法的組網(wǎng)覆蓋率對比[6],如表1所示。
表1 雷達組網(wǎng)覆蓋率對比
3臺雷達
4臺雷達
5臺雷達
傳統(tǒng)遺傳算法覆蓋率
81.5%
86.1%
88.2%
免疫遺傳算法覆蓋率
87.84%
91.53%
94.27%
由表1可以看出,基于免疫遺傳算法雷達組網(wǎng)的覆蓋率要高很多,將免疫遺傳算法應(yīng)用于雷達組網(wǎng)是一種有效的方法。
5 結(jié) 論
雷達組網(wǎng)能夠有效地提高雷達系統(tǒng)的整體性能,更好的適應(yīng)高科技電子戰(zhàn)、信息戰(zhàn)。將免疫遺傳算法應(yīng)用于雷達組網(wǎng)系統(tǒng),能夠有效地提高雷達網(wǎng)的面積覆蓋率。免疫遺傳算法與傳統(tǒng)的遺傳算法相比,能夠克服傳統(tǒng)遺傳算法的未成熟收斂,提高全局搜索能力。
參考文獻
[1] 彭獲由. 區(qū)域性雷達組網(wǎng)[J]. 電子科學(xué)技術(shù)評論, 1992(3) :1- 6
[2] LUO WENJIAN, CAO XIANBIN,WANG XUFA. An Immune Genetic Algorithm Based on Immune Regulation[A]. Proceedings of the 2002 Congress on Evolutionary ComputationCEC2002. May 12-17, 2002. Honolulu, USA. Vol.2: 801-806.
[3] 段玉波, 任偉建, 霍鳳財, 等. 一種新的免疫遺傳算法及其應(yīng)用[J], 2005, 20(10):1185-1188.
[4] 謝剛, 武斌, 謝克明. 基于免疫遺傳算法的TSP優(yōu)化問題求解[J]. 太原理工大學(xué)學(xué)報, 2007, 38(3): 199-201.
[5] 劉旭旺, 齊微. 免疫遺傳算法的研究及其在TSP中的應(yīng)用[EB/OL]. 中國科技論文在線,http://www.paper.edu.cn/releasepaper/content/200706-566, 2007-06-28
[6] 莫建文, 韓傳久, 張彤. 遺傳算法在雷達組網(wǎng)最優(yōu)化布陣中算法研究[J]. 信息技術(shù), 2005(3): 66-68.
X11Y11X12Y12X13Y13X14Y14X15Y15第1種雷達組網(wǎng)方案
X21Y21X22Y22X23Y23X24Y24X25Y25
第2種雷達組網(wǎng)方案
……
Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5
第n種雷達組網(wǎng)方案
其中,Xn1Yn1Xn2Yn2Xn3Yn3Xn4Yn4Xn5Yn5(n表示雷達組網(wǎng)的任意一種方案)分別是第1臺、第2臺、第3臺、第4臺和第5臺雷達的橫坐標(biāo)和縱坐標(biāo)。雷達布陣的每一種方案對應(yīng)為免疫遺傳算法中抗體種群中的每一個個體。
3.3 親和力計算
親和力是指發(fā)生免疫系統(tǒng)識別的抗體和抗原的結(jié)構(gòu)越互補,結(jié)合越可能發(fā)生,而把彼此結(jié)合的強度稱之為親和力??乖枰M可能好的與入侵的抗體相結(jié)合。二者匹配程度越好,結(jié)合就越好,抗原和抗體親和力就越大。對于雷達組網(wǎng)問題,抗原對應(yīng)的是雷達組網(wǎng)的最大覆蓋系數(shù),由于雷達組網(wǎng)的總責(zé)任區(qū)域的面積是一定的,可以定義抗體Ab與抗原Ag的之間的親和力為:
(Ab,Ag)=1Tb-Tg(3)
式中:Tb,Ta分別為抗體Ab,Ag對應(yīng)的雷達網(wǎng)的覆蓋系數(shù),其中Tg也是所求的最大的雷達網(wǎng)覆蓋系數(shù)。
定義抗體Ab1與抗體Ab2之間的親和力為:
App(Ab1,Ab2)=1|Tb1-Tb2|(4)
式中:Tb1、Tb2分別為抗體Ab1與抗體Ab2對應(yīng)的雷達網(wǎng)覆蓋系數(shù)[4]。
3.4 算法選擇、交叉、變異算子
免疫遺傳算法能使抗體保持多樣性并且最終能夠收斂到最優(yōu)解的主要操作,就是因為在算法中有選擇、交叉和變異算子的存在,從而使整個抗體群沿著適應(yīng)度較好的方向搜索。
1)選擇算子
選擇是根據(jù)適者生存原則選擇下一代抗體,在基于免疫遺傳算法的雷達組網(wǎng)中選擇的是下一代的組網(wǎng)方案。采用如下的選擇算子:
PS(xi)=αρ(xi)∑ni=1ρ(xi)+(1-α)1Ne-ciβ(5)
式中:ρ(xi)是以式(1)為適應(yīng)度函數(shù)的雷達組網(wǎng)覆蓋系數(shù);ci是抗體xi的濃度,即群體中相似抗體所占的比重[5],即
ci=和抗體i親和度大于λ的抗體數(shù)N(6)
其中λ為相似常數(shù),一般取為0.9~1;
α和β是常數(shù)調(diào)節(jié)因子; N是該種群內(nèi)抗體的總數(shù)。
2)交叉
交叉是在選中的抗體中,對兩個不同的個體按交叉概率Pc隨機選中相同的位置進行基因交換(一般交叉概率取值為0.15~0.75),從而產(chǎn)生新的抗體,也就是新的雷達組網(wǎng)方案。如果對于5臺組網(wǎng)雷達,隨機選擇的是抗體1和抗體2進行交叉,抗體1和抗體2的編碼如下:
抗體1:X11Y11X12Y12X13Y13X14Y14X15Y15
抗體2:X21Y21X22Y22X23Y23X24Y24X25Y25
按照交叉概率產(chǎn)生的交叉點為4,則交叉以后的抗體1和抗體2的編碼分別是:
抗體1:X11Y11X12Y12X23Y23X24Y24X25Y25
抗體2:X21Y21X22Y22X13Y13X14Y14X15Y15
3)變異
變異是在選中的抗體中,對個體中的某些基因以一定的變異概率對某些抗體的某些位執(zhí)行異向轉(zhuǎn)化(一般變異概率的取值為0.01~0.2)。在變異的時候,對于交叉過程中產(chǎn)生的抗體方案,隨機產(chǎn)生一個介于[-MAX/2, +MAX/2]的數(shù)字rand,其中MAX為橫坐標(biāo)和縱坐標(biāo)可取的最大值。變異以后坐標(biāo)值x與原來值x之間的關(guān)系如下:
x,=x+rand(7)
變異以后如果x′>MAX,則取x′=MAX;如果x,<0,則取x′= 0。
單靠變異不能在求解中得到好處,但是,它能保證算法過程不會產(chǎn)生無法進化的單一群體。因為在所有的個體一樣的時候,選擇和交叉無法產(chǎn)生新的基因,這時只能靠變異產(chǎn)生新的基因,即變異增加了全局優(yōu)化的特質(zhì)。
3.5 基于濃度的種群更新
由式(5)可見,本文中的選擇算子是采用基于抗體的雷達組網(wǎng)覆蓋系數(shù)以及抗體濃度的選擇概率Ps(xi),從而有效地保證了抗體的多樣性,避免了算法的早熟問題,能夠獲得更好的覆蓋系數(shù)。
4 仿真實驗
根據(jù)前面的算法流程,選取目標(biāo)函數(shù)為N臺雷達所覆蓋的有效責(zé)任區(qū)域占總責(zé)任區(qū)域的比重。雷達的臺數(shù)分別為3臺,4臺,5臺,初始抗體的數(shù)量為N=50,即一共有50種方案。利用matlab仿真的結(jié)果如圖2所示。圖中矩形框表示的是組網(wǎng)雷達陣地的有效范圍(長為200公里,寬為100公里),雷達的類型一共有兩種,即探測范圍為50公里的雷達和探測范圍為30公里的雷達。參照雷達組網(wǎng)的實際要求,在仿真的過程中,最多只提供2臺性能優(yōu)越的雷達,效能較差的雷達數(shù)量不限。
(a)3臺雷達的面積覆蓋率
(b)4臺雷達的面積覆蓋率
(c)5臺雷達的面積覆蓋率
圖2 matlab仿真結(jié)果
在圖2中當(dāng)雷達數(shù)量為3臺時,組網(wǎng)雷達的面積覆蓋率為87.84%(圖2(a));當(dāng)為4臺時,組網(wǎng)雷達的面積覆蓋率為91.53%(圖2(b));當(dāng)為5臺時,組網(wǎng)雷達的覆蓋率為94.27%(圖2(c))。將仿真的結(jié)果與傳統(tǒng)遺傳算法的組網(wǎng)覆蓋率對比[6],如表1所示。
表1 雷達組網(wǎng)覆蓋率對比
3臺雷達
4臺雷達
5臺雷達
傳統(tǒng)遺傳算法覆蓋率
81.5%
86.1%
88.2%
免疫遺傳算法覆蓋率
87.84%
91.53%
94.27%
由表1可以看出,基于免疫遺傳算法雷達組網(wǎng)的覆蓋率要高很多,將免疫遺傳算法應(yīng)用于雷達組網(wǎng)是一種有效的方法。
5 結(jié) 論
雷達組網(wǎng)能夠有效地提高雷達系統(tǒng)的整體性能,更好的適應(yīng)高科技電子戰(zhàn)、信息戰(zhàn)。將免疫遺傳算法應(yīng)用于雷達組網(wǎng)系統(tǒng),能夠有效地提高雷達網(wǎng)的面積覆蓋率。免疫遺傳算法與傳統(tǒng)的遺傳算法相比,能夠克服傳統(tǒng)遺傳算法的未成熟收斂,提高全局搜索能力。
參考文獻
[1] 彭獲由. 區(qū)域性雷達組網(wǎng)[J]. 電子科學(xué)技術(shù)評論, 1992(3) :1- 6
[2] LUO WENJIAN, CAO XIANBIN,WANG XUFA. An Immune Genetic Algorithm Based on Immune Regulation[A]. Proceedings of the 2002 Congress on Evolutionary ComputationCEC2002. May 12-17, 2002. Honolulu, USA. Vol.2: 801-806.
[3] 段玉波, 任偉建, 霍鳳財, 等. 一種新的免疫遺傳算法及其應(yīng)用[J], 2005, 20(10):1185-1188.
[4] 謝剛, 武斌, 謝克明. 基于免疫遺傳算法的TSP優(yōu)化問題求解[J]. 太原理工大學(xué)學(xué)報, 2007, 38(3): 199-201.
[5] 劉旭旺, 齊微. 免疫遺傳算法的研究及其在TSP中的應(yīng)用[EB/OL]. 中國科技論文在線,http://www.paper.edu.cn/releasepaper/content/200706-566, 2007-06-28
[6] 莫建文, 韓傳久, 張彤. 遺傳算法在雷達組網(wǎng)最優(yōu)化布陣中算法研究[J]. 信息技術(shù), 2005(3): 66-68.