非協(xié)同方式多微小型飛行器區(qū)域覆蓋多目標(biāo)的彈道建模與設(shè)計(jì)

都業(yè)宏，郁浩，張軍，李杰，徐興建

（1.中國(guó)白城兵器試驗(yàn)中心，吉林白城137001；2.北京理工大學(xué)機(jī)電學(xué)院，北京100081）

非協(xié)同方式多微小型飛行器區(qū)域覆蓋多目標(biāo)的彈道建模與設(shè)計(jì)

都業(yè)宏1，2，郁浩1，張軍1，李杰2，徐興建2

（1.中國(guó)白城兵器試驗(yàn)中心，吉林白城137001；2.北京理工大學(xué)機(jī)電學(xué)院，北京100081）

為提高微小型飛行器（MAV）群對(duì)作戰(zhàn)目標(biāo)偵察攻擊的適應(yīng)能力，建立多MAV非協(xié)同方式區(qū)域覆蓋多目標(biāo)問(wèn)題模型。應(yīng)用云模型模擬MAV的起飛點(diǎn)坐標(biāo)值，應(yīng)用蒙特卡洛方法模擬起飛方向角度值，建立與問(wèn)題模型相關(guān)的目標(biāo)區(qū)域模型、MAV探測(cè)模型、起飛點(diǎn)模型、初始飛行方向模型和彈道模型。為提高優(yōu)化設(shè)計(jì)效率，在對(duì)期望覆蓋率指標(biāo)進(jìn)行了穩(wěn)定性分析的基礎(chǔ)上，確定計(jì)算期望覆蓋率的最小模擬次數(shù)。應(yīng)用多目標(biāo)遺傳算法對(duì)建立的多目標(biāo)問(wèn)題模型進(jìn)行了求解，得到了非協(xié)同方式條件下多MAV區(qū)域覆蓋的彈道設(shè)計(jì)方案。從MAV的探測(cè)半徑與目標(biāo)區(qū)域半徑的比例關(guān)系出發(fā)，對(duì)這些方案適于解決的問(wèn)題進(jìn)行了分析。該彈道降低了MAV控制難度，兼顧了覆蓋率、MAV數(shù)量和航程3個(gè)目標(biāo)。

兵器科學(xué)與技術(shù)；非協(xié)同方式；微小型飛行器；區(qū)域覆蓋；彈道設(shè)計(jì)；多目標(biāo)優(yōu)化

0 引言

為激發(fā)出單獨(dú)無(wú)人飛行器（UAV）所不具備的作戰(zhàn)效能或完成單獨(dú)UAV無(wú)法完成的任務(wù)，國(guó)內(nèi)外很多學(xué)者提出多UAV協(xié)同作戰(zhàn)方法，以衍生出UAV群系統(tǒng)的整體作戰(zhàn)效能，完成單獨(dú)UAV無(wú)法完成的任務(wù)［1-2］。借鑒UAV協(xié)同作戰(zhàn)研究成果，國(guó)內(nèi)外學(xué)者在微小型飛行器（MAV）協(xié)同作戰(zhàn)領(lǐng)域也取得了一些研究成果［3-7］，但這些成果大都是在進(jìn)行了一些假設(shè)條件下完成的。常見(jiàn)的或默認(rèn)的假設(shè)有：

1）通信鏈路完全可靠；

2）可在最小或近乎最小轉(zhuǎn)彎半徑下進(jìn)行滾轉(zhuǎn)和偏航協(xié)同轉(zhuǎn)彎；

3）定速飛行；

4）完全抗干擾、抗電子誘騙，或無(wú)干擾、無(wú)電子誘騙。

MAV作為一種小型UAV，它比常見(jiàn)的UAV更小、更輕、更具機(jī)動(dòng)能力。MAV與體積較大的UAV相比有較多的共性，但也有很多不同，并且與當(dāng)前某些研究成果的假設(shè)相矛盾，主要有：

1）MAV采用的通信硬件相對(duì)于UAV采用的通信硬件性能有很大差距，導(dǎo)致MAV間通信效能明顯低于UAV間通信效能；

2）MAV采用的主控硬件相對(duì)于UAV采用的主控硬件性能有很大差距，加之MAV強(qiáng)不確定性、強(qiáng)耦合性和強(qiáng)擾動(dòng)特性，MAV飛控系統(tǒng)承擔(dān)的任務(wù)更加艱巨；

3）通常MAV采用電動(dòng)力驅(qū)動(dòng)飛行方式，導(dǎo)致其承載能力有限、航程有限，由于其抗風(fēng)能力有限，故僅能在較低空域（通常不大于500 m）低速（通常不大于70 m/s）飛行；

4）MAV抗電磁干擾和電子誘騙能力相對(duì)于UAV有很大差距；

5）MAV實(shí)時(shí)處理MAV間作戰(zhàn)信息能力相對(duì)于UAV差別巨大。

Kress等［8］針對(duì)不同的目標(biāo)識(shí)別率條件下進(jìn)行了UCAV的協(xié)同作戰(zhàn)模擬，從統(tǒng)計(jì)結(jié)果來(lái)看，在目標(biāo)識(shí)別率較低的情況下進(jìn)行協(xié)同作戰(zhàn)是沒(méi)有意義的。MAV群組的非協(xié)同方式作戰(zhàn)與協(xié)同作戰(zhàn)的相同之處：均具有群組作戰(zhàn)特征，即能夠勝任單架MAV無(wú)法完成的任務(wù)。MAV群組的非協(xié)同方式作戰(zhàn)與協(xié)同作戰(zhàn)的主要不同之處在于：

1）MAV間無(wú)通信數(shù)據(jù)交互，可以避免外界電子干擾與欺騙，可規(guī)避 MAV群組作戰(zhàn)效能減弱或喪失；

2）MAV群組采用非協(xié)同作戰(zhàn)方式時(shí)，MAV單機(jī)按預(yù)先設(shè)定作戰(zhàn)計(jì)劃工作，不需要兼顧MAV間的信息互通、MAV與MAV群組間的信息共享，降低了MAV控制系統(tǒng)的負(fù)擔(dān)，使MAV飛行控制系統(tǒng)設(shè)計(jì)更加簡(jiǎn)潔、可靠；

3）MAV采用離線(xiàn)工作模式，不需要安裝無(wú)線(xiàn)通信硬件設(shè)備，大幅降低了MAV研制成本；

4）協(xié)同能夠?qū)?dòng)態(tài)威脅環(huán)境進(jìn)行一定程度的反應(yīng)，而非協(xié)同方式則難以做出相應(yīng)的反應(yīng)。

而且，MAV群組采用非協(xié)同方式作戰(zhàn)時(shí)，需要MAV群組在某一空間點(diǎn)投放并完成彈機(jī)轉(zhuǎn)換和初始轉(zhuǎn)平飛過(guò)程，MAV的平飛起點(diǎn)和初始航向具有隨機(jī)性［9］；載荷有限，完成任務(wù)不需要回收，“打出去不管”［10］。

通常情況下，在建立彈道的多目標(biāo)模型時(shí)，UAV介入目標(biāo)區(qū)域的位置和飛行方向是確定的，完成覆蓋任務(wù)的時(shí)間也是按照這樣的條件進(jìn)行計(jì)算的。然而，由于MAV具有前述的兩個(gè)特點(diǎn)，就需要布撒器遠(yuǎn)距離拋撒MAV，MAV完成彈機(jī)轉(zhuǎn)換和初始轉(zhuǎn)平飛過(guò)程后，經(jīng)過(guò)一段dubins路徑［11］飛行至為其設(shè)計(jì)的覆蓋彈道起始位置，這樣就增加了完成覆蓋任務(wù)的時(shí)間。

從戰(zhàn)場(chǎng)環(huán)境方面考慮，高新技術(shù)應(yīng)用于現(xiàn)代戰(zhàn)爭(zhēng)的水平差異，已導(dǎo)致非對(duì)稱(chēng)作戰(zhàn)在一定程度、某些方面有所呈現(xiàn)，軍事博弈致使戰(zhàn)場(chǎng)環(huán)境日趨非線(xiàn)性和復(fù)雜化，并具有隨機(jī)性、模糊性、強(qiáng)對(duì)抗性、欺騙性等不確定性特點(diǎn)。在當(dāng)前復(fù)雜的戰(zhàn)場(chǎng)環(huán)境對(duì)抗中，強(qiáng)電磁干擾和欺騙成為應(yīng)對(duì)MAV群組協(xié)同作戰(zhàn)的一種有力武器；受自身載荷和航程的限制，需要MAV配備小體積、輕重量、低功耗、高可靠、抗干擾、低成本的通信器件，這類(lèi)通信器件仍亟待研制。

因此，為提高M(jìn)AV群對(duì)目標(biāo)區(qū)域覆蓋的速度，及其在復(fù)雜戰(zhàn)場(chǎng)環(huán)境下的適應(yīng)能力，降低MAV其控制系統(tǒng)負(fù)擔(dān)，本文充分利用MAV的工作特點(diǎn)，設(shè)計(jì)了非協(xié)同方式的 MAV群組區(qū)域覆蓋的彈道方案。

1 云模型

為對(duì)MAV起飛點(diǎn)和初始航向進(jìn)行模擬，本文引入了云模型方法模擬MAV的起飛點(diǎn)，引入了蒙特卡洛方法［12］模擬MAV的初始航向。下面重點(diǎn)對(duì)云模型方法進(jìn)行簡(jiǎn)要介紹。

云模型［13］是李德毅院士提出的定性和定量轉(zhuǎn)換模型，文獻(xiàn)［13-14］中詳細(xì)介紹了云模型基礎(chǔ)理論。文獻(xiàn)［15］應(yīng)用云模型對(duì)武器系統(tǒng)故障進(jìn)行了預(yù)測(cè)，對(duì)正向云發(fā)生器和逆向云發(fā)生器算法進(jìn)行了詳細(xì)敘述，并用Mat1ab語(yǔ)言對(duì)這兩種算法進(jìn)行了實(shí)現(xiàn)，因此本文不再贅述。文獻(xiàn)［16］應(yīng)用二維云模型對(duì)射擊精度合格率進(jìn)行了評(píng)價(jià)。

設(shè)U是一個(gè)精確數(shù)值集合的定量論域，C是U上的定性概念。若定量值x?U，且x是定性概念C的一次隨機(jī)實(shí)現(xiàn)，x對(duì)C的隸屬度μi?［0，1］是具有穩(wěn)定傾向的隨機(jī)數(shù)，則x在論域U上的分布成為云，每個(gè)x成為一個(gè)云滴x?U.

云模型的數(shù)字特征是描述云模型、產(chǎn)生虛擬云、實(shí)現(xiàn)云計(jì)算、完成云變換的數(shù)值基礎(chǔ)，也是利用云技術(shù)從含有不確定性的空間數(shù)據(jù)庫(kù)或空間數(shù)據(jù)倉(cāng)庫(kù)中發(fā)現(xiàn)空間知識(shí)的基礎(chǔ)。云的數(shù)字特征用期望值Ex、熵En、超熵He 3個(gè)數(shù)值表示，如圖1所示，其中：縱坐標(biāo)為云滴的隸屬度MG，橫坐標(biāo)為云滴的數(shù)值V；Ex反映了定性概念的云滴群的中心值；En在云模型中用來(lái)綜合度量定性概念的模糊度和概率，En的大小直接決定論域中可被模糊概念接受的范圍；He 是En的不確定度量，其大小間接地表示了云的離散程度和厚度。

取三維云模型的數(shù)字特征為

其1 000個(gè)云滴的效果在三維上展開(kāi)如圖2所示。

2 涉及到的模型

MAV彈道的設(shè)計(jì)涉及到目標(biāo)區(qū)域模型、MAV探測(cè)模型、起飛點(diǎn)模型、初始航向模型和彈道模型的建立。

圖1　云模型的數(shù)字特征Fig.1 Numerica1 characteristics of c1oud mode1

圖2　三維云模型示意圖Fig.2 Schematic diagram of three-dimensiona1 c1oud mode1

2.1目標(biāo)區(qū)域模型

目標(biāo)區(qū)域如圖3所示，O（XO，YO）為目標(biāo)區(qū)域的中心位置，R為目標(biāo)區(qū)域半徑。

圖3　圓形作戰(zhàn)區(qū)域Fig.3 Round combat zone

對(duì)戰(zhàn)場(chǎng)平面進(jìn)行柵格化處理，表示成節(jié)點(diǎn)（X，Y，p（X，Y）），｛（X，Y）|（X-XO）2+（Y-YO）2≤R2｝為戰(zhàn)場(chǎng)的柵格點(diǎn)。令p為戰(zhàn)場(chǎng)目標(biāo)域（X，Y）柵格處的探測(cè)率。p的更新公式為

式中：NM為MAV數(shù)量；pi（X，Y）為第i枚MAV對(duì)（X，Y）柵格處的探測(cè)率。

2.2MAV探測(cè)模型

MAV探測(cè)器探測(cè)范圍如圖4所示。

圖4　MAV探測(cè)示意圖Fig.4 Schematic diagram of MAV detection

2.3起飛點(diǎn)

設(shè)定每枚布撒器額定攜帶12枚MAV，介入方向相互垂直，起飛點(diǎn)坐標(biāo)如圖5所示，MAV經(jīng)彈機(jī)轉(zhuǎn)換和初始轉(zhuǎn)平飛過(guò)程后即開(kāi)始遂行區(qū)域覆蓋任務(wù)。第1枚布撒器拋灑的MAV起飛點(diǎn)以黑色圓點(diǎn)標(biāo)記，第2枚布撒器拋灑的MAV以黑色方塊標(biāo)記。

MAVi的起飛點(diǎn)坐標(biāo)（SXi，SYi，SZi）服從三維獨(dú)立正態(tài)分布，則

圖5　起飛點(diǎn)分布示意圖Fig.5 Schematic diagram of take-off point distribution

式中：CG為正向云發(fā)生器［13-16］；（ExXi，ExYi，ExZi）為MAVi的期望起飛點(diǎn)，

｛EnXi，EnYi，EnZi｝和｛HeXi，HeYi，HeZi｝為應(yīng)用逆向云發(fā)生器算法［13-16］和大量MAV初始轉(zhuǎn)平飛試驗(yàn)數(shù)據(jù)計(jì)算所得到的起飛點(diǎn)熵值和超熵，

2.4起飛航向

MAV初始轉(zhuǎn)平飛存在平飛航向的不確定性［9］。起飛航向?yàn)橐痪S均勻分布，即

2.5彈道

本文的彈道設(shè)計(jì)需要有兩點(diǎn)考慮：1）能夠覆蓋圓形作戰(zhàn)區(qū)域的邊界區(qū)域；2）便于MAV在離線(xiàn)條件下跟蹤控制。覆蓋作戰(zhàn)區(qū)域的邊界區(qū)域，MAV要在近邊界區(qū)域飛行一段距離；設(shè)計(jì)的彈道便于跟蹤控制，盡量不采用小的轉(zhuǎn)彎半徑，到達(dá)指定的目標(biāo)區(qū)域盡量采用大的轉(zhuǎn)彎半徑，同時(shí)分段航線(xiàn)在銜接點(diǎn)處相切。

圖6　MAV起飛航向示意圖Fig.6 Take-off direction of MAV

MAV的彈道可以參數(shù)化地表示為

式中：θSp為MAV的掃掠角。為了減少控制系統(tǒng)的負(fù)擔(dān)，本文借鑒dubins路徑設(shè)計(jì)方法，采用3段相切圓弧彈道拼接成MAV飛行彈道，令MAV在距地面100 m空域飛行。起飛點(diǎn)在（0，0，0），起飛航向?yàn)?，θSp為0.5πrad的彈道f（X，Y，100，0，0），如圖7所示。

圖7　彈道f（X，Y，100，0，0）示意圖Fig.7 Schematic diagram of trajectory f（X，Y，100，0，0）

對(duì)圖7中的彈道按照第1節(jié)和第2節(jié)中的方法賦予起飛點(diǎn)和起飛航向隨機(jī)變量，彈道形狀變?yōu)槿鐖D8所示的形式。

同時(shí)，由（2）式得

圖8　賦予隨機(jī)變量值后的彈道Fig.8 Trajectory after giving a random variab1e

3 目標(biāo)函數(shù)

NT為仿真次數(shù)為第j次仿真（X，Y）柵格點(diǎn)處的覆蓋率，battle為作戰(zhàn)域柵格點(diǎn)集是平均MAV數(shù)量，

由于J3是關(guān)于θSp的正比例函數(shù)，因此可以將J3進(jìn)一步簡(jiǎn)化為

實(shí)際作戰(zhàn)條件下，對(duì)目標(biāo)區(qū)域覆蓋率有明確的要求，因此將目標(biāo)轉(zhuǎn)化為約束條件，本文設(shè)定為覆蓋率不小于0.90.

因此目標(biāo)函數(shù)變?yōu)?/p>

4 仿真示例

4.1單次模擬飛行試驗(yàn)

起飛點(diǎn)按布撒器拋灑排列，引入起飛點(diǎn)、起飛航向隨機(jī)變量的條件下，某次模擬的24枚MAV的彈道如圖9所示。24枚MAV的探測(cè)軌跡疊加后的覆蓋效果如圖10所示。隨著MAV數(shù)量的增加，覆蓋率的變化如圖11所示。

圖9　單次模擬時(shí)MAV1～MAV24的彈道Fig.9 Trajectories of MAV1～MAV24in a simu1ation

圖10　24枚MAV彈道疊加后的覆蓋效果Fig.10 Coverage effect after superimposition of 24 MAV trajectories

圖11　單次模擬時(shí)隨著MAV數(shù)量增加的覆蓋率變化Fig.11 Vairiation of coverage rate when MAV number increases in the case of sing1e simu1ation

圖12　不同數(shù)量MAV條件下的覆蓋率均值隨模擬次數(shù)的變化Fig.12 The mean va1ue of coverage rate when number of simu1ations increases in the cases of different number of MAVs

圖13　不同數(shù)量MAV條件下的覆蓋率方差隨模擬次數(shù)的變化Fig.13 Variation of coverage variance when the number of simu1ations increases in the cases of different number of MAVs

4.2試驗(yàn)次數(shù)的確定

先排除起飛角度θS影響，將其設(shè)置為定值。不同數(shù)量MAV的條件下，覆蓋率均值和方差隨模擬次數(shù)的變化規(guī)律分別如圖12和圖13所示。

從圖13中可以看出，不同的MAV數(shù)量的作戰(zhàn)區(qū)域有效覆蓋率，在模擬次數(shù)5次以上時(shí)趨于穩(wěn)定。因此，對(duì)同一自變量進(jìn)行5次模擬以計(jì)算平均覆蓋率目標(biāo)。同時(shí)，對(duì)40次模擬的覆蓋率數(shù)據(jù)進(jìn)行正態(tài)分布假設(shè)檢驗(yàn)，確定覆蓋率服從正態(tài)分布。

處理多目標(biāo)問(wèn)題，可以將各個(gè)目標(biāo)進(jìn)行歸一化，并賦予不同的權(quán)重后，應(yīng)用單目標(biāo)智能優(yōu)化算法進(jìn)行求解；也可以應(yīng)用基于非劣解集的方法，對(duì)單目標(biāo)智能優(yōu)化算法進(jìn)行改進(jìn)，以適應(yīng)求解多目標(biāo)優(yōu)化問(wèn)題，這類(lèi)改進(jìn)后的算法常被稱(chēng)為多目標(biāo)優(yōu)化算法。當(dāng)前，求解多目標(biāo)問(wèn)題的智能算法獲得了很大的發(fā)展，且應(yīng)用領(lǐng)域也較為廣泛，如多目標(biāo)粒子群算法（MOPSO）［17］、多目標(biāo)遺傳算法（MOGA）［18］、多目標(biāo)蟻群算法（MOACO）等，也有較為年輕的多目標(biāo)智能優(yōu)化算法，如多目標(biāo)煙花算法（MOFWA）［19-20］等。

本文的彈道設(shè)計(jì)是一個(gè)離線(xiàn)問(wèn)題，對(duì)算法求解的實(shí)時(shí)性要求不高，只要能收斂到近似最優(yōu)解即可，因此采用MOGA對(duì)第3節(jié)中的目標(biāo)函數(shù)（17）式進(jìn)行優(yōu)化，優(yōu)化結(jié)果與分析見(jiàn)4.3節(jié)。

4.3結(jié)果與分析

優(yōu)化結(jié)果如圖14和表1所示。

由于多枚MAV在布撒器中為對(duì)稱(chēng)排列，MAV群組數(shù)量為偶數(shù)，故排除第4組、第5組和第6組解，第1組、第2組和第3組解可供決策者依據(jù)攻擊成本、攻擊耗時(shí)進(jìn)行決策。

從圖14和表1可以看出，MAV群組數(shù)量和覆蓋耗時(shí)之間在一定范圍內(nèi)存在此消彼漲的關(guān)系。MAV群組數(shù)量不大于12時(shí)，隨著MAV群組數(shù)量增加，覆蓋耗時(shí)明顯減少；但MAV群組數(shù)量大于12時(shí)，隨著MAV群組數(shù)量增加，覆蓋耗時(shí)不再顯著減少，MAV群組數(shù)量存在冗余。這是因?yàn)楸疚脑O(shè)計(jì)的彈道是基于2.5節(jié)中的兩點(diǎn)考慮，即便于MAV彈道跟蹤控制和覆蓋近邊界區(qū)域，所設(shè)計(jì)的彈道（見(jiàn)圖6）僅有1個(gè)自變量。

本文中MAV的探測(cè)半徑與目標(biāo)區(qū)域的半徑之比為1∶5，比值較大，對(duì)于較大目標(biāo)區(qū)域的覆蓋則可以考慮應(yīng)用文獻(xiàn)［17］中的方法進(jìn)行覆蓋；在起飛點(diǎn)和初始轉(zhuǎn)平飛方向存在不確定性的條件下，可以考慮重新設(shè)計(jì)彈道形狀，如阿基米德螺旋線(xiàn)或組合彈道等彈道形狀。而更大目標(biāo)區(qū)域的覆蓋則不應(yīng)是MAV群組的工作范疇。

圖14　多目標(biāo)遺傳算法優(yōu)化結(jié)果Fig.14 The optimized resu1ts of MOGA

表1　多目標(biāo)遺傳算法優(yōu)化結(jié)果Tab.1 The optimized resu1ts of MOGA

5 結(jié)論

本文提出的非協(xié)同方式多MAV區(qū)域覆蓋彈道設(shè)計(jì)方案，通過(guò)第1～4節(jié)的研究得出如下結(jié)論：

1）本文設(shè)計(jì)的多MAV區(qū)域覆蓋彈道方案可以提升MAV群對(duì)作戰(zhàn)目標(biāo)的偵查攻擊適應(yīng)能力，是在強(qiáng)電磁干擾、強(qiáng)電子誘騙條件下對(duì)協(xié)同攻擊的有力補(bǔ)充。

2）設(shè)計(jì)的彈道為3段圓弧拼接而成，其半徑均比MAV最小轉(zhuǎn)彎半徑大，降低了控制難度。

3）分析了覆蓋率數(shù)據(jù)的分布規(guī)律，確定了計(jì)算單組目標(biāo)函數(shù)的最小模擬次數(shù)，降低了計(jì)算成本。

4）設(shè)計(jì)的彈道兼顧了覆蓋率、MAV數(shù)量和航程3個(gè)目標(biāo)。

（References）

［2］沈林成，牛軼峰，朱華勇.多無(wú)人機(jī)自主協(xié)同控制理論與方法［M］.北京：國(guó)防工業(yè)出版社，2013. SHFN Lin-cheng，NIU Yi-feng，ZHU Hua-yong，et a1.Theories and methods of autonomous cooperative contro1 for mu1tip1e UAVs ［M］.Beijing：Nationa1 Defense Industry Press，2013.（in Chinese）

［2］Gu X Q，Shen L C，Chen J，et a1.A virtua1 motion camouf1age approach for cooperative trajectory p1anning of mu1tip1e UCAVs ［J］.Mathematica1 Prob1ems in Fngineering，2014（2）：1-15.

［3］王玥，馮順山，曹紅松，等.基于遺傳算法的多枚巡飛子彈探測(cè)階段協(xié)同彈道規(guī)劃［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2006（5）：80-82. WANG Yue，F(xiàn)FNG Shun-shan，CAO Hong-song，et a1.Coordinated trajectory p1anning of 1oitering sub-munition on detective stage based on genetic a1gorithm［J］.Tactica1 Missi1e Techno1ogy，2006（5）：80-82.（in Chinese）

［4］吳蔚楠.巡飛彈協(xié)同攻擊任務(wù)規(guī)劃問(wèn)題建模與方法研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2013. WU Wei-nan.Research on mode1ing and method for cooperative combat task p1anning of mu1tip1e cruise f1ight missi1es［D］.Harbin：Harbin Institute of Techno1ogy，2013.（in Chinese）

［5］李杰，游寧，李兵，等.云彈藥系統(tǒng)［J］.兵工學(xué)報(bào)，2015，36（2）：250-254. LI Jie，YOU Ning，LI Bing，et a1.C1oud ammunition system［J］. Acta Armamentarii，2015，36（2）：250-254.（in Chinese）

［6］Li B，Li J，Huang K W.Mode1ing and f1ocking consensus ana1ysis for 1arge-sca1e UAV swarms［J］.Mathematica1 Prob1ems in Fngineering，2013（6）：1-9.

［7］Li B，Li J，He G L，et a1.Research on cooperative combat for integrated reconnaissance-attack-BDA of group LAVs［J］.Mathematica1 Prob1ems in Fngineering，2014（1）：1-6.

［8］Kress M，Baggesen A，Gofer F.Probabi1ity mode1ing of autonomous unmanned combat aeria1 vehic1es（UCAVs）［J］.Mi1itary Operations Research，2006，11（4）：5-24.

［9］劉暢.傘降式小型飛行器初始轉(zhuǎn)平飛控制技術(shù)研究［D］.北京：北京理工大學(xué)，2014. LIU Chang.Technica1 research of initia1 state to 1eve1 f1ight contro1 for parachute-type sma11 aircraft［D］.Beijing：Beijing Institute of Techno1ogy，2014.（in Chinese）

［10］苗昊春，楊栓虎.智能化彈藥［M］.北京：國(guó)防工業(yè)出版社，2014. MIAO Hao-chun，YANG Shuan-hu.Inte11igent1ized ammunition ［M］.Beijing：Nationa1 Defense Industry Press，2014.（in Chinese）

［11］關(guān)震宇，楊東曉，李杰，等.基于Dubins路徑的無(wú)人機(jī)避障規(guī)劃算法［J］.北京理工大學(xué)學(xué)報(bào)，2014，34（6）：571-575. GUAN Zhen-yu，YANG Dong-xiao，LI Jie，et a1.Obstac1e avoidance p1anning a1gorithm for UAV based on Dubins path［J］. Transactions of Beijing Institute of Techno1ogy，2014，34（6）：571-575.（in Chinese）

［12］康崇祿.蒙特卡羅方法理論和應(yīng)用［M］.北京：科學(xué)出版社，2014. TANG Chong-1u.Theories and methods of mente-car1o［M］.Beijing：Science Press，2014.（in Chinese）

［13］李德毅，杜鹢.不確定性人工智能［M］.第2版.北京：國(guó)防工業(yè)出版社，2014. LI De-yi，DU Yi.Artificia1 inte11igence with uncertainty［M］. 2nd ed.Beijing：Nationa1 Defense Industry Press，2014.（in Chinese）

［14］苗奪謙，王國(guó)胤，姚一豫，等.云模型與粒計(jì)算［M］.北京：科學(xué)出版社，2012. MIAO DUO-qian，WANG GUO-yin，YAO Yi-yu，et a1.C1oud mode1 and granu1ar computing［M］.Beijing：Science Press，2012.（in Chinese）

［15］都業(yè)宏，連素杰，趙學(xué)剛.基于云模型和SWRL預(yù)測(cè)兵器試驗(yàn)故障的研究［J］.科學(xué)技術(shù)與工程，2012，20（5）：1175-1180. DU Ye-hong，LIAN Jie-su，ZHAO Xue-gang.Research on c1oud mode1 and swr1-based in weapon test fau1t prediction［J］. Science Techno1ogy and Fngineering，2012，20（5）：1175-1180.（in Chinese）

［16］都業(yè)宏，郁浩，李國(guó)富，等.基于云模型的射擊精度合格率及可靠性評(píng)估方法［J］.彈道學(xué)報(bào)，2013，25（3）：44-49. DU Ye-hong，YU Hao，LI Guo-fu，et a1.Fva1uation method of shooting accuracy rate and re1iabi1ity based on c1oud mode1［J］. Journa1 of Ba11istics，2013，25（3）：44-49.（in Chinese）

［17］Du Y H，Cui H，Li B，et a1.Reasearch on regiona1 coverage with LAVs on MOPSOA［J］.Journa1 of Information&Computiona1 Science，2015，12（6）：2081-2092.

［18］都業(yè)宏，郁浩，趙靜，等.武器系統(tǒng)預(yù)防性維修間隔期的多目標(biāo)決策研究［J］.兵工學(xué)報(bào)，2015，36（6）：1898-1905. DU Ye-hong，YU Hao，ZHAO Jing，et a1.Mu1ti-objective decision of preventive maintenance interva1s of weapon systems［J］. Acta Armamentarii，2015，36（6）：1898-1905.（in Chinese）

［19］Zheng Y J，Song Q，Chen S Y.Mu1tiobjective fireworks optimization for variab1e-rate ferti1ization in oi1 crop production［J］. App1ied Soft Computing，2013，13（11）：4253-4263.

［20］譚營(yíng).煙花算法引論［M］.北京：科學(xué)出版社，2015. TAN Ying.Introduction of fireworks a1gorithm［M］.Beijing：Science Press，2015.（in Chinese）

Modeling and Design of Trajectory of Multi-MAVs for Coveraging Multi-targets in Non-cooperative Manner

DU Ye-hong1，2，YU Hao1，ZHANG Jun1，LI Jie2，XU Xing-jian2
（1.Baicheng Ordnance Test Center of China，Baicheng 137001，Ji1in，China；2.Schoo1 of Mechatronica1 Fngineering，Beijing Institute of Techno1ogy，Beijing 100081，China）

To improve the adaptabi1ities of reconnaissance and attack of mu1ti-MAVs，the trajectory of mu1ti-MAVs for coveraging the mu1ti-targets in non-cooperative manner is mode1ed.The take-off point coordinates of MAVs are simu1ated using c1oud mode1，and the taking-off direction is simu1ated by Monte Car1o method.The target zone mode1，MAV detection mode1，take-off point mode1，f1ight direction mode1 and trajectory mode1 are estab1ished.To improve the efficiency of the optimization design，the minimum number of simu1ations is determined based on the stabi1ity ana1ysis of the expected coverage.The mu1titarget prob1em is so1ved by MOGA，and the regiona1 coverage trajectory of mu1ti-MAVs in non-cooperative manner is determined.The prob1ems that can be so1ved by the schemes are ana1yzed according to the ratio of detecting radius to target zone radius.The designed trajectories reduce difficu1ty in contro1，and take coverage，the number of MAVs and f1ying range of MAV into account.

ordnance science and techno1ogy；non-cooperative manner；MAV；regiona1 coverage；trajectory design；mu1ti-target optimization

TJ012.3；V249.1

A

1000-1093（2016）05-0936-09

10.3969/j.issn.1000-1093.2016.05.023

2015-07-16

總裝備部預(yù)先研究項(xiàng)目（2014年）

都業(yè)宏（1984—），男，工程師。F-mai1：yehong_du@126.com；郁浩（1971—），男，高級(jí)工程師。F-mai1：yuhao1i1ing@sina.com