CN114444203A - 基于战场态势变权的空中集群威胁评估方法 - Google Patents

基于战场态势变权的空中集群威胁评估方法 Download PDF

Info

Publication number
CN114444203A
CN114444203A CN202210084034.8A CN202210084034A CN114444203A CN 114444203 A CN114444203 A CN 114444203A CN 202210084034 A CN202210084034 A CN 202210084034A CN 114444203 A CN114444203 A CN 114444203A
Authority
CN
China
Prior art keywords
cluster
threat
index
weight
calculating
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
Application number
CN202210084034.8A
Other languages
English (en)
Other versions
CN114444203B (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.)
Air Force Engineering University of PLA
Original Assignee
Air Force Engineering University of 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 Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN202210084034.8A priority Critical patent/CN114444203B/zh
Publication of CN114444203A publication Critical patent/CN114444203A/zh
Application granted granted Critical
Publication of CN114444203B publication Critical patent/CN114444203B/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/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

本发明提出一种基于战场态势变权的空中集群威胁评估方法,包括下列步骤:建立空中集群威胁评估指标体系;构建评估指标的量化模型;基于层次分析法确定评估指标的主观权重;基于熵权法确定客观权重;基于博弈组合确定常权权重;确定指标变权权重;确定威胁隶属度矩阵的正负理想点;计算感知效用值;确定感知效用矩阵的正负理想解;计算与正负理想解的关联系数;计算与正负理想解的加权灰色关联度;计算与正负理想解的加权欧式距离;基于综合贴近度的待评估空中集群威胁评估。本发明能够解决传统赋权方法指标权重反映战场态势特点和变化不足的问题,以及传统多属性决策方法没有考虑决策者心理因素以及灰色关联分析和TOPSIS法评估结果不一致的问题。

Description

基于战场态势变权的空中集群威胁评估方法
技术领域
本发明涉及防空作战指挥控制技术领域,具体涉及一种基于战场态势变权的空中集群威胁评估方法。
背景技术
空中集群威胁评估是防空作战指挥控制领域的核心问题。由于空中目标集群的组成复杂多变,空袭目标特征属性众多,对空中集群进行威胁评估需综合目标的特征属性、集群的组成样式和战场态势的动态变化,其实质是一类不确定动态多属性群决策问题。
在联合空中作战中,为提高指挥协同的效率,提升空中作战效能,集群作战已经成为现代空中作战的主要样式。威胁评估作为作战过程的重要环节,是后续兵力部署和指挥决策的重要依据。但目前对于空中目标威胁评估的研究大都停留在战术层次,即仅考虑目标数量较少的情况下,一般采用多属性决策、直觉模糊集、神经网络、支持向量机和贝叶斯网络等方法,选取目标类型、航路捷径、飞行速度和飞行高度等目标属性,对单个来袭目标进行威胁评估和排序。而在联合空中作战中,以多种不同类型的作战飞机组成的空中集群已经成为基本作战单元。因此,从战役层次对空中集群的威胁度进行合理有效的评估,对于联合空中作战的指挥决策具有重要意义。但目前对于空中集群的威胁评估研究较少,方法有限。
发明内容
针对以上不足,本发明提出一种基于战场态势变权的空中集群威胁评估方法,具体包括下列步骤:
Step 1建立空中集群威胁评估指标体系
为充分反映空中集群的威胁程度,选取火力打击能力、指控信息能力、战场机动能力、生存防护能力、综合保障能力作为空中集群威胁评估的指标,建立空中集群威胁评估指标体系;
Step 2构建评估指标的量化模型
根据联合空中作战的原则和空中集群的构成特点,将空中集群按照单机、编组和集群三个层次,构建从单机到集群的层次聚合威胁指标量化模型,进而得到空中集群的威胁属性值;
Step 2.1计算单机威胁度
构建单架飞机的威胁度模型为:
fj=-ln(1-p(j)) (1)
式中:fj表示单架飞机在第j个评估指标下的威胁度,p(j)表示单架飞机取得第j个评估指标下作战效果的概率,j为正整数;
Step 2.2计算编组威胁度
根据作战实际,将空中集群中的同类飞机作为一个编组进行处理,编组威胁度通过单机威胁指数线性相加的聚合方式得到,计算方法为:
Ibj=Nb·fj (2)
式中:Ibj为第b个编组在第j个指标下的威胁度,Nb为集群中第b个编组内飞机的数量,b为正整数;
Step 2.3计算集群威胁度
通过幂指数模型将不同编组的威胁度进行聚合,空中集群在第j个评估指标威胁度的计算方法为:
Figure BSA0000264684700000021
式中:K为调整系数,Ej为空中集群在第j个指标的威胁度,B为集群中编组数量,B为正整数;αbj为幂指数,反映第b个编组对于第j个指标的重要性;
Step 2.4引入协同系数修正空中集群威胁度
引入协同系数对空中集群威胁度进行修正,协同系数计算方法为:
Figure BSA0000264684700000031
式中:ebc为第b个编组中第c项协同内容的协同度值,c为正整数,C为总的协同内容项数,C为正整数,c=1,2,…,C;αbc和βbc分别表示ebc的最大值和最小值;λb为第b个编组在空中集群的权重;
通过协同系数修正得到最终的空中集群在第j个评估指标威胁量化值为:
Zj=ρ·Ej (5)
Step 3基于层次分析法确定评估指标的主观权重
Step 3.1构建判断矩阵
假设有n个威胁评估指标,采用Saaty九标度对威胁评估指标采用两两比较的方式构建判断矩阵;
Figure BSA0000264684700000032
其中,aij为第i个评估指标相对于第j个评估指标的重要程度,i,j=1,2,…,n;
Step 3.2对判断矩阵进行一致性检验
一致性比例的计算方法为:
Figure BSA0000264684700000033
式中:λmax为判断矩阵的最大特征值,n为指标数量,RI为平均随机一致性指标;
若一致性比例CR小于0.1,则认为判断矩阵的一致性能够接受,继续下一步;否则按照Step3.1重新构造判断矩阵,直到一致性比例满足要求为止;
Step 3.3计算指标主观权重
根据矩阵论的基础知识,求出判断矩阵的最大特征值及其对应的特征向量,对特征向量进行归一化处理,得到评估指标的主观权重向量w=(w1,w2,…,wn),wj(j=1,2…,n)为第j个评估指标的权重;
Step 4基于熵权法确定客观权重
Step 4.1构建威胁属性值矩阵
假设待评估空中集群有m个,评估指标共有n个,则威胁属性值矩阵为:
Figure BSA0000264684700000041
其中,zkj为第k个待评估空中集群在第j个评估指标下的威胁属性值,k=1,2,…,m,j=1,2,…,n;
Step 4.2计算评估指标的信息熵
第j个评估指标的信息熵为:
Figure BSA0000264684700000042
Step 4.3计算评估指标的熵权
第j个评估指标的熵权为:
Figure BSA0000264684700000043
得到评估指标的客观权重向量v=(v1,v2,…vn);
Step 5基于博弈组合确定常权权重
Step 5.1将主观权重向量w和客观权重向量v线性组合得到的指标常权权重向量为:
Figure BSA0000264684700000051
式中:λ1,λ2为第一和第二线性组合系数;
Step 5.2根据博弈论思想,建立η与w和v离差和最小的目标函数为:
Figure BSA0000264684700000052
Step 5.3根据微分原理,使上述目标函数最小需要满足的一阶导数条件为:
Figure BSA0000264684700000053
标准化处理得:
Figure BSA0000264684700000054
式中,
Figure BSA0000264684700000055
为第一和第二标准化线性组合系数;
得到最终指标常权权重向量为:
Figure BSA0000264684700000056
Step 6确定指标变权权重
Step 6.1计算状态变权参数
第i个待评估空中集群在第j个评估指标的状态变权参数计算方法为:
Figure BSA0000264684700000057
式中:K1为激励幅度系数;K2为惩罚幅度系数;zij为第i个待评估空中集群在第j个评估指标下的威胁属性值;s1,s2,s3分别为激励变权指标集合、常权指标集合和惩罚变权指标集合,s1,s2,s3根据战场态势进行确定;
Step 6.2计算指标变权值
第i个待评估空中集群在第j个评估指标变权值为:
Figure BSA0000264684700000061
式中:о为哈达玛(Hadamard)乘积;
Step 7确定威胁隶属度矩阵的正负理想点
从威胁属性值矩阵Z中取
Figure BSA0000264684700000062
为正理想点,取
Figure BSA0000264684700000063
为负理想点,其中:
Figure BSA0000264684700000064
Figure BSA0000264684700000065
Step 8计算感知效用值
第i个待评估空中集群在第j个评估指标的感知效用表达式为:
Figure BSA0000264684700000066
式中:
Figure BSA0000264684700000067
为后悔值,
Figure BSA0000264684700000068
为欣喜值,R(·)为后悔欣喜函数,为单调递增的凹函数,且满足R′(·)>0,R″(·)<0和R(0)=0;
以正理想点作为参考,计算待评估空中集群的后悔值为:
Figure BSA0000264684700000069
以负理想点作为参考,计算待评估空中集群的欣喜值为:
Figure BSA00002646847000000610
式中:δ(δ>0)为后悔规避系数,δ越大则决策者的后悔规避系数越大;
根据后悔和欣喜值得到感知效用矩阵为:
Figure BSA00002646847000000611
Step 9确定感知效用矩阵的正负理想解
从感知效用矩阵U中取
Figure BSA00002646847000000612
为正理想解,取
Figure BSA00002646847000000613
作为负理想解,其中:
Figure BSA0000264684700000071
Figure BSA0000264684700000072
Step 10计算与正负理想解的关联系数
计算第i个待评估空中集群第j个评估指标与正、负理想解的关联系数
Figure BSA0000264684700000073
Figure BSA0000264684700000074
Figure BSA0000264684700000075
Figure BSA0000264684700000076
式中:ρ为分辨系数,通常取0.5;
Step11计算与正负理想解的加权灰色关联度
计算第i个待评估评估集群与正、负理想解的加权灰色关联度
Figure BSA0000264684700000077
Figure BSA0000264684700000078
Figure BSA0000264684700000079
Figure BSA00002646847000000710
Step12计算与正负理想解的加权欧式距离
计算第i个待评估评估集群与正负理想解的加权欧式距离
Figure BSA00002646847000000711
Figure BSA00002646847000000712
Figure BSA00002646847000000713
Figure BSA00002646847000000714
Tep13基于综合贴近度的待评估空中集群威胁评估
第i个待评估空中集群与正、负理想解的贴近度分别为:
Figure BSA00002646847000000715
Figure BSA00002646847000000716
其中,ω为权重偏好因子,ω∈[0,1];
得到第i个待评估空中集群的综合贴近度为:
Figure BSA0000264684700000081
根据综合贴近度对空中集群进行排序,综合贴近度越大,说明集群的威胁程度越大,反之说明威胁程度越小。
本发明基于战场态势变权的空中集群威胁评估方法。采用层次分析法和熵权法得到主客观权重,通过博弈组合得到常权权重,引入变权理论确定变权权重,解决了传统赋权方法指标权重反映战场态势特点和变化不足的问题;构建基于后悔理论的灰色TOPSIS威胁评估模型,通过计算综合贴近度得到最终威胁评估结果,解决了传统多属性决策方法没有考虑决策者心理因素以及灰色关联分析和TOPSIS法评估结果不一致的问题。
附图说明
图1示出本发明一种基于战场态势变权的空中集群威胁评估方法处理流程;
图2示出空中集群威胁评估指标体系。
具体实施方式
为使本发明的目的、技术路线及优点更加清楚,下面结合附图及实施例,对本发明进行进一步详细说明。
如图1所示,一种基于战场态势变权的空中集群威胁评估方法包含以下步骤:
Step 1建立空中集群威胁评估指标体系
根据联合空中作战中空中集群的作战特点,对于空中集群的威胁评估主要针对其作战效能,即空中集群的作战能够对我方产生多大的威胁。因此,为了充分反映空中集群的威胁程度,选取火力打击能力、指控信息能力、战场机动能力、生存防护能力、综合保障能力作为空中集群威胁评估的指标,建立空中集群威胁评估指标体系如图2所示。
Step 2构建评估指标的量化模型
根据联合空中作战的原则和空中集群的构成特点,将空中集群按照单机、编组和集群三个层次,构建从单机到集群的层次聚合威胁指标量化模型,进而得到空中集群的威胁属性值。
Step 2.1计算单机威胁度
单机的威胁度可通过作战效能进行反映,使用作战概率度量武器装备的作战效能,武器装备达成某一作战效果的概率越大,则该作战效果对应的作战效能就越大,对应的威胁程度也越大。构建单架飞机的威胁度模型为:
fj=-ln(1-p(j)) (1)
式中:fj表示单架飞机在第j个评估指标下的威胁度,p(j)表示单架飞机取得第j个评估指标下作战效果的概率,j为正整数。
Step 2.2计算编组威胁度
根据作战实际,将空中集群中的同类飞机作为一个编组进行处理,编组威胁度通过单机威胁指数线性相加的聚合方式得到,计算方法为:
Ibj=Nb·fj (2)
式中:Ibj为第b个编组在第j个指标下的威胁度,Nb为集群中第b个编组内飞机的数量,b为正整数。
Step 2.3计算集群威胁度
一个空中集群往往由多个编组所组成,不同编组包含的飞机类型不同,因此对指标威胁值的重要性不同,通过幂指数模型将不同编组的威胁度进行聚合,空中集群在第j个评估指标威胁度的计算方法为:
Figure BSA0000264684700000101
式中:K为调整系数,Ej为空中集群在第j个指标的威胁度,B为集群中编组数量,B为正整数。αbj为幂指数,反映第b个编组对于第j个指标的重要性。
Step 2.4引入协同系数修正空中集群威胁度
在联合空中作战中,空中集群往往是由多军兵种的不同类型飞机所组成,通过协同遂行各种作战任务。协同系数作为联合作战中反映各编组之间协同能力的数量表示,对空中集群的威胁程度有着直接的影响,因此引入协同系数对空中集群威胁度进行修正,协同系数计算方法为:
Figure BSA0000264684700000102
式中:ebc为第b个编组中第c项协同内容的协同度值,c为正整数,C为总的协同内容项数,C为正整数,c=1,2,…,C,一般协同内容包括火力协同、时间协同、空间协同和信息协同等;αbc和βbc分别表示ebc的最大值和最小值;λb为第b个编组在空中集群的权重。
通过协同系数修正得到最终的空中集群在第j个评估指标威胁量化值为:
Zj=ρ·Ej (5)
Step 3基于层次分析法确定评估指标的主观权重
Step 3.1构建判断矩阵
假设有n个威胁评估指标,采用Saaty九标度(Saaty九标度为本领域技术人员熟知,不再赘述)对威胁评估指标采用两两比较的方式构建判断矩阵;
Figure BSA0000264684700000103
其中,aij为第i个评估指标相对于第j个评估指标的重要程度,i,j=1,2,…,n。
Step 3.2对判断矩阵进行一致性检验
一致性比例的计算方法为:
Figure BSA0000264684700000111
式中:λmax为判断矩阵的最大特征值,n为指标数量,RI为平均随机一致性指标,可以通过查表获得(查表获得RI值的方法为本领域技术人员熟知,不再赘述)。
若一致性比例CR小于0.1,则认为判断矩阵的一致性可以接受,继续下一步;否则按照Step3.1重新构造判断矩阵,直到一致性比例满足要求为止。
Step 3.3计算指标主观权重
根据矩阵论的基础知识,求出判断矩阵的最大特征值及其对应的特征向量,对特征向量进行归一化处理,得到评估指标的主观权重向量w=(w1,w2,…,wn),wj(j=1,2…,n)为第j个评估指标的权重。该步骤的实施办法,为本领域技术人员熟知,不再累述。
Step 4基于熵权法确定客观权重
Step 4.1构建威胁属性值矩阵
假设待评估空中集群有m个,评估指标共有n个,则威胁属性值矩阵为:
Figure BSA0000264684700000112
其中,zkj为第k个待评估空中集群在第j个评估指标下的威胁属性值,k=1,2,…,m,j=1,2,…,n。
Step 4.2计算评估指标的信息熵
第j个评估指标的信息熵为:
Figure BSA0000264684700000113
Step 4.3计算评估指标的熵权
第j个评估指标的熵权为:
Figure BSA0000264684700000121
可得评估指标的客观权重向量v=(v1,v2,…vn)。
Step 5基于博弈组合确定常权权重
通过博弈论的思想进行组合赋权,将主观权重作为双方博弈的一方,将客观权重作为博弈的另一方,当博弈双方达到纳什均衡状态时得到的指标权重最合理,且主观权重和客观权重的离差之和最小。
Step 5.1将主观权重向量w和客观权重向量v线性组合得到的指标常权权重向量为:
Figure BSA0000264684700000122
式中:λ1,λ2为第一和第二线性组合系数。
Step 5.2根据博弈论思想,建立η与w和v离差和最小的目标函数为:
Figure BSA0000264684700000123
Step 5.3根据微分原理,使上述目标函数最小需要满足的一阶导数条件为:
Figure BSA0000264684700000124
标准化处理得:
Figure BSA0000264684700000125
式中,
Figure BSA0000264684700000126
为第一和第二标准化线性组合系数。
得到最终指标常权权重向量为:
Figure BSA0000264684700000127
Step 6确定指标变权权重
博弈组合权重虽然能够兼顾专家主观经验和指标客观数据,得到较为合理的权重,但当空中战场态势发生变化时,常权评估无法反映战场态势的特点和动态变化,容易造成评估结果的不合理。通过变权理论构造状态变权向量,计算相应战场态势下不同待评估空中集群的指标变权权重向量,能够得到更加科学的威胁评估结果。
Step 6.1计算状态变权参数
不同的作战态势会直接影响不同评估指标对于威胁评估的重要性程度,比如在敌方进行战略空袭的作战态势下,火力打击能力和战场机动能力的作用更加突出,对我方威胁的影响也越大,应做激励性变权处理;而在空袭态势下指控信息能力和生存防护能力对威胁程度的影响相比就会下降,应做惩罚变权处理;综合保障能力的作用在空袭态势下变化不大,可做常权处理处理,但由于需要满足权重归一化条件,因此也会相应变化。另外,不同待评估空中集群由于在同一评估指标下的威胁值不同,惩罚和激励的幅度也需要与威胁值的大小相适应。
第i个待评估空中集群在第j个评估指标的状态变权参数计算方法为:
Figure BSA0000264684700000131
式中:K1为激励幅度系数;K2为惩罚幅度系数;zij为第i个待评估空中集群在第j个评估指标下的威胁属性值;s1,s2,s3分别为激励变权指标集合、常权指标集合和惩罚变权指标集合,s1,s2,s3根据战场态势进行确定。
Step 6.2计算指标变权值
第i个待评估空中集群在第j个评估指标变权值为:
Figure BSA0000264684700000132
式中:о为哈达玛(Hadamard)乘积。
Step 7确定威胁隶属度矩阵的正负理想点
从威胁属性值矩阵Z中取
Figure BSA0000264684700000141
为正理想点,取
Figure BSA0000264684700000142
为负理想点,其中:
Figure BSA0000264684700000143
Figure BSA0000264684700000144
Step 8计算感知效用值
根据期望效用理论,决策者的感知效用值随着欣喜值和后悔值的变化产生波动,由本身的效用值、后悔值和欣喜值三部分组成,第i个待评估空中集群在第j个评估指标的感知效用表达式为:
Figure BSA0000264684700000145
式中:
Figure BSA0000264684700000146
为后悔值,
Figure BSA0000264684700000147
为欣喜值,R(·)为后悔欣喜函数,为单调递增的凹函数,且满足R′(·)>0,R″(·)<0和R(0)=0。
以正理想点作为参考,计算待评估空中集群的后悔值为:
Figure BSA0000264684700000148
以负理想点作为参考,计算待评估空中集群的欣喜值为:
Figure BSA0000264684700000149
式中:δ(δ>0)为后悔规避系数,δ越大则决策者的后悔规避系数越大。
根据后悔和欣喜值得到感知效用矩阵为:
Figure BSA00002646847000001410
Step 9确定感知效用矩阵的正负理想解
从感知效用矩阵U中取
Figure BSA00002646847000001411
为正理想解,取
Figure BSA00002646847000001412
作为负理想解,其中:
Figure BSA00002646847000001413
Figure BSA00002646847000001414
Step 10计算与正负理想解的关联系数
计算第i个待评估空中集群第j个评估指标与正、负理想解的关联系数
Figure BSA00002646847000001415
Figure BSA00002646847000001416
Figure BSA0000264684700000151
Figure BSA0000264684700000152
式中:ρ为分辨系数,通常取0.5。
Step11计算与正负理想解的加权灰色关联度
计算第i个待评估评估集群与正、负理想解的加权灰色关联度
Figure BSA0000264684700000153
Figure BSA0000264684700000154
Figure BSA0000264684700000155
Figure BSA0000264684700000156
Step12计算与正负理想解的加权欧式距离
计算第i个待评估评估集群与正负理想解的加权欧式距离
Figure BSA0000264684700000157
Figure BSA0000264684700000158
Figure BSA0000264684700000159
Figure BSA00002646847000001510
Tep13基于综合贴近度的待评估空中集群威胁评估
第i个待评估空中集群与正、负理想解的贴近度分别为:
Figure BSA00002646847000001511
Figure BSA00002646847000001512
其中,ω为权重偏好因子,ω∈[0,1]。
得到第i个待评估空中集群的综合贴近度为:
Figure BSA00002646847000001513
根据综合贴近度对空中集群进行排序,综合贴近度越大,说明集群的威胁程度越大,反之说明威胁程度越小。
具体实施例
为验证本发明的可行性和有效性,假设在联合空中作战中发现了大规模的来袭目标,经过处理得到6个空中集群及其具体信息;受篇幅所限,仅显示第一个空中集群的兵力组成及作战效果概率情况如表1所示
表1空中集群信息
Figure BSA0000264684700000161
由式(1-2)得到编组威胁度矩阵为:
Figure BSA0000264684700000162
根据特征向量法得到幂指数矩阵为:
Figure BSA0000264684700000163
由式(3-5)得到该集群的威胁属性值向量为:
(5.5704,4.2102,5.0522,2.3917,3.2779)
同理可计算其他5个空中集群的威胁属性值向量,得到威胁属性值矩阵为:
Figure BSA0000264684700000171
根据层次分析法得到主观权重为:
(0.398,0.263,0.136,0.125,0.078)
由式(7-8)得到客观权重为:
(0.378,0.096,0.212,0.183,0.131)
由式(9-12)得到线性组合系数为:
Figure BSA0000264684700000172
由式(13)得到指标常权权重为:
(0.393,0.219,0.156,0.140,0.092)
由式(14-15)得到各空中集群的指标变权权重如表2所示
表2指标变权权重
Figure BSA0000264684700000173
取后悔规避系数为0.5,由式(16-18)得到感知效用矩阵为
Figure BSA0000264684700000174
从感知效用矩阵U中取出正、负理想解分别为
Figure BSA0000264684700000175
Figure BSA0000264684700000176
得到与正、负理想解的关联系数矩阵分别为
Figure BSA0000264684700000181
Figure BSA0000264684700000182
由式(19-24)得到加权灰色关联度和加权欧式距离如表3所示
表3加权灰色关联度和加权欧式距离
Figure BSA0000264684700000183
取偏好系数为0.5,表示欧式距离与曲线关联度同样重要,由式(25-27)得到综合贴近度为:
(0.6855,0.3881,0.5670,0.5212,0.4817,0.3695)
根据综合贴近度得到威胁评估结果,6个空中集群的威胁度排序为:
集群1>集群3>集群4>集群5>集群2>集群6。
一种基于战场态势变权的空中集群威胁评估方法。首先,根据联合空中作战实际,建立了空中集群的威胁评估指标体系,根据离散信源与武器装备的相似原理构造了基于层次聚合的威胁指标量化模型,构建了一种从单机到集群的层次聚合模型对空中集群威胁度进行量化;接下来,引入变权理论根据战场态势确定指标变权权重,采用层次分析法和熵权法得到主客观权重,通过博弈组合得到常权权重,引入变权理论确定变权权重,解决了传统赋权方法指标权重反映战场态势特点和变化不足的问题;最后,构建基于后悔理论的灰色TOPSIS威胁评估模型,通过计算综合贴近度得到最终威胁评估结果,解决了传统多属性决策方法没有考虑决策者心理因素以及灰色关联分析和TOPSIS法评估结果不一致的问题;本发明解决了现有威胁评估方法忽略空中集群威胁和战场态势变化而造成评估精度不高的问题,为新一代防空智能指挥控制系统的研制提供方法支持。
应当理解,上述具体实施例仅仅用于解释本发明,并不用于限定本发明。

Claims (1)

1.一种基于战场态势变权的空中集群威胁评估方法,其特征在于,具体包括下列步骤:
Step 1建立空中集群威胁评估指标体系
为充分反映空中集群的威胁程度,选取火力打击能力、指控信息能力、战场机动能力、生存防护能力、综合保障能力作为空中集群威胁评估的指标,建立空中集群威胁评估指标体系;
Step 2构建评估指标的量化模型
根据联合空中作战的原则和空中集群的构成特点,将空中集群按照单机、编组和集群三个层次,构建从单机到集群的层次聚合威胁指标量化模型,进而得到空中集群的威胁属性值;
Step2.1计算单机威胁度
构建单架飞机的威胁度模型为:
fj=-ln(1-p(j)) (1)
式中:fj表示单架飞机在第j个评估指标下的威胁度,p(j)表示单架飞机取得第j个评估指标下作战效果的概率,j为正整数;
Step 2.2计算编组威胁度
根据作战实际,将空中集群中的同类飞机作为一个编组进行处理,编组威胁度通过单机威胁指数线性相加的聚合方式得到,计算方法为:
Ibj=Nb·fj (2)
式中:Ibj为第b个编组在第j个指标下的威胁度,Nb为集群中第b个编组内飞机的数量,b为正整数;
Step 2.3计算集群威胁度
通过幂指数模型将不同编组的威胁度进行聚合,空中集群在第j个评估指标威胁度的计算方法为:
Figure FSA0000264684690000021
式中:K为调整系数,Ej为空中集群在第j个指标的威胁度,B为集群中编组数量,B为正整数;αbj为幂指数,反映第b个编组对于第j个指标的重要性;
Step 2.4引入协同系数修正空中集群威胁度
引入协同系数对空中集群威胁度进行修正,协同系数计算方法为:
Figure FSA0000264684690000022
式中:ebc为第b个编组中第c项协同内容的协同度值,c为正整数,C为总的协同内容项数,C为正整数,c=1,2,…,C;αbc和βbc分别表示ebc的最大值和最小值;λb为第b个编组在空中集群的权重;
通过协同系数修正得到最终的空中集群在第j个评估指标威胁量化值为:
Zj=ρ·Ej (5)
Step 3基于层次分析法确定评估指标的主观权重
Step 3.1构建判断矩阵
假设有n个威胁评估指标,采用Saaty九标度对威胁评估指标采用两两比较的方式构建判断矩阵;
Figure FSA0000264684690000023
其中,aij为第i个评估指标相对于第j个评估指标的重要程度,i,j=1,2,…,n;
Step 3.2对判断矩阵进行一致性检验
一致性比例的计算方法为:
Figure FSA0000264684690000024
式中:λmax为判断矩阵的最大特征值,n为指标数量,RI为平均随机一致性指标;
若一致性比例CR小于0.1,则认为判断矩阵的一致性能够接受,继续下一步;否则按照Step3.1重新构造判断矩阵,直到一致性比例满足要求为止;
Step 3.3计算指标主观权重
根据矩阵论的基础知识,求出判断矩阵的最大特征值及其对应的特征向量,对特征向量进行归一化处理,得到评估指标的主观权重向量w=(w1,w2,…,wn),wj(j=1,2…,n)为第j个评估指标的权重;
Step 4基于熵权法确定客观权重
Step 4.1构建威胁属性值矩阵
假设待评估空中集群有m个,评估指标共有n个,则威胁属性值矩阵为:
Figure FSA0000264684690000031
其中,zkj为第k个待评估空中集群在第j个评估指标下的威胁属性值,k=1,2,…,m,j=1,2,…,n;
Step 4.2计算评估指标的信息熵
第j个评估指标的信息熵为:
Figure FSA0000264684690000032
Step 4.3计算评估指标的熵权
第j个评估指标的熵权为:
Figure FSA0000264684690000033
得到评估指标的客观权重向量v=(v1,v2,…vn);
Step 5基于博弈组合确定常权权重
Step 5.1将主观权重向量w和客观权重向量v线性组合得到的指标常权权重向量为:
Figure FSA0000264684690000041
式中:λ1,λ2为第一和第二线性组合系数;
Step 5.2根据博弈论思想,建立η与w和v离差和最小的目标函数为:
Figure FSA0000264684690000042
Step 5.3根据微分原理,使上述目标函数最小需要满足的一阶导数条件为:
Figure FSA0000264684690000043
标准化处理得:
Figure FSA0000264684690000044
式中,
Figure FSA0000264684690000045
为第一和第二标准化线性组合系数;
得到最终指标常权权重向量为:
Figure FSA0000264684690000046
Step 6确定指标变权权重
Step 6.1计算状态变权参数
第i个待评估空中集群在第j个评估指标的状态变权参数计算方法为:
Figure FSA0000264684690000047
式中:K1为激励幅度系数;K2为惩罚幅度系数;zij为第i个待评估空中集群在第j个评估指标下的威胁属性值;s1,s2,s3分别为激励变权指标集合、常权指标集合和惩罚变权指标集合,s1,s2,s3根据战场态势进行确定;
Step 6.2计算指标变权值
第i个待评估空中集群在第j个评估指标变权值为:
Figure FSA0000264684690000051
式中:
Figure FSA0000264684690000052
为哈达玛(Hadamard)乘积;
Step 7确定威胁隶属度矩阵的正负理想点
从威胁属性值矩阵Z中取
Figure FSA0000264684690000053
为正理想点,取
Figure FSA0000264684690000054
为负理想点,其中:
Figure FSA0000264684690000055
Figure FSA0000264684690000056
Step 8计算感知效用值
第i个待评估空中集群在第j个评估指标的感知效用表达式为:
Figure FSA0000264684690000057
式中:
Figure FSA0000264684690000058
为后悔值,
Figure FSA0000264684690000059
为欣喜值,R(·)为后悔欣喜函数,为单调递增的凹函数,且满足R′(·)>0,R″(·)<0和R(0)=0;
以正理想点作为参考,计算待评估空中集群的后悔值为:
Figure FSA00002646846900000510
以负理想点作为参考,计算待评估空中集群的欣喜值为:
Figure FSA00002646846900000511
式中:δ(δ>0)为后悔规避系数,δ越大则决策者的后悔规避系数越大;
根据后悔和欣喜值得到感知效用矩阵为:
Figure FSA00002646846900000512
Step 9确定感知效用矩阵的正负理想解
从感知效用矩阵U中取
Figure FSA00002646846900000513
为正理想解,取
Figure FSA00002646846900000514
作为负理想解,其中:
Figure FSA0000264684690000061
Figure FSA0000264684690000062
Step 10计算与正负理想解的关联系数
计算第i个待评估空中集群第j个评估指标与正、负理想解的关联系数
Figure FSA0000264684690000063
Figure FSA0000264684690000064
Figure FSA0000264684690000065
Figure FSA0000264684690000066
式中:ρ为分辨系数,通常取0.5;
Step11计算与正负理想解的加权灰色关联度
计算第i个待评估评估集群与正、负理想解的加权灰色关联度
Figure FSA0000264684690000067
Figure FSA0000264684690000068
Figure FSA0000264684690000069
Figure FSA00002646846900000610
Step12计算与正负理想解的加权欧式距离
计算第i个待评估评估集群与正负理想解的加权欧式距离
Figure FSA00002646846900000611
Figure FSA00002646846900000612
Figure FSA00002646846900000613
Figure FSA00002646846900000614
Tep13基于综合贴近度的待评估空中集群威胁评估
第i个待评估空中集群与正、负理想解的贴近度分别为:
Figure FSA00002646846900000615
Figure FSA00002646846900000616
其中,ω为权重偏好因子,ω∈[0,1];
得到第i个待评估空中集群的综合贴近度为:
Figure FSA0000264684690000071
根据综合贴近度对空中集群进行排序,综合贴近度越大,说明集群的威胁程度越大,反之说明威胁程度越小。
CN202210084034.8A 2022-01-25 2022-01-25 基于战场态势变权的空中集群威胁评估方法 Active CN114444203B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210084034.8A CN114444203B (zh) 2022-01-25 2022-01-25 基于战场态势变权的空中集群威胁评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210084034.8A CN114444203B (zh) 2022-01-25 2022-01-25 基于战场态势变权的空中集群威胁评估方法

Publications (2)

Publication Number Publication Date
CN114444203A true CN114444203A (zh) 2022-05-06
CN114444203B CN114444203B (zh) 2024-04-05

Family

ID=81369154

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210084034.8A Active CN114444203B (zh) 2022-01-25 2022-01-25 基于战场态势变权的空中集群威胁评估方法

Country Status (1)

Country Link
CN (1) CN114444203B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545584A (zh) * 2022-12-02 2022-12-30 中国兵器科学研究院 防空威胁控制方法、装置、电子设备及存储介质
CN116739431A (zh) * 2023-06-28 2023-09-12 北京流体动力科学研究中心 基于层次分析法的飞行器实时威胁评估方法
CN117408493A (zh) * 2023-12-08 2024-01-16 中国人民解放军海军航空大学 一种要地一体化防空平台协同方法、系统及介质
CN117436282A (zh) * 2023-12-07 2024-01-23 南京信息工程大学 基于双重变权和topsis-灰关联的架构优选方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005124580A1 (en) * 2004-06-15 2005-12-29 The University Of Melbourne A threat assessment system and process
CN109146334A (zh) * 2018-10-09 2019-01-04 中国人民解放军海军航空大学 一种适用于高动态环境的多目标威胁评估方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005124580A1 (en) * 2004-06-15 2005-12-29 The University Of Melbourne A threat assessment system and process
CN109146334A (zh) * 2018-10-09 2019-01-04 中国人民解放军海军航空大学 一种适用于高动态环境的多目标威胁评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LONGFEI YUE等: "Air Target Threat Assessment Based on Improved Moth Flame Optimization-Gray Neural Network Model", MATH- EMATICALPROBLEMSINENGINEERING, 10 October 2019 (2019-10-10) *
杜鹏枭;卢盈齐;: "基于惩罚与激励变权和TOPSIS的空中目标威胁评估", 军械工程学院学报, no. 03, 15 June 2017 (2017-06-15) *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545584A (zh) * 2022-12-02 2022-12-30 中国兵器科学研究院 防空威胁控制方法、装置、电子设备及存储介质
CN116739431A (zh) * 2023-06-28 2023-09-12 北京流体动力科学研究中心 基于层次分析法的飞行器实时威胁评估方法
CN116739431B (zh) * 2023-06-28 2024-03-26 北京流体动力科学研究中心 基于层次分析法的飞行器实时威胁评估方法
CN117436282A (zh) * 2023-12-07 2024-01-23 南京信息工程大学 基于双重变权和topsis-灰关联的架构优选方法
CN117436282B (zh) * 2023-12-07 2024-04-02 南京信息工程大学 基于双重变权和topsis-灰关联的架构优选方法
CN117408493A (zh) * 2023-12-08 2024-01-16 中国人民解放军海军航空大学 一种要地一体化防空平台协同方法、系统及介质
CN117408493B (zh) * 2023-12-08 2024-03-01 中国人民解放军海军航空大学 一种要地一体化防空平台协同方法、系统及介质

Also Published As

Publication number Publication date
CN114444203B (zh) 2024-04-05

Similar Documents

Publication Publication Date Title
CN114444203A (zh) 基于战场态势变权的空中集群威胁评估方法
CN107844835B (zh) 基于动态权重m-topsis多属性决策的多目标优化改进遗传算法
CN110502031A (zh) 基于任务需求的异构无人机集群协同最优配置方法
CN108564592A (zh) 基于动态多种群集成差分进化算法的图像分割方法
CN108983800B (zh) 一种基于深度学习的飞机姿态控制方法
CN110490422A (zh) 一种基于博弈云模型的目标作战效能态势评估方法
CN106934225A (zh) 一种电力巡检用多旋翼无人机系统效能评估方法
CN112163719B (zh) 基于指派模型的导弹装备系统弹种优化设计方法及装置
CN112541273A (zh) 一种基于改进pcm-幂指数法的非致命武器效能评估方法
CN113177583A (zh) 一种空中目标聚类分群方法
CN116739428A (zh) 一种主客观赋权的变权topsis炮兵火力打击目标价值分析方法
CN116934148A (zh) 一种基于聚类组合赋权的空袭目标威胁评估方法
CN110782062A (zh) 一种用于防空系统的多对多分组拦截目标分配方法及系统
CN115186378A (zh) 一种空战模拟环境中的战术控制距离实时解算方法
CN115661576A (zh) 一种样本不平衡下的飞机群组意图识别方法
Yangyang et al. A new method of multi-target threat assessment for air combat
CN114065386A (zh) 一种快速预计整机mtbf的多元线性回归方程方法
CN113792984A (zh) 基于云模型的防空反导指挥控制模型能力评估方法
CN116739431B (zh) 基于层次分析法的飞行器实时威胁评估方法
CN109866931B (zh) 一种基于自编码器的飞机油门控制方法
CN117436282B (zh) 基于双重变权和topsis-灰关联的架构优选方法
Wu et al. Fuzzy Dynamic Interval AHP Threat Assessment Method for Low Observable Targets in Air and Space
Wu et al. Multi-sensor management based on information increment and covariance of target priority discrimination
CN117590757B (zh) 基于高斯分布海鸥优化算法的多无人机协同任务分配方法
CN115964640B (zh) 一种改进的基于模板匹配的二次目标分群方法

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