甄昕，丁力軍，辛慶偉

（1.92941部隊；2.92493部隊，遼寧葫蘆島125001；3.海軍航空工程學(xué)院接改裝訓(xùn)練大隊，山東煙臺264001）

簡化導(dǎo)航算法在導(dǎo)彈飛行試驗安全分析中的應(yīng)用

甄昕1，丁力軍2，辛慶偉3

（1.92941部隊；2.92493部隊，遼寧葫蘆島125001；3.海軍航空工程學(xué)院接改裝訓(xùn)練大隊，山東煙臺264001）

安全分析是反艦導(dǎo)彈飛行試驗前必須開展的一項重要工作，安全分析以導(dǎo)彈彈道數(shù)據(jù)為基礎(chǔ)，而彈道仿真模型獲取不易?；诖?，提出了一種可用于反艦導(dǎo)彈飛行試驗安全分析并能取代彈道仿真的簡化導(dǎo)航算法，對其實現(xiàn)原理和實現(xiàn)方法進行了闡述。經(jīng)過實踐檢驗，此種算法計算簡單、效率高，其精度滿足導(dǎo)彈飛行試驗安全分析要求。

反艦導(dǎo)彈；飛行試驗；安全分析；威脅概率

靶場進行反艦導(dǎo)彈飛行試驗，首先要選擇試驗航區(qū)，并且要對試驗航區(qū)內(nèi)或周邊的石油鉆井平臺、大型島礁以及岸上目標可能造成的安全威脅作出評估，也就是要進行所謂的安全分析。

進行安全分析必須明確導(dǎo)彈飛行軌跡，也就是要確知導(dǎo)彈的相關(guān)彈道數(shù)據(jù)。目前，從解算導(dǎo)彈彈道角度而言，反艦導(dǎo)彈飛行試驗安全分析可分為2種：一種是從研制單位獲取彈道方程及相關(guān)氣動參數(shù)進行彈道仿真，來得到導(dǎo)彈飛行軌跡的安全分析方法。此種方法的優(yōu)點是精度高，缺點是彈道仿真模型難以獲取，即使獲取了仿真模型，編程計算也相當繁瑣，不是短時間能夠完成的[1]，而飛行試驗安全分析是開展飛行試驗必不可少的前提性工作，時間緊迫；另外一種方法即工程近似法，就是本文提出的基于簡化導(dǎo)航算法來模擬導(dǎo)彈飛行軌跡的工程近似法。實踐表明，此種方法易于理解，計算簡單，通用性強，效率高且精度滿足要求。

1 基本原理

導(dǎo)彈只有在其制導(dǎo)系統(tǒng)的導(dǎo)引下，才有可能飛向目標，一般反艦導(dǎo)彈在導(dǎo)航坐標系確立后，通過解算導(dǎo)航方程并進行相應(yīng)的誤差補償，得出導(dǎo)彈飛行所需的導(dǎo)航信息，如導(dǎo)彈經(jīng)度、緯度、目標方位角、導(dǎo)彈姿態(tài)角以及沿3個方向的加速度、速度等[2-3]。解算導(dǎo)航方程的過程就是在進行彈道計算，其模型復(fù)雜，計算繁瑣，輸入、輸出參數(shù)眾多，而在導(dǎo)彈飛行試驗安全分析中，一般只需要知道導(dǎo)彈飛行過程中的位置信息（任一點的經(jīng)緯度和方位角），且輸入?yún)?shù)很少，因而有必要尋找一種能夠較精確解算導(dǎo)彈飛行過程中位置信息的簡化導(dǎo)航算法，以提高試驗效率。

反艦導(dǎo)彈在地球表面附近運動，實際上是沿橢球面運動（地球近似橢球），由于反艦導(dǎo)彈射程較小，飛行高度低，其在空中的運動軌跡完全可以用參考橢球面上的大地線來近似[4-5]。

所謂大地線，就是橢球面上2點間的最短線，而且它是唯一的[6]。在理想狀態(tài)下，不進行扇面發(fā)射和航路規(guī)劃的反艦導(dǎo)彈在空中的飛行軌跡在橢球面上的投影，應(yīng)該和橢球面上連接導(dǎo)彈發(fā)射點和目標點的大地線近似重合[7]。因此，在對反艦導(dǎo)彈飛行試驗進行安全分析時，用參考橢球面上的大地線來近似代替其在空中的運動軌跡，就是簡化導(dǎo)航原理的出發(fā)點。

簡化的導(dǎo)航算法正是基于簡化導(dǎo)航原理，其核心就是貝塞爾大地問題正解公式[8]。

如圖1所示，已知P1點大地坐標(L1,B1)及其大地方位角α1，P1點至P2點的大地線長度為s，N為真北方向。計算P2點的大地坐標(L2,B2)和大地線在P2點的大地反方位角α2，稱為大地問題正解。

由球面三角學(xué)可知，球面三角公式都是用角度的三角函數(shù)表示的，球面上距離的大小與解算球面三角形的精度無關(guān)。又因為地球橢球的扁率很小，若將其上的諸元素（經(jīng)度、緯度、邊長及方位角）化為球面上對應(yīng)元素并均用角度表出時，其對應(yīng)的用角度表出的改正數(shù)一定是不大的，而且也和距離幾無關(guān)系。因此，貝塞爾公式解算大地問題的基本原理就是按一定的投影條件，建立橢球面上的元素和球面上對應(yīng)元素間的投影關(guān)系，通過解算球面三角形的方法而使問題得解。用此方法求解，其終點經(jīng)緯度及大地方位角可精確至0.003″，而且適用于長距離解算，具體計算公式詳見參考文獻[8]。

圖1 大地問題正解示意圖Fig.1 Positive solution schematic diagram ground question

在反艦導(dǎo)彈飛行試驗中，知道了導(dǎo)彈發(fā)射點O的經(jīng)緯度(L0,B0)，導(dǎo)彈射向α1，取解算步長h為100 m，運用貝塞爾大地問題正解公式不斷進行迭代計算，就可給出導(dǎo)彈飛行軌跡任一位置的經(jīng)緯度(L2,B2)和大地線反方位角α2（迭代時要將α2換算為大地線方位角），這樣就可以近似模擬導(dǎo)彈飛行過程中的位置信息。在實際應(yīng)用過程中，導(dǎo)彈射程數(shù)百千米，運用貝塞爾大地問題正解公式經(jīng)過迭代（h=100 m）解算得到的靶船經(jīng)緯度，與通過彈道仿真給出的靶船經(jīng)緯度相比，經(jīng)度相差不到2'、緯度相差不到1'（等同于在以研制單位給出的靶船位置為中心、半徑4 km的圓內(nèi)），在精度上能夠滿足導(dǎo)彈飛行試驗安全分析的要求。

2 算例分析

假設(shè)某型反艦導(dǎo)彈在O(L0,B0)發(fā)射，射向為α1，如圖2所示，已知導(dǎo)彈射程dk，動力航程為dl，平均速度為v，試計算導(dǎo)彈對陸地威脅概率。

圖2 導(dǎo)彈飛行試驗區(qū)域示意圖Fig.2 Regional schematic diagram of missile flying test

假設(shè)反艦導(dǎo)彈只有在導(dǎo)彈故障且安控失效的情況下，才可能對陸地構(gòu)成威脅。又假設(shè)導(dǎo)彈有2種主要故障模式即失控無偏墜毀和失控偏航墜毀，只有失控偏航墜毀才能對陸地構(gòu)成威脅。如果導(dǎo)彈在飛行彈道某點出現(xiàn)故障，失控偏航墜毀導(dǎo)彈的散布服從正態(tài)分布，散布范圍是以導(dǎo)彈射向為中心線的前半圓，不會對后方構(gòu)成威脅；從對陸地最大威脅情況考慮，導(dǎo)彈發(fā)生故障后以100%的概率飛完剩余航程（動力航程—已飛距離）。

2.1 建立導(dǎo)航坐標系

首先，將導(dǎo)彈彈道上某點的大地坐標(L,B,H)轉(zhuǎn)換為地心空間直角坐標(x,y,z)；再通過坐標平移和轉(zhuǎn)換，可將地心空間直角坐標(x,y,z)轉(zhuǎn)換為以此點為原點的導(dǎo)航坐標系(X,Y,Z)，此坐標系X軸與真北所成角度為α，方向指向目標一側(cè)，Y軸與水平面垂直、向上，Z軸在水平面內(nèi)與X、Y軸成右手法則[9]。它可以近似模擬反艦導(dǎo)彈真實飛行中的游動方位坐標系，其在發(fā)射點與導(dǎo)彈發(fā)射坐標系重合[10]。

2.2 威脅原理推導(dǎo)

在圖2上沿海岸線上取2點M(LM,BM)，N(LN,BN)，經(jīng)坐標平移和轉(zhuǎn)換后得到M(Z1,X1)，N(Z2,X2)，見圖3，使直線MN近似代替海岸線。取步長h=100 m（也可取其他值），則導(dǎo)彈發(fā)射后，彈道上間隔100 m的任一遞推點Ki的大地坐標和大地反方位角可由貝塞爾大地問題正解公式得到。再經(jīng)過坐標平移和轉(zhuǎn)換，可建立以Ki點為原點的導(dǎo)航坐標系。以此坐標系原點為圓心，r為半徑畫圓（r=動力航程-已飛距離），如果此圓與直線MN有2個交點，那么導(dǎo)彈可能對陸有威脅，否則沒有威脅。

圖3 對陸威脅概率計算示意圖Fig.3 Calculation schematic diagram of land threat probability

如圖3所示，經(jīng)推導(dǎo)，計算圓Ki與直線MN交點P1、P2及角度φ1、φ2的公式如下。

直線MN方程：

圓Ki方程：

解聯(lián)立方程（1）、（2），可得直線MN和圓Ki的交點P1(F,AF+B)、P2(G,AG+B)，且有坐標參數(shù)值如下：

在滿足D2＞4CE條件下，有式（3）~（6）成立。

當AF+B＞0與AG+B＞0同時成立時，有：

當AF+B＞0與AG+B≤0同時成立時，有：

當AG+B＞0與AF+B≤0同時成立時，有

當AG+B≤0與AF+B≤0同時成立時，Wi=0（導(dǎo)彈不會對陸地造成威脅）。

在D2≤4CE時，圓與直線MN有交點個數(shù)不超過1個，Wi=0（導(dǎo)彈不會對陸地造成威脅）。

運用彈道逐點積分原理（具體原理及若干計算公式詳見參考文獻[11]），可以得到導(dǎo)彈在Ki點故障落入到陸地的概率Pi。再由導(dǎo)彈在Ki點故障的概率Fi1和安控系統(tǒng)在Ki點失效概率Fi2，則Ki點對陸地的威脅概率為：Wi=PiFi1Fi2。導(dǎo)彈飛完全程對陸地的總威脅概率為

計算導(dǎo)彈對海上鉆井平臺和島礁的威脅概率比計算對陸威脅概率要簡單，只要給出鉆井平臺或島礁的大地坐標和它們的半徑數(shù)值，運用彈道逐點積分原理就可求得，其坐標系的建立和總威脅概率W的計算公式與對陸威脅概率計算一致。

2.3 編程計算

計算對陸威脅概率W的程序框圖見圖4。

圖4 計算導(dǎo)彈對陸威脅概率程序框圖Fig.4 Program diagram of calculation willile land threat probability

以上算例采用Matlab語言只需近百行源程序，且編程計算方便快捷，結(jié)果精度滿足試驗要求，如采用全彈道仿真方法（在獲得彈道仿真模型的前提下）進行對陸地威脅概率計算，編程計算就繁瑣得多，至少需近千行Matlab源程序[12]，花費大量的時間和人力，其效費比要低得多。假設(shè)發(fā)射點經(jīng)緯度為（121.532 7°,41.091 9°），射向為156°，步長100 m，射程290 000 m，，用以上方法編程計算可得目標經(jīng)緯度為（122.876 2°,38.708 9°）；用彈道仿真方法得到的目標經(jīng)緯度為（122.869 3°,38.692 1°）。經(jīng)計算2個目標點距離為4 000 m，滿足導(dǎo)彈飛行試驗安全分析要求。

3 結(jié)論

如反艦導(dǎo)彈有扇面發(fā)射或航路規(guī)劃，也可用貝塞爾大地問題正解公式求解，只是在導(dǎo)彈飛行彈道上選幾個已知經(jīng)緯度和方位角的轉(zhuǎn)折點，將原來一條大地線變?yōu)閹讞l首尾相連的大地線就可以了。采用以貝塞爾大地問題正解公式為核心的簡化導(dǎo)航算法，只要知道導(dǎo)彈發(fā)射點經(jīng)緯度、導(dǎo)彈射向、目標點經(jīng)緯度就可以進行導(dǎo)彈飛行位置信息模擬，并通過彈道逐點積分方法進行導(dǎo)彈威脅概率解算，可得到精確滿足要求安全分析結(jié)果，這是一種可行的工程近似方法。

[1]王精業(yè).仿真科學(xué)與技術(shù)原理[M].北京：電子工業(yè)出版社，2012：6-56. WANG JINGYE.Science and Technology Theory of Simulation[M].Beijing：Electronic Industry Press，2012：6-56.（in Chinese）

[2]楊軍.導(dǎo)彈控制原理[M].北京：國防工業(yè)出版社，2010：2-15 YANG JUN.Missile control theory[M].Beijing：The National Defense and Industry Press，2010：2-15.（in Chinese）

[3]潘榮霖.飛航導(dǎo)彈自動控制系統(tǒng)[M].北京：宇航出版社，1991：56-59. PAN RONGLIN.Missile auto control system[M].Beijing：Aero space Press，1991：56-59.（in Chinese）

[4]德米特里耶夫斯基.外彈道學(xué)[M].韓子鵬，譯.北京：國防工業(yè)出版社，2000：9-35. DMITRIYEVSKIY.The outer ballistic theory[M].Traslated by Han Zipeng.Beijing：The National Defense and Industry Press，2000：9-35.（in Chinese）

[5]孔祥云，郭際明，劉宗泉.大地測量學(xué)基礎(chǔ)[M].武漢：武漢大學(xué)出版社，2010：6-69. KONG XIANGYUN，GUO JIMING，LIU ZONGQUAN. The base theory of the ground measuring[M].Wuhan：Wuhan University Press，2010：6-69.（in Chinese）

[6]胡明城.現(xiàn)代大地測量學(xué)理論及其應(yīng)用[M].北京：測繪出版社，2003：12-23. HU MINGCHENG.The modern theory and application of the ground measuring[M].Beijing：Measuring Press，2003：12-23.（in Chinese）

[7]韓子鵬.彈箭外彈道學(xué)[M].北京：北京理工大學(xué)出版社，2008：4-41. HAN ZIPENG.Outer ballistic theory of missile and rocket [M].Beijing：Beijing Theory and Industry University Press，2008：4-41.（in Chinese）

[8]朱華統(tǒng)，黃繼文.橢球大地計算[M].北京：解放軍出版社，2001：79-84. ZHU HUATONG，HUANG JIWEN.The computation of the ground of the elliptic sphere[M].Beijing：PLA Press，2001：79-84.（in Chinese）

[9]羅建軍.組合導(dǎo)航原理及應(yīng)用[M].西安：西北工業(yè)大學(xué)出版社，2012：3-11. LUO JIANJUN.The Theory and Application of Composed Navigation[M].Xi'an：North West Industry University Press，2012：3-11.（in Chinese）

[10]孟秀云.導(dǎo)彈制導(dǎo)與控制系統(tǒng)原理[M].北京：北京理工大學(xué)出版社，2003：3-43. MENG XIUYUN.The missile guidance and control system theory[M].Beijing：Beijing Theory and Industry University Press，2003：3-43.（in Chinese）

[11]曲寶中，孫曉峰，李守秀.海軍戰(zhàn)術(shù)導(dǎo)彈試驗與鑒定[M].北京：國防工業(yè)出版社，2006：223-227. QU BAOZHONG，SUN XIAOFENG，LI SHOUXIU.Navy tactical missile test and evaluation[M].Beijing：The National Defence and Industry Press，2006：223-227.（in Chinese）

[12]DAVID ALLERTON.飛行仿真原理[M].劉興科，譯.北京：電子工業(yè)出版社，2012：8-26. DAVID ALLERTON.Flying simulation theory[M].Translated by LIU XINGKE.Beijing：Electronic Industry Press，2012：8-26.（in Chinese）

Application of Simplified Navigation Algorithm for the Safety Analyss During the Missile Flying Test

ZHEN Xin1,DING Li-jun2,XIN Qing-wei3
（1.The 92941stUnit of PLA;2.The 92493rdUnit of PLA,Huludao Liaoning 125001,China;3.Training Brigade of Equipment Acceptance and Modification,NAAU,Yantai Shandong 264001,China）

Safety analysis is an important task that must be finished before launching anti-ship missiles.The missile ballistic simulation data is the basic part for the safety analysis,but the missile ballistic simulation is difficult for us to complete. Because of this,a kind of simplified navigation algorithm for the safety analysis in the anti-ship missile flying test was put forword,and its principle and method was elaborated in this paper.This algorithm had some advantages such as simple computation,high efficiency.The result of its actual application was satisfied with the safety analysis in the anti-ship missile flying test.

anti-ship missile;flying test;safety analysis;threat probability

TJ760.6+28

A

1673-1522（2014）01-0005-04

10.7682/j.issn.1673-1522.2014.01.002

2013-10-10；

2013-11-25

甄昕（1970-），男，高工，大學(xué)。