CN107292014A - 一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 - Google Patents
一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 Download PDFInfo
- Publication number
- CN107292014A CN107292014A CN201710450604.XA CN201710450604A CN107292014A CN 107292014 A CN107292014 A CN 107292014A CN 201710450604 A CN201710450604 A CN 201710450604A CN 107292014 A CN107292014 A CN 107292014A
- Authority
- CN
- China
- Prior art keywords
- guiding
- radar
- msub
- target
- mrow
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种定量评估雷达站引导能力的方法,属于预警探测和军事运筹交叉领域。基于雷达站引导运作流程,根据排队论理论,将雷达站视为一个随机服务系统。通过构建空袭兵器模型、己方兵力模型和动态引导模型等,结合计算机模拟技术,采用事件驱动模拟方法,通过蒙特卡洛试验,定量分析、评估雷达站的引导能力,解决了雷达站动态引导过程缺乏量化评估的难题,有助于有针对性地加强和完善雷达站建设,优化引导资源配置,提升空防作战能力。
Description
技术领域
本发明属于预警探测和军事运筹交叉领域,尤其涉及一种定量评估雷达站引导能力的方法。
背景技术
担负引导任务的航空兵地面雷达站是空防作战指挥的核心组成要素,其引导能力直接关系空防作战的成效。目前,国内对雷达站引导能力的评估多是针对平时的演习演练,设置的空中目标数量非常有限,评估的结论主要是针对既定的、已知的任务而言,定性的给出雷达站的引导能力是否满足航空兵实施拦截的保障要求。
战时,在各种复杂空情和复杂电磁环境条件下,尤其是空袭兵器和电磁环境具有较大的不确定性,我担负引导任务的雷达站的“能力底数”到底是多少、我方的引导方案该如何拟制等问题是空防作战准备急待研究和急需解决的重难点问题,具有更强的现实意义。
发明内容
本发明的目的是提供一种定量评估雷达站引导能力的方法,它充分考虑了雷达站引导运作流程的特点,尤其针对战时不确定的空袭兵器和电磁环境,通过蒙特卡洛模拟试验和统计分析,为量化评估雷达站的引导能力提供了一种有效、可行的方法。本发明的基于蒙特卡洛模拟的定量评估雷达站引导能力的方法主要包括空袭兵力模拟、己方兵力模拟、雷达站动态引导过程建模和蒙特卡洛模拟试验与统计分析等四个部分。
本发明的技术方案是提供了一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:
步骤1、进行空袭兵力模拟,即对空袭兵器和为空袭兵器提供掩护的远距离电子干扰构建模型,具体包括对空袭兵器的飞行速度、到达责任空域时间、在责任空域滞留时间以及远距离电子干扰对我方雷达探测能力的压制影响的参数进行数学描述;
步骤2、进行己方兵力模拟,即完成己方引导能力建模和己方拦截兵力建模,具体包括对己方引导站雷达对空袭兵器的最远探测距离、干扰压制系数、引导时间以及己方担负拦截任务的歼击机数量、歼击机毁伤目标概率的建模;
步骤3、进行雷达站动态引导过程建模,即确定雷达站动态引导模拟的时间逻辑关系,按照事件发生的先后顺序,运用事件驱动方法模拟雷达站动态引导过程;
步骤4、完成蒙特卡洛模拟试验与统计分析,得到雷达站引导能力的量化评估指标结果。
本发明还提供了一种用于定量评估雷达站引导能力的雷达站动态引导方法,其特征在于:
雷达站按照先到先引导规则依次实施指挥引导;若目标进入引导责任空域时,引导通道正在实施引导,目标排队等待;目标在责任空域滞留期间,若引导通道释放,则对其实施引导;若引导通道一直未释放,则不经引导飞出引导责任空域;
首先,根据包括空袭兵器的数量、飞行高度、飞行速度和目标平均到达率在内的参数,逐一批次模拟空袭兵器进入我方责任空域,计算当前批次空袭兵器在责任空域的滞留时间;
然后,模拟电子干扰压制,分别计算在有/无电子干扰条件下,我方雷达对空袭兵器的最大探测距离,得出干扰压制系数;
检查我方引导通道是否空闲,如果空闲,我方雷达开始引导,记录开始引导时刻,按照引导时间的随机分布,计算完成引导时刻;如果引导通道被占用,则计算引导通道释放时刻和等待引导时间,判断引导通道释放时空袭兵器是否已离开责任空域,若具备引导条件,模拟引导过程,否则进入下一批次的模拟;
引导完成后,判断当前引导是否为有效引导,如果是有效引导,一方面计算有效引导时间、有效引导保障距离和有效引导保障中断距离等指标结果;另一方面根据歼击机毁伤目标概率,模拟我方歼击机拦截空袭兵器;否则,进行下一批次的模拟;
若模拟的空袭兵器数量小于空袭兵器的总数量,则进行下一批次空袭目标进入责任空域、我方雷达实施引导以及歼击机拦截等事件的模拟;否则,结束当前模拟试验。
本发明的有益效果在于:本发明将引导雷达站视为一个随机服务系统,通过建立空袭兵器、己方兵力以及和动态引导等模型,结合计算机模拟技术,采用事件驱动模拟方法,通过蒙特卡洛试验,实现雷达站引导能力的量化评估,为引导站平时建设、战时引导方案的拟制以及作战决策提供了定量的参考依据。
附图说明
图1是引导模拟模型的时间逻辑关系图。
图2是引导过程动态模拟示意图。
具体实施方式
以下结合附图1-2对本发明的技术方案进行详细说明。
雷达站指挥引导的流程是,引导员利用雷达信息和引导指挥设施,对进入责任空域的空袭兵器,派出己方飞机,使之占据空中有利位置实施拦截。基于这一特点,根据排队论理论,可以将雷达站视为一个有限等待随机服务系统。其中,空袭兵器进入责任空域视为“顾客到达”;引导通道(包括引导员、雷达以及引导指挥设施等)视为“服务设施”;为空袭兵器分配引导通道视为“为顾客提供服务”;将己方飞机引导至有利拦截位置的时间视为“服务时间”;进入责任空域,但是未被引导的空袭兵器视为等待服务的“队列”。
下面结合附图1-2对本发明方法进一步说明。
本发明提供了一种定量评估雷达站引导能力的方法,其具体包括以下步骤:
步骤1、进行空袭兵力模拟,即对空袭兵器和为空袭兵器提供掩护的远距离电子干扰构建模型,具体包括对空袭兵器的飞行速度、到达责任空域时间、在责任空域滞留时间以及远距离电子干扰对我方雷达探测能力压制影响的参数进行数学描述。
空袭兵力模拟具体通过以下方式实现:
步骤1-1、描述空袭兵器的参数主要有:空袭持续时间Td、空袭目标总数Nr、目标高度h、目标速度v、目标平均到达率λ以及目标平均滞留时间τz等。根据经验假设:
空袭目标速度是一个随机值,服从均匀分布,均值为v;
空袭兵器到达责任空域是一个随机值,服从平均到达率为λ的泊松分布,λ=Nr/Td;
空袭兵器在责任空域滞留时间是一个随机值,服从平均滞留率为ε的指数分布,ε=1/τz。
步骤1-2、建立远距离电子干扰对我雷达探测能力压制模型,得到干扰压制系数;通过干扰压制系数,描述因远距离电子干扰而导致的引导延迟及引导中断的程度。
远距离电子干扰主要对雷达站引导的及时性和连续性产生影响,对其模拟具体表现在延迟引导和引导不连续两个方面。假设无电子干扰时,雷达最大探测距离为Rm;存在电子干扰时,雷达最大探测距离为Rj,则电子干扰导致雷达延迟引导的时间Δt为:
其中,v表示空袭目标速度。
设干扰压制系数则Δt也可表示为:
雷达引导保障距离设定为Ls:
Ls=vτs,
式中,τs为雷达站引导时间;
假设雷达因电子干扰而中断引导的距离,即雷达引导保障中断距离为Lb,它服从平均引导保障中断距离为La的指数分布,La=0.1(2-Cj)Ls。
步骤2、进行己方兵力模拟,包括己方引导能力建模和己方拦截兵力建模两个部分,具体包括己方引导站雷达对空袭兵器的最远探测距离、干扰压制系数、引导时间以及己方担负拦截任务的歼击机数量、歼击机毁伤目标概率等。
步骤2-1、己方引导能力建模是根据雷达装备参数和电子干扰参数,基于雷达方程,计算雷达对空袭兵器的最大探测距离和干扰压制系数。
雷达对空袭兵器的最大探测距离Rm:
其中,Pt为雷达的发射功率,Gt为发射天线增益,Gr为接收天线增益,λ为雷达工作波长,σ为目标的散射截面积,k是波耳兹曼常数,T0为接收机参考温度,Br为雷达工作带宽,F为噪声系数,Vmin为系统检测因子,Lr为雷达系统总的损耗。
有源干扰下,雷达最大探测距离Rj:
其中,Rj为干扰机相对雷达的距离,Lj为干扰的综合损耗,Pj为干扰机的发射功率,为干扰机天线在雷达方向上的增益,Gr(θ)为雷达天线在干扰机方向上的增益,rj为干扰信号对雷达天线的极化损失,S为雷达接收目标回波功率,J为干扰机进入雷达的功率。
干扰压制系数Cj:
雷达站引导时间为τs,假设它是一个随机值,服从平均引导率为μ的指数分布,μ=1/τs。
步骤2-2、建立己方拦截兵力建模;即针对己方担负拦截任务的歼击机构建模型,主要有歼击机数量、一架歼击机毁伤目标概率。其中,假设一架歼击机毁伤目标概率是一个随机值,服从均值为P1的均匀分布。
步骤3、进行雷达站动态引导过程建模,即确定雷达站动态引导模拟的时间逻辑关系,按照事件发生的先后顺序,运用事件驱动方法模拟雷达站动态引导过程。
雷达站引导过程建模具体通过以下方式实现:
步骤3-1、确定雷达站动态模拟模型的时间逻辑。
为了建立雷达站动态引导模拟模型,首先描述模型的时间逻辑,如附图1所示。其中:
(1)目标等待时间wq
当ts≤ta,目标等待时间wq=0;
当ta<ts<td,目标等待时间wq=ts-ta;
当ts≥td,目标等待时间wq=τz。
其中,ts表示引导通道空闲时刻及雷达站开始引导时刻,ta表示目标到达责任空域时刻。
(2)雷达站开始引导时刻ts
当ts≤ta,开始引导时刻ts=ta;
当ta<ts<td,开始引导时刻ts=ts;
当ts≥td,不具备开始引导时刻。
其中,td表示目标离开责任空域时刻。
(3)雷达站有效引导时间τs
当ts≤ta且tc≤td,视为有效引导,引导时间τs=tc-ta;
当ta<ts<td且tc≤td,视为有效引导,引导时间τs=tc-ts;
否则,有效引导时间τs=0。
其中,tc表示雷达站引导结束时刻。
(4)雷达站引导结束时刻tc
当ts≤ta且tc≤td,引导结束时刻tc=ta+τs;
当ta<ts<td且tc≤td,引导结束时刻tc=ts+τs;
当ts≥td,引导结束时刻等于上次引导结束时刻;
否则,引导结束时刻tc=td。
步骤3-2、建立雷达站动态引导模拟模型。运用事件驱动的方法模拟雷达站动态引导过程,如附图2所示。
将雷达站动态引导过程分解为空袭兵器进入、电子干扰压制、指挥引导、航空兵拦截等若干个关键事件,按照时间顺序,通过事件驱动,完成引导过程模拟。
雷达站按照先到先引导规则依次实施指挥引导;若目标进入引导责任空域时,引导通道正在实施引导,目标排队等待;目标在责任空域滞留期间,若引导通道释放,则对其实施引导;若引导通道一直未释放,则不经引导飞出引导责任空域。
首先,根据空袭兵器的数量、飞行高度、飞行速度和目标平均到达率等参数,逐一批次模拟空袭兵器进入我方责任空域,计算当前批次空袭兵器在责任空域的滞留时间。
然后,模拟电子干扰压制,分别计算在有/无电子干扰条件下,我方雷达对空袭兵器的最大探测距离,得出干扰压制系数。
检查我方引导通道是否空闲,如果空闲,我方雷达开始引导,记录开始引导时刻,按照引导时间的随机分布,计算完成引导时刻;如果引导通道被占用,则计算引导通道释放时刻和等待引导时间,判断引导通道释放时空袭兵器是否已离开责任空域,若具备引导条件,模拟引导过程,否则进入下一批次的模拟。
引导完成后,判断当前引导是否为有效引导,如果是有效引导,一方面计算有效引导时间、有效引导保障距离和有效引导保障中断距离等指标结果;另一方面根据歼击机毁伤目标概率,模拟我方歼击机拦截空袭兵器;否则,进行下一批次的模拟。
若模拟的空袭兵器数量小于空袭兵器的总数量,则进行下一批次空袭目标进入责任空域、我方雷达实施引导以及歼击机拦截等事件的模拟;否则,结束当前模拟试验。
步骤4、完成蒙特卡洛模拟试验与统计分析,得到雷达站引导能力的量化评估指标结果。
由于描述空袭兵器、电子干扰压制、己方兵力等的参数属于随机变量,具有不确定性,对雷达站动态引导过程的模拟需进行满足置信水平和精度条件要求的足够次数的蒙特卡洛模拟试验。对模拟试验结果进行统计分析,得到雷达站引导能力的量化评估指标结果。
每次模拟试验均按照已建立的动态引导模拟模型执行,模拟结束后获得大量的模拟统计数据,用于评估雷达站的引导能力。
雷达站引导能力的评估指标和计算方法如下:
(1)引导目标概率Pg
Pg=N/Nr
其中,N为有效引导目标数量,Nr为空袭目标总数。
(2)平均每分钟引导目标数量Sa
Sa=N/(to-ti)
其中,to为所有空袭兵器最晚离开时刻,ti为所有空袭兵器最早到达时刻。
(3)队中平均目标等待时间
其中,wqi为第i个空袭目标等待时间。
(4)歼击机拦截作战效能E
假定己方派给每个空袭目标的歼击机数量为n,歼击机拦截作战效能为:
E=PgPk
其中,Pg为引导目标概率;Pk为派出n架歼击机拦截毁伤目标的概率,且Pk=1-(1-Pl)n;Pl为一架歼击机毁伤空袭目标概率。
每次蒙特卡洛模拟试验结束后,按照雷达站引导能力评估指标和计算方法,可以得到一组描述雷达站引导能力的指标结果。对多次蒙特卡洛试验的模拟结果进行统计、分析、计算,便可得到雷达站引导能力的量化评估指标结果。
Claims (6)
1.一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:
步骤1、进行空袭兵力模拟,即对空袭兵器和为空袭兵器提供掩护的远距离电子干扰构建模型,具体包括对空袭兵器的飞行速度、到达责任空域时间、在责任空域滞留时间以及远距离电子干扰对我方雷达探测能力的压制影响的参数进行数学描述;
步骤2、进行己方兵力模拟,即完成己方引导能力建模和己方拦截兵力建模,具体包括对己方引导站雷达对空袭兵器的最远探测距离、干扰压制系数、引导时间以及己方担负拦截任务的歼击机数量、歼击机毁伤目标概率的建模;
步骤3、进行雷达站动态引导过程建模,即确定雷达站动态引导模拟的时间逻辑关系,按照事件发生的先后顺序,运用事件驱动方法模拟雷达站动态引导过程;
步骤4、完成蒙特卡洛模拟试验与统计分析,得到雷达站引导能力的量化评估指标结果。
2.根据权利要求1所述的基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:步骤1中,空袭兵力模拟具体通过以下方式实现:
步骤1-1、描述空袭兵器的参数包括:空袭持续时间Td、空袭目标总数Nr、目标高度h、目标速度v、目标平均到达率λ以及目标平均滞留时间τz;根据经验假设:
空袭目标速度是一个随机值,服从均匀分布,均值为v;
空袭兵器到达责任空域是一个随机值,服从平均到达率为λ的泊松分布,λ=Nr/Td;
空袭兵器在责任空域滞留时间是一个随机值,服从平均滞留率为ε的指数分布,ε=1/τz;
步骤1-2、建立远距离电子干扰对我雷达探测能力压制模型,得到干扰压制系数;通过干扰压制系数,描述因远距离电子干扰而导致的引导延迟及引导中断的程度;
远距离电子干扰主要对雷达站引导的及时性和连续性产生影响,对其模拟具体表现在延迟引导和引导不连续两个方面;假设无电子干扰时,雷达最大探测距离为Rm;存在电子干扰时,雷达最大探测距离为Rj,则电子干扰导致雷达延迟引导的时间Δt为:
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>j</mi>
</msub>
</mrow>
<mi>V</mi>
</mfrac>
<mo>,</mo>
</mrow>
其中,v表示空袭目标速度。
设干扰压制系数则Δt也可表示为:
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mi>V</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
雷达引导保障距离设定为Ls:
Ls=vτs,
式中,τs为雷达站引导时间;
假设雷达因电子干扰而中断引导的距离,即雷达引导保障中断距离为Lb,它服从平均引导保障中断距离为La的指数分布,La=0.1(2-Cj)Ls。
3.根据权利要求1所述的基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:步骤2中包括:
步骤2-1、己方引导能力建模是根据雷达装备参数和电子干扰参数,基于雷达方程,计算雷达对空袭兵器的最大探测距离和干扰压制系数;
雷达对空袭兵器的最大探测距离Rm:
<mrow>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mi>t</mi>
</msub>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<msub>
<mi>G</mi>
<mi>r</mi>
</msub>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<mi>&sigma;</mi>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<msub>
<mi>kT</mi>
<mn>0</mn>
</msub>
<msub>
<mi>B</mi>
<mi>r</mi>
</msub>
<msub>
<mi>FV</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mi>r</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>4</mn>
</mrow>
</msup>
</mrow>
其中,Pt为雷达的发射功率,Gt为发射天线增益,Gr为接收天线增益,λ为雷达工作波长,σ为目标的散射截面积,k是波耳兹曼常数,T0为接收机参考温度,Br为雷达工作带宽,F为噪声系数,Vmin为系统检测因子,Lr为雷达系统总的损耗;
有源干扰下,雷达最大探测距离Rj:
其中,Rj为干扰机相对雷达的距离,Lj为干扰的综合损耗,Pj为干扰机的发射功率,为干扰机天线在雷达方向上的增益,Gr(θ)为雷达天线在干扰机方向上的增益,rj为干扰信号对雷达天线的极化损失,S为雷达接收目标回波功率,J为干扰机进入雷达的功率;
干扰压制系数Cj:
<mrow>
<msub>
<mi>C</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>j</mi>
</msub>
<msub>
<mi>R</mi>
<mi>m</mi>
</msub>
</mfrac>
</mrow>
雷达站引导时间为τs,假设它是一个随机值,服从平均引导率为μ的指数分布,μ=1/τs;
步骤2-2、建立己方拦截兵力建模;即针对己方担负拦截任务的歼击机构建模型,包括歼击机数量、一架歼击机毁伤目标概率;其中,假设一架歼击机毁伤目标概率是一个随机值,服从均值为P1的均匀分布。
4.根据权利要求1所述的基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:步骤3中,雷达站引导过程建模通过以下方式实现:
步骤3-1、确定雷达站动态模拟模型的时间逻辑;
(1)目标等待时间wq
当ts≤ta,目标等待时间wq=0;
当ta<ts<td,目标等待时间wq=ts-ta;
当ts≥td,目标等待时间wq=τz;
其中,ts表示引导通道空闲时刻及雷达站开始引导时刻,ta表示目标到达责任空域时刻,τz是目标平均滞留时间;
(2)雷达站开始引导时刻ts
当ts≤ta,开始引导时刻ts=ta;
当ta<ts<td,开始引导时刻ts=ts;
当ts≥td,不具备开始引导时刻;
其中,td表示目标离开责任空域时刻;ta表示目标到达责任空域时刻;
(3)雷达站有效引导时间τs
当ts≤ta且tc≤td,视为有效引导,引导时间τs=tc-ta;
当ta<ts<td且tc≤td,视为有效引导,引导时间τs=tc-ts;
否则,有效引导时间τs=0;
其中,tc表示雷达站引导结束时刻;
(4)雷达站引导结束时刻tc
当ts≤ta且tc≤td,引导结束时刻tc=ta+τs;
当ta<ts<td且tc≤td,引导结束时刻tc=ts+τs;
当ts≥td,引导结束时刻等于上次引导结束时刻;
否则,引导结束时刻tc=td;
步骤3-2、建立雷达站动态引导模拟模型,将雷达站动态引导过程分解为包括空袭兵器进入、电子干扰压制、指挥引导、航空兵拦截在内的关键事件,按照时间顺序,通过事件驱动,完成引导过程模拟。
5.根据权利要求1所述的基于蒙特卡洛模拟的定量评估雷达站引导能力的方法,其特征在于:步骤4中,雷达站引导能力的评估指标和计算方法如下:
(1)引导目标概率Pg
Pg=N/Nr,其中,N为有效引导目标数量,Nr为空袭目标总数;
(2)平均每分钟引导目标数量Sa
Sa=N/(to-ti)
其中,to为所有空袭兵器最晚离开时刻,ti为所有空袭兵器最早到达时刻;
(3)队中平均目标等待时间
<mrow>
<msub>
<mover>
<mi>w</mi>
<mo>&OverBar;</mo>
</mover>
<mi>q</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>r</mi>
</msub>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
其中,wqi为第i个空袭目标等待时间;
(4)歼击机拦截作战效能E
假定己方派给每个空袭目标的歼击机数量为n,歼击机拦截作战效能为:
E=PgPk
其中,Pg为引导目标概率;Pk为派出n架歼击机拦截毁伤目标的概率,且Pk=1-(1-Pl)n;Pl为一架歼击机毁伤空袭目标概率;
每次蒙特卡洛模拟试验结束后,按照雷达站引导能力评估指标和计算方法,得到一组描述雷达站引导能力的指标结果。
6.一种用于定量评估雷达站引导能力的雷达站动态引导方法,其特征在于:
雷达站按照先到先引导规则依次实施指挥引导;若目标进入引导责任空域时,引导通道正在实施引导,目标排队等待;目标在责任空域滞留期间,若引导通道释放,则对其实施引导;若引导通道一直未释放,则不经引导飞出引导责任空域;
首先,根据包括空袭兵器的数量、飞行高度、飞行速度和目标平均到达率在内的参数,逐一批次模拟空袭兵器进入我方责任空域,计算当前批次空袭兵器在责任空域的滞留时间;
然后,模拟电子干扰压制,分别计算在有/无电子干扰条件下,我方雷达对空袭兵器的最大探测距离,得出干扰压制系数;
检查我方引导通道是否空闲,如果空闲,我方雷达开始引导,记录开始引导时刻,按照引导时间的随机分布,计算完成引导时刻;如果引导通道被占用,则计算引导通道释放时刻和等待引导时间,判断引导通道释放时空袭兵器是否已离开责任空域,若具备引导条件,模拟引导过程,否则进入下一批次的模拟;
引导完成后,判断当前引导是否为有效引导,如果是有效引导,一方面计算有效引导时间、有效引导保障距离和有效引导保障中断距离等指标结果;另一方面根据歼击机毁伤目标概率,模拟我方歼击机拦截空袭兵器;否则,进行下一批次的模拟;
若模拟的空袭兵器数量小于空袭兵器的总数量,则进行下一批次空袭目标进入责任空域、我方雷达实施引导以及歼击机拦截等事件的模拟;否则,结束当前模拟试验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710450604.XA CN107292014B (zh) | 2017-06-15 | 2017-06-15 | 一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710450604.XA CN107292014B (zh) | 2017-06-15 | 2017-06-15 | 一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292014A true CN107292014A (zh) | 2017-10-24 |
CN107292014B CN107292014B (zh) | 2018-03-30 |
Family
ID=60096594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710450604.XA Active CN107292014B (zh) | 2017-06-15 | 2017-06-15 | 一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292014B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109324527A (zh) * | 2018-09-04 | 2019-02-12 | 西北工业大学 | 基于空天地海信息网支撑的雷达资源管理方法 |
CN112834987A (zh) * | 2021-01-18 | 2021-05-25 | 成都老鹰信息技术有限公司 | 一种飞机擦尾毫米波测距系统 |
CN113641191A (zh) * | 2021-10-14 | 2021-11-12 | 中国人民解放军空军预警学院 | 预警机和干扰机协同作业的空域配置方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183754A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种鲁棒智能雷达海上目标检测系统及方法 |
CN102183751A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海上目标检测系统及方法 |
CN102183752A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种自适应雷达海杂波预报系统及方法 |
CN103870654A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | 基于并行矩量法与物理光学混合的电磁散射仿真方法 |
-
2017
- 2017-06-15 CN CN201710450604.XA patent/CN107292014B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183754A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种鲁棒智能雷达海上目标检测系统及方法 |
CN102183751A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种智能雷达海上目标检测系统及方法 |
CN102183752A (zh) * | 2011-03-03 | 2011-09-14 | 浙江大学 | 一种自适应雷达海杂波预报系统及方法 |
CN103870654A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | 基于并行矩量法与物理光学混合的电磁散射仿真方法 |
Non-Patent Citations (4)
Title |
---|
刘兴堂等: "精确制导武器与仿真技术研究", 《系统仿真学报》 * |
姜林等: "武装直升机拦截反辐射无人机过程建模及效能评估", 《指挥与控制学报》 * |
陈永光: "基于电子对抗环境的地空 导弹系统仿真技术研究", 《电子对抗技术》 * |
陈永光等: "综合电子战仿真系统中关键技术的研究", 《系统工程与电子技术》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109324527A (zh) * | 2018-09-04 | 2019-02-12 | 西北工业大学 | 基于空天地海信息网支撑的雷达资源管理方法 |
CN109324527B (zh) * | 2018-09-04 | 2021-04-02 | 西北工业大学 | 基于空天地海信息网支撑的雷达资源管理方法 |
CN112834987A (zh) * | 2021-01-18 | 2021-05-25 | 成都老鹰信息技术有限公司 | 一种飞机擦尾毫米波测距系统 |
CN113641191A (zh) * | 2021-10-14 | 2021-11-12 | 中国人民解放军空军预警学院 | 预警机和干扰机协同作业的空域配置方法及设备 |
CN113641191B (zh) * | 2021-10-14 | 2022-01-14 | 中国人民解放军空军预警学院 | 预警机和干扰机协同作业的空域配置方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107292014B (zh) | 2018-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Aircraft re-routing optimization and performance assessment under uncertainty | |
CN107292014B (zh) | 一种基于蒙特卡洛模拟的定量评估雷达站引导能力的方法 | |
Svyd et al. | Method of increasing the identification friend or foe systems information security | |
CN103093649B (zh) | 用于推断飞行器参数的方法和系统 | |
CN110380897A (zh) | 基于改进bp神经网络的网络安全态势感知模型和方法 | |
CN103345587B (zh) | Ads‑b监视数据与雷达航迹的直觉模糊关联方法、装置 | |
CN109214419B (zh) | 一种基于空间聚类的鸟击事件风险综合防控方法 | |
CN109975774A (zh) | 基于算子的雷达交接效能柔性评估建模方法 | |
CN111612673A (zh) | 一种无人驾驶飞行器对多要地威胁程度确认方法及系统 | |
CN111583083A (zh) | 低空飞行的非合作目标对要地威胁程度的确定方法及系统 | |
CN113570910B (zh) | 空中交通流量管理方法、装置和电子设备 | |
Rivas et al. | Uncertainty | |
Lee et al. | Evaluating alerting and guidance performance of a UAS detect-and-avoid system | |
CN112396298B (zh) | 一种无人直升机多机协同任务规划方法 | |
Nikoleris et al. | Effect of trajectory prediction and stochastic runway occupancy times on aircraft delays | |
Punitha | Design and evaluation of traffic delays in toll plaza using combination of queueing and simulation | |
Zhao et al. | Study on combat effectiveness of air defense missile weapon system based on queuing theory | |
Atkinson et al. | Spatial queueing analysis of an interdiction system to protect cities from a nuclear terrorist attack | |
Dmochowski et al. | Air traffic smoothness as a universal measure for air traffic quality assessment | |
Zheng et al. | Research on channel optimization of Ads‐B aviation target surveillance radar based on improved filtering algorithm | |
Tian et al. | Probabilistic demand prediction model for en-route sector | |
He et al. | Constrained multiple model bayesian filtering for target tracking in cluttered environment | |
Li et al. | Distributed interference optimization method of large-scale UAV based on tabu search artificial bee colony algorithm | |
CN110378388A (zh) | 飞行物识别方法、装置、计算机可读存储介质及终端设备 | |
Volkov et al. | The protection of critical infrastructure facilities from air strikes due to compatible use of various forces and means |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210429 Address after: No. 303, group 3, liangfengding village, Zhengxing street, Tianfu New District, Chengdu, Sichuan 610218 Patentee after: Chengdu Furui Kongtian Technology Co.,Ltd. Address before: 100085 No. 11, Anning Road, Haidian District, Beijing Patentee before: BEIHANG University |