基于反預警的反攔截中段規(guī)避突防策略

周啟航，劉延芳，＊，齊乃明，閆俊豐，2

基于反預警的反攔截中段規(guī)避突防策略

周啟航1，劉延芳1，＊，齊乃明1，閆俊豐1，2

1．哈爾濱工業(yè)大學 航天學院，哈爾濱 150006
2．61906部隊，廊坊 065000

為了充分挖掘多脈沖機動突防的潛力，提高彈道導彈中段機動突防能力，提出一種基于反預警的反攔截中段規(guī)避突防策略。首先，分別以脈沖點火點及點火間隔為搜索節(jié)點與步長，視敵方攔截系統(tǒng)為威脅源并轉(zhuǎn)化為航跡規(guī)劃中的約束條件，將多脈沖彈道的設計問題轉(zhuǎn)化成無人機避障航跡的規(guī)劃問題。然后，在綜合考慮敵方探測系統(tǒng)延遲和攔截系統(tǒng)部署及性能的基礎上，選取最優(yōu)的脈沖點火參數(shù)。最后，研究了脈沖增量和探測系統(tǒng)延遲對多脈沖規(guī)劃的影響。仿真結(jié)果表明：所提出的突防策略通過多脈沖機動增加敵方預測誤差，使其不能收斂與準確預警，當敵方的預測誤差收斂且準確預警時，導彈卻處于其攔截范圍之外，無法進行攔截；能根據(jù)獲取的最新戰(zhàn)場信息，在線下60s內(nèi)設計出一條多脈沖突防彈道。此外，還能根據(jù)任務的需求，提高（降低）導彈的飛行高度、平均飛行速度，縮短（增加）20%～35%的飛行時間。

彈道導彈；多脈沖；威脅源；避障航跡規(guī)劃；彈道設計；機動突防

隨著導彈突防技術的發(fā)展與進步［1－2］，出現(xiàn)了跳躍突防［3］、滑翔彈突防［4］、變射面突防［5］以及多脈沖機動突防［6］等多種中段機動突防策略。多脈沖機動突防能快速地改變導彈的彈道與射程，大幅提高敵方防御系統(tǒng)的探測與攔截難度，有效彌補彈道導彈彈道容易被探測以致被攔截的缺點。相比其他突防策略，多脈沖機動突防對彈頭的氣動外形、再入角度和防護等要求低且工程上相對容易實現(xiàn)。鑒于此，此項技術受到了國內(nèi)外學者的廣泛關注，并取得了一系列的研究成果。

由于導彈對時間、突防和打擊精度的特殊要求，在其他航天器上已經(jīng)成熟應用的多脈沖變軌原理及方法對其并不適用，在此不再作過多的闡述［7－8］。導彈多脈沖機動突防的關鍵是在保證導彈能擊中目標的前提下選取合適的脈沖點火參數(shù)，以使其突防能力最大。文獻［9］提出了在彈道最高點通過單脈沖改變導彈射程與剩余飛行時間來進行突防的策略；文獻［10］則將其推廣到任意點火時刻，并基于遍歷法對點火參數(shù)進行優(yōu)化。但文獻［9－10］只考慮了單次脈沖作用的情況，且采用的遍歷法的計算量很大。文獻［11－12］提出了利用多脈沖機動增加敵方預測誤差的突防策略，但沒有考慮機動后如何保證導彈打擊精度的問題。文獻［13－15］提出了利用敵方防御系統(tǒng)延遲進行突防的策略，并從同面及異面變軌的角度進行了分析。文獻［16］研究了點火時刻對該突防策略的影響，分析了其與導彈射程、飛行時間之間的關系。文獻［17］提出了利用攔截器點火空檔進行智能規(guī)避的策略。但以上多脈沖機動突防策略還存在以下局限性：

1）導彈的突防策略可以分為反探測和反攔截兩類。以上突防策略都只是針對其中某一方面，沒有將兩者綜合起來考慮導彈的突防問題，無法充分挖掘多脈沖機動突防的潛力。

2）采用的彈道設計方法必須預先選取中間彈道的落地點，使彈道方程組封閉，但選取規(guī)則并沒有給出。

3）不能依據(jù)戰(zhàn)場上敵方防御系統(tǒng)的最新情報信息，快速且有針對性地設計突防彈道。導彈的突防潛力沒能完全發(fā)揮且打擊任務的靈活性不足。

針對以上問題，本文借鑒路徑規(guī)劃中節(jié)點搜索及避障規(guī)劃的思想［18］，結(jié)合基于防御系統(tǒng)延遲的突防策略，提出了一種基于反預警的反攔截中段規(guī)避突防策略。綜合考慮探測系統(tǒng)的延遲與攔截系統(tǒng)的部署、攔截能力，通過“欺騙＋躲避”的方式完成多脈沖機動突防。本文提出的突防策略主要有以下4個創(chuàng)新點：

1）研究導彈突防問題時，綜合考慮敵方探測系統(tǒng)的延遲及攔截系統(tǒng)的性能與部署。使導彈具備反預警與反攔截的突防能力，提高其生存概率。

2）將攔截陣地視為障礙物，根據(jù)最新情報信息更新其位置、大小及威脅強度，最后將其變成路徑規(guī)劃中的多約束。導彈不是無目的程序式機動，而是有針對性地采用規(guī)避策略進行突防，以進一步提高其突防能力。

3）將規(guī)避彈道的設計問題轉(zhuǎn)化為避障路徑的規(guī)劃問題，采用路徑規(guī)劃的方法設計彈道，不再需要預先指定中間彈道的落地點。該方法能根據(jù)最新的戰(zhàn)場信息快速地設計一條突防能力高的彈道，提高導彈作戰(zhàn)的靈活性。

4）將三維空間中彈道到攔截陣地的距離問題轉(zhuǎn)換成二維平面內(nèi)點到橢圓的距離問題，將問題予以簡化。

1 多脈沖點火模型

1．1 基本定義

本文將地球視為不轉(zhuǎn)動的勻質(zhì)球體，采用軌道要素描述彈道。導彈點火時認為其直接獲取一個固定幅值的速度增量ΔV，ΔV在彈道面內(nèi)的投影與導彈速度vi（i＝0，1，…，n）（n為點火次數(shù)）的夾角為σi，與彈道平面的夾角為γi，點火點到地心的距離（即點火點的極半徑）為ri。

防御系統(tǒng)延遲：導彈防御系統(tǒng)中的探測系統(tǒng)能夠準確確定彈道所需的最短時間稱為防御系統(tǒng)延遲τ［13］，文獻［19］中為4～6min。

基于防御系統(tǒng)延遲的突防策略：導彈飛行過程中，每隔時間Ti（兩次點火間隔間導彈的飛行時間）將機動變軌一次，且Ti≤τ。由于防御系統(tǒng)延遲，敵方無法在Ti時間內(nèi)準確確定導彈的彈道并進行有效的攔截，使導彈突防成功［15］。

如圖1（a）所示，定義如下坐標系：

參考坐標系（OXYZ）：以地心O為原點，地心與參考彈道的關機點的連線為Y軸，在參考彈道面內(nèi)過球心且垂直Y軸指向?qū)嶋H目標點的直線為X軸，按照右手定則確定Z軸。

地心慣性坐標系（OXiYiZi）：以地心O為坐標原點，地心與第i次點火點連線為Yi軸，在當前彈道平面內(nèi)過地心并且垂直Yi軸指向當前虛擬目標的直線為Xi軸，按照右手定則確定Zi軸；OX0Y0Z0的Y0軸為脈沖彈道的初始關機點與地心的連線。

參考彈道：借鑒文獻［9－10］，人為給定一組符合實際情況的關機點參數(shù)，其射程大于、小于或等于打擊目標所需的射程。通常為了充分利用多脈沖提高導彈射程的能力，選取射程小于打擊目標所需射程的彈道作為參考彈道。不失一般性，本文選取由關機點到真實目標點之間的能量最省彈道為參考彈道，如圖1（a）中虛線所示。其關機點的極半徑、速度大小、彈道傾角和方位角記為（rr，vr，θr，ζr），同理多脈沖彈道的關機點參數(shù)記為（r0，v0，θ0，ζ0）。圖1（a）中的實線部分為多脈沖彈道，1，2，…，n為脈沖點火點，。令二者在關機點處的極半徑和速度等參數(shù)相同、彈道傾角與方位角等參數(shù)不同，即

式中：Δθ0和Δη0分別為多脈沖彈道與參考彈道在關機點處的彈道傾角差值和方位角差值。

多脈沖彈道初次關機時，其落地點為與真實目標有較大差距的虛擬目標；直到最后一次點火結(jié)束，導彈打擊的目標才為真實目標。虛擬目標與真實目標之間的差距由Δθ0和Δη0決定?；谧兩涿娴乃枷雽Ζう?和Δη0進行優(yōu)化，使虛擬目標與真實目標差距最大［5］，從而縮短敵方預警反應時間，提高導彈的突防能力。

圖1（b）為子彈道示意圖。如圖1（b）中所示，相鄰兩次點火點（A、B分別為i次和第i＋1此點火點）之間的彈道為子彈道。導彈脈沖點火n次，則有n＋1段子彈道。導彈在i次點火后，以點火點A作為關機點，在坐標系OXiYiZi中按照開普勒軌道定律飛行。飛行到B點時，進行第i＋1次點火，此時點火點的極半徑為ri＋1，該段彈道的射程角為Δfi＋1，飛行時間為Ti＋1，點火前后彈道面的夾角為Δηi＋1（即點火前后方位角差值）。從坐標系OXiYiZi到坐標系OXi＋1Yi＋1Zi＋1的轉(zhuǎn)換矩陣［20］為

同理，可求得參考坐標系OXYZ到坐標系OX0Y0Z0的轉(zhuǎn)換矩陣L0（Δ f0，Δη0），Δf0＝0°。

1．2 第i＋1次脈沖點火

本文設計的多脈沖彈道是由位于n＋1個共原點坐標系里的子彈道組成。在每個坐標系中將點火后的彈道參數(shù)作為關機點參數(shù)計算下次點火前的彈道參數(shù)。圖2為導彈第i次點火后在坐標系OXiYiZi中的子彈道，其在XiOYi面內(nèi)飛行，沒有Zi方向上的速度與位移。直線CD表示點火點當?shù)氐乃骄€。令ri、vi和θi分別表示在該子彈道上關機點的極半徑、速度大小和彈道傾角，（x，y，z）為導彈的空間位置，vx、vy和vz為沿各方向的速度分量，各變量的上標“－”和“＋”分別表示點火前和點火后，下標i表示在坐標系OXiYiZi中（下文同），則導彈i次點火后的彈道參數(shù)為

根據(jù)式（3）及文獻［21］可得第i段彈道的軌道參數(shù)為

式中：μ 為地球引力系數(shù)；υk，i、vi、ei、pi、θi、ai、bi和ci分別為該軌道的能量參數(shù)、導彈速度大小、離心率、半通徑、彈道傾角、長半軸、短半軸和焦距。

計算i次點火時的偏近點角Ei0及對應的真近點角

第i次點火完成并飛行Ti＋1時間后，導彈偏近點角為真近點角為則有

第i＋1次點火前的軌道參數(shù)［21］為

導彈的位置與速度大小分別為

式中：x－i＋1、y－i＋1、z－i＋1、v－x，i＋1、v－y，i＋1和v－z，i＋1為導彈第i＋1次點火前在坐標系OXiYiZi中的彈道參數(shù)。

第i＋1次點火時，脈沖速度增量為ΔVi＋1（大小為ΔVi＋1），方向角為σi＋1和γi＋1，則有

由式（10）可得

坐標系OXiYiZi到坐標系OXi＋1Yi＋1Zi＋1的轉(zhuǎn)換矩陣為Li＋1（Δfi＋1，Δηi＋1），點火后導彈在坐標系OXi＋1Yi＋1Zi＋1中的位置與速度參數(shù)組成向量Si＋1。坐標系OXi＋1Yi＋1Zi＋1中的Yi＋1軸為點火點與地心的連線，Yi＋1軸垂直當?shù)厮矫婕碭i＋1軸與當?shù)厮矫嫫叫?，所以點火后的彈道傾角為

在坐標系 OXi＋1Yi＋1Zi＋1中將Si＋1和θi＋1作為關機點參數(shù)計算第i＋1次點火前的彈道參數(shù)。

參考坐標系OXYZ 到坐標系OXi＋1Yi＋1Zi＋1的轉(zhuǎn)換矩陣Ci＋1為

經(jīng)過n次點火，導彈在坐標系OXYZ中的落地點坐標為（xtf，ytf，ztf）。

2 多脈沖機動彈道突防策略分析

2．1 突防策略分析

圖3為節(jié)點航跡搜索示意圖，圖中從關機點到目標點之間有3條彈道，其中彈道1和彈道2為多脈沖彈道，彈道3為參考彈道。圖3中有個威脅源；oj（j＝1，2，…，） 為威脅源的中心；Dj為每個威脅源的攔截范圍；tmi（m＝1，2）為第m條彈道的第i次脈沖點火時刻；tmi＋τ為第m 條彈道在第i次點火后敵方探測系統(tǒng)能夠準確確定彈道的最早時刻；dij（j＝1，2，…，） 為第i段子彈道與第j個威脅源中心的最短距離。此外，定義ti為任意彈道的第i次點火時刻。

彈道1為傳統(tǒng)的無人機避障航跡。導彈沒有穿過威脅源，而是從威脅源的邊緣繞飛過去。最后從突防概率最高的方向進入威脅源，擊中目標，稱為機動方式1。

彈道2不同于傳統(tǒng)無人機避障航跡規(guī)劃，規(guī)劃中考慮了導彈防御系統(tǒng)延遲τ這個因素。在滿足一定條件的情況下，導彈可以穿過威脅源而不被攔截，稱為機動方式2。具體情況如下：

1）ti＋1－ti≤τ

當子彈道的飛行時長小于防御系統(tǒng)延遲時，導彈無法被敵方防御系統(tǒng)確定并被攔截，因此是安全的。

2）ti＋1－ti＞τ

當子彈道的飛行時長大于防御系統(tǒng)延遲時，導彈在ti時刻點火，從ti＋τ到ti＋1之間的彈道到各個威脅源的最短距離dij大于對應威脅源的攔截范圍Dj，那么這段子彈道是安全的。

當子彈道的飛行時長大于防御系統(tǒng)的延遲時，敵方防御系統(tǒng)有足夠的時間來確定并攔截導彈。假如導彈在ti時刻機動，彈道發(fā)生了改變，由于系統(tǒng)延遲，敵方探測系統(tǒng)能夠完全確定彈道的最早時刻是ti＋τ。假如ti＋τ時刻導彈在威脅源之外，那么導彈此時無法被攔截。假如從ti＋τ到ti＋1之間彈道處于所有威脅源之外，那么認為這段彈道是安全的。對于其他子彈道，可依此類推。以圖3中的彈道2為例，t21＋τ與t22之間的彈道到威脅源j的最短距離為d11，威脅源1的攔截范圍為D1（d11＞D1）。該段彈道與威脅源之間關系依次為：d12＞D2，…，d1j＞Dj，所以該段彈道是安全的。

通過以上兩種機動方式，即可實現(xiàn)導彈中段的智能規(guī)避突防。

2．2 彈道設計問題轉(zhuǎn)化

圖4 中，（x0，y0，z0）和 （xt，yt，zt）分 別 為 給定的起點與終點位置。圖4（a）表示無人機避障航跡規(guī)劃中一條航跡鏈表，每個節(jié)點記錄了航跡節(jié)點的空間坐標 （ux，i，uy，i，uz，i）（i＝0，1，…，n）及對應的狀態(tài)變量Bi；（ux，t，uy，t，uz，t，Bt）為最終狀態(tài)信息。Bi含有當前節(jié)點的以下信息：①節(jié)點滿足航跡約束條件的情況；②該節(jié)點與下一個節(jié)點之間的航跡段是否可行；當且僅當航跡鏈表中所有節(jié)點與航跡段都可行時，航跡才可行。所以航跡集合中存在可行航跡及不可行航跡。在此基礎上，利用評價函數(shù)對每一條航跡的進行評價，搜索得到最優(yōu)航跡。

圖4 （b）所展示的為多脈沖彈道Pj（j＝1，2，…）。將點火點作為規(guī)劃的節(jié)點，第1個節(jié)點為初始關機點，最后1個節(jié)點為真實的目標點。節(jié)點個數(shù)固定，節(jié)點與節(jié)點之間的航跡不再是直線段，而是滿足開普勒軌道方程的橢圓彈道。圖4（b）用點火間隔Ti代替圖4（a）中節(jié)點在空間中的位置。因為在已知Ti的情況下，通過開普勒方程可得到Ti時間后導彈的位置，即當前節(jié)點的空間坐標。在多脈沖彈道設計中，彈道鏈表的節(jié)點參數(shù)為 （Ti，σi，γi，Ni），當前節(jié)點的彈道參數(shù)Si、θi及Ti＋1共同決定了下個節(jié)點的空間位置。Ni為狀態(tài)變量（N0和Nt分別為初態(tài)和終態(tài)），記錄以下信息：①脈沖點滿足彈道約束條件的情況；②當前點火點到下次點火點之間的子彈道是否可行；③該段子彈道與威脅源的關系。當Pj中所有節(jié)點及子彈道都可行時，彈道Pj才是可行的。

通過上述轉(zhuǎn)換，可將多脈沖中段規(guī)避突防彈道的設計問題轉(zhuǎn)化為無人機避障航跡的規(guī)劃問題。在以有分析基礎上，基于評價函數(shù)對彈道進行評價，再利用搜索算法可找到適應值最高的彈道。

2．3 子彈道與威脅源距離計算

計算子彈道與威脅源的距離采用以下步驟：

步驟1 將坐標系OXYZ中的威脅源Dj轉(zhuǎn)換到子彈道所在坐標系OXiYiZi中，并投影到當前彈道平面內(nèi)，得到投影g。

步驟2 視第i段子彈道為標準橢圓的一段，將其與g轉(zhuǎn)換到標準橢圓所在的坐標系中。

步驟3 在標準橢圓坐標系中，利用橢圓參數(shù)方程計算g與子彈道的最短距離，得到威脅源j與第i段子彈道的距離。

通過上述方法，將三維空間中點與彈道的距離問題轉(zhuǎn)換成二維平面內(nèi)點到橢圓的距離問題。

如圖5所示，過地心O以第i段橢圓彈道的近地點指向遠地點的直線為Yc軸，在彈道平面內(nèi)以過地心且垂直Yc軸的直線為Xc軸，建立直角坐標系XcΩOYc。以第i段橢圓彈道所在的標準橢圓的中心i（i＝1，2，…，n＋1）為Ω原點建立焦點在Y0軸上的直角坐標系X0iY0并建立輔助圓。

根據(jù)式（7）的定義可知，ti＋1、ti＋τ時刻導彈的偏近點角分別為和，根據(jù)式（2）和式（8）可得到從XcOYc到XiOYi的轉(zhuǎn)換矩陣l （）為

在Ω標準橢圓參Ω數(shù)方程中，定義：ti＋τ時刻，導彈和Ωi的連線與ΩiX0的夾角為αi；ti＋1時刻，導彈與i的連線與iX0的夾角為βi。則有

根據(jù)式（4）給出Ω的橢圓參數(shù)，可知兩時刻之間的彈道在平面X0iY0中用參數(shù)方程表示為

1）ti＋1－ti≤τ

第i段子彈道與第j個威脅源的距離為

式中：k為威脅源的距離常數(shù)。

為了使優(yōu)化算法更快地收斂，根據(jù)威脅源攔截范圍的不同對k作如下定義：

2）ti＋Ω1－ti＞τ

假設j在坐標系OXYZ 中的坐標為），在坐標系 OXiYiZi中的坐標為），在 當Ω 前 彈 道 平 面 內(nèi) 的 投 影 為，則在 X0iY0面中的投影為

此時，子彈道到威脅源的距離為

3 彈道的搜索與評價

3．1 威脅源性能參數(shù)

本文選取SM－3作為威脅源分析。SM－3導彈包括Block 0、Block I和Block II系列。查閱文獻［17］可知現(xiàn)役SM－3的最大射程為600km，最大攔截高度為160km。未來部署的Block II2A的射程可達1 200km，攔截高度為500km。多個SM－3攔截陣地可以相互通信，聯(lián)合在一起。通過這種方式組成的“威脅網(wǎng)”的威脅程度大于導彈各自威脅程度之和。

3．2 彈道約束表達與分析

下文在討論約束時，若無單獨說明，均視為在當前子彈道所在的坐標系OXiYiZi中。

1）最小及最大航跡段長度

通過點火間隔Ti來控制子彈道的長短。

a）最小子彈道長度：導彈點火前需要一定的時間來調(diào)整彈體姿態(tài)及發(fā)動機噴嘴方向，假設其所需的調(diào)整時間為Tmin，則可將約束表示為

最后一次點火結(jié)束后，導彈按照傳統(tǒng)彈道落地。如果打擊目標在威脅源內(nèi)，Tn＋1必須小于τ時才能突防成功，則有

b）最大子彈道長度：最大子彈道長度為導彈不點火，從關機點到落地點的飛行時間，設其為Tmax，約束表示為

2）最大拐彎角及最大爬升／下滑角

a）在彈道設計中，通過控制脈沖作用的角度來限制導彈最大（最小）的拐彎角及爬升角。為了使脈沖點火增加導彈的射程，要求點火時脈沖與導彈速度的夾角為銳角，即

b）為了迷惑敵方，由初始關機點決定的虛擬目標與實際目標的距離越遠越好。在擊中目標、滿足距離與時間要求的前提下，多脈沖彈道與參考彈道在關機點的彈道傾角與方位角差值應滿足：

3）最低飛行高度約束限制

本文中的最小高度約束為地球半徑R，根據(jù)式（3）可得

4）端點約束

彈道設計中的端點約束為初始與終端約束。初始約束與參考彈道的關機點參數(shù)一致，

為了能擊中目標，導彈需要滿足的終端約束為

為了保證導彈的攻擊效果，終端彈道參數(shù)還需要滿足速度和落角約束，即

式中：vtf和θtf分別為導彈的落點速度與彈道傾角；vf為落點速度；θf為落地點彈道傾角；下標min和max分別表示相應物理量的最小值與最大值。

5）威脅指數(shù)

第i段子彈道受到第j個威脅源的威脅指數(shù)FTAJ由Ti和dij決定，即

3．3 性能指標的選取

本文設計多脈沖彈道是為了提高彈道導彈的突防能力。因此，除了考慮飛行時間最優(yōu)外，還需考慮多脈沖彈道對敵方探測及防御系統(tǒng)的影響。

“虛擬落點”表示敵方探測系統(tǒng)依據(jù)點火前的彈道數(shù)據(jù)預測的導彈落地點。圖6為多脈沖對彈道落點的影響示意圖。如圖所示，1，2，…，n為脈沖點火點，n＋1為最終落地點；2′，…，i′，n′為“虛擬落點”；R1，R2，…，Rn為相鄰“虛擬落點”之間的距離。Ri的存在使得敵方在導彈點火后需要重新確定其落地點位置，其值越大對敵方的誤導性也越大即導彈突防能力更高；所以在優(yōu)化過程中，也將Ri作為優(yōu)化指標之一。

考慮到以上因素，本文采用如下優(yōu)化函數(shù)：

式中：Tt為總的飛行時間；w1，w2，m1，…，mn為權值，且都大于零；Rref為彈道落點差值的參考值，Tref為參考飛行時間，均用來對指標函數(shù)作無量綱處理。在實際彈道設計中，權值w1，w2，m1，…，mn依據(jù)敵方防御及探測系統(tǒng)的部署情況與能力的強弱來選取。

3．4 基于遺傳算法的彈道設計

為了提高計算效率，采用實數(shù)編碼的方式編寫染色體｛Δθ，Δη，T1，…，Tn，γ1，…，γn，σ1，…，σn｝。

假設一共點火n次且有個威脅源，則每條染色體含有變量個數(shù)為3n＋2，每一條染色體Pj代表一條彈道。

優(yōu)化中需要滿足的不等式約束Ωi為

不考慮彈道終端的速度與角度約束，只考慮其位置約束。因此優(yōu)化中需要滿足的等式約束i：

基于以上分析，航跡評價函數(shù)為

式中：U為懲罰系數(shù)。

4 仿真實驗算例

4．1 仿真參數(shù)

采用文獻［15］中的仿真例子，假設敵方防御系統(tǒng)的延遲時間為4min，導彈需要4次點火，打擊距離關機點4 666km外的目標。仿真中的參考彈道及攔截陣地的參數(shù)如表1和表2所示，表3仿真硬件環(huán)境參量。

表1 仿真參數(shù)Table 1 Parameters of simulation

表2 威脅源參數(shù)Table 2 Parameters of threats

表3 硬件環(huán)境Table 3 Hardware environment

4．2 仿真分析

分析過程中，導彈點火時獲取的速度增量大小ΔV 的取值范圍為100～500m／s，每隔10m／s取一個值進行仿真分析。對于每個ΔV，成功優(yōu)化30次后終止，取JGA值最小的一組彈道。選取ΔV＝250m／s、300m／s的例子對機動方式1進行仿真分析，優(yōu)化結(jié)果見表4與表5。選取ΔV＝360m／s、410m／s例子對機動方式2進行仿真分析，其優(yōu)化結(jié)果如表6與表7所示。表中：Tcost為優(yōu)化所需的時間；ΔVx，ΔVy和ΔVz分別表示速度增量在3個方向上的分量；i表示第幾次點火。

圖7和圖8分別為采用兩種機動方式的多脈沖彈道與參考彈道在三維彈道、飛行高度h、飛行速度大小v 3個方面的對比（均在參考坐標系OXYZ中，下文同）。圖7中，5條彈道從相同的關機點出發(fā)，經(jīng)過若干時間的飛行，擊中相同的目標。參考彈道不采取任何機動突防措施，將在威脅源1或3被敵方攔截。除參考彈道外，其余4條彈道都是由差異較大的5段異面子彈道組成。

表4 機動方式1的仿真結(jié)果（V＝250m／s）Table 4 Simulation results of the first maneuver mode（V＝250m／s）

表5 機動方式1的仿真結(jié)果（V＝300m／s）Table 5 Simulation results of the first maneuver mode（V＝300m／s）

表6 機動方式2的仿真結(jié)果（V＝360m／s）Table 6 Simulation results of the second maneuver mode（V＝360m／s）

表7 機動方式2的仿真結(jié)果（V＝410m／s）Table 7 Simulation results of the second maneuver mode（V＝410m／s）

依據(jù)表4和表5中的數(shù)據(jù)，結(jié)合圖7分析有如下結(jié)論：

1）ΔV＝250m／s、300m／s的彈道采用第1種機動方式，其與傳統(tǒng)無人機避障航跡類似，飛越威脅源，最后擊中目標，全程都處于威脅源外，敵方攔截系統(tǒng)無法進行攔截。

2）ΔV＝360m／s、410m／s的彈道采用第2種機動方式，彈道穿越威脅源。在敵方探測系統(tǒng)確定導彈彈道之前飛出威脅源，使敵方攔截系統(tǒng)無法攔截。

3）ΔR的范圍為229～343km，說明敵方探測系統(tǒng)必須重新采集導彈點火后的數(shù)據(jù)并確定其新的落地點。所以，直到導彈最后一次點火結(jié)束，敵方才能確定其真實落點即敵方預警時間被大大縮短。

4）由表4和表5可知，設計及優(yōu)化彈道所需的時間不超過60s。同時在蒙特卡羅打靶試驗中，設計彈道所需時間也不超過60s。這說明本文提出的中段規(guī)避突防策略可以根據(jù)瞬息變化的戰(zhàn)場情況，迅速地做出有針對性的反應，在很短時間內(nèi)得到一條線下規(guī)劃彈道，從而進一步提高導彈的突防能力。此外，導彈可以隨時根據(jù)新的目標信息，在距發(fā)射前很短的時間內(nèi)規(guī)劃出最優(yōu)彈道，提高了導彈打擊任務的靈活性。

5）所有的多脈沖彈道最后都擊中了目標，說明本文提出的中段突防策略能保證導彈的打擊精度。

以上結(jié)論說明采用兩種機動方式得到的突防彈道都可以使導彈成功規(guī)避敵方攔截系統(tǒng)，這也說明本文提出的中段規(guī)避突防策略是可行且有效的。

依據(jù)表4和表5中的數(shù)據(jù)，結(jié)合圖8分析有如下結(jié)論：

1）在飛行高度、時間方面機動方式1大于參考彈道，參考彈道大于機動方式2。根據(jù)二體軌道理論，導彈要維持橢圓彈道，導彈的平均速度增加，會導致其飛行高度的降低，導彈的平均速度降低，其飛行高度會增加。

2）相比參考彈道：機動方式1的平均飛行速度降低，飛行時間延長了21．25%～36．67%；機動方式2平均飛行速度提高，飛行時間縮短了30．14%～37．24%。通過多脈沖機動可以調(diào)節(jié)導彈的飛行速度、時間、彈道高度，根據(jù)突防需求設計滿足要求的彈道。

以上仿真分析結(jié)果說明，本文提出的中段規(guī)避突防策略可以有效地改變彈道參數(shù)，使敵方探測系統(tǒng)的探測數(shù)據(jù)失真，誤導及迷惑敵方的探測系統(tǒng)，提高導彈的突防能力。

4．3 防御系統(tǒng)延遲對彈道規(guī)劃的影響

為了研究系統(tǒng)延遲τ對彈道規(guī)劃的影響，在其他參數(shù)不變的情況下，改變其大小，觀察多脈沖突防彈道的變化情況。

仿真中，τ＝120～300s，間隔時間為5s，ΔV＝250m／s，其他仿真條件與4．1節(jié)一致。以τ＝180s和τ＝230s的情況進行舉例分析，圖9為突防彈道、導彈飛行高度及速度隨防御系統(tǒng)延遲的變化而變化得關系，圖10為穿越彈道的總飛行時間隨系統(tǒng)延遲的變化關系。1）采用本文提出突防策略可以得到兩種彈道：一種是從威脅源攔截范圍上限飛過，為飛越彈道；另一種是從威脅源中穿過或從其攔截下限飛過，表現(xiàn)形式為穿越威脅源，為穿越彈道。

2）τ＝180s時，只有飛越威脅源的彈道，無穿越彈道。τ＝230s時，飛越及穿越威脅源的彈道同時存在，且都滿足導彈的突防要求。但是兩種彈道的飛行高度及飛行速度差別很大。

3）根據(jù)圖9可知，飛越彈道的飛行高度高，飛行速度慢，飛行時間長；穿越彈道的飛行高度低，飛行速度快，飛行時間短。

4）當τ≥210s時，飛越及穿越彈道會同時存在，當τ＜210s時，只有飛越彈道。具體變化關系如圖10所示，隨著系統(tǒng)延遲的增加，導彈的飛行時間被縮短。所以，可以根據(jù)時間、距離及打擊任務的要求增加約束，對彈道進行取舍。

5 結(jié) 論

1）多脈沖彈道設計問題可以轉(zhuǎn)變成無人機避障航跡規(guī)劃問題，并使問題得以簡化。從而可以利用已經(jīng)比較成熟的無人機避障航跡規(guī)劃方法來解決多脈沖彈道設計問題。這樣能在保證導彈打擊精度前體下實現(xiàn)導彈的成功突防。

2）綜合考慮敵方探測及防御系統(tǒng)，使導彈獲得反預警及反攔截能力，確實能提高其突防能力。將敵方攔截系統(tǒng)視為障礙物，并處理成跡規(guī)劃中的多約束是可行的。根據(jù)實時獲取的敵方防御系統(tǒng)信息及時地改變規(guī)劃中的約束條件，得到最優(yōu)的線下規(guī)劃彈道，從而提高導彈打擊任務的靈活性。

3）本文提出的突防策略即可以獲得平均飛行速度大、飛行時間短且彈道高度低的彈道，也可以獲得平均飛行速度小、飛行時間長且彈道飛行高度高的彈道。根據(jù)具體突防需求選擇最合適的機動方式。

4）防御系統(tǒng)延遲會影響多脈沖彈道的規(guī)劃，要根據(jù)實際需要合理規(guī)劃彈道，使彈道的突防能力最大。

5）在其他條件固定時，如何選取點火次數(shù)、點火次數(shù)增減的規(guī)則、點火次數(shù)對彈道規(guī)劃的影響需要進一步進行深入的研究，也將是下一步研究的重點。

［1］ 鮮勇，李少朋，雷剛，等．彈道導彈中段機動突防技術研究綜述［J］．飛航導彈，2015，18（9）：43－46．XIAN Y，LI S P，LEI G，et al．Discussion of ballistic missile penetrating missile defense system［J］．Aerodynamic Missile Journal，2015，18（9）：43－46（in Chinese）．

［2］ 周啟航，張劉，霍明英，等．彈道導彈中段突防彈道設計與驗證［J］．光學精密工程，2015，36（9）：2645－2655．ZHOU Q H，ZHANG L，HUO M Y．et al．Ballistic missile midcourse penetration trajectory design and validation［J］．Optics and Precision Engineering，2015，36（9）：2645－2655（in Chinese）．

［3］ 鮮勇，李少朋，李振華，等．基于梯度粒子群算法的縱橫向機動跳躍彈道設計及優(yōu)化［J］．彈道學報，2015，27（3）：1－6．XIAN Y，LI S P，LI Z H，et al．Design and optimization of vertical wavy and crosswise maneuvering trajectory based on grads particle swarm algorithm［J］．Journal of Ballistics，2015，27（3）：1－6（in Chinese）．

［4］ 李廣華，張洪波，湯國建．高超聲速滑翔飛行器典型彈道特性分析［J］．宇航學報，2015，46（4）：397－403．LI G H，ZHANG H B，TANG G J．Typical trajectory characteristics of hypersonic glide vehicle［J］．Journal of Astronautics，2015，46（4）：397－403（in Chinese）．

［5］ 鮮勇，劉炳琪，王安民，等．變射面彈道橫向轉(zhuǎn)彎飛行程序的設計與優(yōu)化［J］．彈道學報，2014，26（4）：30－35．XIAN Y，LIU B Q，WANG A M，et al．Design and optimization of lateral turning flight programon changeable launching plane trajectory［J］．Journal of Ballistics，2014，26（4）：30－35（in Chinese）．

［6］ 明超，孫瑞勝，白宏陽，等．基于hp自適應偽譜法的多脈沖導彈彈道優(yōu)化設計［J］．固體火箭技術，2015，38（2）：151－155．MING C，SUN R S，BAI H Y，et al．Optimizing design of trajectory for multiple－pulse missiles based on hp－adaptive pseudo－spectral method［J］．Journal of Solid Rocket Technology，2015，38（2）：151－155（in Chinese）．

［7］ 鄧逸凡，李超兵，王志剛．一種基于軌道要素形式終端約束的航天器空間變軌迭代制導算法［J］．航空學報，2015，36（6）：1975－1982．DENG Y F，LI C B，WANG Z G．An iterative guidance algorithm using orbit element as terminal constraints for spacecraft orbit transfer［J］．Acta Aeronautica et Astronautica Sinica，2015，36（6）：1975－1982（in Chinese）．

［8］ 馮令兵，葛致磊，李憲強，等．彈道導彈被動段最優(yōu)脈沖變軌方法［J］．指揮控制與仿真，2014，36（2）：114－117．FENG L B，GE Z L，LI X Q，et al．A study on the ballistic missile’s optimal impulse orbit transfer in the passive ballistic curve［J］．Command Control Simulation，2014，36（2）：114－117（in Chinese）．

［9］ GILLIS M P I．Optimal mid－course modication of ballistic missile trajectories：Air Force Institute of Technology Re－port［R］．Washington，D．C．：Wright－Patterson Air Force Base，1975．

［10］ DRUCKMANN E，BENASHER J Z．Optimal in－flight trajectory modifications for ballistic missiles and free rockets［J］．Journal of Guidance，Control，and Dynamics，2012，35（2）：462－470．

［11］ GUO Y，YAO Y，WANG S，et al．Maneuver control strategies to maximize prediction errors in ballistic middle phase［J］．Journal of Guidance，Control，and Dynamics，2013，36（4）：1225－1234．

［12］ GUO Y，YAO Y，WANG S，et al．Maneuver strategy of evader considering detection system［C］／Proceedings of AIAA Guidance，Navigation，and Control Conference．Reston：AIAA，2011．

［13］ 張廣明，高普云，唐乾剛．導彈在被動段突防的脈沖式異面變軌方法［J］．宇航學報，2008，29（1）：89－94．ZHANG G M，GAO P Y，TANG Q G．The method of the impulse trajectory transfer in a different plane for the ballistic missile penetrating missile defense system in the passive ballistic curve［J］．Journal of Astronautics，2008，29（1）：89－94（in Chinese）．

［14］ 高普云，張廣明，唐乾剛，等．導彈在被動段突防的脈沖式軸 向 加 速 方 法 ［J］． 宇 航 學 報，2008，29（4）：1126－1130．GAO P Y，ZHANG G M，TANG Q G，et al．The method of increasing velocity along its original velocity direction by the impulse firing for the ballistic missile to penetrate the anti－missile system in the passive ballistic curve［J］．Journal of Astronautics，2008，29（4）：1126－1130（in Chinese）．

［15］ 高普云，張廣明，馮志剛，等．基于脈沖點火彈頭空間機動的被動段彈道設計方法和數(shù)值仿真［J］．宇航學報，2009，30（3）：884－889．GAO P Y，ZHANG G M，F(xiàn)ENG Z G，et al．The method of the trajectory design and numerical stimulation of the missile maneuvers in space based on the impulse firing in the passive curve［J］．Journal of Astronautics，2009，30（3）：884－889（in Chinese）．

［16］ 張斌，任軍學，王一白，等．導彈在被動段軸向加速變軌彈道設計與優(yōu)化［J］．北京航空航天大學學報，2010，21（6）：650－653．ZHANG B，REN J X，WANG Y B．Design and optimization on trajectory of increasing velocity along its original velocity direction by impulse firing for ballistic missile to penetrate anti－missile system in passive ballistic curve［J］．Journal of Beijing University of Aeronautics and Astronautics，2010，21（6）：650－653（in Chinese）．

［17］ 趙秀娜．機動彈頭的智能規(guī)避策略研究［D］．長沙：國防科學技術大學，2006：10－25．ZHAO X N．The intelligent evasive strategy of the ma－neuver missiles［D］．Changsha：National University of Defense Technology，2006：10－25（in Chinese）．

［18］ 丁明躍，鄭昌文，周成平，等．無人飛行器航跡規(guī)劃［M］．北京：電子工業(yè)出版社，2009：49－50．DING M Y，ZHENG C W，ZHOU C P，et al．Flight path planning for unmanned aerial vehicles［M］．Beijing：Electronics Industry Press，2009：49－50（in Chinese）．

［19］ 張廣明，高普云，唐乾剛，等．彈道導彈在被動段突防的平面脈沖式變軌方法［J］．彈道學報，2008，34（3）：83－86．ZHANG G M，GAO P Y，TANG Q G，et al．The method of the impulse trajectory transfer in a plane for the ballistic missile breaking through a defense line in the passive ballistic curve［J］．Journal of Ballistics，2008，34（3）：83－86（in Chinese）．

［20］ 賈沛然，陳克俊，何力．遠程火箭彈道學［M］．長沙：國防科技大學出版社，1993：173－199．JIA P R，CHEN K J，HE L．Long－range rocket ballistics［M］．Changsha：National University of Defense Technology Press，1993：173－199（in Chinese）．

［21］ 劉暾，趙鈞．空間飛行器動力學［M］．哈爾濱：哈爾濱工業(yè)大學出版社，2003：62－68．LIU D，ZHAO J．Aero spacecraft dynamics［M］．Harbin：Harbin Institute of Technology Press，2003：62－68 （in Chinese）．

［22］ 趙鈞．航天器軌道動力學［M］．哈爾濱：哈爾濱工業(yè)大學出版社，2011：42－87．ZHAO J．Orbital dynamics ofspacecraft［M］．Harbin：Harbin Institute of Technology Press，2011：42－87 （in Chinese）．

Anti－warning－based anti－interception avoiding penetration strategy in midcourse

ZHOU Qihang1，LIU Yanfang1，＊ ，QI Naiming1，YAN Junfeng1，2
1．School of Astronautics，Harbin Institute of Technology，Harbin 150006，China
2．61906Forces，Langfang 065000，China

An anti－warning－based anti－interception avoiding penetration strategy is proposed for fully exploiting potential of the multi－pulse maneuvering and improving the midcourse penetration ability of ballistic missiles．With the interception system being regarded as threats and transformed into the constraint conditions，the multi－pulse trajectory design is turned into an obstacle avoidance path planning for the unmanned aerial vehicle，wherein the pulse ignition point and interval are regarded as the search nodes and step，respectively．The optimal ignition parameters are selected based on a comprehensive consideration of the detection system delay and the performance and deployment of the interception system．Influence of pulse increment and system delay on multi－pulse trajectory is studied Simulation results show that the proposed strategy can damage the convergence of enemy’s prediction algorithm by multi－pulse maneuvering，and the accuracy of hostile warning results is thus reduced．When the enemy’s prediction error sufficiently converges and the hostile warning is accurate enough，the missile is already being out of the range of interception．According to the latest information in battlefield，an optimal multi－pulse trajectory can be designed offline in 60s，which can increase（decrease）the flight altitude and the average flight speed，and shorten（prolong）the flight time by 20%－35%．

ballistic missile；multi－pulse；threat；obstacle avoidance；ballistic design；maneuvering penetration

2015－12－03；Revised：2016－01－26；Accepted：2016－04－07；Published online：2016－04－13 15：35

URL：www．cnki．net／kcms／detail／11．1929．V．20160413．1535．004．html

Shanghai Aerospace Science and Technology Innovation Fund（SAST201312）

V448．131

A

1000－6893（2017）01－319922－13

http：／hkxb．buaa．edu．cn hkxb＠buaa．edu．cn

10．7527／S1000－6893．2016．0116

2015－12－03；退修日期：2016－01－26；錄用日期：2016－04－07；網(wǎng)絡出版時間：2016－04－13 15：35

www．cnki．net／kcms／detail／11．1929．V．20160413．1535．004．html

上海航天科技創(chuàng)新基金 （SAST201312）

＊通訊作者 ．E－mail：lyf04025121＠126．com

周啟航，劉延芳，齊乃明，等．基于反預警的反攔截中段規(guī)避突防策略［J］．航空學報，2017，38（1）：319922．ZHOU Q H，LIU Y F，QI N M，et al．Anti－warning－based anti－interception avoiding penetration strategy in midcourse［J］．Acta Aeronautica et Astronautica Sinica，2017，38（1）：319922．

（責任編輯：張玉）

＊Corresponding author．E－mail：lyf04025121＠126．com