一种基于D数和威胁性优先的弹炮结合火力分配方法

摘要

本发明基于D数，提供一种在威胁度优先条件下的弹炮结合火力分配的方法，涉及弹炮结合的火力分配优化技术领域。本发明对空中目标选取威胁因素根据威胁函数生成威胁值，对威胁程度建立梯形模糊数，根据各因素威胁值生成基本概率分配函数，利用邓熵生成各因素权重后加权平均，用D数实现数据融合，生成目标威胁性矩阵，并结合命中率矩阵、火力单元矩阵和火力效能矩阵，根据威胁性对火力进行分配。本发明提出的D数和梯形模糊数相结合的方法，计算简单；本发明用梯形模糊数对战场信息建模，解决了模糊信息的表示；本发明提出的火力效能矩阵，为火力分配提供了依据；本发明提出的威胁度优先的火力分配原则，对低空突袭作战模式效果良好。

权利要求

1.一种基于D数和威胁性优先的弹炮结合火力分配方法，其特征在于包括下述步骤：步骤一：对空中目标威胁因素评估的方法为：

输入n个所探测空中目标，记为F1,F2,...,Fn，5种威胁因素的Dij，其中i＝1,2,…,n,j＝

1,2,…,5，对于每一个目标，计算各个因素威胁值，所述威胁因素为速度v、航路捷径p、高度h、干扰r和目标类型t，所述威胁值用威胁评估函数计算，分别为y(r)＝r和y(t)＝t；

步骤二：建立威胁度评估因子的威胁程度模糊数的方法为：输入空中目标各因素威胁值，根据梯形模糊数建立模型生成威胁程度的基本概率分配函数m(Fjk)，所述威胁程度k分为低l、中低ml、中m、中高mh和高h五个等级，所述模糊数以梯形模糊数表示，对应隶属度函数分别为(0,0,0.1,0.2)、(0.1,0.2,0.3,0.4)、(0.3,0.4,

0.6,0.7)、(0.6,0.7,0.8,0.9)和(0.8,0.9,1,1)；所述梯形模糊数是给定论域U上的一个模糊集，是指对任何x∈U，都有一个数μ(x)∈[0,1]与之对应，μ(x)称为x对U的隶属度，μ称为x的隶属函数，步骤三：生成威胁因素权重W1,W2,…W5，权重生成方法为：输入基本概率分配函数，对于同一目标不同威胁因素的分配函数，使用邓熵计算威胁因素权重，所述生成权重公式为 对权重归一化，所述归一化方法为：

步骤四：将生成的基本概率分配函数用特征权重W1,W2,…W5加权平均得到D(k)方法为：对于一个空中目标，对应威胁因素权重与基本概率分配函数m(Fjk)相乘，加权方法公式为步骤五：使用D数融合方法为：将加权平均得到的D(k)与自身使用D数融合j-1次得到D(A),D数融合规则为：其中

D1(B)和D2(C)是两组待组合的BPA，A、B、C指的是威胁等级，AreaB指的是威胁等级B的梯形模糊数面积，AreaB∩C指威胁等级B和威胁等级C的梯形模糊数相交面积；

步骤六：生成威胁度的方法为：

将融合后结果使用公式:S＝D(l)×0.1+D(ml)×0.3+D(m)×0.5+D(mh)×0.7+D(h)×

0.9，得到目标威胁度；所有目标按照步骤一至步骤六生成威胁度矩阵；

步骤七：生成命中率矩阵的方法为：

计算当前位置导弹和火炮对目标的命中率，

导弹命中率公式为：

火炮命中率公式为：

其中，p为航路捷径；O为命中率；

步骤八：生成火力单元矩阵的方法为：

火力单元矩阵表示当前火力单元是否可以工作，当该火力单元不能发射，即视为不可工作；

火力单元矩阵需根据实际情况生成；

步骤九：生成打击效能矩阵的方法为：

将步骤六的威胁度矩阵、步骤七的命中率矩阵和步骤八所产生的火力单元矩阵对应相乘得到打击效能矩阵，方法公式为：E＝S×O×D；E代表打击效能；

步骤十：根据威胁性进行火力分配的方法为：

对威胁性最大的目标优先进行火力分配，并分配打击效能最大的火力单元，然后将该火力单元和该目标对应火力单元置0，选取下一个威胁性最大的目标按上述方法进行火力分配，当空中目标多于火力单元时，将火力单元全部分配后，对打击效能矩阵进行更新，再进行火力分配，依次迭代。

说明书

一种基于D数和威胁性优先的弹炮结合火力分配方法

技术领域

[0001] 本发明涉及弹炮结合的火力分配优化技术领域，特别涉及一种基于D数和威胁性优先的弹炮结合火力分配方法。

背景技术

[0002] 弹炮结合防空武器系统是一种高性能的中近程防空武器系统，它整合了高射炮和地空导弹，同时具备高炮反应速度快、火力密级度好、近距离毁伤概率大，以及防空导弹射程远、精度高、单发杀伤概率大的优点，是未来反空袭作战的重要方向。火力分配是弹炮结合的关键因素，对其充分发挥作战效能起着至关重要的作用。

[0003] 针对复杂战场态势和高度不确定性环境下的弹炮系统火力分配模型，可以更加合理地协调和控制各子武器系统的行为，实现火力转换时机的优选和火力分配方案的优化，从而有效提升弹炮系统的作战效能。

[0004] 信息融合技术是协同利用多源信息，以获得对事物或目标更客观、更本质认识的信息综合处理技术，是智能科学研究的关键技术之一。在诸多的融合模型和方法中，D-S证据理论算法是最为有效的算法之一。证据理论把概率论中的基本事件空间拓宽为基本事件的幂集，又称为辨识框架，在辨识框架上建立了基本概率分配函数(Basic Probability Assignment，BPA)。此外，证据理论还提供了一个Dempster组合规则，该规则可以在没有先验信息的情况下实现证据的融合。特别地，当BPA只在辨识框架的单子集命题上进行分配时，BPA就转换为概率论中的概率，而组合规则的融合结果与概率论中的Bayes公式相同。D数是一种最新的不确定信息建模方法，是DS证据理论的推广，具备比D-S证据更好的不确定性表示和处理能力，使其在信息融合、决策分析、分析评估、模式识别等领域得到了广泛的应用。

[0005] 而针对低空突袭的作战模式，突袭目标需要突破层层防线，应该优先对威胁性大的目标进行打击。

发明内容

[0006] 为了实现弹炮结合的火力分配，本发明基于D数和威胁性优先原则，提供一种火力分配方法。使用该方法实现的火力分配方案，在低空突防等防空作战任务中，有明显效果。

[0007] 本发明解决其技术问题所采用的技术方案包括如下步骤：

[0008] 步骤一：对空中目标威胁因素评估的方法为：

[0009] 输入n个所探测空中目标，记为F1,F2,...,Fn，5种威胁因素的Dij，其中i＝1,2,…,n,j＝1,2,…,5，对于每一个目标，计算各个因素威胁值，所述威胁因素为速度v、航路捷径p、高度h、干扰r和目标类型t，所述威胁值用威胁评估函数计算，分别为y(r)＝r和y(t)＝t；

[0010] 步骤二：建立威胁度评估因子的威胁程度模糊数的方法为：

[0011] 输入空中目标各因素威胁值，根据梯形模糊数建立模型生成威胁程度的基本概率分配函数m(Fjk)，所述威胁程度k分为低l、中低ml、中m、中高mh和高h五个等级，所述模糊数以梯形模糊数表示，对应隶属度函数分别为(0,0,0.1,0.2)、(0.1,0.2,0.3,0.4)、(0.3,0.4,0.6,0.7)、(0.6,0.7,0.8,0.9)和(0.8,0.9,1,1)；所述梯形模糊数是给定论域U上的一个模糊集，是指对任何x∈U，都有一个数μ(x)∈[0,1]与之对应，μ(x)称为x对U的隶属度，μ称为x的隶属函数，

[0012] 步骤三：生成威胁因素权重W1,W2,…W5，权重生成方法为：

[0013] 输入基本概率分配函数，对于同一目标不同威胁因素的分配函数，使用邓熵计算威胁因素权重，所述生成权重公式为 对权重归一化，所述归一化方法为：

[0014] 步骤四：将生成的基本概率分配函数用特征权重W1,W2,…W5加权平均得到D(k)方法为：

[0015] 对于一个空中目标，对应威胁因素权重与基本概率分配函数m(Fjk)相乘，加权方法公式为

[0016] 步骤五：使用D数融合方法为：

[0017] 将加权平均得到的D(k)与自身使用D数融合j-1次得到D(A),D数融合规则为：

[0018]

[0019] 其中

[0020]

[0021] D1(B)和D2(C)是两组待组合的BPA，A、B、C指的是威胁等级，AreaB指的是威胁等级B的梯形模糊数面积，AreaB∩C指威胁等级B和威胁等级C的梯形模糊数相交面积；

[0022] 步骤六：生成威胁度的方法为：

[0023] 将融合后结果使用公式:S＝D(l)×0.1+D(ml)×0.3+D(m)×0.5+D(mh)×0.7+D(h)×0.9，得到目标威胁度；所有目标按照步骤一至步骤五生成威胁度矩阵；

[0024] 步骤七：生成命中率矩阵的方法为：

[0025] 计算当前位置导弹和火炮对目标的命中率，

[0026] 导弹命中率公式为：

[0027]

[0028] 火炮命中率公式为：

[0029]

[0030] 其中，p为航路捷径；O为命中率；

[0031] 步骤八：生成火力单元矩阵的方法为：

[0032] 火力单元矩阵表示当前火力单元是否可以工作，当该火力单元不能发射，即视为不可工作；

[0033]

[0034] 火力单元矩阵需根据实际情况生成；

[0035] 步骤九：生成打击效能矩阵的方法为：

[0036] 将步骤六的威胁度矩阵、步骤七的命中率矩阵和步骤八所产生的火力单元矩阵对应相乘得到打击效能矩阵，方法公式为

[0037] E＝S×O×D；E代表打击效能；

[0038] 步骤十：根据威胁性进行火力分配的方法为：

[0039] 对威胁性最大的目标优先进行火力分配，并分配打击效能最大的火力单元，然后将该火力单元和该目标对应火力单元置0，选取下一个威胁性最大的目标按上述方法进行火力分配，当空中目标多于火力单元时，将火力单元全部分配后，对打击效能矩阵进行更新，再进行火力分配，依次迭代。

[0040] 本发明的有益效果在于本发明采用D数和梯形模糊数相结合的方法，具有计算简单的特点；本发明用梯形模糊数对战场信息建模，解决了模糊信息的表示；本发明提出的特征权重加权融合方法，取得了较好的融合效果；本发明提出的火力效能矩阵，为火力分配提供了依据；本发明提出的威胁度优先的火力分配原则，对低空突袭作战模式效果良好。

附图说明

[0041] 图1本发明实现的总流程图。

[0042] 图2是样本数据图。

[0043] 图3是目标1各威胁因素威胁值。

[0044] 图4梯形模糊数隶属度图像

[0045] 图5威胁程度的概率分布

[0046] 图6威胁因素权重

[0047] 图7加权平均融合后的BPA。

[0048] 图8命中率矩阵

[0049] 图9火力效能矩阵

[0050] 图10火力分配矩阵

具体实施方式

[0051] 下面结合附图和实例对本发明进一步说明，此处给出火力分配的实际案例。共给出8个空中目标以及每个目标有5个威胁因素的数据，以目标F1为例，说明所提的威胁度计算具体步骤，以8个目标为例，设置2个导弹火力单元，2个火炮火力单元，说明所提出火力分配方法的实施步骤。

[0052] 步骤一：对空中目标威胁因素评估的方法为：

[0053] 输入8个所探测空中目标，记为F1,F2,...,F8，5种威胁因素的Dij，其中i＝1,2,…,8，j＝1,2,…,5，对于每一个目标，计算各个因素威胁值，所述威胁因素为速度v、航路捷径p、高度h、干扰r和目标类型t，所述威胁值用威胁评估函数计算，分别为
y(r)＝r和y(t)＝t；

[0054] 我们在本例子中选取8个空中目标，数据如图2，以F1为例子计算各个因素威胁值结果如图3。

[0055] 步骤二：建立威胁度评估因子的威胁程度模糊数的方法为：

[0056] 输入空中目标各因素威胁值，根据梯形模糊数建立模型生成威胁程度的基本概率分配函数m(Fjk)，所述威胁程度k分为低l、中低ml、中m、中高mh和高h五个等级，所述模糊数以梯形模糊数表示，对应隶属度函数分别为(0,0,0.1,0.2)、(0.1,0.2,0.3,0.4)、(0.3,0.4,0.6,0.7)、(0.6,0.7,0.8,0.9)和(0.8,0.9,1,1)；所述梯形模糊数是给定论域U上的一个模糊集，是指对任何x∈U，都有一个数μ(x)∈[0,1]与之对应，μ(x)称为x对U的隶属度，μ称为x的隶属函数，

[0057] 隶属度函数图像如图4，F1各个因素对于威胁程度隶属度如图5；

[0058] 步骤三：生成威胁因素权重W1,W2,…W5，权重生成方法为：

[0059] 输入基本概率分配函数，对于同一目标不同威胁因素的分配函数，使用邓熵计算威胁因素权重，所述生成权重公式为 对权重归一化，所述归一化方法为：

[0060] 将图5数据带入上述公式，归一化得到各威胁因素权重如图6；

[0061] 步骤四：将生成的基本概率分配函数用特征权重W1,W2,…W5加权平均得到D(k)方法为：

[0062] 对于一个空中目标，对应威胁因素权重与基本概率分配函数m(Fjk)相乘，加权所述方法公式为

[0063] 根据图5和图6数，使图5每行数据与图6对于威胁因素权重相乘，然后每列相加，即得到加权平均后的数值，称为关于威胁程度等级的概率分布值，如图7；

[0064] 步骤五：使用D数融合方法为：

[0065] 将加权平均得到的D(k)与自身使用D数融合j-1次得到D(A),D数融合规则为[0066]

[0067] 其中

[0068]

[0069] D1(B)和D2(C)是两组待组合的BPA，A、B、C指的是威胁等级，AreaB指的是威胁等级B的梯形模糊数面积，AreaB∩C指威胁等级B和威胁等级C的梯形模糊数相交面积；

[0070] 以图7数据为例，说明融合方法；

[0071] 本方法采用对同一数据融合四次，根据图4，Areamh∩h为mh与h隶属度函数相交部分的面积0.025，Areah为h隶属度函数面积0.15。融合四次可得D(mh)＝0.49，D(h)＝0.51；

[0072] 步骤六：生成威胁度的方法为：

[0073] 将融合后结果使用公式:S＝D(l)×0.1+D(ml)×0.3+D(m)×0.5+D(mh)×0.7+D(h)×0.9，F1的威胁度为：S＝0.49×0.7+0.51×0.9＝0.802；

[0074] 8个空中目标的威胁度计算方法与F1相同，数值如8；

[0075] 步骤七：生成命中率矩阵的方法为：计算当前位置导弹和火炮对目标的命中率；

[0076] 导弹命中率公式为：

[0077]

[0078] 火炮命中率公式为：

[0079] 其中，p为航路捷径；O为命中率

[0080] 将8个目标的航路捷径带入公式，得到4个火力单元对目标的命中率矩阵，如图9；

[0081] 步骤八：生成火力单元矩阵的方法为：

[0082] 火力单元矩阵表示当前火力单元是否可以工作，当该火力单元不能发射，即视为不可工作；

[0083]

[0084] 本例子假设火力单元均可工作；

[0085] 步骤九：生成打击效能矩阵的方法为：

[0086] 将步骤六的威胁度矩阵、步骤七的命中率矩阵和步骤八所产生的火力单元矩阵对应相乘得到打击效能矩阵，方法公式为：

[0087] E＝S×O×D；E代表打击效能；

[0088] 步骤十：根据威胁性进行火力分配的方法为：

[0089] 对威胁性最大的目标优先进行火力分配，并分配打击效能最大的火力单元，然后将该火力单元和该目标对应火力单元置0，选取下一个威胁性最大的目标按上述方法进行火力分配。当空中目标多于火力单元时，将火力单元全部分配后，对打击效能矩阵进行更新，再进行火力分配，依次迭代；

[0090] 火力分配方案如图9，前两行为第一次火力分配，后两列为第二次火力分配。