面向臨近空間高機(jī)動(dòng)目標(biāo)的改進(jìn)預(yù)測(cè)命中點(diǎn)規(guī)劃方法

尹中杰，劉 凱

(大連理工大學(xué)航空航天學(xué)院，大連 116024)

0 引言

近年來(lái)，隨著高超聲速技術(shù)的發(fā)展，以美國(guó)為代表的軍事大國(guó)正在大力發(fā)展臨近空間高超聲速武器，并相繼完成了一系列的飛行試驗(yàn)[1-7]，如X-43A、X-51A、HTV-2等，同時(shí)還提出了武器化的HSSW、AHW等項(xiàng)目，以及SR-72高超聲速平臺(tái)項(xiàng)目。2018年10月，美國(guó)陸軍太空與導(dǎo)彈防御司令(SMDC)披露，將與美國(guó)空軍、海軍和導(dǎo)彈防御局合作研制一型高超聲速滑翔飛行器，其中美國(guó)空軍和陸軍分別計(jì)劃在2022年前列裝“高超聲速常規(guī)打擊武器”(HCSW)和“遠(yuǎn)程高超聲速武器”(LRHW)。

由于高超武器帶來(lái)的威脅，相應(yīng)的攔截方案研究受到了各國(guó)的重視，而預(yù)測(cè)命中點(diǎn)設(shè)計(jì)更是因?yàn)橛?jì)算時(shí)間長(zhǎng)受到廣泛關(guān)注。預(yù)測(cè)命中點(diǎn)是攻擊彈、攔截彈同時(shí)間到達(dá)的一點(diǎn)，預(yù)測(cè)命中規(guī)劃是指在攔截彈發(fā)射前，通過(guò)算法計(jì)算出合理的攔截點(diǎn)，并為攔截彈裝訂相應(yīng)的諸元，使得攔截彈準(zhǔn)時(shí)精確地飛抵預(yù)測(cè)命中點(diǎn)的工作。為了完成預(yù)測(cè)命中點(diǎn)的設(shè)計(jì)，國(guó)內(nèi)外學(xué)者普遍采用基于標(biāo)準(zhǔn)彈道[8-16]的預(yù)測(cè)命中點(diǎn)規(guī)劃。

基于基本彈道的預(yù)測(cè)命中點(diǎn)設(shè)計(jì)是指攔截方提前通過(guò)對(duì)所有可以覆蓋的區(qū)域，繪制標(biāo)準(zhǔn)彈道族，并通過(guò)調(diào)整攔截彈發(fā)射方位角形成等時(shí)間攔截曲面，并結(jié)合攻擊方彈道尋找時(shí)間相同、位置重合的相遇點(diǎn)。該計(jì)算方法可以針對(duì)不同類(lèi)型的目標(biāo)飛行器快速有效地尋找準(zhǔn)確的預(yù)測(cè)命中點(diǎn)。

由于需要攔截彈與目標(biāo)飛行器同一時(shí)間到達(dá)相同地點(diǎn)，所以解算攔截時(shí)間就是問(wèn)題的關(guān)鍵。李轅等[8]、張建祥等[9]、尤劉球[10]提出了以航程或攔截彈狀態(tài)參數(shù)作為輸入迭代相遇時(shí)間的方法。在此基礎(chǔ)上， 張榮升等[11]通過(guò)以斜距為輸入進(jìn)行時(shí)間迭代，獲得攔截彈與目標(biāo)斜距相同的時(shí)間點(diǎn)，從而獲得相應(yīng)的攔截彈彈道。但是高超武器機(jī)動(dòng)性強(qiáng)，彈道形式復(fù)雜，預(yù)測(cè)彈道上可能存在多個(gè)斜距相同的點(diǎn)，僅通過(guò)斜距無(wú)法描述預(yù)測(cè)命中點(diǎn)的信息，可能導(dǎo)致攔截失敗。而且由于算法原因，標(biāo)準(zhǔn)彈道族要逐一時(shí)間、逐一彈道與目標(biāo)進(jìn)行位置比對(duì)，耗時(shí)較多。

隨著攔截任務(wù)的明確與優(yōu)化算法的應(yīng)用，標(biāo)準(zhǔn)彈道族也可以根據(jù)所需的性能指標(biāo)，進(jìn)行有針對(duì)的設(shè)計(jì)。盛永智等[13]首先對(duì)攔截任務(wù)進(jìn)行分析，并針對(duì)不同戰(zhàn)場(chǎng)情況提出了相應(yīng)的性能指標(biāo)。在此基礎(chǔ)上，謝經(jīng)緯等[14]針對(duì)不同作戰(zhàn)需求設(shè)計(jì)了不同類(lèi)型的標(biāo)準(zhǔn)彈道族，大體分為考慮攔截時(shí)間需求的標(biāo)準(zhǔn)彈道族與考慮效費(fèi)比的標(biāo)準(zhǔn)彈道族，并根據(jù)防御方探測(cè)系統(tǒng)首次捕獲到目標(biāo)的時(shí)間選取針對(duì)相應(yīng)需求的標(biāo)準(zhǔn)彈道族，從而提高攔截效率。但是上述算法均只考慮了當(dāng)前時(shí)刻的預(yù)測(cè)命中點(diǎn)規(guī)劃問(wèn)題，并沒(méi)有對(duì)發(fā)射時(shí)間窗口計(jì)算問(wèn)題進(jìn)行深入研究，具有一定的局限性。

針對(duì)以上問(wèn)題，本文采用建立標(biāo)準(zhǔn)彈道族結(jié)合斜距單調(diào)化與時(shí)間搜索算法的方式解決因機(jī)動(dòng)造成的目標(biāo)飛行器與攔截陣地斜距非線性存在多個(gè)斜距相同點(diǎn)，導(dǎo)致獲得錯(cuò)誤預(yù)測(cè)命中點(diǎn)的問(wèn)題。利用攔截彈標(biāo)準(zhǔn)彈道族的物理特性增加判斷條件，降低計(jì)算量，解決運(yùn)算耗時(shí)長(zhǎng)的問(wèn)題，利用狀態(tài)轉(zhuǎn)換的思想設(shè)計(jì)發(fā)射時(shí)間窗口計(jì)算方法，從而快速準(zhǔn)確完成預(yù)測(cè)命中點(diǎn)設(shè)計(jì)與發(fā)射時(shí)間窗口計(jì)算工作。

1 攔截彈動(dòng)力學(xué)建模與問(wèn)題描述

1.1 質(zhì)點(diǎn)動(dòng)力學(xué)模型

模型考慮地球自轉(zhuǎn)、球形大地假設(shè)，在發(fā)射坐標(biāo)系下建立攔截彈質(zhì)點(diǎn)動(dòng)力學(xué)模型，根據(jù)攔截彈的特性完成坐標(biāo)系的轉(zhuǎn)化、受力分析并建立動(dòng)力學(xué)方程。

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

h=r-R

(9)

(10)

式中，x、y、z為發(fā)射系x、y、z軸位置，Vx、Vy、Vz為發(fā)射系x、y、z軸速度，θ、σ為發(fā)射系下彈道傾角、彈道偏角，D、L、Z、P分別為阻力、升力、側(cè)向力、推力，Mfv為發(fā)射系與速度坐標(biāo)系的轉(zhuǎn)換矩陣，h、r分別為高度、地心距，m為攔截彈質(zhì)量，g0為地表重力加速度，Re為地球半徑，ωe為地球自轉(zhuǎn)角速度，ωex、ωey、ωez為發(fā)射系下地球自轉(zhuǎn)角速度分量。

1.2 問(wèn)題描述

預(yù)測(cè)命中點(diǎn)是預(yù)先計(jì)算來(lái)襲目標(biāo)和攔截彈的飛行彈道, 確定出瞬時(shí)遭遇點(diǎn), 即預(yù)測(cè)命中點(diǎn), 當(dāng)預(yù)測(cè)命中點(diǎn)確定后，就導(dǎo)引攔截彈向其接近。

基于此，預(yù)測(cè)命中點(diǎn)規(guī)劃問(wèn)題可以描述為：結(jié)合預(yù)報(bào)彈道與標(biāo)準(zhǔn)彈道族，通過(guò)搜索算法尋找滿(mǎn)足時(shí)間約束、位置約束的攔截點(diǎn)。并為攔截彈裝訂相應(yīng)的諸元，確保其準(zhǔn)時(shí)精確到達(dá)攔截點(diǎn)，如圖1所示。

圖1 攔截示意圖Fig.1 Interception schematic

約束條件如下：

x(tf)=xtf
y(tf)=ytf
z(tf)=ztf

(11)

其中，tf為預(yù)測(cè)飛行時(shí)間，xtf、ytf、ztf為預(yù)測(cè)命中點(diǎn)在發(fā)射坐標(biāo)系下的位置。

考慮到攔截過(guò)程有延遲發(fā)射的需求，需要計(jì)算發(fā)射時(shí)間窗口。發(fā)射窗口是指在基于當(dāng)前標(biāo)準(zhǔn)彈道族，即使推遲發(fā)射，依舊存在預(yù)測(cè)命中點(diǎn)可行解的時(shí)間區(qū)間。

基于此，時(shí)間窗口計(jì)算問(wèn)題可以描述為：結(jié)合目標(biāo)預(yù)報(bào)彈道與標(biāo)準(zhǔn)彈道族，通過(guò)搜索算法獲得可發(fā)射攔截彈至預(yù)測(cè)命中點(diǎn)的延遲時(shí)間區(qū)間。

2 預(yù)測(cè)命中點(diǎn)規(guī)劃

2.1 現(xiàn)有方法的問(wèn)題

現(xiàn)有方法普遍基于攔截陣地與目標(biāo)之間的斜距進(jìn)行時(shí)間的迭代，通過(guò)假設(shè)目標(biāo)飛行器與發(fā)射陣地斜距是不斷縮小，選取時(shí)間t0，將目標(biāo)飛行器與發(fā)射陣地斜距ra與攔截彈斜距rd相對(duì)比，若ra>rd則t0增加，反之減小，通過(guò)迭代的方式獲得ra=rd的時(shí)間t1，從而獲得相應(yīng)的攔截彈彈道。而高超聲速飛行器會(huì)因?yàn)槎喾N因素，導(dǎo)致距離攔截陣地的斜距不單調(diào)，即彈目距離不斷縮小的假設(shè)不成立。此時(shí)基于上述原則的時(shí)間迭代方向可能存在錯(cuò)誤，這使得依靠現(xiàn)有方法解算預(yù)測(cè)命中點(diǎn)時(shí)，存在即使有可行解卻搜索不到的問(wèn)題。

當(dāng)目標(biāo)飛行器基于感知信息發(fā)現(xiàn)攔截陣地，會(huì)通過(guò)機(jī)動(dòng)遠(yuǎn)離攔截陣地，當(dāng)離開(kāi)攔截陣地一定距離后再次修正航向。這就導(dǎo)致預(yù)測(cè)彈道上存在若干斜距相同的點(diǎn)的情況，使得目標(biāo)飛行器與攔截陣地的斜距時(shí)而降低，時(shí)而增加。此時(shí)如果僅僅考慮斜距ra=rd，依舊采用ra>rd則t0增加，反之減小的迭代原則，可能會(huì)導(dǎo)致在時(shí)間迭代跳入斜距單調(diào)遞增的區(qū)間后，被提供了錯(cuò)誤收斂方向，在迭代達(dá)到上限后，輸出錯(cuò)誤的預(yù)測(cè)命中點(diǎn)，導(dǎo)致攔截失敗，如圖2所示。

圖2 斜距非單調(diào)示意圖Fig.2 Diagonal non-monotonic sketch

除此之外，目標(biāo)飛行器的打擊目標(biāo)不在攔截陣地的覆蓋范圍內(nèi)，目標(biāo)飛行器彈道穿過(guò)攔截陣地覆蓋范圍。此時(shí)必然存在兩個(gè)以上的斜距相同的點(diǎn)，致使斜距不單調(diào)，這同樣會(huì)導(dǎo)致預(yù)測(cè)命中點(diǎn)計(jì)算錯(cuò)誤。

基于以上問(wèn)題，這里提出了攔截時(shí)間區(qū)間的思想，將目標(biāo)與攔截陣地之間的斜距進(jìn)行單調(diào)化處理，確保在攔截區(qū)間內(nèi)任意特定斜距只存在唯一解，來(lái)解決斜距解算無(wú)法收斂到正確解的問(wèn)題。

2.2 標(biāo)準(zhǔn)彈道族設(shè)計(jì)

標(biāo)準(zhǔn)彈道設(shè)計(jì)采用目前工程規(guī)劃方法中典型的多級(jí)運(yùn)載火箭上升段指令生成規(guī)律進(jìn)行設(shè)計(jì)。該方法將火箭上升過(guò)程分為若干飛行階段，每段設(shè)定固定的程序角變化率指令形式，而這些指令由有限個(gè)參數(shù)決定，根據(jù)任務(wù)需求預(yù)先進(jìn)行優(yōu)化得到這些參數(shù)及其對(duì)應(yīng)程序指令，裝訂在助推器系統(tǒng)中，上升段飛行過(guò)程中利用姿態(tài)控制系統(tǒng)實(shí)現(xiàn)程序指令即可。為尋求彈道軌跡的快速生成，將整個(gè)上升段分為3段，包括垂直上升段(0-T11)、負(fù)攻角轉(zhuǎn)彎段(T11-T12)、重力轉(zhuǎn)彎段(T12-T13)。

(1)垂直上升段

(12)

式中，1/v0為推重比。

(2)負(fù)攻角轉(zhuǎn)彎段

在轉(zhuǎn)彎段中，T11到T12為負(fù)攻角轉(zhuǎn)彎段，如圖3所示。在之后的大動(dòng)壓段一般依靠重力的法向分量緩慢的轉(zhuǎn)彎，即重力轉(zhuǎn)彎段(T12～T13)。

圖3 攻角變化規(guī)律曲線Fig.3 Angle of attack variation with time

在負(fù)攻角轉(zhuǎn)彎段中，根據(jù)實(shí)際要求，攻角變化規(guī)律可以由經(jīng)驗(yàn)公式給出

α(t)=-αm·sin2f(t)

(13)

(3)重力轉(zhuǎn)彎段

此段內(nèi)程序攻角為0°，攔截彈全程飛行俯仰角近似公式為

(14)

當(dāng)確定了T11、T12、T13，所選取的攔截彈彈道控制變量就可以確定為助推段負(fù)攻角大小和發(fā)射方位角。

由此，以攔截陣地為原點(diǎn)(0，0，0)在最大負(fù)攻角區(qū)間αm∈[αmmin,αmmax]中，每隔Δα積分標(biāo)準(zhǔn)彈道，再由tmin開(kāi)始，將隔條攔截彈道中時(shí)間為tmin的點(diǎn)連接，構(gòu)成等時(shí)間線，等時(shí)間線遵守tmin至tmax每間隔Δt連接一條的原則。由此可以獲得完整的標(biāo)準(zhǔn)彈道族，如圖4所示。

圖4 標(biāo)準(zhǔn)彈道族Fig.4 Standard trajectories collection

2.3 目標(biāo)飛行器攔截時(shí)間段劃分

為解決不同時(shí)間、斜距相同的問(wèn)題，需要對(duì)目標(biāo)飛行器進(jìn)行攔截區(qū)域劃分，并分段討論。

首先，獲得斜距與時(shí)間的關(guān)系，確定一定時(shí)間內(nèi)的斜距極值點(diǎn)個(gè)數(shù)n與相對(duì)應(yīng)的時(shí)間tn，由此將攔截時(shí)間分為n+1段，時(shí)間段為t∈[ti-1,ti],i∈[0,n+1]，由此確保在規(guī)定時(shí)間段內(nèi)，斜距變化是單調(diào)的，不存在斜距相同的可能，如圖5所示。

圖5 目標(biāo)飛行器與攔截陣地斜距-時(shí)間圖Fig.5 Range-time diagram of target vehicle and interception position

2.4 基于搜索算法的預(yù)測(cè)命中點(diǎn)計(jì)算

通過(guò)對(duì)攔截時(shí)間的劃分可以確保在對(duì)應(yīng)的時(shí)間段內(nèi)，目標(biāo)飛行器與攔截陣地的斜距是單調(diào)遞增(遞減)的，可以通過(guò)t時(shí)刻，目標(biāo)飛行器的經(jīng)緯度與高度與攔截彈的經(jīng)緯度與高度是否相同來(lái)迭代預(yù)測(cè)命中點(diǎn)時(shí)間t與發(fā)射方位角af。但是由于攔截彈的發(fā)射方位角度af∈[0,360]，負(fù)攻角α∈[αmin,αmax]，產(chǎn)生的標(biāo)準(zhǔn)彈道族巨大，所以逐一時(shí)間、逐一標(biāo)準(zhǔn)彈道進(jìn)行對(duì)比是不現(xiàn)實(shí)的，需要利用標(biāo)準(zhǔn)彈道族的物理特性基于新的搜索算法進(jìn)行預(yù)測(cè)命中點(diǎn)搜索，提高運(yùn)算效率。

攔截彈的可控變量有3個(gè)，分別是攔截彈的最大負(fù)攻角、發(fā)射方位角、飛行時(shí)間。最大負(fù)攻角與飛行時(shí)間共同決定了攔截彈的飛行高度、飛行距離，發(fā)射方位角決定攔截彈的飛行方向。其中飛行方向只與發(fā)射方位角有關(guān)，但是剩余兩項(xiàng)存在耦合現(xiàn)象。

因此，搜索算法可以在標(biāo)準(zhǔn)彈道族中先篩選滿(mǎn)足水平航程的攔截時(shí)間區(qū)間[tmin，tmax]，再根據(jù)目標(biāo)飛行器的飛行高度在攔截時(shí)間區(qū)間[tmin，tmax]，及負(fù)攻角區(qū)間[αmin,αmax]中搜索滿(mǎn)足高度條件的飛行時(shí)間tf及負(fù)攻角α，最后利用tf時(shí)刻的目標(biāo)飛行器位置確定發(fā)射方位角af，確定攔截彈的飛行軌跡，從而裝訂發(fā)射諸元。

根據(jù)攔截條件可知，只有全部滿(mǎn)足高度、水平航程、飛行方向條件的攔截彈才具有攔截的能力，缺少任意條件均無(wú)法有效攔截。但是一次性完成3個(gè)條件的篩選，算法過(guò)于復(fù)雜，所以可以率先篩選出滿(mǎn)足水平航程的攔截時(shí)間區(qū)間作為攔截時(shí)間估計(jì)區(qū)域，并在后續(xù)步驟中進(jìn)一步篩選符合高度條件的時(shí)間點(diǎn)，確定最大負(fù)攻角，最后計(jì)算發(fā)射方位角。具體步驟如下：

(1)初步估計(jì)時(shí)間區(qū)間[tmin，tmax]

根據(jù)攔截彈的物理特性可知，相同時(shí)間內(nèi)，最大負(fù)攻角最大的攔截彈軌跡水平航程最大，則通過(guò)搜索最大負(fù)攻角最大的攔截彈道與目標(biāo)彈道中同一時(shí)刻水平航程誤差最小的點(diǎn)，就可以獲得tmin。

假設(shè)初始飛行時(shí)間為t1，并通過(guò)目標(biāo)飛行器預(yù)報(bào)信息計(jì)算出該時(shí)刻的發(fā)射方位角，同時(shí)獲得t1時(shí)刻的目標(biāo)飛行器與攔截陣地的水平距離La1、t1時(shí)刻最大負(fù)攻角最大的攔截彈水平距離Ld1，誤差為ΔL1；獲得攔截彈攔截窗口終端t3時(shí)刻的目標(biāo)飛行器與攔截陣地的水平距離La3、t3時(shí)刻最大負(fù)攻角最大的攔截彈水平距離Ld3，誤差為ΔL3；獲得中間時(shí)刻t2的目標(biāo)飛行器與攔截陣地的水平距離La2、t2時(shí)刻與最大負(fù)攻角最大的攔截彈水平距離Ld2，誤差為ΔL2。

這里分為4種情況，并設(shè)計(jì)搜索算法：

① ΔL1>0, ΔL3>0時(shí)，說(shuō)明該時(shí)間區(qū)域內(nèi)沒(méi)有攔截的可能，搜索結(jié)束。

② ΔL1<0, ΔL2<0時(shí)，說(shuō)明該時(shí)間區(qū)域內(nèi)沒(méi)有攔截的可能，搜索結(jié)束。

圖6 時(shí)間搜索示意圖Fig.6 Time search diagram

確定了tmin，即可利用tmin，縮小tmax的搜索范圍。同理，相同時(shí)間內(nèi)，最大負(fù)攻角最小的攔截彈軌跡水平航程最小，則通過(guò)搜索最大負(fù)攻角最小的攔截彈道與目標(biāo)彈道中同一時(shí)刻水平航程誤差最小的點(diǎn)，就可以獲得tmax。

假設(shè)初始飛行時(shí)間t1=tmin，并通過(guò)目標(biāo)飛行器預(yù)報(bào)信息計(jì)算出該時(shí)刻的發(fā)射方位角，同時(shí)獲得t1時(shí)刻的目標(biāo)飛行器與攔截陣地的水平距離La1、t1時(shí)刻最大發(fā)射傾角的攔截彈水平距離Ld1，兩者誤差為ΔL1；獲得攔截彈攔截時(shí)間區(qū)域終端t3時(shí)刻的目標(biāo)飛行器與攔截陣地的水平距離La3、t3時(shí)刻最大負(fù)攻角最小的攔截彈水平距離Ld3，兩者誤差為ΔL3；獲得中間時(shí)刻t2的目標(biāo)飛行器與攔截陣地的水平距離La2、t2時(shí)刻與最大負(fù)攻角最小的攔截彈水平距離Ld2，誤差為ΔL2。并再次利用二分法確定最大值tmax。

(2)標(biāo)準(zhǔn)彈道選取

由于前期的發(fā)射時(shí)間估計(jì)區(qū)域的確立，根據(jù)t∈[tmin,tmax]的目標(biāo)飛行器的位置信息確定t時(shí)刻相對(duì)應(yīng)的攔截彈發(fā)射方位角af∈[afmin,afmax]，根據(jù)目標(biāo)飛行器的高度，確定與目標(biāo)飛行器距離誤差最小的負(fù)攻角α，存儲(chǔ)t時(shí)刻ΔLt，最后選取t∈[tmin,tmax]中ΔLt最小的彈道作為攔截彈道。

具體流程圖如圖7、圖8所示。

圖7 篩選流程圖Fig.7 Screening flow

2.5 發(fā)射時(shí)間窗口計(jì)算

由于攔截任務(wù)的需要，攔截彈可能需要延時(shí)發(fā)射，所以有必要計(jì)算出可行的攔截彈延時(shí)發(fā)射時(shí)間區(qū)間以及相關(guān)諸元。

由于攻擊彈與攔截陣地的斜距可以轉(zhuǎn)化為距離攔截陣地的水平距離與垂直距離，同理攔截彈與攔截陣地的斜距也可以轉(zhuǎn)化為距離攔截陣地的水平距離與垂直距離。以水平距離為橫坐標(biāo)，垂直距離為縱坐標(biāo)，不同標(biāo)準(zhǔn)彈道與預(yù)報(bào)彈道的交點(diǎn)就是可能的預(yù)測(cè)命中點(diǎn)。由于攻擊彈到達(dá)可能的預(yù)測(cè)命中點(diǎn)的時(shí)間t4與攔截彈到達(dá)可能的預(yù)測(cè)命中點(diǎn)的時(shí)間t5是已知的，只要滿(mǎn)足t4>t5就認(rèn)為是可選的預(yù)測(cè)命中點(diǎn)。通過(guò)比較所有的可選預(yù)測(cè)命中點(diǎn)，得到最大的延遲時(shí)間Δt=max(t4-t5)，則獲得攔截時(shí)間窗口t∈[0,Δt]，并通過(guò)t4時(shí)刻的目標(biāo)位置計(jì)算出相對(duì)應(yīng)的發(fā)射方位角af。

3 仿真分析

3.1 仿真輸入

設(shè)目標(biāo)飛行器初始飛行速度Ma=10，初始為平衡滑翔狀態(tài)，60s～80s只進(jìn)行縱向機(jī)動(dòng)，80s～120s只進(jìn)行橫向機(jī)動(dòng)。攔截陣地經(jīng)緯度為(-3.9°，79.8°)，目標(biāo)飛行器發(fā)現(xiàn)點(diǎn)經(jīng)緯度為(-2.10°，80.58°)，滿(mǎn)足預(yù)報(bào)精度誤差的預(yù)報(bào)時(shí)間為120s，如圖9所示。

圖8 預(yù)測(cè)命中點(diǎn)及發(fā)射諸元解算流程圖Fig.8 Flow charts for predicted impact point and launch data calculation

圖9 目標(biāo)軌跡圖Fig.9 Trajectory diagram

3.2 預(yù)測(cè)命中點(diǎn)規(guī)劃仿真分析

由于目標(biāo)飛行器可能存在機(jī)動(dòng)，所以斜距與飛行時(shí)間不是簡(jiǎn)單的線性關(guān)系，需要對(duì)目標(biāo)預(yù)報(bào)彈道進(jìn)行單調(diào)化處理，具體關(guān)系如圖10所示。

圖10 目標(biāo)飛行器、攔截陣地斜距—時(shí)間圖Fig.10 Range-time diagram of target vehicle and interception position

由圖10可見(jiàn)，目標(biāo)飛行器與攔截陣地斜距先降低后增加，所以發(fā)射窗口被分為兩個(gè)區(qū)間，分別為t∈[0,90],t∈[90,120]。

在t∈[0,90]區(qū)間內(nèi)，通過(guò)迭代發(fā)現(xiàn)在89s時(shí)，發(fā)射方位角為-35.73°，負(fù)攻角為-10.26°，預(yù)測(cè)命中點(diǎn)為(-3.6°，79.7°)，此時(shí)采用標(biāo)稱(chēng)彈道的攔截彈與目標(biāo)的預(yù)報(bào)位置誤差為0.52km，可以通過(guò)后續(xù)的末制導(dǎo)環(huán)節(jié)進(jìn)行修正，認(rèn)為符合攔截條件，如圖11所示。

圖11 攔截示意圖Fig.11 Interception schematic

在t∈[90,120]區(qū)間內(nèi)，通過(guò)迭代沒(méi)有發(fā)現(xiàn)滿(mǎn)足當(dāng)前時(shí)刻發(fā)射的攔截點(diǎn)，所以符合命中條件的預(yù)測(cè)命中點(diǎn)只有一個(gè)，為(-3.6°，79.7°)，計(jì)算用時(shí)0.2s。

3.3 攔截時(shí)間窗口計(jì)算仿真分析

通過(guò)對(duì)發(fā)射窗口計(jì)算，發(fā)現(xiàn)最大延遲時(shí)間為11s，發(fā)射方位角為-82.19°，負(fù)攻角為-12.25°,預(yù)測(cè)命中點(diǎn)(-3.9°，79.4°)，此時(shí)采用標(biāo)稱(chēng)彈道的攔截彈與目標(biāo)的預(yù)報(bào)位置誤差為0.34km，同樣可以通過(guò)后續(xù)的末制導(dǎo)環(huán)節(jié)進(jìn)行修正，計(jì)算用時(shí)0.3s，如圖12所示。

圖12 最大延遲時(shí)間攔截示意圖Fig.12 Maximum delay time interception diagram

針對(duì)拋物線彈道延遲發(fā)射諸元如表1所示。

由表1可知，由雷達(dá)首次發(fā)現(xiàn)目標(biāo)的時(shí)刻至隨后的11s內(nèi)，隨著攔截高度的降低，最大負(fù)攻角幅值變大。同時(shí)所有選取的時(shí)間節(jié)點(diǎn)均有攔截彈道可以準(zhǔn)時(shí)飛抵預(yù)測(cè)命中點(diǎn)，可以認(rèn)為在0s～11s這個(gè)攔截窗口內(nèi)攔截窗口時(shí)間連續(xù)，所有的預(yù)測(cè)命中點(diǎn)均有標(biāo)準(zhǔn)彈道與其相對(duì)應(yīng)。

表1 延時(shí)發(fā)射諸元表

4 結(jié)論

本文針對(duì)臨近空間高機(jī)動(dòng)飛行器機(jī)動(dòng)能力強(qiáng)，所帶來(lái)的彈道中存在多個(gè)斜距相同的點(diǎn)，導(dǎo)致基于斜距迭代預(yù)測(cè)命中點(diǎn)的方法不適用的問(wèn)題，通過(guò)彈道劃分的方法，保證了每個(gè)搜索區(qū)間的斜距具有單調(diào)性。在此基礎(chǔ)上利用改進(jìn)的快速搜索算法，準(zhǔn)確有效地完成預(yù)測(cè)命中點(diǎn)規(guī)劃與發(fā)射時(shí)間窗口計(jì)算?？紤]高超聲速飛行器預(yù)報(bào)精度差，所帶來(lái)的更新后的預(yù)測(cè)命中點(diǎn)與原有預(yù)測(cè)命中點(diǎn)位置偏差較大的問(wèn)題，未來(lái)將開(kāi)展攔截彈道在線快速糾偏方法的研究，來(lái)彌補(bǔ)此方面的不足。