基于動態規劃的反艦導彈火力分配優化方法

粘松雷，嚴建鋼，陳 榕

(海軍航空工程學院，山東 煙臺 264001)

基于動態規劃的反艦導彈火力分配優化方法

粘松雷，嚴建鋼，陳 榕

(海軍航空工程學院，山東 煙臺 264001)

遠程精確打擊是反艦導彈攻擊的必然趨勢和高技術戰場的主要特點。導彈攻擊是打擊敵水面艦船的主要手段，如何充分發揮反艦導彈的作戰效能，減少不必要的浪費就成為戰前制定作戰方案以及實施作戰指揮時必須解決的首要問題，而火力分配則是其中的一個重要環節。

反艦導彈;火力分配;動態規劃

0 引言

“反艦導彈火力分配”是指根據導彈攻擊的任務，合理進行目標配對決策，并向導彈裝定相應的決策指令，依靠導彈自身的目標選擇能力選擇火分目標，實現上級的作戰意圖，保障作戰協同［1］。其實質是按照指揮員的作戰意圖，合理分配導彈火力，打擊預定目標。其實現形式主要體現在導彈的火分目標選擇能力上，依靠導彈自主完成對目標的搜索、識別、選擇和攻擊，提高導彈的智能化作戰水平。

1 火力分配數學模型的建立

1.1 作戰想定

紅方對藍方艦艇編隊實施超視距多彈飽和攻擊，藍方艦艇對來襲反艦導彈實施硬攔截和軟對抗;在攻防過程中不考慮藍方對紅方的導彈反擊［2］。

紅方編隊對藍方編隊導彈攻擊時的火力分配，需解決以下2個問題：

1)為了達到對藍方編隊預定的毀傷指標，如何合理部署紅方導彈兵力，才能使紅方消耗的導彈價值最小;

2)紅方如何合理分配導彈，才能使其對藍方編隊造成的毀傷程度最大，即使紅方導彈的作戰效能最大化［3］。

可見，紅方與藍方的攻防對抗問題是具有目標價值收益最大、耗彈量成本最少的約束條件眾多而復雜的優化問題，該問題是多目標和多約束的優化問題，即多目標優化。MO問題需要優化一組費用函數，其解不是單一點，而是一組點的集合，稱之為Pareto 最優集［4－5］。

1.2 數學模型的建立

紅方多種型號的反艦導彈對藍方艦艇編隊實施飽和攻擊。假設i為導彈的型號;λi為1枚第i類導彈的價值系數;j和r為目標的批號;Ni為第i類導彈的數量;xij為第i類導彈用于射擊第j批目標的數量;xj為射擊第j批目標的導彈總數;vr為第r批目標的價值系數;ωir為第i類導彈按照毀傷指標，毀傷第r批目標所需的導彈數;Pijr為第i類導彈對第j批目標射擊時對第r批目標的命中概率;ki為第i類導彈的價值系數。

根據作戰想定可以建立導彈攻擊火力分配問題的規劃模型如下：

式中：Pi，可靠為第i類導彈的可靠飛行概率;Pijr捕捉為第i類導彈對第j批目標射擊時，對第r批目標的捕捉概率;Pir自導為第i類導彈對第r批目標射擊時的自導命中概率;Qir，突防為第i類導彈對第r批目標射擊時的突防概率。ωir，輕傷為第i類導彈輕傷第r批目標平均所需的導彈數;ωir，重創為第i類導彈重創第r批目標平均所需的導彈數;ωir，擊沉為第i類導彈擊沉第r批目標平均所需的導彈數。

因為反艦導彈末制導雷達對某批目標的捕捉概率具有很大的不確定性。為此，在不考慮反艦導彈的目標捕捉概率Pijr捕捉對作戰效能影響的情況下，可將式(1)進一步簡化為：

2 基于動態規劃算法的火力分配優化方法

動態規劃是運籌學的一個分支，是求解決策過程最優化的數學方法［4］。20世紀50年代初，美國數學家R．E．Bellman等在研究多階段決策過程的優化問題時，提出了著名的最優化原理，把多階段過程轉化為一系列單階段問題，逐個求解，創立了解決這類過程優化問題的新方法——動態規劃［6］。

多階段決策過程是指這樣一類特殊的活動過程，它們可以按時間順序分解成若干相互聯系的階段，在每個階段都要做出決策，全部過程的決策是一個決策序列［7］。要使整個活動的總體效果達到最優的問題，稱為多階段決策問題。雖然動態規劃主要用于求解以時間劃分階段的動態過程的優化問題，但一些與時間無關的靜態規劃(如線性規劃、非線性規劃)，只要人為地引進時間因素，把它視為多階段決策過程，也可以用動態規劃方法方便地求解。

2.1 基本假設

在反艦導彈的數量較多、型號單一時，可采用動態規劃算法優化火力分配，使導彈群的整體作戰效能最大化。采用動態規劃算法進行火力分配優化，需基于如下3點假設：

1)對每個目標射擊時，不存在火力轉移問題;2)參與攻擊的反艦導彈均為同一型號;

3)反艦導彈總數量不足以對所有被攻擊目標艦艇同時實施飽和攻擊。

這3點假設對應的數學描述為

其中：i=1，2，…，k;j=1，2，…，m;r=1，2，…，m。雖然導彈火力分配問題是個多目標優化問題，但依據假設將參與攻擊的反艦導彈總數量常數化了，因此多目標優化問題即轉化為單目標函數的非線性規劃問題，如式(4)所示。

2.2 模型解算過程

動態規劃問題的求解過程共分5步：

1)根據所給問題，將過程恰當地分為若干階段(一般根據時間和空間來劃分)，即把導彈的火力分配過程看作是有次序地將導彈分配給第1、第2，…，第m個目標的過程，且把分配給第j個目標的過程作為第j個階段，這樣就把導彈火力分配過程轉化為相互聯系的m個階段。

2)正確選擇狀態變量來描述過程及其在各個階段的狀態。動態規劃中的狀態概念，必須滿足“無后效性這樣一個條件”。所謂無后效性是指，如果給定了某一階段的某一狀態，則其后續過程的發展不受該狀態以前各階段狀態的影響。這也就是說，過程的狀態不影響它未來的發展。

設狀態變量Sj為分配給第j個目標到第m個目標的導彈總數，即j階段初所擁有的導彈數。顯然，0≤Sj≤N。

3)確定決策變量及每段允許決策集合。

決策變量xj為分配給第j個目標的導彈數，0≤xj≤Sj(即允許決策集合)。

4)建立狀態轉移方程。

一般來說，對于j階段某一給定的狀態Sj，只要決策變量xj一經確定，其下一階段的狀態變量值Sj+1也就完全確定。因此，可建立Sj+1隨Sj和xj變化的對應關系，這種狀態轉移的規律，可用式(5)表示Sj+1與Sj和xj之間的函數關系，稱為狀態轉移方程，即

則該導彈火力分配問題的狀態轉移方程為

5)根據題意列出指標函數，該指標函數要滿足遞推性。

用階段毀傷效能函數gj(Sj，xj)表示對第j個目標分配xj枚導彈后所得到的毀傷效能，最優效能指標函數fj(Sj)為將Sj枚導彈分配給第j個目標到第m個目標后，所獲得的最優毀傷效能。其動態規劃的遞推關系為

按照以上5個步驟，用遞推方法，先求后部最優子策略，然后逐步向前遞推得全過程的最優策略。

3 實例分析

設紅方飛機編隊在一次突襲中使用8枚某型反艦導彈對藍方水面艦艇編隊實施攻擊，藍方編隊由1艘大型補給船T1、1艘驅逐艦T2、1艘護衛艦T3和另外1艘護衛艦T4組成。藍方艦艇的價值系數為vj，單枚XX型導彈對單艘目標艦艇的毀傷概率為ej，如表1所示。

依據遞推關系式(7)及表1的數據，求解最優火力分配方案，其中，導彈數N=8;目標數m=4。

表1 目標艦艇的價值系數和毀傷概率Tab．1 Important coefficients and kill probability of targets

1)當j=4時，S4=0，1，…8;X4=0，1，…8;該階段最大毀傷效能為

2)當j=3時，S3=0，1，…8;即S3枚導彈分配給第3和第4個目標，最大毀傷效能為

3)當j=2時，S2=0，1，…8;S2枚導彈分配給第2、第3和第4個目標，最大毀傷效能為

4)當j=1時，只有S1=8，即8枚導彈分配給第1、第2、第3和第4個目標，最大毀傷效能為

鑒于文章篇幅問題，這里僅給出j=1時的計算結果，如表2所示。

表2 j=1時的階段效能Tab．2 Stage efficiency when j=1

由表2可知，最大毀傷效能指標為f1(8)=3.231，此時決策變量x1=2，即對第1個目標發射2枚反艦導彈。然后，按j=2，j=3，j=4的順序逆推，從而求得導彈最優分配方案為：X1=2，X2=2，X3=2，X4=2。導彈的火力分配矩陣為［2，2，2，2］，即對大型補給船T1、驅逐艦T2、護衛艦T3和護衛艦T4分別發射2枚導彈，實施攻擊所獲得的綜合毀傷效能最大。

4 結語

本文在作戰想定的基礎上，建立反艦導彈群攻擊火力分配問題的完整數學模型;并運用動態規劃算法對火力分配的數學模型進行了優化研究，在實例中證明了數學模型和優化算法的正確和可行性。

［1］楊祖快，劉鼎臣．反艦導彈自主式火力分配［J］．火力與指揮控制，2005，30(2)：1 －4．

YANG Zu-kuai，LIU Ding-chen．Study on anti-ship missile firepower distribution［J］．Fire Control and Command Control，2005，30(2)：1 －4．

［2］滕克難．艦空導彈系統抗飽和攻擊反導作戰效能評估［J］．彈箭與制導學報，2005，25(3)：37 －39．

TENG Ken-nan．Effectiveness evalaution of anti-missile saturation attack for the ship-to － air missile［J］．Arrow and Guidance Transaction，2005，25(3)：37 －39．

［3］張年春，沈培華．編隊導彈攻擊的火力分配研究［J］．戰術導彈技術，2007，(3)：36 －39．

ZHANG Nian-chun， SHEN Pei-hua．Researchonfire disposition model for formated anti-ship missile［J］．Tactical Missile Technolog，2007，(3)：36 －39．

［4］王小平，曹立明．遺傳算法－理論、應用與軟件實現［M］．西安：西安交通大學出版社，2002.115－121．

WANG Xiao-ping，CHAO Li-ming．Inheritance algorithmtheory，application and software realization［M］．Xi'an：Xi'an Traffic University Publishing House，2002：115 －121．

［5］崔遜學．多目標進化算法及其應用［M］．北京：國防工業出版社，2006.40－50．

CUI Xun-xue．Many objectives evolution algorithm and application［M］．Beijing：National Defense Industry Publishing House，2006：40－50．

［6］張最良，李長生，等．軍事運籌學［M］．北京：軍事科學出版社，1999．50－56．

ZHANG Zui-liang， LIChang-sheng， etal． Military operations research［M］．Beijing：Military Science Publishing House，1999．50 －56．

［7］謝春燕，李為民，劉付顯．整數規劃與動態規劃在防空導彈火力分配中的應用研究［J］．彈箭與制導學報，2003，23(4)：72－75．

XIE Chun-yan， LIWei-min， LIU Fu-xian．Integer programming and dynamic programming to airforce missile firepower allocation［J］．Arrow and Guidance Transaction，2003，23(4)：72－75．

The optimization method of antiship missile fire distribution based on dynamic programming

NIAN Song-lei，YAN Jian-gang，CHEN Rong
(Naval Aeronautical Engineering Institute，Yantai 264001，China)

The long-range accurate stike is the certain trend of antiship missile attack and the major characteristic of high technology battlefield，the missile attack is the major means hitting enemy's surface ships，how to sufficiently elaborate the operation efficacy of antiship missileand reduce unnecessary waste is the first importance question which must be solved in course of action and operational command，while the fire distribution is the important link．

antiship missile;fire distribution;dynamic programming

E927

A

1672－7649(2012)07－0110－04

10.3404/j.issn.1672－7649.2012.07.025

2011－08－01;

2011－10－12

粘松雷(1982－)，男，博士研究生，研究方向為海軍兵種作戰數理戰術分析。