CN112464451B - 基于作战仿真系统的防空导弹武器命中概率修正方法 - Google Patents

基于作战仿真系统的防空导弹武器命中概率修正方法 Download PDF

Info

Publication number
CN112464451B
CN112464451B CN202011278680.5A CN202011278680A CN112464451B CN 112464451 B CN112464451 B CN 112464451B CN 202011278680 A CN202011278680 A CN 202011278680A CN 112464451 B CN112464451 B CN 112464451B
Authority
CN
China
Prior art keywords
target
aircraft
correction
hit probability
airplane
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011278680.5A
Other languages
English (en)
Other versions
CN112464451A (zh
Inventor
黄隽
刘方
李晓宝
刘玥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN202011278680.5A priority Critical patent/CN112464451B/zh
Publication of CN112464451A publication Critical patent/CN112464451A/zh
Application granted granted Critical
Publication of CN112464451B publication Critical patent/CN112464451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Abstract

本发明所设计的基于作战仿真系统的防空导弹武器命中概率修正方法,包括如下步骤:步骤1:计算抗击飞机的抗击距离修正后的命中概率;计算抗击飞机目标速度修正后的命中概率;计算抗击飞机目标机动性修正后的命中概率;步骤2:计算防空导弹武器抗击飞机目标或抗击制导目标修正后的命中概率;步骤3:计算抗击制导目标方位角修正的抗击制导目标命中概率。本发明基于作战仿真系统,采用数学解析模型和数值拟合相结合的方法实时在线修正防空导弹武器命中概率,提高效率,增强工程操作性。

Description

基于作战仿真系统的防空导弹武器命中概率修正方法
技术领域
本发明涉及导弹武器系统运用技术领域,具体地指一种基于作战仿真系统的防空导弹武器命中概率修正方法。
背景技术
防空导弹命中概率(hit probability)属于导弹作战效能评估指标。导弹的命中包括导弹的正确定位、破片及爆轰产物覆盖目标,破片及爆轰产物命中目标三个贯序发生的事件,目标运动及几何特征等特性,不同战斗部形成的威力环和威力圆的几何特征两个基本要素。命中概率的估计属于导弹试验鉴定与效能评估问题。
防空导弹命中概率建模关键技术一般分为解析法、模拟法、拟合法和Bayes方法四种基本方法类型,或者基本类型结合的方法。事件概率分析法是一种常用解析方法,首先将一项任务划分为若干环节,对各个环节都进行衡量并定出标准,最后通过评价各个环节的完成概率得出命中概率。防空导弹从发射到引爆,经历初制导段、中制导段、末制导段、截击目标段。
防空导弹命中概率解析模型通常表示为制导误差概率分布密度在散布圆内的积分。制导误差概率分布密度参数包括导弹随机落点脱靶量距离的期望和随机落点对平均落点分布距离的均方差。导弹随机落点脱靶量距离的期望和随机落点对平均落点分布距离的均方差可根据验前试验数据和现场试验数据利用Bayes方法获取。根据有无系统误差,制导误差概率分布密度分为瑞利分布和莱斯分布两类,莱斯分布可用常用近似表示,也可用正交多项式精确逼近。
空空导弹命中概率解析模型也可由一般中制导的捕获概率、末制导的追踪概率以及导引头发现目标距离计算。而捕获概率为制导误差的积分,追踪概率由拦截角(目标速度矢量方向与防空导弹速度矢量方向的夹角)分布密度函数积分。研究表明,空空导弹的攻击飞机目标进入角(目标速度矢量与空空导弹-飞机目标相对距离矢量之间的夹角,也叫目标攻击角)对目标命中概率的影响,迎头和尾追时命中概率最大,前后半球的命中概率较大,越靠近两侧离轴方向命中概率越小。目标速度越大、目标过载越大,命中概率较高的目标进入角区域就越小。
空空导弹命中概率解析模型还可以表示为干扰条件下目标探测概率、自控终点对目标的截获跟踪概率、中末制导交接班捕获概率、末制导精度的乘积。其中,中末制导交接班捕获概率表示为中制导末端导弹位置、速度大小、速度方向实际值/理论值比值的乘积,末制导精度表示为脱靶量、末速度、入射角的实际值/理论值比值的乘积。
模拟计算法需要构建和解算由质量方程、动力学方程、运动学方程、制导方程、控制方程组成的导弹模拟微分弹道方程组,进行拦截仿真或者半实物仿真,获得防空导弹的导弹模拟打靶试验数据,得到脱靶量范围,然后修正得到计算的理论命中概率,该方法作为理论研究固然可以探讨,但工程上只能由导弹研制部门实现,未见公开文献报道的可完全求解的6自由度导弹模拟全弹道方程组。
数值拟合方法,存在函数难确定、函数分段范围难把握的问题;拟合函数的形式多种多样。对于数据的一维拟合问题,可通过分段处理,将节点光滑的连接起来,可以做到拟合误差任意小。但是对于多维问题,若进行分段,就非常难于处理,因为多维问题分段处将会出现一个曲面或者一条曲线甚至会是一个更高维的情况。对于高维情况的出现,拟合分段就很不好处理了。
Bayes方法需要验前试验数据和现场试验数据,主要用于防空导弹命中概率的参数估计和鉴定,并不适用于作战仿真推演。
因此,研究防空导弹武器命中概率修正方法对于防空作战模拟具有重要意义。
发明内容
本发明的目的就是要提供一种基于作战仿真系统的防空导弹武器命中概率修正方法,本发明基于作战仿真系统,采用数学解析模型和数值拟合相结合的方法实时在线修正防空导弹武器命中概率,提高效率,增强工程操作性。
为实现此目的,本发明所设计的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于,它包括如下步骤:
步骤1:根据抗击飞机目标与防空导弹武器的距离R、防空导弹武器的发射远界R发射远界和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正
根据抗击飞机目标当前速度v、可抗击飞机目标最大速度vmax和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机目标速度修正后的命中概率P目标速度修正
根据抗击飞机目标基准机动系数A分别计算飞机巡航高度修正机动性系数A1,预设飞行员训练水平修正机动性系数A2,根据抗击飞机目标载荷系数fG和飞行员训练水平修正机动性系数A2计算抗击飞机目标载荷修正机动性系数A3,根据抗击飞机目标毁伤比例fD和飞机载荷修正机动性系数A3计算抗击飞机目标毁伤修正机动性系数 A4,根据抗击飞机目标毁伤修正机动性系数A4计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5,根据防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5计算抗击飞机目标机动性修正后的命中概率P机动系数修正
步骤2:根据抗击飞机目标处于的突防海拔高度H抗击飞机目标突防海拔和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算防空导弹武器抗击飞机目标修正后的命中概率 P抗击飞机目标突防海拔高度修正
根据抗击制导目标处于的突防海拔高度H抗击制导目标突防海拔和作战仿真系统数据库提供的抗击制导目标基准命中概率P抗击制导目标基准计算防空导弹武器抗击制导目标修正后的命中概率 P抗击制导目标突防海拔高度修正
步骤3:根据防空导弹抗击制导目标的方位角q计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正
本发明提出了一种防空导弹武器命中概率工程修正方法,高效利用获取的抗击距离、目标速度、机动系数、突防海拔高度、目标信号特征,基于作战仿真系统,采用数学模型和数值拟合相结合的方法实现了飞机目标和制导目标基准命中概率的快速准确修正。
本发明与现有技术相比的优点在于:传统的命中概率的分析方法需要大量实际数据,采用离线处理方式,但是多数情况下并不能得到足够的实际数据,特别是对一些引进的或较新的型号。本发明基于作战仿真系统,采用数学解析模型和数值拟合相结合的方法实时在线修正防空导弹武器命中概率,提高效率,增强工程操作性。
附图说明
图1为抗击方位角修正飞机机动性,飞机T为圆心,防空导弹为M,飞机速度矢量
Figure GDA0003116267620000041
方向为0°方向,顺时针为正向,抗击方位角指飞机速度矢量
Figure GDA0003116267620000042
与飞机与防空导弹相对距离矢量
Figure GDA0003116267620000043
之间的夹角q;
图2为抗击方位角修正制导目标,制导目标T为圆心,防空导弹为M,制导目标矢量
Figure GDA0003116267620000044
方向为0°方向,顺时针为正向,抗击方位角指制导目标速度矢量
Figure GDA0003116267620000045
与制导目标与防空导弹相对距离矢量
Figure GDA0003116267620000046
之间的夹角q。
图3为目标螺旋机动方式弹道。
其中,xL表示目标螺旋机动的螺旋轴心线,OL为目标开始螺旋机动在螺旋轴心线上对应的起点,zL、yL表示OL为原点、平行于地面坐标轴z、y的坐标轴,vs为目标沿xL轴运动的线速度,ω为目标绕xL轴的角速度,v为目标绕xL轴的线速度,LJ为目标沿xL轴的位移,h表示螺距,R表示螺旋半径,H0为螺旋轨迹最低点的高度。
图4为目标蛇形机动方式弹道。
图5为目标跃升俯冲方式弹道。
具体实施方式
以下结合附图和具体实施例对本发明作进一步的详细说明:
本发明设计的一种基于作战仿真系统的防空导弹武器命中概率修正方法,如图1和2所示,它包括如下步骤:
步骤1:根据抗击飞机目标与防空导弹武器的距离R、防空导弹武器的发射远界R发射远界和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正
根据抗击飞机目标当前速度v、可抗击飞机目标最大速度vmax和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机目标速度修正后的命中概率P目标速度修正
根据抗击飞机目标基准机动系数A分别计算飞机巡航高度修正机动性系数A1,预设飞行员训练水平修正机动性系数A2,根据抗击飞机目标载荷系数fG和飞行员训练水平修正机动性系数A2计算抗击飞机目标载荷修正机动性系数A3,根据抗击飞机目标毁伤比例fD和飞机载荷修正机动性系数A3计算抗击飞机目标毁伤修正机动性系数 A4,根据抗击飞机目标毁伤修正机动性系数A4计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5,根据防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5计算抗击飞机目标机动性修正后的命中概率P机动系数修正
步骤2:根据抗击飞机目标处于的突防海拔高度H抗击飞机目标突防海拔和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算防空导弹武器抗击飞机目标修正后的命中概率 P抗击飞机目标突防海拔高度修正
根据抗击制导目标处于的突防海拔高度H抗击制导目标突防海拔和作战仿真系统数据库提供的抗击制导目标基准命中概率P抗击制导目标基准计算防空导弹武器抗击制导目标修正后的命中概率 P抗击制导目标突防海拔高度修正
步骤3:根据防空导弹抗击制导目标的方位角q计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正
上述技术方案中,抗击制导目标为飞航导弹、再入飞行器RV、高超音速飞行器等制导武器。
上述技术方案中,抗击飞机目标和抗击制导目标均为防空导弹武器抗击的对象。
上述技术方案的步骤1中,防空导弹武器发射距离大于等于防空导弹武器发射远界的二分之一时,进行P抗击距离修正的计算;
所述步骤1中,根据抗击飞机目标与防空导弹武器的距离R,防空导弹武器的发射远界R发射远界计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正的具体方法为:
当PR≥Pf时,P抗击距离修正=P抗击飞机目标基准*(1+Pf-PR);
当R发射近界/R发射远界≤PR<Pf时,P抗击距离修正=P抗击飞机目标基准
其中,Pf是中间变量,Pf=0.5,PR=R/R发射远界,R为防空导弹武器与抗击飞机目标的距离,R发射远界为防空导弹武器的发射远界,P 抗击飞机目标基准为作战仿真系统数据库提供的抗击飞机目标基准命中概率。
上述技术方案的步骤1中,所述步骤1中,根据抗击飞机目标当前速度v、可抗击飞机目标最大速度vmax和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机目标速度修正后的命中概率P目标速度修正的具体方法为:
Figure GDA0003116267620000071
上述技术方案的步骤1中,所述步骤1中,根据抗击飞机目标基准机动系数A计算飞机巡航高度修正机动性系数A1的方法为:
Figure GDA0003116267620000072
其中,HEmax为抗击飞机目标发动机最大工作高度, H抗击飞机目标海拔为抗击飞机目标当前海拔高度。超机动性是超音速机动和过失速机动的综合,一般有标志性动作:回形针动作、眼镜蛇机动、尾冲机动、落叶飘等,一般只有四代半和五代战机具备超机动性。
上述技术方案中,预设飞行员训练水平修正机动性系数A2为:
新手A2=A1*0.3;实习A2=A1*0.5;普通A2=A1*0.8;老手 A2=A1*1.0;顶级A2=A1*1.2。
上述技术方案的步骤1中,所述步骤1中,根据抗击飞机目标载荷系数fG和飞行员训练水平修正机动性系数A2计算抗击飞机目标载荷修正机动性系数A3的方法为:
如果抗击飞机目标的挂载和油料为0就不修正,否则修正,根据以下公式计算抗击飞机目标的载荷系数fG
抗击飞机目标重量Gc由空重GEmpty、挂载载荷重量GPayload和燃料重量GFuel三部分组成,即
Gc=GEmpty+GPayload+GFuel
设抗击飞机目标最大载荷重量为Gmax,则抗击飞机目标有效载荷重量Gvalid
Gvalid=Gmax-(GEmpty+0.6GFuel)
则抗击飞机目标的载荷系数fG为:
Figure GDA0003116267620000081
抗击飞机目标载荷修正机动性系数A3
A3=[0.4+0.6(1-fG)]A2
上述技术方案的步骤1中,所述步骤1中,根据抗击飞机目标毁伤比例fD(飞机被抗击下的完好程度)和飞机载荷修正机动性系数 A3计算抗击飞机目标毁伤修正机动性系数A4的方法为:
A4=(1-fD)A3
根据抗击飞机目标毁伤修正机动性系数A4计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5的方法为:
当防空导弹抗击飞机的方位角范围在345°~360°或0°~15°内时, A5=0.6*A4,即得到防空导弹武器迎头攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在15°~60°或300°~345°内时, A5=0.7*A4,即得到防空导弹武器前向攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在60°~110°或250°~300°内时,A5=1*A4,此时防空导弹武器的侧向攻击对抗击飞机目标机动性系数无影响;
当防空导弹抗击飞机的方位角范围在110°~165°或195°~250°内时,A5=0.85*A4,即得到防空导弹武器后向攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在165°~195°内时,A5=0.5*A4,即得到防空导弹武器尾追攻击效应修正后的抗击飞机目标机动性系数;
根据防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5计算抗击飞机目标机动性修正后的命中概率P机动系数修正的方法为:
P机动系数修正=P抗击飞机目标基准–10%*A5
防空导弹抗击飞机目标方位角为飞机速度矢量与飞机-防空导弹相对距离矢量之间的角度。
上述技术方案的步骤2中,当抗击目标处于不同突防海拔高度 H海拔(单位m),且抗击武器不具备抗击掠海目标能力时,距离海平面越近,命中概率下降越大,根据抗击飞机目标处于的突防海拔高度H抗击飞机目标突防海拔和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算防空导弹武器抗击飞机目标修正后的命中概率P抗击飞机目标突防海拔高度修正的方法为:
Figure GDA0003116267620000091
根据抗击制导目标处于的突防海拔高度H抗击制导目标突防海拔和作战仿真系统数据库提供的抗击制导目标基准命中概率P抗击制导目标基准计算防空导弹武器抗击制导目标修正后的命中概率 P抗击制导目标突防海拔高度修正的方法为:
Figure GDA0003116267620000092
上述技术方案的步骤3中,所述步骤3中,根据防空导弹抗击制导目标的方位角q计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正的方法为:
Figure GDA0003116267620000101
上述技术方案的步骤3中,还包括根据抗击制导目标末端机动能力,结合预测命中点与防空导弹武器发射点的水平距离R预测计算抗击制导目标末端机动突防修正的抗击制导目标命中概率 P目标末端机动突防修正
Figure GDA0003116267620000102
所述步骤3中,如果防空导弹武器为雷达制导方式,根据抗击制导目标雷达截射面积RCS,计算目标信号特性修正后的抗击制导目标命中概率P目标信号特征修正的具体方法为:
Figure GDA0003116267620000103
如果防空导弹武器为红外制导方式,根据抗击制导目标的红外探测距离RD,计算目标信号特性修正后的抗击制导目标命中概率 P目标信号特征修正的具体方法为:
Figure GDA0003116267620000104
实施例:抗击飞机目标时,由上述技术方案的步骤1分别计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正、抗击飞机目标速度修正后的命中概率P目标速度修正、计算抗击飞机目标机动性修正后的命中概率P机动系数修正,由上述技术方案的步骤2计算防空导弹武器抗击飞机目标修正后的命中概率P抗击飞机目标突防海拔高度修正,相乘得到抗击飞机目标的修正命中概率。
抗击制导目标时,由上述技术方案的步骤2计算防空导弹武器抗击制导目标修正后的命中概率P抗击制导目标突防海拔高度修正,由上述技术方案的步骤3计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正,相乘得到抗击制导目标的修正命中概率。
本实施例中,作战仿真系统采用基于HLA体系架构的模拟海空军联合行动的智能化联合作战推演系统。
计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正
某型空对空导弹的推进系统为火箭发动机,因此修正适用条件为当发射距离大于等于发射远界的二分之一,则进行距离修正。
在此次抗击空中突击的过程中,某型空对空导弹的基准命中概率为90%,发射远界为92.6km,发射远界的二分之一为46.3km,而发射单元(纬度latitude=9.58145963313289,经度 longitude=112.808972261327)与目标(latitude=9.73616712165069,longitude=112.276743341381)的距离为32.85864nm,即60.8542km,大于发射远界的二分之一,因此满足修正条件,根据距离修正公式:P抗击距离修正=P抗击飞机目标基准*(1+Pf-PR)
可知,P抗击距离修正=75.85449%;
计算抗击飞机目标速度修正后的命中概率P目标速度修正
目标当前速度为1179.724km/h,大于等于其最大速度 (1713.1km/h)的60%(68.8%),因此根据修正条件:
P目标速度修正=P抗击飞机目标基准–0.10(V≥0.6Vmax)
可知,P目标速度修正=命中概率为75.75%;
计算飞机巡航高度修正机动性系数A1
目标苏-27SK型战斗机的基准机动系数A为4.5,目标的飞行高度为24.384m,低于海拔3000m,因此机动系数A1无需高度修正,为4.5;
预设飞行员训练水平修正机动性系数A2为:
训练水平为普通,因此根据修正公式:
A2=A1*0.8
训练水平修正后的机动系数A2为3.6;
计算抗击飞机目标载荷修正机动性系数A3
从输出消息中已知目标的重量分数fG为0.32,因此根据修正公式:
A3=[0.4+0.6(1-fG)]A2
飞机重量修正后的机动系数A3为2.92;
计算抗击飞机目标毁伤修正机动性系数A4
目标在攻击前的毁伤为0,因此飞机毁伤修正后机动系数A4保持不变,为2.92;
计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5
目标相对武器的方位角的126°,因此根据修正公式:
A5=0.85*A4
后向攻击效应修正后的机动系数A5为2.5;
计算抗击飞机目标机动性修正后的命中概率P机动系数修正
P机动系数修正=P-A5*10%
可知,经机动系数修正后的命中概率为51%。
对掠海目标抗击能力修正:
由于该型防空武器具有抗击掠海目标的能力,因此不适用本项修正条件,即武器的命中概率保持不变,为51%。
在某次抗击飞机的过程中,系统输出消息如下:
某型空对空导弹正在攻击目标(苏-27SK型战斗机,基准命中概率:90%,命中概率(经距离修正):76%,命中概率(经目标速度修正): 75.8%,苏27目标基准机动系数:4.5,高度修正:4.5,训练水平修正(普通):3.6,飞机重量分数为0.32,机动系数修正为2.92,机动系数(后向攻击效应修正):2.5,命中概率(机动系数修正量):-25%,最终命中概率:51%。
对制导武器命中概率的修正:
经查询,某型防空导弹具有对掠海目标的抗击能力,因此不适用本项修正。目标为AS-18A型“芦笛”反舰导弹,是制导武器,因此适用本项修正;
武器与目标之间角度修正:
武器某型防空导弹的基准命中概率是80%,武器与目标之间的方位角是0°,因此根据修正公式:
P抗击方位角修正=P抗击制导目标基准*(1.0–q*0.5/90.0)
可知,经武器与目标之间角度修正后的命中概率为80%。
武器突防措施修正:
AS-18A型“芦笛”反舰导弹采用跃升俯冲机动方式,预测命中点与某型防空导弹武器发射点的水平距离R预测为13.85264nm,大于 4.5km,因此不适用本项修正。
目标信号特征修正:
武器某型防空导弹的制导方式为半主动雷达寻的,因此需根据目标雷达信号强度对命中概率进行修正,由输出消息可知,武器与目标之间的角度为0°,目标信号特征取0.12m2,根据修正公式:
P目标信号特征修正=P抗击制导目标基准–10%(0.1<RCS≤1)
可知,经目标信号特征修正的命中概率为70%。
在某次抗击制导武器的过程中,系统输出消息如下:
武器:某型防空导弹正在攻击目标AS-18A型“芦笛”反舰导弹,基准命中概率:80%,武器与目标之间角度为0°,概率修正为80%,目标信号特征修正:-10%,最终命中概率:70%。
本说明书未作详细描述的内容属于本领域专业技术人员公知的现有技术。

Claims (10)

1.一种基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于,它包括如下步骤:
步骤1:根据抗击飞机目标与防空导弹武器的距离R、防空导弹武器的发射远界R发射远界和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正
根据抗击飞机目标当前速度v、可抗击飞机目标最大速度vmax和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机目标速度修正后的命中概率P目标速度修正
根据抗击飞机目标基准机动系数A分别计算飞机巡航高度修正机动性系数A1,预设飞行员训练水平修正机动性系数A2,根据抗击飞机目标载荷系数fG和飞行员训练水平修正机动性系数A2计算抗击飞机目标载荷修正机动性系数A3,根据抗击飞机目标毁伤比例fD和飞机载荷修正机动性系数A3计算抗击飞机目标毁伤修正机动性系数A4,根据抗击飞机目标毁伤修正机动性系数A4计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5,根据防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5计算抗击飞机目标机动性修正后的命中概率P机动系数修正
步骤2:根据抗击飞机目标处于的突防海拔高度H抗击飞机目标突防海拔和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算防空导弹武器抗击飞机目标修正后的命中概率P抗击飞机目标突防海拔高度修正
根据抗击制导目标处于的突防海拔高度H抗击制导目标突防海拔和作战仿真系统数据库提供的抗击制导目标基准命中概率P抗击制导目标基准计算防空导弹武器抗击制导目标修正后的命中概率P抗击制导目标突防海拔高度修正
步骤3:根据防空导弹抗击制导目标的方位角q计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正
2.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,根据抗击飞机目标与防空导弹武器的距离R,防空导弹武器的发射远界R发射远界计算抗击飞机的抗击距离修正后的命中概率P抗击距离修正的具体方法为:
当PR≥Pf时,P抗击距离修正=P抗击飞机目标基准*(1+Pf-PR);
当R发射近界/R发射远界≤PR<Pf时,P抗击距离修正=P抗击飞机目标基准
其中,Pf是中间变量,Pf=0.5,PR=R/R发射远界,R为防空导弹武器与抗击飞机目标的距离,R发射远界为防空导弹武器的发射远界,P抗击飞机目标基准为作战仿真系统数据库提供的抗击飞机目标基准命中概率。
3.根据权利要求2所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,防空导弹武器发射距离大于等于防空导弹武器发射远界的二分之一时,进行P抗击距离修正的计算。
4.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,根据抗击飞机目标当前速度v、可抗击飞机目标最大速度vmax和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算抗击飞机目标速度修正后的命中概率P目标速度修正的具体方法为:
Figure FDA0003116267610000021
5.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,根据抗击飞机目标基准机动系数A计算飞机巡航高度修正机动性系数A1的方法为:
Figure FDA0003116267610000031
其中,HEmax为抗击飞机目标发动机最大工作高度,H抗击飞机目标海拔为抗击飞机目标当前海拔高度。
6.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,根据抗击飞机目标载荷系数fG和飞行员训练水平修正机动性系数A2计算抗击飞机目标载荷修正机动性系数A3的方法为:
如果抗击飞机目标的挂载和油料为0就不修正,否则修正,根据以下公式计算抗击飞机目标的载荷系数fG
抗击飞机目标重量Gc由空重GEmpty、挂载载荷重量GPayload和燃料重量GFuel三部分组成,即
Gc=GEmpty+GPayload+GFuel
设抗击飞机目标最大载荷重量为Gmax,则抗击飞机目标有效载荷重量Gvalid
Gvalid=Gmax-(GEmpty+0.6GFuel)
则抗击飞机目标的载荷系数fG为:
Figure FDA0003116267610000032
抗击飞机目标载荷修正机动性系数A3
A3=[0.4+0.6(1-fG)]A2
7.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤1中,根据抗击飞机目标毁伤比例fD和飞机载荷修正机动性系数A3计算抗击飞机目标毁伤修正机动性系数A4的方法为:
A4=(1-fD)A3
根据抗击飞机目标毁伤修正机动性系数A4计算防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5的方法为:
当防空导弹抗击飞机的方位角范围在345°~360°或0°~15°内时,A5=0.6*A4,即得到防空导弹武器迎头攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在15°~60°或300°~345°内时,A5=0.7*A4,即得到防空导弹武器前向攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在60°~110°或250°~300°内时,A5=1*A4,此时防空导弹武器的侧向攻击对抗击飞机目标机动性系数无影响;
当防空导弹抗击飞机的方位角范围在110°~165°或195°~250°内时,A5=0.85*A4,即得到防空导弹武器后向攻击效应修正后的抗击飞机目标机动性系数;
当防空导弹抗击飞机的方位角范围在165°~195°内时,A5=0.5*A4,即得到防空导弹武器尾追攻击效应修正后的抗击飞机目标机动性系数;
根据防空导弹抗击飞机目标方位角修正的抗击飞机目标机动性系数A5计算抗击飞机目标机动性修正后的命中概率P机动系数修正的方法为:
P机动系数修正=P抗击飞机目标基准–10%*A5
8.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤2中,根据抗击飞机目标处于的突防海拔高度H抗击飞机目标突防海拔和作战仿真系统数据库提供的抗击飞机目标基准命中概率P抗击飞机目标基准计算防空导弹武器抗击飞机目标修正后的命中概率P抗击飞机目标突防海拔高度修正的方法为:
Figure FDA0003116267610000051
根据抗击制导目标处于的突防海拔高度H抗击制导目标突防海拔和作战仿真系统数据库提供的抗击制导目标基准命中概率P抗击制导目标基准计算防空导弹武器抗击制导目标修正后的命中概率P抗击制导目标突防海拔高度修正的方法为:
Figure FDA0003116267610000052
9.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤3中,根据防空导弹抗击制导目标的方位角q计算抗击制导目标方位角修正的抗击制导目标命中概率P抗击方位角修正的方法为:
Figure FDA0003116267610000053
10.根据权利要求1所述的基于作战仿真系统的防空导弹武器命中概率修正方法,其特征在于:所述步骤3中,还包括根据抗击制导目标末端机动能力,结合预测命中点与防空导弹武器发射点的水平距离R预测计算抗击制导目标末端机动突防修正的抗击制导目标命中概率P目标末端机动突防修正
Figure FDA0003116267610000054
Figure FDA0003116267610000061
所述步骤3中,如果防空导弹武器为雷达制导方式,根据抗击制导目标雷达截射面积RCS,计算目标信号特性修正后的抗击制导目标命中概率P目标信号特征修正的具体方法为:
Figure FDA0003116267610000062
如果防空导弹武器为红外制导方式,根据抗击制导目标的红外探测距离RD,计算目标信号特性修正后的抗击制导目标命中概率P目标信号特征修正的具体方法为:
Figure FDA0003116267610000063
CN202011278680.5A 2020-11-16 2020-11-16 基于作战仿真系统的防空导弹武器命中概率修正方法 Active CN112464451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011278680.5A CN112464451B (zh) 2020-11-16 2020-11-16 基于作战仿真系统的防空导弹武器命中概率修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011278680.5A CN112464451B (zh) 2020-11-16 2020-11-16 基于作战仿真系统的防空导弹武器命中概率修正方法

Publications (2)

Publication Number Publication Date
CN112464451A CN112464451A (zh) 2021-03-09
CN112464451B true CN112464451B (zh) 2021-08-13

Family

ID=74836908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011278680.5A Active CN112464451B (zh) 2020-11-16 2020-11-16 基于作战仿真系统的防空导弹武器命中概率修正方法

Country Status (1)

Country Link
CN (1) CN112464451B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117390830A (zh) * 2023-09-25 2024-01-12 中国人民解放军海军航空大学 一种无人飞行器跨平台制导仿真方法、系统及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE8704096L (sv) * 1987-10-22 1989-04-23 Bofors Ab Saett att oeka traeffsannolikheten
US5836540A (en) * 1994-03-25 1998-11-17 Rheinmetall W & M Gmbh Projectile having an apparatus for flight-path correction
CN105910495A (zh) * 2016-05-09 2016-08-31 哈尔滨工业大学 基于性能指标的面向效能的导弹武器系统设计方法
CN110414061A (zh) * 2019-06-29 2019-11-05 西南电子技术研究所(中国电子科技集团公司第十研究所) 建立飞机平台系统对地攻击生存能力概率模型的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1005871B (zh) * 1987-07-25 1989-11-22 空军第一航空技术专科学校 飞机激光空中射击模拟器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE8704096L (sv) * 1987-10-22 1989-04-23 Bofors Ab Saett att oeka traeffsannolikheten
US5836540A (en) * 1994-03-25 1998-11-17 Rheinmetall W & M Gmbh Projectile having an apparatus for flight-path correction
CN105910495A (zh) * 2016-05-09 2016-08-31 哈尔滨工业大学 基于性能指标的面向效能的导弹武器系统设计方法
CN110414061A (zh) * 2019-06-29 2019-11-05 西南电子技术研究所(中国电子科技集团公司第十研究所) 建立飞机平台系统对地攻击生存能力概率模型的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HLA海上编队对抗仿真系统构建;刘方 等;《火力与指挥控制》;20080915;第33卷(第9期);第29-32页 *
Simulation of covariance analysis describing equation technique (CADET) in missile hit probability calculation;Xinsheng Wang 等;《2010 Sixth International Conference on Natural Computation》;20200923;第4282-4285页 *

Also Published As

Publication number Publication date
CN112464451A (zh) 2021-03-09

Similar Documents

Publication Publication Date Title
CN111649624B (zh) 一种空间微型精确制导武器控制方法
CN109902933B (zh) 一种反舰战斗部打击水面舰船毁伤等级的评价方法
CN113359813B (zh) 一种多领弹多从弹群组协同制导方法
CN109857145B (zh) 一种基于迭代预测命中点的增程型拦截弹预测制导方法
CN111027862B (zh) 一种基于多维度的层次聚合作战模拟训练评估方法
CN102314537A (zh) 一种掠飞击顶灵巧弹药毁伤概率计算方法
CN108073742B (zh) 基于改进粒子滤波算法的拦截导弹末段飞行状态估计方法
CN112464451B (zh) 基于作战仿真系统的防空导弹武器命中概率修正方法
CN112836383B (zh) 一种基于靶式数据修正的武器效能建模方法
Liang et al. Multi-attacks effectiveness evaluation of UUV based on wake guidance
CN116360489A (zh) 一种基于成型制导律的协同中制导律设计方法
CN114548674B (zh) 面向多智能体对抗场景的威胁态势评估方法、装置及设备
CN112396288B (zh) 一种装备体系火力打击能力评估方法、装置、设备和介质
RU2602162C2 (ru) Способ стрельбы реактивными снарядами реактивной системы залпового огня в условиях контрбатарейной борьбы
CN114819055A (zh) 一种基于lstm网络的导弹拦截点预测方法
CN111832193A (zh) 一种超视距空战中空中目标的威胁评估方法
Liu et al. Research and analysis on firing accuracy of naval gun
Zuoe et al. Study on vertical attack penetration probability of anti-ship missile
CN112729015B (zh) 一种基于拦截弹导引律辨识的机动突防方法
KR102489645B1 (ko) 탄착군 확률 모델 기반 대공표적 실시간 명중평가 방법 및 장치
RU2794300C2 (ru) Способ получения вероятностной оценки возможности преодоления зон поражения зенитных управляемых ракет маневрирующим беспилотным летательным аппаратом
Herrmann Air-to-air missile engagement analysis using the USAF Trajectory Analysis Program (TRAP)
Wei et al. The fire control calculation of new shipborne gun to sea target based on terminal impact angle constraint
Lee et al. A METHOD OF CALCULATING WEAPON EFFECTIVENESS INDEX USING MUNITION TRAJECTORY AND VULNERABILITY DATA
CN112861330A (zh) 一种基于matlab的导弹杀伤效果计算和可视化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant