吳 濤 王 迅
(91404部隊93分隊 秦皇島 066001)
干擾資源優化分配在艦艇光電對抗作戰中有非常重要的作用。它是根據裝備的干擾資源以及戰術要求,合理分配干擾資源,以達到最佳的干擾效果,干擾資源優化分配問題可借助于目標優化理論,將問題抽象為尋求系統最優的干擾資源優化分配模型。文中利用光電信息系統作戰運籌中目標分配原則,建立了干擾資源優化分配模型,并應用量子遺傳算法對該模型進行了解算研究。
1)最先告警,最先干擾
對于最先出現告警的目標,盡量在遠界進行致盲、致眩等干擾,減少對己方的威脅。
2)目標最近,最先干擾
對于距離我防區最近目標,要優先干擾。
3)重點目標,優先干擾
對于上級指定的重點目標或者威脅程度大的目標要優先進行干擾。
4)整體最優分配
第一,系統干擾效能最大;第二,對目標的總威脅最大;第三,干擾目標數最多;第四,己方受損最小。
干擾資源分配是把空中的威脅目標分配到某個干擾單元的某個目標通道。由決策變量及各影響因素來決定資源分配優化目標函數及各項約束條件,進行干擾資源優化分配解算,得到分配結果,為干擾決策提供依據。假定威脅目標是由多種類型目標群組成,包括飛機、導彈等共r批目標。光電對抗系統共有n個干擾措施單元,在動態分配過程中,若某個干擾措施單元正在對某批目標進行干擾,則此干擾單元不能馬上使用。干擾資源與目標不是一對一得關系,有可能一個目標分配有多個干擾單元,也可能沒有;可能一個干擾單元同時對抗多批目標,也可能只能對抗一批目標。選擇xij作為決策變量,其含義為:xij=1,表示第i個干擾單元對第j個目標進行干擾;xij=0,表示第i個干擾單元不對第j個目標進行干擾;依據光電信息系統干擾資源整體最優分配原則,取系統總的干擾效能作為目標函數,目的是使系統總的干擾效能達到最大,干擾目標的總威脅最大,干擾目標數最多。建立的干擾資源優化分配數學模型如下[2,3,5~7]:
目標函數:

式中:C為光電對抗系統總的干擾效能;K為干擾目標數;pij為第i個干擾單元對第j個目標干擾有效率;Rj為第j個目標威脅程度權值;Sj為第j批目標是否受到干擾,若受到干擾則Sj=1,否則Sj=0。Sj=x1j∩x2j∩…∩xnj。
量子遺傳算法(quantum genetic algorithm,QGA)是量子計算理論和遺傳算法原理相結合的產物。主要以量子理論和量子計算為基礎,采用量子比特實現染色體編碼,通過量子門對其進行更新,產生種群的多樣性。QGA具有種群規模小、尋優能力強、收斂速度快和計算時間短的特點。
1)量子比特編碼
在量子信息論中,信息的載體不再是經典的比特,而是量子比特或量子位。量子比特可以處于0和1這兩個基態的任意疊加狀態。一個量子計算比特可以表示為:

其中,α和β是兩個復數,分別表示狀態|0〉和狀態|1〉的概率幅。|α|2和|β|2別表示量子比特處于|0〉和|1〉的概率。
一個m位量子比特的編碼形式如下:

2)量子旋轉門
量子旋轉門是演化操作的執行機構,其調整操作如下式:

后的形式,θ為量子門的旋轉角。
1)染色體編碼
由干擾資源優化分配模型解xij的含義,對染色體進行合理編碼,采用二進制字符串表示。個體的串長為n×m(n個干擾單元,m個威脅目標),用xij表示串中的項,代表第i個干擾單元對第j個目標的分配情況。染色體可表示為:x11x12…x1mx21x22…x2m…xn1xn2…xnm。
2)種群初始化
針對光電信息作戰系統中的n個干擾單元和m個威脅目標,根據系統的先驗知識及約束條件,分析最優解子空間在整個問題空間中的分布情況,將最優解所在子空間設定為初始群體生產區域,或者將整個可行域設定為初始群體的生成區域,再隨機生成一定數目的個體,從中挑選出較好的個體構成QGA的初始種群。
3)適應度函數
干擾資源優化分配的目的是最大限度地干擾空中威脅目標,將系統多目標函數進行組合,組合后函數可作為適應度函數,即:

4)量子旋轉門更新策略
通過量子旋轉門,根據式(7)更新種群的量子位,采用一種通用的旋轉角調整策略,如下式所示:

式中:s(αi,βi)是搜索方向函數,Δθ為旋轉角度,其大小影響算法的收斂速度;如果太大,算法易出現早熟現象而陷入局部最優解;如果太小,可能出現停滯狀態,因此,需要自適應調整搜索。本文Δθ取10e-t/maxt,t為進化代數,maxt為最大進化代數,主要使算法向最優解方向進行搜索。
5)終止條件
QGA是一種反復迭代的搜索算法,通過多次進化逐漸逼近最優解。可采用規定遺傳代數與個體適應度相結合的方法作為終止條件,即經過t′次迭代后,所選出的最優個體的適應度不再增加,表明全局最優個體從第t′代開始出現,則終止遺傳搜索。該個體即為所獲得的最優個體,亦即干擾資源優化分配的最優方案。
假設某系統內共有4個干擾單元,10批威脅目標。其中,干擾單元1只能分配給目標3或5或6,干擾單元3只能分配目標2或8或9,干擾單元2或4只能最多分配給序號連續的3批目標。干擾資源分配,就是使整個系統獲得總的最佳干擾效果。干擾有效率(pij)及威脅程度權值(Rj)分別如表2、表3所示,N=60,t=100。

表2 干擾有效率(pij)

表3 威脅程度權值(Rj)
經過Matlab調試和運行,得到局部最優解(x11x12…x1m x21…x2m…xn1xn2…xnm)=[0,0,0,0,1,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,1,1,0,0],對應的適應度值為9.995,每代種群中個體的適應度與迭代次數之間的關系如圖1所示。

圖1 每一代最優適應度和迭代次數之間的關系曲線
從圖中可以看出,當迭代到60后,搜索到局部最優解,文獻[4]中在180次迭代后得到最優解,且適應度為9.941,比較文中的基于量子遺傳算法得到的最優解而言,在搜索速度和適應度兩方面具備一定優勢。故采用量子遺傳算法的明顯優于文獻中的遺傳算法。
文中通過建立光電干擾資源優化分配模型,提出了將量子遺傳算法應用于優化模型的求解,為解決光電信息系統作戰決策提供了新的思路,仿真結果證明該方法是有效、可行的。
[1]ZHANG Ye-peng.Military Operations Basis[M].Beijing:Higher Education Publishing Company,2006
[2]Liu Ming,Li Wei-min,Wang Ying-long,et al.Optimization of the regional air defense disposition based on genetic algorithms[J].system engineering and electronic technology,2003,25(2):191~193
[3]ZHOU Ming,SUN Shu-dong.Genetic Algorithms Principle and Aplication[M].Beijing:National Defense Industry Press,2000
[4]王鐵紅,呂躍廣,劉珍陽,等.基于遺傳算法的光電干擾資源優化分配研究[J].紅外與激光工程,2007,36(6):963~967
[5]Xin Qing-wu,Yang Can-jun.A method of optimizing artillery firing based on fuzzy multiobjective decision making[J].Operations Research and Management Science,2005,14(4):29~34
[6]Yao Xiao-bai,Zhao Xiao-zhe.Weapon-target assignment model based on the decision-aid system for the ship's air defense[J].Fire control and Command control,2005,30(4):10~13
[7]Zhang Jing-yuan.Application of genetic algorithm in parameter identification based on fuzzy model[J].Computer Engineering and Design,2006,27(2):262~264
[8]吳濤,葉曉慧,王紅霞,等.基于量子遺傳算法測試選擇問題的研究[J].計算機測量與控制,2010,18(11)
[9]趙知勁,彭振,鄭仕鏈,等.基于量子遺傳算法的認知無線電頻譜分配[J].物理學報,2009,2(2):1358~1359
[10]吳平健,李智勇.目標分配的遺傳算法改進研究[J].計算機與數字工程,2009,37(9)
[11]羅紅明,王家映,朱培民,等.量子遺傳算法在大地電磁反演中的應用[J].地球物理學報