CN112083393B - 基于谱图平均时间特征的间歇采样转发干扰识别方法 - Google Patents

基于谱图平均时间特征的间歇采样转发干扰识别方法 Download PDF

Info

Publication number
CN112083393B
CN112083393B CN202011160606.3A CN202011160606A CN112083393B CN 112083393 B CN112083393 B CN 112083393B CN 202011160606 A CN202011160606 A CN 202011160606A CN 112083393 B CN112083393 B CN 112083393B
Authority
CN
China
Prior art keywords
isdj
signal
isrj
iscj
intermittent sampling
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.)
Active
Application number
CN202011160606.3A
Other languages
English (en)
Other versions
CN112083393A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202011160606.3A priority Critical patent/CN112083393B/zh
Publication of CN112083393A publication Critical patent/CN112083393A/zh
Application granted granted Critical
Publication of CN112083393B publication Critical patent/CN112083393B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提出了一种基于谱图平均时间特征的间歇采样转发干扰识别方法，旨在提高低干噪比下间歇采样转发干扰的识别概率，实现步骤为：获取时域雷达检测信号集；提取谱图平均时间特征作为时域雷达检测信号集的特征数据集；构造训练数据集和测试数据集；训练随机森林模型；获取间歇采样转发干扰识别结果。本发明在进行间歇采样转发干扰的识别时，所提取的谱图平均时间特征在低干噪比下对噪声的敏感性低，在各类信号之间的区分性大，具有良好的干扰识别效果，同时，仅使用谱图平均时间特征用于干扰信号的分类，所采用的随机森林训练方法的学习过程较为快速且能产生高准确度的分类器具有实时性高，识别效率高的特点。

Description

基于谱图平均时间特征的间歇采样转发干扰识别方法
技术领域
本发明属于雷达对抗干扰技术领域，涉及一种间歇采样转发干扰识别方法，特别涉及一种对间歇采样转发干扰信号的雷达干扰识别方法，可用于雷达对间歇采样转发干扰信号的对抗识别。
背景技术
基于数字射频存储(Digital Radio Frequency Memory，DRFM)技术的间歇采样转发干扰(Interrupted Sampling and Repeater Jamming)是针对线性调频脉冲压缩雷达而提出的，干扰机通过对接收到的雷达信号进行间歇性“欠采样”处理，基于天线收发分时体制，巧妙地利用脉冲压缩雷达的脉内相干性，复制出与雷达发射信号具有相当大的相干处理增益的干扰信号，可以产生多个逼真假目标，在一定条件下还可以产生压制干扰的效果，为雷达的检测与跟踪带来巨大的挑战。随着应用逐渐广泛，间歇采样转发干扰得到不断改进，在间歇采样直接转发干扰(Interrupted Sampling and Direct Repeater Jamming,ISDJ)的基础上，又发展出间歇采样重复转发干扰(Interrupted Sampling and RepeatRepeater Jamming,ISRJ)和间歇采样循环转发干扰(Interrupted Sampling and CircleRepeater Jamming,ISCJ)样式，可以实现对多种新体制雷达的干扰。间歇采样转发干扰是目前一种主流的干扰样式，一方面可以对雷达进行噪声压制，同时也可以产生多假目标进行速度、距离欺骗。
间歇采样转发干扰包括间歇采样直接转发干扰、间歇采样重复转发干扰和间歇采样循环转发干扰三种样式。间歇采样直接转发干扰是指干扰机截获到雷达信号后，高保真采样其中一小段信号并马上进行转发，然后再采样、转发下一段，采样与转发分时交替工作。间歇采样重复转发干扰是指干扰机采样到雷达发射的一小段信号后，按照干扰机设定的次数重复读出当前采样信号并转发，然后再采样一小段信号并重复进行转发，如此重复上述过程，直到雷达脉冲结束。间歇采样循环转发干扰是指第1次采样后，转发第1次采样信号，接着进行第2次采样，转发第2次采样信号后再转发存储的第1次采样信号，第3次采样后转发第3次采样信号后再把存储的前2次信号按采样顺序进行逆序转发，依此类推，直到雷达脉冲结束。
干扰识别是抗干扰流程中的关键环节，采取有效抗干扰措施的前提是对干扰类型进行正确的分类识别。对干扰类型进行分类识别主要从识别准确率和识别效率两个方面考虑，一方面要提高干扰识别准确率，只有对干扰类型进行准确的识别，后续才能提供精准有效的干扰对抗方式，另一方面要提高干扰识别效率，干扰识别过程中必须能够快速准确的识别干扰类型，具有实时性，以应对瞬息万变的电磁干扰环境。
目前关于雷达干扰识别的已有方法大多是从信号处理的角度出发，提取不同的干扰信号在多个变换域上的特征参数，得到干扰信号幅度、相位、频率、能量、波形等多方面的信息，这些信息反映了不同种类信号之间的差异性。一般随着噪声功率增大，代表信号特点的特征参数越不容易提取，当噪声功率增强到一定程度时，信号将淹没在噪声之中，特征参数不能被提取。雷达一般工作在复杂的电磁环境中，因此，为了提高干扰信号的识别准确率，特征参数应该一方面要对噪声的敏感性要低，即特征参数应该随干噪比的增大变化不大，另一方面特征参数在各类信号之间的区分性大。
比如，周超等人于2017年在《信号处理》上发表了一篇题为“DRFM间歇采样转发干扰辨识方法研究”的文章，公开了一种基于滑动截断匹配滤波的间歇采样转发干扰识别方法，该方法首先通过对匹配滤波器的参考窗宽度和延时进行二维搜索，输出脉压后的二维幅度分布；然后基于该幅度分布对干扰切片宽度和转发周期进行估计；并通过分析切片宽度与转发周期之间的关系，实现典型转发式干扰的辨识。但是，该方法仅在干噪比大于5dB时，正确识别概率不低于90％，而在干噪比较小时识别准确率较低。
发明内容
本发明的目的在于克服上述现有技术存在的缺陷，提出了一种基于谱图平均时间特征的间歇采样转发干扰识别方法，用于解决现有技术中存在的低干噪比下识别准确率较低的技术问题。
为实现上述目的，本发明采取的技术方案包括如下步骤：
(1)获取时域雷达检测信号集X：
(1a)设雷达发射机发射的L个信号为s0(t)＝{s0(t1)，s0(t2)，…,s0(tl)，…,s0(tL)}，雷达接收机接收s0(t)经目标反射后包含有高斯白噪声的目标信号为S(t)＝{S(t′1),S(t'2),…,S(t'l),…,S(t'L)}，其中，s0(tl)表示第l个发射信号，S(t′l)表示s0(tl)反射的目标信号：
Figure BDA0002744077760000031
S(t′l)＝n(t′l)+s0(t′l)exp(jwdt′l)
其中，L≥1000，tl表示s0(tl)的发射时间，T表示s0(tl)的脉冲宽度，
Figure BDA0002744077760000032
表示矩形窗函数，/>
Figure BDA0002744077760000033
exp[·]表示以e为底的指数函数，j表示虚数单位，ws表示载频，π表示圆周率，k表示脉冲调频斜率，/>
Figure BDA0002744077760000034
B表示s0(tl)的带宽，t′l＝tl-Δt表示S(t′l)的时间，Δt表示时延，/>
Figure BDA0002744077760000035
R0表示目标与雷达之间的距离，c表示光速，n(t′l)表示高斯白噪声，wd表示由目标运动速度引起的多普勒频移；
(1b)设雷达接收机接收到的含有高斯白噪声的间歇采样直接转发干扰信号为ISDJ(t)＝{ISDJ(t1),ISDJ(t2),…,ISDJ(tl),…,ISDJ(tL)}、间歇采样重复转发干扰信号为ISRJ(t)＝{ISRJ(t1),ISRJ(t2),…,ISRJ(tl),…,ISRJ(tL)}和间歇采样循环转发干扰信号为ISCJ(t)＝{ISCJ(t1),ISCJ(t2),…,ISCJ(tl),…,ISCJ(tL)},其中，ISDJ(tl)、ISRJ(tl)、ISCJ(tl)分别表示第l个间歇采样直接转发干扰信号、间歇采样重复转发干扰信号和间歇采样循环转发干扰信号：
ISDJ(tl)＝n(tl)+xs(tl-τ)
Figure BDA0002744077760000036
Figure BDA0002744077760000037
其中，xs(tl)表示目标所携带干扰机对接收到的第l个信号s0'(tl)进行间歇采样所获取的间歇采样信号，xs(tl)＝xs1(tl)+xs2(tl)+…+xsi(tl)+…+xsN(tl)，s0'(tl)＝s0(tl)exp(jwdtl)，xsi(tl)表示第i次间歇采样得到的子脉冲信号，i＝0,1,…,N，N表示T内的采样次数，
Figure BDA0002744077760000041
Figure BDA0002744077760000042
表示向下取整，Ts表示间歇采样的周期，τ表示间歇采样的时长，Σ表示求和，M表示Ts内对xsi(tl)转发的次数，/>
Figure BDA0002744077760000043
R＝min{M,N}；
(1c)对步骤(1a)得到的S(t)和步骤(1b)得到的ISDJ(t)、ISCJ(t)、ISRJ(t)进行合并，得到包含四类共4×L个信号的时域雷达检测信号集X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}；
(2)提取时域雷达检测信号集X的特征数据集D：
(2a)设时域雷达检测信号集X中的任意一个时域雷达检测信号为z(t)，并对每个z(t)进行短时傅里叶变换，得到4×L个频域信号，其中z(t)的频域信号STFTz(t)(t,w)表达式为：
Figure BDA0002744077760000044
其中w表示短时傅里叶变化的频率频点，g(u-t)表示短时傅里叶变换中对z(t)加的窗函数，
Figure BDA0002744077760000045
表示短时傅里叶变换；
(2b)对频域信号STFTz(t)(t,w)取模值后取平方，得到4×L个信号频谱图，其中z(t)的信号频谱图SPECz(t)(t,w)的表达式为：
SPECz(t)(t,w)＝|STFTz(t)(t,w)|2
其中|·|表示取模值；
(2c)根据信号频谱图SPECz(t)(t,w)在(-∞,+∞)的整体平均时间，计算z(t)的谱图平均时间特征tSP_z(t)，得到X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}对应的谱图平均时间特征集tSP＝{tSP_Sl,tSP_ISDJl,tSP_ISCJl,tSP_ISRJl}，其中：
Figure BDA0002744077760000051
其中∫∫表示二重积分；
(2d)以
Figure BDA0002744077760000052
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，得到时域雷达检测信号集X的特征数据集D，D＝{tySP_Sl,tySP_ISDJl,tySP_ISCJl,tySP_ISRJl}；
以
Figure BDA0002744077760000053
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，其中的真实标签值为任意实数，且同一类型干扰的标签值相同，不同类型干扰的标签值不同。
(3)构造训练数据集E和测试数据集T：
将从特征数据集D中随机选取的每类信号的V个特征数据组成训练数据集E＝{ty_Sv,ty_ISDJv,ty_ISCJv,ty_ISRJv}，每类信号其余的L-V个特征数据组成包含4×(L-V)个特征数据的测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}，其中，1≤v≤V，1≤x≤L-V；
(4)训练随机森林模型：
(4a)以决策树为基学习器，并使用Bagging集合策略构建包含NUM棵决策树的随机森林模型RF＝{Tree1,Tree2,…,Treenum,…,TreeNUM}，其中，Treenum表示第num棵决策树，NUM≥50，1≤num≤NUM；
(4b)令num＝1；
(4c)从训练数据集E中独立地、随机地、有放回地抽取出p个特征数据作为随机森林模型的输入，对决策树Treenum的根节点进行节点分裂，得到2个分支节点，直到每个分支节点的所有训练数据都属于同一类型，得到训练好的决策树Tree'num，其中p＜P，P＝4×V表示训练数据集E中的特征数据总数；
(4d)判断num＝NUM是否成立，若是，得到训练好的随机森林模型，否则，令num＝num+1，并执行步骤(4c)；
(5)获取间歇采样转发干扰识别结果：
(5a)将测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}中的4×(L-V)个特征数据去除真实类别标签
Figure BDA0002744077760000061
得到用于分类识别的数据，并将其输入训练好的随机森林分类模型进行测试，得到对应的4×(L-V)个预测类别标签
Figure BDA0002744077760000062
(5b)将4×(L-V)个预测类别标签
Figure BDA0002744077760000063
与真实类别标签
Figure BDA0002744077760000064
进行对比，得到目标信号S(t)的识别准确率/>
Figure BDA0002744077760000065
间歇采样直接转发干扰信号ISDJ(t)的识别准确率/>
Figure BDA0002744077760000066
间歇采样重复转发干扰信号ISCJ(t)的识别准确率/>
Figure BDA0002744077760000067
和间歇采样循环转发干扰信号ISRJ(t)的识别准确率
Figure BDA0002744077760000068
其中RS、RISDJ、RISCJ、RISRJ分别表示S(t)、ISDJ(t)、ISCJ(t)、ISRJ(t)的特征数据中预测类别标签和真实类别标签一致的特征数据的个数。
本发明与现有技术相比，具有以下优点：
1、本发明在对间歇采样转发干扰的识别过程中，首先使用海明窗对信号进行短时傅里叶变换得到信号谱图，提取谱图平均时间作为信号的特征，该谱图平均时间特征在低干噪比下对噪声的敏感性低，在各类信号之间的区分性大，因此具有良好的干扰识别效果，解决了现有技术中低干噪比下信号识别率低的问题，提高了间歇采样转发干扰信号的识别准确率。
2、本发明仅使用一种特征参数用于干扰信号的分类识别，所采用的随机森林模型训练学习过程较为快速且能产生高准确度的分类结果，因此本发明具有实时性高，识别效率高的特点。
附图说明
图1是本发明的实现流程图；
图2是本发明谱图平均时间随干噪比变化曲线仿真图；
图3是用本发明中4种信号的识别准确率随干噪比变化曲线仿真图；
图4是现有技术与本发明平均识别准确率的仿真对比图。
具体实施方式
下面结合附图和具体实施例，对本发明作进一步详细描述：
参照图1，本发明包括如下步骤：
步骤1)获取时域雷达检测信号集X：
步骤1a)设雷达发射机发射的L个信号为s0(t)＝{s0(t1)，s0(t2)，…,s0(tl)，…,s0(tL)}，雷达接收机接收s0(t)经目标反射后包含有高斯白噪声的目标信号为S(t)＝{S(t′1),S(t'2),…,S(t'l),…,S(t'L)}，其中，L≥1000，s0(tl)表示第l个发射信号，S(t′l)表示s0(tl)反射的目标信号，本实施例中L＝2000：
Figure BDA0002744077760000071
S(t′l)＝n(t′l)+s0(t′l)exp(jwdt′l)
其中，tl表示s0(tl)的发射时间，T表示s0(tl)的脉冲宽度，
Figure BDA0002744077760000072
表示矩形窗函数，/>
Figure BDA0002744077760000073
exp[·]表示以e为底的指数函数，j表示虚数单位，ws表示载频，π表示圆周率，k表示脉冲调频斜率，/>
Figure BDA0002744077760000074
B表示s0(tl)的带宽，t′l＝tl-Δt表示S(t′l)的时间，Δt表示时延，/>
Figure BDA0002744077760000075
R0表示目标与雷达之间的距离，c表示光速，n(t′l)表示高斯白噪声，wd表示由目标运动速度引起的多普勒频移；
步骤1b)设雷达接收机接收到的含有高斯白噪声的间歇采样直接转发干扰信号为ISDJ(t)＝{ISDJ(t1),ISDJ(t2),…,ISDJ(tl),…,ISDJ(tL)}、间歇采样重复转发干扰信号为ISRJ(t)＝{ISRJ(t1),ISRJ(t2),…,ISRJ(tl),…,ISRJ(tL)}和间歇采样循环转发干扰信号为ISCJ(t)＝{ISCJ(t1),ISCJ(t2),…,ISCJ(tl),…,ISCJ(tL)},其中，ISDJ(tl)、ISRJ(tl)、ISCJ(tl)分别表示第l个间歇采样直接转发干扰信号、间歇采样重复转发干扰信号和间歇采样循环转发干扰信号：
ISDJ(tl)＝n(tl)+xs(tl-τ)
Figure BDA0002744077760000081
Figure BDA0002744077760000082
其中，xs(tl)表示目标所携带干扰机对接收到的第l个信号s0'(tl)进行间歇采样所获取的间歇采样信号，xs(tl)＝xs1(tl)+xs2(tl)+…+xsi(tl)+…+xsN(tl)，
s0'(tl)＝s0(tl)exp(jwdtl)，xsi(tl)表示第i次间歇采样得到的子脉冲信号，i＝0,1,…,N，N表示T内的采样次数，
Figure BDA0002744077760000083
Figure BDA0002744077760000084
表示向下取整，Ts表示间歇采样的周期，Σ表示求和，M表示Ts内对xsi(tl)转发的次数，/>
Figure BDA0002744077760000085
τ表示Ts的间歇采样时长，R＝min{M,N}；
间歇采样直接转发干扰是指目标所携带干扰机接收到的第l个信号s0'(tl)后，在间歇采样周期Ts内，采样一段时长为τ的信号并马上进行转发，然后在下一个间歇采样周期Ts内再采样、转发，采样与转发分时交替工作；间歇采样重复转发干扰是指目标所携带干扰机采样到一段时长为τ的信号后，按照干扰机设定的次数M重复读出当前采样信号并转发，然后在下一个间歇采样周期Ts内再采样一段时长为τ的信号信号并重复进行转发，如此重复上述过程，直到雷达脉冲结束；间歇采样循环转发干扰是指第1次采样后，转发第1次采样信号，接着进行第2次采样，转发第2次采样信号后再转发存储的第1次采样信号，第3次采样后转发第3次采样信号后再把存储的前2次信号按采样顺序进行逆序转发，依此类推，直到雷达脉冲结束。
步骤1c)对步骤1a)得到的S(t)和步骤1b)得到的ISDJ(t)、ISCJ(t)、ISRJ(t)进行合并，得到包含四类共4×L个信号的时域雷达检测信号集X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}；
步骤2)提取时域雷达检测信号集X的特征数据集D：
步骤2a)设时域雷达检测信号集X中的任意一个时域雷达检测信号为z(t)，并对每个z(t)进行短时傅里叶变换，得到4×L个频域信号，其中z(t)的频域信号STFTz(t)(t,w)表达式为：
Figure BDA0002744077760000091
其中w表示短时傅里叶变化的频率频点，g(u-t)表示短时傅里叶变换中对z(t)加的窗函数，
Figure BDA0002744077760000092
表示短时傅里叶变换；
时域雷达检测信号集X中的任意一个时域雷达检测信号为z(t)，都是包含高斯白噪声的时域信号，属于非平稳信号；
通过短时傅里叶变换，可以将含有高斯白噪声的非平稳信号z(t)看成是由一系列短时平稳信号构成的，对窗内信号做傅立叶变换，得到非平稳信号的时变频谱图；
当选择的窗函数在时间上对称时，就能根据频谱图正确给出信号的波形中心和信号频率，本实施例中选择海明窗作为短时傅里叶变换中的窗函数；
步骤2b)对频域信号STFTz(t)(t,w)取模值后取平方，得到4×L个信号频谱图，其中z(t)的信号频谱图SPECz(t)(t,w)的表达式为：
SPECz(t)(t,w)＝|STFTz(t)(t,w)|2
其中|·|表示取模值；
步骤2c)根据信号频谱图SPECz(t)(t,w)在(-∞,+∞)的整体平均时间，计算z(t)的谱图平均时间特征tSP_z(t)，得到X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}对应的谱图平均时间特征集tSP＝{tSP_Sl,tSP_ISDJl,tSP_ISCJl,tSP_ISRJl}，其中：
Figure BDA0002744077760000093
其中∫∫表示二重积分；
步骤2d)以
Figure BDA0002744077760000101
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，得到时域雷达检测信号集X的特征数据集D，D＝{tySP_Sl,tySP_ISDJl,tySP_ISCJl,tySP_ISRJl}；
以
Figure BDA0002744077760000102
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，其中的真实标签值为任意实数，且同一类型干扰的L个信号的谱图平均时间特征标签值相同，不同类型干扰的标签值不同；
本实施例设置真实类别标签
Figure BDA0002744077760000103
步骤3)构造训练数据集E和测试数据集T：
将从特征数据集D中随机选取的每类信号的V个特征数据组成训练数据集E＝{ty_Sv,ty_ISDJv,ty_ISCJv,ty_ISRJv}，每类信号其余的L-V个特征数据组成包含4×(L-V)个特征数据的测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}，其中，1≤v≤V，1≤x≤L-V，本实施例中V＝500；
步骤4)训练随机森林模型：
步骤4a)以决策树为基学习器，并使用Bagging集合策略构建包含NUM棵决策树的随机森林模型RF＝{Tree1,Tree2,…,Treenum,…,TreeNUM}，其中，Treenum表示第num棵决策树，NUM≥50，1≤num≤NUM，本实施例设置NUM＝50；
步骤4b)令num＝1；
步骤4c)从训练数据集E中独立地、随机地、有放回地抽取出p个特征数据作为随机森林模型的输入，对决策树Treenum的根节点进行节点分裂，得到2个分支节点，直到每个分支节点的所有训练数据都属于同一类型，得到训练好的决策树Tree'num，其中p＜P，P＝4×V表示训练数据集E中的特征数据总数；
步骤4d)判断num＝NUM是否成立，若是，得到训练好的随机森林模型，否则，令num＝num+1，并执行步骤(4c)；
步骤5)获取间歇采样转发干扰识别结果：
步骤5a)将测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}中的4×(L-V)个特征数据去除真实类别标签
Figure BDA0002744077760000111
得到用于分类识别的数据，并将其输入训练好的随机森林分类模型进行测试，得到对应的4×(L-V)个预测类别标签
Figure BDA0002744077760000112
步骤5b)将4×(L-V)个预测类别标签
Figure BDA0002744077760000113
与真实类别标签
Figure BDA0002744077760000114
进行对比，得到目标信号S(t)的识别准确率/>
Figure BDA0002744077760000115
间歇采样直接转发干扰信号ISDJ(t)的识别准确率/>
Figure BDA0002744077760000116
间歇采样重复转发干扰信号ISCJ(t)的识别准确率/>
Figure BDA0002744077760000117
和间歇采样循环转发干扰信号ISRJ(t)的识别准确率/>
Figure BDA0002744077760000118
其中RS、RISDJ、RISCJ、RISRJ分别表示S(t)、ISDJ(t)、ISCJ(t)、ISRJ(t)的特征数据中预测类别标签和真实类别标签一致的特征数据的个数。
下面结合仿真实验，对本发明的技术效果作以说明。
1、仿真条件和内容：
雷达发射信号为线性调频信号，发射L＝2000个信号，脉冲宽度T＝100μs，载频ws＝12e6Hz，带宽B＝10e6Hz，雷达采样频率fs＝30e6Hz，目标与雷达之间的距离R0＝10km，由目标运动速度引起的多普勒频移wd＝10Hz；目标所携带干扰机进行间歇采样时间歇采样周期Ts＝10μs，间歇采样时长τ＝2.5μs，采样次数N＝10，在Ts内对子脉冲转发的次数M＝3，R＝3；高斯白噪声的噪声功率与干扰信号的功率的比值为干噪比，令干噪比在-6dB到21dB之间，间隔3dB进行仿真实验；在提取信号的谱图平均时间特征时，使用长度为128个采样点海明窗，重叠点数为50％。
仿真1，在-6dB到21dB之间，间隔3dB提取L＝2000个目标回波信号、间歇采样直接转发信号、间歇采样重复转发信号和间歇采样循环转发信号得谱图平均时间特征，将这2000个谱图平均时间特征取平均值，得到图2谱图平均时间随干噪比变化仿真曲线。
仿真2，在-6dB到21dB之间，间隔3dB提取L＝2000个目标回波信号、间歇采样直接转发信号、间歇采样重复转发信号和间歇采样循环转发信号得谱图平均时间特征，将每种信号500个特征加上真实类型标签作为随机森林模型的输入数据训练随机森林，将剩下的数据用来测试，得到图3为用本发明中4种信号的识别准确率随干噪比变化仿真曲线。
仿真3，在-6dB到21dB之间，间隔3dB提取L＝2000个目标回波信号、间歇采样直接转发信号、间歇采样重复转发信号和间歇采样循环转发信号得谱图平均时间特征，将每种信号500个特征加上真实类型标签作为随机森林模型的输入数据训练随机森林，将剩下的数据用来测试，将得到的4种型号在不同干噪比下的识别准确率取平均值，得到本发明方法的识别准确率随干噪比变化的曲线，与现有方法相比得到图4是现有技术与本发明平均识别准确率的仿真对比图。
仿真过程中软硬件环境:
硬件环境：CPU为Inter Core i7-8750H，主频为2.2Ghz，主存为8GB。
软件环境：Windows 10企业版，MATLAB 2020a仿真软件。
2、结果分析：
图2为谱图平均时间随干噪比变化仿真曲线，从图2可以看出，在不同的干噪比下，几种信号的谱图平均时间特征在数值上明显区分，所以此特征是一个噪声不敏感的特征且对间歇采样转发干扰的识别效果会更好。
图3为用本发明中4种信号的识别准确率随干噪比变化仿真曲线，在0dB时四种信号的识别准确率均能达到95％以上，在-6dB时，识别准确率最低的间歇采样循环转发干扰也能达到68％左右的识别准确率。
图4是现有技术与本发明平均识别准确率的仿真对比图，从现有方式和本发明的识别准确率对比曲线来看，在低干噪比下，本发明的平均识别率要高于现有方式。

Claims (2)

1.一种基于谱图平均时间特征的间歇采样转发干扰识别方法，其特征在于，包括以下步骤：
(1)获取时域雷达检测信号集X：
(1a)设雷达发射机发射的L个信号为s0(t)＝{s0(t1)，s0(t2)，…,s0(tl)，…,s0(tL)}，雷达接收机接收s0(t)经目标反射后包含有高斯白噪声的目标信号为S(t)＝{S(t′1),S(t'2),…,S(t'l),…,S(t'L)}，其中，s0(tl)表示第l个发射信号，S(t′l)表示s0(tl)反射的目标信号：
Figure FDA0002744077750000011
S(t′l)＝n(t′l)+s0(t′l)exp(jwdt′l)
其中，L≥1000，tl表示s0(tl)的发射时间，T表示s0(tl)的脉冲宽度，
Figure FDA0002744077750000012
表示矩形窗函数，/>
Figure FDA0002744077750000013
exp[·]表示以e为底的指数函数，j表示虚数单位，ws表示载频，π表示圆周率，k表示脉冲调频斜率，/>
Figure FDA0002744077750000014
B表示s0(tl)的带宽，t′l＝tl-Δt表示S(t′l)的时间，Δt表示时延，/>
Figure FDA0002744077750000015
R0表示目标与雷达之间的距离，c表示光速，n(t′l)表示高斯白噪声，wd表示由目标运动速度引起的多普勒频移；
(1b)设雷达接收机接收到的含有高斯白噪声的间歇采样直接转发干扰信号为ISDJ(t)＝{ISDJ(t1),ISDJ(t2),…,ISDJ(tl),…,ISDJ(tL)}、间歇采样重复转发干扰信号为ISRJ(t)＝{ISRJ(t1),ISRJ(t2),…,ISRJ(tl),…,ISRJ(tL)}和间歇采样循环转发干扰信号为ISCJ(t)＝{ISCJ(t1),ISCJ(t2),…,ISCJ(tl),…,ISCJ(tL)},其中，ISDJ(tl)、ISRJ(tl)、ISCJ(tl)分别表示第l个间歇采样直接转发干扰信号、间歇采样重复转发干扰信号和间歇采样循环转发干扰信号：
ISDJ(tl)＝n(tl)+xs(tl-τ)
Figure FDA0002744077750000021
Figure FDA0002744077750000022
其中，xs(tl)表示目标所携带干扰机对接收到的第l个信号s0'(tl)进行间歇采样所获取的间歇采样信号，xs(tl)＝xs1(tl)+xs2(tl)+…+xsi(tl)+…+xsN(tl)，s0'(tl)＝s0(tl)exp(jwdtl)，xsi(tl)表示第i次间歇采样得到的子脉冲信号，i＝0,1,…,N，N表示T内的采样次数，
Figure FDA0002744077750000023
Figure FDA0002744077750000024
表示向下取整，Ts表示间歇采样的周期，τ表示间歇采样的时长，Σ表示求和，M表示Ts内对xsi(tl)转发的次数，/>
Figure FDA0002744077750000025
R＝min{M,N}；
(1c)对步骤(1a)得到的S(t)和步骤(1b)得到的ISDJ(t)、ISCJ(t)、ISRJ(t)进行合并，得到包含四类共4×L个信号的时域雷达检测信号集X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}；
(2)提取时域雷达检测信号集X的特征数据集D：
(2a)设时域雷达检测信号集X中的任意一个时域雷达检测信号为z(t)，并对每个z(t)进行短时傅里叶变换，得到4×L个频域信号，其中z(t)的频域信号STFTz(t)(t,w)表达式为：
Figure FDA0002744077750000026
其中w表示短时傅里叶变化的频率频点，g(u-t)表示短时傅里叶变换中对z(t)加的窗函数，
Figure FDA0002744077750000027
表示短时傅里叶变换；
(2b)对频域信号STFTz(t)(t,w)取模值后取平方，得到4×L个信号频谱图，其中z(t)的信号频谱图SPECz(t)(t,w)的表达式为：
SPECz(t)(t,w)＝|STFTz(t)(t,w)|2
其中|·|表示取模值；
(2c)根据信号频谱图SPECz(t)(t,w)在(-∞,+∞)的整体平均时间，计算z(t)的谱图平均时间特征tSP_z(t)，得到X＝{S(t),ISDJ(t),ISCJ(t),ISRJ(t)}对应的谱图平均时间特征集tSP＝{tSP_Sl,tSP_ISDJl,tSP_ISCJl,tSP_ISRJl}，其中：
Figure FDA0002744077750000031
其中∫∫表示二重积分；
(2d)以
Figure FDA0002744077750000032
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，得到时域雷达检测信号集X的特征数据集D，D＝{tySP_Sl,tySP_ISDJl,tySP_ISCJl,tySP_ISRJl}；
(3)构造训练数据集E和测试数据集T：
将从特征数据集D中随机选取的每类信号的V个特征数据组成训练数据集E＝{ty_Sv,ty_ISDJv,ty_ISCJv,ty_ISRJv}，每类信号其余的L-V个特征数据组成包含4×(L-V)个特征数据的测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}，其中，1≤v≤V，1≤x≤L-V；
(4)训练随机森林模型：
(4a)以决策树为基学习器，并使用Bagging集合策略构建包含NUM棵决策树的随机森林模型RF＝{Tree1,Tree2,…,Treenum,…,TreeNUM}，其中，Treenum表示第num棵决策树，NUM≥50，1≤num≤NUM；
(4b)令num＝1；
(4c)从训练数据集E中独立地、随机地、有放回地抽取出p个特征数据作为随机森林模型的输入，对决策树Treenum的根节点进行节点分裂，得到2个分支节点，直到每个分支节点的所有训练数据都属于同一类型，得到训练好的决策树Tree'num，其中p＜P，P＝4×V表示训练数据集E中的特征数据总数；
(4d)判断num＝NUM是否成立，若是，得到训练好的随机森林模型，否则，令num＝num+1，并执行步骤(4c)；
(5)获取间歇采样转发干扰识别结果：
(5a)将测试数据集T＝{ty_Sx,ty_ISDJx,ty_ISCJx,ty_ISRJx}中的4×(L-V)个特征数据去除真实类别标签
Figure FDA0002744077750000041
得到用于分类识别的数据，并将其输入训练好的随机森林分类模型进行测试，得到对应的4×(L-V)个预测类别标签
Figure FDA0002744077750000042
(5b)将4×(L-V)个预测类别标签
Figure FDA0002744077750000043
与真实类别标签
Figure FDA0002744077750000044
进行对比，得到目标信号S(t)的识别准确率/>
Figure FDA0002744077750000045
间歇采样直接转发干扰信号ISDJ(t)的识别准确率/>
Figure FDA0002744077750000046
间歇采样重复转发干扰信号ISCJ(t)的识别准确率/>
Figure FDA0002744077750000047
和间歇采样循环转发干扰信号ISRJ(t)的识别准确率
Figure FDA0002744077750000048
其中RS、RISDJ、RISCJ、RISRJ分别表示S(t)、ISDJ(t)、ISCJ(t)、ISRJ(t)的特征数据中预测类别标签和真实类别标签一致的特征数据的个数。
2.根据权利要求1所述的基于谱图平均时间特征的间歇采样转发干扰识别方法，其特征在于，步骤(2d)中所述的以
Figure FDA0002744077750000049
为真实类别标签分别对tSP_Sl、tSP_ISDJl、tSP_ISCJl、tSP_ISRJl进行标注，其中的真实标签值为任意实数，且同一类型干扰的标签值相同，不同类型干扰的标签值不同。
CN202011160606.3A 2020-10-27 2020-10-27 基于谱图平均时间特征的间歇采样转发干扰识别方法 Active CN112083393B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011160606.3A CN112083393B (zh) 2020-10-27 2020-10-27 基于谱图平均时间特征的间歇采样转发干扰识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011160606.3A CN112083393B (zh) 2020-10-27 2020-10-27 基于谱图平均时间特征的间歇采样转发干扰识别方法

Publications (2)

Publication Number Publication Date
CN112083393A CN112083393A (zh) 2020-12-15
CN112083393B true CN112083393B (zh) 2023-06-30

Family

ID=73729907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011160606.3A Active CN112083393B (zh) 2020-10-27 2020-10-27 基于谱图平均时间特征的间歇采样转发干扰识别方法

Country Status (1)

Country Link
CN (1) CN112083393B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113030880B (zh) * 2021-03-17 2024-05-10 南京航空航天大学 基于admm的间歇采样转发式干扰的信号估计方法
CN113075698B (zh) * 2021-03-24 2023-07-07 重庆大学 卫星导航接收机中欺骗式干扰抑制方法
CN113269135A (zh) * 2021-06-17 2021-08-17 中国空间技术研究院 一种卫星转发器识别模型及其训练方法和使用方法
CN116609732B (zh) * 2023-07-21 2023-12-08 江西红声技术有限公司 反间歇采样转发干扰雷达方法、系统、存储介质及计算机

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN110659656A (zh) * 2018-06-29 2020-01-07 中兴通讯股份有限公司 一种检测干扰的方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014048193A1 (zh) * 2012-09-28 2014-04-03 北京理工大学 一种用于舰艇编队情况下同型雷达同频干扰抑制方法
CN110659656A (zh) * 2018-06-29 2020-01-07 中兴通讯股份有限公司 一种检测干扰的方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于频域稀疏性的间歇采样转发干扰识别;蒋莹;何明浩;郁春来;王冰切;;空军预警学院学报(第05期);全文 *

Also Published As

Publication number Publication date
CN112083393A (zh) 2020-12-15

Similar Documents

Publication Publication Date Title
CN112083393B (zh) 基于谱图平均时间特征的间歇采样转发干扰识别方法
CN106330385B (zh) 一种干扰类型识别方法
CN111399002B (zh) 一种基于两级神经网络的gnss接收机组合干扰分类识别方法
CN106125053B (zh) 脉冲多普勒雷达极化抗干扰方法
CN111722199A (zh) 一种基于卷积神经网络的雷达信号检测方法
US20210399817A1 (en) Wireless channel scenario identification method and system
CN113962151B (zh) 基于深度卷积迁移学习的距离假目标智能识别方法
CN113221631B (zh) 基于卷积神经网络的序列脉冲抗干扰目标检测方法
CN112859012B (zh) 一种基于级联卷积神经网络的雷达欺骗干扰识别方法
CN113472390B (zh) 一种基于深度学习的跳频信号参数估计方法
CN109782261A (zh) 对水下入侵目标引起接收信号变化的无监督学习检测方法
CN109061632A (zh) 一种无人机识别方法
CN114114166A (zh) 一种基于dtm算法的雷达脉冲去交错方法
CN111985349A (zh) 一种雷达接收信号类型分类识别方法及系统
CN113608193A (zh) 一种基于UNet的雷达多目标距离和速度估计方法
Rigling et al. Acf-based classification of phase modulated waveforms
Shapero et al. Identifying agile waveforms with neural networks
CN110515044A (zh) 低慢小雷达空频域四维杂波抑制方法
CN116359854A (zh) 一种基于YOLOv5的防空警戒雷达复合干扰参数估计方法
Gao et al. Specific emitter identification based on instantaneous frequency characteristics
KR102407834B1 (ko) 긴 펄스 레이더 신호의 속성 기반 분류 방법 및 장치
Ma et al. Radar working mode recognition based on hierarchical feature representation and clustering
CN115755020A (zh) 一种机械扫描雷达天线扫描位置的跟踪方法
CN109031212A (zh) 一种雷达跟踪状态下的工作频率优化方法
Wu et al. Using Range-Doppler Spectrum-Based Deep Learning Method to Detect Radar Target in Interrupted Sampling Repeater Jamming

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