基于信息维数的欺骗干扰存在性检测
祝宏1,2,3, 张海1,3, 唐高弟1, 李中云1, 刘以农2
1.中国工程物理研究院 电子工程研究所,四川 绵阳 621999
2.清华大学 工程物理系,北京 100084
3.中国工程物理研究院 中物院复杂电磁环境科学与技术重点实验室,四川 绵阳 621999
通讯作者:张海(1968-),男,研究员.研究方向:雷达系统及信号处理.E-mail:e2zhanghai@aliyun.com

作者简介:祝宏(1987-),男,博士研究生.研究方向:雷达信号处理,雷达抗干扰技术.E-mail:zhu-h06@163.com

摘要

利用欺骗干扰信号与目标回波包络上的不同,基于分形理论中的信息维数提出了一种欺骗干扰的存在性检测方法.当信号本身的信息维数与信号包络的信息维数之差超过一定门限时,可判断存在欺骗干扰,同时可获知欺骗干扰的起始时刻.仿真结果表明:本文算法在干噪比不低于6 dB时,检测概率在90%以上.

关键词: 雷达工程; 抗干扰; 欺骗干扰; 分形维数; 信息维数
中图分类号:TN974 文献标志码:A 文章编号:1671-5497(2016)02-0616-05
Detection of presence of deception jamming based on information dimension
ZHU Hong1,2,3, ZHANG Hai1,3, TANG Gao-di1, LI Zhong-yun1, LIU Yi-nong2
1.Institute of Electronic Engineering,CAEP,Mianyang 621999,China
2.Department of Engineering Physics,Tsinghua University,Beijing 100084,China
3.Key Lab of Science and Technology in Complex Electromagnetic Environment,CAEP,Mianyang 621999,China
Abstract

A detection algorithm of the presence of deception jamming is proposed based on information dimension with the fact that the signal envelopes of deception jamming and target return are different. The algorithm can detect the existence or not of a deception jamming and find its initial time after calculating the difference between the signal's information dimension and the signal envelope's dimension. Simulation results show that this approach can detect deception jamming effectively when JNR is no less than 6 dB, and the detection probability is higher than 90%.

Keyword: radar engineering; anti-jamming; deception jamming; fractal dimension; information dimension
0 引 言

抗欺骗式干扰技术一直是雷达研究者的关注重点.目前, 关于欺骗式干扰的研究主要集中在原理分析, 性能仿真和效果评估方面, 关于抗干扰的研究则相对较少[1, 2].文献[3]分析了数字射频存储器(DRFM)产生的距离波门拖引(RGPO)干扰信号由于DRFM相位量化误差在频域上产生的细微变化, 但这个细微变化在DRFM阶数较高时不明显, 且不适用于目标回波和干扰同时存在的情况.文献[4]分析了目标回波与欺骗干扰信号脉内幅度起伏的差异; 文献[5]研究了两者在幅度起伏, 高阶统计量, 双谱特性上的差异, 但均分析简单, 没有研究目标回波与欺骗干扰部分重叠的情况.

自分形理论被提出以来, 分形理论已成功地应用于图像处理, 语音识别, 雷达信号处理等领域.文献[6, 7, 8]分别从不同角度研究了分形理论在雷达目标检测中的应用, 但均不适用于存在欺骗干扰的情况.对此, 本文利用欺骗式干扰与目标回波在信息维数上的差异, 提出了一种新的欺骗干扰存在性检测方法.本文采用实测的大地回波数据进行仿真, 通过分析大地回波和欺骗干扰的信息维数特征, 得出判定欺骗干扰存在与否的依据, 为后续的抗欺骗干扰算法设计打下了良好基础.

1 信息维数

分形维数是分形理论的中心概念, 它有多种定义和计算方法, 常用的是Hausdorff维数, 盒维数, 信息维数, 相似维数等[9, 10].雷达信号作为一种时间序列, 可以用分形维数对其进行刻画, 比如盒维数用以描述信号的复杂度, 信息维数用以描述信号在二维平面上的分布信息.本文采用后者对雷达信号进行研究.

(X, d)为一个度量空间, HX的非空紧集族, ε为一个非负实数, 令 B(x, ε)表示一个中心在 x, 半径为 ε的闭球.设A是X中的一个非空紧集, {Aj}(j=1, 2, , et al., K)是集合A的一个有限 ε-格形覆盖, pj表示 A的元素在集合 Aj中的概率.令信息熵为I(ε )=- j=1Kpjlgpj, 如果 I(ε)满足: I(ε)~-lgεDI(f), 则称 DI为集合 A的信息维数, 定义为[11]:

DI=-limε0I(ε)lgε(1)

实际计算中通常通过粗视化变换得到信息维数:

DI=[I(λ)-I2λ)]/lg2

基于信息维数的定义, 本文设计了以下算法步骤进行计算:

Step1 数据预处理.首先, 对信号数据线性插值, 使得数据长度为2的整数次幂, 设插值后的数据长度为Length; 然后, 将每个数据点都减去这段信号的最小值, 使得这段信号数据的最小值为0, 设处理后的数据最大值为Data_max; 最后, 将每个数据同乘以Length/Data_max, 使得这段信号数据的最大值等于数据长度.以盒子覆盖信号时, 计算落入每个盒子的信号曲线长度比较麻烦, 所以本文采用计算落入每个盒子的数据点数来替代, 故而需要对信号数据进行插值, 以提高计算精度.使数据长度为2的整数次幂的目的是为了便于设定盒子尺寸, 以及便于运算.对信号的最小值和最大值预处理的目的是为了便于摆放盒子, 且能够满足分形维数与信号幅度无关的特点.

Step2 计算信息熵.依照信息维数的定义, 计算不同盒子尺度下的信息熵.盒子尺度依次选取2的整数次幂, 且不大于数据长度.

Step3 选取无标度区间.作 I(ε)~-lgε曲线, 选取线性度较好的部分作为无标度区间.

Step4 计算信息维数.对无标度区间内的曲线进行线性拟合, 得到的斜率即为信息维数.

2 检测算法设计
2.1 检测原理

欺骗式干扰主要是通过DRFM对接收到的雷达发射脉冲进行采集, 存储, 调制放大, 转发来产生的, 从而逼真地模拟目标的相关特性, 使得二者在时域, 频域上的差异很小.随着现代数字处理软硬件技术的发展, 一些新型的干扰机在转发信号的同时能够结合脉内调制, RCS调制等复合手段, 实现欺骗信号与目标回波的高度相似.然而这些新型干扰机目前大多还处于研制和实验阶段, 并未大量装备[2].因此, 本文从实际需要出发, 只针对传统的转发式欺骗干扰, 对于新型欺骗干扰则没有过多考虑.由于雷达的发射脉冲往往是恒包络的, 所以转发式欺骗干扰的包络也是恒定的, 而目标回波的幅度受到目标散射面的调制, 其包络是不规则的, 如图1所示.

图1 大地回波和欺骗干扰的时域波形Fig.1 Time-domain signal of terrain return and deception jamming

图1可知:大地回波的包络与欺骗干扰的包络具有明显差异, 所以它们的包络的分形特征也会有明显差异, 而一般的目标分形检测并没有对信号包络的分形维数进行分析.对于存在欺骗干扰的数据段, 由于干扰信号的包络是个直流电平, 由分形理论的基本性质和信息维数的计算流程均不难看出该数据段的信息维数与噪声比较接近.在信噪比较大时, 存在目标回波的数据段信息维数将小于噪声的信息维数.经过研究发现, 如果仅仅利用对信号包络求信息维数并取最小值来检测目标回波, 其性能很差, 一是因为对信噪比的要求很高; 二是因为当欺骗干扰较强时, 在目标回波和欺骗干扰之间会存在极小值, 影响检测的准确性.为此, 本文先逐段分别计算信号包络和信号本身的信息维数, 然后求两个信息维数的差.前面的分析已经指出, 欺骗干扰的两个信息维数差值比较大, 而噪声, 目标回波的则相对比较小, 故可以据此判断是否存在欺骗干扰, 并且检测出干扰的存在位置.

2.2 算法设计

依次以每个信号数据点为起点, 划分长度为发射脉冲宽度的数据段, 计算每个数据段 Dataj的信息维数 Aj, 信号包络的信息维数 Bj, 两者之差为 Lj根据前文的分析, 干扰检测算法如下:

Lj=Bj-Aj> ε, 受到干扰ε, 未受干扰(2)

式中: ε为干扰检测门限. ε应大于目标回波单独存在时 Lj能取到的最大值, 避免目标回波被误判为欺骗干扰.

(1)当 Ljε时, 说明不存在欺骗干扰, 或者欺骗干扰不够强.此时, 可按常规的匹配滤波法进行检测, 也可以按文献[6, 7, 8]研究的分形检测法进行后续检测.

(2)当 Lj> ε时, 说明雷达受到干扰. Lj取到最大值的时刻即为欺骗式干扰的起始时刻.知道了欺骗干扰的存在位置, 可以后续采用很多抗干扰算法来抑制干扰, 比如简单的置零处理.

3 仿真及分析

本文采用雷达高度表实测的大地回波数据作为目标回波.选取信噪比很好的某个大地回波, 将其下变频到零中频(消除载波影响), 作为干净的目标回波, 其发射脉冲是线性调频信号, 带宽为10 MHz, 脉宽为5 μ s.由于测量时挂飞的直升机垂直速度很小, 且多普勒频率对分形维数影响很小, 故本文不考虑目标回波的多普勒频率.欺骗式干扰的参数与发射脉冲参数相同, 多普勒频率也假设为零.构造一个长度为20 μ s的雷达接收数据, 其中目标回波起始时刻是1 μ s, 欺骗干扰的起始时刻是4 μ s(与目标回波部分重叠)或者8 μ s(与目标回波分离), 噪声为高斯白噪声.

分形理论所研究的系统是建立在自相似基础之上的, 对于实际的系统, 这只可能在一定的区间内成立, 即无标度区间, 超出这个区间的自相似性就不存在了, 分形也就失去了意义.图2给出了大地回波的信息熵与尺度的单对数坐标曲线, 由图可知, 应当选取曲线的中间部分进行拟合来求出信息维数.欺骗干扰和噪声的曲线与图2类似, 只是斜率不同.

图2 信息熵与尺度的单对数坐标曲线Fig.2 Semi-logarithmic curve of scale and comentropy

分别对欺骗干扰信号与目标信号分离和部分重叠两种情况进行仿真, 如图3所示, 信噪比(SNR)为5 dB, 干噪比(JNR)为8 dB.两种情况下的曲线峰值都在欺骗干扰的起始位置, 图3(a)的峰值在8.01 μ s, 图3(b)的峰值在4.06 μ s, 可见本文算法对欺骗干扰的检测精度很高, 且对欺骗干扰的存在位置不敏感.

图3 L值的时域分布图Fig.3 Temporal distribution of L

图4 L值随信噪比/干噪比的变化曲线Fig.4 Curve of L changing with SNR or JNR

图4仿真的是目标回波和欺骗干扰的L值与信噪比, 干噪比的关系, 其光滑曲线是由50次仿真平均得到的.由图4可知:目标回波的L值在SNR为8 dB时达到最大值0.141.考虑到这是50次仿真的平均结果, 单次计算会存在波动性, 通过仿真发现目标回波的L值在SNR为8 dB时的100次仿真中最大可达0.1674, 故将干扰检测门限 ε设为0.1675.在该检测门限下, 欺骗式干扰的检测概率曲线如图5所示.

图5 欺骗干扰的检测概率曲线Fig.5 Curve of detection probability of repeat jamming

图5可以看出:在JNR为6 dB时, 欺骗干扰的检测概率即可达到0.9.故当干噪比强于6 dB时, 本文方法可以有效地检测出欺骗干扰, 并易于进一步抑制干扰.还可以看出, 当干噪比弱于6 dB时, 只要目标回波的信噪比满足不低于6 dB, 则欺骗干扰不能构成威胁.故在信噪比不低于6 dB时, 本文算法可有效抗欺骗干扰, 而这个条件在实际应用环境中可以利用相参积累得到满足.

4 结束语

基于DRFM技术的欺骗式干扰信号在时, 频域均可与目标信号高度相似, 对雷达系统具有较强欺骗性.针对这种情况, 本文基于欺骗干扰与目标信号在信息维数上的不同, 提出了一种新的欺骗干扰存在性检测方法.通过仿真验证了本文算法的有效性, 仿真结果表明, 在干噪比不弱于6 dB时, 本文方法能有效检测欺骗干扰.本文并未研究后续的抗干扰算法, 而简单地置零处理有时也会造成目标信号的能量损失, 故下一步将研究在抑制欺骗干扰的同时不损失目标信号能量的抗干扰算法.

The authors have declared that no competing interests exist.

参考文献
[1] 卢刚. 雷达有源假目标抑制方法研究[D]. 成都: 电子科技大学电子工程学院, 2012.
Lu Gang. Study of algorithms on suppression of radar active false targets[D]. Chengdu: School of Electronic Engineering, University of Electronic Science and Technology of China, 2012. [本文引用:1]
[2] 宋志勇, 肖怀铁. 基于回波幅度特征的拖曳式诱饵存在性检测[J]. 电子与信息学报, 2011, 33(6): 1515-1519.
Song Zhi-yong, Xiao Huai-tie. Detection of presence of towed radar active decoy based on amplitude characteristic[J]. Journal of Electronics & Information Technology, 2011, 33(6): 1515-1519. [本文引用:2]
[3] Greco M, Gini F, Farina A, et al. Effect of phase and range gate pull-off delay quantization on jammer signal[J]. IEE Proceeding on Radar Sonar Navigation, 2006, 153(5): 454-459. [本文引用:1]
[4] 韩伟, 汤子跃, 朱振波, . 基于脉内特征的雷达抗欺骗干扰方法[J]. 雷达科学与技术, 2012, 10(5): 481-485.
Han Wei, Tang Zi-yue, Zhu Zhen-bo, et al. A method against radar's deception-jamming based on intra-pulse feature[J]. Radar Science and Technology, 2012, 10(5): 481-485. [本文引用:1]
[5] Li Jian-xun, Shen Qi, Yan Hai. Signal feature analysis and experimental verification of radar deception jamming[C]//Proceedings of 2011 IEEE CIE International Conference on Radar, Chengdu, 2011: 230-233. [本文引用:1]
[6] Cherouat S, Soltani F. Radar signal detection using fractal analysis in K-distributed clutter[C]//2008 International Conference on Signals, Monastir, 2008: 1-4. [本文引用:2]
[7] Guan J, Liu N B, Huang Y, et al. Fractal characteristic in frequency domain for target detection within sea clutter[J]. IET Radar, Sonar and Navigation, 2012, 6(5): 293-306. [本文引用:2]
[8] Zhang Kun-hua, Yang Xuan, Zhang Li. Target detection in complex background based on fuzzy enhanced fractal feature[C]//2009 IEEE International Conference on Information and Automation, Zhuhai, 2009: 408-411. [本文引用:2]
[9] 庞新良, 李鹏, 范江兵, . 小波计盒分形维数核爆电磁脉冲信号模糊识别[J]. 强激光与粒子束, 2012, 24(9): 2259-2263.
Pang Xin-liang, Li Peng, Fan Jiang-bing, et al. Fuzzy identification of NEMP signal based on wavelet and box-counting dimension eigenvalue[J]. High Power Laser and Particle Beams, 2012, 24(9): 2259-2263. [本文引用:1]
[10] 黄晚晴, 张颖, 王文义, . 利用分形法表征光学元件中高频相位畸变[J]. 强激光与粒子束, 2013, 25(5): 1171-1175.
Huang Wan-qing, Zhang Ying, Wang Wen-yi, et al. Evaluating middle-high frequency phase distortion of optics by fractal method[J]. High Power Laser and Particle Beams, 2013, 25(5): 1171-1175. [本文引用:1]
[11] Ye Fei, Yu Zhi-fu, Luo Jing-qing. Analysis of radar emitter signal feature based on mutifractal theory[C]//The Eighth International Conference on Electronic Measurement and Instruments, Xi'an, 2007: 114-117. [本文引用:1]