CN109239670B - 基于结构嵌入和深度神经网络的雷达hrrp识别方法 - Google Patents

基于结构嵌入和深度神经网络的雷达hrrp识别方法 Download PDF

Info

Publication number
CN109239670B
CN109239670B CN201810998871.5A CN201810998871A CN109239670B CN 109239670 B CN109239670 B CN 109239670B CN 201810998871 A CN201810998871 A CN 201810998871A CN 109239670 B CN109239670 B CN 109239670B
Authority
CN
China
Prior art keywords
neural network
quantization
sample
data
layer
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
CN201810998871.5A
Other languages
English (en)
Other versions
CN109239670A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201810998871.5A priority Critical patent/CN109239670B/zh
Publication of CN109239670A publication Critical patent/CN109239670A/zh
Application granted granted Critical
Publication of CN109239670B publication Critical patent/CN109239670B/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于结构嵌入和深度神经网络的雷达HRRP识别方法，首先提取数据的时域特征并对提取出的时域特征进行切分和非均匀量化编码得到其局部结构的编码，然后利用局部结构和其周围数个局部结构的关系得到其编码之间的共现矩阵，进而通过该共现矩阵得到数据的结构嵌入特征，然后再将提取出的嵌入特征送入由全连接层和基于LSTM的双向循环神经网络构成的深度神经网络进行训练，最后利用softmax层对HRRP进行识别并得到识别结果。

Description

基于结构嵌入和深度神经网络的雷达HRRP识别方法
技术领域
本发明属于雷达目标识别领域，具体的讲，涉及一种基于高分辨距离像结构嵌入特征和深度神经网络的雷达高分辨距离像识别方法。
背景技术
雷达高分辨距离像(HRRP)是用宽带雷达信号获取的目标散射点子回波在雷达视线方向投影的矢量和，它包含了目标的尺寸和散射点分布等重要结构信息，并且易于获取和快速处理。如今已经成为了雷达自动目标识别(RATR)领域的一个重要的研究方向，因此利用HRRP进行目标识别已经成为当前研究的热点。针对雷达HRRP识别的一些经典的模型如基于HRRP的一阶统计特性的模板匹配分类器(MCC-TMM)和基于融合一阶和二阶HRRP特性的自适应高斯分类器(AGC)，基于主成分分析(PCA)子空间投影模型等都是将HRRP样本作为一个整体来考虑的，未考虑到样本内部结构的相关性。鉴于这一特性，后来又发展出了基于HRRP样本内部结构特性的隐马尔可夫模型(HMM)和基于注意力机制的循环神经网络模型，它们均利用了HRRP背后所隐含的固定空间结构特性来进行识别。目前已有不少文献应用基于HRRP样本内部结构特性的方法对雷达HRRP的时域特征以及雷达整体HRRP序列进行实验，并且取得了很好的识别效果，然而在实际应用中，还存在以下几个主要的问题：(1)HMM在模型的建立过程中，假设信号符合一阶马尔可夫性，即当前时间点上的信号仅和前一个时间点信号相关，应用该方法对雷达HRRP的时域特征、谱图特征识别时，雷达HRRP特征的当前局部特征结构仅依赖于上一个局部特征结构，未能进一步发掘样本各局部特征结构之间隐含的整体物理结构相关性，尚存在很大改进空间；(2)基于注意力机制的循环神经网络模型为单向的循环神经网络模型，只能朝着一个方向建立依赖性，未能很好的利用HRRP整体的结构信息；(3)基于注意力机制的循环神经网络模型中所采用的循环神经网络过于简单，未使用长短时记忆(LSTM)单元，其循环神经网络的输出会过度的依赖于局部的结构，而很难顾全全局性的信息。由于这些问题的存在，大大限制了它们在实际工程中的应用，因此解决以上存在的问题在实际工程应用中是十分有意义的。
发明内容
本发明针对上述提到的应用HMM模型和基于注意力机制的循环神经网络模型进行雷达HRRP目标识别中存在的几点不足，提出了一种基于高分辨距离像结构嵌入特征与深度神经网络的雷达高分辨距离像目别方法，首先提取数据的时域特征并对提取出的时域特征进行切分和非均匀量化编码得到其局部结构的编码，然后利用局部结构和其周围数个局部结构的关系得到其编码之间的共现矩阵，进而通过该共现矩阵得到数据的结构嵌入特征，然后再将提取出的嵌入特征送入由全连接层和基于LSTM的双向循环神经网络构成的深度神经网络进行训练，最后利用softmax层对HRRP进行识别并得到识别结果。
为实现上述目的，本发明的技术方案为一种基于结构嵌入和深度神经网络的雷达HRRP识别方法，包括以下步骤；
S1:采集数据集，将Q个雷达采集到的HRRP数据集依据目标的种类进行合并，每种种类的数据分别在不同数据段里选择训练样本和测试样本，在训练集和测试集样本的选取过程中，保证所选取训练集的数据与雷达所成姿态涵盖测试数据集与雷达所成的姿态。
S2:对经过S1所挑选数据集中的每一个样本使用补偿对齐法进行对齐，得到对齐后的数据；
S3:对S2处理后的每一个样本进行能量归一化处理，得到归一化后的数据；
S4：对经过S3处理后的雷达训练目标数据，提取时域特征作为识别特征，并用于训练；
S5:循环切分经过处理后的数据集，设定切分的向量长度为p，使p为2的幂级数，如4；之后对数据集内的每一个HRRP样本，将其移位p-1次并进行连接；
S6:提取样本的嵌入表示特征，将每一个经过S5处理后的HRRP样本按10级进行非均匀量化；设定窗口值，初始化一个空矩阵，并依次统计量化向量在窗长度下共现的次数，对矩阵中的对应元素赋值；最后对矩阵中的每个元素限制其动态范围，根据共现矩阵得到HRRP训练样本的结构嵌入特征；
S7:构建深度神经网络模型，并利用HRRP训练样本的结构嵌入特征训练该模型。该深度神经网络模型由三部分组成：首先是一个全连接网络层，用于嵌入特征降维，然后将降维后的数据输入至基于LSTM(长短时记忆网络)的双向循环神经网络，最后采用softmax层对双向循环神经网络的输出进行识别；
测试阶段：
S8:对由S1采集到的测试数据进行训练阶段的步骤S2、S3、S4、S5和S6操作；
S9:将经过S8所提取的嵌入表示特征送入S7构建的深度神经网络得到测试样本所对应的目标。
进一步的，步骤S1中各类目标训练集和测试集样本数的比例为8:2，将挑选的数据集记作
Figure BDA0001782160770000031
其中Xi表示第i个样本，yk表示第k种目标，一共采集了5种目标，i0表示样本总数。
进一步的，步骤S2具体操作步骤如下：
S2.1:由S1所选的数据集T中的每个样本皆为256维度，原始数据样本表示为：Xi＝[x1,x2,...,x256]；
S2.2:将三个原始样本按照横向粘贴在一起构成一个256×3新样本数据样本，表示为X′i＝abs([Xi,Xi,Xi])＝[x′1,x′2,…,x′768]，其中abs(·)为对其中每个元素的值取绝对值；
S2.3:构建一个256维的序列L＝[l1,l2,…,l256]，其中li＝exp(-k|i-128.5|)，k表示控制因子，k＞0，|·|表示求绝对值操作；
S2.4:求
Figure BDA0001782160770000032
得到d之后，从X'中取出第d+1到第d+256个连续的点作为对齐后的数据X”＝[x'd+1,x'd+2,…,x'd+256]；
S2.5:将T数据集中的每一个样本重复步骤S2.2、S2.3以及S2.4操作得到对齐后的样本数据，记为T1。
进一步的，步骤S3中，归一化公式如下：
Figure BDA0001782160770000041
其中||·||2表示求2范数。
进一步的，步骤S4中，从X″′i内提取每一帧的时域特征作为识别特征,记作
Figure BDA0001782160770000042
其中x(c)表示的是第c个目标的训练样本的时域特征集，其中提取时域特征之后的第i个样本记作
Figure BDA0001782160770000049
C为目标总数.
进一步的，步骤S5具体操作步骤如下：
S5.1:
Figure BDA0001782160770000043
其中i指示第i个目标，k指示的是第k个HRRP；
S5.2:设定切分的向量的长度为p,使得为p为2的幂级数；
S5.3:对数据集内的每一个样本，将其移位p-1次并进行连接，
Figure BDA0001782160770000044
若
Figure BDA0001782160770000045
则
Figure BDA0001782160770000046
进一步的，步骤S6具体操作步骤如下：
S6.1:对经过S5处理后的每一个
Figure BDA0001782160770000047
按10级进行非均匀量化，得到其量化编码特征
Figure BDA0001782160770000048
其中对于每个元素：
若y＜e-3.6，则量化为0
若e-3.6≤y＜e-3.2，则量化为1
若e-3.2≤y＜e-2.8，则量化为2
若e-2.8≤y＜e-2.4，则量化为3
若e-2.4≤y＜e-2.0，则量化为4
若e-2.0≤y＜e-1.6，则量化为5
若e-1.6≤y＜e-1.2，则量化为6
若e-1.2≤y＜e-0.8，则量化为7
若e-0.8≤y＜e-0.4，则量化为8
若e-0.4≤y＜e0，则量化为9
S6.2:设定窗口值K，计算向量与向量之间共同出现的矩阵；
S6.3:首先初始化一个空矩阵，记作V＝{Vij}i＝0,1…9999,j＝0,1…9999，即初始化的空矩阵有10000行和10000列，其中0000代表取出来的列向量[yi,yi+1,yi+2,yi+3]T的量化编码为[0,0,0,0]T，其余的也是一样，矩阵中的每一个元素Vij为量化编码i和量化编码j之间的共现次数；
S6.4:依次统计列向量之间在窗长度下共现的次数，如
Figure BDA0001782160770000051
若窗K＝3(为保证左右对称，K一般取奇数)，则当循环至向量
Figure BDA0001782160770000052
的时候，统计其左侧
Figure BDA0001782160770000053
和其右侧
Figure BDA0001782160770000054
不同向量共同出现的次数，我们遍历所有的量化编码特征，将不同量化编码特征共同出现的总次数对应地记录在空矩阵V中，其中Vij为量化编码特征
Figure BDA0001782160770000055
与量化编码特征
Figure BDA0001782160770000056
共同出现的次数，其中i1、i2、i3和i4分别代表i的千位、百位、十位以及个位数字，j1、j2、j3和j4分别代表j的千位、百位、十位以及个位数字；
S6.4:对矩阵V中的每个元素限制其动态范围，公式如下：V′ij＝10ln(1+Vij)，并将更新后的矩阵记作V'＝{V″′ij}；
S6.5:根据循环切分的结果去V'中查找
Figure BDA0001782160770000057
所对应的嵌入表示，如
Figure BDA0001782160770000058
中的第m列为[0,1,3,8]T，即V'中行标号为第0138所对应的行向量即为该列的嵌入表示，该嵌入表示为10000×1维的，依次查找
Figure BDA0001782160770000059
中每一列所对应的嵌入表示并记录，遍历
Figure BDA00017821607700000510
中所有列之后即可得到
Figure BDA00017821607700000511
的嵌入表示矩阵A，其中A＝[a1(10000×1)，a2(10000×1)，…，a256(10000×1)]，其中ai(10000×1)代表
Figure BDA00017821607700000512
第i列所对应的嵌入表示，ai为10000×1维的向量。
进一步的，步骤S7具体操作步骤如下：
S7.1:构建一层全连接神经网络，设置第一层(输入层)神经元数目为10000，第二层(输出层)神经元数目为300，将经过步骤S6.5处理后的A按列切分依次送入全连接神经网络，得到输出
Z＝[z1＝f(Wija1+b),z2＝f(Wija2+b),…,z256＝f(Wija256+b)]，其中zi代表第i个输出，为300×1维的，Wij代表输入层第i个神经元到输出层第j个神经元的权重，b代表输入层神经元到输出层神经元的偏置，所有的神经元共享偏置；
S7.2:构建基于LSTM(长短时记忆网络)的双向循环神经网络，具体操作步骤如下：
S7.2.1:双向循环神经网络为两个单向其方向相反的循环神经网络(前向网络和反向网络)的结合，在每一个时刻t，网络的输出由这两个方向相反的单向循环神经网络共同决定，设置n_in＝256，n_step＝300，n_nums＝128，n_cls＝5，其中n_in＝256基于LSTM的双向循环神经网络的输入层神经元数目，n_step＝300代表每一个神经元的输入向量，为一个300×1的向量，n_nums＝128代表LSTM隐藏层神经元的数目，n_cls＝5代表输出层神经元数目，即代表所挑选的目标类别数目为5类；
S7.2.2:对双向循环神经网络的权重和偏置进行初始化，输入层和隐层之间的权重记为W1，其维度为shape1＝(n_in,n_nums),隐层和输出层之间的权重记为W2，其维度为shape2＝(n_nums,n_cls)，权重初始化采用服从均值为0，方差为1的高斯分布，输入层和隐层之间的偏置记为bias1，其维度为shape3＝(n_nums)，隐层和输出层之间的偏置记为bias2，其维度为shape4＝(n_cls)，均初始化为常量1，双向的初始状态皆记为state＝0，前向循环神经网络的所有时刻共享前向的权重，偏置和初始，反向也是如此；
S7.2.3:由S7.1可知，全连接神经网络的输出Z作为循环神经网络的输入，则循环神经网络第一层的输出为X1_out＝tanh(Z×W1+bias1)，同时作为隐层的输入，其中tanh(·)是激活函数，根据tensorflow中的计算公式
lstm_fw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)
lstm_bw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)
output,_,_＝tf.contrib.rnn.static_bidirectional_rnn(lstm_fw_cell,lstm_bw_cell,X1_out)
其中，lstm_fw_cell代表前向的循环神经网络，lstm_bw_cell代表反向循环神经网络，forget_bias＝1.0代表LSTM不会忘记任何信息，output,_,_分别代表X1_out经过隐层后的输出，前向循环神经网络和反向循环神经网络被更新后的状态，取出序列中的最后一个X2_out＝output[-1]作为输出层的输入，X3＝tanh(X2_out×W2+bias2)，为双向循环神经网络每一个时刻256个输出的组合，256个输出中每一个序列的维度是n_cls，X3＝{X31,X32,X33,X34,X35}；
S7.3:构建softmax层，将X3作为softmax分类器的输入特征，进行分类，由S1可知雷达HRRP样本来自5个不同的目标，训练得到softmax回归模型的参数集θ之后，雷达HRRP样本的特征X3对应于第i个目标的概率值p(i|X3)定义为：
Figure BDA0001782160770000071
其中，i表示目标的标号，i＝1,2,…,5，5为总目标数，θi表示对应于目标i的参数，
Figure BDA0001782160770000072
表示归一化向量，使雷达HRRP样本对应于所有目标模型的概率和为1。
本发明的有益效果是：
1、本发明加入了基于LSTM的双向循环神经网络，该网络不仅可以较好的从两个方向上利用HRRP样本所蕴含的局部物理结构信息的依赖关系，并且基于LSTM结构的循环神经网络较原始神经网络可更好的存储和利用全局结构信息。
2、现有技术仅利用了HRRP中的一阶马尔可夫性，它只注重于当前距离单元和之前距离单元的局部结构关系而忽略了各距离单元之间整体所存在的物理关系，在本发明中，通过对HRRP的局部结构提取其嵌入式特征表示可以更好地反映该局部结构之间存在的关联。
3、本发明所采用的双向循环神经网络方法可以朝着两个方向建立局部机构的依赖关系，相比单向循环神经网络可以更好地利用HRRP整体的结构信息。
附图说明
图1：基于高分辨距离像结构嵌入特征与深度神经网络的雷达高分辨距离像识别方法的步骤流程图。
具体实施方式
参照图1，为本发明的一种基于高分辨距离像结构嵌入特征与深度神经网络的雷达高分辨距离像识别技术流程图，具体实施步骤如下：
S1:采集数据集，将Q个雷达采集到的HRRP数据集依据目标的种类进行合并，每种种类的数据分别在不同的数据段里选择训练样本和测试样本，在训练集和测试集样本的选取过程中，保证所选取的训练集的数据与雷达所成姿态涵盖测试数据集与雷达所成的姿态。各类目标训练集和测试集样本数的比例为8:2，将挑选的数据集记作
Figure BDA0001782160770000081
其中Xi表示第i个样本，yk表示第k种目标，一共采集了5种目标，i0表示样本总数。
S2:对经过S1所提取数据集中的每一个样本做对齐预处理，具体操作步骤如下：
S2.1:由S1所选的数据集T中的每个样本皆为256维度，原始数据样本表示为：Xi＝[x1,x2,...,x256]
S2.2:将三个原始样本按照横向粘贴在一起构成一个256×3新样本数据样本，表示为X′i＝abs([Xi,Xi,Xi])＝[x′1,x'2,…,x'768]，其中abs(·)为对其中每个元素的值取绝对值；
S2.3:构建一个256维的序列L＝[l1,l2,…,l256]，其中li＝exp(-k|i-128.5|)，k表示控制因子，k＞0，|·|表示求绝对值操作；
S2.4:求
Figure BDA0001782160770000082
得到d之后，从X'中取出第d+1到第d+256个连续的点作为对齐后的数据X”＝[x'd+1,x'd+2,…,x'd+256]；
S2.5:将T数据集中的每一个样本重复步骤S3.2、S3.3以及S3.4操作得到对齐后的样本数据，记为T1；
S3:对T1中的每一个样本进行能量归一化，得到归一化后的数据，归一化公式如下：
Figure BDA0001782160770000091
其中||·||2表示求2范数；
S4:对经过S3处理后的雷达训练目标数据，从Xi”'内提取每一帧的时域特征作为识别特征,记作
Figure BDA0001782160770000092
其中x(c)表示的是第c个目标的训练样本的时域特征集，其中提取时域特征之后的第i个样本记作
Figure BDA0001782160770000093
C为目标总数；
S5:循环切分经过处理后的数据集，具体操作步骤如下：
S5.1:
Figure BDA0001782160770000094
其中i指示第i个目标，k指示的是第k个HRRP；
S5.2:设定切分的向量的长度为p,使得为p为2的幂级数，如4；
S5.3:对数据集内的每一个样本，将其移位p-1次并进行连接，
Figure BDA0001782160770000095
若
Figure BDA0001782160770000096
则
Figure BDA0001782160770000097
举例说明切分：
Figure BDA0001782160770000098
切分：
Figure BDA0001782160770000099
切分后的数据集可以表示为
Figure BDA00017821607700000910
S6:提取样本的嵌入表示特征，具体操作步骤如下：
S6.1:对经过S5处理后的每一个
Figure BDA00017821607700000911
按10级进行非均匀量化，得到其量化编码特征
Figure BDA00017821607700000912
其中对于每个元素：
若y＜e-3.6，则量化为0
若e-3.6≤y＜e-3.2，则量化为1
若e-3.2≤y＜e-2.8，则量化为2
若e-2.8≤y＜e-2.4，则量化为3
若e-2.4≤y＜e-2.0，则量化为4
若e-2.0≤y＜e-1.6，则量化为5
若e-1.6≤y＜e-1.2，则量化为6
若e-1.2≤y＜e-0.8，则量化为7
若e-0.8≤y＜e-0.4，则量化为8
若e-0.4≤y＜e0，则量化为9
S6.2:设定窗口值K，计算向量与向量之间共同出现的矩阵；
S6.3:首先初始化一个空矩阵，记作V＝{Vij}i＝0,1…9999,j＝0,1…9999，即初始化的空矩阵有10000行和10000列，其中0000代表取出来的列向量[yi,yi+1,yi+2,yi+3]T的量化编码为[0,0,0,0]T，其余的也是一样，矩阵中的每一个元素Vij为量化编码i和量化编码j之间的共现次数；
S6.4:依次统计列向量之间在窗长度下共现的次数，如
Figure BDA0001782160770000101
若窗K＝3(为保证左右对称，K一般取奇数)，则当循环至向量
Figure BDA0001782160770000102
的时候，统计其左侧
Figure BDA0001782160770000103
和其右侧
Figure BDA0001782160770000104
不同向量共同出现的次数，我们遍历所有的量化编码特征，将不同量化编码特征共同出现的总次数对应地记录在空矩阵V中，其中Vij为量化编码特征
Figure BDA0001782160770000105
与量化编码特征
Figure BDA0001782160770000106
共同出现的次数，其中i1、i2、i3和i4分别代表i的千位、百位、十位以及个位数字，j1、j2、j3和j4分别代表j的千位、百位、十位以及个位数字；
S6.4:对矩阵V中的每个元素限制其动态范围，公式如下：V′ij＝10ln(1+Vij)，并将更新后的矩阵记作V'＝{V′ij}；
S6.5:根据循环切分的结果去V'中查找
Figure BDA0001782160770000111
所对应的嵌入表示，如
Figure BDA0001782160770000112
中的第m列为[0,1,3,8]T，即V'中行标号为第0138所对应的行向量即为该列的嵌入表示，该嵌入表示为10000×1维的，依次查找
Figure BDA0001782160770000113
中每一列所对应的嵌入表示并记录，遍历
Figure BDA0001782160770000114
中所有列之后即可得到
Figure BDA0001782160770000115
的嵌入表示矩阵A，其中A＝[a1(10000×1)，a2(10000×1)，…，a256(10000×1)]，其中ai(10000×1)代表
Figure BDA0001782160770000116
第i列所对应的嵌入表示，ai为10000×1维的向量；
S7:构建深度神经网络，该深度神经网络由三部分组成：首先是一个全连接网络层用于嵌入特征降维，然后降维后的数据输入至基于LSTM(的双向循环神经网络，最后采用softmax层对双向循环神经网络的输出进行识别，具体操作步骤如下：
S7.1:构建一层全连接神经网络，设置第一层(输入层)神经元数目为10000，第二层(输出层)神经元数目为300，将经过步骤S6.5处理后的A按列切分依次送入全连接神经网络，得到输出
Z＝[z1＝f(Wija1+b)，z2＝f(Wija2+b)，…，z256＝f(Wija256+b)]，其中zi代表第i个输出，为300×1维的，Wij代表输入层第i个神经元到输出层第j个神经元的权重，b代表输入层神经元到输出层神经元的偏置，所有的神经元共享偏置；
S7.2:构建基于LSTM(长短时记忆网络)的双向循环神经网络，具体操作步骤如下：
S7.2.1:双向循环神经网络为两个单向其方向相反的循环神经网络(前向网络和反向网络)的结合，在每一个时刻t，网络的输出由这两个方向相反的单向循环神经网络共同决定，设置n_in＝256，n_step＝300，n_nums＝128，n_cls＝5，其中n_in＝256基于LSTM的双向循环神经网络的输入层神经元数目，n_step＝300代表每一个神经元的输入向量，为一个300×1的向量，n_nums＝128代表LSTM隐藏层神经元的数目，n_cls＝5代表输出层神经元数目，即代表所挑选的目标类别数目为5类；
S7.2.2:对双向循环神经网络的权重和偏置进行初始化，输入层和隐层之间的权重记为W1，其维度为shape1＝(n_in,n_nums),隐层和输出层之间的权重记为W2，其维度为shape2＝(n_nums,n_cls)，权重初始化采用服从均值为0，方差为1的高斯分布，输入层和隐层之间的偏置记为bias1，其维度为shape3＝(n_nums)，隐层和输出层之间的偏置记为bias2，其维度为shape4＝(n_cls)，均初始化为常量1，双向的初始状态皆记为state＝0，前向循环神经网络的所有时刻共享前向的权重，偏置和初始，反向也是如此；
S7.2.3:由S7.1可知，全连接神经网络的输出Z作为循环神经网络的输入，则循环神经网络第一层的输出为X1_out＝tanh(Z×W1+bias1)，同时作为隐层的输入，其中tanh(·)是激活函数，根据tensorflow中的计算公式
lstm_fw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)
lstm_bw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)
output,_,_＝tf.contrib.rnn.static_bidirectional_rnn(lstm_fw_cell,lstm_bw_cell,X1_out)
其中，lstm_fw_cell代表前向的循环神经网络，lstm_bw_cell代表反向循环神经网络，forget_bias＝1.0代表LSTM不会忘记任何信息，output,_,_分别代表X1_out经过隐层后的输出，前向循环神经网络和反向循环神经网络被更新后的状态，取出序列中的最后一个X2_out＝output[-1]作为输出层的输入，X3＝tanh(X2_out×W2+bias2)，为双向循环神经网络每一个时刻256个输出的组合，256个输出中每一个序列的维度是n_cls，X3＝{X31,X32,X33,X34,X35}；
S7.3:构建softmax层，将X3作为softmax分类器的输入特征，进行分类，由S1可知雷达HRRP样本来自5个不同的目标，训练得到softmax回归模型的参数集θ之后，雷达HRRP样本的特征X3对应于第i个目标的概率值p(i|X3)定义为：
Figure BDA0001782160770000131
其中，i表示目标的标号，i＝1,2,…,5，5为总目标数，θi表示对应于目标i的参数，
Figure BDA0001782160770000132
表示归一化向量，使雷达HRRP样本对应于所有目标模型的概率和为1；
测试阶段：
S8:对由S1采集到的测试数据进行训练阶段的步骤S3、S4、S5、S6和S7操作；
S9:将经过S8所提取的嵌入表示特征送入S7构建的基于深度神经网络的模型进行测试求得最后的结果,并通过最大后验概率来对样本进行判决，将测试样本分类到最大目标条件似然值所对应的目标
Figure BDA0001782160770000133
其中e代表模型判断X3属于5个类别分别会有一个对应的概率，从判断的5个概率中挑选最大概率所属的类别赋给e，即为所判定的雷达目标。
经过上述9个步骤，即可得到本发明所提出的一种基于高分辨距离像结构嵌入与深度神经网络的雷达高分辨距离像识别模型。

Claims (7)

1.一种基于结构嵌入和深度神经网络的雷达HRRP识别方法，包括以下步骤；
S1:采集数据集，将Q个雷达采集到的HRRP数据集依据目标的种类进行合并，每种种类的数据分别在不同数据段里选择训练样本和测试样本，在训练集和测试集样本的选取过程中，保证所选取训练集的数据与雷达所成姿态涵盖测试数据集与雷达所成的姿态；
S2:对经过S1所挑选数据集中的每一个样本使用补偿对齐法进行对齐，得到对齐后的数据；
S3:对S2处理后的每一个样本进行能量归一化处理，得到归一化后的数据；
S4：对经过S3处理后的雷达训练目标数据，提取时域特征作为识别特征，并用于训练；
S5:循环切分经过处理后的数据集，设定切分的向量长度为p，使p为2的幂级数；之后对数据集内的每一个HRRP样本，将其移位p-1次并进行连接；
S6:提取样本的嵌入表示特征，将每一个经过S5处理后的HRRP样本按10级进行非均匀量化；设定窗口值，初始化一个空矩阵，并依次统计量化向量在窗长度下共现的次数，对矩阵中的对应元素赋值；最后对矩阵中的每个元素限制其动态范围，根据共现矩阵得到HRRP训练样本的结构嵌入特征；
S7:构建深度神经网络模型，并利用HRRP训练样本的结构嵌入特征训练该模型；该深度神经网络模型由三部分组成：首先是一个全连接网络层，用于嵌入特征降维，然后将降维后的数据输入至基于长短时记忆网络LSTM的双向循环神经网络，最后采用softmax层对双向循环神经网络的输出进行识别；
测试阶段：
S8:对由S1采集到的测试数据进行训练阶段的步骤S2、S3、S4、S5和S6操作；
S9:将经过S8所提取的嵌入表示特征送入S7构建的深度神经网络得到测试样本所对应的目标；
步骤S7具体操作步骤如下：
S7.1:构建一层全连接神经网络，设置第一层输入层神经元数目为10000，第二层输出层神经元数目为300，将经过步骤S6.5处理后的A按列切分依次送入全连接神经网络，得到输出
Z＝[z1＝f(Wija1+b),z2＝f(Wija2+b),…,z256＝f(Wija256+b)]，其中zi代表第i个输出，为300×1维的，Wij代表输入层第i个神经元到输出层第j个神经元的权重，b代表输入层神经元到输出层神经元的偏置，所有的神经元共享偏置；
S7.2:构建基于长短时记忆网络LSTM的双向循环神经网络，具体操作步骤如下：
S7.2.1:双向循环神经网络为两个单向其方向相反的循环神经网络，即前向网络和反向网络的结合，在每一个时刻t，网络的输出由这两个方向相反的单向循环神经网络共同决定，设置n_in＝256，n_step＝300，n_nums＝128，n_cls＝5，其中n_in＝256代表基于LSTM的双向循环神经网络的输入层神经元数目，n_step＝300代表每一个神经元的输入向量，为一个300×1的向量，n_nums＝128代表LSTM隐藏层神经元的数目，n_cls＝5代表输出层神经元数目，即代表所挑选的目标类别数目为5类；
S7.2.2:对双向循环神经网络的权重和偏置进行初始化，输入层和隐层之间的权重记为W1，其维度为shape1＝(n_in,n_nums),隐层和输出层之间的权重记为W2，其维度为shape2＝(n_nums,n_cls)，权重初始化采用服从均值为0，方差为1的高斯分布，输入层和隐层之间的偏置记为bias1，其维度为shape3＝(n_nums)，隐层和输出层之间的偏置记为bias2，其维度为shape4＝(n_cls)，均初始化为常量1，双向的初始状态皆记为state＝0，前向循环神经网络的所有时刻共享前向的权重，偏置和初始，反向也是如此；
S7.2.3:由S7.1可知，全连接神经网络的输出Z作为循环神经网络的输入，则循环神经网络第一层的输出为X1_out＝tanh(Z×W1+bias1)，同时作为隐层的输入，其中tanh(·)是激活函数，根据tensorflow中的计算公式lstm_fw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)lstm_bw_cell＝tf.contrib.rnn.BasicLSTMCell(n_nums,forget_bias＝1.0)output,_,_＝tf.contrib.rnn.static_bidirectional_rnn(lstm_fw_cell,lstm_bw_cell,X1_out)
其中，lstm_fw_cell代表前向的循环神经网络，lstm_bw_cell代表反向循环神经网络，forget_bias＝1.0代表LSTM不会忘记任何信息，output,_,_分别代表X1_out经过隐层后的输出，前向循环神经网络和反向循环神经网络被更新后的状态，取出序列中的最后一个X2_out＝output[-1]作为输出层的输入，X3＝tanh(X2_out×W2+bias2)，为双向循环神经网络每一个时刻256个输出的组合，256个输出中每一个序列的维度是n_cls，
X3＝{X31,X32,X33,X34,X35}；
S7.3:构建softmax层，将X3作为softmax分类器的输入特征，进行分类，由S1可知雷达HRRP样本来自5个不同的目标，训练得到softmax回归模型的参数集θ之后，雷达HRRP样本的特征X3对应于第i个目标的概率值p(i|X3)定义为：
Figure FDA0002493840190000031
其中，i表示目标的标号，i＝1,2,…,5，5为总目标数，θi表示对应于目标i的参数，
Figure FDA0002493840190000032
表示归一化向量，使雷达HRRP样本对应于所有目标模型的概率和为1。
2.如权利要求1所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S1中各类目标训练集和测试集样本数的比例为8:2，将挑选的数据集记作
Figure FDA0002493840190000033
其中Xi表示第i个样本，yk表示第k种目标，一共采集了5种目标，i0表示样本总数。
3.如权利要求2所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S2具体操作步骤如下：
S2.1:由S1所选的数据集T中的每个样本皆为256维度，原始数据样本表示为：Xi＝[x1,x2,...,x256]；
S2.2:将三个原始样本按照横向粘贴在一起构成一个256×3新样本数据样本，表示为X′i＝abs([Xi,Xi,Xi])＝[x′1,x′2,…,x′768]，其中abs(·)为对其中每个元素的值取绝对值；
S2.3:构建一个256维的序列L＝[l1,l2,…,l256]，其中li＝exp(-k|i-128.5|)，k表示控制因子，k＞0，|·|表示求绝对值操作；
S2.4:求
Figure FDA0002493840190000041
得到d之后，从X'中取出第d+1到第d+256个连续的点作为对齐后的数据X”＝[x′d+1,x′d+2,…,x′d+256]；
S2.5:将T数据集中的每一个样本重复步骤S2.2、S2.3以及S2.4操作得到对齐后的样本数据，记为T1。
4.如权利要求3所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S3中，归一化公式如下：
Figure FDA0002493840190000042
其中||·||2表示求2范数。
5.如权利要求4所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S4中，从X″i内提取每一帧的时域特征作为识别特征,记作
Figure FDA0002493840190000043
其中x(c)表示的是第c个目标的训练样本的时域特征集，其中提取时域特征之后的第i个样本记作
Figure FDA0002493840190000044
C为目标总数。
6.如权利要求5所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S5具体操作步骤如下：
S5.1:
Figure FDA0002493840190000045
其中i指示第i个目标，k指示的是第k个HRRP；
S5.2:设定切分的向量的长度为p,使得为p为2的幂级数；
S5.3:对数据集内的每一个样本，将其移位p-1次并进行连接，
Figure FDA0002493840190000046
若
Figure FDA0002493840190000051
则
Figure FDA0002493840190000052
7.如权利要求5所述的基于结构嵌入和深度神经网络的雷达HRRP识别方法，其特征在于：步骤S6具体操作步骤如下：
S6.1:对经过S5处理后的每一个
Figure FDA0002493840190000053
按10级进行非均匀量化，得到其量化编码特征
Figure FDA0002493840190000054
其中对于每个元素：
若y＜e-3.6，则量化为0
若e-3.6≤y＜e-3.2，则量化为1
若e-3.2≤y＜e-2.8，则量化为2
若e-2.8≤y＜e-2.4，则量化为3
若e-2.4≤y＜e-2.0，则量化为4
若e-2.0≤y＜e-1.6，则量化为5
若e-1.6≤y＜e-1.2，则量化为6
若e-1.2≤y＜e-0.8，则量化为7
若e-0.8≤y＜e-0.4，则量化为8
若e-0.4≤y＜e0，则量化为9
S6.2:设定窗口值K，计算向量与向量之间共同出现的矩阵；
S6.3:首先初始化一个空矩阵，记作V＝{Vij}i＝0,1…9999,j＝0,1…9999，即初始化的空矩阵有10000行和10000列，其中0000代表取出来的列向量[yi,yi+1,yi+2,yi+3]T的量化编码为[0,0,0,0]T，其余的也是一样，矩阵中的每一个元素Vij为量化编码i和量化编码j之间的共现次数；
S6.4:依次统计列向量之间在窗长度下共现的次数，如
Figure FDA0002493840190000055
为保证左右对称，K取奇数，若窗K＝3，则当循环至向量
Figure FDA0002493840190000056
的时候，统计其左侧
Figure FDA0002493840190000057
和其右侧
Figure FDA0002493840190000058
不同向量共同出现的次数，我们遍历所有的量化编码特征，将不同量化编码特征共同出现的总次数对应地记录在空矩阵V中，其中Vij为量化编码特征
Figure FDA0002493840190000061
与量化编码特征
Figure FDA0002493840190000062
共同出现的次数，其中i1、i2、i3和i4分别代表i的千位、百位、十位以及个位数字，j1、j2、j3和j4分别代表j的千位、百位、十位以及个位数字；
S6.4:对矩阵V中的每个元素限制其动态范围，公式如下：V′ij＝10ln(1+Vij)，并将更新后的矩阵记作V'＝{V′ij}；
S6.5:根据循环切分的结果去V'中查找
Figure FDA0002493840190000063
所对应的嵌入表示，如
Figure FDA0002493840190000064
中的第m列为[0,1,3,8]T，即V'中行标号为第0138所对应的行向量即为该列的嵌入表示，该嵌入表示为10000×1维的，依次查找
Figure FDA0002493840190000065
中每一列所对应的嵌入表示并记录，遍历
Figure FDA0002493840190000066
中所有列之后即可得到
Figure FDA0002493840190000067
的嵌入表示矩阵A，其中A＝[a1(10000×1),a2(10000×1),…,a256(10000×1)]，其中ai(10000×1)代表
Figure FDA0002493840190000068
第i列所对应的嵌入表示，ai为10000×1维的向量。
CN201810998871.5A 2018-08-29 2018-08-29 基于结构嵌入和深度神经网络的雷达hrrp识别方法 Active CN109239670B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810998871.5A CN109239670B (zh) 2018-08-29 2018-08-29 基于结构嵌入和深度神经网络的雷达hrrp识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810998871.5A CN109239670B (zh) 2018-08-29 2018-08-29 基于结构嵌入和深度神经网络的雷达hrrp识别方法

Publications (2)

Publication Number Publication Date
CN109239670A CN109239670A (zh) 2019-01-18
CN109239670B true CN109239670B (zh) 2020-08-04

Family

ID=65068801

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810998871.5A Active CN109239670B (zh) 2018-08-29 2018-08-29 基于结构嵌入和深度神经网络的雷达hrrp识别方法

Country Status (1)

Country Link
CN (1) CN109239670B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110045362B (zh) * 2019-03-01 2022-02-22 合肥工业大学 基于rnn的雷达目标航迹识别方法和装置
CN110109109B (zh) * 2019-04-26 2023-03-21 西安电子科技大学 基于多分辨率注意力卷积网络的hrrp目标识别方法
CN111444790B (zh) * 2020-03-13 2022-07-01 北京理工大学 一种对多功能雷达工作模式序列的脉冲级智能识别方法
CN112732871B (zh) * 2021-01-12 2023-04-28 上海畅圣计算机科技有限公司 一种机器人催收获取客户意向标签的多标签分类方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8872693B1 (en) * 2009-04-29 2014-10-28 The United States of America as respresented by the Secretary of the Air Force Radar signature database validation for automatic target recognition
CN107121679A (zh) * 2017-06-08 2017-09-01 湖南师范大学 用于雷达回波外推的循环神经网络预测法及记忆单元结构
CN107463966B (zh) * 2017-08-17 2019-06-18 电子科技大学 基于双深度神经网络的雷达一维距离像目标识别方法
CN108133232B (zh) * 2017-12-15 2021-09-17 南京航空航天大学 一种基于统计字典学习的雷达高分辨距离像目标识别方法
CN108009154B (zh) * 2017-12-20 2021-01-05 哈尔滨理工大学 一种基于深度学习模型的图像中文描述方法
CN108254741B (zh) * 2018-01-16 2021-02-09 中国人民解放军海军航空大学 基于循环神经网络的目标航迹预测方法

Also Published As

Publication number Publication date
CN109239670A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN109214452B (zh) 基于注意深度双向循环神经网络的hrrp目标识别方法
CN110443143B (zh) 多分支卷积神经网络融合的遥感图像场景分类方法
CN109239670B (zh) 基于结构嵌入和深度神经网络的雷达hrrp识别方法
CN109086700B (zh) 基于深度卷积神经网络的雷达一维距离像目标识别方法
CN111583263B (zh) 一种基于联合动态图卷积的点云分割方法
CN108090472B (zh) 基于多通道一致性特征的行人重识别方法及其系统
CN112116563A (zh) 一种基于谱维与空间协作邻域注意力的高光谱图像目标检测方法与系统
CN111736125B (zh) 基于注意力机制和双向堆叠循环神经网络的雷达目标识别方法
CN109743642B (zh) 基于分层循环神经网络的视频摘要生成方法
CN107169117B (zh) 一种基于自动编码器和dtw的手绘图人体运动检索方法
CN113297936B (zh) 一种基于局部图卷积网络的排球群体行为识别方法
CN111814685A (zh) 基于双支路卷积自编码器的高光谱图像分类方法
CN109446897B (zh) 基于图像上下文信息的场景识别方法及装置
CN114926693A (zh) 基于加权距离的sar图像小样本识别方法及装置
CN114937173A (zh) 一种基于动态图卷积网络的高光谱图像快速分类方法
CN115048870A (zh) 基于残差网络和注意力机制的目标轨迹识别方法
CN111596292A (zh) 基于重要性网络和双向堆叠循环神经网络的雷达目标识别方法
CN112132257A (zh) 基于金字塔池化及长期记忆结构的神经网络模型训练方法
CN116206158A (zh) 基于双超图神经网络的场景图像分类方法及系统
Lv et al. Region-based adaptive association learning for robust image scene recognition
CN109871907B (zh) 基于sae-hmm模型的雷达目标高分辨率距离像识别方法
CN113344069B (zh) 一种基于多维关系对齐的无监督视觉表征学习的图像分类方法
CN114187506A (zh) 视点意识的动态路由胶囊网络的遥感图像场景分类方法
Kumar et al. Region driven remote sensing image captioning
CN112836736B (zh) 一种基于深度自编码器构图的高光谱图像半监督分类方法

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190118

Assignee: Hangzhou kunjun Ocean Engineering Technology Co., Ltd

Assignor: HANGZHOU DIANZI University

Contract record no.: X2020330000110

Denomination of invention: Radar HRRP recognition method based on structure embedding and depth neural network

Granted publication date: 20200804

License type: Common License

Record date: 20201203