周榮坤，張永利，計(jì)文平

（中國(guó)電子科技集團(tuán)公司電子科學(xué)研究院，北京100041）

0 引 言

對(duì)空目標(biāo)威脅估計(jì)是艦艇編隊(duì)防空作戰(zhàn)過(guò)程進(jìn)行火力分配的重要理論依據(jù)。現(xiàn)代海戰(zhàn)中，空襲的主要模式是多層次、多批次、多類型的飽和攻擊，所以威脅估計(jì)需要針對(duì)不同類型的來(lái)襲目標(biāo)，并把來(lái)襲目標(biāo)各評(píng)估要素關(guān)聯(lián)起來(lái)。本文建立了適于防空系統(tǒng)的威脅評(píng)估模型，得到目標(biāo)的威脅度并實(shí)現(xiàn)評(píng)估結(jié)果排序，同時(shí)根據(jù)評(píng)估結(jié)果，結(jié)合艦空導(dǎo)彈毀傷概率對(duì)目標(biāo)實(shí)行火力分配。

1 威脅估計(jì)

防空作戰(zhàn)中，及時(shí)準(zhǔn)確地判斷來(lái)襲目標(biāo)的威脅程度，是合理進(jìn)行防空火力分配的重要前提。關(guān)于目標(biāo)威脅評(píng)估的研究有很多［1－6］。通常情況下，根據(jù)來(lái)襲目標(biāo)的屬性，建立威脅評(píng)估模型可以分為4個(gè)層次：

（1）確定威脅評(píng)估要素；

（2）建立各評(píng)估要素的數(shù)學(xué)表達(dá)式；

（3）綜合各因素求取目標(biāo)威脅度；

（4）確定來(lái)襲目標(biāo)威脅等級(jí)。

1.1 空中目標(biāo)威脅評(píng)估影響因素

選擇來(lái)襲目標(biāo)威脅評(píng)估因素不僅要反映目標(biāo)的威脅程度，而且要便于量化。威脅評(píng)估影響因素既包含定性因子，也包含定量因子。

本文將來(lái)襲目標(biāo)類型、高度、距離、速度作為威脅評(píng)估的主要影響因素，其中目標(biāo)類型為定性指標(biāo)，具有模糊集特征，由專家根據(jù)作戰(zhàn)經(jīng)驗(yàn)等實(shí)際情況用語(yǔ)言變量描述給出；高度、距離、速度為定量指標(biāo)，由傳感器測(cè)得參數(shù)并進(jìn)行數(shù)據(jù)處理得到。

1.2 目標(biāo)威脅評(píng)估模型

1.2.1 定性指標(biāo)的Vague值表示方法

目標(biāo)類型對(duì)于目標(biāo)威脅度的影響，很難通過(guò)具體公式進(jìn)行處理。本文通過(guò)Vague值11級(jí)語(yǔ)言變量表示，將指標(biāo)變量轉(zhuǎn)換成Vague值，取值在［0，1］之間，如表1所示。

表1 Vague值表示的11級(jí)語(yǔ)言變量

1.2.2 定量指標(biāo)的數(shù)學(xué)表達(dá)式

1.2.2.1 目標(biāo)與艦艇編隊(duì)的距離

目標(biāo)距離是指來(lái)襲目標(biāo)到達(dá)艦艇編隊(duì)的距離。按距離因素計(jì)算目標(biāo)威脅度的數(shù)學(xué)表達(dá)式為：

式中：r1，r2為目標(biāo)威脅的最近距離和最遠(yuǎn)距離。

當(dāng)目標(biāo)與艦艇編隊(duì)的距離逼近到一定范圍內(nèi)，即r≤r1時(shí)，來(lái)襲目標(biāo)對(duì)艦艇編隊(duì)造成毀滅性的打擊，即其威脅計(jì)算結(jié)果為最大值1；當(dāng)來(lái)襲目標(biāo)與艦艇編隊(duì)遠(yuǎn)離到一定程度，即r≥r2時(shí)，對(duì)艦艇編隊(duì)基本不造成威脅，可認(rèn)為其威脅計(jì)算結(jié)果為最小值0。

1.2.2.2 目標(biāo)距艦艇編隊(duì)的高度

來(lái)襲目標(biāo)的高度是指目標(biāo)所在位置到水平航路的距離。威脅目標(biāo)離艦艇編隊(duì)的高度越低，攻擊意圖越明顯，對(duì)艦艇編隊(duì)的威脅也越大；目標(biāo)高度越高，威脅程度越小。按高度因素計(jì)算目標(biāo)威脅度的數(shù)學(xué)表達(dá)式為：

式中：h1，h2分別為目標(biāo)威脅的最低高度和最高高度。

1.2.2.3 目標(biāo)的速度

當(dāng)目標(biāo)速度達(dá)到其最大值的時(shí)候，認(rèn)為該目標(biāo)對(duì)艦艇編隊(duì)的威脅達(dá)到最大值1，當(dāng)目標(biāo)以最小速度飛行時(shí)，認(rèn)為其對(duì)艦艇編隊(duì)的威脅為較弱的值，可定義為0.3。目標(biāo)威脅度與目標(biāo)速度的數(shù)學(xué)表達(dá)式為：

式中：vmax，vmin為目標(biāo)飛行速度的最大值和最小值。

1.3 艦艇編隊(duì)目標(biāo)威脅估計(jì)

1.3.1 數(shù)據(jù)源描述數(shù)據(jù)源描述如表2所示。

表2 數(shù)據(jù)源描述

根據(jù)Vague值表示以及公式（1）～（3），得到目 標(biāo)定性和定量指標(biāo)威脅屬性值μij，如表3所示。

表3 目標(biāo)影響因素的威脅屬性值

1.3.2 各威脅評(píng)估因素的權(quán)重設(shè)置

對(duì)來(lái)襲目標(biāo)進(jìn)行威脅評(píng)估，不同的評(píng)估影響因子具有不同的權(quán)重系數(shù)。文中利用層次分析法確定威脅目標(biāo)影響因子的權(quán)重值。

計(jì)算判斷矩陣每行所有元素的幾何平均值：

于是得到W ＝ （ω1，ω2，…，ωi）T。

計(jì)算判斷矩陣的最大特征根：

為了避免由于估計(jì)誤差引起特征值和特征向量產(chǎn)生偏差，導(dǎo)致判斷不一致，需要求出一個(gè)判斷矩陣一致性的指標(biāo)IC，即：

隨著n的增加，判斷誤差就會(huì)增加，因此應(yīng)當(dāng)考慮到n的影響，使用隨機(jī)性一致性比值：

式中：IR為平均隨機(jī)一致性指標(biāo)，如表4所示。

當(dāng)RC＜0.1時(shí)，判斷矩陣的一致性是可以接受的。

表4 平均隨機(jī)一致性指標(biāo)

采用層次分析法確定來(lái)襲目標(biāo)影響因子的權(quán)重向量為：W＝（0.387 3，0.139 7，0.274 8，0.198 1）T，其中λmax＝4.121 7，IC＝0.040 6，IR＝0.89，RC＝0.045 6＜0.1，判斷矩陣的一致性是可以接受的。

表5 GAMILLER量化理論

將各目標(biāo)的威脅度轉(zhuǎn)換為5個(gè)威脅等級(jí)，如表6所示。

表6 各目標(biāo)相當(dāng)于艦艇編隊(duì)威脅計(jì)算結(jié)果

通過(guò)表6可知，敵巡航導(dǎo)彈的威脅程度最大，威脅等級(jí)最高，應(yīng)首先對(duì)其實(shí)施射擊。其他依次為F－35低空、F－18低空、F－35高空和F－18高空。

2 火力分配的數(shù)學(xué)模型

假設(shè)一個(gè)編隊(duì)具有m組不同類型的武器系統(tǒng)（也可以理解為m個(gè)導(dǎo)彈武器發(fā)控中心），編隊(duì)通過(guò)其偵察預(yù)警系統(tǒng)發(fā)現(xiàn)有n個(gè)敵來(lái)襲目標(biāo)。第i組武器系統(tǒng)由mi（i＝1，2，…，m）個(gè)單位組成?？罩衼?lái)襲目標(biāo)的威脅度系數(shù)為ωj，易損性系數(shù)為eij。其中eij（i＝1，2，…，m；j＝1，2，…，n）為第i組武器系統(tǒng)對(duì)第j個(gè)目標(biāo)的單位毀傷概率。導(dǎo)彈武器的火力分配可用矩陣X＝（xij）m×n來(lái)描述，其中xij是用于第j個(gè)目標(biāo)的第i組導(dǎo)彈武器的數(shù)量。

最優(yōu)分配的目標(biāo)函數(shù)是使對(duì)目標(biāo)的毀傷效能指標(biāo)達(dá)到最大。通常一定數(shù)量的第i類武器毀傷第j個(gè)目標(biāo)的概率可表示為：

因而，所有m類武器對(duì)目標(biāo)j的毀傷概率：

毀傷目標(biāo)數(shù)的數(shù)學(xué)期望值M為：

所以，最優(yōu)化火力分配模型為：

將上述非線性的整數(shù)規(guī)劃問(wèn)題轉(zhuǎn)化為能夠用匈牙利算法求解的線性整數(shù)規(guī)劃問(wèn)題：

式中：cij＝ωjeij。

3 艦空導(dǎo)彈火力分配

假設(shè)我艦艇編隊(duì)某型導(dǎo)彈驅(qū)逐艦具有2種不同類型的艦空導(dǎo)彈，攻擊敵威脅目標(biāo)，第1種導(dǎo)彈具有3個(gè)單位，第2種導(dǎo)彈的數(shù)量為2個(gè)單位。目標(biāo)的威脅度和每種導(dǎo)彈毀傷各個(gè)目標(biāo)的概率如表7所示，要求攻擊每個(gè)目標(biāo)的導(dǎo)彈數(shù)不超過(guò)2個(gè)，研究艦空導(dǎo)彈分配方案，使得目標(biāo)毀傷效能指標(biāo)達(dá)到最大。

根據(jù)上節(jié)可得各威脅目標(biāo)的威脅度，結(jié)合2型艦空導(dǎo)彈對(duì)各威脅目標(biāo)的毀傷概率（見(jiàn)表7），設(shè)xij為第i（i＝1，2）種導(dǎo)彈分配給第j（j＝1，2，…，5）個(gè)目標(biāo)的數(shù)量，則目標(biāo)函數(shù)為：

表7 目標(biāo)的重要程度和導(dǎo)彈毀傷概率

同時(shí)，還應(yīng)滿足下列約束條件：

目標(biāo)威脅度ωj＝（0.756 8，0.372 2，0.508 3，0.474 4，0.610 4）

矩陣R運(yùn)用匈牙利算法［10－12］得到最優(yōu)解矩陣：

艦艇編隊(duì)艦空導(dǎo)彈最優(yōu)火力分配方案見(jiàn)表8。

表8 艦空導(dǎo)彈火力分配方案

4 結(jié)束語(yǔ)

空中多目標(biāo)對(duì)艦艇編隊(duì)攻擊的威脅評(píng)估和排序是目標(biāo)分配和火力分配的基礎(chǔ)。本文針對(duì)艦艇編隊(duì)的空中多目標(biāo)威脅因素進(jìn)行分析，對(duì)來(lái)襲目標(biāo)威脅程度進(jìn)行了評(píng)估，并根據(jù)目標(biāo)的威脅度和艦空導(dǎo)彈對(duì)各目標(biāo)的毀傷概率，得到艦空導(dǎo)彈對(duì)威脅目標(biāo)的最優(yōu)火力分配方案。

［1］蘇曼，康鳳舉，謝攀.艦艇多目標(biāo)威脅評(píng)估及可視化方法綜合 研 究 ［J］.科 學(xué) 技 術(shù) 與 工 程，2010，10（8）：1893－1897.

［2］曹睿婷，章衛(wèi)國(guó)，李廣文.基于多層次分析方法的目標(biāo)威脅估計(jì)研究［J］.彈箭與制導(dǎo)學(xué)報(bào)，2011，31（2）：30－32.

［3］郭輝，徐浩軍，周莉.基于區(qū)間數(shù)TOPSIS法的空襲目標(biāo)威脅評(píng)估［J］.空軍工程大學(xué)學(xué)報(bào)，2011，12（1）：40－45.

［4］王小藝，劉載文，侯朝楨.基于模糊多屬性決策的目標(biāo)威脅估計(jì)方法［J］.控制與決策，2007，22（8）：859－32.

［5］耿濤，盧廣山，張安.基于Vague集的空中目標(biāo)威脅評(píng)估群決策方法［J］.系統(tǒng)工程與電子技術(shù)，2011，33（12）：2686－2690.

［6］尹高揚(yáng)，周紹磊，張文廣.艦艇對(duì)來(lái)襲反艦導(dǎo)彈威脅評(píng)估算法［J］.兵工自動(dòng)化，2011，30（12）：1－4.

［7］張年春，沈培華.編隊(duì)導(dǎo)彈攻擊的火力分配研究［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2007（3）：36－39.

［8］依鳴.中國(guó)艦艇編隊(duì)編隊(duì)中的驅(qū)逐艦［J］.兵工科技，2013（S0）：27－35.

［9］張洋.多目標(biāo)防空火力分配技術(shù)研究［D］.大連：大連大學(xué)，2012.

［10］黃力偉，許品剛，王勤.基于匈牙利算法求解的火力分配問(wèn)題［J］.火力與指揮控制，2007，32（6）：25－28.

［11］柳毅，佟明安.匈牙利算法在多目標(biāo)分配中的應(yīng)用［J］.火力與指揮控制，2002，27（4）：34－37.

［12］王勇，吳梅，侯俊，潘泉.基于匈牙利算法的協(xié)同作戰(zhàn)排序問(wèn)題最優(yōu)解［J］.彈箭與制導(dǎo)學(xué)報(bào)，2005，25（3）：209－211.