導引頭空域搜索算法研究

栗高峰，董帥君

（上海機電工程研究所，上海 200233）

0 引言

從近代多次局部戰爭可發現，空襲兵器呈防區外發射、多批次、飽和攻擊態勢，要求防空導彈具有發射后不管及超視距攔截能力。隨著主動導引頭的研制成功，使防空導彈實現這一目標成為可能。為減小對地（艦）面武器系統的依賴，防空導彈武器系統常采取指令＋末端主／被動的復合制導模式，該制導模式的關鍵是中末制導交班的制導精度。受末端導引頭威力限制，系統的交班誤差與交班時導彈的飛行時間（距離）成正比。隨著防空導彈攔截目標空域縱深的擴展，會出現中末制導交班時因系統精度誤差的放大而無法滿足交班要求，這時常需要導引頭能實現快速空域搜索并鎖定目標。

當導引頭波速寬度和交班誤差已知時，為實現預定空域的覆蓋需要的最少搜索波位是系統必須求解的問題。誤差較小時，可粗略估算，誤差較大時，因目前尚無成熟的計算方法，難以估算搜索波位個數。目前對中末制導交班導引頭空域搜索算法研究較少，為此本文對導引頭空域搜索算法進行了研究。

1 空域搜索波位算法

1.1 空域搜索方式

1.1.1 圓掃描

圓掃描即按圓環形式進行掃描（即環掃描方式），導引頭在一周掃描過程中，波位中心移動軌跡呈圓形。表現為一層掃描波位的圓心在同一個圓上，各層掃描波位圓心對應的圓具有同一個圓心。如圖1所示：虛線為波位軌跡。

1.1.2 行掃描

行掃描即按行進行掃描，導引頭在一次掃描過程中，波位中心移動軌跡呈線性。根據交班時系統誤差散布往往呈圓概率散布的特點，用最少數量的圓覆蓋最大的面積，行掃描時采取由中間向兩端逐層遞減的方式。如圖2所示。

1.2 空域搜索波位算法

1.2.1 圓掃描波位

圖1 圓掃描Fig.1 Circular scanning

圖2 行掃描Fig.2 Line scanning

用圓掃描方式實現空域搜索時，采取兩層覆蓋法，即中心1個波位＋外圍一圈多個波位的方式，各波位兩兩間有一定的重疊區域。在波位寬度和波位數已知條件下，求所能覆蓋的最大掃描空域的半徑及各波位的位置信息。設導引頭方位和俯仰波束寬度一致，均為r；交班時系統精度誤差服從瑞利分布。圓掃描時，為保證波位間無縫交替，則要求相鄰掃描波位覆蓋空域有一定的重疊，設內外圈間的重疊區域為r－d，如圖3所示。

圖3中：待掃空域－大圓（圓1）方程為

中心波位－小圓（圓2）方程為

圓3方程為

圓4方程為

式中：θ為圓2的圓心至原點連線與X軸正方向夾角，且θ＝2π／N。此處：N為第二圈波位個數。圓掃描所需波位數n＝N＋1。

當搜索波位個數和波位半徑固定時，為獲得最大的搜索空域，且不能有搜索空隙，最理想的狀態是滿足外圈相鄰波位與中心圓有一個共同的交點，且外圈各相鄰波位的交點在中心圓上分布均勻。

a）求圓3、4交點坐標 （x1，y1），（x2，y2）

由式（3）、（4）可得

將式（5）代入（3），有

解方程（6）得

同時，保證該兩圓有兩個交點的條件

將點 （x1，y1）代入x2＋y2＝r2，有

b）計算能覆蓋的最大空域的半徑并確定各波位圓心坐標

第二層各掃描波位圓心坐標為

1.2.2 7波位行掃描波位

7波位行掃描如圖4所示。共用7個波位進行掃描，中間是3個波位的橫向掃描，由于其波位分布與7波位圓掃描相同，故其最大覆蓋空域半徑及波位位置信息求法同7波位圓掃描。

圖4 7波位行掃描Fig.4 Line scanning with 7beam positions

1.2.3 14波位行掃描波位

14波位行掃描如圖5所示。共用14個波位進行掃描，中間是4個波位的橫向掃描，其他層波位數由中間向兩端逐層遞減。波位采用圍繞圓1～4分別用六波位圓掃描的分布形式。

其中：第一行波位圓心坐標

第二行波位圓心坐標

圖5 14波位行掃描Fig.5 Line scanning with 14beam positions

第三行波位圓心坐標

第四行波位圓心坐標

第五行波位圓心坐標

式中：dmax為六圓形掃描的最大掃描半徑。

本文提出改進算法如下。

a）第一、五層的波位可向內移動一定位置，以減少溢出率并增加重疊面積；

b）去掉第一、五層的圓，并將第二、四層的圓向外移動一頂距離，減少掃描次數（即減少掃描時間），但會損失一定覆蓋率。

1.2.4 19波位行掃描波位

19波位行掃描如圖6所示。總共用19個波位進行掃描，中間是5個波位的橫向掃描，其他層波位數由中間向兩端逐層遞減。用圍繞波位圓0采用六波位圓掃描的分布得到波位1，2，3，4，5，6的位置，再分別圍繞該6個波位采用六波位圓掃描方式得到第三層波位的分布位置。

其中：波位1～6的圓心坐標

圖6 19波位行掃描Fig.6 Line scanning with 16beam positions

波位7的圓心坐標

波位8的圓心坐標

波位9的圓心坐

波位10的圓心坐標

波位11的圓心坐標

波位12的圓心坐標

波位13的圓心坐標

波位14的圓心坐標

波位15的圓心坐標

波位16的圓心坐標

波位17的圓心坐標

波位18的圓心坐標

式中：n為波位號。

1.2.5 6波位擴展行掃描方式波位

采用圓掃描的波位分布形式，第二層采用6個波位覆蓋后，針對第二層每一個波位，對外層每個波位依次6波位擴展；針對第三層每一個波位，對外層每個波位依次6波位擴展；其他層依次類推。得到如圖7所示的波位形式，最終使用的掃描方式為行掃描方式。

圖7 6波位擴展行掃描Fig.7 Extended line scanning with 6beam positions

2 算例

以一固定波束寬度r＝20，用上述各種波位個數搜索法求得能覆蓋的最大空域面積，結果見表1。

由表1可知：14波位行掃描方式下，覆蓋空域的面積大于14，15，16，17，18波位圓掃描的面積，且有更大的重疊距離；19波位行掃描方式較20波位以上的掃描方式有較大的覆蓋面積，且重疊距離更大。因此，在固定r＝20的條件下，若待掃描空域半徑范圍為54.641 0～57.298 9時，則優先使用14波位行掃描方式；若待掃描空域半徑大于57.298 9時，則優先使用19波位行掃描方式。這時，行掃描方式的覆蓋空域面積較相應圓掃描方式的更大，相鄰掃描波位間的重疊度更高。

表1 以半徑r＝20為例得到覆蓋的最大空域半徑結果Tab.1 Maxiuim coverage gained by r＝20

2.1 搜索波位個數計算

搜索波位個數計算方法的原則是在波速寬度和待掃描空域半徑固定條件下，得到搜索區域最大、重疊度高，且波束個數最小（即搜索時間最短）的方案。

假設待掃描空域半徑為R，固定波束寬度為r，令x＝R／r，y為要完全覆蓋半徑為R的待掃描空域而所需波位的最少個數，則

a）當x∈（1，1.618］時，y＝6，采用6波位圓覆蓋算法。

b）當x∈（1.618，2.732］時，

式中：［］表示向下取整數。

c）當x∈（2.732，2.866］時，y＝14，采用14波位行掃描算法。

d）當x∈（2.866，3.732］時，y＝19，采用19波位行掃描算法。

e）當x∈（3.732，＋∞）時，采用6波位擴展行掃描算法。

2.2 空域搜索波位算法

本文的搜索波位算法為：確定待掃描空域需要的波位個數及掃描方式，求出掃描波位的位置坐標及波位重疊度，調整相鄰掃描波位間的重疊度，獲得最佳空域掃描方式。算法流程如圖8所示。

圖8 算法流程Fig.8 Flowchart of algorithm

2.3 待掃描空域中心坐標確定

由本文算法算得的波位中心坐標是在待掃描空域中心坐標為中心的坐標系中的坐標。因導彈和目標均處于運動中，故待掃描空域中心的位置不斷變化，這時待掃描空域中心坐標的確定（即目標相對導彈的位置）顯得非常重要。目標和導彈均為勻速運動條件下的待掃描空域中心坐標確定方法如下：對如圖9所示的空情，設目標（t）和導彈（m）都為勻速直線運動，在時刻t＝0目標的起始位置和速度分別為 （xt0，yt0，zt0），（vtx，vty，vtz）；導彈的起始位置和速度分別為 （xm0，ym0，zm0），（vmx，vmy，vmz），則

a）目標在時刻t的當前位置為 （xt（t），yt（t），zt（t））。此處：xt（t）＝xt0＋vtxt；yt（t）＝yt0＋vtyt；zt（t）＝zt0＋vtzt。

b）導彈在時刻t的當前位置為 （xm（t），ym（t），zm（t））。此 處：xm（t）＝xm0＋vmxt；ym（t）＝ym0＋vmyt；zm（t）＝zm0＋vmzt。

目標相對導彈的坐標，即在以導彈為原點的坐標系中的掃描空域圓心坐標為

圖9 空情Fig.9 Air situation

3 結束語

本文對導引頭空域搜索算法進行了研究。給出了圓掃描和行掃描方式下的搜索波位個數計算方法。提出了待搜索空域半徑和搜索波位寬度已知條件下搜索方式的確定方法，以及相應的空域搜索波位算法。考慮目標與導彈的相對運動，給出了待掃描空域中心坐標的求法。該方法在導引頭搜索方式的設計中，有一定的參考作用。