張 會(huì), 王 安, 程 健, 屈澤均
(1.海軍潛艇學(xué)院,山東 青島 266199;2.青島大學(xué) 自動(dòng)化學(xué)院, 山東 青島 266071)
搜索行動(dòng)決策是艦船海上作戰(zhàn)行動(dòng)決策的重要內(nèi)容。根據(jù)任務(wù)背景和行動(dòng)策略的不同,艦船海上搜索行動(dòng)主要有檢查搜索、應(yīng)召搜索和巡邏搜索3種。直線往返搜索作為巡邏搜索的一種基本樣式,因?yàn)榻M織簡(jiǎn)單,在部隊(duì)演習(xí)訓(xùn)練和海上封控行動(dòng)籌劃中運(yùn)用廣泛。而作為作戰(zhàn)雙方博弈的重點(diǎn),更復(fù)雜的巡邏搜索樣式和搜索策略則是作戰(zhàn)決策理論研究的重點(diǎn)。
可靠的搜索效能評(píng)估是搜索行動(dòng)優(yōu)化決策的基礎(chǔ)。艦船搜索效能評(píng)估方法主要包括仿真方法和解析方法兩類。已有的艦船搜索行動(dòng)效能分析方法包括單艦船往返搜索發(fā)現(xiàn)概率計(jì)算的仿真方法[1]、單艦船往返搜索發(fā)現(xiàn)概率計(jì)算的解析方法[2],以及其他樣式的搜索行動(dòng)效能評(píng)估的仿真和解析方法[3-4]。
隨著協(xié)同作戰(zhàn)和無人集群系統(tǒng)技術(shù)的發(fā)展,編隊(duì)和協(xié)同搜索行動(dòng)決策成為搜索行動(dòng)決策技術(shù)研究的熱點(diǎn)。為降低研究難度,一般使用仿真方法構(gòu)建編隊(duì)和協(xié)同搜索效能評(píng)估系統(tǒng)[5-6],采用這種方法可實(shí)現(xiàn)對(duì)編隊(duì)樣式和協(xié)同策略的靈活設(shè)計(jì)。由于仿真方法計(jì)算復(fù)雜度高,當(dāng)搜索平臺(tái)數(shù)量較多時(shí)只能通過人為的模式將搜索優(yōu)化限制在相對(duì)很小的空間內(nèi),難以保證方案優(yōu)化結(jié)果的可靠。
為實(shí)現(xiàn)高效的編隊(duì)和協(xié)同搜索行動(dòng)效能評(píng)估,出現(xiàn)了以單平臺(tái)搜索效能計(jì)算解析方法為基礎(chǔ)的多平臺(tái)協(xié)同搜索效能評(píng)估和方案優(yōu)化方法[7],這種方法的研究對(duì)單平臺(tái)解析計(jì)算方法的可靠性提出了更高的要求。已有單平臺(tái)搜索行動(dòng)效能解析計(jì)算方法多為近似方法,計(jì)算結(jié)果的可靠性沒有保證,因此限制了其在編隊(duì)和協(xié)同搜索行動(dòng)效能評(píng)估和方案優(yōu)化方法研究中的應(yīng)用。
作為多平臺(tái)搜索效能評(píng)估方法研究的基礎(chǔ),本文建立模型進(jìn)行單艦船往返攔截搜索發(fā)現(xiàn)概率計(jì)算公式推導(dǎo),并使用仿真方法對(duì)其正確性進(jìn)行檢驗(yàn)。
攔截搜索,也稱線搜索,即在作戰(zhàn)海區(qū)設(shè)置攔截陣地,艦船在陣地內(nèi)按照一定的樣式機(jī)動(dòng)搜索,意圖發(fā)現(xiàn)從該作戰(zhàn)海區(qū)經(jīng)過的目標(biāo)。為便于分析,進(jìn)行攔截搜索效能分析時(shí)可將搜索者的機(jī)動(dòng)過程簡(jiǎn)化為分段的勻速直線運(yùn)動(dòng)。
攔截搜索是一種常見的對(duì)運(yùn)動(dòng)目標(biāo)的搜索行動(dòng)。如果事先能根據(jù)目標(biāo)任務(wù)方向、類型和作戰(zhàn)海區(qū)水文條件等確定目標(biāo)的機(jī)動(dòng)方向、機(jī)動(dòng)速度以及目標(biāo)的通過位置和通過時(shí)間范圍,則可將攔截搜索效能分析問題抽象為一種典型的攔截搜索效能分析問題,即目標(biāo)航向和速度已知,目標(biāo)通過位置和通過時(shí)間均在一定區(qū)間內(nèi)均勻分布情形下的攔截搜索發(fā)現(xiàn)概率計(jì)算問題,稱為典型攔截搜索發(fā)現(xiàn)概率計(jì)算問題。針對(duì)典型攔截搜索發(fā)現(xiàn)概率計(jì)算問題所進(jìn)行的研究,可在一定程度上揭示攔截搜索行動(dòng)的一般規(guī)律。
對(duì)于典型攔截搜索,當(dāng)搜索樣式為單兵力往返搜索時(shí),有一個(gè)常用的效能計(jì)算公式[2],本文使用“基于有效相對(duì)搜索區(qū)分析” 的方法對(duì)該公式的計(jì)算原理進(jìn)行分析,并編制仿真計(jì)算程序?qū)ζ溥M(jìn)行檢驗(yàn),發(fā)現(xiàn)其計(jì)算結(jié)果誤差較大。為消除或減小計(jì)算誤差,通過進(jìn)一步分析對(duì)計(jì)算公式進(jìn)行修正,使用仿真方法驗(yàn)證修正后公式的正確性。最后對(duì)使用“基于有效相對(duì)搜索區(qū)分析” 的方法進(jìn)行一般樣式或多艦船協(xié)同的攔截搜索發(fā)現(xiàn)概率計(jì)算的難點(diǎn)進(jìn)行分析,給出了開展研究工作的思路。
所謂經(jīng)典往返攔截搜索,指的是搜索者的搜索線為一條與目標(biāo)機(jī)動(dòng)方向相垂直的直線段,搜索者在搜索線上往返機(jī)動(dòng),搜索者搜索線覆蓋范圍與目標(biāo)通過搜索線的位置范圍一致,搜索者搜索時(shí)間范圍覆蓋目標(biāo)通過搜索線時(shí)間,目標(biāo)通過位置和通過時(shí)間均在各自可能取值范圍內(nèi)均勻分布的一種攔截搜索情形。
經(jīng)典往返攔截搜索的搜索過程如圖1所示,目標(biāo)以某速度自上向下機(jī)動(dòng),設(shè)置一條與目標(biāo)機(jī)動(dòng)方向垂直、長(zhǎng)度為L(zhǎng)的搜索線AB,搜索艦船以某速度在搜索線AB上往返機(jī)動(dòng),搜索艦船發(fā)現(xiàn)目標(biāo)當(dāng)且僅當(dāng)目標(biāo)在某個(gè)時(shí)刻進(jìn)入艦船探測(cè)扇面。作為經(jīng)典搜索,假設(shè)目標(biāo)通過搜索線的位置在搜索線AB上均勻分布,通過時(shí)間在區(qū)間[Tx,Ts]內(nèi)均勻分布,其中Tx、Ts分別為目標(biāo)通過搜索線的最早和最晚時(shí)間。
圖1 往返巡邏攔截搜索過程示意圖
運(yùn)用相對(duì)運(yùn)動(dòng)原理,可以將典型情形下的攔截搜索發(fā)現(xiàn)概率計(jì)算轉(zhuǎn)化為求“搜索艦船有效相對(duì)搜索區(qū)面積與目標(biāo)可能位置區(qū)域面積的比值”,然后通過對(duì)搜索者有效搜索區(qū)分析得到攔截搜索發(fā)現(xiàn)概率的解析計(jì)算公式。
將搜索者的探測(cè)范圍簡(jiǎn)化為一個(gè)圓形區(qū)域,則搜索者機(jī)動(dòng)一個(gè)直航段,其探測(cè)圓掃過的區(qū)域?yàn)橐粋€(gè)膠囊形區(qū)域,稱為搜索者該直航段的搜索區(qū),如圖2所示。圖2中,搜索者對(duì)目標(biāo)的探測(cè)距離為ds,直線段P1P2為搜索者的一個(gè)直線航跡段。搜索者航段P1P2的搜索區(qū)為一個(gè)半徑為ds的半圓Q1W1Q2、一個(gè)寬度為2ds的矩形Q1Q2Q3Q4和一個(gè)半徑為ds的半圓Q3W2Q4組成的膠囊型區(qū)域。
圖2 單航段搜索區(qū)示意圖
搜索艦船機(jī)動(dòng)多個(gè)直航段,所形成的搜索區(qū)就是所有其單航段搜索區(qū)的并,如果是多艦船聯(lián)合搜索,則搜索區(qū)就是全部艦船的所有單航段搜索區(qū)的并。
顯然,不論是單艦船搜索還是多艦船聯(lián)合搜索,如果目標(biāo)不動(dòng),搜索者能發(fā)現(xiàn)目標(biāo)當(dāng)且僅當(dāng)目標(biāo)位于搜索者搜索區(qū)內(nèi)。如果目標(biāo)靜止不動(dòng)且位置在某個(gè)區(qū)域內(nèi)均勻分布,則搜索者發(fā)現(xiàn)目標(biāo)的概率等于“搜索者搜索區(qū)位于目標(biāo)可能位置區(qū)域內(nèi)部分的面積”與“目標(biāo)可能位置區(qū)域面積”的比值。為簡(jiǎn)潔,以下稱“搜索者搜索區(qū)位于目標(biāo)可能位置區(qū)域內(nèi)部分”為“搜索者的有效搜索區(qū)”。
為了將搜索區(qū)概念用于攔截搜索發(fā)現(xiàn)概率計(jì)算,運(yùn)用相對(duì)運(yùn)動(dòng)原理對(duì)攔截搜索過程進(jìn)行分析。具體就是“將目標(biāo)看成不動(dòng),分析搜索者的相對(duì)運(yùn)動(dòng)所形成的搜索區(qū),即相對(duì)搜索區(qū)”。
(1)往返攔截搜索相對(duì)搜索區(qū)和目標(biāo)可能位置區(qū)域
以圖1中情形為例,將搜索者某次位于搜索線左端點(diǎn)時(shí)刻定義為0時(shí)刻,將目標(biāo)看成在其0時(shí)刻位置不動(dòng),繪制出搜索者的相對(duì)機(jī)動(dòng)航跡段及各航段搜索區(qū)如圖3所示。圖3左側(cè)從下向上的兩個(gè)相對(duì)運(yùn)動(dòng)速度合成三角形,分別為搜索者從左向右機(jī)動(dòng)和從右向左機(jī)動(dòng)時(shí)的搜索者相對(duì)運(yùn)動(dòng)速度合成三角形,其中,Vw、Vm、Vx分別為搜索者運(yùn)動(dòng)速度向量、目標(biāo)運(yùn)動(dòng)速度向量和搜索者的相對(duì)運(yùn)動(dòng)速度向量。
從0時(shí)刻開始,搜索者從左向右機(jī)動(dòng)時(shí),由相對(duì)運(yùn)動(dòng)速度合成三角形可以得到,其相對(duì)運(yùn)動(dòng)方向?yàn)閛A1方向,到達(dá)右端點(diǎn)后,搜索者改為由右向左機(jī)動(dòng),相對(duì)運(yùn)動(dòng)方向變?yōu)锳1A2方向,然后搜索者再?gòu)淖笙蛴覚C(jī)動(dòng),……,如此往復(fù),搜索者的相對(duì)搜索區(qū)為由這些直航段所形成的膠囊型區(qū)域組成的多個(gè)膠囊型區(qū)域的并。從0時(shí)刻開始的3個(gè)搜索單程,搜索者的相對(duì)運(yùn)動(dòng)航跡為折線oA1A2A3,相對(duì)搜索區(qū)域則為圖3中的3個(gè)膠囊型區(qū)域的并。
圖3 相對(duì)有效搜索區(qū)和目標(biāo)可能位置區(qū)域示意圖
如圖3所示,以搜索線左端點(diǎn)為坐標(biāo)原點(diǎn),以水平向右的方向?yàn)閤軸正向,建立平面直角坐標(biāo)系,則目標(biāo)0時(shí)刻位置橫坐標(biāo)取值范圍為[0,L],再由目標(biāo)通過巡邏線時(shí)刻范圍[Tx,Ts]推算目標(biāo)0時(shí)刻位置縱坐標(biāo)范圍為[VmTx,VmTs],故目標(biāo)0時(shí)刻的可能位置區(qū)域?yàn)橐粋€(gè)矩形,如圖3中矩形oE1E3E2所示,該矩形面積為
S=L(Ts-Tx)Vm
(1)
由對(duì)典型攔截搜索目標(biāo)通過位置和通過時(shí)間的假設(shè)可知,目標(biāo)0時(shí)刻位置在矩形oE1E3E2內(nèi)均勻分布。
如果目標(biāo)通過時(shí)間范圍區(qū)間長(zhǎng)度為其一個(gè)巡邏單程時(shí)間的m倍,則搜索艦船一個(gè)巡邏單程對(duì)應(yīng)的目標(biāo)可能位置區(qū)域?yàn)?/p>
(2)
式中T為搜索艦船一個(gè)巡邏單程時(shí)間。
(2)往返攔截搜索發(fā)現(xiàn)概率計(jì)算公式推導(dǎo)
下面分析搜索者的有效搜索區(qū)面積,也就是相對(duì)搜索區(qū)位于目標(biāo)可能位置區(qū)域內(nèi)的面積。觀察搜索者0時(shí)刻后的第一個(gè)搜索單程的相對(duì)搜索區(qū),即相對(duì)運(yùn)動(dòng)航跡線段oA1的搜索區(qū),其位于目標(biāo)可能位置區(qū)域內(nèi)的部分為圖4中的曲邊多邊形區(qū)域D1C1B2D2C4B3。該圖形面積近似等于矩形B1B2B4B3的面積,事實(shí)上,矩形B1B2B4B3比有效搜索區(qū)D1C1B2D2C4B3多了兩個(gè)面積相等的三角形oB1D1和A1D2B4,少了兩個(gè)面積相等的扇形oC1B2和A1C4B3,鑒于三角形的面積和扇形的面積相差不大,將偏差忽略。
圖4 單航段有效相對(duì)搜索區(qū)面積估算原理示意圖
矩形B1B2B4B3的面積為直線段oA1的搜索區(qū)的矩形部分,即
(3)
綜合以上分析可得,當(dāng)m為整數(shù)時(shí)典型攔截搜索發(fā)現(xiàn)概率的解析計(jì)算公式為
(4)
整理得
(5)
這是一個(gè)經(jīng)常被使用的攔截搜索發(fā)現(xiàn)概率計(jì)算公式[2]。該公式被用于進(jìn)行多平臺(tái)線式搜索發(fā)現(xiàn)概率估算和搜索方案優(yōu)化決策[5]。
因?yàn)榫幊淘砗?jiǎn)單、不易出錯(cuò),仿真方法經(jīng)常被用于計(jì)算效率要求不高情形下各種搜索行動(dòng)的效能分析以及其他方法的正確性檢驗(yàn)[1,3,6-8]。為檢驗(yàn)式(5)計(jì)算結(jié)果的可靠性,編制計(jì)算機(jī)程序進(jìn)行攔截搜索發(fā)現(xiàn)概率仿真計(jì)算對(duì)比。
(1)仿真檢驗(yàn)
取目標(biāo)速度18節(jié)(≈33.3 km/h),搜索艦船速度8節(jié)(≈14.8 km/h)、搜索線長(zhǎng)度12海里(≈22.2 km)、對(duì)目標(biāo)探測(cè)距離3海里(≈5.6 km)。取目標(biāo)可能通過位置范圍與艦船巡邏線覆蓋位置范圍一致,用式(5)計(jì)算得到的艦船發(fā)現(xiàn)目標(biāo)的概率為0.547 2。取目標(biāo)通過時(shí)間范圍為0~72 h,樣本容量10 000,仿真時(shí)間步長(zhǎng)10 s,仿真結(jié)果為0.471 0,二者偏差較大,式(5)的計(jì)算結(jié)果偏大。
取目標(biāo)速度8節(jié)(≈14.8 km/h),其他參數(shù)值保持不變,用式(5)計(jì)算得到的發(fā)現(xiàn)概率為0.707 1,而仿真結(jié)果為0.582 2,式(5)的計(jì)算結(jié)果還是偏大。
(2)計(jì)算公式修正
為找出式(5)計(jì)算結(jié)果偏大的原因,對(duì)圖3中搜索者有效搜索區(qū)進(jìn)行仔細(xì)觀察發(fā)現(xiàn),除了前面提及的計(jì)算單航段有效搜索區(qū)面積所用矩形面積實(shí)際上比搜索者單航段有效搜索面積略大外,搜索者相鄰兩個(gè)單程搜索區(qū)之間還有重疊,進(jìn)行多階段搜索有效搜索區(qū)面積計(jì)算時(shí),沒有將重疊部分的面積減去??紤]這兩點(diǎn),對(duì)單航程有效搜索區(qū)面積計(jì)算公式進(jìn)行以下兩項(xiàng)修正:
最后得到,一個(gè)巡邏單程搜索艦船的有效搜索面積實(shí)際為
(6)
目標(biāo)可能位置區(qū)域不變,故攔截搜索發(fā)現(xiàn)概率應(yīng)該為
(7)
針對(duì)前述的計(jì)算情形1,使用式(7)得到的發(fā)現(xiàn)概率計(jì)算結(jié)果為0.472 3,與仿真結(jié)果相吻合。針對(duì)計(jì)算情形2,使用式(7)得到的發(fā)現(xiàn)概率計(jì)算結(jié)果0.582 2,與仿真結(jié)果同樣吻合。
取目標(biāo)速度18節(jié)(≈33.3 km/h),搜索者搜索速度8節(jié)(≈14.8 km/h)、對(duì)目標(biāo)探測(cè)距離3海里(≈5.6 km)、巡邏線長(zhǎng)度范圍10~20海里(≈18.52~37.04 km),使用式(5)、式(7)以及仿真分別計(jì)算得到不同巡邏線長(zhǎng)度時(shí)搜索者發(fā)現(xiàn)目標(biāo)的概率值,并將計(jì)算結(jié)果繪制成曲線,如圖5所示。圖5中,橫坐標(biāo)為巡邏線長(zhǎng)度,單位海里(1海里≈1.852 km),縱坐標(biāo)為搜索者發(fā)現(xiàn)目標(biāo)的概率,取值較大的紅色光滑曲線為式(5)的計(jì)算結(jié)果,值較小的綠色光滑曲線為式(7)的計(jì)算結(jié)果,值較小的藍(lán)色非光滑曲線為仿真法所得結(jié)果。由圖5可以看出,修正后的式(7)計(jì)算結(jié)果與仿真結(jié)果吻合,說明其是正確的。
圖5 不同方法計(jì)算所得發(fā)現(xiàn)概率曲線比較圖
式(7)是假設(shè)目標(biāo)速度和航向確定情況下的往返攔截搜索發(fā)現(xiàn)概率計(jì)算公式。實(shí)際的作戰(zhàn)行動(dòng)決策過程中,會(huì)遇到目標(biāo)速度和目標(biāo)航向同樣具有不確定性以及平臺(tái)搜索樣式更加靈活的情形,能否將該公式用于這些情況下的發(fā)現(xiàn)概率估算,或者是否可將該公式的研究方法推廣應(yīng)用于這些情形下的發(fā)現(xiàn)概率計(jì)算公式推導(dǎo)呢?為回答這個(gè)問題,進(jìn)行計(jì)算公式對(duì)目標(biāo)信息不確定性的適應(yīng)性和研究思路對(duì)攔截搜索樣式變化適應(yīng)性的分析。
(1)計(jì)算公式對(duì)目標(biāo)信息不確定性的適應(yīng)性分析
假設(shè)除了目標(biāo)通過位置和通過時(shí)間的不確定性外,目標(biāo)速度和航向也是不能完全確定的[9-10]。將目標(biāo)速度和目標(biāo)航向均描述為隨機(jī)變量,使用仿真確定搜索行動(dòng)發(fā)現(xiàn)概率,與使用式(7)計(jì)算所得結(jié)果進(jìn)行比較。
取巡邏線長(zhǎng)度10海里(≈18.52 km),其他參數(shù)同前,假設(shè)目標(biāo)速度和航向均服從均勻分布,取不同的目標(biāo)速度分布均方差和目標(biāo)航向分布均方差,得到的發(fā)現(xiàn)概率計(jì)算偏差隨目標(biāo)速度散布均方差變化曲線如圖6所示,發(fā)現(xiàn)概率計(jì)算偏差隨目標(biāo)航向散布均方差變化曲線如圖7所示。
由圖6和圖7可以看出,當(dāng)目標(biāo)速度均方差和目標(biāo)航向散布均方差較小時(shí),計(jì)算結(jié)果偏差較小,故一定情況下仍可使用式(7)進(jìn)行發(fā)現(xiàn)概率計(jì)算。
比較圖6和圖7還可發(fā)現(xiàn),對(duì)比計(jì)算偏差隨目標(biāo)速度均方差的變化,計(jì)算偏差隨目標(biāo)航向均方差的變化率更大,且偏差值為正值,即目標(biāo)航向不確定情形下式(7)的計(jì)算結(jié)果偏大,分析其原因,是因?yàn)椴糠謽颖緦?duì)應(yīng)的目標(biāo)0時(shí)刻位置超出理論分析所使用的目標(biāo)可能位置區(qū)域,理論計(jì)算所使用的目標(biāo)可能位置區(qū)域小于實(shí)際的目標(biāo)位置區(qū)域,因此所得概率值偏大。
圖6 發(fā)現(xiàn)概率計(jì)算偏差隨目標(biāo)速度散布均方差變化曲線
圖7 發(fā)現(xiàn)概率計(jì)算偏差隨目標(biāo)航向散布均方差變化曲線
(2)研究思路對(duì)攔截搜索樣式變化的適應(yīng)性
為滿足戰(zhàn)法設(shè)計(jì)靈活性的需要,需要給出任意樣式的攔截搜索發(fā)現(xiàn)概率計(jì)算的方法,以及更多情形下的多實(shí)體協(xié)同搜索效能分析方法[11-13]。若使用本文所用“基于有效相對(duì)搜索區(qū)分析” 的方法進(jìn)行任意樣式的攔截搜索發(fā)現(xiàn)概率計(jì)算,需要進(jìn)行數(shù)量較多且相互關(guān)系多樣的單航段搜索區(qū)合并圖形面積的計(jì)算。
由本文對(duì)單艦船往返攔截搜索發(fā)現(xiàn)概率公式的推導(dǎo)和修正過程可以看出,由于涉及不同航段搜索區(qū)合并圖形面積的計(jì)算,即便是對(duì)單艦船的非常簡(jiǎn)單的搜索樣式,攔截搜索效能計(jì)算公式的推導(dǎo)也是比較困難且極易出錯(cuò)的。相比往返搜索,其他樣式搜索的搜索區(qū)更為復(fù)雜,而多艦船聯(lián)合搜索時(shí),包含的不同直線航段數(shù)更多,形成的聯(lián)合搜索區(qū)更難以直接給出面積計(jì)算公式。
要實(shí)現(xiàn)對(duì)任意樣式攔截搜索發(fā)現(xiàn)概率的高效計(jì)算,需要研究專門的聯(lián)合搜索區(qū)面積計(jì)算算法或設(shè)法提高仿真方法的計(jì)算效率。任意樣式攔截搜索發(fā)現(xiàn)概率計(jì)算方法研究是下一步的一項(xiàng)重點(diǎn)工作。
攔截搜索是反潛戰(zhàn)中重要的搜索樣式,本文使用“基于有效相對(duì)搜索區(qū)分析” 的方法對(duì)常用的單艦船往返攔截搜索發(fā)現(xiàn)概率計(jì)算公式進(jìn)行了理論推導(dǎo)和仿真檢驗(yàn),并根據(jù)檢驗(yàn)結(jié)果和進(jìn)一步的理論分析對(duì)其進(jìn)行修正,得到更為精確的計(jì)算公式,為任意樣式攔截搜索發(fā)現(xiàn)概率計(jì)算公式推導(dǎo)和算法設(shè)計(jì)提供了思路。根據(jù)戰(zhàn)場(chǎng)實(shí)際,除了典型攔截搜索行動(dòng)效能計(jì)算外,還需要研究目標(biāo)航向和速度不是確定值情形下的發(fā)現(xiàn)概率計(jì)算,這是下一步工作的另一個(gè)重點(diǎn)。