Research on the Collision Probability of Missile Group Launching

WANG Zheng-wei; LIANG Bo; XIE Peng

doi:10.3969/j.issn.1673-3819.2020.04.021

Command Control and Simulation >

2020 , Vol. 42 >Issue 4: 107 - 110

DOI: https://doi.org/10.3969/j.issn.1673-3819.2020.04.021

Engineering & Application

Research on the Collision Probability of Missile Group Launching

WANG Zheng-wei ,
LIANG Bo ,
XIE Peng

Expand

The 28th Research Institute of China Electronics Technology Group Corporation,Nanjing 210007, China

Received date: 2020-03-26

Revised date: 2020-05-08

Online published: 2022-05-10

Fold

Abstract

In order to solve the minimum launch-timing sequence problem of ballistic missile cluster launching, the dynamic model of ballistic missile flight and the error covariance matrix model of missile position are established and the collision probability in 3D-space has been simulated in different time intervals. The computing of minimum launch timing and accurate solution of trajectory intersection problem is introduced, which provides a theoretical method for the optimization of fire plan of cluster launching.

Key words： group launch; trajectory intersection; collision probability; launching operation; mission scheduling

Cite this article

WANG Zheng-wei , LIANG Bo , XIE Peng . Research on the Collision Probability of Missile Group Launching[J]. Command Control and Simulation, 2020 , 42(4) : 107 -110 . DOI: 10.3969/j.issn.1673-3819.2020.04.021

鉴于美国的“发射即摧毁”能力,在弹道导弹大规模集群作战过程中,如图1所示,部队需要在最短的时间内完成导弹的发射,然后,快速进行撤收和隐蔽,使作战部队暴露时间最少,提升其作战生存概率。但是,大规模的导弹集中发射时并不能在同一时间进行点火发射,必须要通过保持一定的发射间隔来避免发射与飞行过程中弹间的相互气动影响和交叉碰撞问题^[1-2]。目前,在大规模火力分配过程中,通常只是将武器飞行弹道进行平面投影^[3],进行平面交叉分析,来确定最小发射时间间隔,这种方法忽略了导弹飞行高度影响,无法准确计算飞行器间的真实距离,会导致火力分配后的发射间隔大大增加。在导弹飞行过程中,还存在由于各种随机误差引起飞行理论偏差^[4⇓⇓-7],导弹飞行实际位置与理论弹道轨迹位置存在概率分布误差,理论弹道位置不能作为实际作战过程中弹道交叉的直接判据。本文利用导弹间各自的位置矢量和速度矢量,并引入了导弹位置误差协方差矩阵计算方法^[8⇓⇓-11],将问题转换为求解导弹位置概率分布下弹间的空间交叉碰撞概率问题,并对三维空间下弹间交叉碰撞概率的求解过程进行了推导,为弹道导弹集群发射时序确定问题提供了最优化解决方法。

显示原图|下载原图ZIP|生成PPT

图1 导弹集群发射中的弹道交叉

1 弹道交叉问题分析与建模

1.1 问题分析

弹道交叉是指导弹在飞行过程中距离过近时容易相互产生气动、电磁干扰甚至物理碰撞的情况,严重影响了导弹飞行安全。本次进行三维空间建模时遵循以下假设:

1)在导弹飞行主动段和再入段主要考虑气动、电磁干扰,导弹间相互干扰简化为弹间需保持一定的安全距离来避免干扰和碰撞。

2)导弹被动段飞行一般处于真空环境,主要考虑弹间的物理碰撞问题。

3)导弹飞行过程中,其位置误差是以理论弹道位置为中心的球状空间分布。

1.2 导弹运动模型

已知导弹的理论轨道方程^[9]可表示如下。

d X d t = F (t, X) d Y d t = F (t, Y) d Z d t = F (t, Z)

其中,t为导弹飞行时间,X、Y、Z分别为导弹在空间中的位置变量。

而在实际导弹的飞行过程中,由于天气、导弹本身控制的累积误差等方面的因素导致与理论计算得到的轨道方程有一定的偏差,为了更好地模拟实际飞行情况,引入随机误差变量

R i

,得到实际导弹的轨道方程为

d (X + R x (t)) d t = F (t, X + R x (t)) d (Y + R y (t)) d t = F (t, Y + R y (t)) d (Z + R z (t)) d t = F (t, Z + R z (t))

其中,

R x

、

R y

、

R z

为由于随机误差引起的轨道偏差,且轨道偏差为时间t的函数。在某一个时刻

t 0

,轨道偏差服从正态分布,即

R x (t 0) ~ N (X (t 0), σ x 2) R y (t 0) ~ N (Y (t 0), σ y 2) R z (t 0) ~ N (Z (t 0), σ z 2)

在任意时刻,且三个方向上的轨道偏差

R x

、

R y

、

R z

均为0的时候,导弹实际的轨道方程会退化为导弹理论的轨道方程。

1.3 弹间最小距离计算方法

设两条空间飞行弹道表达式^[12⇓-14]

r → 1 = x 1 (t 1) ι → + y 1 (t 1) J → + z 1 (t 1) k → r → 2 = x 2 (t 2) ι → + y 2 (t 2) J → + z 2 (t 2) k →

为了计算两条飞行弹道间的最短间距,建立变量

t → = [t 1, t 2]

和目标函数f(t)使得

f (t →) = ‖ r → 1 - r → 2 ‖ 2 = [x 1 (t 1) - x 2 (t 2)] 2 + [y 1 (t 1) - y 2 (t 2)] 2 + [z 1 (t 1) - z 2 (t 2)] 2

分别求导

∂ f ∂ t 1 = 2 [x 1 (t 1) - x 2 (t 2)] · x' 1 (t 1) + 2 [y 1 (t 1) - y 2 (t 2)] ·

y' 1 (t 1) + 2 [z 1 (t 1) - z 2 (t 2)] · z' 1 (t 1)

要求解目标函数的极小值,则需要满足:

∂ f ∂ t 1 = ∂ f ∂ t 2 = 0

[x 1 (t 1) - x 2 (t 2)] · x' 1 (t 1) + [y 1 (t 1) - y 2 (t 2)] ·

y' 1 (t 1) + [z 1 (t 1) - z 2 (t 2)] · z' 1 (t 1) = 0

[x 1 (t 1) - x 2 (t 2)] · x' 2 (t 2) - [y 1 (t 1) - y 2 (t 2)] ·

y' 2 (t 2) + [z 1 (t 1) - z 2 (t 2)] · z' 2 (t 2) = 0

则飞行弹道间最小距离求解转换为非线性方程组的求解问题,求解即可得到导弹飞行弹道表征的异面曲线之间的最小距离。

1.4 弹道交叉碰撞概率计算方法

通过对大量导弹飞行的理论与精确弹道数据间的误差进行统计,分析导弹飞行位置误差分布,进而得到导弹飞行的位置误差和状态矢量的协方差矩阵,并将上述协方差矩阵转换至同一坐标系中,可以得出弹间碰撞概率为三维高斯密度分布概率函数在总包络球中的积分。然后,对三维高斯密度分布函数在相对速度方向进行数值积分,将碰撞概率的计算公式简化成二维高斯密度分布概率函数,并求解其在总包络球投影圆域中的积分。

假设已知两枚导弹在各自发射坐标系中的位置误差协方差矩阵分别记为

C S

与

C R

。根据在地心坐标系中两枚导弹位置矢量和速度矢量,进行坐标换算,可以求解发射坐标系到地心坐标系的转置矩阵。这里将由发射坐标系到地心坐标系的转置矩阵表示为P,由地心坐标系到发射坐标系的转置矩阵表示为Q,如下式所示。

C S 1 = P C s P - 1 C R l = Q C R Q - 1

两导弹在速度坐标系中位置误差协方差矩阵均是对角矩阵。

C s = σ s 2 x 00 0 σ s 2 y 0 00 σ s 2 z C R = σ R 2 x 00 0 σ R 2 y 0 00 σ R 2 z

上述两协方差矩阵是不相关的,相加即可得到相对位置误差的协方差矩阵。

C c = C S 1 + C R 1

相对位置矢量和相对速度矢量的表达式为:

X c = r S - r R V c = v s - v R

定义地心坐标系到相遇坐标系的转换矩阵为U,在相遇坐标系中,相对位置矢量与协方差矩阵的转换公式如下所示^[8,11-12]。

X = U X e C = U C e U - 1

在完成相关系数的转换后,即可用三维高斯分布概率函数来表示两个相遇物体间的不确定性关系。

f (X) = 1 (2 π) 3 / 2 | C | 1 / 2 e x p - 12 X' C - 1 X

式中,|C|是协方差矩阵的行列式的值。

将上述方程沿着其相对运动方向进行数值积分,即可消掉与速度平行方向上的一维,完成三维概率分布函数向二维的简化。

已知相遇坐标系中,导弹的位置矢量为

X = x y z

概率函数密度的表达式为

f (X) = 1 (2 π) 3 / 2 | C | 1 / 2 e x p - 12 (c 11 x 2 + c 22 y 2 + c 33 z 2 + 2 c 12 x y + 2 c 13 x z + 2 c 23 y z)

随后我们对z方向积分,将三维积分化为二维;此时将x,y看作已知量,对后半部分做进一步简化,令

h x, y = ∫ - ∞ + ∞ e x p - 12 c 33 z 2 + 2 c 13 x z + 2 c 23 y z d z

将z分离,配方得到

h (x, y) = e x p (c 13 x + c 23 y) 2 2 c 33 ·

∫ - ∞ + ∞ e x p - c 33 2 z + c 13 x + c 23 y c 33 2 d z

f (x, y) = 1 2 3 / 2 π | C | 1 / 2 a 1 / 2 · e x p [- e x 2 - f y 2 - g x y]

其中,

a = c 33 2 e = c 11 2 - c 132 2 c 33 f = c 22 2 - c 232 2 c 33 g = c 12 - c 13 c 23 c 33

。

直接对函数f(x,y)在投影区域内进行数值积分,即可求出弹道交叉碰撞概率,如下式所示。

P = ∬ f (x, y) d x d y

由于式中含有交叉项,为简化计算,进行坐标旋转可得

f (x, y) = 1 2 3 / 2 π | C | 1 / 2 a 1 / 2 · e x p [- α x 2 - β y 2]

其中,

α = e + f - g 2 + (f - e) 2 2 β = e + f + g 2 + (f - e) 2 2

此时,经坐标系旋转后得到的坐标系表示为O-XY,其圆心O位于导弹的质心,X轴与Y轴分别对应于坐标系旋转之前相遇平面坐标系的X轴与Y轴。积分区域则为总包络球投影区域经过坐标旋转后的区域。

假设积分圆域中心为

(μ x, μ y)

μ x 2 + μ y 2 = μ r 2

,则

μ r 2

为导弹间间距。在坐标系转换过程中,只有导弹位置进行了投影变化,进行积分的区域形状没有变化,所以,积分区域的圆形半径R不变,仍为总球的半径。

所以其碰撞概率计算公式为

P = ∫ ∫ (x - μ x) 2 + (y - μ y) 2 ≤ R 2 1 2 3 / 2 π | C | 1 / 2 a 1 / 2 · e x p [- α x 2 - β y 2]

对O-XY进行坐标平移,让积分区域中心移动到坐标系原点,则碰撞概率求解公式可表示为

P = ∫ ∫ x 2 + y 2 ≤ R 2 1 2 3 / 2 π | C | 1 / 2 a 1 / 2 · e x p [- α (x + μ x) 2 - β (y + μ y) 2]

2 碰撞概率计算流程

碰撞概率计算流程如图2所示。

显示原图|下载原图ZIP|生成PPT

图2 碰撞概率计算流程

3 结束语

本文针对弹道导弹弹道交叉碰撞问题,建立了空间概率分布模型,推导求解了弹道导弹在三维空间中的空间碰撞概率,为确定弹道导弹集群发射时最小时间确定提供了计算方法,为发射阵地和发射时序的最优化确定提供了决策依据。

References

Publishing order | Descend order by publishing year | Descend order by cited within

[1]	张家铭, 刘忠. 考虑飞行初段弹道交叉的多导弹协同发射时序规划方法[J]. 系统工程与电子技术, 2016, 39(2): 316-321.

[2]	孟伦, 赵冠光. 巡航导弹飞行初段弹道交叉近似计算方法研究[J]. 装备制造技术, 2010(12): 61-62,75.

[3]	杨其, 李军良, 江光德. 弹道平面交叉算法[J]. 火力与指挥控制, 2009, 34(4): 54-56.

[4]	李泽秀, 朱昱, 鲜勇, 等. 弹道导弹飞行时间偏差修正方法[J]. 现代防御技术, 2017, 45(3):54-60,103.

[5]	魏宗康, 高荣荣, 周姣, 等. 基于椭圆概率误差的落点精度评定方法[J]. 导弹与航天运载技术, 2020(2):58-63.

[6]	杨云刚, 刘钧圣, 王琨, 等. 主要扰动因素对飞行弹道参数的影响分析[J]. 指挥控制与仿真, 2018, 40(3):127-131.

[7]	王彪, 吴楠, 孟凡坤, 等. 基于恩克法的雷达交接时刻预报与误差分析[J]. 指挥控制与仿真, 2019, 41(5):47-53.

[8]	张明选. 航天器碰撞概率的计算方法研究[D]. 哈尔滨:哈尔滨工业大学, 2010.

[9]	陈磊, 白显宗, 梁彦刚. 空间目标轨道数据应用[M]. 北京: 国防工业出版社, 2015.

[10]	何金铭, 刘江龙, 李瑞君. 一种计算空间异面高压线间最短距离的方法[J]. 电力勘测设计, 2015(z1):265-267,279.

[11]	马清华, 王明海. 弹道导弹常用坐标系及其转换[J]. 火力与指挥控制, 2004, 29(S1):89-91.

[12]	郝辉, 李雪瑞, 舒健生, 等. 导弹常用空间直角坐标系间转换方法[J]. 兵器装备工程学报, 2013, 34(2):18-20.

[13]	冯燕来, 王红杰. 基于外弹道修正理论的导弹落点预测精度评估方法[J]. 指挥信息系统与技术, 2017, 8(4):48-52.

[14]	曹静, 周茜, 基于机动检测的目标轨迹快速提取算法[J]. 指挥信息系统与技术, 2018, 9(3):70-74.

Options

Outlines

模态框（Modal）标题