摘要： 傳統(tǒng)多目標(biāo)優(yōu)化算法在解決多于兩個(gè)目標(biāo)函數(shù)的火力分配問題時(shí)收斂效果不佳，多樣性差，耗時(shí)過大?；诖?，提出了一種自適應(yīng)網(wǎng)格多目標(biāo)鯨魚優(yōu)化算法（AG-MOWOA）來解決以震塌比例、彈藥成本和自身剩余價(jià)值為目標(biāo)函數(shù)的火力分配問題。該算法引入混沌映射和外部Pareto存檔進(jìn)化策略提高了種群的多樣性，通過自適應(yīng)網(wǎng)格選取最優(yōu)個(gè)體的方法極大地減少了算法運(yùn)行時(shí)間。仿真實(shí)驗(yàn)結(jié)果表明，該算法較其他算法收斂速度更快、收斂質(zhì)量更高、解集分布更多樣，能夠有效解決火力分配問題。

關(guān)鍵詞： 火力分配； 混沌映射； 自適應(yīng)網(wǎng)格劃分； 多目標(biāo)優(yōu)化； 鯨魚優(yōu)化算法

中圖分類號(hào)： TP301.6;E91

文獻(xiàn)標(biāo)志碼： A

文章編號(hào)： 1671-6841（2024）06-0017-08

DOI： 10.13705/j.issn.1671-6841.2023010

Research on Firepower Assignment Problem Based on Adaptive

Grid Multi-objective Whale Algorithm

SHE  Wei1，2，3， WANG Yeteng1，3， KONG Defeng4， LIU Wei1，2，3， LI Yinghao1，2， TIAN Zhao1，3

（1.School of Cyber Science and Engineering， Zhengzhou University， Zhengzhou 450002， China;

2.Henan Collaborative Innovation Center for Internet Medical and Health Services， Zhengzhou University，

Zhengzhou 450052， China; 3.Zhengzhou Key Laboratory of Blockchain and Data Intelligence，

Zhengzhou 450002， China; 4.Institute of Engineering Protection， National Defense Engineering

Research Institute， The Academy of Military Sciences， Luoyang 471023， China）

Abstract： The traditional multi-objective optimization algorithm had poor convergence effect， bad diversity and serious time-consuming problems when solving the firepower assignment problem with more than two objective functions. Based on this situation， an adaptive grid multi-objective whale optimization algorithm （AG-MOWOA） was proposed. The algorithm was to solve the firepower assignment optimization problem in which the collapse ratio， ammunition cost， and the own surplus value were taken as the objective functions. Besides， the algorithm increased the population diversity by introducing chaotic mapping and using external Pareto archival evolutionary strategy， and reduced the algorithm running time greatly by selecting the optimal individuals through adaptive grid. The results of the simulation experiments showed that the algorithm had better convergence speed and stability， and more diverse solution set distributions than other traditional algorithms in solving the firepower assignment problem.

Key words： firepower assignment; chaotic mapping; adaptive grid generation; multi-objective optimization; whale optimization algorithm

0 引言

隨著現(xiàn)代戰(zhàn)場(chǎng)智能化程度的不斷提高，對(duì)于火力籌劃方法的要求也越來越高，傳統(tǒng)的根據(jù)某項(xiàng)指標(biāo)對(duì)火力分配方案進(jìn)行篩選的方法已經(jīng)無法適應(yīng)現(xiàn)代戰(zhàn)場(chǎng)中復(fù)雜多變的環(huán)境，在如今仿真作戰(zhàn)場(chǎng)景下需要更加高效準(zhǔn)確的方法來處理大量戰(zhàn)場(chǎng)信息，對(duì)敵方目標(biāo)進(jìn)行火力分配方案的優(yōu)選是火力籌劃方法研究中的關(guān)鍵問題之一。

關(guān)于火力分配方案的選擇已經(jīng)不單單是基于彈藥對(duì)目標(biāo)最大毀傷效果的篩選，而是有選擇地根據(jù)現(xiàn)有資源合理分配，在對(duì)敵方造成盡可能大的毀傷效果的情況下使自身損耗盡可能?。?］。文獻(xiàn)［2］以作戰(zhàn)資源消耗最少、近距離作戰(zhàn)武器損失價(jià)值最小、敵方戰(zhàn)場(chǎng)剩余價(jià)值最小為優(yōu)化目標(biāo)構(gòu)建防御場(chǎng)景下的動(dòng)態(tài)火力分配模型，并利用A-NSGA-GKM算法對(duì)模型進(jìn)行求解，通過仿真實(shí)驗(yàn)證明了算法解決動(dòng)態(tài)火力分配問題的有效性。文獻(xiàn)［3］以實(shí)際毀傷與期望毀傷的差值平方和為目標(biāo)函數(shù)構(gòu)建火力分配模型，利用改進(jìn)遺傳算法對(duì)模型進(jìn)行求解，在達(dá)到預(yù)期毀傷要求的基礎(chǔ)上盡可能減少火力資源浪費(fèi)，有效增加了火力分配方案的靈活性。這些研究均根據(jù)不同的作戰(zhàn)場(chǎng)景設(shè)計(jì)了相適應(yīng)的優(yōu)化目標(biāo)。

火力分配方案的求解是一種多目標(biāo)優(yōu)化問題，如何在多個(gè)優(yōu)化目標(biāo)中尋求一種平衡是多目標(biāo)優(yōu)化要解決的主要問題，求解此類問題的本質(zhì)就在于找到問題模型的Pareto前沿解集［4］。文獻(xiàn)［5］提出一種基于協(xié)同進(jìn)化的混合變量多目標(biāo)粒子群優(yōu)化算法，利用協(xié)同進(jìn)化策略提高搜索效率，引入基于結(jié)構(gòu)學(xué)習(xí)的重組策略提高種群多樣性，解決了無人機(jī)協(xié)同多任務(wù)分配問題。文獻(xiàn)［6］提出了一種改進(jìn)NSGA-Ⅲ算法，利用連續(xù)編碼方式增加初始種群的多樣性，采用隨機(jī)選擇，模擬二進(jìn)制交叉和多項(xiàng)式變異作為種群進(jìn)化策略，有效解決了集群目標(biāo)火力分配問題。

上述方法雖然在算法收斂速度與效果上較一些經(jīng)典算法有所提升，但還是普遍存在難以跳出局部最優(yōu)解、全局搜索能力差的問題。鯨魚優(yōu)化算法（whale optimization algorithm， WOA）是一種元啟發(fā)式優(yōu)化算法，在2016年由Mirjalili等首次提出［7］，它通過模仿鯨魚的狩獵行為實(shí)現(xiàn)對(duì)目標(biāo)問題的求解，具有控制參數(shù)少、可擴(kuò)展性強(qiáng)、跳出局部最優(yōu)能力強(qiáng)的特點(diǎn)。在傳統(tǒng)的WOA基礎(chǔ)上，文獻(xiàn)［8］利用多種群探索機(jī)制增加種群多樣性、引入二次插值法提高算法收斂速度、通過控制參數(shù)設(shè)置平衡算法的開發(fā)與探索能力。文獻(xiàn)［9］使用柯西逆累積分布函數(shù)提高WOA的全局搜索能力，并利用改變權(quán)重的方法提高了算法的局部搜索能力和收斂精度。文獻(xiàn)［10］通過非線性種群初始化、增加邊界優(yōu)化條件和引入精英個(gè)體引導(dǎo)機(jī)制對(duì)WOA改進(jìn)，有效提高了算法的多樣性與收斂速度。

傳統(tǒng)的單目標(biāo)算法在求解復(fù)雜問題上存在較大的局限性，因此多目標(biāo)優(yōu)化算法的設(shè)計(jì)與研究具有重大的意義，文獻(xiàn)［11］提出了一種非支配排序的WOA，通過非支配排序擁擠度計(jì)算策略與狩獵行為對(duì)種群非支配解集的迭代更新實(shí)現(xiàn)多目標(biāo)優(yōu)化。文獻(xiàn)［12］利用反向精英保留策略和種群引導(dǎo)策略加快了算法收斂速度，引入種群引導(dǎo)策略和Levy變異策略，提高了算法的全局搜索能力。上述方法都利用計(jì)算個(gè)體擁擠度的方式在新組成的種群中挑選最優(yōu)個(gè)體，不但會(huì)增加算法的計(jì)算量，還會(huì)因?yàn)榉N群數(shù)量的限制而降低種群的多樣性。

本文針對(duì)火力分配提出一種改進(jìn)的多目標(biāo)鯨魚優(yōu)化算法，該方法利用混沌映射初始化種群，進(jìn)行自適應(yīng)網(wǎng)格劃分選擇最優(yōu)個(gè)體，進(jìn)而結(jié)合外部Pareto檔案進(jìn)化種群。仿真實(shí)驗(yàn)表明，本文方法較之同類方法具有更好的搜索能力與收斂性。

1 相關(guān)技術(shù)

多目標(biāo)優(yōu)化算法的重要環(huán)節(jié)就是尋找解空間中的非支配解集，本節(jié)對(duì)多目標(biāo)優(yōu)化、非支配排序以及鯨魚優(yōu)化算法進(jìn)行介紹。

1.1 多目標(biāo)優(yōu)化

最小化的多目標(biāo)優(yōu)化問題的數(shù)學(xué)模型表示為

minF（x）=（f1（x），f2（x），…，fm（x）），x∈S，（1）

其中：S為解空間；m為目標(biāo)函數(shù)的個(gè)數(shù)。多目標(biāo)優(yōu)化的目的就是找到x*=

argminx*∈S（f1（x），f2（x），…，fm（x））的解x*。在實(shí)際應(yīng)用中，這些目標(biāo)函數(shù)都是互相矛盾的，即優(yōu)化一個(gè)目標(biāo)函數(shù)的代價(jià)是劣化其他目標(biāo)函數(shù)，很難達(dá)到一種資源分配的理想狀態(tài)，并且這種理想狀態(tài)往往都存在多種解［13］，如何在解空間內(nèi)找到所有這樣的解集是多目標(biāo)優(yōu)化算法要解決的主要問題。

1.2 非支配排序

解決多目標(biāo)優(yōu)化問題首先要引入Pareto支配關(guān)系的概念，對(duì)于m個(gè)目標(biāo)函數(shù)fi（x），i=1，2，…，m，給定任意兩個(gè)決策變量x1和x2，如果滿足式（2）條件，則可稱解x1支配x2，

fi（x1）≤fi（x2）， i∈1，2，…，m，

fi（x1）<fi（x2）， i∈1，2，…，m。（2）

如圖1，設(shè)目標(biāo)函數(shù)f1和f2為越小越好，A、B、C、D、E和F為存在于目標(biāo)空間中的6個(gè)解。其中B在兩個(gè)目標(biāo)函數(shù)上都要優(yōu)于E，根據(jù)式（2）就可以稱B支配E，同理也可以得出C支配F，而對(duì)于A與D，由于解空間中不存在支配它們的解，因此A、B、C、D所在的超曲面就構(gòu)成了目前解空間中如圖1中曲線所示的Pareto前沿。

因此如果對(duì)于一個(gè)決策變量來說，沒有其他任何決策變量能夠支配它，那這個(gè)決策變量就稱為非支配解。非支配排序就是將所有決策變量中滿足此條件的解放入一個(gè)集合中，構(gòu)成Pareto前沿解集，標(biāo)記為第一層，之后排除找到的所有非支配解并對(duì)剩余個(gè)體再次進(jìn)行非支配排序，得到的Pareto前沿解集標(biāo)記為第二層，以此類推直至處理完種群中的所有個(gè)體［14］。

1.3 鯨魚優(yōu)化算法

鯨魚優(yōu)化算法模擬了自然界中鯨魚的捕食過程：鯨魚在發(fā)現(xiàn)獵物之后會(huì)成群地將獵物圍住，螺旋運(yùn)動(dòng)游向目標(biāo)并在這個(gè)過程中不斷吐出氣泡，形成一種柱狀“氣泡網(wǎng)”將獵物緊緊包圍在中心，直至鯨魚靠近一口吞下。鯨魚的捕食策略可以將其看作一個(gè)不斷逼近優(yōu)化問題最優(yōu)解的過程，其中每頭鯨魚可以看作是優(yōu)化問題中的一個(gè)解，獵物可以看作是需要求得的最優(yōu)解，因此可以模仿鯨魚的捕食過程來解決優(yōu)化問題［15］。

1） 環(huán)繞獵物

鯨魚是群居動(dòng)物，只要種群中某個(gè)個(gè)體發(fā)現(xiàn)獵物，則其他個(gè)體就會(huì)向它移動(dòng)爭(zhēng)搶獵物。在優(yōu)化問題中把具有最佳適應(yīng)度函數(shù)的個(gè)體視為種群的最優(yōu)個(gè)體，種群中其他個(gè)體朝著最優(yōu)個(gè)體變化位置的過程用數(shù)學(xué)模型表示為

X（t+1）=X*（t）-AD，

D=CX*（t）-X（t），（3）

其中：t為當(dāng)前的迭代次數(shù);X（t）為當(dāng)前鯨魚的位置;X*（t）為目前得到的最優(yōu)鯨魚位置;X（t+1）為更新后的鯨魚位置。A和C為相關(guān)系數(shù)，如式（4）所示，

A=2aR1-a，C=2R2，a=2（tmax-t）/tmax，（4）

其中：R1和R2為［0，1］上的隨機(jī)浮點(diǎn)數(shù)；tmax為設(shè)定的最大迭代次數(shù)；a為收斂因子，會(huì)隨著迭代次數(shù)的增加逐漸遞減至0。

2） “氣泡網(wǎng)”狩獵

鯨魚在狩獵的過程中會(huì)吐出“氣泡網(wǎng)”包圍獵物，因此設(shè)計(jì)出兩種數(shù)學(xué)模型來表示這種捕獵行為。

收縮包圍。通過公式（3）不斷收縮接近獵物，與公式（3）不同之處在于A的取值范圍，由于鯨魚在接近獵物的時(shí)候才會(huì)發(fā)出“氣泡網(wǎng)”困住獵物，而A隨著迭代次數(shù)的不斷增加會(huì)逐漸減小，所以這里設(shè)置在A<1時(shí)，鯨魚更新后的位置可以是它與最優(yōu)個(gè)體之間任意位置。

螺旋游動(dòng)。鯨魚會(huì)螺旋游向獵物，可以用數(shù)學(xué)模型表示為

X（t+1）=X*（t）+Dpeblcos（2πl(wèi)），

Dp=X*（t）-X（t），（5）

其中：Dp代表鯨魚與最優(yōu)個(gè)體之間的距離；b是一個(gè)用來定義螺旋形狀的常數(shù)，通常取1；l為（-1，1）上的隨機(jī)浮點(diǎn)數(shù)。

以上介紹的兩種狩獵行為發(fā)生在不同鯨魚的捕食過程中，為了更好地模擬行為的異步性，引入行動(dòng)概率P來決定鯨魚是選擇收縮包圍還是螺旋游動(dòng)接近獵物，鯨魚個(gè)體在生成的隨機(jī)數(shù)小于P時(shí)進(jìn)行收縮包圍，不小于P時(shí)進(jìn)行螺旋游動(dòng)，因此鯨魚“氣泡網(wǎng)”狩獵的數(shù)學(xué)模型為

X（t+1）=

X*（t）-ACX*（t）-X（t）， p<P，

X*（t）+X*（t）-X（t）eblcos（2πl(wèi)）， p≥P，（6）

其中：p為區(qū)間［0，1］上的隨機(jī)浮點(diǎn)數(shù)；P為行動(dòng)概率常數(shù)，通常設(shè)置為0.5。

3） 搜索狩獵

鯨魚還可以在全局搜索獵物，其數(shù)學(xué)模型表示為

X（t+1）=Xrand（t）-AD，

D=CXrand（t）-X（t），（7）

其中：Xrand（t）為在種群隨機(jī)選擇的鯨魚位置，當(dāng)A≥1時(shí)，鯨魚會(huì)隨機(jī)選擇種群中的個(gè)體進(jìn)行位置的更新，增強(qiáng)算法的全局搜索能力。

2 火力分配問題描述

假設(shè)有m種彈藥資源用于打擊n個(gè)地面目標(biāo)，彈藥總資源集合為W=｛w1，w2，…，wm｝，打擊目標(biāo)集合為T=｛t1，t2，…，tn｝，彈藥成本向量為C=｛c1，c2，…，cm｝，決策變量為彈藥數(shù)量集合X=｛x1，x2，…，xm｝。

針對(duì)問題描述選定火力分配的目標(biāo)函數(shù)為以下3個(gè)。

1） MTxT8Grxns+81l1zsWDPyPv/k93dpRKASKTtoe/ztgM=震塌比例。在目標(biāo)毀傷效果選取方面，目前大多數(shù)的相關(guān)研究均通過設(shè)定毀傷概率來模擬某彈藥打擊某目標(biāo)［16］，此類方法雖然能夠一定程度上反映火力分配過程中彈藥對(duì)目標(biāo)的毀傷效果，但是存在較大的主觀性與隨機(jī)性，并不能真實(shí)反映具體種類彈藥對(duì)于具體種類目標(biāo)的毀傷效果。因此本文以有限元仿真獲得的彈藥對(duì)特定目標(biāo)的毀傷效應(yīng)數(shù)據(jù)，在提高毀傷效果真實(shí)程度的同時(shí)，縮短仿真時(shí)間。

ye=Ansys（x1，x2，…，xm），（8）

其中：Ansys表示有限元仿真的過程；x1，x2，…，xm表示輸入的m種彈藥種類。

2） 彈藥成本。彈藥資源的消耗量也是火力分配方案優(yōu)劣性評(píng)估的重要指標(biāo)，表示為

yc=∑mi=1cixi。（9）

3） 自身剩余價(jià)值?？紤]在實(shí)際作戰(zhàn)中敵方作戰(zhàn)單元遭受火力打擊后的反擊行為和我方作戰(zhàn)單元作戰(zhàn)后對(duì)自身價(jià)值的消耗，這里將敵我雙方的對(duì)抗性行為和自身價(jià)值的消耗用自身剩余價(jià)值來描述，

yt=∑mi=1cixi（1-hiε（xi）），（10）

其中：hi為武器打擊目標(biāo)后自身的損傷概率；ε（xi）為階躍函數(shù)。

綜上所述，本文提出的火力分配問題可以描述為

max（ye），

min（yc），

max（yt），（11）

s.t. minBdi≤xi≤maxBdi，xi∈N，i=1，2，…，m，（12）

其中：式（11）為火力分配問題的優(yōu)化目標(biāo);式（12）為約束條件;maxBdi和minBdi分別為彈藥數(shù)量的上、下界。

3 改進(jìn)的多目標(biāo)鯨魚優(yōu)化算法

對(duì)所提出的改進(jìn)的自適應(yīng)網(wǎng)格多目標(biāo)鯨魚優(yōu)化算法（adaptive grid multi-objective WOA，AG-MOWOA）進(jìn)行介紹，包括種群初始化、最優(yōu)個(gè)體選擇、Pareto存檔進(jìn)化等過程。算法流程如圖2所示。

3.1 種群初始化

采用混沌映射的方法來進(jìn)行種群的初始化，Logistic方程是一個(gè)典型的混沌映射系統(tǒng)［17］，

Zn+1=μZn（1-Zn），n=0，1，2，…，（13）

其中：Z為［0，1］上的隨機(jī)浮點(diǎn)數(shù);μ被稱為L(zhǎng)ogistic控制參量，是［0，4］上的隨機(jī)浮點(diǎn)數(shù)，研究表明μ在（3.567，4］時(shí)整個(gè)混沌系統(tǒng)會(huì)呈現(xiàn)出一種偽隨機(jī)分布的狀態(tài)，且在靠近4時(shí)可以通過迭代不重復(fù)地經(jīng)歷空間中所有位置。所以用這種方法來初始化種群可以極大地提高初始種群在解的空間內(nèi)的分布質(zhì)量，具體操作步驟如下。

1） 初始化一個(gè)規(guī)模為m的種群，在［0，1］內(nèi)生成隨機(jī)解向量z0=（z0，1，z0，2，…，z0，m）。

2） 利用Logistic方程進(jìn)行迭代產(chǎn)生混沌序列。

3） 將得到的混沌序列映射到輸入解空間中，

Zi，j=zi，j（Cmax-Cmin）+Cmin，

Zi=（Zi，1，Zi，2，…，Zi，m），（14）

其中：Zi為混沌序列向量zi對(duì)應(yīng)的輸入解空間向量;Cmax和Cmin分別為輸入解空間的最大與最小值向量。

4） 假設(shè)種群大小為N，計(jì)算所有輸入解空間向量的目標(biāo)函數(shù)，并根據(jù)非支配關(guān)系選出其中最優(yōu)的N個(gè)個(gè)體作為初始種群。

3.2 最優(yōu)個(gè)體選擇

傳統(tǒng)的單目標(biāo)鯨魚優(yōu)化算法可以根據(jù)適應(yīng)函數(shù)值的大小選擇最優(yōu)個(gè)體，但是擁有多個(gè)目標(biāo)函數(shù)的種群往往有多個(gè)最優(yōu)個(gè)體，以什么樣的標(biāo)準(zhǔn)來選取領(lǐng)導(dǎo)鯨魚群更新位置是設(shè)計(jì)多目標(biāo)鯨魚優(yōu)化算法的重要問題之一。

目前大多數(shù)多目標(biāo)優(yōu)化算法都是通過計(jì)算個(gè)體間的擁擠度來實(shí)現(xiàn)非支配解集優(yōu)劣的排序，但是這個(gè)計(jì)算過程比較復(fù)雜且不利于算法的全局搜索，因此本文選用基于自適應(yīng)網(wǎng)格劃分［18］的方法來選擇最優(yōu)個(gè)體。

網(wǎng)格密度可以用來描述目標(biāo)空間中每個(gè)區(qū)域所包含的個(gè)體數(shù)量，區(qū)域內(nèi)的個(gè)體越多，網(wǎng)格密度就越大；個(gè)體越少，網(wǎng)格密度就越小。為了維持種群個(gè)體的多樣性，網(wǎng)格密度越小的區(qū)域內(nèi)的個(gè)體越有機(jī)會(huì)被選為最優(yōu)個(gè)體。以目標(biāo)函數(shù)一為例，關(guān)于網(wǎng)格信息計(jì)算方法的數(shù)學(xué)模型為

M=［Grid×yi（t）-minY（t）maxY（t）-minY（y）］，（15）

其中：M代表個(gè)體所處的網(wǎng)格編號(hào)；［·］表示括號(hào)內(nèi)數(shù)值向下取整;取小于括號(hào)內(nèi)數(shù)的最大整數(shù)；

Grid為網(wǎng)格要?jiǎng)澐值木W(wǎng)格數(shù)；t為迭代次數(shù)；yi（t）為第t代個(gè)體i的目標(biāo)函數(shù)值；maxY（t）為第t代種群所有個(gè)體中對(duì)應(yīng)目標(biāo)函數(shù)的最大值；minY（t）為第t代種群所有個(gè)體中對(duì)應(yīng)目標(biāo)函數(shù)的最小值。

此時(shí)通過式（15）得到所有個(gè)體的網(wǎng)格編號(hào)元組，利用輪盤賭法選擇最優(yōu)個(gè)體，處于同一個(gè)編號(hào)組的個(gè)體越多，該網(wǎng)格的密度就越大，選中這個(gè)網(wǎng)格中個(gè)體作為最優(yōu)個(gè)體的概率就越小。當(dāng)選中的網(wǎng)格有多個(gè)個(gè)體時(shí)，隨機(jī)選擇其中一個(gè)個(gè)體作為最優(yōu)個(gè)體。

3.3 Pareto存檔進(jìn)化

在目前許多群智能算法的研究中都會(huì)引入Pareto外部檔案用來存儲(chǔ)種群迭代過程中產(chǎn)生的非支配解，當(dāng)外部檔案中的個(gè)體超過提前設(shè)定好的值之后，就需要根據(jù)相應(yīng)的策略來刪掉多余的個(gè)體，保持外部檔案中個(gè)體的優(yōu)越性［19］。

按照這種思路，當(dāng)外部檔案中存放的非支配解超過設(shè)定值時(shí)，可以以不同個(gè)體的網(wǎng)格密度為評(píng)價(jià)標(biāo)準(zhǔn)對(duì)外部檔案進(jìn)行處理，將網(wǎng)格密度大的個(gè)體從檔案中刪除，直到外部檔案中的個(gè)體數(shù)量滿足要求。

4 仿真實(shí)驗(yàn)與結(jié)果分析

4.1 實(shí)驗(yàn)環(huán)境與參數(shù)設(shè)置

針對(duì)本文提出的地面目標(biāo)火力分配問題，設(shè)彈藥資源種類為W={w1，w2，w3，w4}，彈藥資源的取值范圍上、下區(qū)間分別為［0，0，0，0］和［100，100，50，100］，彈藥資源的成本向量為C=（20，50，30，15），hi為［0，0.3］區(qū)間內(nèi)的隨機(jī)值。仿真實(shí)驗(yàn)所使用的操作系統(tǒng)為Windows10，環(huán)境為Python3.7，PyCharm版本為2021.1.3 x64。

為了驗(yàn)證AG-MOWOA求解地面目標(biāo)火力分配方案的有效性，在實(shí)驗(yàn)?zāi)K分別采用NSGA-Ⅱ［20］、MOPSO［21］、NSWOA［11］以及AG-MOWOA四種算法對(duì)火力分配問題進(jìn)行對(duì)比分析，關(guān)于算法相關(guān)參數(shù)設(shè)置如下。

NSGA-Ⅱ：種群大小為100，迭代次數(shù)為50，存檔規(guī)模為100，交叉概率為0.8，變異概率為0.05。

MOPSO：種群大小為100，迭代次數(shù)為50，存檔規(guī)模為100，網(wǎng)格劃分?jǐn)?shù)為10。

NSWOA：種群大小為100，迭代次數(shù)為50，存檔規(guī)模為100。

AG-MOWOA：種群大小為100，迭代次數(shù)為50，存檔規(guī)模為100，網(wǎng)格劃分?jǐn)?shù)為10，Logistic控制參量為4。

4.2 實(shí)驗(yàn)結(jié)果分析

為了盡可能減少誤差，使實(shí)驗(yàn)結(jié)果更加準(zhǔn)確，這里獨(dú)立進(jìn)行三組實(shí)驗(yàn)，每組實(shí)驗(yàn)進(jìn)行10次，取這10次實(shí)驗(yàn)的平均值為實(shí)驗(yàn)結(jié)果，三種算法得到的最優(yōu)目標(biāo)函數(shù)值如表1所示。

從表1可以看出，在相同的迭代次數(shù)下，MOPSO、NSGA-Ⅱ和NSWOA得到的最優(yōu)解質(zhì)量不相上下，而AG-MOWOA所求得的三個(gè)最優(yōu)目標(biāo)函數(shù)值的質(zhì)量和穩(wěn)定性均優(yōu)于其他三種算法，說明該算法擁有更優(yōu)的全局搜索能力。

以第一組算法得到的解集分布為例，初始的解集分布如圖3所示，AG-MOWOA的分布狀況如圖4所示。從圖中可以看出，相較于初始解集， AG-MOWOA的分布更加集中，但是仍無法準(zhǔn)確判斷是否收斂。因此，需要用更加科學(xué)的方法來對(duì)比這四種算法的效果。

4.3 算法性能分析

為了測(cè)試各算法求解火力分配問題的性能，這里用三個(gè)指標(biāo)對(duì)求得的解集進(jìn)行評(píng)估，其中真實(shí)Pareto前沿取在范圍區(qū)間內(nèi)均勻分布的500條真實(shí)數(shù)據(jù)。

1） 反轉(zhuǎn)世代距離（IGD）

反轉(zhuǎn)世代距離計(jì)算的是每個(gè)參考點(diǎn)到最近解的歐氏距離［22］，如式（16）所示，可以用來評(píng)估算法的收斂性與多樣性，IGD的值越小，算法的綜合性能就越好，

IGD（P，P*）=∑x∈P*miny∈Pdis（x，y）P*，（16）

其中：P為算法求得的Pareto前沿解集；P*為真實(shí)Pareto前沿上均勻的參考點(diǎn)；P*為參考點(diǎn)的數(shù)量；dis（x，y）為點(diǎn)x與點(diǎn)y之間的歐氏距離。

通過計(jì)算得到的最小歐氏距離分布如圖5所示，IGD值大小如表2所示。

從圖5中可以看出AG-MOWOA求得的解集距離參考Pareto前沿的距離分布更加均勻，結(jié)合表2可以看出AG-MOWOA求得的IGD值更小，表明該算法獲得的解集更加均勻，并且在收斂性和多樣性上都占明顯優(yōu)勢(shì)。

2） C-metric解集覆蓋率

解集覆蓋率可以表示不同算法得到解集之間的互相支配關(guān)系，如式（17）所示，分子為滿足存在解集A中的點(diǎn)支配解集B中的任意條件的點(diǎn)的個(gè)數(shù)，分母為解集B中解的個(gè)數(shù)，如果解集覆蓋率C（A，B）大于C（B，A），則說明解集A的個(gè)體更接近真實(shí)Pareto前沿。

C（A，B）={x∈By∈A：yx}B，（17）

其中：A和B表示需要比較的兩個(gè)解集；yx表示y支配x；·表示滿足括號(hào)內(nèi)條件的個(gè)數(shù)。

記NSGA-Ⅱ?yàn)樗惴╝，MOPSO為算法b，NSWOA為算法c，AG-MOWOA為算法d，計(jì)算三組實(shí)驗(yàn)得到的各算法解集覆蓋率如表3所示。

從表3中可以看到，在三組實(shí)驗(yàn)中AG-MOWOA的解集覆蓋率遠(yuǎn)大于其他三種算法，表明該算法求得的解集相較于其他三種算法質(zhì)量更高，更接近真實(shí)Pareto前沿。

3） 算法消耗時(shí)間

每組10次實(shí)驗(yàn)的平均算法消耗時(shí)間統(tǒng)計(jì)得到的數(shù)據(jù)如表4所示。從表4中可以看出AG-MOWOA在消耗時(shí)間上也遠(yuǎn)遠(yuǎn)小于其他三種算法。

通過以上評(píng)估，可以看出本文提出的AG-MOWOA在收斂性、多樣性、均勻性、消耗時(shí)間上較其他算法均有提升，可以有效解決本文提出的火力分配問題。

5 結(jié)語

本文針對(duì)火力分配問題構(gòu)建了以震塌比例、彈藥成本和自身剩余價(jià)值為目標(biāo)函數(shù)的火力分配模型，并提出了一種自適應(yīng)網(wǎng)格多目標(biāo)鯨魚優(yōu)化算法，該方法采用混沌映射初始化種群，使初始種群在可行解空間內(nèi)分布更加均勻，利用自適應(yīng)網(wǎng)格劃分法為個(gè)體編號(hào)，不但可以挑選出最優(yōu)個(gè)體，還能減少算法計(jì)算量，最后引入Pareto存檔進(jìn)化策略有效增加了種群的多樣性。仿真實(shí)驗(yàn)結(jié)果表明，本文提出的算法在收斂性、多樣性、運(yùn)行速度等多方面性能較其他算法都有提升。針對(duì)多于三個(gè)目標(biāo)函數(shù)的火力分配問題的求解是本文所提算法今后需要改進(jìn)的方向。

參考文獻(xiàn)：

［1］ 劉志超， 石章松， 姜濤， 等. 基于最小資源損耗的火力分配研究［J］. 火力與指揮控制， 2018， 43（6）： 167-170.

LIU Z C， SHI Z S， JIANG T， et al. Research on method of weapon-target assignment based on the minimal resource loss［J］. Fire control & command control， 2018， 43（6）：167-170.

［2］ 于博文， 呂明. 基于改進(jìn)NSGA-III算法的動(dòng)態(tài)武器協(xié)同火力分配方法［J］. 火力與指揮控制， 2021， 46（8）：71-77， 82.

YU B W， LYU M. Method for dynamic weapon coordinative firepower distribution based on improved NSGA-Ⅲ algorithm［J］. Fire control & command control， 2021， 46（8）：71-77， 82.

［3］ 張明雙， 徐克虎. 基于最小火力浪費(fèi)的火力優(yōu)化分配［J］. 電光與控制， 2020， 27（9）： 55-59.

ZHANG M S， XU K H. Optimal firepower distribution based on minimum firepower waste［J］. Electronics optics & control， 2020， 27（9）： 55-59.

［4］ KONG L R， WANG J Z， ZHAO P. Solving the dynamic weapon target assignment problem by an improved multiobjective particle swarm optimization algorithm［J］. Applied sciences， 2021， 11（19）： 9254.

［5］ 王峰， 張衡， 韓孟臣， 等. 基于協(xié)同進(jìn)化的混合變量多目標(biāo)粒子群優(yōu)化算法求解無人機(jī)協(xié)同多任務(wù)分配問題［J］. 計(jì)算機(jī)學(xué)報(bào)， 2021， 44（10）： 1967-1983.

WANG F， ZHANG H， HAN M C， et al. Co-evolution based mixed-variable multi-objective particle swarm optimization for UAV cooperative multi-task allocation problem［J］. Chinese journal of computers， 2021， 44（10）： 1967-1983.

［6］ 聶俊峰， 陳行軍， 蘇琦. 基于NSGA-Ⅲ算法的集群目標(biāo)來襲火力分配建模與優(yōu)化［J］. 兵工學(xué)報(bào)， 2021， 42（8）：1771-1779.

NIE J F， CHEN X J， SU Q. Modeling and optimization of weapon-target assignment for group targets defense based on NSGA-Ⅲ algorithm［J］. Acta armamentarii， 2021， 42（8）：1771-1779.

［7］ MIRJALILI S， LEWIS A. The whale optimization algorithm［J］. Advances in engineering software， 2016， 95： 51-67.

［8］ SUN Y J， CHEN Y. Multi-population improved whale optimization algorithm for high dimensional optimization［J］. Applied soft computing， 2021， 112： 107854.

［9］ 郭振洲， 王平， 馬云峰， 等. 基于自適應(yīng)權(quán)重和柯西變異的鯨魚優(yōu)化算法［J］. 微電子學(xué)與計(jì)算機(jī)， 2017， 34（9）： 20-25.

GUO Z Z， WANG P， MA Y F， et al. Whale optimization algorithm based on adaptive weight and cauchy mutation［J］. Microelectronics & computer， 2017， 34（9）： 20-25.

［10］REN G， YANG R H， YANG R Y， et al. A parameter estimation method for fractional-order nonlinear systems based on improved whale optimization algorithm［J］. Modern physics letters B， 2019， 33（7）： 1950075.

［11］JANGIR P， JANGIR N. Non-dominated sorting whale optimization algorithm （NSWOA）： a multi-objective optimization algorithm for solving engineering design problems［J］. Global journals of research in engineering，2017： 17： 15-42.

［12］梁倩. 基于反向精英保留和Levy變異的多目標(biāo)鯨魚優(yōu)化算法［J］. 現(xiàn)代計(jì)算機(jī)， 2021（18）： 25-31.

LIANG Q. Multi-objective whale optimization algorithm based on reverse elite retention and Levy mutation［J］. Modern computer， 2021（18）： 25-31.

［13］SUN H， CAO A R， HU Z Y， et al. A novel quantile-guided dual prediction strategies for dynamic multi-objective optimization［J］. Information sciences， 2021， 579： 751-775.

［14］SONG Z H， LIU T， LIN Q Z. Multi-objective optimization of a solar hybrid CCHP system based on different operation modes［J］. Energy， 2020， 206： 118125.

［15］WANG Y K， WANG S L， LI D， et al. An improved multi-objective whale optimization algorithm for the hybrid flow shop scheduling problem considering device dynamic reconfiguration processes［J］. Expert systems with applications， 2021， 174： 114793.

［16］孫海文， 謝曉方， 孫濤， 等. 改進(jìn)型布谷鳥搜索算法的防空火力優(yōu)化分配模型求解［J］. 兵工學(xué)報(bào)， 2019， 40（1）： 189-197.

SUN H W， XIE X F， SUN T， et al. Improved cuckoo search algorithm for solving antiaircraft weapon-target optimal assignment model［J］. Acta armamentarii， 2019， 40（1）： 189-197.

［17］WANG L， LIU X Y， SUN M H， et al. A new chaotic starling particle swarm optimization algorithm for clustering problems［J］. Mathematical problems in engineering， 2018， 2018： 1-14.

［18］方冰， 李太勇， 吳江. 一種基于網(wǎng)格劃分的自適應(yīng)粒子群優(yōu)化算法［J］. 計(jì)算機(jī)應(yīng)用研究， 2010， 27（11）：4136-4139.

FANG B， LI T Y， WU J. Grid-based adaptive particle swarm optimization［J］. Application research of computers， 2010， 27（11）：4136-4139.

［19］KNOWLES J D， CORNE D W. Approximating the nondominated front using the Pareto archived evolution strategy［J］. Evolutionary computation， 2000， 8（2）： 149-172.

［20］DEB K， PRATAP A， AGARWAL S， et al. A fast and elitist multiobjective genetic algorithm： NSGA-II［J］. IEEE transactions on evolutionary computation， 2002， 6（2）： 182-197.

［21］COELLO C A， LECHUGA M S. MOPSO： a proposal for multiple objective particle swarm optimization［C］∥Proceedings of the 2002 Congress on Evolutionary Computation. Piscataway： IEEE Press， 2002： 1051-1056.

［22］ISHIBUCHI H， MASUDA H， TANIGAKI Y， et al. Difficulties in specifying reference points to calculate the inverted generational distance for many-objective optimization problems［C］∥IEEE Symposium on Computational Intelligence in Multi-criteria Decision-making. Piscataway： IEEE Press， 2015： 170-177.