CN104007431B - 基于dpLVSVM模型的雷达HRRP的目标识别方法 - Google Patents

基于dpLVSVM模型的雷达HRRP的目标识别方法 Download PDF

Info

Publication number
CN104007431B
CN104007431B CN201410234677.1A CN201410234677A CN104007431B CN 104007431 B CN104007431 B CN 104007431B CN 201410234677 A CN201410234677 A CN 201410234677A CN 104007431 B CN104007431 B CN 104007431B
Authority
CN
China
Prior art keywords
power spectrum
resolution range
range profile
radar
cluster
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
CN201410234677.1A
Other languages
English (en)
Other versions
CN104007431A (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 CN201410234677.1A priority Critical patent/CN104007431B/zh
Publication of CN104007431A publication Critical patent/CN104007431A/zh
Application granted granted Critical
Publication of CN104007431B publication Critical patent/CN104007431B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于dpLVSVM模型的雷达HRRP的目标识别方法，其步骤为：步骤1，对雷达HRRP数据进行特征提取得到功率谱特征集X；步骤2，构建dpLVSVM模型并得出功率谱特征的概率密度函数及各个参数的联合条件后验分布；步骤3，推导各个参数的条件后验分布；步骤4，对各个参数进行循环采样I次；步骤5，保存T0次测试阶段所需参数的采样结果；步骤6，判断测试雷达HRRP是否为库外样本，若为库外样本则拒判；否则步骤7；步骤7，采样得到测试雷达HRRP的功率谱特征的聚类标号；步骤8，输出测试雷达HRRP的目标类别标号本发明具有分类器设计复杂度小，识别性能高以及拒判性能好的优点，可用于对雷达目标的识别。

Description

基于dpLVSVM模型的雷达HRRP的目标识别方法
技术领域
本发明属于雷达技术领域，涉及雷达目标识别方法，尤其涉及一种基于dpLVSVM(Dirichlet process latent variable support vector machine，Dirichlet过程隐变量支撑向量机)模型的雷达高分辨距离像HRRP的目标识别方法，用于对飞机，车辆等目标进行识别。
背景技术
雷达目标识别就是利用目标的雷达回波信号，实现对目标类型的判定。宽带雷达通常工作在光学区，此时目标可以看作是由大量强度不同的散射点构成。高分辨距离像(High-resolution range profile,HRRP)是用宽带雷达信号获取的目标体上各散射点回波的矢量和。它反映了目标体上散射点沿雷达视线的分布情况，包含了目标重要的结构特征，被广泛应用于雷达目标识别领域。由于目标具有姿态敏感性，同一目标的HRRP具有多模分布特性，尤其是随着目标库的增大，训练样本个数也会随之增加，数据分布也变得更加复杂。多模分布数据的分类界面往往是高度非线性的，需要采用非线性分类器对其分类。
作为一种常用的非线性分类器，核方法分类器是将原始空间线性不可分的数据映射成为高维空间中线性可分的数据，然后进行线性分类。然而核方法分类器面临核函数选择以及核参数选择的问题，且当训练样本数过大时，核方法分类器计算困难。另外，若使用所有雷达高分辨距离像数据来训练一个分类器会增加分类器的训练复杂度，而且容易忽视样本的内在结构，不利于分类。混合专家模型提出后，避免了复杂分类器设计，从而大大简化分类器设计的复杂度。
混合专家模型将数据集划分成若干子集，然后在各个子集上分别训练简单的分类器，最终构造全局非线性的复杂分类器，称为有限混合专家模型模型。这类模型存在两个缺点：一是模型选择问题，即如何选择样本子集(聚类)个数；二是样本集的聚类过程是无监督的，独立于后端的分类器任务，因此较难保证每个聚类中数据的可分性，从而影响全局的分类性能。
发明内容
为了克服以上困难，本发明提出一种基于dpLVSVM模型的雷达HRRP的目标识别方法，用于提高分类性能，降低模型求解复杂度。
为达到上述目的，本发明采用以下技术方案予以实现：
一种基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，包括以下步骤：
步骤1，雷达接收M个类别的目标的高分辨距离像HRRP；再对每个高分辨距离像进行特征提取，得到雷达高分辨距离像的功率谱特征xn，将M个类别的目标的高分辨距离像的功率谱特征组成功率谱特征集X；不属于该M个类别的目标的雷达高分辨距离像为库外样本；属于该M个类别的目标的雷达高分辨距离像为库内样本；
步骤2，利用功率谱特征集X，将LVSVM分类器和TSB-DPM模型结合构建dpLVSVM模型；根据dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数，以及dpLVSVM模型各个参数的联合条件后验分布；
dpLVSVM模型各个参数的联合条件后验分布即：TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布；
步骤3，通过步骤2的dpLVSVM模型各个参数的联合条件后验分布，推导各个参数的条件后验分布，即聚类的高斯分布参数的条件后验分布、雷达高分辨距离像的功率谱特征的聚类标号Z的条件后验分布、TSB-DPM模型截棍参数υ的条件后验分布，以及LVSVM分类器系数的条件后验分布、雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的条件后验分布；
步骤4，设定聚类高斯分布参数的初始值、雷达高分辨距离像的功率谱特征的聚类标号Z的初始值、TSB-DPM模型的截棍参数υ的初始值、LVSVM分类器系数的初始值以及雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的初始值；
在设定初始值之后，根据设定初始值的参数在步骤3中对应的条件后验分布，按照Gibbs采样技术对设定初始值的参数依次进行采样，总共对设定初始值的参数循环采样I0次，I0为自然数；
步骤5，在对设定初始值的参数循环采样I0次之后，从第I0+1次开始每间隔Sp次保存聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ，以及LVSVM分类器系数总共保存T0次参数的采样结果；
在保存T0次参数的采样结果后完成高分辨距离像数据HRRP的样本训练阶段，同时得到训练的LVSVM分类器和训练的TSB-DPM模型；
步骤6，对测试雷达高分辨距离像进行特征提取得到测试雷达高分辨距离像的功率谱特征计算测试雷达高分辨距离像的功率谱特征的概率密度函数值，并预先设定拒判门限Th，再与预先设定的拒判门限Th比较，根据比较结果判断测试雷达高分辨距离像是否为库外样本；否则继续步骤7；
步骤7，将保存的T0次参数的采样结果中的聚类的高斯分布参数截棍参数υ代入测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布，得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布；
在得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布之后，从测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布中采样得到测试雷达高分辨距离像的功率谱特征的聚类标号
步骤8，根据测试雷达高分辨距离像的功率谱特征所属的聚类标号将测试雷达高分辨距离像的功率谱特征依次输入其所属聚类对应的M个训练的LVSVM分类器中，即，将步骤7得到的测试雷达高分辨距离像的功率谱特征所属的聚类标号以及步骤5中保存的LVSVM分类器的系数代入到训练的LVSVM分类器的判别公式中，输出测试雷达高分辨距离像的目标类别标号
上述技术方案的特点和进一步改进在于：
(1)步骤2包括以下子步骤：
2a)利用TSB-DPM模型将功率谱特征集X进行聚类，包括以下2a1)、2a2)和2a3)：
2a1)在TSB-DPM模型中设定功率谱特征集X的最大聚类个数为C，每个聚类中雷达高分辨距离像的功率谱特征服从高斯分布；
2a2)设定TSB-DPM模型中的基分布G0采用Normal-Wishart分布其中，μ表示高斯分布的均值，Σ表示高斯分布的协方差矩阵，μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；
2a3)将以上2a1)和2a2)中的设定代入TSB-DPM模型得到下式(1-a)；
2b)利用LVSVM分类器对每个聚类中的雷达高分辨距离像的功率谱特征进行分类，包括以下2b1)、2b2)和2b3)：
2b1)设定每个LVSVM分类器系数先验分布为高斯分布表示高斯分布，I表示单位矩阵；
2b2)根据TSB-DPM模型中的功率谱特征集X的最大聚类个数为C个聚类且目标类别个数为M，采用一对多策略，即分别将M个类别中的一类目标看作正类目标，其它类别看作负类目标，分别训练LVSVM分类器，则需训练C×M个LVSVM分类器；
2b3)将LVSVM分类器系数ωcm的先验分布代入到训练的C×M个LVSVM分类器，得到以下式(1-b)；
2c)通过公式(1-a)和(1-b)共同构建dpLVSVM模型；
其中，υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；表示Normal-Wishart分布；{μc,Σc}表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵；G0表示基分布；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中雷达高分辨距离像的功率谱特征个数；π＝[π1,π2,...,πc...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；Mult(·)表示多项分布；ωc表示第c个聚类中所有M个LVSVM分类器的系数m＝1,2,...,M，M表示目标类别个数；ωcm表示第c个聚类中第m个LVSVM分类器的系数；λcm表示第c个聚类中的雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像功率谱特征对应的第m个LVSVM分类器的隐变量且有ym表示雷达高分辨距离像的功率谱特征对应于第m个LVSVM的类别标号；ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，且有：若雷达高分辨距离像的功率谱特征xn属于第m类目标则ynm＝+1，否则ynm＝-1；
表示第n个雷达高分辨距离像的功率谱特征的增广向量；γ表示调和系数；I表示单位矩阵；表示高斯分布；(·)T表示转置操作；
2d)由dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数和dpLVSVM模型各个参数的联合后验分布；dpLVSVM模型各个参数即TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布；
雷达高分辨距离像的功率谱特征的概率密度函数见以下公式(2)：
其中，π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数；表示均值为μc协方差矩阵为Σc的高斯分布，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；
dpLVSVM模型各个参数的联合后验分布见以下公式(3)：
其中，表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；表示雷达高分辨距离像的功率谱特征的聚类标号，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；表示LVSVM分类器系数，ωc表示第c个聚类中所有M个LVSVM分类器的系数m＝1,2,...,M，M表示目标类别个数；λ表示雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量；y表示雷达高分辨距离像的类别标号，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的标号；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；γ表示调和系数；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；I表示单位矩阵。
(2)步骤3包括以下子步骤：
3a)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类的高斯分布参数{μc,Σc}条件后验分布为：
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；X表示功率谱特征集；μ0为第c个聚类的高斯分布参数{μc,Σc}的先验分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；第c个聚类的高斯分布参数{μc,Σc}的条件后验分布为Normal-Wishart分布，其均值为尺度矩阵W为尺度因子β＝β0+Nc，尺度因子υ＝υ0+Nc，Nc表示功率谱特征集X中属于第c个聚类的雷达高分辨距离像的功率谱特征的个数；zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；
3b)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征的聚类标号zn的条件后验分布为：
zn～Mult(κn),κn＝[κn1,...,κnC]；
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；κnc表示第n个雷达高分辨距离像的功率谱特征属于第c个聚类的概率；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵；γ表示调和系数；ωcm表示第c个聚类中第m个LVSVM分类器的系数；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，m＝1,2,...,M，M表示目标类别个数；表示第n个雷达高分辨距离像的功率谱特征的增广向量，(·)T表示转置操作；Mult(·)表示多项分布；
3c)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类中第m个LVSVM分类器的系数ωcm条件后验分布为：
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；Z表示雷功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号；高斯分布的均值协方差矩阵为γ表示调和系数；表示第n个雷达高分辨距离像的功率谱特征的增广向量；(·)T表示转置操作；
3d)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量λnm条件后验分布为：
其中，m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；表示第zn个聚类中第m个LVSVM分类器的系数，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号；ynm表示第n个雷达高分辨距离像的功率谱特征xn对应于第m个LVSVM分类器的标号；表示第n个雷达高分辨距离像的功率谱特征的增广向量；表示逆高斯分布；
3e)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到STB-DPM模型的截棍参数υ中第c个变量vc的条件后验分布为：
p(vc|Z,α)＝Beta(vc；a,b) (7)
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Z表示功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；a＝1+Nc，Nc表示属于第c个聚类的雷达高分辨距离像的功率谱特征的个数，Nk表示属于第k个聚类的雷达高分辨距离像的功率谱特征的个数，k＝c+1,c+2,...,C，α表示TSB-DPM模型的截棍参数υ的先验分布的参数。
(3)步骤6包括以下子步骤：
6a)将保存的采样结果中的聚类的高斯分布参数{μc,Σc}和截棍参数υ代入步骤2得到的雷达高分辨距离像的功率谱特征的概率密度函数公式(2)计算测试雷达高分辨距离像的功率谱特征的概率密度函数值；
6b)给定预先设定的拒判门限Th；将测试雷达高分辨距离像的功率谱特征x·的概率密度函数值与拒判门限Th比较，判断测试雷达高分辨距离像是否为库外样本；
6c)根据保存T0次参数的采样结果得到测试雷达高分辨距离像的T0个判断结果；对T0个判断结果采用投票法则，即采用出现比例大于等于50％的判断结果，判断测试雷达高分辨距离像是否为库外样本；若为库外样本则将测试雷达高分辨距离像拒判，即不给定目标类别号并结束测试阶段；否则继续步骤7。
(4)步骤7中测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布公式见公式(9)：
其中，t表示保存的第t次采样，t＝1,2,...,T0，T0表示步骤5中设定的保存参数采样的个数；表示根据保存的第t次采样参数所确定的测试雷达高分辨距离像的功率谱特征属于第c个聚类的概率；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；表示根据保存的第t次采样参数所确定的第c个聚类的权系数，j＝1,2,...,c-1；{μc,Σc}t表示保存的第t次采样的第c个聚类的均值和协方差矩阵；c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；(·)T表示转置操作；Mult(·)表示多项分布。
(5)步骤8中LVSVM分类器的判别公式为下式(10)：
其中，表示第zt个聚类所对应的第m个LVSVM分类器的系数，m＝1,2,...,M，M表示目标类别个数，zt表示步骤5中得到的测试雷达高分辨距离像的功率谱特征的聚类标号，t＝1,2,...,T0，T0表示步骤5中设定的保存参数采样的个数；表示测试雷达高分辨距离像的功率谱特征的增广向量；ρm表示第m个LVSVM分类器的平均输出；表示求解最大值对应的m的值，(·)T表示转置操作。
本发明中提出dpLVSVM是一种无限混合专家模型，它是通过LVSVM(Latentvariable SVM，将隐变量SVM)分类器，见[Polson N.G.,Scott S.L..Data augmentationfor support vector machines[J].Bayesian Analysis,2011,vol.6(1),1-24]，引入TSB-DPM模型(截断Stick-breaking构造的Dirichlet process mixture model)，见[BleiD.M.and Jordan M.I..Variational inference for Dirichlet process mixtures[J].Bayesian Analysis,2006,vol.1(1),121-144]，来实现。
本发明与现有方法相比具有以下优点：
(1)与单分类器方法相比，本发明将数据划分为若干个聚类，可以通过多个简单的分类器实现复杂的全局分类，由于每个聚类中样本个数较少从而降分类器设计复杂度。
(2)与有限混合专家模型相比，本发明采用TSB-DPM模型可以自动选择数据的聚类个数并能将TSB-DPM模型与LVSVM分类器联合求解，能够保证每个聚类中的样本都具有良好的可分性，从而取得了更好的识别性能。
(3)本发明采用TSB-DPM模型，可以得到雷达高分辨距离像的功率谱特征的概率密度函数，由此可以描述数据的整体分布。而且通过根据雷达高分辨距离像的功率谱特征的概率密度值与预先设定的拒判门限进行比较从而实现对库外目标的拒判。
(4)与现有方法相比，本发明采用LVSVM作为分类器，可以通过Gibbs采样算法对参数进行估计，大大简化了求解复杂度。
本发明dpLVSVM模型利用TSB-DPM模型来自动将数据划分为多个具有高斯分布的聚类且不需事先确定样本聚类个数；同时在每个子集上训练一个形式简单的线性LVSVM分类器。由于该模型将聚类过程和分类器的训练过程进行联合优化，在一定程度上保证了各个聚类在分布上一致而且有一定可分性。dpLVSVM模型通过对数据潜在结构的挖掘，将非线性分类问题分解为多个线性可分的子问题，从而实现对整个数据的非线性分类并提高识别性能。本发明将LVSVM与DPM模型统一在一个框架下，可以采用Gibbs采样技术对参数进行简单有效的估计。另外，当观测目标不属于模板库内的任一目标类别时，需要能够对该库外目标进行拒判。dpLVSVM通过采用TSB-DPM模型对对数据进行描述，可以实现对库外样本的拒判。可用于处理大规模多模分布数据，将非线性分类问题分解为多个线性可分的子问题，从而实现对整个数据的非线性分类。
附图说明
下面结合附图和具体实施方式对本发明做进一步说明。
图1是基于本发明的目标识别算法流程图；
图2是本发明与现有技术的三种方法在不同特征维数下三类飞机的识别结果图；
图3是本发明与现有技术的三种方法的ROC曲线对比图。
具体实施方式
参照图1，说明本发明的一种基于dpLVSVM模型的雷达HRRP的目标识别方法，其具体步骤如下：
图1给出了整个识别***的流程，可以看出整个***包括两部分：训练阶段(左边部分)和测试阶段(右边部分)。其中，训练阶段的任务是对dpLVSVM模型进行参数估计，在训练阶段之后，测试阶段的任务是先执行拒判任务，然后根据训练得到参数计算样本所属聚类，最后输出目标的类别标号从而完成识别任务。
步骤1，雷达接收M个类别的目标的高分辨距离像HRRP；再对每个高分辨距离像进行特征提取，得到雷达高分辨距离像的功率谱特征xn，将M个类别的目标的高分辨距离像的功率谱特征组成功率谱特征集X；不属于该M个类别的目标的雷达高分辨距离像为库外样本；属于该M个类别的目标的雷达高分辨距离像为库内样本。
步骤2，利用功率谱特征集X，将LVSVM分类器和TSB-DPM模型结合构建dpLVSVM模型；根据dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数，以及dpLVSVM模型各个参数的联合条件后验分布；
dpLVSVM模型各个参数的联合条件后验分布即：TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布。
步骤2包括以下子步骤：
2a)利用TSB-DPM模型将功率谱特征集X进行聚类，包括以下2a1)、2a2)和2a3)：
2a1)在TSB-DPM模型中设定功率谱特征集X的最大聚类个数为C，每个聚类中雷达高分辨距离像的功率谱特征服从高斯分布；
2a2)设定TSB-DPM模型中的基分布G0采用Normal-Wishart分布其中，μ表示高斯分布的均值，Σ表示高斯分布的协方差矩阵，μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；
2a3)将以上2a1)和2a2)中的设定代入TSB-DPM模型得到下式(1-a)；
2b)利用LVSVM分类器对每个聚类中的雷达高分辨距离像的功率谱特征进行分类，包括以下2b1)、2b2)和2b3)：
2b1)设定每个LVSVM分类器系数先验分布为高斯分布表示高斯分布，I表示单位矩阵；
2b2)根据TSB-DPM模型中的功率谱特征集X的最大聚类个数为C个聚类且目标类别个数为M，采用一对多策略，即分别将M个类别中的一类目标看作正类目标，其它类别看作负类目标，分别训练LVSVM分类器，则需训练C×M个LVSVM分类器；
2b3)将LVSVM分类器系数ωcm的先验分布代入到训练的C×M个LVSVM分类器，得到以下式(1-b)；
2c)通过公式(1-a)和(1-b)共同构建dpLVSVM模型；
其中，υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；表示Normal-Wishart分布；{μc,Σc}表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵；G0表示基分布；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中雷达高分辨距离像的功率谱特征个数；π＝[π1,π2,...,πc...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；Mult(·)表示多项分布；ωc表示第c个聚类中所有M个LVSVM分类器的系数m＝1,2,...,M，M表示目标类别个数；ωcm表示第c个聚类中第m个LVSVM分类器的系数；λcm表示第c个聚类中的雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像功率谱特征对应的第m个LVSVM分类器的隐变量且有ym表示雷达高分辨距离像的功率谱特征对应于第m个LVSVM的类别标号；ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，且有：若雷达高分辨距离像的功率谱特征xn属于第m类目标则ynm＝+1，否则ynm＝-1；
表示第n个雷达高分辨距离像的功率谱特征的增广向量；γ表示调和系数；I表示单位矩阵；表示高斯分布；(·)T表示转置操作。
2d)由dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数和dpLVSVM模型各个参数的联合后验分布；dpLVSVM模型各个参数即TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布；
雷达高分辨距离像的功率谱特征的概率密度函数见以下公式(2)：
其中，π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数；表示均值为μc协方差矩阵为Σc的高斯分布，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数。
dpLVSVM模型各个参数的联合后验分布见以下公式(3)：
其中，表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；表示雷达高分辨距离像的功率谱特征的聚类标号，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；表示LVSVM分类器系数，ωc表示第c个聚类中所有M个LVSVM分类器的系数m＝1,2,...,M，M表示目标类别个数；λ表示雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量；y表示雷达高分辨距离像的类别标号，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的标号；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；γ表示调和系数；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；I表示单位矩阵。
步骤3，通过步骤2的dpLVSVM模型各个参数的联合条件后验分布，推导各个参数的条件后验分布，即聚类的高斯分布参数的条件后验分布、雷达高分辨距离像的功率谱特征的聚类标号Z的条件后验分布、TSB-DPM模型截棍参数υ的条件后验分布，以及LVSVM分类器系数的条件后验分布、雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的条件后验分布。
步骤3包括以下子步骤：
3a)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类的高斯分布参数{μc,Σc}条件后验分布为：
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；X表示功率谱特征集；μ0为第c个聚类的高斯分布参数{μc,Σc}的先验分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；第c个聚类的高斯分布参数{μc,Σc}的条件后验分布为Normal-Wishart分布，其均值为尺度矩阵W为尺度因子β＝β0+Nc，尺度因子υ＝υ0+Nc，Nc表示功率谱特征集X中属于第c个聚类的雷达高分辨距离像的功率谱特征的个数；zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数。
3b)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征的聚类标号zn的条件后验分布为：
zn～Mult(κn),κn＝[κn1,...,κnC]；
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；κnc表示第n个雷达高分辨距离像的功率谱特征属于第c个聚类的概率；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有j＝1,2,...,c-1；μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵；γ表示调和系数；ωcm表示第c个聚类中第m个LVSVM分类器的系数；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，m＝1,2,...,M，M表示目标类别个数；表示第n个雷达高分辨距离像的功率谱特征的增广向量，(·)T表示转置操作；Mult(·)表示多项分布。
3c)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类中第m个LVSVM分类器的系数ωcm条件后验分布为：
其中，c＝1,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；Z表示雷功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号；高斯分布的均值协方差矩阵为γ表示调和系数；表示第n个雷达高分辨距离像的功率谱特征的增广向量；(·)T表示转置操作。
3d)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量λnm条件后验分布为：
其中，m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；表示第zn个聚类中第m个LVSVM分类器的系数，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号；ynm表示第n个雷达高分辨距离像的功率谱特征xn对应于第m个LVSVM分类器的标号；表示第n个雷达高分辨距离像的功率谱特征的增广向量；表示逆高斯分布。
3e)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到STB-DPM模型的截棍参数υ中第c个变量vc的条件后验分布为：
p(vc|Z,α)＝Beta(vc；a,b) (16)
其中c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Z表示功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；a＝1+Nc，Nc表示属于第c个聚类的雷达高分辨距离像的功率谱特征的个数，Nk表示属于第k个聚类的雷达高分辨距离像的功率谱特征的个数，k＝c+1,c+2,...,C，α表示TSB-DPM模型的截棍参数υ的先验分布的参数。
步骤4，设定聚类高斯分布参数的初始值、雷达高分辨距离像的功率谱特征的聚类标号Z的初始值、TSB-DPM模型的截棍参数υ的初始值、LVSVM分类器系数的初始值以及雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的初始值；
在设定初始值之后，根据设定初始值的参数在步骤3中对应的条件后验分布，按照Gibbs采样技术对设定初始值的参数依次进行采样，总共对设定初始值的参数循环采样I0次，I0为自然数。
Gibbs采样技术可见于[Casella G.,George E.I..Explaining the Gibbssampler[J].The American Statisticain,1992,vol.46(3),167-174.]。
本发明中，使用Gibbs采样技术的原因为：dpLVSVM模型采用LVSVM作为分类器，整个模型可以用概率框架进行描述，见公式(1-a)和(1-b)，因而可以通过Gibbs采样算法对参数进行估计，可以大大简化求解复杂度。
步骤5，在对设定初始值的参数循环采样I0次之后，从第I0+1次开始每间隔Sp次保存聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ，以及LVSVM分类器系数总共保存T0次参数的采样结果；
在保存T0次参数的采样结果后完成高分辨距离像数据HRRP的样本训练阶段，同时得到训练的LVSVM分类器和训练的TSB-DPM模型。
本发明中步骤1到步骤5完成训练阶段。在执行完成步骤5之后，进入以下的测试阶段(目标识别阶段)即判断测试雷达高分辨距离像的功率谱特征的功率谱特征的目标类别号
步骤6，对测试雷达高分辨距离像进行特征提取得到测试雷达高分辨距离像的功率谱特征计算测试雷达高分辨距离像的功率谱特征的概率密度函数值，并预先设定拒判门限Th，再与预先设定的拒判门限Th比较，根据比较结果判断测试雷达高分辨距离像是否为库外样本；否则继续步骤7。
步骤6包括以下子步骤：
6a)将保存的采样结果中的聚类的高斯分布参数{μc,Σc}和截棍参数υ代入步骤2得到的雷达高分辨距离像的功率谱特征的概率密度函数公式(2)计算测试雷达高分辨距离像的功率谱特征的概率密度函数值；
6b)给定预先设定的拒判门限Th；将测试雷达高分辨距离像的功率谱特征的概率密度函数值与拒判门限Th比较，判断测试雷达高分辨距离像是否为库外样本；
6c)根据保存T0次参数的采样结果得到测试雷达高分辨距离像的T0个判断结果；对T0个判断结果采用投票法则，即采用出现比例大于等于50％的判断结果，判断测试雷达高分辨距离像是否为库外样本；若为库外样本则将测试雷达高分辨距离像拒判，即不给定目标类别号并结束测试阶段；否则继续步骤7。
步骤7，将保存的T0次参数的采样结果中的聚类的高斯分布参数截棍参数υ代入测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布，得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布；
在得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布之后，从测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布中采样得到测试雷达高分辨距离像的功率谱特征的聚类标号
具体的，测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布公式见公式(9)：
其中，t表示保存的第t次采样，t＝1,2,...,T0，T0表示步骤5中设定的保存参数采样的个数；表示根据保存的第t次采样参数所确定的测试雷达高分辨距离像的功率谱特征属于第c个聚类的概率；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；表示根据保存的第t次采样参数所确定的第c个聚类的权系数，j＝1,2,...,c-1；{μc,Σc}t表示保存的第t次采样的第c个聚类的均值和协方差矩阵；c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；(·)T表示转置操作；Mult(·)表示多项分布。
步骤8，根据测试雷达高分辨距离像的功率谱特征所属的聚类标号将测试雷达高分辨距离像的功率谱特征依次输入其所属聚类对应的M个训练的LVSVM分类器中，即，将步骤7得到的测试雷达高分辨距离像的功率谱特征所属的聚类标号以及步骤5中保存的LVSVM分类器的系数代入到训练的LVSVM分类器的判别公式中，输出测试雷达高分辨距离像的目标类别标号
具体的，LVSVM分类器的判别公式为下式(10)：
其中，表示第zt个聚类所对应的第m个LVSVM分类器的系数，m＝1,2,...,M，M表示目标类别个数，zt表示步骤5中得到的测试雷达高分辨距离像的功率谱特征的聚类标号，t＝1,2,...,T0，T0表示步骤5中设定的保存参数采样的个数；表示测试雷达高分辨距离像的功率谱特征的增广向量；ρm表示第m个LVSVM分类器的平均输出；表示求解最大值对应的m的值，(·)T表示转置操作。
本发明中，dpLVSVM模型中采用TSB-DPM模型将功率谱特征集划分为至多C个聚类，同时在每个聚类上训练LVSVM分类器，可以通过多个简单的分类器实现复杂的全局分类。由于每个聚类中的样本数远小于总体样本数(也就是，功率谱特征集X中属于第c个聚类的功率谱特征的个数Nc远小于功率谱特征集X中功率谱特征个数N)，从而降底了LVSVM分类器的设计复杂度。
与有限混合专家模型相比，dpLVSVM由于采用了TSB-DPM模型，因而可以自动选择数据的聚类个数。C为最大的聚类个数，在实际中，dpLVSVM模型根据雷达高分辨距离像的实际分布自动确定聚类的个数，最终聚类个数要小于C。
不同于现有技术中的聚类和分类相互独立的方法，dpLVSVM构建了TSB-DPM模型聚类与LVSVM分类器分类的相关性，见公式(1-a)和(1-b)，即在对功率谱特征集进行聚类的同时，在每个聚类上训练分类器。将整个模型联合求解，从而能够保证每个聚类中的样本都具有良好的可分性。
根据dpLVSVM模型，推导出了雷达高分辨距离像的功率谱特征的概率密度函数，通过比较雷达高分辨距离像的概率密度值与预先设定的拒判门限实现对库外目标的拒判。
下面结合仿真实验对本发明的效果做进一步说明。
(1)实验条件
本实验采用维数较高且分布相对复杂的实测雷达HRRP数据。该数据为某院的Ｃ波段雷达实测飞机的一维HRRP数据。数据中包含三类飞机目标(雅-42、奖状、安-26)。雷达参数和三类飞机目标的参数如表1所示。
表1
三类飞机的HRRP数据均被划分成了若干段。分别选择“雅-42”的第2、5段，“奖状”的第6、7段以及“安-26”的第5、6段共600个样本作为训练数据集，选择其余段中2400个雷达高分辨距离像样本作为测试数据集。
预处理：采用幅度2范数归一的方法，对HRRP信号进行归一。然后提取功率谱特征。原始HRRP维数为256，由于功率谱具有称性，只需取128维作为特征。为了提高计算效率，采用PCA算法对数据进行降维，并比较了不同维数下各个分类器的识别性能。
dpLVSVM模型参数设置如下：γ＝1，W0＝1e-6Iq，β0＝0.01，υ0＝q，其中q为样本维数，α＝0.1，I＝1000，Sp＝10，T＝100。
(2)实验内容
(2a)为了进一步说明本发明dpLVSVM模型在识别性能上的优势，与下面现有技术中的三种模型进行对比：线性SVM(LSVM)，Km+SVM，dp+SVM。其中Km+SVM表示先采用K-means算法将训练样本进行聚类，然后每个聚类分别训练一个SVM分类器，两个过程是分离的而且其聚类个数同过交叉验证的方法确定；dp+SVM表示先采用DPM模型对样本进行聚类，然后每个聚类分别训练一个SVM分类器，两个过程是分离的，其聚类个数不用事先确定。
(2b)选取了4类其它飞机目标作为库外目标，每个目标等间隔抽取200个样本(共800个样本)作为库外目标样本。实验中比较了现有技术中的SVDD、Km+SVM、dp+SVM以及本发明dpLVSVM模型四种方法的拒判性能。
(3)实验结果分析
图2给出了不同方法在不同特征维数下的识别结果，其中图2横坐标为雷达高分辨距离像功率谱特征的维数，纵坐标为识别率。参照图2，本发明dpLVSVM模型在各个特征维度下的性能优于现有技术中的三种模型(LSVM，Km+SVM，dp+SVM)，特别是当特征维数为15时平均正确识别率达到了最高的0.930。图2同时表明维数对识别率产生一定的影响：当特征维数较小时由于损失了较多的信息，识别率较低；特征维数较大时，特征中会包含一定的冗余信息，对识别有一定干扰作用，识别率有所降低。
分类器的拒判性能通常由接收机工作特性(Receiver operatingcharacteristic，ROC)曲线来衡量。ROC曲线的横轴为虚警概率，纵轴为检测概率。ROC曲线下的面积AUC(Area under an ROC curve)越大，说明分类器的拒判性能越好。表2示出本发明dpLVSVM模型和现有技术中的SVDD、Km+SVM、dp+SVM方法的AUC值比较结果。图3为本发明dpLVSVM模型和现有技术中的SVDD、Km+SVM、dp+SVM方法的ROC曲线图，其中横坐标为虚警概率纵坐标为检测概率。由图3和表2可知，采用DPM模型的分类模型(dp+SVM以及dpLVSVM)均可以对数据的分布进行较好的描述，其拒判性能要强于现有技术的SVDD方法以及K-means方法。
表2
方法 SVDD Km+SVM dp+SVM dpLVSVM
AUC值 0.630 0.648 0.881 0.882
综合识别和拒判结果可见，本发明dpLVSVM模型既能提高分类性能又有良好的拒判性能。

Claims (6)

1.一种基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，包括以下步骤：
步骤1，雷达接收M个类别的目标的高分辨距离像HRRP；再对每个高分辨距离像进行特征提取，得到雷达高分辨距离像的功率谱特征xn，将M个类别的目标的高分辨距离像的功率谱特征组成功率谱特征集X；不属于该M个类别的目标的雷达高分辨距离像为库外样本；属于该M个类别的目标的雷达高分辨距离像为库内样本；
步骤2，利用功率谱特征集X，将LVSVM分类器和TSB-DPM模型结合构建dpLVSVM模型；即通过TSB-DPM模型公式(1-a)和LVSVM分类器公式(1-b)共同构建dpLVSVM模型为：
其中，υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；表示Normal-Wishart分布；{μc,Σc}表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，Σc表示第c个聚类的协方差矩阵；G0表示基分布；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中雷达高分辨距离像的功率谱特征个数；π＝[π1,π2,...,πc...,πC]表示每个聚类的权系数且有Mult(·)表示多项分布；wc表示第c个聚类中所有M个LVSVM分类器的系数M表示目标类别个数；wcm表示第c个聚类中第m个LVSVM分类器的系数；λcm表示第c个聚类中的雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像功率谱特征对应的第m个LVSVM分类器的隐变量且有ym表示雷达高 分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号；ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，且有：若雷达高分辨距离像的功率谱特征xn属于第m类目标则ynm＝+1，否则ynm＝-1；
表示第n个雷达高分辨距离像的功率谱特征的增广向量；γ表示调和系数；I表示单位矩阵；
表示高斯分布；(·)T表示转置操作；
根据dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数，以及dpLVSVM模型各个参数的联合条件后验分布；
dpLVSVM模型各个参数的联合条件后验分布即：TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布；
步骤3，通过步骤2的dpLVSVM模型各个参数的联合条件后验分布，推导各个参数的条件后验分布，即聚类的高斯分布参数的条件后验分布、雷达高分辨距离像的功率谱特征的聚类标号Z的条件后验分布、TSB-DPM模型截棍参数υ的条件后验分布，以及LVSVM分类器系数的条件后验分布、雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的条件后验分布；
步骤4，设定聚类高斯分布参数的初始值、雷达高分辨距离像的功率谱特征的聚类标号Z的初始值、TSB-DPM模型的截棍参数υ的初始值、LVSVM分类器系数的初始值以及雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的初始值；
在设定初始值之后，根据设定初始值的参数在步骤3中对应的条件后验分布，按照Gibbs采样技术对设定初始值的参数依次进行采样，总共对设定初始值的参数循环采样I0次，I0为自然数；
步骤5，在对设定初始值的参数循环采样I0次之后，从第I0+1次开始每间隔Sp次保存 聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ，以及LVSVM分类器系数总共保存T0次参数的采样结果；
在保存T0次参数的采样结果后完成高分辨距离像HRRP的样本训练阶段，同时得到训练的LVSVM分类器和训练的TSB-DPM模型；
步骤6，对测试雷达高分辨距离像进行特征提取得到测试雷达高分辨距离像的功率谱特征计算测试雷达高分辨距离像的功率谱特征的概率密度函数值，并预先设定拒判门限Th，将所述测试雷达高分辨距离像的功率谱特征的概率密度函数值与预先设定的拒判门限Th比较，根据比较结果判断测试雷达高分辨距离像是否为库外样本；若为库内样本则继续步骤7；
步骤7，将保存的T0次参数的采样结果中的聚类的高斯分布参数截棍参数υ代入测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布，得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布，为根据保存的第t次采样的参数得到的测试雷达高分辨距离像的功率谱特征的聚类标号，T0表示步骤5中设定的保存参数的采样结果的次数；
在得到测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布之后，从测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布中采样得到测试雷达高分辨距离像的功率谱特征所属的聚类标号
步骤8，根据测试雷达高分辨距离像的功率谱特征所属的聚类标号将测试雷达高分辨距离像的功率谱特征依次输入其所属聚类对应的M个训练的LVSVM分类器中，即，将步骤7得到的测试雷达高分辨距离像的功率谱特征所属的聚类标号以及步骤5中保存的LVSVM分类器的系数代入到训练的LVSVM分类器的判别公式中，输出测试雷达高分辨距离像的目标类别标号
2.根据权利要求1所述的基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，步骤2包括以下子步骤：
2a)利用TSB-DPM模型将功率谱特征集X进行聚类，包括以下步骤2a1)、2a2)和 2a3)：
2a1)在TSB-DPM模型中设定功率谱特征集X的最大聚类个数为C，每个聚类中雷达高分辨距离像的功率谱特征服从高斯分布；
2a2)设定TSB-DPM模型中的基分布G0采用Normal-Wishart分布其中，μ表示高斯分布的均值，Σ表示高斯分布的协方差矩阵，μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；
2a3)将以上步骤2a1)和2a2)中的设定代入TSB-DPM模型得到所述公式(1-a)；
2b)利用LVSVM分类器对每个聚类中的雷达高分辨距离像的功率谱特征进行分类，包括以下步骤2b1)、2b2)和2b3)：
2b1)设定每个LVSVM分类器系数先验分布为高斯分布 表示高斯分布，I表示单位矩阵；
2b2)根据TSB-DPM模型中的功率谱特征集X的最大聚类个数为C个聚类且目标类别个数为M，采用一对多策略，即分别将M个类别中的一类目标看作正类目标，其它类别看作负类目标，分别训练LVSVM分类器，则需训练C×M个LVSVM分类器；
2b3)将LVSVM分类器系数wcm的先验分布代入到训练的C×M个LVSVM分类器，得到所述公式(1-b)；
2c)通过公式(1-a)和(1-b)共同构建所述dpLVSVM模型；
2d)由dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数和dpLVSVM模型各个参数的联合后验分布；dpLVSVM模型各个参数即TSB-DPM模型中聚类的高斯分布参数雷达高分辨距离像的功率谱特征的聚类标号Z、TSB-DPM模型的截棍参数υ、LVSVM分类器系数雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ；
雷达高分辨距离像的功率谱特征的概率密度函数见以下公式(2)：
其中，π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数；表示均值为μc协方差矩阵为∑c的高斯分布，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚 类个数；
dpLVSVM模型各个参数的联合后验分布见以下公式(3)：
其中，表示第c个聚类的高斯分布参数，μc表示第c个聚类的均值，∑c表示第c个聚类的协方差矩阵，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；表示雷达高分辨距离像的功率谱特征的聚类标号，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有 表示LVSVM分类器系数，wc表示第c个聚类中所有M个LVSVM分类器的系数M表示目标类别个数；λ表示雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量，λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量；y表示雷达高分辨距离像的类别标号，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号；Beta(·)表示Beta分布；α表示TSB-DPM模型的截棍参数υ的先验分布的参数；γ表示调和系数；μ0为Normal-Wishart分布的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；I表示单位矩阵。
3.根据权利要求2所述的基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，步骤3包括以下子步骤：
3a)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类的高斯分布参数{μc,∑c}条件后验分布为：
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；X表示功率谱特征集；μc为第c个聚类的均值，W0为尺度矩阵，β0、υ0为两个尺度因子；第c个聚类的高斯分布参数{μc,∑c}的条件后验分布为Normal-Wishart分布，其均值为尺度矩阵W为尺度因子β＝β0+Nc，尺度因子υ＝υ0+Nc，Nc表示功率谱特征集X中属于第c个聚类的雷达高分辨距离像的功率谱特征的个数；zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号，n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；
3b)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征的聚类标号zn的条件后验分布为：
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；κnc表示第n个雷达高分辨距离像的功率谱特征属于第c个聚类的概率；π＝[π1,π2,...,πc,...,πC]表示每个聚类的权系数且有μc表示第c个聚类的均值，∑c表示第c个聚类的协方差矩阵；γ表示调和系数；wcm表示第c个聚类中第m个LVSVM分类器的系数；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号，m＝1,2,...,M，M表示目标类别个数；表示第n个雷达高分辨距离像的功率谱特征的增广向量， (·)T表示转置操作；Mult(·)表示多项分布；
3c)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第c个聚类中第m个LVSVM分类器的系数wcm条件后验分布为：
其中，c＝1,2,...C,，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；Z表示功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；zn表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号；λnm表示第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量，ynm表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号；高斯分布的均值 协方差矩阵为γ表示调和系数； 表示第n个雷达高分辨距离像的功率谱特征的增广向量；
3d)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到第n个雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量λnm条件后验分布为：
其中，m＝1,2,...,M，M表示目标类别个数；n＝1,2,...,N，N表示功率谱特征集X中功率谱特征个数；wznm表示第zn个聚类中第m个LVSVM分类器的系数，zn表示第n个雷达高分辨距离像的功率谱特征的聚类标号；ynm表示第n个雷达高分辨距离像的功率谱特征xn对应于第m个LVSVM分类器的类别标号；表示第n个雷达高分辨距离像的功率谱特征的增广向量；表示逆高斯分布；
3e)根据贝叶斯公式和dpLVSVM模型各个参数的联合条件后验分布，得到TSB-DPM模型的截棍参数υ中第c个变量vc的条件后验分布为：
p(vc|Z,α)＝Beta(vc；a,b) (8)
其中，c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；Z表示功率谱特征集X中雷达高分辨距离像的功率谱特征的聚类标号；
a＝1+Nc，Nc表示属于第c个聚类的雷达高分辨距离像的功率谱特征的个数，Nk表示属于第k个聚类的雷达高分辨距离像的功率谱特征的个数，k＝c+1,c+2,...,C，α表示TSB-DPM模型的截棍参数υ的先验分布的参数。
4.根据权利要求2所述的基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，步骤6包括以下子步骤：
6a)将保存的采样结果中的聚类的高斯分布参数{μc,Σc}和截棍参数υ代入步骤2得到的雷达高分辨距离像的功率谱特征的概率密度函数公式(2)计算测试雷达高分辨距离像的功率谱特征的概率密度函数值；
6b)给定预先设定的拒判门限Th；将测试雷达高分辨距离像的功率谱特征的概率密度函数值与拒判门限Th比较，判断测试雷达高分辨距离像是否为库外样本；
6c)根据保存T0次参数的采样结果得到测试雷达高分辨距离像的T0个判断结果；对T0个判断结果采用投票法则，即采用出现比例大于等于50％的判断结果，判断测试雷达高分辨距离像是否为库外样本；若为库外样本则将测试雷达高分辨距离像拒判，即不给定目标类别号并结束测试阶段；否则继续步骤7。
5.根据权利要求1所述的基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，
步骤7中测试雷达高分辨距离像的功率谱特征的聚类标号的条件后验分布公式见公式(9)：
其中，表示高斯分布，t表示保存的第t次采样，t＝1,2,...,T0，T0表示步骤5中设定的保存参数的采样结果的次数；表示根据保存的第t次采样参数所确定的测试雷达高分辨距离像的功率谱特征属于第c个聚类的概率；υ＝[v1,v2,...,vc,...,vC]表示TSB-DPM模型的截棍参数；表示根据保存的第t次采样参数所确定的第c个聚类的权系数，{μc,∑c}t表示保存的第t次采样的第c个聚类的均值和协方差矩阵；c＝1,2,...,C，C为TSB-DPM模型设定的功率谱特征集X的最大聚类个数；(·)T表示转置操作；Mult(·)表示多项分布。
6.根据权利要求1所述的基于dpLVSVM模型的雷达HRRP的目标识别方法，其特征在于，
步骤8中LVSVM分类器的判别公式为下式(10)：
其中，表示第zt个聚类所对应的第m个LVSVM分类器的系数，m＝1,2,...,M，M表示目标类别个数，zt表示步骤7中得到的测试雷达高分辨距离像的功率谱特征所属的聚类标号，t＝1,2,...,T0，T0表示步骤5中设定的保存参数的采样结果的次数；表示测试雷达高分辨距离像的功率谱特征的增广向量；ρm表示第m个LVSVM分类器的平均输出；表示求解最大值对应的m的值，(·)T表示转置操作。
CN201410234677.1A 2014-05-29 2014-05-29 基于dpLVSVM模型的雷达HRRP的目标识别方法 Active CN104007431B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410234677.1A CN104007431B (zh) 2014-05-29 2014-05-29 基于dpLVSVM模型的雷达HRRP的目标识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410234677.1A CN104007431B (zh) 2014-05-29 2014-05-29 基于dpLVSVM模型的雷达HRRP的目标识别方法

Publications (2)

Publication Number Publication Date
CN104007431A CN104007431A (zh) 2014-08-27
CN104007431B true CN104007431B (zh) 2017-03-29

Family

ID=51368175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410234677.1A Active CN104007431B (zh) 2014-05-29 2014-05-29 基于dpLVSVM模型的雷达HRRP的目标识别方法

Country Status (1)

Country Link
CN (1) CN104007431B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016080156A1 (ja) * 2014-11-19 2016-05-26 三菱電機株式会社 レーダ装置
CN105116397B (zh) * 2015-08-25 2017-09-26 西安电子科技大学 基于mmfa模型的雷达高分辨率距离像目标识别方法
CN105116400B (zh) * 2015-08-31 2018-04-10 西安电子科技大学 基于iMMFA模型的雷达高分辨距离像的目标识别方法
CN105334504B (zh) * 2015-09-08 2017-11-21 西安电子科技大学 基于大边界的非线性判别投影模型的雷达目标识别方法
CN105223552B (zh) * 2015-09-15 2018-03-06 上海无线电设备研究所 基于无线电谱纹识别的干扰辨识方法
CN105975994B (zh) * 2016-05-19 2019-05-21 西安电子科技大学 基于非相似性变换一类svm模型的sar目标鉴别方法
CN106951915B (zh) * 2017-02-23 2020-02-21 南京航空航天大学 一种基于类别置信度的一维距离像多分类器融合识别法
CN107219510B (zh) * 2017-05-18 2019-11-01 西安电子科技大学 基于无限最大间隔线性判别投影模型的雷达目标识别方法
CN109409117B (zh) * 2017-08-15 2021-10-22 创新先进技术有限公司 敏感数据的差分隐私保护方法和装置
CN107576948B (zh) * 2017-08-15 2020-09-25 电子科技大学 一种基于高分辨距离像imf特征的雷达目标识别方法
CN108416393A (zh) * 2018-03-16 2018-08-17 西安电子科技大学 一种基于高分辨一维距离像的目标识别优化方法
CN111208483B (zh) * 2020-01-03 2022-04-19 西安电子科技大学 基于贝叶斯支持向量数据描述的雷达库外目标识别方法
CN111401470B (zh) * 2020-03-31 2023-03-10 西安电子科技大学 基于特征空间分布的Fisher特征选择方法
CN112198488B (zh) * 2020-09-08 2021-07-20 北京理工大学 一种极窄脉冲雷达角度关联目标回波空间聚合方法
CN113361472B (zh) * 2021-07-01 2023-04-18 西安电子科技大学 基于ILFAcs模型的雷达HRRP目标识别方法

Also Published As

Publication number Publication date
CN104007431A (zh) 2014-08-27

Similar Documents

Publication Publication Date Title
CN104007431B (zh) 基于dpLVSVM模型的雷达HRRP的目标识别方法
CN105116397B (zh) 基于mmfa模型的雷达高分辨率距离像目标识别方法
CN113378632B (zh) 一种基于伪标签优化的无监督域适应行人重识别方法
CN109190524B (zh) 一种基于生成对抗网络的人体动作识别方法
Chang et al. Learning representations of emotional speech with deep convolutional generative adversarial networks
CN105116400B (zh) 基于iMMFA模型的雷达高分辨距离像的目标识别方法
CN109117826B (zh) 一种多特征融合的车辆识别方法
CN107507197B (zh) 一种基于聚类算法和卷积神经网络的肺实质提取方法
Cateni et al. A hybrid feature selection method for classification purposes
CN111000553B (zh) 一种基于投票集成学习的心电数据智能分类方法
CN106951915B (zh) 一种基于类别置信度的一维距离像多分类器融合识别法
CN106980848A (zh) 基于曲波变换和稀疏学习的人脸表情识别方法
CN105389583A (zh) 图像分类器的生成方法、图像分类方法和装置
CN109359525B (zh) 基于稀疏低秩的判别谱聚类的极化sar图像分类方法
Wan et al. Recognizing the HRRP by combining CNN and BiRNN with attention mechanism
CN103886335B (zh) 基于模糊粒子群和散射熵的极化sar图像分类方法
CN109840518B (zh) 一种结合分类与域适应的视觉追踪方法
CN107219510B (zh) 基于无限最大间隔线性判别投影模型的雷达目标识别方法
Bahari Speaker age estimation using Hidden Markov Model weight supervectors
CN111273288A (zh) 一种基于长短期记忆网络的雷达未知目标识别方法
CN106951822B (zh) 一种基于多尺度稀疏保持投影一维距离像融合识别方法
CN109871907B (zh) 基于sae-hmm模型的雷达目标高分辨率距离像识别方法
CN108805181B (zh) 一种基于多分类模型的图像分类装置及分类方法
Xin et al. Random part localization model for fine grained image classification
CN113239199A (zh) 一种基于多方数据集的信用分类方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant