CN110414061A - 建立飞机平台系统对地攻击生存能力概率模型的方法 - Google Patents
建立飞机平台系统对地攻击生存能力概率模型的方法 Download PDFInfo
- Publication number
- CN110414061A CN110414061A CN201910580684.XA CN201910580684A CN110414061A CN 110414061 A CN110414061 A CN 110414061A CN 201910580684 A CN201910580684 A CN 201910580684A CN 110414061 A CN110414061 A CN 110414061A
- Authority
- CN
- China
- Prior art keywords
- probability
- aircraft
- indicate
- receiver
- operational
- 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.)
- Granted
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开的一种建立飞机平台系统对地攻击生存能力概率模型的方法,旨在提供一种实用性强,能够清楚表达飞机生存力和突防成功的生存概率的建模方法。本发明通过下述技术方案予以实现:判断飞机是否落入对方探测系统有效搜索空域,将飞机所在轨迹的地理坐标转入对方每个地面防空探测系统的探测坐标系,解出作战飞机对探测装置的相对位置,判断出当前时刻作战飞机是否落入对方探测系统有效搜索空域,这个事件出现概率就是飞机进入对方探测系统有效搜索空域的概率,结合飞机在执行任务过程中出现故障的概率,计算作战活动的每个阶段的每个时刻生存概率,从而构建飞机平台系统生存能力概率模型。
Description
技术领域
本发明涉及一种建立飞机平台系统对地攻击生存能力概率模型的方法,也适用于生存效能、突防效能、杀伤效能、任务完成效能的评估。
背景技术
按照评估的方式和手段,可以将作战飞机对地攻击作战效能的评估方法分为五大类:飞行性能对比法、解析计算法、空战仿真法、专家评估法和实验统计法等。解析计算法的特点是根据描述效能指标与给定条件之间的函数关系的解析表达式计算指标值。这个解析表达式可以是直接根据军事运筹理论建立的,也可以是用数学方法求解所建立的效能方程得到的。解析计算法按其计算特点和评估方式可分为参数计算法、概率分析法和需要量估算法3类。参数计算法是最常用的效能评估方法。它还可分为对数法、顺序评估法、相对值评估法、相对指数法和多参数(品质)分析法等。概率分析法是利用概率的基本原理和方法,对武器装备在现代战争中使用的效果进行分析和预测,概率法是当前在武器装备论证及作战效能评估中常用的方法。一般来说,凡是涉及到处理不确定因素的问题都可采用概率法。
生存性是作战飞机在战场上“保存自己”的能力体现,也是进一步消灭敌人的前提和基础,它直接影响到作战平台持续战斗能力。随着装备设计要求的不断提高和现代空战的不断发展,生存力已经越来越引起人们的重视。生存力是现代飞机作战效能评估与平台设计的重要指标之一,可用生存概率来度量,系飞机躲避和(或)抵抗人为敌对环境的能力。要想击败对方,飞机不但要求具备承受损伤的能力,而且要有避免被威胁目标命中的能力,它属于特定的作战环境下飞机敏感性和易损性的综合。其中敏感性指飞机的特性,在对方武器威胁的环境下,表示飞机难于躲避雷达、威胁传播物等一切人为敌对环境因素的能力。易损性是关于飞机承受打击的能力。飞机的敏感性、易损性越高,则飞机的生存力越低。在对地攻击作战样式下,飞机被发现概率与飞机进入对方探测系统有效搜索空域后对方的探测手段有关,典型的探测手段有截获我机辐射(或散射)的电磁波和红外直接探测两种。于是,飞机被发现概率等于飞机进入对方搜索空域的概率与对方截获我机电磁波概率或者是红外感知我方概率的乘积;飞机被识别的概率是飞机进入对方识别系统作用空域的概率与飞机被对方识别系统判别为对方威胁的概率的乘积;飞机被跟踪的概率是以下三类概率的乘积,它们分别是飞机进入对方跟踪系统作用空域的概率、飞机被捕获的概率以及飞机被确认跟踪的概率;目前运用概率分析方法建立的飞机生存效能模型文献不多,有飞机在电子对抗情况下的生存概率计算模型,它表示为几类事件概率的乘积,包括邻近飞机的威胁武器系统活动概率、威胁处于活动状态时飞机被探测的条件概率、在探测到飞机情况下飞机被追踪并向飞机发射导弹或枪弹的条件概率和传播物已经截击飞机下击中飞机的概率乘积;有突防飞机在雷达网检测和防空武器系统攻击情况下进行突防作战的生存概率模型,通过1减去飞机被击毁的概率得到飞机生存概率;对对方空防系统实施电子干扰时飞机生存概率是1减去飞机被击毁的概率与对方雷达一次扫描对飞机的发现概率的乘积。所有模型及其参量的计算只是按典型情况给出的,但都具有一个共同点,就是在确定了击毁概率以后才计算飞机的生存概率,且生存概率计算不准确,生存概率模型构建过于理想,以及对生存概率影响因素考虑不完整等不足。
发明内容
本发明的目的是针对现有技术存在的不足之处,提供一种方法科学、实用性强,可以进行对抗环境、对抗双方变量间关系分析,能够清楚表达飞机生存力和突防成功的生存概率的建模方法。
本发明的上述目的可以通过下述技术方案予以实现:一种建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于包括如下步骤,根据地面防空探测区域为防空探测系统每个探测单元同时探测区域Xi(i=1,…,m)的并集(X1∪X2∪…∪Xm)构成的包络区域,分析地面固定目标所在防空探测系统对我飞机的发现概率、识别概率、跟踪概率和毁伤概率,判断飞机是否落入对方探测系统有效搜索空域,判断当前时刻作战飞机是否存在包络空域区域内,将作战飞机所在轨迹的地理坐标分别转入对方每个地面防空探测系统的探测坐标系,解出作战飞机对探测装置的相对位置,判断出当前时刻作战飞机是否落入对方探测系统有效搜索空域;然后根据飞机在执行任务过程中影响生存概率主要因素的基础上,计算出飞机在执行对地攻击任务过程中出现故障的概率、被发现概率,被识别概率和被跟踪的概率,计算这些阶段事件出现概率的乘积,结合飞机在执行任务过程中出现故障的概率,计算作战活动的每个阶段的每个时刻生存概率,给出飞机生存概率的计算过程,得到飞机生存概率,依据飞机生存概率,建立飞机平台系统生存能力概率模型。
本发明相比于现有技术具有如下有益效果,
本发明根据飞机在执行任务过程中影响生存概率主要因素的基础上,给出飞机生存概率的计算过程和分析方法,建立飞机平台系统生存能力概率模型,弥补了以往文献中装备生存概率计算不准确,生存概率模型构建过于理想,以及对生存概率影响因素不完整等缺陷。
本发明通过研究作战飞机对地面固定目标进行攻击事件过程分析,综合考虑了在攻防对抗过程中作战飞机生存能力的影响因素,提出了飞机平台系统生存能力概率模型。通过计算作战活动的每个阶段的每个时刻生存概率,可以帮助机载任务系统实时地调整任务规划,使电子信息系统和武器系统最大效用地发挥作用,有效回避各种威胁和障碍,提高了飞机生存能力和突防任务的成功率。
本发明针对作战飞机对地面固定目标进行攻击时的效能评估问题,在分析了地面固定目标所在的防空作战系统对我作战飞机的发现概率、跟踪概率、识别概率和毁伤概率的基础上,结合我作战飞机在执行任务过程中出现致命故障的概率,提出单平台作战飞机对地攻击生存概率概念以及计算模型,该计算模型适用于常规战斗机对地攻击时的生存效能的评估,也适用于生存效能、突防效能、杀伤效能、任务完成效能的评估。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明的建立飞机平台系统生存能力概率模型的流程图。
具体实施方式
参阅图1。根据本发明,飞机进入对方的防空探测系统有效搜索空域的概率与对方探测设备的引导概率、导航系统精度、地面防空探测设备特性和目标区地形/气象等条件都有关系。地面防空探测区域为防空探测系统每个探测单元同时探测区域Xi(i=1,…,m)的并集(X1∪X2∪…∪Xm)构成的包络区域。判断飞机是否落入对方探测系统有效搜索空域也即判断当前时刻作战飞机是否存在于上述包络空域区域内部,其方法是将作战飞机所在轨迹的地理坐标分别转入对方每个地面防空探测系统的探测坐标系,解出作战飞机对探测装置的相对位置,从而可判断该时刻作战飞机是否落入对方探测系统有效搜索空域。
分析地面固定目标所在防空探测系统对我飞机的发现概率、识别概率、跟踪概率和毁伤概率;根据飞机在执行对地攻击任务过程中影响生存概率主要因素的基础上,生存概率模型子系统POSM计算出飞机在执行任务过程中被发现概率、被识别概率、被跟踪的概率和被毁伤概率,计算这些阶段事件出现概率的乘积,结合飞机在执行任务过程中出现故障的概率,计算作战活动的每个阶段的每个时刻生存概率,给出飞机生存概率的计算过程,得到飞机生存概率;依据飞机生存概率,建立飞机平台系统生存能力概率模型。在计算生存概率中,首先计算得到被发现概率、被识别的概率、被跟踪的概率和被毁伤概率,然后计算得到飞机在当前时刻的生存概率,使用飞机被毁伤概率计算公式计算我作战飞机被对方毁伤概率Pk,
飞机在t0时刻生存概率的计算公式为:P=(1-Pf)·(1-Pr·Pi·Pt·Pk)
式中,n是对方防空武器发射武器总数目,Pka(i)是我作战飞机进入对方第i种防空武器作用空域的概率,Pht(i)是我作战飞机被对方第i种防空武器击中的概率,Pmr(i)是我作战飞机被对方第i种防空武器击中条件下的损毁概率;P表示飞机生存概率,Pf表示飞机在执行任务过程中出现致命故障的概率,Pr表示飞机被对方发现的概率,Pi表示飞机被对方发现条件下被对方识别的概率,Pt表示飞机被对方发现条件下被对方跟踪的概率,Pk表示飞机被对方击毁的概率。
S1:生存概率模型子系统POSM计算飞机故障概率Pf,
并且
其中,Pf表示飞机故障概率,MTBFf表示飞机在执行任务的时候出现故障的平均间隔时间,表示系统处于正常工作状态的数量特征;MTTRf表示飞机在执行任务的时候故障平均修复时间,表征系统出于故障状态的数量特征,(MTTRf)i表示构成飞机的第i个可更换的部件LRU完成一次修复性维修所需要的时间;λi是构成飞机的第i个可更换的部件LRU的故障率;n表示构成飞机的第i个可更换的部件LRU的总数。
S2:POSM根据被发现概率是作战飞机的电磁信号特征或辐射散射特征被对方防空探测系统中所有侦察探测传感器单元联合探测的发现概率,POSM计算飞机被发现概率Pr,Pr=Pea·Pei
式中:Pea是飞机进入对方探测系统有效搜索空域的概率,Pei是飞机辐射或散射的电磁波被对方探测系统截获的概率或对方红外探测系统对作战飞机红外辐射探测概率。
S21:POSM计算飞机进入对方探测系统有效搜索空域的概率Pea,假设作战飞机在预警机数据链导引下飞往对方的防空探测区域s,在时刻t(t∈T)进入目标防空探测区s的概率为Pea,设对方的防空探测传感器互相独立,互不影响,在时间t(t∈T)进入目标防空探测区s的概率为Pea,则
式中,P(X1∪X2∪…∪Xm)为飞机落入防空探测系统每个探测单元在同一时刻的探测区域Xi(i=1,…,m)的并集(X1∪X2∪…∪Xm)构成的包络区域的概率。飞机在预警机数据链导引下飞往对方的目标防空探测区s的第i个探测设备的探测区域Xi中出现的联合概率密度,假设对方探测设备在二维空域进行搜索,即在方位角和俯仰角范围内搜索,则
其中,是飞机在目标防空探测区s的第i个探测设备的探测坐标系的真实相对位置,是对方的目标防空探测区s的第i个探测设备对我飞机方位角和俯仰角的估计值,为对方的目标防空探测区s的第i个探测设备对我飞机探测方位角和俯仰角的误差的标准方差,为我飞机落入对方第i个探测设备的探测区域Xi内的概率,此时,Xi(i=1,…,m)为对方第i个探测设备的 为对方第i个探测设备的方位覆盖范围,为对方第i个探测设备的俯仰覆盖范围。
S22:POSM计算飞机辐射或反射的电磁波被对方探测系统截获的概率或对方红外探测系统对我机红外辐射探测概率Pei,假设对方防空区域内对我飞机探测的雷达数估值为并赋予定义并设每部雷达的探测概率相同;对方防空区域内对我飞机电子侦察传感器数估值为并赋予定义并设每部电子侦察传感器的侦察概率相同;对方防空区域内对我飞机红外探测传感器数估值为并赋予定义并设每部红外探测传感器的探测概率相同,则我飞机被对方防空探测传感器探测概率为:
S221:POSM计算散射电磁波被对方探测系统截获概率P1,假设Z事件“雷达截获到我方飞机RCS波瓣σRCS(θ)≠0”,H1表示事件“雷达发射波束对准目标”;POSM根据对我方飞机RCS截获概率P(Z|Pr≥Smin,H1)和对我方飞机的检测概率Pd,计算我机散射的电磁波被对方探测系统截获的概率P1,且
P1=P(Pr≥Smin|H1)×P(Z|Pr≥Smin,H1)≈Pd·P(Z|Pr≥Smin,H1),其中,Pr为对方雷达发射信号被我飞机反射后到达对方截获接收机的信号功率,Smin为对方接收机识别信号的灵敏度门限。
S222:POSM根据第k台接收机发生一次截获的截获概率Pi_k=Pd·Ps·Pt·Pf计算辐射电磁波被对方探测系统截获的概率P2,
其中,n表示辐射波束区域内存在n部接收器;Pi_cum为整个波束覆盖区内发生一次截获的截获概率;Pd表示对方接收机的检测概率,即到达对方接收机的信号功率超过敌接收机最小可检测的门限功率时,对方接收机探测到我发射机波束能量的概率;Ps表示单个对方接收机的空域截获概率,即我发射机波束照射到对方飞机所在位置的概率;Pt表示单个对方接收机的时间截获概率,表示我方发射机设备正在辐射信号的时间内,对方接收设备处于开机状态的概率;Pf表示单个对方接收机的频率截获概率,即我发射机的当前工作频率正好落入对方接收机瞬时测频带宽,且满足对信号的测频条件的概率。
POSM将对方接收机的检测概率Pd、对方接收机的空域截获概率Ps、对方接收机的时间截获概率Pt和对方接收机的频率截获概率Pf综合起来,得到第k台接收机发生一次截获的截获概率计算公式为:
式中,MF表示主瓣覆盖面积(3dB);DI表示截获接收机密度;TOT表示我辐射源照射对方接收机时间;TI表示截获接收机的扫描时间。
S223:POSM计算对方红外探测系统对我飞机红外辐射探测概率P3,
并且η=1-exp(-λT),当对方红外系统为光机扫描型时,红外系统扫描周期T=1/τd,为凝视型时,T=1/tmin,无干扰条件下红外探测系统对目标的发现概率按下式计算第i次单次发现概率Pdi:TNR可以由虚警概率来获得:
式中,Ph为红外探测系统N次观察中最大的单次发现概率,λ表示红外探测系统工作波长,N表示红外探测系统的观测次数,T为红外系统扫描周期,τd为扫描单帧的时间;tmin为探测器积分时间,SNR对方红外探测器的信号噪声比,TNR为阈值信噪比,也称最小可检测信噪比,Pf为虚警概率,一般取Pf=10-6。
S3:POSM计算飞机被对方识别的概率Pi,Pi=Pra·Ptt
式中,Pra表示飞机进入对方识别系统作用空域的概率;Ptt表示飞机被对方识别系统判别为对方威胁的概率。
S31:飞机进入对方识别系统作用空域的概率Pra。
POSM判断我飞机是否落入对方综合识别系统有效空域的方法是判断当前时刻的我飞机是否存在于构成综合识别系统的所有传感器当前时刻探测侦察空域内,其方法是根据我作战飞机所在的轨迹的地理坐标独自转化为每个构成敌我识别系统的每个传感器坐标系中我飞机对每个敌传感器的相对位置,即可判断该时刻我飞机是否落入对方综合识别系统有效空域。
我机在时间t(t∈T)处于对方综合识别系统有效空域的概率为Pra,并假设敌构成综合识别系统传感器互相独立,互不影响,则定义
定义为我作战飞机落在Xi内的概率,
其中,
是我飞机在对方的构成综合识别系统第i个传感器真实相对位置,是对方的目标防空探测区s的第i个探测设备在t时刻的对我飞机方位角和俯仰角的估计值。此时,Xi(i=1,…,m)为对方第i个探测设备的为对方第i个探测设备的方位覆盖范围,为对方第i个探测设备的俯仰覆盖范围。
S32:飞机被对方识别系统判别为对方威胁的概率Ptt。假设对方防空区域内有n个传感器,包括雷达、敌我识别器、雷达信号侦察传感器、通信信号侦察传感器等,对我飞机的敌我属性进行识别,设定敌我识别分别为我方、对方,有K种敌我属性Hj(j=1,2,…,K);对我飞机的机型属性进行识别,有K种机型Hj(j=1,2,…,K),如小型机、大型机、民航机等;对我飞机的类型属性进行识别,有K种类型Hj(j=1,2,…,K),如直升机、轰炸机、歼击机、预警机、电子干扰机等;对我作战飞机的型号属性进行识别,有K种型号Hj(j=1,2,…,K),如F35、F22、E2C、E2D、Su是27、B是52等。先验分布概率Pr(x/Hj),x为敌、我(j=1,2,…,K),作为证据n个传感器证据的目标特征信息为Ei(i=1,2,…,n),这些特征信息来自相互独立传感器。
(1)POSM计算对我飞机机型估计的后验概率
j=1,2,…,K;Pr(Hj/E1,E2,…,En)为n个传感器观测发生情况下,目标属于属性类Hj概率;当n个传感器观测证据互相独立时,即E1,E2,…,En相互独立时,
P((E1,E2,…,En)/Hj)=P(E1/Hj)×P(E2/Hj)×…×P(En/Hj);
P(Ei/Hj),i=1.2,…,n表示第i个传感器观测证据在固定的机型参数Hj下的似然函数;
Pr(Hj)是假设Hj的先验概率,则采用极大后验判定逻辑,选择
(2)计算对我作战飞机敌我识别的后验概率
(3)POSM根据类型属性与底层机型的隶属关系,若给出机型Hj为真条件下类型属性为Ai的先验概率,则可以计算出该机类型属性的后验估计概率:
采用极大后验判定逻辑,选择
(4)根据类型属性与底层机型的隶属关系,若给出类型Ai为真条件下型号属性为Bm的先验概率,则可以计算出该机型号属性的后验估计概率:
采用极大后验判定逻辑,选择;
(5)飞机被对方识别系统判别为对方威胁的概率Ptt。如果仅给出我方飞机的机型属性Hj的先验概率Pr(Hj)和Pr(我(或敌)/Hj),则Ptt=Pr(敌/E1,E2,…,En);如果给出我方飞机的机型属性Hj的先验概率Pr(Hj)、Pr(我(或敌)/Hj)和Pr(Ai/Hj),则Ptt=Pr(A/E1,E2,…,En);如果给出我方飞机的机型属性Hj的先验概率Pr(Hj)、Pr(我(或敌)/Hj)、Pr(Ai/Hj)和Pr(Bm/Ai),则Ptt=Pr(B/E1,E2,…,En);
(6)飞机被对方识别系统判别为对方威胁的概率Ptt。POSM利用Bayes算法,由观测得到的条件概率和指定的先验概率,得到对我方飞机的后验概率;然后在新观测的下一次循环中用后验概率更新先验概率,并最终得到多次循环后收敛的对我方的识别概率。
S4:POSM计算飞机被对方跟踪的概率Pt。根据飞机进入对方跟踪系统作用空域的概率Pta,在对方跟踪系统作用空域内飞机被捕获的概率Pct和对方跟踪系统作用空域内飞机摆脱跟踪的概率Pst,计算飞机被对方跟踪的概率Pt=Pta·Pct·(1-Pst)。
S41:POSM计算飞机进入对方跟踪系统作用空域的概率Pta。POSM判断飞机是否落入对方跟踪系统有效空域的方法是判断当前时刻的作战飞机是否存在于构成对方跟踪系统的所有传感器作用空域内,其方法是根据飞机所在的轨迹的地理坐标独自转化为每个构成对方跟踪系统的每个传感器坐标系中飞机对对方每个传感器的相对位置,即可判断该时刻我飞机是否落入对方跟踪系统有效空域。
我机在时间t(t∈T)处于对方跟踪系统有效空域的概率为Pra,并假设对方跟踪系统传感器互相独立,互不影响,则定义
为我作战飞机落在Xi内的概率,
其中,
是我飞机在构成对方跟踪系统第i个传感器真实相对位置,是作战飞机在对方引导传感器引导下以及在对方导航系统作用下在t时刻的对我飞机方位角和俯仰角的估计值。此时,Xi(i=1,…,m)为对方第i个探测设备的[αi min,αi max]为对方第i个探测设备的方位覆盖范围,为对方第i个探测设备的俯仰覆盖范围。
S42:在对方跟踪系统作用空域内POSM计算飞机被捕获的概率Pct。假设虚警观测与实际观测无关,Z为多传感器对多目标观测集,W为虚警或杂波观测集,X为多目标状态集,观测形式为Z=T∪W,T=T(X)是实际目标产生的观测,f(Z|X)是在整个传感器组合获取的所有目标数据的似然函数,根据递归最佳多目标跟踪的Bayesian递归非线性滤波方程为:
式中,k表示第k时刻,Xk表示k时刻多个对方目标的状态集,Zk表示k时刻多传感器对多目标的观测集,Z1:k是在时-阶k时的累积证据集。
假设对我方目标的初始状态先验分布f(X0)已知,经过非线性滤波递推可以获得关于当前状态分布f(Xk|Z1:k)的完整信息,则飞机被捕获的概率,
S43:在对方跟踪系统作用空域内,飞机摆脱跟踪的概率Pst,作战飞机为了提髙自身的生存能力,通常在被对方雷达锁定情况下会作规避的机动动作、释放干扰力图摆脱跟踪,前者主要通过自身运动状态的快速变化导致雷达跟踪器精度变差甚至丢失跟踪目标,后者则通过制造假目标掩护自身,即为雷达探测到目标的虚假位置或者由雷达探测到的数据无法通过跟踪模型对其进行精确跟踪。以作战飞机机动逃逸建立飞机摆脱跟踪的概率Pst为例。如果作战飞机要逃逸雷达跟踪,则需要作战飞机在一个对方雷达扫描周期内,机动出对方雷达扫面波束范围,假设作战飞机在下一个雷达扫描周期时刻的位置为点B,则点B与雷达波束中轴线的最短距离为d,与对方雷达波束中轴线垂直且经过点B的法平面所截的横截面的扫描半径为r,为了摆脱雷达的跟踪,则必须满足的条件为d>r,d-r差值越大,下一时刻雷达重新捕获我作战飞机的难度越大,我机越容易逃逸。逃逸成功概率与我机加速度大小、逃逸时刻的速度大小和方向有关系,逃逸模型的最佳选取策略选择逃逸加速度方向与逃逸时刻的速度方向之间角尽可能小的情况下d-r尽可能大,在此策略下我作战飞机摆脱跟踪的概率为
k为修正因子。
S5:飞机被对方防空武器击毁的概率Pk,地面火力防区一般是整个防空体系中各类防空武器对各类目标的综合作战空域。由于地面防空火控系统要按照空中目标所处位置的不同,需要控制不同的武器来杀伤目标。当发现目标后能够迅速对目标实施有效的射击,射出的弹丸能够对目标造成致命的杀伤。毁伤概率是在给定射击条件下火控系统击毁目标可能性度量。
我作战飞机被对方毁伤概率
式中:n是对方防空武器发射武器总数目,Pka(i)是我作战飞机进入对方第i种防空武器作用空域的概率,Pht(i)是我作战飞机被对方第i种防空武器击中的概率,Pmr(i)是我作战飞机被对方第i种防空武器击中条件下的损毁概率。
S51:我作战飞机进入对方第i种防空武器杀伤空域的概率Pka(i),我飞机进入对方第i种防空武器杀伤空域的方法是根据我飞机所在的轨迹的地理坐标独自转化为每个防空武器坐标系中我飞机对每个防空武器的相对位置,即可判断该时刻我飞机是否落入防空武器系统有效空域。这里简化计算为:Pka(i)=第i种武器攻击区角度范围/360。
S52:POSM计算击中概率与击中条件下杀伤概率的乘积Pht(i)·Pmr(i),其杀伤概率算可用如下表达式来进行解算:Pht(i)·Pmr(i)=1-(1-P1/r)nk
式中,P1为防空导弹单发杀伤概率,n为每次攻击投射的弹数,k为敌发动攻击波次数,r为平均必需命中弹数,是指摧毁某个目标所需弹数的平均值.各类目标受到不同武器攻击时的平均必需命中弹数不同,各种武器对不同目标的平均必需命中弹数可由靶场实验数据得到。
S521:POSM计算防空导弹单发杀伤概率P1,P1=PK/H×PH,并且
式中:PH代表敏感性,对目标的杀伤度,是指导弹在完成任务过程中命中目标的攻击效,果,RK为导弹战斗部杀伤半径,Rδ为圆概率误差,RT为脱靶量,PK/H代表易损性,是飞机被击中后的杀伤概率,是一种条件概率。地空导弹战斗部大多数属于高爆破片式战斗部,其击中目标机后,破片对目标的杀伤概率可以表示为PK/H=1-exp(-Nv×Pai)
式中:Nv为落入易损面积上的破片数,Pai为单枚破片的杀伤概率。
S6:POSM计算飞机的生存概率P。作战飞机在t0时刻生存概率的计算公式为:
P=(1-Pf)·(1-Pr·Pi·Pt·Pk)
式中:P表示飞机生存概率;Pf表示飞机在执行任务过程中出现致命故障的概率;Pr表示飞机被对方发现的概率;Pi表示飞机被对方发现条件下被对方识别的概率;Pt表示飞机被对方发现条件下被对方跟踪的概率;Pk表示飞机被对方击毁的概率。
尽管上述已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于包括如下步骤,根据地面防空探测区域为防空探测系统每个探测单元同时探测区域Xi(i=1,…,m)的并集(X1∪X2∪…∪Xm)构成的包络区域,判断飞机是否落入对方探测系统有效搜索空域,判断当前时刻作战飞机是否存在包络空域区域内,将作战飞机所在轨迹的地理坐标分别转入对方每个地面防空探测系统的探测坐标系,解出作战飞机对探测装置的相对位置,判断出当前时刻作战飞机是否落入对方探测系统有效搜索空域,确定飞机进入对方探测系统有效搜索空域的概率;在飞机在执行任务过程中影响生存概率因素的基础上,生存概率模型子系统POSM计算出飞机在执行对地攻击任务过程中出现故障的概率、被发现概率、被识别概率、被跟踪的概率和被毁伤概率,计算这些阶段事件出现概率的乘积和作战活动的每个阶段的每个时刻生存概率,给出飞机生存概率的计算过程,得到飞机生存概率,然后依据飞机生存概率,建立飞机平台系统生存能力概率模型。
2.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:在计算生存概率中,首先计算得到被发现概率、被识别的概率、被跟踪的概率和被毁伤概率,然后计算得到飞机在当前时刻的生存概率;使用飞机被敌方防空探测系统发现概率公式Pr=Pea·Pei,计算我作战飞机被对方发现概率Pr,使用飞机被对方识别概率公式Pi=Pra·Ptt计算我作战飞机被对方识别概率Pi;使用飞机被对方跟踪概率公式Pt=Pta·Pct·(1-Pst)计算我作战飞机被对方跟踪概率Pt;使用飞机被毁伤概率计算公式计算我作战飞机被对方毁伤概率Pk;
飞机在t0时刻生存概率的计算公式为:P=(1-Pf)·(1-Pr·Pi·Pt·Pk)
其中,Pea是飞机进入敌方探表示飞机进入对方识别系统作用空域的概率,Ptt表示飞机被对方识别系统判别为对方威胁的概率测系统有效搜索空域的概率,Pei是飞机辐Pra射率;或散射的电磁波被敌方探测系统截获的概率或敌方红外探测系统对作战飞机红外辐射探测概,n是对方防空武器发射武器总数目,Pka(i)是我作战飞机进入对方第i种防空武器作用空域的概率,Pht(i)是我作战飞机被对方第i种防空武器击中的概率,Pmr(i)是我作战飞机被对方第i种防空武器击中条件下的损毁概率;P表示飞机生存概率,Pf表示飞机在执行任务过程中出现致命故障的概率,Pr表示飞机被对方发现的概率,Pi表示飞机被对方发现条件下被对方识别的概率,Pt表示飞机被对方发现条件下被对方跟踪的概率,Pk表示飞机被对方击毁的概率。
3.如权利要求2所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:生存概率模型子系统POSM计算飞机故障概率Pf,
并且
其中,Pf表示飞机故障概率,MTBFf表示飞机在执行任务的时候出现故障的平均间隔时间,表示系统处于正常工作状态的数量特征;MTTRf表示飞机在执行任务的时候故障平均修复时间,表征系统出于故障状态的数量特征,(MTTRf)i表示构成飞机的第i个可更换的部件LRU完成一次修复性维修所需要的时间;λi是构成飞机的第i个可更换的部件LRU的故障率;n表示构成飞机的第i个可更换的部件LRU的总数。
4.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:POSM计算飞机进入对方探测系统有效搜索空域的概率Pea,作战飞机在预警机数据链导引下飞往对方的防空探测区域s,在时刻t(t∈T)进入目标防空探测区s的概率为Pea,则
式中,P(X1∪X2∪…∪Xm)为飞机落入防空探测系统每个探测单元在同一时刻的探测区域Xi(i=1,…,m)的并集(X1∪X2∪…∪Xm)构成的包络区域的概率。
5.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:飞机在预警机数据链导引下飞往对方的目标防空探测区s的第i个探测设备的探测区域Xi中出现的联合概率密度,假设对方探测设备在二维空域进行搜索,即在方位角和俯仰角范围内搜索,则
其中,是飞机在目标防空探测区s的第i个探测设备的探测坐标系的真实相对位置,是对方的目标防空探测区s的第i个探测设备在t时刻对我飞机方位角和俯仰角的估计值,为对方的目标防空探测区s的第i个探测设备对我飞机探测方位角和俯仰角的误差的标准方差。
6.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:POSM根据对方雷达对我方飞机雷达散射截面积RCS截获概率P(Z|Pr≥Smin,H1)和对我方飞机的检测概率Pd,计算我机散射的电磁波被对方探测系统截获概率P1,且P1=P(Pr≥Smin|H1)×P(Z|Pr≥Smin,H1)≈Pd·P(Z|Pr≥Smin,H1),其中,Pr为对方雷达发射信号被我飞机反射后到达对方截获接收机的信号功率,Smin为对方接收机识别信号的灵敏度门限,H1表示事件“雷达发射波束对准目标”,Z事件表示“雷达截获到我方飞机RCS波瓣σRCS(θ)≠0”事件。
7.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:POSM根据第k台接收机发生一次截获的截获概率Pi_k=Pd·Ps·Pt·Pf计算辐射电磁波被对方探测系统截获的概率P2,
其中,n表示辐射波束区域内存在n部接收器;Pi_cum为整个波束覆盖区内发生一次截获的截获概率;Pd表示对方接收机的检测概率,即到达对方接收机的信号功率超过敌接收机最小可检测的门限功率时,对方接收机探测到我发射机波束能量的概率;Ps表示单个对方接收机的空域截获概率,即我发射机波束照射到对方飞机所在位置的概率;Pt表示单个对方接收机的时间截获概率,表示我方发射机设备正在辐射信号的时间内,对方接收设备处于开机状态的概率;Pf表示单个对方接收机的频率截获概率,即我发射机的当前工作频率正好落入对方接收机瞬时测频带宽,且满足对信号的测频条件的概率。
8.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:POSM将对方接收机的检测概率Pd、对方接收机的空域截获概率Ps、对方接收机的时间截获概率Pt和对方接收机的频率截获概率Pf综合起来,得到第k台接收机发生一次截获的截获概率计算公式为:
式中,MF表示主瓣覆盖面积(3dB);DI表示截获接收机密度;TOT表示我辐射源照射对方接收机时间;TI表示截获接收机的扫描时间。
9.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:飞机被对方识别系统判别为对方威胁的概率Ptt,如果仅给出我方飞机的机型属性Hj的先验概率Pr(Hj)和Pr(我(或敌)/Hj),则如果给出我方飞机的机型属性Hj的先验概率Pr(Hj)、Pr(我(或敌)/Hj)和Pr(Ai/Hj),则如果给出我方飞机的机型属性Hj的先验概率Pr(Hj)、Pr(我(或敌)/Hj)、Pr(Ai/Hj)和Pr(Bm/Ai),则
10.如权利要求1所述的建立飞机平台系统对地攻击生存能力概率模型的方法,其特征在于:
POSM在对方跟踪系统作用空域内,,假设虚警观测与实际观测无关,同时假设对我方目标的初始状态先验分布f(X0)已知,则根据多传感器对多目标观测集Z和经过非线性滤波递推获得关于当前状态分布f(Xk|Z1:k)的完整信息,Z为多传感器对多目标观测集,和递归最佳多目标跟踪的Bayesian递归非线性滤波方程:
计算出飞机被捕获的概率其中,X为多目标状态集,k表示第k时刻,Xk表示k时刻多个对方目标的状态集,Zk表示k时刻多传感器对多目标的观测集,Z1:k是在时-阶k时的累积证据集,f(Z|X)是在整个传感器组合获取的所有目标数据的似然函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910580684.XA CN110414061B (zh) | 2019-06-29 | 2019-06-29 | 建立飞机平台系统对地攻击生存能力概率模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910580684.XA CN110414061B (zh) | 2019-06-29 | 2019-06-29 | 建立飞机平台系统对地攻击生存能力概率模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110414061A true CN110414061A (zh) | 2019-11-05 |
CN110414061B CN110414061B (zh) | 2022-10-28 |
Family
ID=68360028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910580684.XA Active CN110414061B (zh) | 2019-06-29 | 2019-06-29 | 建立飞机平台系统对地攻击生存能力概率模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110414061B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112464451A (zh) * | 2020-11-16 | 2021-03-09 | 中国人民解放军海军工程大学 | 基于作战仿真系统的防空导弹武器命中概率修正方法 |
CN113190989A (zh) * | 2021-04-26 | 2021-07-30 | 中国人民解放军海军大连舰艇学院 | 复杂电磁环境下复合制导舰空导弹截获目标概率仿真方法 |
CN113255118A (zh) * | 2021-05-11 | 2021-08-13 | 上海机电工程研究所 | 基于杀伤链的武器装备体系优化方法和系统 |
CN113852972A (zh) * | 2021-09-13 | 2021-12-28 | 金华航大北斗应用技术有限公司 | 一种基于波束共享的高速移动终端波束调度方法 |
CN114638298A (zh) * | 2022-03-15 | 2022-06-17 | 天津大学 | 飞行器攻击行为预测方法及电子设备 |
CN115051766A (zh) * | 2022-06-02 | 2022-09-13 | 中国电子科技集团公司第十研究所 | 一种通信侦察检测设备及检测方法 |
CN115858403A (zh) * | 2023-03-01 | 2023-03-28 | 中国电子科技集团公司第十研究所 | 一种电子系统的虚警率预计方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122759A1 (en) * | 1994-03-15 | 2004-06-24 | Spivey Loehr H. | Computer-implemented process and mechanism for implementing an employee stock ownership plan |
CN101136081A (zh) * | 2007-09-13 | 2008-03-05 | 北京航空航天大学 | 基于蚁群智能的无人作战飞机多机协同任务分配方法 |
CN102323963A (zh) * | 2011-06-16 | 2012-01-18 | 北京航空航天大学 | 一种飞行器动态综合隐身性能评估软件仿真平台及构建方法 |
CN102902977A (zh) * | 2012-06-21 | 2013-01-30 | 中国人民解放军电子工程学院 | 基于风场扰动特征的空中目标分类方法 |
US9104915B2 (en) * | 2008-08-19 | 2015-08-11 | Digimarc Corporation | Methods and systems for content processing |
US9462801B2 (en) * | 2006-07-05 | 2016-10-11 | Nipro Corporation | Method for preserving pancreatic islet, container for preserving pancreatic islet, and kit for transplanting pancreatic islet |
CN108318032A (zh) * | 2017-12-21 | 2018-07-24 | 上海交通大学 | 一种考虑攻防对抗的无人机航迹智能规划方法 |
CN108509728A (zh) * | 2018-04-04 | 2018-09-07 | 西北工业大学 | 一种对海协同作战飞机生存力权衡设计方法 |
CN208798093U (zh) * | 2018-09-21 | 2019-04-26 | 湖北蔚蓝国际航空学校股份有限公司 | 一种低空空域飞行器监视装置 |
-
2019
- 2019-06-29 CN CN201910580684.XA patent/CN110414061B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040122759A1 (en) * | 1994-03-15 | 2004-06-24 | Spivey Loehr H. | Computer-implemented process and mechanism for implementing an employee stock ownership plan |
US9462801B2 (en) * | 2006-07-05 | 2016-10-11 | Nipro Corporation | Method for preserving pancreatic islet, container for preserving pancreatic islet, and kit for transplanting pancreatic islet |
CN101136081A (zh) * | 2007-09-13 | 2008-03-05 | 北京航空航天大学 | 基于蚁群智能的无人作战飞机多机协同任务分配方法 |
US9104915B2 (en) * | 2008-08-19 | 2015-08-11 | Digimarc Corporation | Methods and systems for content processing |
CN102323963A (zh) * | 2011-06-16 | 2012-01-18 | 北京航空航天大学 | 一种飞行器动态综合隐身性能评估软件仿真平台及构建方法 |
CN102902977A (zh) * | 2012-06-21 | 2013-01-30 | 中国人民解放军电子工程学院 | 基于风场扰动特征的空中目标分类方法 |
CN108318032A (zh) * | 2017-12-21 | 2018-07-24 | 上海交通大学 | 一种考虑攻防对抗的无人机航迹智能规划方法 |
CN108509728A (zh) * | 2018-04-04 | 2018-09-07 | 西北工业大学 | 一种对海协同作战飞机生存力权衡设计方法 |
CN208798093U (zh) * | 2018-09-21 | 2019-04-26 | 湖北蔚蓝国际航空学校股份有限公司 | 一种低空空域飞行器监视装置 |
Non-Patent Citations (4)
Title |
---|
ANDRES SAN ANTONIO 等: ""Using simulation to estimate critical paths and survival functions in aircraft turnaround processes"", 《 2017 WINTER SIMULATION CONFERENCE》, 8 January 2018 (2018-01-08), pages 3394 - 3403 * |
宋文彬 等: ""基于四面体六棱长的三站测时差立体定位方法"", 《火力与指挥控制》, vol. 42, no. 6, 15 June 2017 (2017-06-15), pages 150 - 153 * |
李琨: ""国外ADS-B应用现状及趋势研究"", 《飞航导弹》, no. 12, 15 December 2018 (2018-12-15), pages 60 - 66 * |
谭勇 等: ""临近空间飞艇与防空识别区隐身飞机探测"", 《现代军事》, no. 7, 5 July 2015 (2015-07-05), pages 56 - 59 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112464451A (zh) * | 2020-11-16 | 2021-03-09 | 中国人民解放军海军工程大学 | 基于作战仿真系统的防空导弹武器命中概率修正方法 |
CN112464451B (zh) * | 2020-11-16 | 2021-08-13 | 中国人民解放军海军工程大学 | 基于作战仿真系统的防空导弹武器命中概率修正方法 |
CN113190989A (zh) * | 2021-04-26 | 2021-07-30 | 中国人民解放军海军大连舰艇学院 | 复杂电磁环境下复合制导舰空导弹截获目标概率仿真方法 |
CN113190989B (zh) * | 2021-04-26 | 2024-03-12 | 中国人民解放军海军大连舰艇学院 | 复杂电磁环境下复合制导舰空导弹截获目标概率仿真方法 |
CN113255118A (zh) * | 2021-05-11 | 2021-08-13 | 上海机电工程研究所 | 基于杀伤链的武器装备体系优化方法和系统 |
CN113255118B (zh) * | 2021-05-11 | 2023-05-30 | 上海机电工程研究所 | 基于杀伤链的武器装备体系优化方法和系统 |
CN113852972A (zh) * | 2021-09-13 | 2021-12-28 | 金华航大北斗应用技术有限公司 | 一种基于波束共享的高速移动终端波束调度方法 |
CN113852972B (zh) * | 2021-09-13 | 2023-10-10 | 金华航大北斗应用技术有限公司 | 一种基于波束共享的高速移动终端波束调度方法 |
CN114638298A (zh) * | 2022-03-15 | 2022-06-17 | 天津大学 | 飞行器攻击行为预测方法及电子设备 |
CN115051766A (zh) * | 2022-06-02 | 2022-09-13 | 中国电子科技集团公司第十研究所 | 一种通信侦察检测设备及检测方法 |
CN115858403A (zh) * | 2023-03-01 | 2023-03-28 | 中国电子科技集团公司第十研究所 | 一种电子系统的虚警率预计方法 |
CN115858403B (zh) * | 2023-03-01 | 2023-06-02 | 中国电子科技集团公司第十研究所 | 一种电子系统的虚警率预计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110414061B (zh) | 2022-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110414061A (zh) | 建立飞机平台系统对地攻击生存能力概率模型的方法 | |
Kang et al. | Protect your sky: A survey of counter unmanned aerial vehicle systems | |
Neri | Introduction to electronic defense systems | |
Sharma et al. | Artificial intelligence aided electronic warfare systems-recent trends and evolving applications | |
CN109740876A (zh) | 目标威胁判断方法 | |
Tyurin et al. | General approach to counter unmanned aerial vehicles | |
KR101378887B1 (ko) | 피아 식별 장치 및 시스템 | |
CN107300691B (zh) | 一种针对反舰导弹末制导雷达的分布式无源协同干扰方法 | |
CN109981212A (zh) | 一种基于数据链检测及反制的低慢小防控系统及方法 | |
CN108507577B (zh) | 一种基于数据链和飞机敏感性的任务规划方法 | |
CN105426634B (zh) | 一种预警机战场生存力评估方法 | |
Zohuri et al. | Electronic countermeasure and electronic counter-countermeasure | |
Popescu | The existing technologies on anti-drone systems | |
MacDonald | Hide and Seek: Remote Sensing and Strategic Stability | |
Lehto et al. | Small Drones’ Swarms and Military Paradigm Change | |
CN110261835A (zh) | 基于探测效率最大化的机载雷达协同探测工作方法 | |
Yildirim | Self-defense of large aircraft | |
CN109724472A (zh) | 一种车载要地区域综合防抗系统及方法 | |
Jiang et al. | Introduction to Network Radar Countermeasure Systems | |
Tang et al. | Evaluation of Radar Jamming Effect | |
Kostis et al. | Stealth aircraft tactical assessment using stealth entropy and digital steganography | |
Gonzalez-Jorge et al. | Counter Drone Technology: A Review | |
Jiang et al. | Four Countermeasure Capacity Analysis of Network Radar Countermeasure Systems | |
Yildiz et al. | Electronic attack and sensor fusion techniques for boot-phase defense against multiple ballistic threat missiles | |
Bukovetskiy et al. | Modeling of processes of security assurance from threats caused by unmanned air vehicles |
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 |