作者简介:祝宏(1987-),男,博士研究生.研究方向:雷达信号处理,雷达抗干扰技术.E-mail:zhu-h06@163.com
利用欺骗干扰信号与目标回波包络上的不同,基于分形理论中的信息维数提出了一种欺骗干扰的存在性检测方法.当信号本身的信息维数与信号包络的信息维数之差超过一定门限时,可判断存在欺骗干扰,同时可获知欺骗干扰的起始时刻.仿真结果表明:本文算法在干噪比不低于6 dB时,检测概率在90%以上.
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%.
抗欺骗式干扰技术一直是雷达研究者的关注重点.目前, 关于欺骗式干扰的研究主要集中在原理分析, 性能仿真和效果评估方面, 关于抗干扰的研究则相对较少[1, 2].文献[3]分析了数字射频存储器(DRFM)产生的距离波门拖引(RGPO)干扰信号由于DRFM相位量化误差在频域上产生的细微变化, 但这个细微变化在DRFM阶数较高时不明显, 且不适用于目标回波和干扰同时存在的情况.文献[4]分析了目标回波与欺骗干扰信号脉内幅度起伏的差异; 文献[5]研究了两者在幅度起伏, 高阶统计量, 双谱特性上的差异, 但均分析简单, 没有研究目标回波与欺骗干扰部分重叠的情况.
自分形理论被提出以来, 分形理论已成功地应用于图像处理, 语音识别, 雷达信号处理等领域.文献[6, 7, 8]分别从不同角度研究了分形理论在雷达目标检测中的应用, 但均不适用于存在欺骗干扰的情况.对此, 本文利用欺骗式干扰与目标回波在信息维数上的差异, 提出了一种新的欺骗干扰存在性检测方法.本文采用实测的大地回波数据进行仿真, 通过分析大地回波和欺骗干扰的信息维数特征, 得出判定欺骗干扰存在与否的依据, 为后续的抗欺骗干扰算法设计打下了良好基础.
分形维数是分形理论的中心概念, 它有多种定义和计算方法, 常用的是Hausdorff维数, 盒维数, 信息维数, 相似维数等[9, 10].雷达信号作为一种时间序列, 可以用分形维数对其进行刻画, 比如盒维数用以描述信号的复杂度, 信息维数用以描述信号在二维平面上的分布信息.本文采用后者对雷达信号进行研究.
设
实际计算中通常通过粗视化变换得到信息维数:
基于信息维数的定义, 本文设计了以下算法步骤进行计算:
Step1 数据预处理.首先, 对信号数据线性插值, 使得数据长度为2的整数次幂, 设插值后的数据长度为Length; 然后, 将每个数据点都减去这段信号的最小值, 使得这段信号数据的最小值为0, 设处理后的数据最大值为Data_max; 最后, 将每个数据同乘以Length/Data_max, 使得这段信号数据的最大值等于数据长度.以盒子覆盖信号时, 计算落入每个盒子的信号曲线长度比较麻烦, 所以本文采用计算落入每个盒子的数据点数来替代, 故而需要对信号数据进行插值, 以提高计算精度.使数据长度为2的整数次幂的目的是为了便于设定盒子尺寸, 以及便于运算.对信号的最小值和最大值预处理的目的是为了便于摆放盒子, 且能够满足分形维数与信号幅度无关的特点.
Step2 计算信息熵.依照信息维数的定义, 计算不同盒子尺度下的信息熵.盒子尺度依次选取2的整数次幂, 且不大于数据长度.
Step3 选取无标度区间.作
Step4 计算信息维数.对无标度区间内的曲线进行线性拟合, 得到的斜率即为信息维数.
欺骗式干扰主要是通过DRFM对接收到的雷达发射脉冲进行采集, 存储, 调制放大, 转发来产生的, 从而逼真地模拟目标的相关特性, 使得二者在时域, 频域上的差异很小.随着现代数字处理软硬件技术的发展, 一些新型的干扰机在转发信号的同时能够结合脉内调制, RCS调制等复合手段, 实现欺骗信号与目标回波的高度相似.然而这些新型干扰机目前大多还处于研制和实验阶段, 并未大量装备[2].因此, 本文从实际需要出发, 只针对传统的转发式欺骗干扰, 对于新型欺骗干扰则没有过多考虑.由于雷达的发射脉冲往往是恒包络的, 所以转发式欺骗干扰的包络也是恒定的, 而目标回波的幅度受到目标散射面的调制, 其包络是不规则的, 如图1所示.
由图1可知:大地回波的包络与欺骗干扰的包络具有明显差异, 所以它们的包络的分形特征也会有明显差异, 而一般的目标分形检测并没有对信号包络的分形维数进行分析.对于存在欺骗干扰的数据段, 由于干扰信号的包络是个直流电平, 由分形理论的基本性质和信息维数的计算流程均不难看出该数据段的信息维数与噪声比较接近.在信噪比较大时, 存在目标回波的数据段信息维数将小于噪声的信息维数.经过研究发现, 如果仅仅利用对信号包络求信息维数并取最小值来检测目标回波, 其性能很差, 一是因为对信噪比的要求很高; 二是因为当欺骗干扰较强时, 在目标回波和欺骗干扰之间会存在极小值, 影响检测的准确性.为此, 本文先逐段分别计算信号包络和信号本身的信息维数, 然后求两个信息维数的差.前面的分析已经指出, 欺骗干扰的两个信息维数差值比较大, 而噪声, 目标回波的则相对比较小, 故可以据此判断是否存在欺骗干扰, 并且检测出干扰的存在位置.
依次以每个信号数据点为起点, 划分长度为发射脉冲宽度的数据段, 计算每个数据段
式中:
(1)当
(2)当
本文采用雷达高度表实测的大地回波数据作为目标回波.选取信噪比很好的某个大地回波, 将其下变频到零中频(消除载波影响), 作为干净的目标回波, 其发射脉冲是线性调频信号, 带宽为10 MHz, 脉宽为5 μ s.由于测量时挂飞的直升机垂直速度很小, 且多普勒频率对分形维数影响很小, 故本文不考虑目标回波的多普勒频率.欺骗式干扰的参数与发射脉冲参数相同, 多普勒频率也假设为零.构造一个长度为20 μ s的雷达接收数据, 其中目标回波起始时刻是1 μ s, 欺骗干扰的起始时刻是4 μ s(与目标回波部分重叠)或者8 μ s(与目标回波分离), 噪声为高斯白噪声.
分形理论所研究的系统是建立在自相似基础之上的, 对于实际的系统, 这只可能在一定的区间内成立, 即无标度区间, 超出这个区间的自相似性就不存在了, 分形也就失去了意义.图2给出了大地回波的信息熵与尺度的单对数坐标曲线, 由图可知, 应当选取曲线的中间部分进行拟合来求出信息维数.欺骗干扰和噪声的曲线与图2类似, 只是斜率不同.
分别对欺骗干扰信号与目标信号分离和部分重叠两种情况进行仿真, 如图3所示, 信噪比(SNR)为5 dB, 干噪比(JNR)为8 dB.两种情况下的曲线峰值都在欺骗干扰的起始位置, 图3(a)的峰值在8.01 μ s, 图3(b)的峰值在4.06 μ s, 可见本文算法对欺骗干扰的检测精度很高, 且对欺骗干扰的存在位置不敏感.
图4仿真的是目标回波和欺骗干扰的L值与信噪比, 干噪比的关系, 其光滑曲线是由50次仿真平均得到的.由图4可知:目标回波的L值在SNR为8 dB时达到最大值0.141.考虑到这是50次仿真的平均结果, 单次计算会存在波动性, 通过仿真发现目标回波的L值在SNR为8 dB时的100次仿真中最大可达0.1674, 故将干扰检测门限
由图5可以看出:在JNR为6 dB时, 欺骗干扰的检测概率即可达到0.9.故当干噪比强于6 dB时, 本文方法可以有效地检测出欺骗干扰, 并易于进一步抑制干扰.还可以看出, 当干噪比弱于6 dB时, 只要目标回波的信噪比满足不低于6 dB, 则欺骗干扰不能构成威胁.故在信噪比不低于6 dB时, 本文算法可有效抗欺骗干扰, 而这个条件在实际应用环境中可以利用相参积累得到满足.
基于DRFM技术的欺骗式干扰信号在时, 频域均可与目标信号高度相似, 对雷达系统具有较强欺骗性.针对这种情况, 本文基于欺骗干扰与目标信号在信息维数上的不同, 提出了一种新的欺骗干扰存在性检测方法.通过仿真验证了本文算法的有效性, 仿真结果表明, 在干噪比不弱于6 dB时, 本文方法能有效检测欺骗干扰.本文并未研究后续的抗干扰算法, 而简单地置零处理有时也会造成目标信号的能量损失, 故下一步将研究在抑制欺骗干扰的同时不损失目标信号能量的抗干扰算法.
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|