劉 振,徐學文,劉 勇

(海軍航空大學 岸防兵學院,山東 煙臺 264001)

空戰(zhàn)決策問題本質(zhì)上是一個目標分配問題,其目的在于將目標合理地分配給攻擊飛機,使得機群的整體效能最大。當前有眾多的研究者利用各種的智能優(yōu)化算法求解目標分配問題,進化算法在尋優(yōu)的簡便性、通用性以及全局尋優(yōu)性上的優(yōu)勢,使得其應用廣泛,例如文獻[1-2]就分別利用遺傳算法和文化基因算法等優(yōu)化方法進行求解,也逐步提出一些改進類型的進化算法應用于目標分配問題,例如文獻[3-4]分別利用精英改選機制粒子群算法和啟發(fā)式粒子群算法求解空戰(zhàn)決策問題?，F(xiàn)在的作戰(zhàn)環(huán)境逐步趨向于分布式網(wǎng)絡化,因此諸多文獻都采用了分布式的求解方法,文獻[5]利用協(xié)同拍賣求解目標分配問題。合同網(wǎng)也逐步應用到目標和任務的分配中[6],基本思想都是利用合同網(wǎng)機制中的拍賣合同、競標合同以及交換和置換合同等方法處理靜態(tài)和動態(tài)的分配問題。

當前已經(jīng)提出了諸多方法用于求解空戰(zhàn)決策和目標分配問題,隨著網(wǎng)絡化作戰(zhàn)環(huán)境的逐步建立和完善,多作戰(zhàn)平臺在戰(zhàn)場環(huán)境下能夠?qū)崿F(xiàn)態(tài)勢的共享。因此,當某架飛機發(fā)射導彈后,若該機失去制導能力或者需要做機動規(guī)避的情形下,如何使得另一架飛機有效地實現(xiàn)對導彈的協(xié)同制導是一個需要解決的問題。當前較多文獻關注多機協(xié)同制導下的制導律設計[7],關于協(xié)同制導下的目標分配和任務分配研究的相對較少,文獻[8]較早關注到了該問題,并提出了協(xié)同制導決策方法,文獻[9-10]進一步拓展推廣了制導優(yōu)勢模型,以此為基礎建立了考慮協(xié)同制導的任務分配模型,所考慮的制導交接因素過多注重于當前空戰(zhàn)態(tài)勢因素。

本文針對多機協(xié)同制導情況下的決策模型和求解方法進行探討和分析。在本文的分析過程中,將發(fā)射導彈的飛機稱為交班平臺,將協(xié)同制導的飛機稱為接班平臺。在危險的戰(zhàn)場環(huán)境下,交班平臺在發(fā)射完導彈后可以進行有效的規(guī)避,而制導的任務則可以交給有制導能力的接班平臺,從而避免人員的傷亡,提高整體作戰(zhàn)效能。通過查閱以及分析國內(nèi)外文獻,本文設計了一種協(xié)同制導交接方案,以敵方毀傷最大以及我方損失最小,建立了考慮協(xié)同制導下的空戰(zhàn)決策模型,該模型不僅充分考慮態(tài)勢優(yōu)勢,同時也將截獲概率和交接成功概率作為候選制導平臺的重要因素,為提高求解該模型的收斂速度和精度,提出一種分布估計免疫算法(estimation distribution immune algorithm,EDIA)來求解該模型。

1 協(xié)同制導空戰(zhàn)決策問題分析及建模

雙機或多機協(xié)同制導就是在分布式網(wǎng)絡化的作戰(zhàn)條件下,當交班平臺無法有效地對所發(fā)射的導彈進行制導時,將空空導彈的中制導權交給制導平臺,協(xié)同制導示意圖如圖1所示。

從圖1可以看到,交班平臺在發(fā)射完導彈后,由于有其他作戰(zhàn)任務或者本身受損,無法繼續(xù)完成所發(fā)射導彈的制導任務,因此將中制導任務進行交接,接班平臺可利用中繼制導平臺,如預警機提供的探測信息,完成對導彈的制導任務。

美國于二十世紀七八十年代就開始對協(xié)同作戰(zhàn)能力系統(tǒng)進行研究,但在國外公開發(fā)表的資料中尚未看到空戰(zhàn)中相關的研究。因此,研究協(xié)同制導下的空戰(zhàn)決策問題具有重要的現(xiàn)實意義。本文對這個問題進行了初步探討分析,提出了協(xié)同制導下的決策模型和求解方法。

1.1 協(xié)同制導接班平臺確定方法

首先進行如下的假設:

①假設1。 交班平臺和接班平臺距離小于接班平臺的最大制導距離,即接班平臺能夠?qū)话嗥脚_發(fā)射的導彈進行有效的制導;

②假設2。 導彈在發(fā)射后,直至命中目標,都處于接班平臺的制導范圍內(nèi)。

文獻[8]將接班平臺優(yōu)勢函數(shù)作為制導交接的依據(jù),但優(yōu)勢函數(shù)最大的接班平臺,并一定能夠有效地對導彈進行制導。因此,本文給出一種多機協(xié)同制導情況下接班平臺的選擇方法,在考慮制導優(yōu)勢外,將導彈導引頭截獲目標概率和交接成功概率也作為接班平臺的選擇依據(jù)。

1.1.1 交接制導優(yōu)勢

交接制導優(yōu)勢主要由態(tài)勢優(yōu)勢和效能優(yōu)勢構(gòu)成。

1)態(tài)勢優(yōu)勢。

空戰(zhàn)態(tài)勢優(yōu)勢可以表示為

TSI=β1Ta+β2Td+β3Tv

(1)

式中:Ta,Td和Tv和分別為角度、距離和速度優(yōu)勢;β1,β2和β3分別為3種態(tài)勢優(yōu)勢的權重系數(shù),并且β1+β2+β3=1,一般可取0.2,0.4和0.4。

2)空戰(zhàn)效能優(yōu)勢。

(2)

式中:Ci和Cj分別為交班平臺和候選接班平臺的空對空作戰(zhàn)能力。

在得到空戰(zhàn)態(tài)勢優(yōu)勢和空戰(zhàn)效能優(yōu)勢后,總的制導優(yōu)勢為

w=a1TSI+a2TC

(3)

式中:a1和a2通常由專家評估或者根據(jù)經(jīng)驗知識得到,分別為態(tài)勢優(yōu)勢和效能優(yōu)勢對制導優(yōu)勢的影響權重系數(shù),本文設定a1=0.6,a2=0.4。

1.1.2 導引頭截獲目標概率

協(xié)同制導的目標就是使導彈能夠命中目標,提高整體效能,因此制導交接需要考慮導彈的導引頭截獲目標的概率。只有成功截獲目標,才能實現(xiàn)對目標的有效制導[11]。導引頭截獲目標概率可按下式獲得:

Pm=PdPaPv

(4)

式中:Pm為導引頭對目標的正確截獲概率;Pd,Pa和Pv分別為距離、角度和速度截獲概率。

1.1.3 交接成功概率

在滿足協(xié)同制導的前提條件下,除了候選制導平臺自身的優(yōu)勢以外,還必須考慮交班平臺與接班平臺的交接成功概率,優(yōu)勢較大的飛機如果不能有效地捕獲導彈,優(yōu)勢即使再大也是沒有意義的。交接成功概率PH由3個方面的因素構(gòu)成[12],分別為:導彈落入概率PD、對導彈的發(fā)現(xiàn)概率PF以及對導彈的鎖定概率PL。不考慮交班設備的可靠性,交接成功概率為

PH=PDPFPL

(5)

通過上述的分析,接班平臺的選擇概率為

PSE=PmPHw

(6)

1.2 協(xié)同制導交接方案設計

在確定了接班平臺以后,針對目前廣泛使用的慣導+指令修正中制導和主動雷達末制導的空空導彈,給出一種協(xié)同制導交接方案:

①導彈發(fā)射后轉(zhuǎn)入中制導段,交班平臺對導彈進行穩(wěn)定跟蹤,同時接收接班平臺傳送的目標信息,經(jīng)時空一致性轉(zhuǎn)換后向?qū)棸l(fā)送中制導指令;

②交班平臺交接形勢判斷模塊在導彈飛行過程中實時進行判斷計算,當滿足制導交接條件時啟動交接程序。交班平臺向接班平臺發(fā)送“開始交接”信號和導彈的飛行空域預測信息;

③接班平臺接收“開始交接”信號和報文,控制機載雷達轉(zhuǎn)向?qū)楊A測飛行空域,由接班平臺的接力交接控制模塊控制自身機載雷達向?qū)棸l(fā)送一次詢問信號,并等待接收應答信號;

④導彈分時接收兩平臺的制導指令,并發(fā)送應答信號;

⑤接班平臺交接狀態(tài)判斷模塊依據(jù)接收的應答信號對交接狀態(tài)進行實時判斷,判斷交接成功后,向交班平臺發(fā)送“交接成功”信號,交班平臺接收“交接成功”信號,停止向?qū)棸l(fā)送指令,隨后接班平臺開始在詢問信號后加入制導信息,以固定的頻率向?qū)棸l(fā)送制導信息。

按照協(xié)同制導交接方案,交班平臺就可以把導彈的制導權順利地移交給接班平臺。

1.3 協(xié)同制導空戰(zhàn)決策模型建立

①每個目標至多分配一枚導彈;

②每枚導彈只能用于攻擊一個目標。

如果此時Ri的機載導彈Mk需要進行協(xié)同制導,Ri發(fā)射后并不對Mk進行制導。在協(xié)同制導的情形下Mk擊毀Bj的概率Tkj并不等于Ri對Bj的威脅Tij,需要考慮制導交接的成功概率,故此時Tkj=PHTij。

則藍方損失的威脅可以表示為

(7)

式中:ω為導彈目標分配方案。

滿足約束條件式(8)和式(9):

(8)

(9)

當?shù)趉枚導彈用于攻擊第j個目標時,則Xkj=1。如果在各架飛機掛載的導彈數(shù)目都為S時,通過觀察該方程可以看出,當i=1,2,…,M,式(7)中k的取值范圍只能為Ui=[iS(i+1)S-1],式(7)可以表示為

(10)

(11)

(12)

在考慮消除威脅最大的同時,還必須考慮我方受到的威脅最小,因Ri對Bj的威脅為Tij,則Bj對Ri的威脅為Tji,Tji可以看作Bj對Ri的擊毀概率,則紅方受到的總威脅為

(13)

由于飛機一般都攜帶有欺騙設備和干擾設備,從而可以有效地降低其所受到的威脅,軟殺傷武器系統(tǒng)可以等效為欺騙設備,假定Ri等效后的欺騙設備的數(shù)量為mi,因此考慮欺騙設備和干擾設備后的總威脅為

(14)

因此總的空戰(zhàn)決策評估函數(shù)為

R(ω)=E(ω)-Q(ω)

(15)

2 協(xié)同制導空戰(zhàn)決策模型求解算法

協(xié)同制導空戰(zhàn)決策問題對求解速度和求解精度有較高的要求,智能優(yōu)化方法由于其強大的尋優(yōu)能力以及操作上的方便性,已經(jīng)被廣泛用于求解該類問題。免疫算法是受生物免疫系統(tǒng)啟示的一種尋優(yōu)搜索算法,具有良好的學習記憶能力和優(yōu)良的全局尋優(yōu)性能,已經(jīng)被廣泛研究和應用[13-14],本文考慮利用免疫算法求解空戰(zhàn)決策問題。但通過對免疫算法的研究發(fā)現(xiàn),基本免疫算法難免存在尋優(yōu)效率不高,尋優(yōu)效果不佳等問題。擴展緊致遺傳算法[15]作為一種進化性能優(yōu)越的分布估計算法,利用優(yōu)良解集的概率模型指導進化的方向,通過將染色體基因位劃分為MPM,從而保證優(yōu)良模式不在進化過程中被破壞。因此為提高免疫算法的求解性能,本文將分布估計的思想融入到免疫算法中,提出分布估計免疫算法(EDIA),基于MPM進行交叉和變異操作,可以有效地提高算法的進化效率,避免算法過早收斂。

2.1 編碼及種群初始化方法

1)編碼方法。

編碼方式如圖2所示。

染色體共有M段組成,染色體長度L=MN,染色體中第i(1≤i≤M)部分的第j(1≤j≤N)位為1,則代表第i架紅機能夠用于攻擊第j架藍機。

2)種群初始化方法。

在算法的運行過程中引入啟發(fā)式信息,啟發(fā)式信息是指在算法的運行開始階段,根據(jù)問題的特點、以往的經(jīng)驗以及實際的資源和環(huán)境條件,在算法的初始化階段所構(gòu)造出的啟發(fā)信息,例如當?shù)趇架紅方飛機不能用于攻擊第j架藍方飛機時,則可將該染色體相應的位置設為0。在算法的運行過程中,根據(jù)問題的實際特點,如果產(chǎn)生不可行解,人工修改染色體中相應的位置,調(diào)節(jié)算法的搜索范圍和方向。

隨機的初始化種群不便于種群個體進化,且有陷入局部極值的風險,為了使種群能保持較好的多樣性,文獻[16]利用反向?qū)W習方法進行種群的初始化,證明能夠有效地提高算法的進化效率,本文利用文獻[16]提出的方法進行種群初始化。

2.2 基于MPM的交叉和變異

在進化算法中,利用交叉方式可以有效地提高種群多樣性,獲得全局優(yōu)良解,但交叉方式選擇至關重要,優(yōu)良的交叉方式可以有效地保證積木塊的累積,因此本文采用文獻[15]提出的基于MPM的交叉。將MPM作為一個整體進行交叉,從而可以有效地保證種群個體基因位之間的鏈接關系,使得優(yōu)良模式得到保留和進化。同時本文也采用了MPM變異,即按照變異概率Pu,對某一MPM中的基因位進行變異操作。在種群進化過程中,如果整個進化過程中采用統(tǒng)一的變異概率,收斂過程緩慢,不便于優(yōu)良個體進化,利用自適應變異概率,第i個個體的變異概率設置為

(16)

2.3 免疫選擇方法

為保證進化種群的多樣性,避免單一選擇方式導致的進化停滯和陷入局部極值,利用適應度共享增強種群進化的多樣性,對于進化種群中的抗體ai(t)和aj(t),其共享函數(shù)可以定義為

(17)

式中:σs為共享半徑,dij為抗體之間的距離,ai(t)的適應度可以表示為

(18)

(19)

2.4 算法求解流程

利用本文提出的分布估計免疫算法求解協(xié)同制導下的決策模型,首先需要初始化進化參數(shù)信息,包括種群規(guī)模N,循環(huán)迭代次數(shù)Tmax,依據(jù)敵我雙方飛機數(shù)目確定的染色體長度L,共享半徑σs,變異概率Pu,其主要流程可概述為:

①步驟1。 依據(jù)當前態(tài)勢信息和空戰(zhàn)決策模型,初始化進化算法參數(shù)信息;

②步驟2。 利用反向?qū)W習方法,初始化種群P(t),并根據(jù)初始啟發(fā)式信息剔除不合理個體;

③步驟3。 對種群P(t)中的染色體分別進行交叉和變異操作,其中交叉和變異均按照MPM方式進行,變異概率按照式(16)設置;

④步驟4。 計算抗體和抗原之間的親和力,依據(jù)式(18)進行適應度共享操作,并依據(jù)式(19)的選擇概率進行免疫選擇;

⑤步驟5。 判斷是否滿足結(jié)束條件,滿足則結(jié)束,輸出分配結(jié)果,否則轉(zhuǎn)步驟2。

3 協(xié)同制導空戰(zhàn)決策仿真分析

在標準態(tài)勢和隨機態(tài)勢下進行仿真分析,其示意圖分別如圖3(a)和圖3(b)所示,標準態(tài)勢的數(shù)據(jù)來自文獻[8],用以對比分析本文采用的制導平臺選擇方法與文獻[8]的區(qū)別聯(lián)系,隨機態(tài)勢用以仿真分析本文提出的空戰(zhàn)決策模型及分布估計免疫算法的正確性。

標準空戰(zhàn)態(tài)勢如圖3(a)所示,其中B、C和D3架飛機對攻擊目標T的某枚導彈進行協(xié)同制導,得到3架飛機的態(tài)勢優(yōu)勢分別為0.627 7,0.412 0和0.326 0,根據(jù)文獻[8],效能優(yōu)勢都為0.75,因此3架飛機總的制導優(yōu)勢為0.676 6,0.547 2和0.495 6。依據(jù)文獻[8],若單純考慮制導優(yōu)勢,應該選擇平臺B作為接班平臺,但從實際空戰(zhàn)情況出發(fā),為有效提高作戰(zhàn)效能,確保制導平臺能夠順利截獲目標并完成對導彈的協(xié)同制導,除了考慮制導優(yōu)勢外,還應考慮目標截獲概率和制導交接成功概率。若B,C和D對需要協(xié)同制導的導彈進行制導后,雷達截獲概率分別為0.90,0.99和0.99,交接成功概率分別為0.924,0.975和0.988,則按照本文計算得到的制導平臺選擇概率分別為0.522 0,0.528 2和0.484 8,故此時應該選擇平臺C作為制導接班平臺。

由此可以看出,當考慮了雷達截獲概率和交接成功概率后,與單純依靠制導平臺優(yōu)勢進行選擇發(fā)生了沖突。經(jīng)分析發(fā)現(xiàn),單純考慮制導優(yōu)勢,并不能綜合考慮制導交接后的空戰(zhàn)效能變化情況,由于平臺B的雷達截獲概率及交接成功概率較低,因此不宜作為接班平臺,平臺C的制導優(yōu)勢雖然不如平臺B高,但其雷達截獲概率和交接成功概率有優(yōu)勢,因此應以平臺C為接班平臺。在復雜的戰(zhàn)場環(huán)境下,單純依靠制導優(yōu)勢進行選擇過于片面,依據(jù)本文的方法,當考慮了影響制導交接的多個約束條件后,能夠更加全面地刻畫此時交接的情況。

隨機交戰(zhàn)態(tài)勢如圖3(b)所示,紅方R飛機數(shù)量M=4,每架飛機可以掛載4枚導彈,藍方B飛機數(shù)量N=14。紅方R此時的總導彈數(shù)目為14,R1,R2,R3和R4掛載的導彈序號分別為1～4,5～8,9～12,13～14。R3由于在前一波次的攻擊中發(fā)射了2枚導彈,因此此時只有2枚導彈,并且其制導系統(tǒng)受到損壞,不能對其發(fā)射的導彈進行中制導。

在隨機態(tài)勢下,R4由于戰(zhàn)術原因需要對其掛載的導彈進行協(xié)同制導。依據(jù)1.1節(jié)提出的方法,最終選擇R3為其接班平臺,R4將第13和14枚導彈的中制導權交給R3,最終計算得到的分配方案為

R1{1,2,3,4}→{B13,B7,B12,B9}R2{5,6,7,8}→{B3,B6,B14,B5}R3{9,10,11,12,13,14}→{B2,B8,B10,B11,B1,B4}

其中左邊括號內(nèi)為紅機Ri制導的導彈編號,右邊括號內(nèi)為Ri分配的藍機編號,從最終的分配方案可以看出,R4掛載的第13和14枚導彈被用于攻擊B1和B4。

將采用協(xié)同制導后的效能以及不同求解算法獲得的效能進行對比,仿真結(jié)果如圖4所示。圖4(a)為采用協(xié)同制導前后的效能對比圖,其中橫坐標為迭代次數(shù)(t),縱坐標為效能值(f),兩者均無量綱。從中可以看出采用協(xié)同制導后的效能值有了明顯提高,這是因為若不采用協(xié)同制導,則R3必須退出攻擊,從而造成了資源的浪費和效能的降低。

圖4(b)為在隨機態(tài)勢下,分別用遺傳算法(GA)、免疫遺傳算法(IGA)以及分布估計免疫算法(EDIA)協(xié)同制導決策問題進行對比,從圖中可以發(fā)現(xiàn),本文提出的算法能有效提高整體作戰(zhàn)效能,并且本文提出的EDIA收斂速度也優(yōu)于GA和IGA。

4 結(jié)束語

本文分析協(xié)同制導下的空戰(zhàn)決策方法,給出了協(xié)同制導的接班平臺選擇方法和協(xié)同制導交接方案設計,建立了協(xié)同制導下的空戰(zhàn)決策模型。為提高問題的求解效率,提出利用分布估計免疫算法求解協(xié)同制導下的空戰(zhàn)決策問題,對本文提出的方法和模型進行了仿真分析,證明了本文模型和方法的正確性。

對協(xié)同制導下的空戰(zhàn)決策問題,現(xiàn)在開展的研究還相對較少,因此在下一步的研究中考慮如何在復雜環(huán)境下進行協(xié)同制導決策,特別是動態(tài)目標的協(xié)同制導空戰(zhàn)決策,是進一步研究的方向。

[1] CHEN Jie,XIN Bin,PENG Zhihong,et al. Evolutionary decision-making for the dynamic weapon-target assignment problem[J]. Science China:Series F-information Science,2009,52(11):2006-2018.

[2] 王曉光,章衛(wèi)國,陳偉. 無人機編隊超視距空戰(zhàn)決策及作戰(zhàn)仿真[J]. 控制與決策,2015,28(2):352-357.

WANG Xiaoguang,ZHANG Weiguo,CHEN Wei. BVR air combat decision making and simulation for UAV formation[J]. Control and Decision,2015,28(2):352-357.(in Chinese)

[3] 王昱,章衛(wèi)國,傅莉,等. 基于精英改選機制的粒子群算法的空戰(zhàn)納什均衡策略逼近[J]. 控制理論與應用,2015,32(7):857-865.

WANG Yu,ZHANG WeiGuo,FU Li,et al. Nash equilibrium strategies approach for aerial combat based on elite re-election particle swarm optimization[J]. Control Theory and Application,2015,32(7):857-865.(in Chinese)

[4] LUO Delin,YANG Zhong,DUAN Haibin,et al. Heuristic particle swarm optimization algorithm for air combat decision-making on CMTA[J]. Transactions of Nanjing University of Aeronautics & Astronautics,2006,23(1):20-26.

[5] EDALAT N,XIAO Wendong,ROY N,et al. Combinatorial auction-based task allocation in multi-application wireless sensor networks[C]//The 9th International Conference on Embedded and Ubiquitous Computing. Melbourne:IEEE,2011:3622-3627.

[6] 廖沫,陳宗基. 基于多Agent分布協(xié)同拍賣的動態(tài)目標分配算法[J]. 北京航空航天大學學報,2007,33(2):180-183.

LIAO Mo,CHEN Zongji. Dynamic target assignment method based on multi-agent decentralized cooperative auction[J]. Journal of Beijing University of Aeronautics and Astronautics,2007,33(2):180-183.(in Chinese)

[7] ZHAO Jiang,ZHOU Rui,DONG Zhuoning. Three-dimensional cooperative guidance laws against stationary and maneuvering targets[J]. Chinese Journal of Aeronautics,2015,28(4):1104-1120.

[8] 肖冰松,方洋旺,胡詩國,等. 多機空戰(zhàn)協(xié)同制導決策方法[J]. 系統(tǒng)工程與電子技術,2009,31(3):610-612.

XIAO Bingsong,FANG Yangwang,HU Shiguo,et al. Decision methods for cooperative guidance in multi-aircraft air warfare[J]. Systems Engineering and Electronics,2009,31(3):610-612.(in Chinese)

[9] 冉華明,周銳,吳江. 超視距空戰(zhàn)中多機協(xié)同制導方法[J]. 北京航空航天大學學報,2014,40(10):1457-1462.

RAN Huaming ZHOU Rui,WU Jiang,et al. Cooperative guidance of multi-aircraft in beyond-visual-range air combat[J]. Journal of Beijing University of Aeronautics and Astronautics,2014,40(10):1457-1462.(in Chinese)

[10] 周德云,楊振,張堃. 多UCAV超視距協(xié)同空戰(zhàn)中的交接制導方法[J]. 彈道學報,2017,29(2):1-7.

ZHOU Deyun,YANG Zhen,ZHANG Kun. Method of guidance handover in beyond-visual-range coordinated air-combat for multi-UCAVs[J]. Journal of Ballistics,2017,29(2):1-7.(in Chinese)

[11] 樊會濤. 復合制導空空導彈截獲目標概率研究[J]. 航空學報,2010,31(6):1225-1229.

FAN Huitao. Study on target acquisition probability of air-to-air missiles with combined guidance[J]. Acta Aeronautica et Astronautica Sinica,2010,31(6):1225-1229.(in Chinese)

[12] 劉興堂,周自全,李為民,等. 現(xiàn)代導航、制導與測控技術[M]. 北京:科學出版社,2010:173-176.

LIU Xingtang,ZHOU Ziquan,LI Weimin,et al. Modern navigation,guidance and TT&C[M]. Beijing:Science Press,2010:173-176.(in Chinese)

[13] SIMONE S F. Artificial immune algorithm applied to distribution system reconfiguration with variable demand[J]. Electrical Power and Energy Systems,2016,82:561-568.

[14] SHANG Ronghua,DU Bingqi,MA Hongna,et al. Immune clonal algorithm based on directed evolution for multi-objective capacitated arc routing problem[J]. Applied Soft Computing,2016,49:748-758.

[15] 劉振,胡云安,彭軍. 協(xié)同進化擴展緊致量子進化算法[J]. 控制與決策,2014,29(2):320-326.

LIU Zhen,HU Yunan,PENG Jun. Coevolutionary quantum evolutionary algorithm based on extended compact[J]. Control and Decision,2014,29(2):320-326.(in Chinese)

[16] RAHNAMAYAN S,TIZHOOSH H R. A novel population initialization method for accelerating evolutionary algorithms[J]. Computers and Mathematics with Applications,2007,53(10):1605-1614.

[17] CHENG Lijun,DING Yongsheng,HAO Kongrong. An ensemble kernel classifier with immune clonal selection algorithm for automatic discriminant of primary open-angle glaucoma[J]. Neurocomputing,2012,83(15):1-11.