基于自适应目标函数的协同干扰策略分配

ＤＯＩ：１０．１１９９１／ｙｙｋｊ２０１８１１０１８

应　　　用　　　科　　　技ＡｐｐｌｉｅｄＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ

Ｖｏｌ．４６№．４Ｊｕｌ．２０１９

基于自适应目标函数的协同干扰策略分配

汲清波，杨　帅

哈尔滨工程大学，信息与通信工程学院，黑龙江哈尔滨１５０００１

摘　要：在多机伴随突防组网雷达过程中，会受到不同雷达的威胁。多机伴随协同干扰在保证目标成功突防中起重要作用，而如何提高己方干扰资源的利用效率是协同干扰的核心。文中提出协同干扰效果是关于距离的函数，如果单纯以某一单一目标函数作为干扰效果评估指标，会大大影响干扰机的利用效率，因此，在协同突防编队由远及近的突防过程中，将该过程划分为搜索与定位两个阶段，分别以检测概率和几何精度因子（ＧＤＯＰ）建立各阶段的干扰评估指标，随着距离的变化，干扰效果评估指标也在逐渐变化，通过改进的粒子群算法对每个阶段的目标函数求解寻最优，既解决了传统的动态规划算法计算量过大的问题，又实现了算法的快速收敛和寻优。通过ＭＡＴＬＡＢ仿真结果分析，该算法能够快速地收敛并且提扰机的利用效率。

关键词：组网雷达；协同干扰；检测概率；几何精度因子；自适应目标函数；干扰资源分配

中图分类号：ＴＮ９１．７　　　　　　文献标志码：Ａ　　　　　　文章编号：１００９⁃６７１Ｘ（２０１９）０４⁃０７０⁃０７

Ｃｏｏｐｅｒａｔｉｖｅｉｎｔｅｒｆｅｒｅｎｃｅｓｔｒａｔｅｇｙａｌｌｏｃａｔｉｏｎｂａｓｅｄ

ｏｎａｄａｐｔｉｖｅｏｂｊｅｃｔｉｖｅｆｕｎｃｔｉｏｎ

ＪＩＱｉｎｇｂｏ，ＹＡＮＧＳｈｕａｉ

ＣｏｌｌｅｇｅｏｆＩｎｆｏｒｍａｔｉｏｎａｎｄＣｏｍｍｕｎｉｃａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＨａｒｂｉｎＥｎｇｉｎｅｅｒｉｎｇＵｎｉｖｅｒｓｉｔｙ，Ｈａｒｂｉｎ１５０００１，Ｃｈｉｎａ

Ａｂｓｔｒａｃｔ：Ｉｎｔｈｅｐｒｏｃｅｓｓｏｆｍｕｌｔｉ⁃ｍａｃｈｉｎｅａｄｊｏｉｎｉｎｇｐｅｎｅｔｒａｔｉｏｎｒａｄａｒ，ｒａｄａｒｗｉｌｌｂｅｔｈｒｅａｔｅｎｅｄｂｙｄｉｆｆｅｒｅｎｔｓｙｓｔｅｍｓ．Ｍｕｌｔｉ－ａｉｒｃｒａｆｔｃｏｎｃｏｍｉｔａｎｔｃｏｏｐｅｒａｔｉｖｅｊａｍｍｉｎｇｐｌａｙｓａｎｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｇｕａｒａｎｔｅｅｉｎｇｔｈｅｓｕｃｃｅｓｓｆｕｌｐｅｎｅｔｒａｔｉｏｎｏｆｔａｒｇｅｔｓ，ａｎｄｈｏｗｔｏｅｆｆｅｃｔｉｖｅｌｙｕｔｉｌｉｚｅｉｔｓｏｗｎｊａｍｍｉｎｇｒｅｓｏｕｒｃｅｓｉｓｔｈｅｃｏｒｅｏｆｃｏｏｐｅｒａｔｉｖｅｊａｍｍｉｎｇ．Ｔｈｉｓｐａｐｅｒｐｒｏ⁃ｐｏｓｅｓｔｈａｔｔｈｅｃｏｏｐｅｒａｔｉｖｅｉｎｔｅｆｅｒｅｎｃｅｅｆｆｅｃｔｉｓａｆｕｎｃｔｉｏｎａｂｏｕｔｄｉｓｔａｎｃｅ．Ｉｔｗｉｌｌｈｅａｖｉｌｙｉｎｆｌｕｅｎｃｅｔｈｅｕｔｉｌｉｚａｔｉｏｎｅｆｆｉ⁃ｃｉｅｎｃｙｏｆｔｈｅｉｎｔｅｒｆｅｒｅｃｅｍａｃｈｉｎｅｉｆｕｓｉｎｇａｓｉｎｇｌｅｏｂｊｅｃｔｉｖｅｆｕｎｃｔｉｏｎａｓｔｈｅｅｖａｌｕａｔｉｏｎｉｎｄｅｘｏｆｉｎｔｅｒｆｅｒｅｎｃｅｅｆｆｅｃｔ．Ｔｈｅｒｅｆｏｒｅ，ｉｎｔｈｅｐｒｏｃｅｓｓｏｆｃｏｏｐｅｒａｔｉｖｅｐｅｎｅｔｒａｔｉｏｎｆｏｒｍａｔｉｏｎｆｒｏｍｆａｒｔｏｎｅａｒ，ｔｈｅｐｅｎｅｔｒａｔｉｏｎｐｒｏｃｅｓｓｉｓｄｉｖｉｄｅｄｉｎ⁃ｔｏｔｗｏｓｔａｇｅｓ：ｓｅａｒｃｈｉｎｇａｎｄｌｏｃａｔｉｎｇ．Ｔｈｅｄｅｔｅｃｔｉｏｎｐｒｏｂａｂｉｌｉｔｙａｎｄｇｅｏｍｅｔｒｉｃｐｒｅｃｉｓｉｏｎｆａｃｔｏｒ（ＧＤＯＰ）ａｒｅｕｓｅｄｒｅｓｐｅｃｔｉｖｅｌｙｔｏｅｓｔａｂｌｉｓｈｔｈｅｉｎｔｅｒｆｅｒｅｎｃｅｅｖａｌｕａｔｉｏｎｉｎｄｅｘｏｆｅｖｅｒｙｓｔａｇｅ．Ｔｈｅｉｎｔｅｒｆｅｒｅｎｃｅｅｆｆｅｃｔｅｖａｌｕａｔｉｏｎｉｎｄｅｘｔｈｅｏｐｔｉｍａｌｓｏｌｕｔｉｏｎｔｏｔｈｅｏｂｊｅｃｔｉｖｅｆｕｎｃｔｉｏｎｉｎｅａｃｈｓｔａｇｅ，ｗｈｉｃｈｎｏｔｏｎｌｙｓｏｌｖｅｓｔｈｅｐｒｏｂｌｅｍｏｆｅｘｃｅｓｓｉｖｅｃｏｍｐｕｔａ⁃ｔｉｏｎｏｆｔｈｅｔｒａｄｉｔｉｏｎａｌｄｙｎａｍｉｃｐｒｏｇｒａｍｍｉｎｇａｌｇｏｒｉｔｈｍ，ｂｕｔａｌｓｏｒｅａｌｉｚｅｓｔｈｅｃｏｎｖｅｒｇｅｎｃｅａｎｄｏｐｔｉｍｉｚａｔｉｏｎｏｆｔｈｅａｌ⁃ｑｕｉｃｋｌｙ，ａｎｄｃａｎｉｎｃｒｅａｓｅｕｔｉｌｉｚａｔｉｏｎｅｆｆｉｃｉｅｎｃｙｏｆｔｈｅｉｎｔｅｒｆｅｒｅｎｃｅｍａｃｈｉｎｅ．ｆｕｎｃｔｉｏｎ；ｉｎｔｅｒｆｅｒｅｎｃｅｒｅｓｏｕｒｃｅｓａｌｌｏｃａｔｉｏｎ

ｃｈａｎｇｅｓｗｉｔｈｔｈｅｄｉｓｔａｎｃｅｃｈａｎｇｉｎｇｓｔｅｐｂｙｓｔｅｐ．Ｔｈｅｉｍｐｒｏｖｅｄｐａｒｔｉｃｌｅｓｗａｒｍｏｐｔｉｍｉｚａｔｉｏｎａｌｇｏｒｉｔｈｍｉｓｕｓｅｄｔｏｆｉｎｄｇｏｒｉｔｈｍｍｏｒｅｑｕｉｃｋｌｙ．ＴｈｒｏｕｇｈｔｈｅａｎａｌｙｓｉｓｏｆｓｉｍｕｌａｔｉｏｎｒｅｓｕｌｔｓｏｆＭＡＴＬＡＢ，ｉｔｐｒｏｖｅｓｔｈａｔｔｈｅａｌｇｏｒｉｔｈｍｃｏｎｖｅｒｇｅｓＫｅｙｗｏｒｄｓ：ｃｏｏｐｅｒａｔｉｖｅｎｅｔｔｅｄｒａｄａｒ；ｃｏｏｐｅｒａｔｉｖｅｉｎｔｅｒｆｅｒｅｎｃｅ；ｄｅｔｅｃｔｉｏｎｐｒｏｂａｂｉｌｉｔｙ；ＧＤＯＰ；ａｄａｐｔｉｖｅｏｂｊｅｃｔｉｖｅ　　随着现代电子技术的发展以及组网雷达的应用

收稿日期：２０１９－１１－１８．　　

基金项目：总装预研重点基金项目（６１４０４１５０１０１）；上海航天科技创

新基金项目（ＳＡＳＴ２０１７－０６８）．

作者简介：汲清波，女，副教授，博士；

杨帅，男，硕士研究生．

通信作者：杨帅，Ｅ－ｍａｉｌ：１２８９９０８８５１＠ｑｑ．ｃｏｍ．

单一干扰资源往往无法成功突防雷达网保护的区

域。在现代战场环境中，多干扰机通过频域、极化域等协同组网干扰是保证编队突防成功与否以及打击摧毁目标的重要手段，而如何有效利用己方干扰资源取得最优的干扰效果成为了研究的重要课题，亦

第４期　　　　　　　　　　　　汲清波，等：基于自适应目标函数的协同干扰策略分配　　　　　　　·７１·

是本文研究的目的。

在协同干扰资源分配领域，众多学者根据雷达侦察检测目标的原理提出了不同的干扰评估准则，并且再此基础上提出了相应分配模型，但目前为止，尚未有干扰分配模型得到众多学者的一致认同。文献［２］以干扰暴露区作为干扰评估准则建立目标函数，利用多级动态优化算法寻求干扰效果最优值，但该算法采用穷举法计算函数值，计算量随着干扰资源以及组网雷达数量的增加而增加；文献［３］将博弈论应用于动态干扰评估，但该算法只考虑了一对一干Ｈｏｐｆｉｅｌｄ扰，并未考虑协同干扰；文献［４］首次将域，具有较快的收敛速度神经网络算法应用于干扰资源优化分配领，但如何选择步径仍亟待研究；文献［５－７］分别采用基于改进离散布谷鸟算法、免疫算法、蚁群算法对雷达干扰资源进行分配，尽管分配结果收敛速度快，但仍然只是以单一评估准则作为目标函数和分配依据，未能将飞机编队与攻击目标的距离作为依据而采用不同的目标函数作为分配依据。

在多机编队突防过程中，由于组网雷达工作的多样性［８］，飞机编队与目标的距离由远及近过程中往往会经历搜索同，例如在距离较远时。处在不同阶段的雷达威胁等级并不相、定位、跟踪、火力引导等４个阶段［９］，组网雷达的任务以搜索目标

为主，即雷达的任务以发现目标是否存在为主要目的，此时，干扰的任务分配也应以降低雷达探测目标的发现概率为主，如果在此阶段以定位精度为目标函数并不能使己方的干扰资源达到最优使用效率。因此，本文将距离作为主要考量，并划分不同距离下的雷达工作阶段，在距离较远时，雷达处于警戒搜索阶段，该阶段主要以目标发现概率作为干扰效果评估准则；随着距离的逐渐缩小，雷达发现目标的概率逐渐增大，此时，雷达从搜索阶段转入定位阶段，该阶段以几何精度因子ＧＤＯＰ）主要利用粒子群算法采用自适应目标函数进行资源

为目标函数进行资源分配（ｇｅｏｍａｔｒｉｃｄｉｌｕｔｉｏｎ。上述两个阶段ｐｒｅｃｉｓｉｏｎ，分配，既保证了收敛速度又有效解决了处于不同阶段的雷达因为任务的不同而导致的资源无法充分利用的问题。

１　干扰资源分配模型以及目标函数的

建立

１．１　干扰资源分配是在侦察分选雷达信号之后雷达威胁等级的确定

，根

据侦（ＰＷ）、察得（ＰＲＦ）带到宽的（ＢＷ雷达）、脉载冲频描（述ＲＦ字）、（脉ＰＤＷ），冲重复即脉频宽定不同的威胁等级等信息，以与目标的距离划分的不同阶段制率。在搜索阶段，飞机编队与攻击目标距离较远，远程警戒搜索雷达的威胁较大。远程搜索警戒雷达由于探测距离远，以发现目标的存在与否为主要目的，因此，其载频通常较低，带宽较窄，在该阶段，低载频的雷达威胁程度更高。在雷达网发现目标之后即转入定位阶段，雷达以精确定位为主，此时，ＰＲＦ越高，速度分辨率越高，ＢＷ越宽，距离分辨率越高，因此，具有该指标的雷达威胁等级也越高。由上述分析可知，在同一距离上，具有不同特征的雷达威胁等级不同；在不同距离上，具有相同特征的同一部雷达由于距离的差异威胁等级也不同。

雷达威胁等级的确定对于能否有效分配资源起重要作用。本文根据文献［１０］的熵权法确定ＰＤＷ各指标的权重，最后确定各雷达在不同距离上的威胁等级，从而有效地制定干扰策略，使得干扰效果达１．２　到最优目标函数的确定。

１．２．１　众所周知搜索阶段目标函数，雷达由于信号功率－检测概率

、灵敏度等，

使其工作距离有限。当雷达与目标的距离超出其最大工作距离时，信号的强度远低于雷达工作灵敏度而无法探测到目标，但是随着雷达与目标之间距离减小，雷达开始搜索到目标并逐渐转换为对目标的跟踪阶段，因此，在搜索阶段需要选择合适的目标函数为雷达资源分配提供依据，本文选择目标检测概率作为依据进行分配。

雷达方程与干扰方程可以分别用式（１）和（２）表示：

ＰＰｔＧ２ｔσλ

２ｒｓ＝ＰＰ（４π）３Ｒ４ｔ

ｊＧｊＧｒｊ＝

雷达的信干比为：

（４π）２ｔ（θ）Ｒλ２ξｊ

２ｊ

（１）ＳＮＲ＝

ＰＰＰｒｓ＝

ｒｊ

４πＰｔＧｔσＲ２ｊ

ｊＧｊＧｊ（θ）ξｊＲ４ｔ

（２）

式中：Ｐｒｓ表示雷达接收机接收到的目标回波信号功率；Ｐｒｊ是雷达接收到的干扰信号功率；Ｐｔ和Ｐｊ是雷达接收天线口面接收到的干扰信号功率；Ｐｉｊ是雷达接收天线口面接收到的干扰信号功率；Ｐｔ和Ｐｊ则分别表示雷达发射机和干扰发射机发射的发射信号功

·７２·应　　　用　　　科　　　技　　　　　　　　　　　　　　　　第４６卷

率；Ｇｔ和Ｇｊ分别表示雷达、干扰机的接收天线增益；σ表示目标的平均散射截面积θ）是雷达接收天线在干扰方向的主瓣增益；λ表示信号波长；Ｇｔ（是干扰机与雷达的距离；Ｒｔ是雷达与目标的距离；Ｒｊ；ξｊ是雷达信号与干扰信号之间的极化适配损失系数。

如果雷达设备中目标检测器输入端的信噪比为Ｓ概率为ｎ，目标检测器的平均信噪比为ＰＳｎ，设雷达的虚警ｆａ同公式表示，则不同类型信号的检测概率可以用不１）已知信号的检测概率：

Ｐｄ＝Ｆ[Ｆ－１

：

式中，Ｆ（∗）和Ｆ－１（∗）（定义如下Ｐｆａ）－Ｓ：

ｎ]

Ｆ(ｘ)＝

若Ｐｆａ＝

∫¥－ｔ２／２

∫¥２π１ｅ

－ｔ２／２

ｘ

ｄｔ

限，Ｆ（ｘ）的值难以计算２π１

－１

ｘ

ｅｄｔ，则ｘ＝Ｆ（Ｐｆａ），ｘ为检测门

，但是可以通过查表的方式获取。

２）未知信号的检测概率Ｐ：未知信号的检测是与虚警概率和雷达目标检测ｄ≈Ｆ[２ｌｎ（１／Ｐｆａ）－

Ｓｎ]

　器的输入信噪比有关　，当信噪比越大、虚警概率越小

时，目标检测概率越大３）斯韦林型信号的检测概率，反之，则越小。

ìï－

：

Ｐï（Ｐｆａ－

）１／（１＋Ｓｎ）

ｄ＝íïæＳｎöæïçî

ç

２＋Ｓ－÷ｌｎＰｆａöæｌｎＰｆａè÷çｎøè１＋Ｓ－２

－ｎ１＋Ｓ－ｎ／２÷øｅｘｐçöè１＋Ｓ－ｎ／２÷ø不同目标的起伏类型是不同的，通常采用斯韦林型（３）

目标起伏，斯韦林Ⅰ、Ⅱ起伏型目标如舰船、渔船、汽车等一般采用公式（３）中第１行公式来近似得到该类型的目标检测概率，斯韦林Ⅲ、Ⅳ起伏型目标如战斗机等则通过公式（３）第２行公式近似得到相应目标的检测概率。

通过对比在实施干扰前、后同一目标检测概率的变化可以得到干扰效果，并根据干扰结果进一步优化干扰分配策略。由于雷达在干扰前、后的信噪比与信干比是衡量目标回波信号质量的重要指标，因此利用式（１）和式（２）计算的结果带入式（３），可以得到实施干扰前、Ｐ、后雷达对目标信号的检测概率分别表示为Ｐｄｄｊ。每部雷达在搜索阶段的目标

函数可以定义为

ＥＰｄ＝

ｄ－Ｐｄｊ

　合后的检测概率可以表示为

　根据文献［１１］，协同组网雷达信号通过主站融Ｐｄ

（４）

ＰＤ＝

ｉ

ｉ

（５）

式中：Ｄ表示组网雷达数据融合中心的判决向量∑Ｄ[Ｒ（Ｄ）∏Ｓ０

１－Ｐｄ∏Ｓ１

Ｐｄ]。融合中心决定是否融合的规则为Ｒ，Ｒ（Ｄ）的可以表示为

ìï１，　Ｎ

Ｒ（Ｄ）＝ï

í

∑ｎ＝１

ｄｎ≥ＫïＮ（６）

ïî

０，　ｎ＝１

＜Ｋ

　∑ｄｎ带入式　式中（４）：Ｎ即可得到雷达网的损失表示组网雷达的雷达数量。

。将式（５）根据式（４），可知Ｅｄ与雷达对目标的检测概率成反比，其数值越大，表示干扰效果越好，反之，则干１．２．２　扰效果越差在组网雷达融合中心判断目标存在定位阶段目标函数。

⁃ＧＤＯＰ

，即发现目

标之后，组网雷达即由搜索目标阶段转入定位目标阶段。在精确定位中，雷达网会根据测得与目标的斜距、方位角、俯仰角３个方面来确定目标空间位置，因此，选取ＧＤＯＰ建立目标函数。

根据文献［１２］，假设目标的真实斜距、方位角、

俯仰角分别为Ｒ距、方位角、俯仰角为１、α１和βＲ１，雷达真实测量的目标斜ｙ、ｚ轴的定位误差

２、α２和β２，在直角坐标系ｘ、

ìïｘ－

＝ｘï２－ｘ１

íïｙ－

＝ｙ２－ｙ１ï－

ｚ＝ｚ　　式中

î

２－ｚ１ìïï

ｘ１＝Ｒ１ｃｏｓα１ｓｉｎβ１íïｙïî

ｚ１＝Ｒ１ｓｉｎβα１ｃｏｓβ１１＝Ｒ１ｓｉｎ１ìïï

ｘ２＝Ｒ２ｃｏｓα２ｓｉｎβ２íïｙ２＝Ｒï２ｓｉｎα２ｃｏｓβ２　î

ｚ２＝Ｒ２ｓｉｎβ２角分别为　设某时刻雷达测量的目标斜距Ｒ、α和β，则其在ｘ、ｙ、ｚ轴的定位误差均、方位角和俯仰方差为

第４期　　　　　　　　　　　　汲清波，等：基于自适应目标函数的协同干扰策略分配　　　　　　　

２２２２２２２２２２２２ìïσｘ＝σＲｃｏｓαｃｏｓβ＋ｒσαｓｉｎαｃｏｓβ＋Ｒσβｃｏｓαｓｉｎβï２２２２２２２２２２２２íσｙ＝σＲｓｉｎαｃｏｓβ＋ｒσαｃｏｓαｃｏｓβ＋Ｒσβｓｉｎαｓｉｎβï·７３·

ïσ２

ｚ＝σ２Ｒ雷达在该时刻的ｓｉｎ２β＋Ｒ２σ２

β２　　î

ｃｏｓβ　　对ＧＤＯＰＧＤＯＰ做处理得

＝ＧＤＯＰ（σ２ｘ＋表达式为σ２ｙ＋σ２ｚ）Ｇｊ

ｄ＝

ＧＤＯＰ－ＧＤＯＰ式中：ＧＤＯＰ为无干扰时的定位精度ＧＤＯＰ

，ＧＤＯＰｊ表示在干扰后的定位精度，采用Ｇｄ作为目标函数，该值越大，说明干扰效果越好，反之，则说明干扰效果越差。

１．２．３　在不同阶段采用不同的目标函数评估干扰效果

自适应目标函数的建立

以适应雷达工作任务的转变，从而更好地对干扰方的干扰资源进行分配。将距离均匀划分为每一小段，记为采样点数，本文采用自适应目标函数的粒子群算法对干扰资源进行分配。算法步骤为：成一个种群１）初始化干扰策略Ｓ＝{Ｘ，随机产生ｍ个干扰策略组

或１，Ｘ１２随机数，…，Ｘｉ；

，…，Ｘｍ}，其中每个干扰策略都是的０２）计算该采样点目标与各雷达间的距离，运用粒子群计算

Ｒ；雷达搜索阶段目标函数并计算每一采样点目标函数３）根据距离Ｒ选择目标函数最优值与最优分配策略；

２　４）协同干扰策略分配算法

输出各采样点目标函数最优值。

协同干扰策略分配算法是策略分配的核心。本文采用改进后的粒子群算法分配干扰策略，该算法的好处是显而易见的，比如能够解决传统粒子群算法中容易陷入到局部最优解而无法成功跳出的缺点，又能够解决采用遗传算法不能快速手链的缺陷。在粒子群算法的初期，由于算法本身的设定问题，容易早熟，为解决该问题，本文将模拟退火机制应用于粒子群算法中，模拟退火算法中的核心是通过温度的突变性导致算法跳出局部最优解，而本文将两种算法融合后既能够解决局部最优解问题又能实现算法的快速收敛。基于模拟退火机制的粒子群算法并不是只在当前最优结果附近搜索，而是有一定概率跳出当前最优解寻找全局最优值，因此解决了早熟的问题，退火算法的温度确定公式为

ＴＦ（ｐｉ）＝∑Ｎ

ｅ－（ｆ（ｐｉ）－ｆ（ｐｇ））／ｔ

ｅ－（ｆ（ｐｉ）－ｆ（ｐｇ））／ｔ

（７）

　　基于模拟退火的粒子群算法的步骤为ｉ＝１

速度，１）该值是随机生成的初始化粒子群中各粒子的初始位置和飞行；

中各粒子的位置和飞行速度以及对应的适应度存储２）计算种群中各粒子的适应度函数；将每一代

在该粒子对应的ｐ的粒子存储在全局最优解ｂｅｓｔ，把所有粒子中的适应度最好３）初始化模拟退火机制中的退火温度ｇｂｅｓｔ中；

；及其温度４）根据公式；

（７）确定每一个各粒子适应度函数５）根据随机概率选出ｐｉ作为在当前粒子中确

定的局部最优用于替代当前的全局最优值ｐｇ式（８）、（９），根据＋１）更新位置和速度＝φ｛ｖ（ｔ）＋；

ｖｉ，ｊ（ｔｉ，ｊｃ１ｒ１（ｐｉ，ｊ－ｘｉ，ｊ（ｔ））＋

ｃ２ｒ２（ｐｇ，ｊ（ｔ）－ｘｉ，ｊ（ｔ））｝（８）ｘｉ，ｊ（ｔ＋１）＝ｘｉ，ｊ（ｔ）＋ｖｉ，ｊ（ｔ＋１）

式中φ＝

２－Ｃ－

２Ｃ２－４Ｃ，Ｃ＝ｃ（９）

１＋ｃ２。

骤５）；

６）计算粒子群中的各粒子适应度函数，重复步７）足，则停止重复步骤８）退温判断步骤；

６）的停止条件是否满足６）输出当前全局最优解，如果满，如果不满足９），如何选取退火机制中的初始温度以及退温

继续步骤２）～６）。

公式的选取都对最终输出结果有影响（１０）和（１１）进行初温和退温操作ｔ。

，一般根据式０＝ｆ（ｐｇｔ）λｔ／ｌｎ５（１０）ｋ＋１＝ｋ

３　仿真结果分析

（１１）

假设雷达网由８部雷达组成，其布站采用方形布站，保护以原点为中心边长为４０ｋｍ的正方形区域，多机编队由４部干扰机以及目标组成，目标编队从空间坐标（２００，２００，１０）开始突防，其突防态势图以及多机编队的航迹如图１所示。

·７４·应　　　用　　　科　　　技　　　　　　　　　　　　　　　　第４６卷

　　根据文中１．２．１节所述，在干扰时，由于距离的离上的威胁程度是不同的。如雷达１，在距离较远时，载频低，探测距离远，距离分辨度以及速度分辨

不同，雷达任务不同，具有同一指标的雷达在不同距

度不高，此时，低载频的权重越高；在距离较近时，由于其带宽、重频较低，因此威胁程度相较于雷达８低。搜索及定位阶段指标权重如表２所示。

图１　突防态势图

表２　不同阶段雷达指标的权重载频／ＭＨｚ

搜索阶段定位阶段

带宽／ＭＨｚ

１０５５４０１００２０６０４５８

１３１６４０２５１０２５１８２４脉宽／μｓ

０．３４０．１５

重频／Ｈｚ０．１５０．３７

脉宽／μｓ０．３００．１２

带宽／ＭＨｚ０．２１０．３６

干扰策略分配是处于侦察分选之后，根据侦察得到的雷达息如表１所示。

表１　雷达信息

编号雷达１雷达２雷达３雷达４雷达５雷达６雷达７雷达８

载频／ＭＨｚ２０００９６００１１０００４８００６１００１２０００１００００３０００

重频／Ｈｚ５００１２５０３８５２５００６２５５０００１６００３５７

　　根据文献［１０］方法确定的雷达威胁等级如表３

所示。由表３可以得知，同一部雷达在相同的工作参数下但在不同的工作阶段，其威胁程度是不同的，但是根据雷达工作原理可知，在搜索阶段，雷达低载频、宽脉宽使得其信号损失较小，探测距离更远，而具有更高重频、宽带宽的雷达信号尽管工作距离有限，但是其探测精度更高，因此在定位跟踪阶段威胁更高。

雷达４０．５３０．６１

雷达５０．３９０．３５

雷达６０．２１０．７１

雷达７０．２２０．５５

雷达８０．４５０．２１

表３　雷达威胁等级

所在阶段搜索阶段定位阶段

雷达１０．４２０．２５

雷达２０．２４０．５７

雷达３０．２８０．３４

　　由前文１．２节所述，通过对单部干扰机对单部雷达的干扰前后的发现概率如图２所示。

标，该转入到定位跟踪阶段，如果还以检测概率为指标评估干扰效果，只会是对干扰资源巨大的浪费。

雷达干扰策略分配根据步骤分为两步，首先是干扰对象的分配：需要确定哪部干扰机干扰哪部雷达；然后再分配相应的干扰样式。以４部干扰机需要干扰４部雷达为例说明，在只考虑干扰对象的一对一分配情况下，一共需要考虑的情况为１８８种；而如果要将一对多考虑进去进行策略分配，分配方法采用传统的动态规划算法，则需要考虑的分配情况可以达到１２９６种，计算量的增长带来的分配结果不具有实时性，而随着干扰机可能释放干扰信号的类型的增多，如３种则需要分配的策略可高达进行计算，单次计算最有分配结果的时间将在粒子群算法后寻找全局最优解只需要４～５ｓ，寻优时间大为缩减。

根据文中１．２节目标函数的论述，融合中心判断规则为：当有２部以上雷达（即式（６）中Ｋ＝２）认为目标存在，即开始计算目标检测概率。突防过程中，采用４架干扰机干扰８部雷达，选取采样点８５时，计算此时目标函数的最优值，如图３所示。

图２　干扰前、后雷达检测概率

由图２可以看出，在干扰前、后的雷达检测概率变化较大；图中采样点的数值越大表示离目标越近。在干扰前，组网雷达的检测概率在采样点为７０，即目标的坐标为（６０，６０，４）时能完全发现目标，而干扰后则为９０，即目标坐标为（２０，２０，２）。由此可看出，干扰效果明显，验证了检测概率在搜索阶段作为干扰效果评估准则的可行性。同时，图２也说明在采样点８５之后，无论如何对雷达进行干扰，对雷达的探测性能影响微乎其微，即雷达已经完全发现目

１０４９７６种，假设仍然采用动态规划算法对资源分配２０ｍｉｎ以上，而这段时间战情早已发生变化，而采用

第４期　　　　　　　　　　　　汲清波，等：基于自适应目标函数的协同干扰策略分配　　　　　　　·７５·

同，雷达网的发现概率在采样点为８５时最高可以降低为无干扰时的０．５倍左右，说明在不同距离上如果采用同一分配策略得到的干扰效果是不同的，也不是最优的，因此，要根据距离按照干扰评估准则做符合该采样点的最优分配策略才能使全程的干扰效果最优。

在突防全程只采用检测概率作为目标函数评估干扰效果，得到单一目标函数全程干扰效果如图５所示。图５在全程只采用检测作为目标函数，在采

图３　采样点８５基本粒子群算法寻优效果

样点９０之后，干扰效果相较于采样点９０之前的干扰效果急剧下降并且都基本保持不变，是无干扰时的０．１倍左右，结合图２可以发现，此时目标基本上已经被雷达发现，如果还以检测概率作为目标函数，无论如何改变干扰策略都基本不会影响雷达网的性能。同时，这也说明，在此采样点之后，干扰机编队应该改变干扰评估效果以及干扰方式，从而取得适应工作阶段的变化。以检测概率作为协同干扰评估指标干扰机的利用率如表４所示。

从图４中可以看出用基本粒子群算法寻优时，在１００代内算法都无法收敛。利用基于模拟退火的粒子群算法寻优如图４所示。

图４　基本模拟退火的粒子群算法寻优效果

种群的大小影响着单采样点的最优解，种群数为３２时的收敛速度相比于种群数为１６和８时的收敛速度要更快。同时，可以看到，由于分配策略的不　　

图５　单一目标函数全程干扰效果

表４　单一目标函数干扰机使用效率

种群数

雷达１０．８１

种群数１６

０．９００．８５０．７９０．８８

种群数３２

０．９００．８９０．９００．９１

种群数６４

０．８９０．９１０．９１

雷达２０．１４０．２１０．１９０．１５０．０９０．０９０．１１０．１３０．０６０．０６０．０６０．０５

雷达３０．４６０．４１０．４５０．６１０．５６０．３９０．５２０．６２０．４７０．８６０．６００．２４

雷达４０．８８０．８８０．８９０．９００．８９０．９３０．９５０．９１０．９３０．９３０．９５０．９５

雷达５０．６１０．５３０．５３０．３７０．５３０．６１０．５３０．３８０．６３０．２３０．４８０．８６

雷达６０．０４０．０４０．０７０．０６０．０３０．０４０．０２０．０４０．０２０．０１０．０１０

雷达７０．１９０．１５０．１３０．２３０．１２０．１５０．０８０．０９０．０５０．０６０．０４０．０５

雷达８０．８７０．８８０．８９０．８９０．９００．８９０．９２０．９５０．９１０．９５０．９５０．９３

·７６·应　　　用　　　科　　　技　　　　　　　　　　　　　　　　第４６卷

　　从表４可以看出，如果突防过程全程以检测概２、７等的分配概率极低，这使得即使在雷达发现了突防编队之后，干扰机仍然会给搜索阶段威胁程度大的雷达分配干扰资源，这会造成对干扰资源的巨大浪费。因此，采用本文中自适应目标函数对协同干扰资源进行分配，其干扰效果如图６所示，其干扰效率如表５所示。

从图５、６的对比中可以看出，在采样点为９０之

场实际，更符合实际突防过程中所遭遇的情况。

率作为干扰评估指标对干扰策略进行分配，雷达６、

前，以发现概率为目标函数，干扰效果是随着距离的９０以后采用定位精度作为干扰效果评估可以更好的利用干扰资源，并且雷达网的定位精度是无干扰更有利于充分利用干扰资源，同时也更符合雷达网工作任务的转变。相比较于图５只是单一的采用发现概率作为干扰效果评估准则，图６在采样点为９０之后的分配效果更佳，分配策略也更加接近于战　　

时的几倍，说明在该阶段，以定位精度作为目标函数不断接近呈现先增加后逐渐减少的过程，采样点为

图６　自适应目标函数全程干扰效果

从表４、５中对比来看，干扰机的资源使用效率更为科学合理。在实际突防过程中，突防编队会先后遭遇搜索、定位两个阶段。在搜索阶段，更多的是将干扰资源向探测距离更远的雷达４、１、８倾斜；而在定位阶段，应该将资源更多地分配给重频高、脉宽宽的雷达６、２、７。结合表３可以看出，以自适应目标函数干扰资源策略分配更为科学更符合实际。

表５　自适应目标函数干扰机使用效率

种群数

雷达１０．７５

种群数１６

０．７１０．７２０．７３０．７６

种群数３２

０．７８０．７９０．７９０．７７

种群数６４

０．７８０．７９０．７８

雷达２０．２８０．３３０．３５０．２３０．２４０．２６０．２１０．１９０．２００．１９０．１８０．１９

雷达３０．５００．５６０．４７０．４４０．４７０．４６０．５００．５８０．５３０．３２０．４２０．７７

雷达４０．７１０．７８０．７８０．７７０．７７０．８１０．８００．７７０．８００．８１０．８２０．８２

雷达５０．５７０．４００．４８０．５７０．５５０．６００．５３０．５００．５５０．７７０．６３０．３１

雷达６０．１７０．１７０．１６０．１５０．１７０．１３０．１４０．１６０．１４０．１４０．１５０１５

雷达７０．３００．２８０．２８０．２８０．２６０．１８０．２５０．２２０．２００．２００．２００．１７

雷达８０．７２０．７７０．７６０．７７０．７８０．７８０．７８０．８００．７９０．８２０．８１０．７

４　结论

大功率噪声压制，而在定位阶段，可以采取距离欺骗、角度欺骗等方式来分析不同干扰方式对干扰效果的影响，这也是下一步有待完善的地方。

从理论分析以及仿真中可以看到，在飞机突防

过程中，干扰效果并非是一成不变的，由于雷达所处的工作阶段不同，采用不同的目标函数以适应雷达工作的不同阶段，该方法是可行的。文中分阶段评估雷达的威胁等级，并且采用行之有效的干扰效果评估可以最大化利用己方干扰资源。同时，论文仍然存在一定问题，没有分阶段７虑采用不同的干扰方式对干扰效果的影响，比如在搜索阶段，可以采用

参考文献：

［１］邢美丽，肖军鹏．复杂电磁环境下协同干扰技术研究［Ｊ］．

遥测遥控，２０１０，（４）：３０－３５．

［２］陈雪岩，李宝龙．协同干扰策略分配及评估的研究与实

现［Ｊ］．应用科技，２０１８（４４）：１－５．

（下转第８１页）

第４期　　　　　　　　　　　　　　　韩月，等：先验信息约束ＮＭＦ的高光谱解混　　　　　　　［６］ＭＩＡＯＬ，ＱＩＨ．Ｅｎｄｍｅｍｂｅｒｅｘｔｒａｃｔｉｏｎｆｒｏｍｈｉｇｈｌｙｍｉｘｅｄ

ｄａｔａｕｓｉｎｇｍｉｎｉｍｕｍｖｏｌｕｍｅｃｏｎｓｔｒａｉｎｅｄｎｏｎ⁃ｎｅｇａｔｉｖｅｍａｔｒｉｘｆａｃｔｏｒｉｚａｔｉｏｎ［Ｊ］．ＩＥＥＥｔｒａｎｓａｃｔｉｏｎｓｏｎｇｅｏｓｃｉｅｎｃｅ＆ｒｅｍｏｔｅ［７］ＱＩＡＮＹ，ＪＩＡＳ，ＺＨＯＵＪ．ＨｙｐｅｒｓｐｅｃｔｒａｌｕｎｍｉｘｉｎｇｖｉａＬ１／２

ｓｐａｒｓｉｔｙ⁃ｃｏｎｓｔｒａｉｎｅｄｎｏｎｎｅｇａｔｉｖｅｍａｔｒｉｘｆａｃｔｏｒｉｚａｔｉｏｎ［Ｊ］．４９（１１）：４２８２－４２９７．（４）：１２４７－１２５４．

ＩＥＥＥｔｒａｎｓａｃｔｉｏｎｓｏｎｇｅｏｓｃｉｅｎｃｅ＆ｒｅｍｏｔｅｓｅｎｓｉｎｇ，２０１４，ｓｅｎｓｉｎｇ，２００７，４５（３）：７６５－７７７．

·８１·

［１０］ＰＡＡＴＥＲＯＰ．Ｌｅａｓｔｓｑｕａｒｅｓｆｏｒｍｕｌａｔｉｏｎｏｆｒｏｂｕｓｔｎｏｎｎｅｇａ⁃

ｔｉｖｅｆａｃｔｏｒａｎａｌｙｓｉｓ［Ｊ］．Ｃｈｅｍｏｍｅｔｒｉｃｓ＆ｉｎｔｅｌｌｉｇｅｎｔｌａｂｏｒａ⁃ｔｏｒｙｓｙｓｔｅｍｓ，１９９７，３７（１）：２３－３５．

［１１］方凌江，粘永健，雷树涛．基于顶点成分分析的高光谱［１２］ＨＯＹＥＲＰ．Ｎｏｎ⁃ｎｅｇａｔｉｖｅｍａｔｒｉｘｆａｃｔｏｒｉｚａｔｉｏｎｗｉｔｈｓｐａｒｓｅ⁃［１３］ＩＯＲＤＡＣＨＥＭ．Ｄ，ＢＩＯＵＣＡＳＤＩＡＳＪ．Ｍ，ＰＬＡＺＡＡ．

２０１４－２０３９．

２００４，５：１４５７－１４６９．１５４－１５７，１８１．

图像端元提取算法［Ｊ］．舰船电子工程，２０１４，３４（８）：

［８］魏一苇，黄世奇，王艺婷．基于体积和稀疏约束的高光

ｎｅｓｓｃｏｎｓｔｒａｉｎｔｓ［Ｊ］．Ｊｏｕｒｎａｌｏｆｍａｃｈｉｎｅｌｅａｒｎｉｎｇｒｅｓｅａｒｃｈ，Ｓｐａｒｓｅｕｎｍｉｘｉｎｇｏｆｈｙｐｅｒｓｐｅｃｔｒａｌｄａｔａ［Ｊ］．ＩＥＥＥｔｒａｎｓａｃ⁃

谱混合像元分解算法［Ｊ］．红外与激光工程，２０１４，４３

［９］ＴＡＮＧＡ，ＳＨＩＺ，ＡＮＺ．Ｎｏｎｎｅｇａｔｉｖｅｍａｔｒｉｘｆａｃｔｏｒｉｚａｔｉｏｎｆｏｒ

ｈｙｐｅｒｓｐｅｃｔｒａｌｕｎｍｉｘｉｎｇｕｓｉｎｇｐｒｉｏｒｋｎｏｗｌｅｄｇｅｏｆｓｐｅｃｔｒａｌｓｉｇｎａｔｕｒｅｓ［Ｊ］．Ｏｐｔｉｃａｌｅｎｇｉｎｅｒｉｎｇ，２０１２，５１（８）：７００１．

ｔｉｏｎｓｏｎｇｅｏｓｃｉｅｎｃｅａｎｄｒｅｍｏｔｅｓｅｎｓｉｎｇ，２０１１，４９（６）：

本文引用格式：韩月，康维新，李慧．先验信息约束ＮＭＦ的高光谱解混［Ｊ］．应用科技，２０１９，４６（４）：７７－８１．

ＨＡＮＹｕｅ，ＫＡＮＧＷｅｉｘｉｎ，ＬＩＨｕｉ．Ｈｙｐｅｒｓｐｅｃｔｒａｌｕｎｍｉｘｉｎｇｏｆｎｏｎｎｅｇａｔｉｖｅｍａｔｒｉｘｆａｃｔｏｒｉｚａｔｉｏｎｗｉｔｈｐｒｉｏｒｉｎｆｏｒｍａｔｉｏｎ［Ｊ］．Ａｐｐｌｉｅｄｓｃｉ⁃ｅｎｃｅａｎｄｔｅｃｈｎｏｌｏｇｙ，ＤＯＩ：１０．１１９９１／ｙｙｋｊ．２０１９，４６（４）：７７－８１．

（上接第７６页）

８７９．８０．

［３］黄玉川，饶妮妮，刘永红，等．博弈论应用于干扰效果动态

评估的研究［Ｊ］．电子科技大学学报，２００７，３６（５）：８７６－［４］房明星，王杰贵，朱宇．基于Ｈｏｐｆｉｅｌｄ神经网络算法的雷

达干扰资源分配［Ｊ］．火力与指挥控制，２０１４，３９（２）：７６－［５］李东生，高杨，雍爱霞．基于改进离散布谷鸟算法的干扰［６］张运豪．基于免疫算法的多策略雷达干扰资源分配方法［７］唐建强，李昊．基于蚁群算法的对抗资源多目标分配［Ｊ］．［８］ＢＡＣＨＡＭＡＮＮＤＪ，ＥＶＡＮＳＲＪ，ＭＯＲＡＮＢ．Ｇａｍｅｔｈｅｏ⁃

电子信息对抗技术，２０１８，３３（２）：４８－５２．［Ｄ］．西安：西安电子科技大学，２０１４．８９９－９０５．

ｒｅｔｉｃａｎａｌｙｓｉｓｏｆａｄａｐｔｉｖｅｒａｄａｒｊａｍｍｉｎｇ［Ｊ］．ＩＥＥＥｔｒａｎｓａｃ⁃１０８１－１１００．

ｔｉｏｎｓｏｎａｅｒｏｓｐａｃｅ＆ｅｌｅｃｔｒｏｎｉｃｓｙｓｔｅｍｓ，２０１１，４７（２）：ｊａｍｍｉｎｇａｌｌｏｃａｔｉｏｎ［Ｊ］．ＩＥＴｒａｄａｒ，ｓｏｎａｒ＆ｎａｖｉｇａｔｉｏｎ，

［９］ＮＩＣＨＬＯＲＳＲ，ＷＡＲＲＥＮＰ．Ｐｌｅｓｓｉｓ．Ｔｈｒｅａｔｅｖａｌｕａｔｉｏｎａｎｄ［１０］呙鹏程，王星，周东青，等．基于熵权的灰色关联辐射源［１１］陈彬，童创明，李西敏．基于检测概率的雷达组网反隐身

建模与仿真［Ｊ］．系统仿真学报，２０１７，２９（１２）：３１００－３１０６．

［１２］崔炳福．雷达对抗干扰有效性评估［Ｍ］．北京：电子工业

出版社，２０１６：７０－７１．１７（６）：５８－６３．２０１７，１１（３）：４５９－４６５．

资源算法分配研究［Ｊ］．电子与信息学报，２０１６，３８（４）：威胁评估［Ｊ］．空军工程大学学报（自然科学版），２０１６，

本文引用格式：

汲清波，杨帅．基于自适应目标函数的协同干扰策略分配［Ｊ］．应用科技，２０１９，４６（４）：７０－７６，８１．ｔｅｃｈｎｏｌｏｇｙ，２０１９，４６（４）：７０－７６，８１．

ＪＩＱｉｎｇｂｏ，ＹＡＮＧＳｈｕａｉ．Ｃｏｏｐｅｒａｔｉｖｅｉｎｔｅｒｆｅｒｅｎｃｅｓｔｒａｔｅｇｙａｌｌｏｃａｔｉｏｎｂａｓｅｄｏｎａｄａｐｔｉｖｅｏｂｊｅｃｔｉｖｅｆｕｎｃｔｉｏｎ［Ｊ］．Ａｐｐｌｉｅｄｓｃｉｅｎｃｅａｎｄ