CN113156391A - 一种雷达信号多维特征智能分选方法 - Google Patents

一种雷达信号多维特征智能分选方法 Download PDF

Info

Publication number
CN113156391A
CN113156391A CN202110445439.5A CN202110445439A CN113156391A CN 113156391 A CN113156391 A CN 113156391A CN 202110445439 A CN202110445439 A CN 202110445439A CN 113156391 A CN113156391 A CN 113156391A
Authority
CN
China
Prior art keywords
pulse
sorting
radar
neural network
radar signal
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.)
Granted
Application number
CN202110445439.5A
Other languages
English (en)
Other versions
CN113156391B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110445439.5A priority Critical patent/CN113156391B/zh
Publication of CN113156391A publication Critical patent/CN113156391A/zh
Application granted granted Critical
Publication of CN113156391B publication Critical patent/CN113156391B/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
    • G01S7/417—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 involving the use of neural networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (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

本发明公开一种雷达信号多维特征智能分选方法，应用于雷达信号处理领域，针对现有技术中存在的复杂电磁环境中雷达信号实时分选准确率低的问题，本发明利用空时频域多参数联合聚类方法，实现雷达信号预分选；基于深度卷积神经网络基本架构，结合雷达信号波形调制特点，构建雷达信号脉内调制特征提取深度卷积神经网络，提取雷达信号脉内调制特征，拓展雷达信号分选特征维度，结合雷达信号PRI脉间调制特征，实现信号参数交叠情况下的雷达信号主分选；不仅提升了雷达信号分选准确率，同时实现了雷达信号快速分选。

Description

一种雷达信号多维特征智能分选方法
技术领域
本发明属于雷达信号处理领域，特别涉及一种电子对抗中雷达辐射源信号分选技术。
背景技术
随着现代电子科学技术的迅猛发展，无线电电子技术在军事战场中广泛应用，电子对抗已成为现代战争中的一个重要的高科技分战场。雷达信号分选技术是电子对抗领域一项至关重要的技术，这一过程具体是指利用同一部雷达信号参数的关联性和不同雷达信号参数的差异性，从雷达侦察接收机接收到的多个辐射源脉冲相互交错信号中分离出每个辐射源脉冲序列，它是进行雷达信号特征提取、识别以及威胁评估的前提和基础。
传统信号预分选技术主要利用雷达信号脉冲描述字(PDW)，结合容差进行简单分类。当前主要的单参数(PRI)分选算法有扩展关联法、累积差值直方图法(CDIF)、序列差值直方图法(SDIF)以及PRI变换法等。然而，随着信息化装备平台的广泛运用，现代战场电磁环境中信号密度增大，波形设计复杂多变，工作频段宽且有部分重叠，信号密集而且脉间参数交叠越来越严重，传统的预分选技术与单参数分选算法难以应对日益复杂的电子对抗环境，利用多参数实现对雷达信号的分选方法陆续被提出。
文献“X.Wei,Y.Pan,J.Han and H.He,Radar signal sorting based onresemblance coefficient of Welch power spectrum.2017 4th InternationalConference on Systems and Informatics(ICSAI),Hangzhou,2017,pp.1120-1125.”提出了一种基于Welch功率谱相似系数的雷达信号分选方法，该方法通过分析雷达信号的Welch功率谱及其能量，提取脉冲序列的功率谱相似系数，形成分选特征参数集，利用支持向量机完成信号主分选。该方法在低信噪比区域对频率调制信号可以获取较高的分选准确率，但由于该方法仅提取单特征参数，无法全面表征雷达信号，且其分选性能受限于截断信号的长度，在低信噪比下对于调相信号的分选效果不是很好。
文献“Z.Zheng,C.Qi and X.Duan,Sorting algorithm for pulse radar basedon wavelet transform.2017 IEEE 2nd Information Technology,Networking,Electronic and Automation Control Conference(ITNEC),Chengdu,2017,pp.1166-1169.”提出了一种基于小波变换的雷达信号分选方法，该方法利用小波变换提取时域信号突变点模量与到达时间，对小波模量进行聚类分析实现雷达信号预分选，有效地抑制子谐波现象，可以通过累积差值直方图算法实现主分选；然而，小波变换的效果受限于小波基函数与分解层数的选取，依赖于人工经验，无法满足复杂电磁环境中雷达信号实时分选要求。
发明内容
为解决上述技术问题，本发明提出了一种雷达信号多维特征智能分选方法，预分选与主分选结合，利用多参数联合分选，引入深度学习框架，在提升分选准确率的同时实现了雷达信号快速分选。
本发明采用的技术方案为：一种雷达信号多维特征智能分选方法，包括：
S1、侦收原始雷达信号及脉内参数；
S2、对步骤S1侦收的脉内参数进行预分选；
S3、构建深度卷积神经网络；
S4、对步骤S2预分选结果中对应的雷达信号进行数据格式标准化处理；
S5、将步骤S4得到的数据作为步骤S3所构建深度卷积神经网络的输入，对深度卷积神经网络进行训练；
S6、采用步骤S5训练完成的深度卷积神经网络对步骤S4预处理后的雷达信号进行脉内调制特征提取；
S7、针对步骤S6中具有相同脉内调制特征的雷达信号，利用SDIF算法分析PRI脉间调制特征，提取来自不同雷达的脉冲序列，实现脉冲去交错，完成雷达信号主分选。
步骤S2具体采用K-means算法对步骤S1侦收的脉内参数进行聚类，得到若干簇。
步骤S3的深度卷积神经网络结构，包括：5个卷积层、5个池化层、2个全连接层；5个卷积层的激活函数为Relu激活函数。
所述深度卷积神经网络采用softmax进行多目标分类，输出为相应的调制类型。
步骤S4具体为：对步骤S2得到的每一簇雷达信号，进行数据格式化标准处理，包括：对每一簇雷达信号进行下采样以获取基带信号，对基带信号进行零填充以统一信号长度，然后对零填充后的基带信号进行快速傅里叶变换，并提取基带信号的实部、虚部与频谱信息共同构成标准格式的数据。
步骤S5所述训练具体为：将每一簇雷达信号经步骤S4处理后得到标准格式的数据输入步骤S3构建的深度卷积神经网络进行前向传播，并计算代价函数值；使用基于梯度下降的后向传播算法对深度卷积神经网络参数进行更新；迭代进行前后向传播，直至代价函数收敛，从而得到训练完成的深度卷积神经网络。
步骤S6具体为：将与步骤S4处理后的序列数据具有相同格式的信号作为测试样本，将测试样本输入步骤S5中训练完成的深度卷积神经网络进行前向传播，得到测试样本属于各类别的后验概率，比较各类后验概率大小，取最大者所对应的类别作为最终预测结果。
步骤S7具体为：对信号的脉冲TOA差值作统计直方图，并计算对应的检测门限，将所有超过门限的值看作是可能的PRI值进行序列搜索，进而完成PRI值的确定以及同一雷达脉冲序列的提取。
计算检测门限的函数表达式为：
Figure BDA0003036689750000031
其中，Np是脉冲总数，N是直方图上脉冲间隔的总刻度值，c是差值直方图的级数，X和K微常数。
本发明的有益效果：本发明的一种雷达信号多维特征智能分选方法，利用截获的侦察参数与雷达信号，结合雷达信号空、时、频域特点，建立基于多维调制域特征的雷达信号总体分选模型。利用空、时、频域多参数联合聚类方法，实现雷达信号预分选；基于深度卷积神经网络基本架构，结合雷达信号波形调制特点，构建雷达信号脉内调制特征提取深度卷积神经网络，提取雷达信号脉内调制特征，拓展雷达信号分选特征维度，分析雷达信号PRI脉间调制特征，实现信号参数交叠情况下的雷达信号分选。本发明方法具有灵活、准确、高效和泛化能力强的优点。
附图说明
图1为本发明实施例提供的雷达信号分选方案流程图；
图2为本发明实施例提供的预分选算法流程图；
图3为本发明实施例提供的脉内调制特征提取网络结构示意图；
图4为本发明实施例提供的PRI脉间调制特征提取流程图；
图5为本发明实施例提供的雷达信号分选结果图。
具体实施方式
为便于本领域技术人员理解本发明的技术内容，下面结合附图对本发明内容进一步阐释。
如图1所示为本发明的方案流程图，本发明的技术方案为：一种雷达信号多维特征智能分选方法，包括：
S1、侦收原始雷达信号及脉内参数，原始雷达信号即时域信号，脉内参数具体包括载频(RF)、到达方向角(DOA)、脉宽(PW)、到达时间(TOA)等。
步骤S1侦收的原始雷达信号(脉冲流)在部分信号参数上有交叠，在这种情况下难以直接进行信号主分选。因此，利用K-means算法对空(DOA)、时(PW)、频(RF)域这三种脉内参数进行聚类分析，对脉冲信号形成若干聚类，稀释了脉冲流密度，在后续步骤中再针对不同的聚类，利用经预处理后的信号和TOA参数，分别通过CNN和SDIF算法进行主分选。主分选的任务是：得到各部雷达的脉内调制类型，以及估计各部雷达的脉冲重复间隔(PRI)、提取来自不同雷达的脉冲序列。
S2、利用K-means算法对步骤S1侦收得到的脉内参数进行预分选，具体为：对于给定的样本集，按照样本之间的距离大小，把N个对象分为k个簇，使得簇内相似度高，簇间相似度低；假设簇划分为(C1,C2,……,Ck)，则我们的目标是最小化平方误差
Figure BDA0003036689750000041
其中x表示由步骤S1中侦收的脉内参数共同组成的样本点，
Figure BDA0003036689750000042
表示向量2-范数，
Figure BDA0003036689750000043
是簇Ci的均值向量(质心)。
本发明实施例中，采用RF、DOA、PW组成的脉冲描述字pdwi＝(rfi,doai,pwi)作为待分选数据点，聚类流程如图2所示；具体包括以下步骤：
S21、随机选取k个聚类中心。
S22、计算待分选数据点到各中心点的欧式距离：
Figure BDA0003036689750000044
式中，dj(pdwi,μj)表示第i个数据点到第j个中心点的欧氏距离，其中，j＝1,2,...,k。根据最小距离原则，将全部样本分配到k个类别中。
S23、根据式(2)重新计算聚类中心：
Figure BDA0003036689750000045
式中，μi为第i类的中心点；k为类别数；pdwi为第i类数据点，ni为第i类数据点的个数。
S24、重复步骤S22和S23，直至到达新的聚类中心与前一聚类中心相等时，结束。
S3、构建深度卷积神经网络，具体为：采用包括5个卷积层、5个池化层、2个全连接层的网络结构；将步骤S4数据预处理后得到的序列数据作为深度卷积神经网络的输入；Relu激活函数应用于所有卷积层，利用softmax进行多目标分类，输出为相应的调制类型；另外，由于在全连接的层中有大量可训练的参数，因此本发明实施例还使用随机失活(Dropout)技术来随机减少连接，避免过度拟合，以提高网络的泛化能力。
如图3所示为本发明实施例的具体网络结构，其中N为待分类目标种类数，卷积层在图3中的表示方法为“卷积层(特征图数目)@(卷积核大小)”；池化层在图3中的表示方法为“池化层@(池化窗口大小)”；全连接层在图3中的表示方法为“全连接层(神经元数目)”。
在本实施例中，考虑到雷达信号的特征更多反映在时序关系上，将卷积核设计为一维卷积核的形式，从而保留卷积后雷达信号在各个维度的时序关系不被破坏，因此各个卷积层对应的卷积核大小分别为：1*60，1*50，1*40，1*30，1*20；为减小卷积层输出后训练参数的数量，利用最大池化层来减少可训练参数，减轻网络训练负担；第一层最大池化层的步长设置为4，其余最大池化层的步长均设置为2。
S4、步骤S2对脉内参数进行聚类后会形成若干簇，针对每一簇对应的雷达信号，进行数据格式标准化处理，具体为：对原始雷达信号进行下采样以获取基带信号；对基带信号进行零填充以统一信号长度，然后对零填充后的基带信号进行快速傅里叶变换，并提取基带信号的实部、虚部与频谱信息共同构成后续网络输入的数据格式；具体包括以下步骤：
S41、以fs＝1.2GHz的采样频率对雷达信号进行下采样，以获得基带信号，则采样点数Ns的计算公式为：
Ns＝τfs (3)
其中，τ为雷达信号的脉冲宽度。
S42、对具有较小脉冲宽度的信号执行零填充操作，以满足网络输入的固定大小；零填充的具体操作可以表示为concatenate(Si,Oi)，其中concatenate(·)代表两个向量的串联，Si是基带信号的实部或虚部，Oi是长度为Li＝(τmax-τi)fs的一维零向量，i代表输入信号的阶数。
S43、对基带信号进行快速傅里叶变换，以克服低信噪比情况下特征提取的困难；将频谱信息与步骤S42零填充得到的向量进行整合，作为卷积神经网络的输入，故网络的输入数据格式为3×12000。
S5、训练深度神经网络，具体为：将步骤S4数据预处理后得到的序列数据输入步骤S3构建的深度卷积神经网络进行前向传播，计算代价函数值；使用基于梯度下降的后向传播算法对深度卷积神经网络参数进行更新；迭代进行前后向传播，直至代价函数收敛；每次迭代对应一个簇的数据。具体包括以下步骤：
S51、前向传播，若第
Figure BDA00030366897500000616
层为卷积层，以
Figure BDA0003036689750000061
为第
Figure BDA00030366897500000617
层的第j个特征映射，则
Figure BDA0003036689750000062
其中，Mj为第
Figure BDA00030366897500000618
层中与之相连的特征映射集合，
Figure BDA0003036689750000063
表示连接第i个输入特征图和第j个输出特征图的卷积核，
Figure BDA0003036689750000064
表示偏置项，f(·)表示非线性激活函数，符号“*”表示二维离散卷积；
若第
Figure BDA00030366897500000619
层为池化层，则其输出为
Figure BDA0003036689750000065
其中，down(·)表示下采样函数，
Figure BDA0003036689750000066
分别表示乘性和加性偏置项；
若第
Figure BDA00030366897500000620
层为全连接层，则其输出为
Figure BDA00030366897500000621
其中，
Figure BDA00030366897500000615
表示第
Figure BDA00030366897500000614
层的特征图；
Figure BDA00030366897500000613
表示第
Figure BDA00030366897500000612
层权重，b(l)为该层偏置项；
若第L层为输出层，则当前样本属于第i类的后验概率为
Figure BDA0003036689750000067
其中，zL表示输出层的输入，C表示总类别数。
定义
Figure BDA0003036689750000068
为第
Figure BDA00030366897500000623
层的误差灵敏项，
Figure BDA00030366897500000622
为非线性激活函数的输入，设第L层为输出层，则其灵敏项为：
Figure BDA0003036689750000069
其中，
Figure BDA00030366897500000610
表示点乘。
S52、计算代价函数值，图3所示实施例中以均方误差作为代价函数，则第n个样本对应的代价函数为：
Figure BDA0003036689750000071
其中，
Figure BDA0003036689750000072
表示第n个样本对应的类别标签的第k维，
Figure BDA0003036689750000073
表示第n个样本对应的网络输出的第k维，C表示总类别数。
S53、利用基于梯度下降的后向传播算法对网络参数进行更新，第
Figure BDA00030366897500000711
层的更新变化量为：
Figure BDA0003036689750000074
Figure BDA0003036689750000075
Figure BDA0003036689750000076
其中，η为学习率，
Figure BDA0003036689750000077
表示与
Figure BDA0003036689750000078
中(u,v)处的值相对应的
Figure BDA0003036689750000079
中与
Figure BDA00030366897500000710
进行卷积运算的区域。
S6、利用步骤S5训练完成的网络对雷达信号进行脉内调制特征提取，具体为：产生与步骤S4处理后的序列数据具有相同格式的信号作为测试样本，将测试样本输入步骤S5中训练完成的网络进行前向传播，得到测试样本属于各类别的后验概率，比较各类后验概率大小，取最大者所对应的类别作为最终预测结果，实现雷达信号脉内调制特征聚类分析。
S7、针对步骤S6中具有相同脉内调制特征的雷达信号，利用SDIF算法分析雷达信号PRI脉间调制特征，至此完成信号主分选，具体为：对信号的脉冲TOA差值作统计直方图，并计算对应的检测门限，将超过门限的值看作是可能的PRI值进行序列检索，进而完成PRI值的确定以及同一雷达脉冲序列的提取。序列检索是SDIF算法的关键步骤，一般将窗口大小设置为5，如果在脉冲序列中有5个脉冲按照此PRI值被选中，则认为成功分离出该PRI值对应的脉冲串，定义为序列检索成功；如果每次选中的脉冲数少于5个，则认为该脉冲序列中没有与该PRI值相关联的脉冲，定义为序列检索不成功。
如图5所示，具体包括以下步骤：
S71、对待处理的Np个信号脉冲，计算相邻两个脉冲TOA的差值，按差值大小进行统计并排序，构成脉冲TOA的第一级差值直方图；
S72、任何直方图算法的关键是门限函数，因为它影响检测速度和检测概率，最佳检测门限最接近直方图的峰值，可以使检测速度最快，检测概率最高，而指数函数特性与之最相符。由于直方图的序数与两脉冲之间的间隔相一致，在观察时间一定时，观察的脉冲间隔越大，观察出现的脉冲数量越少，因而门限与序数(τ)成反比，于是门限函数可表示为：
Figure BDA0003036689750000081
其中，Np是观察的脉冲总数，X是小于1的常数。
若观察的脉冲总数足够多，且同时有多个辐射源，则可认为相邻两脉冲的间隔是随机事件，也就可认为脉冲是随机的泊松流。若把一定的观察时间T分为若干个子间隔，则在时间间隔(τ＝t2-t1)内有K个子间隔出现的概率为：
Figure BDA0003036689750000082
其中，λ＝n/T(假定λ为常量)，n是子间隔总数，T是观察时间片段，λ是泊松流参数，表示事件在某个时间间隔里发生的平均数。
相邻两脉冲的时间间隔为τ(即K＝0的情况)的概率为：
P0(τ)＝e-λτ (15)
第一级差值直方图呈(15)式的形式，因为直方图是随机事件概率分布函数的近似值，所以较高级差值直方图也呈指数分布形式。第c级差值直方图的脉冲数量是(Np-c)，观察的时间间隔与直方图总序数成正比，泊松流参数λ＝1/KN，这里K是小于1的正常数。因此，最佳检测门限函数为：
Figure BDA0003036689750000083
其中，τ是直方图的序数，Np是脉冲总个数，N是直方图上脉冲间隔的总刻度值，c是当前差值直方图的级数，常数X和K由经验值或实验确定。实际上，若丢失脉冲发生的概率分布也为泊松分布，处理过程不变，则门限代表直方图的峰值。因此，常数X由假定的脉冲丢失概率决定。
S73、根据式(15)计算检测门限，若只有一个差值大于该门限，则将它默认为可能的PRI值并进行序列检索，转至步骤S75；若有多个差值大于该门限，则不进行序列检索，转至步骤S74；
S74、如果脉冲序列中包含的脉冲数大于5(即设定的窗口大小)，则计算下一级差值直方图，根据式(15)计算检测门限，将大于门限的那些差值作为可能的PRI值，对它们进行排序和序列检索，转至步骤S75；如果脉冲序列中包含的脉冲数小于5，转至步骤S76；
S75、如果序列检索不成功，则转至步骤S74；如果序列检索成功，则在脉冲序列中将此脉冲剔除，并对剩下的脉冲序列从第一级开始形成差值直方图，并重新进行步骤S73；
S76、进行参差分析与处理，算法结束。
最后本发明还包括：
(1)对步骤S5训练完成的网络进行脉内调制类型识别性能测试，具体为：将测试样本输入S5中训练完成的网络进行前向传播，得到测试样本属于各类别的后验概率，比较各类后验概率大小，取最大者所对应的类别作为最终预测结果，验证所提网络对雷达信号脉内调制类型的识别性能。
(2)对所提雷达信号多维特征智能分选方法的分选性能测试，具体为：利用测试样本脉内参数进行多域参数联合聚类实现雷达信号预分选，根据预分选结果利用深度卷积神经网络分析脉内调制特征，利用SDIF算法分析PRI脉间调制特征，完成雷达信号分选。分选成功率定义为正确分选脉冲个数与待分选脉冲个数的比值。
表1为实施例中所用网络训练与测试使用的调制信号数据集，包括不同调制类型和不同信号带宽的27类基带信号，具体包括非线性调频(NLFM)，线性调频(LFM)，单脉冲(MP)，脉内频率捷变(FA)，Costas频率编码，Barker相位编码，Costas/Barker调频调相复合编码，幅度调制(AM)和相位调制(PM)信号。卷积神经网络由6480个训练样本训练，训练样本SNR＝12dB；并由1620个测试样本进行测试，当SNR＝12dB时，得到调制类型识别结果如图5所示，本发明总体识别率达到95.93％。
表2为实施例中完成雷达信号多维特征智能分选方法测试使用的雷达信号数据集，包括不同调制类型、带宽介于100～880MHz的23类雷达辐射源，共639个脉冲；当SNR＝12dB时，正确分选612个脉冲，本发明实施例所提方法的总体分选成功率可达95.77％。
表1调制信号数据集参数对应表
Figure BDA0003036689750000101
表2雷达信号多维特征智能分选方法测试数据集
调制类型 标签 数目 PRI/us 脉宽/us 带宽/MHz
Barker 1 20 1250 5.2 100
Costas 2 32 700/900/1200 2.4 420
LFM 3 22 1100±5％ 2.4 420
Barker 4 41 1150/1201/1050/1601 6.0 620
Costas 5 26 1530±5％ 3.4 780
FA 6 20 1801 6 620
AM 7 26 1350/1750 4.2 550
MP 8 36 1620/750/1250 6.5 560
FA 9 24 1301±5％ 7.4 750
LFM 10 20 1380 5.5 820
Costas/Barker 11 45 950/1601/1220 7.5 630
NLFM 12 25 1560 7.2 530
FM 13 25 1501±5％ 5.8 280
FA 14 32 1150/1800/1070/1350 2.7 770
Costas/Barker 15 36 250/1250/2250 2.7 770
NLFM 16 24 1216/954.5/2460 5.5 750
LFM 17 24 1820/1300/1460 6.1 280
MP 18 20 1750 3.0 280
Costas/Barker 19 24 1300/1500/1700 4.3 880
FA 20 25 1720±5％ 5.5 750
NLFM 21 36 1160/1500/1080 8.1 750
FA 22 32 265.7/1860/812.8/1589 6.5 420
Costas/Barker 23 24 1340±5％ 4.8 840
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。

Claims (9)

1.一种雷达信号多维特征智能分选方法，其特征在于，包括：
S1、侦收原始雷达信号及脉内参数；
S2、对步骤S1侦收的脉内参数进行预分选；
S3、构建深度卷积神经网络；
S4、对步骤S2预分选结果中对应的雷达信号进行数据格式标准化处理；
S5、将步骤S4得到的数据作为步骤S3所构建深度卷积神经网络的输入，对深度卷积神经网络进行训练；
S6、采用步骤S5训练完成的深度卷积神经网络对步骤S4预处理后的雷达信号进行脉内调制特征提取；
S7、针对步骤S6中具有相同脉内调制特征的雷达信号，利用SDIF算法分析PRI脉间调制特征，提取来自不同雷达的脉冲序列，实现脉冲去交错，完成雷达信号主分选。
2.根据权利要求1所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S2具体采用K-means算法对步骤S1侦收的脉内参数进行聚类，得到若干簇。
3.根据权利要求1所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S3的深度卷积神经网络结构，包括：5个卷积层、5个池化层、2个全连接层；5个卷积层的激活函数为Relu激活函数。
4.根据权利要求3所述的一种雷达信号多维特征智能分选方法，其特征在于，所述深度卷积神经网络采用softmax进行多目标分类，输出为相应的调制类型。
5.根据权利要求2所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S4具体为：对步骤S2得到的每一簇雷达信号，进行数据格式化标准处理，包括：对每一簇雷达信号进行下采样以获取基带信号，对基带信号进行零填充以统一信号长度，然后对零填充后的基带信号进行快速傅里叶变换，并提取基带信号的实部、虚部与频谱信息共同构成标准格式的数据。
6.根据权利要求5所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S5所述训练具体为：将每一簇雷达信号经步骤S4处理后得到标准格式的数据输入步骤S3构建的深度卷积神经网络进行前向传播，并计算代价函数值；使用基于梯度下降的后向传播算法对深度卷积神经网络参数进行更新；迭代进行前后向传播，直至代价函数收敛，从而得到训练完成的深度卷积神经网络。
7.根据权利要求6所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S6具体为：将与步骤S4处理后的序列数据具有相同格式的信号作为测试样本，将测试样本输入步骤S5中训练完成的深度卷积神经网络进行前向传播，得到测试样本属于各类别的后验概率，比较各类后验概率大小，取最大者所对应的类别作为最终预测结果。
8.根据权利要求7所述的一种雷达信号多维特征智能分选方法，其特征在于，步骤S7具体为：对信号的脉冲TOA差值作统计直方图，并计算对应的检测门限，将所有超过门限的值看作是可能的PRI值进行序列搜索，进而完成PRI值的确定以及同一雷达脉冲序列的提取。
9.根据权利要求8所述的一种雷达信号多维特征智能分选方法，其特征在于，计算检测门限的函数表达式为：
Figure FDA0003036689740000021
其中，Np是脉冲总数，N是直方图上脉冲间隔的总刻度值，c是差值直方图的级数，X和K为常数。
CN202110445439.5A 2021-04-25 2021-04-25 一种雷达信号多维特征智能分选方法 Active CN113156391B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110445439.5A CN113156391B (zh) 2021-04-25 2021-04-25 一种雷达信号多维特征智能分选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110445439.5A CN113156391B (zh) 2021-04-25 2021-04-25 一种雷达信号多维特征智能分选方法

Publications (2)

Publication Number Publication Date
CN113156391A true CN113156391A (zh) 2021-07-23
CN113156391B CN113156391B (zh) 2022-08-05

Family

ID=76870675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110445439.5A Active CN113156391B (zh) 2021-04-25 2021-04-25 一种雷达信号多维特征智能分选方法

Country Status (1)

Country Link
CN (1) CN113156391B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113721219A (zh) * 2021-10-08 2021-11-30 中国电子科技集团公司第三十八研究所 一种基于多参数聚类的雷达信号分选方法和系统
CN113962261A (zh) * 2021-10-21 2022-01-21 中国人民解放军空军航空大学 一种用于雷达信号分选的深度网络模型
CN114595732A (zh) * 2022-05-10 2022-06-07 西安晟昕科技发展有限公司 基于深度聚类的雷达辐射源分选方法
CN114624645A (zh) * 2022-03-10 2022-06-14 扬州宇安电子科技有限公司 一种基于微型天线阵列的小型旋翼无人机雷达侦察系统
CN115062658A (zh) * 2022-06-14 2022-09-16 电子科技大学 基于自适应门限网络的交叠雷达信号调制类型识别方法
CN115166678A (zh) * 2022-09-06 2022-10-11 中国人民解放军海军工程大学 复杂场景雷达辐射源脉冲信号分选与认知对抗系统
CN115390037A (zh) * 2022-09-06 2022-11-25 中国人民解放军海军工程大学 多类别、未知雷达辐射源脉冲信号分选系统
CN116559817A (zh) * 2023-07-04 2023-08-08 西安晟昕科技股份有限公司 一种用于复杂电磁环境的雷达信号分选方法
CN117129947A (zh) * 2023-10-26 2023-11-28 成都金支点科技有限公司 一种基于miniunet的平面变换法雷达信号识别方法
CN117554919A (zh) * 2024-01-11 2024-02-13 成都金支点科技有限公司 一种基于双向lstm网络的雷达信号分选搜索方法
CN115062658B (zh) * 2022-06-14 2024-05-28 电子科技大学 基于自适应门限网络的交叠雷达信号调制类型识别方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106405518A (zh) * 2016-12-07 2017-02-15 中国船舶重工集团公司第七二四研究所 一种复杂体制雷达信号等级关联聚类分选方法
CN108898157A (zh) * 2018-05-28 2018-11-27 浙江理工大学 基于卷积神经网络的数值型数据的雷达图表示的分类方法
CN109117747A (zh) * 2018-07-23 2019-01-01 电子科技大学 基于循环平稳与深度卷积神经网络的雷达信号分类方法
CN110109060A (zh) * 2019-03-27 2019-08-09 西安电子科技大学 一种基于深度学习网络的雷达辐射源信号分选方法及系统
US20200025877A1 (en) * 2018-07-18 2020-01-23 Qualcomm Incorporated Object verification using radar images
CN110764063A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 一种基于sdif与pri变换法结合的雷达信号分选方法
CN111239692A (zh) * 2020-01-10 2020-06-05 西安电子科技大学 基于深度学习的pri联合脉内信息辐射源信号识别方法
CN111693945A (zh) * 2020-06-19 2020-09-22 清华大学 脉间调制方式识别方法、系统、计算机设备及存储介质
CN111722188A (zh) * 2020-05-18 2020-09-29 中国人民解放军63892部队 基于stft预分选的pri变换雷达信号分选方法
CN112528774A (zh) * 2020-11-27 2021-03-19 中国运载火箭技术研究院 一种复杂电磁环境下未知雷达信号智能分选系统及方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106405518A (zh) * 2016-12-07 2017-02-15 中国船舶重工集团公司第七二四研究所 一种复杂体制雷达信号等级关联聚类分选方法
CN108898157A (zh) * 2018-05-28 2018-11-27 浙江理工大学 基于卷积神经网络的数值型数据的雷达图表示的分类方法
US20200025877A1 (en) * 2018-07-18 2020-01-23 Qualcomm Incorporated Object verification using radar images
CN109117747A (zh) * 2018-07-23 2019-01-01 电子科技大学 基于循环平稳与深度卷积神经网络的雷达信号分类方法
CN110109060A (zh) * 2019-03-27 2019-08-09 西安电子科技大学 一种基于深度学习网络的雷达辐射源信号分选方法及系统
CN110764063A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 一种基于sdif与pri变换法结合的雷达信号分选方法
CN111239692A (zh) * 2020-01-10 2020-06-05 西安电子科技大学 基于深度学习的pri联合脉内信息辐射源信号识别方法
CN111722188A (zh) * 2020-05-18 2020-09-29 中国人民解放军63892部队 基于stft预分选的pri变换雷达信号分选方法
CN111693945A (zh) * 2020-06-19 2020-09-22 清华大学 脉间调制方式识别方法、系统、计算机设备及存储介质
CN112528774A (zh) * 2020-11-27 2021-03-19 中国运载火箭技术研究院 一种复杂电磁环境下未知雷达信号智能分选系统及方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
QIANG LIU等: ""Deep learning and recognition of radar jamming based on CNN"", 《2019 12TH INTERNATIONAL SYMPOSIUM ON COMPUTATIONAL INTELLIGENCE AND DESIGN》 *
史雨璇: "" 认知侦察中的信号分选算法研究"", 《万方学位论文数据库》 *
杨晓红: ""LPI雷达信号智能检测与识别技术"", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *
王益明: ""复杂环境下雷达信号分选算法研究"", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113721219A (zh) * 2021-10-08 2021-11-30 中国电子科技集团公司第三十八研究所 一种基于多参数聚类的雷达信号分选方法和系统
CN113962261A (zh) * 2021-10-21 2022-01-21 中国人民解放军空军航空大学 一种用于雷达信号分选的深度网络模型
CN113962261B (zh) * 2021-10-21 2024-05-14 中国人民解放军空军航空大学 一种用于雷达信号分选的深度网络模型构建方法
CN114624645A (zh) * 2022-03-10 2022-06-14 扬州宇安电子科技有限公司 一种基于微型天线阵列的小型旋翼无人机雷达侦察系统
CN114624645B (zh) * 2022-03-10 2022-09-30 扬州宇安电子科技有限公司 一种基于微型天线阵列的小型旋翼无人机雷达侦察系统
CN114595732A (zh) * 2022-05-10 2022-06-07 西安晟昕科技发展有限公司 基于深度聚类的雷达辐射源分选方法
CN115062658A (zh) * 2022-06-14 2022-09-16 电子科技大学 基于自适应门限网络的交叠雷达信号调制类型识别方法
CN115062658B (zh) * 2022-06-14 2024-05-28 电子科技大学 基于自适应门限网络的交叠雷达信号调制类型识别方法
CN115166678B (zh) * 2022-09-06 2022-12-20 中国人民解放军海军工程大学 复杂场景雷达辐射源脉冲信号分选与认知对抗系统
CN115390037A (zh) * 2022-09-06 2022-11-25 中国人民解放军海军工程大学 多类别、未知雷达辐射源脉冲信号分选系统
CN115166678A (zh) * 2022-09-06 2022-10-11 中国人民解放军海军工程大学 复杂场景雷达辐射源脉冲信号分选与认知对抗系统
CN116559817A (zh) * 2023-07-04 2023-08-08 西安晟昕科技股份有限公司 一种用于复杂电磁环境的雷达信号分选方法
CN116559817B (zh) * 2023-07-04 2023-09-26 西安晟昕科技股份有限公司 一种用于复杂电磁环境的雷达信号分选方法
CN117129947A (zh) * 2023-10-26 2023-11-28 成都金支点科技有限公司 一种基于miniunet的平面变换法雷达信号识别方法
CN117129947B (zh) * 2023-10-26 2023-12-26 成都金支点科技有限公司 一种基于miniunet的平面变换法雷达信号识别方法
CN117554919A (zh) * 2024-01-11 2024-02-13 成都金支点科技有限公司 一种基于双向lstm网络的雷达信号分选搜索方法
CN117554919B (zh) * 2024-01-11 2024-03-29 成都金支点科技有限公司 一种基于双向lstm网络的雷达信号分选搜索方法

Also Published As

Publication number Publication date
CN113156391B (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
CN113156391B (zh) 一种雷达信号多维特征智能分选方法
CN108416290B (zh) 基于残差深度学习的雷达信号特征提取方法
Petrov et al. Radar emitter signals recognition and classification with feedforward networks
Liu Multi-feature fusion for specific emitter identification via deep ensemble learning
CN107301381A (zh) 基于深度学习和多任务学习策略的雷达辐射源识别方法
CN109379153B (zh) 一种频谱感知方法
Kauppi et al. Hierarchical classification of dynamically varying radar pulse repetition interval modulation patterns
CN106483514B (zh) 一种基于eemd和支持向量机的飞机运动模式识别方法
CN108051781A (zh) 一种基于dbn模型的雷达信号工作模式识别方法
Li Research on radar signal recognition based on automatic machine learning
Cain et al. Convolutional neural networks for radar emitter classification
CN107830996B (zh) 一种飞行器舵面系统故障诊断方法
CN112597820A (zh) 一种基于雷达信号分选的目标聚类方法
CN112684427A (zh) 基于串行二次强化训练的雷达目标识别方法
CN113050797A (zh) 一种通过毫米波雷达实现手势识别的方法
Zhao et al. A novel aggregated multipath extreme gradient boosting approach for radar emitter classification
CN114636975A (zh) 基于谱图融合和注意力机制的lpi雷达信号识别方法
CN111208483B (zh) 基于贝叶斯支持向量数据描述的雷达库外目标识别方法
Yar et al. A complete framework of radar pulse detection and modulation classification for cognitive EW
Nuhoglu et al. Image segmentation for radar signal deinterleaving using deep learning
CN114371474A (zh) 基于卷积降噪自编码器的雷达信号智能分选方法及系统
Zhang et al. A hybrid method for classification and identification of emitter signals
CN112990125B (zh) 一种辐射源雷达属于目标平台的判别方法
Willson Radar classification using a neural network
Li et al. Sea/land clutter recognition for over-the-horizon radar via deep CNN

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