基于一维堆叠池化融合卷积自编码器的HRRP目标识别方法

第43卷 第12期系统工程与电子技术Vol.43 No.122021年12月SystemsEngineering                                                  andElectronicsDecember2021
文章编号:1001 506X(2021)12 3533 09 网址:www.sys ele.com收稿日期:20200628;修回日期:20200724;网络优先出版日期:20210315。网络优先出版地址:https:∥kns.cnki.net/kcms/detail/11.2422.TN.20210315.1551.008.html基金项目:国家自然科学基金(61876189,61273275,61503407,61806219,61703426)资助课题 通讯作者.引用格式:张国令,吴崇明,李睿,等.基于一维堆叠池化融合卷积自编码器的HRRP目标识别方法[J].系统工程与电子技术,2021,43(12):3533 3541.犚犲犳犲狉犲狀犮犲犳狅狉犿犪狋:ZHANGGL,WUCM,LIR,etal.HRRPtargetrecognitionmethodbasedonone dimensionalstackedpoolingfusionconvolutionalautoencoder[J].SystemsEngineeringa
卷积编码
ndElectronics,2021,43(12):3533 3541.基于一维堆叠池化融合卷积自编码器的犎犚犚犘目标识别方法
张国令1,吴崇明2, ,李 睿1,来 杰1,向 前1
(1.空军工程大学防空反导学院,陕西西安710051;2.西京学院,陕西西安710123)
  摘 要:针对雷达高分辨距离像(highresolutionrangeprofile,HRRP)目标识别中特征提取及识别问题,提出了一种基于一维堆叠池化融合卷积自编码器(one dimensionalstackedpoolingfusionconvolutionalautoencoder,1DSPF CAE)的识别方法。首先构造一维池化融合卷积自编码器(one dimensionalpoolingfusionconvolutionalauto encoder,1DPF CAE),在编码阶段,采用最大池化和平均池化同时提取不同的编码特征并进行融合来提取HRRP的结构特征;然后堆叠多个1DPF CAE形成1DSPF CAE;最后使用标签数据对网络进行微调,实现HR RP目标识别。并使用AdaBound算法优化网络训练来提高识别性能。基于弹道中段目标仿真数据的实验结果表明,该方法具有较强的特征提取能力,对于HRRP目标识别准确率高、鲁棒性强。关键词:雷达自动目标识别;高分辨距离像;卷积自编码器;特征提取;池化融合中图分类号:TP183    文献标志码:A    犇犗犐:10.12305/j.issn.1001 506X.2021.12.15犎犚犚犘狋犪狉犵犲狋狉犲犮狅犵狀犻狋犻狅狀犿犲狋犺狅犱犫犪狊犲犱狅狀狅狀犲 犱犻犿犲狀狊犻狅狀犪犾狊狋犪犮犽犲犱狆狅狅犾犻狀犵
狌狊犻狅狀犮狅狀狏狅犾狌狋犻狅狀犪犾犪狌狋狅犲狀犮狅犱犲狉ZHANGGuoling1,WUChongming2, ,L
IRui1,LAIJie1,XIANGQian1(1.犃犻狉犪狀犱犕犻狊狊犻犾犲犇犲犳犲狀狊犲犆狅犾犾犲犵犲,犃犻狉犉狅狉犮犲犈狀犵犻狀犲犲狉犻狀犵犝狀犻狏犲狉狊犻狋狔,犡犻’犪狀710051,犆犺犻狀犪;2.犡犻犼犻狀犵犆狅犾犾犲犵犲,犡犻’犪狀710123,犆犺犻狀犪)  犃犫狊狋狉犪犮狋:Aimingattheproblemoffeatureextractionandrecognitioninhighresolutionrangeprofile(HRRP)targetrecognition,arecognitionmethodbasedonone dimensionalstackedpoolingfusionconvolutionalautoencoder(1DSPF CAE)isproposedinthispaper.Firstly,aone dimensionalpoolingfusionconvolutionalautoencoder(1DPF CAE)isconstructed.Intheencodingstage,themaximumpoolingandaveragepoolingareusedtoextractdifferentencodingfeaturesandfusethemtoextractthestructural
featuresofHRRP.Then,multiple1DPF CAEsarestackedtoform1DSPF CAE.Finally,thenetworkisfine tunedusingl
abeldatatorealizeHRRPtargetrecognition.AndtheAdaBoundalgorithmisusedtooptimizenetworktrainingforimprovingtherecognitionperformance.Theexperimentalresultsbasedonthesimulateddataofthetargetinthemiddlepartofthetrajectoryshowthatthemethodhasstrongfeatureextractioncapability,andhashighaccuracyandrobustnessforHRRPtargetrecognition.犓犲狔狑狅狉犱狊:radarautomatictargetrecognition;highresolutionrangeprofile;convolutionautoencoder;featureextraction;poolingf
usion0 引 言雷达高分辨率距离像(highresolutionrang
eprofile,HRRP)不仅能够反映目标散射点沿雷达视线方向的分布情况,而且包含目标的结构特征,为目标识别提供了重要的目标结构特征信息。且相较于二维距离像,一维HRRP更
 ·3534 ·系统工程与电子技术第43卷                                                   易获取、处理以及存储,因此HRRP目标识别成为雷达自动目标识别研究领域的热点[13]。如何提取有效特征是HR RP目标识别的关键,现已有不少相关研究[46]。如基于核主分量相关判别分析的特征提取法[4]、基于字典学习的鲁棒性特征提取法[5]、基于特权信息的特征提取法[6]等。由于这些方法基于浅层线性结构,所以提取的特征表达能力受限,而且特征需要人为设计,限制了识别性能的提升。因此,如何自动提取目标的深层特征成为雷达目标识别的研究热点。深度学习本质上是利用深度神经网络自动提取目标的深层特征,目前已经在HRRP识别领域取得一些研究成果[710]。Pan等[7]提出了一种基于t分布随机邻域嵌入(t distributedstochasticneighborembedding,t SNE)和判别深度置信网络(discriminantdeepbeliefnetwork,DDBN)的HRRP识别模型,该模型利用了HRRP的深层特征进行识别。杨予昊等[8]提出一种基于卷积神经网络(convolutionalneuralnetwork,CNN)的HRRP识别方法,充分发掘HR RP所包含的深层属性特征,并实现了自动特征提取。Feng等[9]提出一种用于HRRP识别的堆叠校正自编码器(stackedcorrectiveautoencoder,SCAE),实验表明该方法具
有较强的特征提取能力。Zhao等[10]利用堆叠自编码器(stackedautoencoder,SAE)提取HRRP高层次特征,利用极限学习机(extremelearningmachine,ELM)代替反向传播(backpropagation,BP)网络,实现了对小样本HRRP数据集的快速准确识别。作为一种深度神经网络,卷积自编码器(convolutionalautoencoder,CAE)因其提取数据深层特征的能力受到了广泛关注[1114]。Chen等[11]利用CAE通过无监督学习提取CT图像特征,成功地应用于医疗领域。Huang等[12]采用卷积核为自编码器的基本单元,提出深度CAE(deepCAE,DCAE),利用CNN的层次特征提取的优势来提取脑功能核磁共振数据中的特征,取得了较好的效果。Luo等[13]提出一种卷积稀疏自编码器(convolutionalsparseautoencoder,CSAE),利用CSAE提取到的图像局部特征初始化CNN进行分类。Du等[14]提出一种堆叠卷积去噪自编码器(stackedconvolutionaldenoisingautoencoder,SCDAE),该网络以卷积方式叠加多层去噪自编码器来实现分层训练优化,每层利用去噪自编码器学习的核函数对下层的特征进行卷积,生成高维特征映射。CAE采用了自编码器的思想进行无监督特征学习,通过卷积和池化操作来提取数据的深层结构特征。从HRRP目标识别的研究现状来看,利用CAE进行HRRP识别的研究很少,本研究旨在设计出性能良好的CAE网络,提取HRRP有效深层特征,实现目标识别。为了充分利用HRRP数据的内部结构信息,进一步提高雷达HRRP目标的识别准确率,本文提出了一种一维池化融合C
AE(one dimensionalpoolingfusionCAE,1DPF CAE)。在编码阶段,利用最大池化和平均池化提取特征并融合,两者提取到的特征互相补充丰富,然后将1DPF CAE堆叠形成一维堆叠池化融合CAE(one dimensionstackedpoolingfusionCAE,1DSPF CAE),从而增强其深层特征提取能力,进一步提出了基于1DSPF CAE的HRRP目标识别方法。同时,考虑到深层网络固有的训练困难问题,使用Ada Bound[15]优化算法优化网络训练,以提高网络的整体性能。实验结果表明,本文提出的基于1DSPF CAE的HRRP目标识别方法取得了良好的效果。1 犆犃犈原理1.1 犆犃犈自动编码器由两部分组成:编码器和解码器,如图1所示
。1 Fig.1 Architectureofautoencoder编码器使用确定性的映射函数将输入狓映射到隐层表示狔,通常映射函数犳是非线性的,狔可表示为狔=犳(犠狓+犫)(1)式中:犠表示输入和隐层表示之间的权重;犫表示偏置。解码器通过隐层表示狔重构输出狕,公式为狕=犳′(犠′狔+犫′)(2)式中:犠′是隐层表示与输出之间的权重;犫′是偏置。自编码器通过最小化重构误差进行训练,假设训练集犇={狓(犻)}犖犻=1,则自编码器的损失函数为犔(
狓,狕)=1犖∑狓∈犇‖狓-狕‖2(3)  CAE将卷积连接和自编码器结合在一起,由卷积编码器和卷积解码器组成[16]。卷积编码器实现输入到特征的卷积转换,卷积解码器实现特征到输出的卷积转换。在CAE中,通过CNN计算提取的特征和重构输出。因此,式(1)和式(2)可改写为狔=ReLU(狑狓+犫)(4)狕=ReLU(狑′狔+犫′)(5)式中:狑为狓和狔之间的卷积核;狑′为狔和狕之间的卷积核。一维堆叠卷积自编码器(one dimensionalstackedCAE,1D SCAE)是一种特殊的卷积自编码器,其原理如下节所述。1.2 1犇 犛犆犃犈1D SCAE由多个一维CAE(one dimensionCAE,1D CAE)分层堆叠而成,前一个1D CAE编码器的输出为后一个1D CAE的输入[17]。1D SCAE能够学习到鲁棒性强、抽象度高的特征,主要有两大原因:一是与单层自动编码器相比,1D SCAE堆叠起来的深层网络结构具有更强的特征表征能力[17];二是与全连接网络相比,1D SCAE网络采用卷积和池化操作,能够在提取局部特征的同时,从全局上把握整体的变化特征,使得提取到的特征更具鲁棒性[18]。图2是由两个1D CAE堆叠而成的1D SCAE。
 第1
2期张国令等:基于一维堆叠池化融合卷积自编码器的HRRP目标识别方法
·3535 ·                                                   图2 1D SCAEFig.2 1D SCAE每层
1D CAE在编码阶段通过一维卷积和池化提取数据的特征,在解码阶段通过上采样和一维反卷积运算重建输入。(1)编码过程由一维卷积和池化组成。其公式为犺犻,1=ReLU(ω犻 狓犻+犫犻)(6)犺犻,2=pool(犺犻,1)(7)式中:狓犻=犺犻-1(犻=1,2,…,狀),犻=1时,狓犻=狓;犺犻,1为第犻个编码器的一维卷积输出;犺犻,2为第犻个编码器输出的特征;ω犻和犫犻分别为第犻个编码器的卷积核和偏置; 表示卷积操作;pool表示池化操作。(2)解码过程由上采样和一维反卷积组成。其公式为犺-犻,1=upsample(犺犻,2)(8)犺-犻,2=ReLU(ω-犻,1 犺-犻,1+犫-犻,1)(9)狓-犻=ReLU(ω-犻,2 犺-犻,2+犫-犻,2)(10)式中:狓-犻是输入数据狓犻的重构;犺-犻,1为第犻个编码器的上采样输出;ω-犻,1和犫-犻,1分别是第犻个1D CAE第一个一维反卷积的卷积核和偏置;ω-犻,2和犫-犻,2分别是第犻个1D CAE第二个一维反卷积的卷积核和偏置,且ω-犻,2卷积核数为1,从而使狓-犻和狓犻保持一致;upsample表示上采样操作。各层1D CAE通过最小化重构误差进行训练,第犻个1D CAE的损失函数为犔(狓犻,狓-犻)=12犖∑犖狀=1‖狓狀犻-狓-狀犻‖2(11)  为进一步提高1D CAE的特征提取能力,在上述经典1D CAE的基础上,本文在1D CAE编码阶段采用最大池化和平均池化同时提取特征,并将提取到的特征进行融合,设计了一种1DPF CAE,进一步将堆叠的1DPF CAE用于HRRP特征的自动提取,从而提出一种基于1DSPF CAE的HRRP识别方法。2 基于1犇犛犘犉 犆犃犈的犎犚犚犘识别方法2.1 1犇犛犘犉 犆犃犈传统的1D CAE因其网络模型采用
了CNN,具有局部连接、参数共享、池化、多卷积核等优势;而且经过池化操作得到的特征对于小范围内的旋转和平移具有较好的鲁棒性,因此具有较强的特征表达能力[19]。为进一步提高网络的表达能力,本文设计了一种1DPF CAE,在编码阶段,同时采用最大池化和平均池化提取不同的编码特征,并将两者提取到的特征进行串联连接为最终的池化输出。在此基础上,将1DPF CAE形成1DSPF CAE,其网络结构如图3所示
。3 1DSPF CAEFig.
3 Networkstructureof1DSPF CAE
 ·3536 ·系统工程与电子技术第43卷                                                     其中,堆叠的三层1DPF CAE主要用于特征提取,后面的网络主要用于分类。1DSPF CAE的每层1DPF CAE都通过无监督学习来提取特征,其无监督学习分为卷积编码和卷积解码两个阶段。编码阶段,首先进行两次一维卷积操作,目的在于充分提取数据的局部结构特征;然后同时进行最大池化和平均池化操作,并将两部分不同特征进行融合作为池化输出。编码的具体过程如下所示:犺犻,1=ReLU(ω犻,1 狓犻+犫犻,1)(12)犺犻,2=ReLU(ω犻,2 犺犻,1+犫犻,2)(13)犺犻,3_max=maxpool(犺犻,2)(14)犺犻,3_avg=avgpool(犺犻,2)(15)犺犻,3=犺犻,3_max 犺犻,3_avg(16)式中:狓犻=犺犻-1(犻=1,2,3),犻=1时,狓犻=狓;犺犻,1为第犻个编码器的第1个一维卷积输出;犺犻,2为第犻个编码器的第2个一维卷积输出;犺犻,3_max为第犻个编码器的最大池化输出;犺犻,3_avg为第犻个编码器的平均池化输出;犺犻,3为第犻个编码器输出的特征;maxpool表示最大池化操作,avgpool表示平均池化操作; 代表串联连接。解码的过程与第1.2节中的解码过程等同,即式(8)~式(10)。网络含有大量可训练的参数θ,容易导致过拟合,因此在全连接层加入Dropout技术[20],以概率随机舍弃部分节点,减少了冗余特征,缓
解了过拟合现象。2.2 训练优化1DPF CAE网络的训练采用小批量训练模式。随机梯度下降(stochasticgradientdescent,SGD)法分批次随机抽取犿个训练样本作为输入,第狋次迭代时,参数更新为θ狋+1=θ狋-α犵狋(17)式中:犵狋为损失函数的梯度;α为学习率。SGD使用固定的学习率进行参数更新,存在学习率选择困难的问题,因此1DPF CAE采用AdaBound算法对参数更新过程进行优化。每次参数更新时,AdaBound利用梯度一阶矩和二阶矩估计分别更新犵狋和α,同时施加动态边界限制了α的变化范围,使得参数更新过程更平稳。基于AdaBound的1DPF CAE训练优化算法见算法1。算法1 基于AdaBound的1DPF CAE训练优化算法输入:训练集、测试集输出:最优模型1参数配置:设置epochs为100,batchsize为200,学习率α0=0.001,α犳=0.1,β1=0.9,β2=0.999,以及小常数γ=10-3,ε=10-7。2模型初始化:设置模型参数θ为近似0的随机值。3参数更新。4初始化梯度一阶矩估计:犿0←05初始化梯度二阶矩估计:狏0←06for狋=1to犜do7根据式(11)计算前向传播误差:犔(θ狋-1)8计算损失函数的梯度:犵狋←Δ
θ狋-1犔(θ狋-1)9更新α:α狋←α01-β狋槡21-β狋110计算α的下界:α犾←α犳(1-1γ狋+1)11计算α的上界:αμ←α犳(1+1γ狋)12更新梯度一阶矩估计:犿狋←β狋犿狋-1+(1-β1)犵狋13更新梯度二阶矩估计:狏狋←β2狏狋-1+(1-β2)犵2狋14更
新θ:θ狋←θ狋-1-犿狋·min{max{α狋狏槡狋+ε,α犾},αμ}15endfor训练好每个1DPF CAE之后将其堆叠构造1DSPF CAE,然后再使用AdaBound算法对整个网络进行微调,最终采用Softmax分类器实现分类。2.3 基于1犇犛犘犉 犆犃犈的犎犚犚犘目标识别方法本节利用提出的1DSPF CAE对雷达HRRP目标进行识别,凭借1DPF CAE提取特征能力强的优势,提取HRRP的深层结构特征;然后堆叠多个1DPF CAE形成1DSPF CAE;最后使用标签数据对网络进行微调,实现HRRP目标识别。基于1DSPF CAE的HRRP目标识别方法主要包括3个阶段:预处理阶段、训练阶段和测试阶段。(1)预处理阶段:对原始的HRRP数据进行去噪处理和归一化处理。(2)训练阶段:训练阶段分为预训练和微调。首先对3个1DPF CAE单独训练,然后将训练好的3个1DPF CAE按照图3进行堆叠,再对Softmax预训练。在微调阶段,对整个网络进行微调,优化1DSPF CAE中的每个1DPF CAE和Softmax分类器的参数,得到最终的训练模型。(3)测试阶段:利用训练好的网络对测试集HRRP数据进行分类得到目标类型。3 实验与结果分析3.1 实验设置为了测试本文提出方法的性能,在此设计了以下实验。实验1 模型的训练实验。实验2 不同池化方式对方法性能的影响分析。实验3 特征可视化分析实验。实验4 参数影响分析实验,分析了不同堆叠层数、不同优化算法以及不同Dropout取值对方法性能的影响。实验5 与其他算法在加噪和未加噪数据集上的性能对比分析。
 第12期张国令等:基于一维堆叠池化融合卷积自编码器的HRRP目标识别方法·3537 ·                                                   3.1.1 实验环境及实验数据本文试验环境为Tensorflow,计算机配置为Intel(R)Core(TM)i7 4790CPU3.60GHz,16GBRAM。本文使用FEKO和Matlab联合仿真得到5类弹道中段目标的HRRP仿真数据。其中,雷达参数设置如下:中心频率10GHz、带宽1GHz、方位角范围0°~180°、方位角间隔0.05°。5类仿真目标及其物理参数如图4所示。图4 5类弹道中段目标Fig.4 Fivetypesofthetargetinthemiddlepartofthetrajectory仿真得到18005个256维的一维HRRP样本数据,每类目标HRRP数量为3601。按照66%、33%的比例随机划分训练集和测试集,构造数据集A;按照50%、50%的比例随机划分训练集和测试集,构造数据集B;按照33%、66%的比例随机划分训练集和测试集,得到数据集C。在不加说明的情况下,以下均将数据集A作为实验数据集。3.1.2 模型参数设置1DSPF CAE训练参数设置如下:学习率为0.001,批量大小为200,Dropout率为0.2。网络结构设置如表1所示。表1 网络参数设置犜犪犫犾犲1 犖犲狋狑狅狉犽狆犪狉犪犿犲狋犲狉狊狊犲狋狋犻狀犵网络层输入核数核/窗口步长输出C1(256×1)×1323×11(256×1)×32C1(256×1)×1323×11(256×1)×32P1_max(256×1)×32—2×12(128×1)×32P1_avg(256×1)×32—2×12(128×1)×32P1————(128×1)×64C2(128×1)×64483×11(128×1)×4
8C2(128×1)×48483×11(128×1)×48P2_max(128×1)×48—2×12(64×1)×48P2_avg(128×1)×48—2×12(64×1)×48P2————(64×1)×96C3(64×1)×96643×11(64×1)×64C3(64×1)×64643×11(64×1)×64P3_max(64×1)×64—2×12(32×1)×64P3_avg(64×1)×64—2×12(32×1)×64P3————(32×1)×128Flatten(32×1)×128———4096×1全连接4096×1———400×1Softmax400×1———5×1
3.2 模型训练在预训练阶段,1DSPF CAE中的每一层1DPF CAE进行无监督式学习,通过对输入数据重构提取数据的特征。在微调阶段,将各层1DPF CAE和Softmax分类器看作一个整体去调整自由参数,参数调优的过程直观表现为损失不断降低,当损失收敛时得到的参数是最优的。实验进行100次迭代,
训练集和测试集的损失变化如图5所示。图5 损失变化Fig.5 Lossvariation由图5看出,在训练初期,损失随着迭代次数的增加而快速下降,训练后期趋于稳定且训练误差较低。说明1DSPF CAE网络具有良好的训练效果,收敛速度快,并且收敛后训练误差小。

本文发布于:2024-09-21 20:30:20,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/1/378044.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:特征   训练   卷积   提取   目标
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议