王， 趙文飛， 滕克難， 周璐， 單鑫

(1.海軍航空大學 岸防兵學院， 山東 煙臺 264001； 2.92199部隊， 山東 青島 290014)

0 引言

海上要地是海軍兵力棲息、休整、待機、補給和維修的根據(jù)地，是海軍兵力生存和作戰(zhàn)的依托，戰(zhàn)時易成為敵空襲首要打擊的目標[1-2]。在海上要地防空作戰(zhàn)時，需要根據(jù)敵空襲兵器的出動規(guī)模以及我方防空武器裝備組成進行合理的火力分配，以實現(xiàn)作戰(zhàn)資源的最優(yōu)配置，對提高防空兵的制空效能有著重要意義[3]。

火力分配(WTA)是將來襲目標分配給可用武器系統(tǒng)，以取得最大目標射擊收益的過程，其中WTA問題又可分為靜態(tài)火力分配(SWTA)與動態(tài)火力分配(DWTA)兩類問題。傳統(tǒng)的SWTA與時間無關(guān)，只針對來襲目標一次攻擊進行一次火力分配。但在實際海上要地防空作戰(zhàn)過程中，原始的WTA方案往往不能有效應對各種突發(fā)狀況，因此需要根據(jù)實際攻防狀態(tài)對目標進行DWTA。

自Hosein等在1990年提出DWTA理念以來，引起越來越多的學者興趣[4]。Ahuja等[5]通過使用網(wǎng)絡流公式描述了WTA問題，并提出了下限解法。Kim等[6]開發(fā)了一個用于將可用武器重新調(diào)度到具有時間相關(guān)殺傷概率目標的模型。Karasakal等[7]提出了一種基于射擊- 觀察- 射擊(SLS)交戰(zhàn)政策和防御單位之間充分協(xié)調(diào)的防空分配模型。Davis等[8]為DWTA問題建模了馬爾可夫決策過程公式，并提出了近似動態(tài)規(guī)劃方法來解決該問題。Summers等[9]提出了一個DWTA 問題，作為來襲導彈分配攔截器的馬爾可夫決策過程。

許多國內(nèi)學者也就DWTA問題展開了深入研究，Wang等[10]提出了含有多種約束的通用型 DWTA 模型。文獻[11]根據(jù)對戰(zhàn)雙方的作戰(zhàn)態(tài)勢，提出了一種DWTA對抗博弈模型。Li等[12]為多階段WTA問題開發(fā)了一種雙目標公式，可最大限度地降低失敗風險和資源成本。雷鳴等[13]建立了以打擊效益最大為優(yōu)化目標的WTA模型，基于云模型的優(yōu)良特性，結(jié)合遺傳算法，提出了一種基于云模型的改進遺傳算法。高志華等[14]運用遺傳算法建立了WTA模型，使防空武器最大限度地發(fā)揮火力單位效能，達到最大毀傷效果。孫海文等[15]將射擊有利度、飛臨時間與威脅度等因素結(jié)合，構(gòu)建一種新的綜合防空WTA模型。

以上國內(nèi)外學者的研究大多數(shù)認為按最初的WTA方案實施作戰(zhàn)，便可獲得最大的來襲目標攔截效率。然而，在實際海上要地防空作戰(zhàn)過程中，原始WTA計劃可能會受到原分配目標被消滅、出現(xiàn)新目標以及某一防空武器發(fā)生故障等不確定因素的影響，而使得初始方案不可行或攔截效率低下。

本文針對海上要地防空火力再分配問題，提出了一種能夠兼顧最大化攔截概率以及最大化武器射擊穩(wěn)定性的雙目標優(yōu)化DWTA模型。其中定義防空武器的射擊穩(wěn)定性為某臺防空武器對空中目標持續(xù)射擊的輪數(shù)。通過將目標函數(shù)線性化的方法進行雙目標函數(shù)求解，對3種不確定突發(fā)情況進行算例仿真求解，驗證了本文提出方法的可行性與有效性。通過與單目標火力規(guī)劃模型的對比驗證了本文算法的優(yōu)越性，通過大規(guī)模樣本的算例仿真驗證了本文算法的穩(wěn)定性。

1 不確定因素下防空火力優(yōu)化模型

1.1 問題描述

假設采用n個武器攔截m個目標，武器和目標之間的交戰(zhàn)組由集合V表示，其中(i,j)∈V表示武器系統(tǒng)j與目標i交戰(zhàn)，i∈M,j∈N，其中M為來襲目標集合，N為防空武器集合。在確定有效的交戰(zhàn)組合時，要綜合考慮武器系統(tǒng)的射程、交戰(zhàn)區(qū)域、攔截概率和火力通道等因素[16-19]。同時假設在開始作戰(zhàn)時按照初始WTA方案進行，但在交戰(zhàn)過程中原方案可能會被以下3種因素影響。

1)原分配方案中的目標已被摧毀。如果一個武器裝備在前幾輪的射擊中摧毀了分配給它的目標，那么原分配方案中之后用來攻擊被摧毀目標的剩余導彈可以用來防御其他目標。

2)武器裝備故障。防空武器裝備由諸多電子和機械部件構(gòu)成，這些組件中的一個或多個由于戰(zhàn)損、電子干擾、自然環(huán)境因素等發(fā)生故障，則可能會導致防空武器在交戰(zhàn)期間無法使用。

3)交戰(zhàn)過程開始后發(fā)現(xiàn)新的來襲目標。交戰(zhàn)過程中的另一個不確定因素是出現(xiàn)了原WTA預案中沒有考慮到的新目標。

1.2 限制條件

當上述突發(fā)情況發(fā)生時，需求重新進行WTA問題，將重新進行WTA的時間間隙設置為重調(diào)度時間點。當前目標和可用武器裝備對的集合會在重新調(diào)度的時間點更新。

針對本文的研究目的提出以下模型假設：

1)在一個時間點，最多發(fā)生一個擾動。

2)如果某武器在交戰(zhàn)過程中發(fā)生故障，則在交戰(zhàn)過程結(jié)束之前無法使用。

3)當發(fā)現(xiàn)到新的空中目標時，能夠檢測出新來襲目標的類型、距離、速度和命中概率。

4)各防空武器之間在WTA方面能夠充分協(xié)調(diào)。

1.3 相關(guān)定義

將整個防空作戰(zhàn)時間范圍劃分為相等、不重疊的時間間隙，每次攔截作戰(zhàn)從這些時間間隙的起始開始，用k∈K表示，其中k表示每次攔截的時間間隙，K表示整個防空作戰(zhàn)時間間隙集[20-22]。假設每次攔截需要恒定的準備時間Δc。防空武器的最大射程和最小射程分別用rmax j和rmin j表示。一個防空武器可以在交戰(zhàn)時段[tfij,tlij]內(nèi)攔截目標，其中tfij是第一次交戰(zhàn)的開始時間，tlij是最后一次交戰(zhàn)的結(jié)束時間。武器j與空中目標i的可交戰(zhàn)時間是由防空裝備的有效射程和速度以及來襲目標的距離和速度計算得到的。

tfij的計算公式如(1)式所示:

(1)

tlij根據(jù)武器最小射程和攔截發(fā)生時目標與防御點距離計算：

(2)

作戰(zhàn)時間為射擊準備時間與敵我目標相遇時間的和，用Δijk表示。飛行時間取決于武器- 目標對的相對速度和交戰(zhàn)的開始時隙ρk。如果交戰(zhàn)在時間間隙k開始，則交戰(zhàn)持續(xù)時間Δijk的公式為

(3)

圖1顯示了時隙、時間范圍和時隙的開始時間。來為每個武器- 目標對計算可交戰(zhàn)時間段。

圖1 防空作戰(zhàn)時間線

根據(jù)SLS防空作戰(zhàn)原則，對每個空中目標進行武器分配。根據(jù)SLS原則，在當前交戰(zhàn)完成之前，不會對空中目標進行新的WTA。定義集合Tij為武器j和目標i的可調(diào)度交戰(zhàn)時間段，因為武器j和目標i對之間作戰(zhàn)[tfij,tlij]內(nèi)實現(xiàn)。

Tij={k∈K:(i,j)∈Vand [ρk,ρk+Δijk]?[tfij,tlij]}

(4)

武器j和目標i之間的最大交戰(zhàn)次數(shù)由fij表示，并根據(jù)SLS原則通過除以可交戰(zhàn)時段計算得出。為了確保根據(jù)SLS原則安排交戰(zhàn)，為每個目標i和時隙k定義特定集，

Lik={(j,τ):(i,j)∈V,τ∈Tijand [ρk,ρk+Δ]?[ρτ,ρτ+Δijτ]}

(5)

如(5)式所示，若目標i和武器j之間的交戰(zhàn)在時隙k之前的時隙開始，并在時隙k結(jié)束時間之后結(jié)束，則(j,τ)對在集合Lik中。

當交戰(zhàn)在時隙k的開頭開始時，武器j對目標i的單發(fā)殺傷概率用pijk表示，并且每個防空武器在重新分配時間點都有可用的次數(shù)uj。

構(gòu)建下列函數(shù)以描述整個防空作戰(zhàn)過程：

(6)

(7)

定義射擊穩(wěn)定性為連續(xù)分配給同一武器的目標總數(shù)。圖2描述了一個防空武器兩個不同射擊方案的示例，箭頭顯示了交戰(zhàn)范圍內(nèi)分配的來襲目標編號。圖2(a)中，防空武器在交戰(zhàn)過程中對同一個空中目標i2進行射擊，并且由于不改變目標而具有更好的射擊穩(wěn)定性。圖2(b)中，武器必須在每次射擊后改變其目標，首先對空中目標2射擊，然后對空中目標1射擊，最后對空中目標3射擊，因此具有較低的射擊穩(wěn)定性。

圖2 防空武器的兩種不同目標分配方案

1.4 WTA優(yōu)化目標函數(shù)

DWTA方案優(yōu)化目標函數(shù)如下：

(8)

(9)

式中：Rj為每型武器最大射擊輪數(shù)。

(10)

(11)

式中：fij為目標交戰(zhàn)次數(shù)。

(12)

目標函數(shù)即攔截概率最大化SC即最大化來襲目標攔截概率；穩(wěn)定性目標函數(shù)WD即最大限度地將相同的空中目標分配到一個武器的作戰(zhàn)安排中。防空武器射擊次數(shù)由(9)式約束，對于每個武器系統(tǒng)j，分配給空中目標i的輪數(shù)必須小于或等于可用輪數(shù)。

約束條件(10)式表示武器射擊方案必須符合SLS原則。根據(jù)集合(j,τ)∈Lik對所有Xijτ變量求和，可確保從不同武器到空中目標的射擊不會發(fā)生時間重疊。約束條件(11)式明確了可以針對每個目標安排的最大武器輪數(shù)。假設根據(jù)SLS原則，在可交戰(zhàn)時限內(nèi)，武器系統(tǒng)j可以分配給空中目標i的輪數(shù)為2，則約束條件(12)式強制規(guī)定，武器系統(tǒng)j最多可以調(diào)度2輪對空目標i進行射擊。

2 考慮射擊穩(wěn)定性的WTA步驟

本節(jié)提出一種考慮防空武器射擊目標交換的WTA優(yōu)化算法。通過在每一步中定義攔截概率的上界來生成最大攔截率與最小攔截率之間的可行解，最后求得各攔截率對應的最佳射擊穩(wěn)定性的目標函數(shù)值[23-25]。

2.1 效率目標函數(shù)的線性化

(13)

(14)

(15)

(16)

(17)

0≤γi1≤C1,?i∈M

(18)

0≤γip≤Cp-Cp-1,?i∈M,p=2,3,…,l

(19)

和(9)式～(12)式。

2.2 生成具有不同攔截概率的所有解決方案

為解決雙目標WTA問題，首先計算在不考慮穩(wěn)定性目標函數(shù)的情況下生成具有不同攔截效率值的所有解。為了獲得具有不同攔截概率的方案，通過在每個步驟中向模型添加以下約束來限制SCl，直到生成具有最小效率值的攔截方案。

(20)

在約束條件(20)式中，將目標函數(shù)值限制為具有不同攔截概率的目標函數(shù)解決方案，SCmax l值是攔截概率的上限。在第一步中SCmax l設置為1，這是攔截概率可以取的最大值。然后，在每個步驟中通過等于先前找到的攔截概率值來更新SCmax l。通過在每一步中更新在約束條件(20)式的右側(cè)，得到一個新的解決方案，該解決方案的攔截概率比之前找到的要小。繼續(xù)迭代，直到生成具有最小攔截概率的WTA方案，本階段不考慮武器系統(tǒng)的射擊順序。

2.3 確定不同攔截率WTA方案穩(wěn)定性值

本節(jié)關(guān)注生成的WTA方案的射擊目標更改總數(shù)，提出一種武器- 目標交換算法，以找到具有不同攔截概率的每個解決方案的最大WD值。交換算法主要包括兩種交換情況。

1)每個目標交換分配的武器組合。在目標的WTA方案表中選擇一對武器組合，并嘗試交換其順序。搜索是否存在可用的時間段來分配兩種武器。如果可以進行交換，將保留新射擊方案作為備選WTA方案，然后為每個目標生成所有交換武器的組合。例如，目標的初始WTA方案如圖3(a)所示，y軸代表時間范圍，矩形表示防空武器。根據(jù)SLS原則，在防空武器攔截目標期間，在同一時間段內(nèi)不得為目標安排其他武器，因此矩形之間沒有重疊。矩形的起點和終點顯示了交戰(zhàn)的起點和終點時間段。矩形的長度代表每次作戰(zhàn)的持續(xù)時間。在交換算法的第一階段，交換武器2和武器1的射擊命令，并為目標1生成新的時間表。圖3(b)顯示了更新后的時間表，新方案對于目標1具有相同的攔截概率。

圖3 武器- 目標分配方案

2)在兩個目標的WTA方案間交換防空武器，保持相同的攔截概率，以在每個可行解決方案中找到最大WD。如果兩個不同的武器對兩個目標具有相同的單發(fā)殺傷概率值，則這些防空武器設計次序可以在兩個目標的WTA方案之間交換，并且整個防空的攔截概率值不會改變。對于每對目標，檢查不同武器是否具有相同的單發(fā)擊殺概率值。如果可以在目標之間交換分配的武器，則將更新當前WTA方案。

交換算法通過更新武器系統(tǒng)的射擊順序，來找到WD值最大且不改變攔截概率的WTA方案，最后一個步驟是去除之前產(chǎn)生的不可行WTA方案。在尋找最大WD值的過程中可能會出現(xiàn)一個武器在一個時間段內(nèi)射擊兩次的情況。若出現(xiàn)此類情況，則必須將其從可行的分配方案集中刪除。計算可行解集中所有解的WD值，并選擇WD值最大的解作為備選方案。

對于每對目標，需要檢查不同防空武器是否存在相同的單發(fā)擊殺概率值。若相同，則交換對原目標射擊的武器，找到所有可能的武器分配方案。例如，假設目標1-武器1、目標1-武器2、目標2-武器1和目標2-武器2之間的單發(fā)殺傷概率相等。圖4(a)顯示了對目標1的初始WTA方案。若交換第二段武器1和武器2射擊目標，則整體的攔截效率不會發(fā)生變化，因為單發(fā)命中概率值相同。另一方面，武器1和武器2的交戰(zhàn)時間表的穩(wěn)定性可以增加。這可以在相同的效率下產(chǎn)生最大的穩(wěn)定性值。圖4(b)顯示了更新后的分配計劃。

圖4 單發(fā)攔截概率相同時可交換對兩個目標射擊的防空武器

圖5表示一個簡單的求WTA方案非支配解的求解過程。首先，在不考慮射擊穩(wěn)定性的條件下，通過求解線性化模型得到具有最大效率目標函數(shù)值的解AS1。其次，通過武器- 目標交換算法在保持相同攔截效率值的同時最大化射擊穩(wěn)定性，并生成解AS2。最后，使用攔截效率值的界限對模型進行求解，以獲得具有不同攔截效率的解，生成解AS3。為了生成此特定攔截效率值的最佳穩(wěn)定性值，再次應用武器- 目標交換算法并生成解AS4。結(jié)果生成了兩個解AS2和AS4。重復進行上述步驟，直到獲得最低攔截效率的解。

圖5 一個WTA方案的生成過程

具體算法步驟如下：

步驟1設置SC=1和WDcur=0。

步驟2進行雙目標武器再分配，尋找滿足SC、WD的武器分配方案。如果不可行，則停止。

步驟3可行設置解cs并將SCcur設置為攔截概率值，計算WD。

步驟4對于當前的方案cs，Ui更改集合Ui中為每個目標i分配的武器順序。找到所有可能的武器序列，并為每個目標i添加設置XLi，XLi為目標i分配的所有可能的武器序列集。

步驟5從集合Ui中為每個目標 選擇一個分配計劃，并找到防御的交戰(zhàn)分配計劃的所有組合并添加到集合JZ(防御交戰(zhàn)分配方案集)。

步驟6判斷集合JZ中的解是否包括武器系統(tǒng)在同一時間段內(nèi)多次射擊，從集合JZ中刪除這些不可行的解決方案。

步驟7對于每個目標對，檢查是否有機會用相同的單發(fā)命中概率對交換分配的武器。如果可以進行交換，則更新當前計劃cs。轉(zhuǎn)步驟2。

步驟8為集合Ui中的每個解決方案計算WD。

步驟9從集合Ui中選擇具有最大WD的解λ。

步驟10如果WDλ則將解WDcur添加到 并設置WDλ=WDcur。將約束SCλ≤SCcur-ε添加到雙目標火力分配方案，轉(zhuǎn)步驟1。

算法實現(xiàn)邏輯圖如圖6所示。

圖6 算法步驟示意圖

3 算例仿真

表1為防空武器和來襲目標性能參數(shù)。根據(jù)海上要地防空作戰(zhàn)的實際情況，從表1的三型防空武器裝備中選取防空武器，對四型空中目標進行防空作戰(zhàn)仿真，驗證本文提出的WTA算法。各武器裝備的參數(shù)，參照各種文獻資料中已經(jīng)公開的參數(shù)進行計算。從每型武器系統(tǒng)中分別選出2臺和3臺武器系統(tǒng)，各自構(gòu)成2×2×2與3×3×3的防空線。仿真環(huán)境為CPU:i7-8850H，16.0GB RAM，操作系統(tǒng)Windows10，仿真實驗工具MATLABR2018a。

表1 防空武器系統(tǒng)和空襲目標參數(shù)

從表1中隨機選擇一個空中目標，并隨機選擇這個空中目標要摧毀的防御單位。重復此步驟，直到選擇夠m個空中目標為止。據(jù)此建立武器與目標(i,j)∈V的有效交戰(zhàn)組合矩陣。隨機設定當前空中目標距離為60～80 km，每個情況中的武器- 目標對都是在表中隨機選擇的。

規(guī)定時隙的持續(xù)時間為1 s，從[0.50，0.80]概率區(qū)間內(nèi)的均勻分布為目標和防空武器創(chuàng)建了一個單次射擊命中概率矩陣樣本。每次射擊的準備時間為9 s，每次可行攔截所花費的時間由恒定的準備時間和目標相遇時間之和計算。

在進行算例仿真時，首先根據(jù)最大攔截概率建立最優(yōu)初始調(diào)度進行防御；然后通過仿真程序模擬1.1節(jié)中考慮的以下3種突發(fā)情況：

1)根據(jù)均勻分布生成一個隨機數(shù)，當隨機數(shù)值小于單發(fā)命中概率時，則武器- 目標對的交戰(zhàn)結(jié)果是擊毀目標。

2)針對武器系統(tǒng)故障的情況，假定其中一臺武器裝備在交戰(zhàn)過程開始后無法發(fā)射，武器的破損和破損時間隨機確定。

3)當有新的來襲目標時，假設雷達發(fā)現(xiàn)一個突然的來襲目標后，交戰(zhàn)過程開始，初始計劃生效。泊松分布隨機確定到達點，新目標的距離和速度是隨機確定的，得到新目標與防空武器之間的有效組合、交戰(zhàn)的可用時隙和最大交戰(zhàn)次數(shù)。

本文定義的相關(guān)屬性參數(shù)符號及含義如表2所示。來襲目標和防空武器從表1中隨機選取，仿真計算得到不同情況下較優(yōu)非支配解的仿真結(jié)果如表3所示。

表2 部分本文定義的屬性符號及含義

表3中，當武器系統(tǒng)數(shù)量為9、來襲目標數(shù)量為3時多次算例仿真的結(jié)果一致，即SC=1,WD=9，因此在表2中只對這一組結(jié)果進行記錄。為更好地評估WTA方案的攔截效率和穩(wěn)定性，本文定義了不同的度量標準，并對結(jié)果進行了分析。具體度量的量及含義如表4所示。計算結(jié)果見表5。由表5可見：當目標數(shù)為6、武器系統(tǒng)數(shù)為6時，平均MIE和APIE值最高；當目標數(shù)為3、武器系統(tǒng)數(shù)為9時，平均MIE和APIE值最低。

表3 各情況下的較優(yōu)非支配解

表4 本文定義的度量的量及含義

表5 不同情況下數(shù)據(jù)處理結(jié)果

在目標被提前消滅的情況下，MIE均值在0～0.285之間。在6個目標-6個武器系統(tǒng)和3個目標-9個武器系統(tǒng)的情況中，MIE的最大值分別為0.282和0.076。結(jié)果表明：目標數(shù)量增加或武器系統(tǒng)數(shù)量減少，平均MIE和APIE結(jié)果增加；在6個目標-6個武器系統(tǒng)，MTC的最大平均值為3.667。在3個目標-9個武器系統(tǒng)的情況中MTC的平均值最低，為0。結(jié)果表明，平均MTC隨著目標數(shù)量的增加而增加，隨著武器系統(tǒng)數(shù)量的增加而減少。在6個目標-6個武器系統(tǒng)的情況中，獲得了最大APIEC值，達到3.367%，因此可以說一次目標轉(zhuǎn)換對該情況下的攔截效率提高了3.367%。當3個目標-9個武器系統(tǒng)時不需要進行射擊次序交換即可完成攔截任務，因此APIEC值無實際意義。

在武器故障情況下，當6個武器系統(tǒng)的情況中目標數(shù)從3個增加到6個時，MIE平均值從0.089增加到0.227。此外，在這些問題集中，MPIE平均值從3.947%增加到47.301%。在6個目標-6個武器系統(tǒng)的情況下，獲得了最高的MIE、MPIE和APIEC平均值。在3個目標-6個武器系統(tǒng)的情況中，獲得了最低的平均MIE、MPIE和APIEC值。平均MTC值隨著武器系統(tǒng)和目標數(shù)量的增加而增加。在6個目標-6個武器系統(tǒng)的情況中，最大平均MTC值分別為1.333。在6個目標-6個武器系統(tǒng)的情況中，最大MPIE值達到49.705%。因此，在該情況下幾乎實現(xiàn)了50%的攔截效率改進增量。

在新目標情況下，在6個目標-6個武器系統(tǒng)的情況中，最高的平均MIE值為0.261；在6個目標-9個武器系統(tǒng)的情況下，它減少到0.228。在6個目標-6個武器系統(tǒng)的情況中，MPIE的平均最高值為40.493%。因此，系統(tǒng)的攔截效率在這種情況下提高了50%。在3個目標-9個武器系統(tǒng)的情況下，MPIE平均值最低的是0.821%。

結(jié)果表明，在各種突發(fā)情況下，當武器系統(tǒng)數(shù)量有限和來襲目標多的情況下，重調(diào)度方法的效果更為顯著。當武器系統(tǒng)數(shù)量增加或空中目標數(shù)量減少時，效率的提高就會降低。MTC值隨著武器系統(tǒng)或目標數(shù)量的增加而增加。

表6顯示了各情況下的仿真運行時間。運行時間取決于武器與目標的有效交戰(zhàn)次數(shù)、各武器系統(tǒng)中武器彈藥的數(shù)量。當問題變得復雜時，運行時間就會增加，這與現(xiàn)實情況基本一致。由于海上要地防空作戰(zhàn)WTA屬于在戰(zhàn)前進行，算法仿真運行時間基本可以滿足要求。

為了驗證本文提出算法的優(yōu)越性，采用盡考慮來襲目標攔截概率的單目標優(yōu)化算法，在相同仿真環(huán)境、相同仿真條件下進行DWAT仿真計算，具體仿真結(jié)果如表7所示。

表6 不同情況下求解的耗時

表7 僅考慮來襲目標攔截概率的單目標優(yōu)化算法仿真結(jié)果

由表6可知，僅考慮來襲目標攔截概率的單目標DWTA算法較本文提出的算法目標攔截概率較高。但單目標算法給出的WTA方案武器射擊穩(wěn)定性過低，在實際防空作戰(zhàn)時由于不同武器系統(tǒng)之間的反復調(diào)度很可能會延誤寶貴的作戰(zhàn)時機、浪費有限的作戰(zhàn)資源，不能很好地滿足海上要地防空作戰(zhàn)對WTA方案時效性的要求。

為測試本文方法的穩(wěn)定性，在相同仿真環(huán)境下，擴大實驗規(guī)模分別采用來襲目標數(shù)9-防空武器數(shù)6；來襲目標數(shù)12-防空武器數(shù)9為初始條件，驗證本文提出的WTA算法。具體仿真結(jié)果見表8。由表8可見，本文提出的模型在來襲目標數(shù)量較大的情況下，依然能夠收斂，給出目標攔截概率及武器系統(tǒng)射擊穩(wěn)定性。

表8 擴大實驗規(guī)模后的仿真結(jié)果

4 結(jié)論

本文提出了一種考慮不確定因素的新的防空武器DWTA模型。將原分配目標被擊毀、防空武器故障和發(fā)現(xiàn)新目標3種情況作為當前分配方案的不確定擾動。建立了一種雙目標模型，將攔截概率作為效率測度，將對同一目標連續(xù)射擊的總次數(shù)作為穩(wěn)定性測度。同時為求解雙目標數(shù)學模型，提出了一種新的求解方法。由于WTA問題是一個非線性整數(shù)規(guī)劃問題，首先對效率目標函數(shù)采用了線性化方法，然后在不考慮防空武器射擊穩(wěn)定性的條件下，生成了具有不同攔截概率的可行解；最后利用這些可行解，在防空武器分配方案中通過交換射擊目標次序生成非支配解。結(jié)果表明，本文提出的方法可用于解決不確定因素下的動態(tài)火力規(guī)劃問題，并能輔助指揮員進行防御決策制定。