-
通过分析红外导引头各部分抗干扰工作原理以及红外弹抗干扰过程,导引头固有属性、抗干扰性能指标以及导弹战技性能指标是影响红外空空导弹综合抗干扰特性的三个方面。结合这三个部分的具体内容构建红外抗干扰综合指标体系。
-
红外成像导引头主要由光学成像系统和随动系统组成,这两个系统性能参数的好坏会直接影响导引头的抗干扰性能。
光学成像系统[11]接收来自目标的红外辐射,形成数字灰度图像并对该图像进行处理,转换为电信号传递至随动系统。它的主要性能参数包括光通透率(Light Transmittance, LT)、焦距(Focal Length, FL)、视场角(Field of View, FOV)、扫描帧速(Scanning Frame Rate, SFR)、探测率(Detection Rate, DR)、图像处理周期(Image Processing Cycle, IPC)。
随动系统[11]接收光学成像系统处理过后的电信号,完成控制与跟踪信息和信号处理,驱动导引头光轴向弹目视线与光轴之间误差减小的方向运动,使其对准目标,实现实时跟踪。其主要性能参数包括最大跟踪角速度(Maximum Tracking Angular Velocity, MTAV)、框架角精度(Frame Angle Accuracy, FAA)和静态跟踪误差(Static Tracking Error, STE)。
-
红外成像导引头的抗干扰性能受自身固有属性以及干扰环境的影响很大,即便是同一导引头在不同干扰环境下其干扰性能也将会不同;即便是同一干扰环境下,用不同的评估指标设计出的导引头其性能也会有天壤之别。通过分析,红外导引头的抗干扰性能主要体现在干扰条件下的识别能力和稳定跟踪能力两方面。
(1)抗干扰识别能力
抗干扰识别能力是指红外成像导引头运用相应的目标识别算法,在打击目标的过程中,红外成像导引头从其所采集到的每一帧图像中识别出正确目标的能力[12]。其性能指标包括检测概率(Detection Probability, DP)、虚警概率(False Alarm Probability, FAP)、单帧识别时间(Single Frame Recognition Time, SFRT)、识别算法成功率(Success Rate of Recognition Algorithm, SRRA)。
检测概率(Pt)是识别算法在干扰环境下可以正确检测出目标的概率,定义为:
$$ {P_t} = {N_1}/N $$ (1) 式中:N1
为干扰环境下检测到目标的次数;N为导引头检测总次数。 虚警概率(Px)是导引头在识别过程中识别错误的概率,可表示为:
$$ {P_x} = {N_2}/N $$ (2) 式中:N2为红外导引头在识别过程中将非目标识别为目标的次数。
单帧识别时间是指导引系统从每一帧图像中提取出目标信息所需要的时间。单帧识别时间越短,留给制导系统的时间越多,导弹的抗干扰性能越好。
识别算法成功率是指在复杂干扰环境下,导引头的识别算法能正常运行并将目标成功识别的概率。
(2)稳定跟踪能力[13]
稳定跟踪能力涵盖了抗干扰跟踪能力和持续抗干扰能力,抗干扰跟踪能力可以用导引头跟踪目标时预测的水平视线角速度与真实水平视线角速度的误差,即目标水平视线角误差(Target Horizontal Line of Sight angle error, THLS)来体现,将其定义为:
$$ \eta = \dfrac{{\dfrac{1}{n}\displaystyle\sum\limits_{i = 1}^n {\left| {{{\dot q}_{HP}} - {{\dot q}_{HR}}} \right|} }}{{\dfrac{1}{n}\displaystyle\sum\limits_{i = 1}^n {\left| {{{\dot q}_{HR}}} \right|} }} $$ (3) 式中:
$\eta $ 为目标水平视线角误差;${{{\dot q}_{HP}}}$ 为跟踪目标时预测的水平视线角速度;${{{\dot q}_{HR}}}$ 为真实水平视线角速度。持续抗干扰能力可用导引头能经受的抗干扰时间占整个弹道时间的比值,即干扰占空比(Duty Ratio of Interference, DRI)来表示,定义为:
$$ \psi = {t_{\rm under}}/{t_{\rm tot}} $$ (4) 式中:
$\psi $ 为干扰占空比;${t_{\rm under}}$ 为导引头能经受的抗干扰时间;${t_{\rm tot}}$ 为弹道总时间。当$\psi $ 为1或接近于1时定义为全程抗干扰。 -
由参考文献[14]可知,弹目距离(Distance Range, DR)、杀伤能力(Kill Capacity, KC)是影响导弹抗干扰效果的战技性能指标。除此之外,导引头抗干扰性能评估还涉及攻击区的划分,根据导弹的进入角划分为尾后攻击区、侧向攻击区、前向攻击区以及禁止攻击区,所以将进入角(Entry Angle, EA)也作为导弹战技性能指标之一[15-17]。
综上所述,建立如图1所示的红外抗干扰综合评估指标体系。
-
首先将指标体系中的各因素划分为不同层次,递归的结构层次一般由三个或三个以上层次组成,红外抗干扰综合性能指标是本问题的预定目标(目标层)为A类;B类与C类指标都是影响目标实现的准则(准则层);D类与E类指标为促使目标实现的措施(措施层)。结构层次划分如图1所示。
-
划分结构层次后,需要对各层元素进行两两比较来构造出判别矩阵。传统判别矩阵的构建主要是基于德尔菲法,但该方法的评价结果在很大程度上取决于评价者的主观因素,主观性过强。故在德尔菲法的基础上引入主成分分析方法,利用试验数据求解指标之间的相关系数并对其排序,根据相关系数的大小来反映指标之间的相关程度,从而对判别矩阵进行修正,图2为PCA修正判别矩阵流程。
德尔菲法构建判别矩阵,对指标之间的重要程度采用九标度表示法进行量化,得到两两比较判别矩阵:
$${\boldsymbol{A}} = {({a_{fg}})_{n \times n}}$$ (5) 式中:
${a_{fg}}$ 表示指标f与指标g相对于上一层级指标重要性的标度量化值,其量化标度如表1所示。表 1 量化标度
Table 1. Quantitative scale
Comparison f and g afg agf f is as important as g 1 1 f is slightly more important than g 3 1/3 f is more important than g 5 1/5 f is obviously more important than g 7 1/7 f is absolutely more important than g 9 1/9 Adjacent judgment median 2,4,6,8 1/2,1/4,1/6,1/8 将构建好的判别矩阵利用主成分分析法[18]进行修正,具体操作步骤如下:
(1)构建原始数据矩阵X。假设有m条试验数据,每条数据都有n个指标,形成
$m \times n$ 的原始数据矩阵X,可表示为:$${\boldsymbol{X}}{\rm{ = }}\left[ {\begin{array}{*{20}{c}} {{x_{11}}}&{{x_{12}}}& \cdots &{{x_{1n}}} \\ {{x_{21}}}&{{x_{22}}}& \cdots &{{x_{2n}}} \\ \vdots & \vdots & \vdots & \vdots \\ {{x_{m1}}}&{{x_{m2}}}& \cdots &{{x_{mn}}} \end{array}} \right]$$ (6) 式中:
$i = 1,2, \cdots ,n$ 为指标数;$j = 1,2, \cdots ,m$ 为试验数据。(2)对原始数据标准化。由于构建的抗干扰评估指标具有不同的量纲与数量级,因此需要进行标准化才能进行比较与分析,故将X中的每一列(代表一个指标字段)进行零均值化处理,表达式为:
$$ \mu _{ji}^* = {x_{ji}} - \frac{1}{m}\sum\limits_{j = 1}^m {{x_{ji}}} \;\;\;\;\;i = 1,2, \cdots ,n;j = 1,2, \cdots ,m $$ (7) (3)求数据零均值化后的相关系数矩阵R,可表示为:
$${\boldsymbol{R}} = \left[ {\begin{array}{*{20}{c}} {{r_{11}}}&{{r_{12}}}& \cdots &{{r_{1n}}} \\ {{r_{21}}}&{{r_{22}}}& \cdots &{{r_{2n}}} \\ \vdots & \vdots & \vdots & \vdots \\ {{r_{n1}}}&{{r_{n2}}}& \cdots &{{r_{nn}}} \end{array}} \right]$$ (8) ${r_{uw}}$ 表示指标u与指标w之间的相关系数$( {u =} $ $ {1,2, \cdots ,n;w = 1,2, \cdots n})$ ,其求解公式为:$$ \left\{ {\begin{array}{*{20}{l}} {{r_{uw}}{\rm{ = }}\dfrac{1}{{m - 1}}\displaystyle\sum\limits_j^m {\left[ {\left( {\mu _{ju}^{*2}/{s_u}} \right) \cdot \left( {\mu _{jw}^{*2}/{s_w}} \right)} \right]} }\\ {{s_k} = \sqrt {\dfrac{1}{{m - 1}}\left[ {\displaystyle\sum\limits_{m = 1}^m {\mu _{mk}^{*2}} } \right]} } \end{array}} \right. $$ (9) (4)求取平均影响度
${R_i}$ [19-20]。由相关系数矩阵可以看出,有${r_{ii}}{\rm{ = }}1,{r_{ij}} = {r_{ji}}\left( {i = 1,2, \cdots n} \right)$ ,故令:$$ {R_i}{\rm{ = }}\sum\limits_{j = 1}^n {{r_{ji}} - 1} $$ (10) ${R_i}$ 反映了指标${X_i}$ 对其他指标的影响程度。其中${r_{ji}}$ 有可能出现负值,即这两指标呈负相关,在求解${R_i}$ 时则取$\left| {{r_{ji}}} \right|$ ,即${R_i}{\rm{ = }}\displaystyle\sum\limits_{j = 1}^n {\left| {{r_{ji}}} \right| - 1} $ 。最后对求解出的${R_i}$ 值进行排序,利用排序结果对判别矩阵的标度值进行观测,若标度值不符合排序结果,则采用如下规则进行修正:$$ \left\{ {\begin{array}{*{20}{c}} {{a_{{f_1}{g_1}}} + 1, {a_{{f_2}{g_1}}} - 1}\\ {\dfrac{1}{{{a_{{f_1}{g_1}}} - 1}}, \dfrac{1}{{{a_{{f_2}{g_1}}} + 1}}} \end{array}} \right.\;\;\;\;\;\;{a_{{f_1}{g_1}}} < {a_{{f_2}{g_1}}},{R_{f1}} > {R_{f2}} $$ (11) 若
${a_{31}}{\rm{ = }}2,{a_{41}}{\rm{ = }}3$ 而${R_3} > {R_4}$ ,则${a_{31}} + 1,{a_{41}} - 1$ 。 -
计算各层次指标权值可以等价为将主成分分析法修改过后的判别矩阵
${{\boldsymbol{A}}^ * }$ ,采用“求和法”来计算其相应的最大特征值与特征向量。(1)按列将判别矩阵
${{\boldsymbol{A}}^ * } = {(a_{fg}^ * )_{n \times n}}$ 进行归一化处理为:$$ {\overline {\left( {a_{fg}^*} \right)} _{n \times n}} \buildrel \Delta \over = {\left[ {\dfrac{{a_{fg}^*}}{{\displaystyle\sum\limits_{f = 1}^n {a_{fg}^*} }}} \right]_{n \times n}} $$ (12) (2)按行取总和:
$$ {\overline \omega _f}{\rm{ = }}\sum\limits_{g = 1}^n {a_{fg}^*} $$ (13) (3)对行和进行归一化处理,得到权重系数:
$$ {\omega _f} = \frac{{{{\overline \omega }_f}}}{{\displaystyle\sum\limits_{f = 1}^n {{{\overline \omega }_f}} }} $$ (14) (4)求最大特征值:
$$ {\lambda _{\max }}{\rm{ = }}\sum\limits_{g = 1}^n {\frac{{\left( {{A^*}\omega } \right)f}}{{n{\omega _f}}}} $$ (15) 则
${{\boldsymbol{A}}^ * }$ 最大特征值${\lambda _{\max }}$ 相应的特征向量W即为指标的权重向量。 -
计算一致性比例CR:
$$ CR = CI/RI $$ (16) 式中:
$CI = \dfrac{{{\lambda _{\max }} - n}}{{n - 1}}$ ;RI值通过查阅平均随机一致性指标表获得,如表2所示。表 2 平均随机一致性指标
Table 2. Mean random consistency index
Order (n) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 RI 0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 当CR<0.1时即可认为判别矩阵满足一致性,否则返回德尔菲法重新建立判别矩阵,循环进行上述步骤。
-
设第k−1层的
${n_{k - 1}}$ 个指标对红外导引系统的综合抗干扰性能的权重向量$ {\boldsymbol{ W}}^{( {k - 1} )} = {\bigg( W_1^{(k - 1)},} $ ${W_2^{(k - 1)},\cdots , } $ $ {W_{{n_k} - 1}^{(k - 1)} \cdots)^{\rm T}}$ ,记第k层上的${n_k}$ 个指标对第(k−1)层各指标的合成权重向量${{\boldsymbol{P}}^{\left( k \right)}} = {({P^{1(k)}},{P^{2(k)}},\cdots ,{P^{{n_{k - 1}}(k)}})_{{n_k} \times {n_{k - 1}}}}$ ,则k层指标对红外导引系统综合抗干扰性能的合成权重向量${{\boldsymbol{W}}^{(k)}}$ 为:$${{\boldsymbol{W}}^{(k)}} = {{\boldsymbol{P}}^{(k)}}{{\boldsymbol{W}}^{(k - 1)}}$$ (17) 递推可得:
$${{\boldsymbol{W}}^{(k)}} = {{\boldsymbol{P}}^{(k)}}{{\boldsymbol{P}}^{(k - 1)}}\cdots{{\boldsymbol{W}}^{(2)}}$$ (18) -
利用红外抗干扰仿真系统获取红外抗干扰数据,构成仿真数据矩阵
${\boldsymbol{B}}{\rm{ = }}\left[ {{B_1},{B_2},\cdots{B_k}} \right]$ ,其中${B_i}$ (i=1,2,···,k)为每一条m维的抗干扰仿真数据,采用Z-Score法对B开展标准化变换得到${{\boldsymbol{B}}^ * }$ ,${{\boldsymbol{B}}^ * }$ 与合成权重向量${{\boldsymbol{W}}^{(k)}}$ 相乘,获得每条抗干扰试验数据的评估结果,计算公式为:$$ b_{ij}^* = \frac{{{b_{ij}} - {{\bar B}_j}}}{{\sqrt {{\mathop{\rm var}} \left( {{B_j}} \right)} }} $$ (19) $$ E = {\boldsymbol{B}^ * } \cdot {\boldsymbol{W}^{\left( k \right)}}$$ (20) 式中:
$i = 1,2, \cdots,k;j = 1,2, \cdots, m$ ;${\bar B_j}$ 和${\mathop{\rm var}} \left( {{B_j}} \right)$ 分别为抗干扰数据j的均值与方差。 -
利用抗干扰评估仿真系统获取10条指标数据构建
$10 \times 18$ 的原始矩阵,指标数据如表3所示,其中D1为视场角(°)、D2为光通透率、D3为焦距(cm)、D4为扫描帧数(fps)、D5为探测率(1011 cm·Hz1/2·W−1)、D6为图像处理周期(ms)、D7为最大跟踪角速度((°)/s)、D8为框架角精度((°)/s)、D9为静态跟踪误差((°)/s)、D10为检测概率、D11为虚警概率、D12为单帧识别时间(s)、D13为识别算法成功率、D14为目标水平视线角误差(°)、D15为干扰占空比、E1为弹目距离(km)、E2为进入角(°)、E3为杀伤能力。对原始矩阵标准化后求取指标之间的相关系数矩阵,利用公式(10)求出指标的平均影响度,如表4所示。由表4可知,水平视线角误差(D14)对其他指标的影响程度最大,其次为图像处理周期(D6)与探测率(D5),视场角(D1)对其他指标的影响程度最小。通过平均影响度值的大小,根据公式(11)的规则修改判别矩阵。
表 3 指标数据
Table 3. Index data
No. D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 E1 E2 E3 1 10 0.965 27.3 2 4.9 15 14 0.3 0.05 0.774 0.075 0.001 0.92 0.18 1 4 10 0.9 2 7 0.987 27.2 2 5.7 16 13 0.2 0.04 0.595 0.016 0.0035 0.95 0.16 0.97 7 30 0.89 3 9 0.967 26.9 2 3.7 15 14 0.3 0.03 0.676 0.044 0.0022 0.91 0.21 0.82 6 60 0.93 4 8 0.975 26.8 2 4.3 15 15 0.4 0.05 0.88 0.059 0.0008 0.9 0.18 0.55 5 40 0.91 5 11 0.99 27.1 4 2.3 17 15 0.2 0.03 0.231 0.056 0.0028 0.93 0.19 0.73 4 90 0.88 6 6 0.985 27 4 2.8 20 14 0.3 0.04 0.702 0.05 0.0014 0.95 0.31 0.64 3 120 0.94 7 9 0.973 27.3 4 1.5 20 12 0.4 0.05 0.336 0.035 0.0045 0.93 0.84 0.46 8 150 0.87 8 12 0.983 26.9 2 3.3 19 13 0.3 0.03 0.498 0.026 0.0024 0.91 0.17 0.39 6 160 0.92 9 7 0.98 26.7 3 5.4 17 13 0.3 0.04 0.159 0.024 0.0032 0.92 0.2 0.84 5 70 0.93 10 8 0.95 27.2 4 4.7 15 14 0.2 0.05 0.967 0.047 0.0038 0.93 0.22 0.78 7 110 0.9 表 4 指标之间的平均影响度
Table 4. Average influence degree between indexes
Index D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 E1 E2 E3 Ri 3.248 3.564 4.019 4.350 5.804 5.847 4.809 3.676 3.829 4.558 5.307 4.955 3.909 6.039 5.009 3.914 5.605 3.893 -
通过德尔菲法构建的判别矩阵经过主成分分析法修正后,求解修正后的判别矩阵的最大特征值以及相对应的特征向量,并进行判别矩阵的一致性检验。
B类指标对总目标的判别矩阵以及相对权重如表5所示。
表 5 A~B判别矩阵及权重
Table 5. A-B discriminant matrix and weight
RA B1 B2 B3 WB B1 1 1/7 2 0.131 B2 7 1 9 0.793 B3 1/2 1/9 1 0.076 C类指标对B类指标的判别矩阵及相对权重如表6所示。
表 6 B1、B2~C判别矩阵及权重
Table 6. B1, B2-C discriminant matrix and weight
RB1 C1 C2 WC1 RB2 C3 C4 WC2 C1 1 2 0.667 C3 1 5 0.833 C2 1/2 1 0.333 C4 1/5 1 0.167 表 7 D~C1判别矩阵及权重
Table 7. D-C1 discriminant matrix and weight
RC1 D1 D2 D3 D4 D5 D6 WD1-6 D1 1 1/2 1/3 1/4 1/5 1/5 0.046 D2 2 1 1/2 1/3 1/4 1/4 0.070 D3 3 2 1 1/2 1/3 1/3 0.112 D4 4 3 2 1 1/2 1/2 0.180 D5 5 4 3 2 1 1 0.296 D6 5 4 3 2 1 1 0.296 表 8 D~C2判别矩阵及权重
Table 8. D-C2 discriminant matrix and weight
RC2 D7 D8 D9 WD7-9 D7 1 3 2 0.540 D8 1/3 1 1/2 0.163 D9 1/2 2 1 0.297 表 9 D~C3判别矩阵及权重
Table 9. D-C3 discriminant matrix and weight
RC3 D10 D11 D12 D13 WD10-13 D10 1 1/4 1/2 2 0.135 D11 4 1 2 7 0.524 D12 2 1/2 1 4 0.271 D13 1/2 1/7 1/4 1 0.070 表 10 D~C4判别矩阵及权重
Table 10. D-C4 discriminant matrix and weight
RC4 D14 D15 WD14-15 D14 1 3 0.75 D15 1/3 1 0.25 E类指标对B类指标的判别矩阵及相对权重如表11所示。
表 11 E~B3判别矩阵及权重
Table 11. E-B3 discriminant matrix and weight
RB3 E1 E2 E3 WE1-3 E1 1 1/3 1 0.2 E2 3 1 3 0.6 E3 1 1/3 1 0.2 对以上判别矩阵进行一致性检验,检验结果均通过利用公 式(17)、(18)计算措施层相对目标层的绝对权重W=[0.004 0.006 0.01 0.016 0.026 0.026 0.024 0.007 0.013 0.089 0.346 0.179 0.046 0.099 0.033 0.015 0.046 0.015]。
-
通过红外抗干扰仿真系统获取15组18维的抗干扰数据,用公式(19)对抗干扰数据标准化,用公式(20)对每条抗干扰数据进行评估打分,原始数据以及最后的评估结果如表12所示。其中抗干扰数据包括诱饵投掷参数、弹目态势参数及目标机动参数,仿真对抗结果为脱靶量MD(m),且以MD<10 m视为导弹命中。
表 12 红外抗干扰数据及评估结果
Table 12. Infrared anti-jamming data and evaluation results
Group Ns Ne Nf Tl/s Ti/s Te/s Vf/m·s-1 ${\theta _l}$/(°) ${\theta _v}$/(°) ${\theta _m}$/(°) Dr/km acc/m·s-2 thro/(°) Ma pitch/(°) yaw/(°) roll/(°) MD E Rank 1 12 1 4 0.5 0.1 0.2 10 0 −90 10 4 0 2.001 0.860 1.453 0 0 2.017 1.804 1 2 24 1 6 0.5 0.4 0.8 10 60 −90 130 6 0 2.001 0.860 1.453 0 0 13.713 0.710 10 3 36 1 4 0.5 0.1 0.4 20 60 −90 70 8 0 2.001 0.860 1.453 0 0 187.404 0.641 13 4 12 2 6 2 0.4 0.8 40 150 90 30 5 15.526 2.019 0.724 2.236 −55.147 15.557 4.560 1.698 2 5 48 2 4 2.5 0.1 0.6 40 30 −60 90 4 15.526 2.019 0.724 2.236 −55.147 15.557 13.546 0.721 9 6 24 2 12 2 0.1 0.8 40 150 90 90 8 15.526 2.019 0.724 2.236 −55.147 15.557 450.930 0.361 15 7 12 1 6 0.5 0.1 0.8 30 60 −90 70 8 −24.523 2.044 0.939 1.141 77.102 −24.551 5.321 1.472 3 8 36 1 4 0.5 0.4 0.8 10 60 −90 10 6 −24.523 2.044 0.939 1.141 77.102 −24.551 12.705 0.798 8 9 24 1 12 1 0.1 0.2 20 120 −90 130 10 −24.523 2.044 0.939 1.141 77.102 −24.551 370.244 0.391 14 10 48 2 4 1.5 0.1 0.8 10 0 60 140 7 −37.325 2.677 0.970 14.495 0.909 −37.308 11.973 0.961 7 11 48 2 12 1.5 0.4 0.6 10 0 60 80 8 −37.325 2.677 0.970 14.495 0.909 −37.308 169.682 0.659 12 12 12 2 4 1.5 0.1 0.2 20 0 60 20 6 −37.325 2.677 0.970 14.495 0.909 −37.308 7.475 1.280 5 13 12 3 12 4 0.3 0.6 30 120 90 40 5 14.294 1.711 1.187 −10.307 −0.005 14.308 6.508 1.342 4 14 36 3 6 4 0.5 1 30 120 90 160 9 14.294 1.711 1.187 −10.307 −0.005 14.308 9.277 1.213 6 15 48 3 6 4 0.5 1 30 120 90 100 8 14.294 1.711 1.187 −10.307 −0.005 14.308 14.220 0.697 11 诱饵投掷参数包括:投掷总数
${N_s}$ (个)、每次投掷个数${N_e}$ (个)、每组投掷个数${N_f}$ (个)、投掷时机${T_l}$ (s)、投掷组内间隔${T_i}$ (s)、投掷组间间隔${T_e}$ (s)、投掷速度${V_f}$ (m/s)、水平投掷角度${\theta _l}$ (°)、垂直投掷角度${\theta _v}$ (°);弹目态势参数包括:进入角${\theta _m}$ (°)、弹目距离${D_r}$ (km);目标机动参数包括:平均加速度aac(m/s2)、平均油门开度thro(°)、平均马赫数Ma、平均俯仰角pitch(°)、平均偏航角yaw(°)、平均滚转角roll(°)。从表12中可以看出,不同的诱饵投掷策略、弹目态势以及目标机动方式所构成的红外干扰环境对红外系统的抗干扰能力会产生不同的影响。将最后的评估结果进行排序可以看出,组别1、4、7、12、13、14的脱靶量均小于10 m,目标完全毁伤或失去作战能力,相对应的评估结果均大于1.0,即在此干扰环境下红外导弹的抗干扰性能较好;而其余组别的脱靶量均大于10 m且评估结果均小于1.0,即可认为在此干扰环境下红外导弹的抗干扰性能较差,导弹脱靶。
将每组红外抗干扰数据的脱靶量与相应评估打分结果绘制如图3所示的对比图,从图3中可以更加直观地看出脱靶量小于10 m,评估打分结果大于1.0;脱靶量大于10 m,打分结果小于1.0。同时还可以看出脱靶量大小与评估打分结果大致呈负相关,即脱靶量越大,相应的评估打分越低。
同时将未用主成分分析法修正的判别矩阵进行正常的层次分析计算权重,用同样的红外抗干扰数据进行评估打分,结果如图4所示。从图4中可以看出,第5、6、15组试验,导弹脱靶量大于10 m,但评分结果均在E=1.0以上;而第7、12组试验,导弹脱靶量小于10 m,但评分却低于E=1.0。
Research on infrared anti-jamming evaluation method based on principal component-analytic hierarchy process
-
摘要: 为了提高红外抗干扰评估的客观性以及准确性,针对层次分析法评估红外抗干扰能力过程中构建判别矩阵主观因素过强的问题,提出了一种基于主成分-层次分析(analytic hierarchy process, AHP)的红外抗干扰评估方法。首先利用主成分分析(principal component analysis, PCA)法求解指标之间的相关系数矩阵,提出平均影响度这一指标;然后通过平均影响度大小,根据既定的调整规则调整判别矩阵;最后利用红外抗干扰仿真数据,结合调整后的判别矩阵进行层次分析,建立主成分-层次分析红外抗干扰评估模型。仿真结果表明:该模型可以较好地评价相应干扰环境下红外对抗的评估结果,使得对抗结果与评估打分呈负相关且打分结果以均线E=1.0为界,导弹命中则评估结果均在均线以上,可以较为直观地观察在该对抗条件下的红外抗干扰评估结果。Abstract: In order to improve the objectivity and accuracy of infrared anti-jamming evaluation, aiming at the problem that the subjective factors of discriminant matrix are too strong in the process of evaluating infrared anti-jamming ability by analytic hierarchy process (AHP), an infrared anti-jamming evaluation method based on principal component-analytic hierarchy process was proposed. Firstly, the correlation coefficient matrix between indexes was solved by principal component analysis (PCA), and the index of average influence degree was put forward. Then, through this index value, the discriminant matrix was adjusted according to the established adjustment rules. Finally, using the infrared anti-jamming simulation data, combined with the adjusted discriminant matrix for AHP, the principal component-analytic hierarchy process infrared anti-jamming evaluation model was established. The simulation results show that the model can better evaluate the evaluation results of infrared countermeasure in the corresponding jamming environment, there is a negative correlation between the confrontation result and the evaluation score. The score is bounded by the E=1.0, and the evaluation results of missile hit are all above the average. The evaluation result of infrared anti-jamming under this countermeasure condition can be observed intuitively.
-
表 1 量化标度
Table 1. Quantitative scale
Comparison f and g afg agf f is as important as g 1 1 f is slightly more important than g 3 1/3 f is more important than g 5 1/5 f is obviously more important than g 7 1/7 f is absolutely more important than g 9 1/9 Adjacent judgment median 2,4,6,8 1/2,1/4,1/6,1/8 表 2 平均随机一致性指标
Table 2. Mean random consistency index
Order (n) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 RI 0 0 0.52 0.89 1.12 1.26 1.36 1.41 1.46 1.49 1.52 1.54 1.56 1.58 表 3 指标数据
Table 3. Index data
No. D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 E1 E2 E3 1 10 0.965 27.3 2 4.9 15 14 0.3 0.05 0.774 0.075 0.001 0.92 0.18 1 4 10 0.9 2 7 0.987 27.2 2 5.7 16 13 0.2 0.04 0.595 0.016 0.0035 0.95 0.16 0.97 7 30 0.89 3 9 0.967 26.9 2 3.7 15 14 0.3 0.03 0.676 0.044 0.0022 0.91 0.21 0.82 6 60 0.93 4 8 0.975 26.8 2 4.3 15 15 0.4 0.05 0.88 0.059 0.0008 0.9 0.18 0.55 5 40 0.91 5 11 0.99 27.1 4 2.3 17 15 0.2 0.03 0.231 0.056 0.0028 0.93 0.19 0.73 4 90 0.88 6 6 0.985 27 4 2.8 20 14 0.3 0.04 0.702 0.05 0.0014 0.95 0.31 0.64 3 120 0.94 7 9 0.973 27.3 4 1.5 20 12 0.4 0.05 0.336 0.035 0.0045 0.93 0.84 0.46 8 150 0.87 8 12 0.983 26.9 2 3.3 19 13 0.3 0.03 0.498 0.026 0.0024 0.91 0.17 0.39 6 160 0.92 9 7 0.98 26.7 3 5.4 17 13 0.3 0.04 0.159 0.024 0.0032 0.92 0.2 0.84 5 70 0.93 10 8 0.95 27.2 4 4.7 15 14 0.2 0.05 0.967 0.047 0.0038 0.93 0.22 0.78 7 110 0.9 表 4 指标之间的平均影响度
Table 4. Average influence degree between indexes
Index D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 E1 E2 E3 Ri 3.248 3.564 4.019 4.350 5.804 5.847 4.809 3.676 3.829 4.558 5.307 4.955 3.909 6.039 5.009 3.914 5.605 3.893 表 5 A~B判别矩阵及权重
Table 5. A-B discriminant matrix and weight
RA B1 B2 B3 WB B1 1 1/7 2 0.131 B2 7 1 9 0.793 B3 1/2 1/9 1 0.076 表 6 B1、B2~C判别矩阵及权重
Table 6. B1, B2-C discriminant matrix and weight
RB1 C1 C2 WC1 RB2 C3 C4 WC2 C1 1 2 0.667 C3 1 5 0.833 C2 1/2 1 0.333 C4 1/5 1 0.167 表 7 D~C1判别矩阵及权重
Table 7. D-C1 discriminant matrix and weight
RC1 D1 D2 D3 D4 D5 D6 WD1-6 D1 1 1/2 1/3 1/4 1/5 1/5 0.046 D2 2 1 1/2 1/3 1/4 1/4 0.070 D3 3 2 1 1/2 1/3 1/3 0.112 D4 4 3 2 1 1/2 1/2 0.180 D5 5 4 3 2 1 1 0.296 D6 5 4 3 2 1 1 0.296 表 8 D~C2判别矩阵及权重
Table 8. D-C2 discriminant matrix and weight
RC2 D7 D8 D9 WD7-9 D7 1 3 2 0.540 D8 1/3 1 1/2 0.163 D9 1/2 2 1 0.297 表 9 D~C3判别矩阵及权重
Table 9. D-C3 discriminant matrix and weight
RC3 D10 D11 D12 D13 WD10-13 D10 1 1/4 1/2 2 0.135 D11 4 1 2 7 0.524 D12 2 1/2 1 4 0.271 D13 1/2 1/7 1/4 1 0.070 表 10 D~C4判别矩阵及权重
Table 10. D-C4 discriminant matrix and weight
RC4 D14 D15 WD14-15 D14 1 3 0.75 D15 1/3 1 0.25 表 11 E~B3判别矩阵及权重
Table 11. E-B3 discriminant matrix and weight
RB3 E1 E2 E3 WE1-3 E1 1 1/3 1 0.2 E2 3 1 3 0.6 E3 1 1/3 1 0.2 表 12 红外抗干扰数据及评估结果
Table 12. Infrared anti-jamming data and evaluation results
Group Ns Ne Nf Tl/s Ti/s Te/s Vf/m·s-1 ${\theta _l}$ /(°)${\theta _v}$ /(°)${\theta _m}$ /(°)Dr/km acc/m·s-2 thro/(°) Ma pitch/(°) yaw/(°) roll/(°) MD E Rank 1 12 1 4 0.5 0.1 0.2 10 0 −90 10 4 0 2.001 0.860 1.453 0 0 2.017 1.804 1 2 24 1 6 0.5 0.4 0.8 10 60 −90 130 6 0 2.001 0.860 1.453 0 0 13.713 0.710 10 3 36 1 4 0.5 0.1 0.4 20 60 −90 70 8 0 2.001 0.860 1.453 0 0 187.404 0.641 13 4 12 2 6 2 0.4 0.8 40 150 90 30 5 15.526 2.019 0.724 2.236 −55.147 15.557 4.560 1.698 2 5 48 2 4 2.5 0.1 0.6 40 30 −60 90 4 15.526 2.019 0.724 2.236 −55.147 15.557 13.546 0.721 9 6 24 2 12 2 0.1 0.8 40 150 90 90 8 15.526 2.019 0.724 2.236 −55.147 15.557 450.930 0.361 15 7 12 1 6 0.5 0.1 0.8 30 60 −90 70 8 −24.523 2.044 0.939 1.141 77.102 −24.551 5.321 1.472 3 8 36 1 4 0.5 0.4 0.8 10 60 −90 10 6 −24.523 2.044 0.939 1.141 77.102 −24.551 12.705 0.798 8 9 24 1 12 1 0.1 0.2 20 120 −90 130 10 −24.523 2.044 0.939 1.141 77.102 −24.551 370.244 0.391 14 10 48 2 4 1.5 0.1 0.8 10 0 60 140 7 −37.325 2.677 0.970 14.495 0.909 −37.308 11.973 0.961 7 11 48 2 12 1.5 0.4 0.6 10 0 60 80 8 −37.325 2.677 0.970 14.495 0.909 −37.308 169.682 0.659 12 12 12 2 4 1.5 0.1 0.2 20 0 60 20 6 −37.325 2.677 0.970 14.495 0.909 −37.308 7.475 1.280 5 13 12 3 12 4 0.3 0.6 30 120 90 40 5 14.294 1.711 1.187 −10.307 −0.005 14.308 6.508 1.342 4 14 36 3 6 4 0.5 1 30 120 90 160 9 14.294 1.711 1.187 −10.307 −0.005 14.308 9.277 1.213 6 15 48 3 6 4 0.5 1 30 120 90 100 8 14.294 1.711 1.187 −10.307 −0.005 14.308 14.220 0.697 11 -
[1] Li Likun. Present situation and development direction of precision guidance technology [J]. Aero Weaponry, 2004(1): 1-4. (in Chinese) doi: 10.3969/j.issn.1673-5048.2004.01.001 [2] 高卫 孙奕帆. 精确制导武器系统电子干扰效果试验与评估[M]. 北京: 国防工业出版社, 2018: 119-123. Gao Wei, Sun Yifan. Test and Evaluation of Electronic Interference Effects of Precision Guided Weapon Systems[M]. Beijing: National Defense Industry Press, 2018: 119-123. (in Chinese) [3] Wang Xiaoming, Wang Mei. ECCM test technology for antiaircraft missile weapon system [J]. Aerospace Shanghai, 2013, 30(2): 34-38. (in Chinese) doi: 10.3969/j.issn.1006-1630.2013.02.008 [4] Forrai D P, Maier J J. Generic models in the advanced IRCM assessment model [C]//Simulation Conference, Proceedings of the Winter, IEEE, 2001: 789-796 [5] Surng-Gabb Jahng, Hong Hyun-Ki, Han Sung-Hyun. Dynamic simulation of the rosette scanning infrared seeker and an infrared counter-countermeasure using the moment technique [J]. Optical Engineering, 1999, 38(5): 921-928. doi: 10.1117/1.602051 [6] Hong Hyun-Ki. Adaptive infrared counter-counterm-easure for two-color spinning concentric-annular-ring reticle seeker [J]. Optical Engineering, 2001, 40(6): 1093-1099. doi: 10.1117/1.1367344 [7] Surng-Gabb Jahng, Hong Hyun-Ki, Seo Dongsun, et al. New infrared counter-counter measure technique using an iterative self-organizing data analysis algorithm for the rosette scanning infrared seeker [J]. Optical Engineering, 2000, 39(9): 2397-2404. [8] 刘泽胤. 基于DODAF的系统效能评估[D]. 哈尔滨工程大学, 2008. Liu Zeyin. System effectiveness evaluation based on DODAF [D]. Harbin: Harbin University of Engineering, 2008. (in Chinese) [9] Wang Quan, Dong Weihao, Liu Xin’ai, et al. Analysis of evaluation for new infrared air-to-air missile countermeasures performance [J]. Aerospace Electronic Warfare, 2019, 35(3): 16-19. (in Chinese) [10] Saaty T L. A scaling method for priorities in hierarchical structures [J]. Journal of Mathematical Psychology, 1977, 15(3): 234-281. [11] 杨宜禾, 岳敏, 周维真. 红外系统(第2版)[M]//红外系统(第2版). 国防工业出版社, 1995: 23-70. Yang Yihe, Yue Min, Zhou Weizhen. Infrared System [M]. 2nd Edition. Beijing: National Defense Industry Press, 1995: 23-70. (in Chinese) [12] Zhou Yangheng, Cui Shaohui, Fang Dan. Construction of antijamming performance index system of infrared imaging seeker [J]. Modern Defence Technology, 2019, 47(3): 175-180. (in Chinese) doi: 10.3969/j.issn.1009-086x.2019.03.26 [13] Tang Shanjun, Wang Feng, Chen Xiaodong. Anti-jamming performance index system and evaluation of infrared missile [J]. Aerospace Shanghai, 2017(4): 144-149. (in Chinese) [14] Li Shenbo, Li Taorui, Tong Zhongxiang, et al. Influence of tactical parameters of missile on jamming effectiveness of surface source infrared decoy [J]. Infrared and Laser Engineering, 2018, 47(7): 0704004. (in Chinese) [15] Dong Minzhou, Yu Yunfeng, Yan Jie. AHP based on exponential scale for better evaluation of IR seeker’s anti-jamming performance during missile research process [J]. Journal of Northwestern Polytechnical University, 2008, 6(1): 69-73. (in Chinese) doi: 10.3969/j.issn.1000-2758.2008.01.014 [16] Xu Yang, Fang Yangwang, Wu Youli, et al. Association rule mining for the infrared countermeasure by the pf-growth algorithm[C]//Proceedings of the 37th China Control Conference, 2018: 8043-8048. [17] 韩培骏. 红外导引系统抗干扰性能评估准则与方法研究[D]. 南京航空航天大学, 2012. Han Peijun. Research on the anti-interference performance evaluation criterion and method of infrared guidance system[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2012. (in Chinese) [18] Liang Shiming. Improvement of Saaty's analytic hierarchy process [J]. Engineering Mathematics, 1995(2): 140-141. (in Chinese) [19] 孙敏. 基于层次分析法的红外制导系统抗干扰性能评估指标分析[D]. 华中科技大学, 2017. Sun Min. Analysis of anti-jamming performance evaluation index of infrared guidance system based on analytic hierarchy process[D]. Wuhan: Huazhong University of Science & Technology, 2017. (in Chinese) [20] Li Lingjuan, Dou Kun. Research on the consistency of the judgment matrix in AHP [J]. Computer Technology and Development, 2009, 19(10): 131-133. (in Chinese) doi: 10.3969/j.issn.1673-629X.2009.10.036 [21] Liu Lijia, Hu Jianwang, Sun Huixian. Adjustment method of comparison matrix in analytic hierarchy process [J]. Journal of Weapons and Equipment Engineering, 2020, 41(2): 221-224. (in Chinese) doi: 10.11809/bqzbgcxb2020.02.044 [22] 周志华. 机器学习[M]. 清华大学出版社, 2016: 229-232. Zhou Zhihua. Machine Learning[M]. Beijing: Tsinghua University Press, 2016: 229-232. (in Chinese) [23] Qu Dongsheng, Liu Lianwei, Fan Hongjie, et al. Evaluation method of infrared system countermeasures performance based on principal components analysis [J]. Laser & Infrared, 2020, 50(6): 712-717. (in Chinese) doi: 10.3969/j.issn.1001-5078.2020.06.011