发明名称 基于dpLVSVM模型的雷达HRRP的目标识别方法
摘要 本发明公开了一种基于dpLVSVM模型的雷达HRRP的目标识别方法,其步骤为:步骤1,对雷达HRRP数据进行特征提取得到功率谱特征集X;步骤2,构建dpLVSVM模型并得出功率谱特征的概率密度函数及各个参数的联合条件后验分布;步骤3,推导各个参数的条件后验分布;步骤4,对各个参数进行循环采样I次;步骤5,保存T<sub>0</sub>次测试阶段所需参数的采样结果;步骤6,判断测试雷达HRRP是否为库外样本,若为库外样本则拒判;否则步骤7;步骤7,采样得到测试雷达HRRP的功率谱特征<img file="DDA0000512982710000011.GIF" wi="42" he="63" />的聚类标号;步骤8,输出测试雷达HRRP的目标类别标号<img file="DDA0000512982710000012.GIF" wi="62" he="57" />本发明具有分类器设计复杂度小,识别性能高以及拒判性能好的优点,可用于对雷达目标的识别。
申请公布号 CN104007431B 申请公布日期 2017.03.29
申请号 CN201410234677.1 申请日期 2014.05.29
申请人 西安电子科技大学 发明人 陈渤;张学峰;陈步华;王鹏辉;刘宏伟
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于dpLVSVM模型的雷达HRRP的目标识别方法,其特征在于,包括以下步骤:步骤1,雷达接收M个类别的目标的高分辨距离像HRRP;再对每个高分辨距离像进行特征提取,得到雷达高分辨距离像的功率谱特征x<sub>n</sub>,将M个类别的目标的高分辨距离像的功率谱特征组成功率谱特征集X;不属于该M个类别的目标的雷达高分辨距离像为库外样本;属于该M个类别的目标的雷达高分辨距离像为库内样本;步骤2,利用功率谱特征集X,将LVSVM分类器和TSB‑DPM模型结合构建dpLVSVM模型;即通过TSB‑DPM模型公式(1‑a)和LVSVM分类器公式(1‑b)共同构建dpLVSVM模型为:<img file="1.GIF" wi="1662" he="324" /><img file="FDA0001114216730000012.GIF" wi="1645" he="83" />其中,υ=[v<sub>1</sub>,v<sub>2</sub>,...,v<sub>c</sub>,...,v<sub>C</sub>]表示TSB‑DPM模型的截棍参数,c=1,2,...,C,C为TSB‑DPM模型设定的功率谱特征集X的最大聚类个数;Beta(·)表示Beta分布;α表示TSB‑DPM模型的截棍参数υ的先验分布的参数;<img file="2.GIF" wi="150" he="60" />表示Normal‑Wishart分布;{μ<sub>c</sub>,Σ<sub>c</sub>}表示第c个聚类的高斯分布参数,μ<sub>c</sub>表示第c个聚类的均值,Σ<sub>c</sub>表示第c个聚类的协方差矩阵;G<sub>0</sub>表示基分布;μ<sub>0</sub>为Normal‑Wishart分布的均值,W<sub>0</sub>为尺度矩阵,β<sub>0</sub>、υ<sub>0</sub>为两个尺度因子;z<sub>n</sub>表示第n个雷达高分辨距离像的功率谱特征所属的聚类标号,n=1,2,...,N,N表示功率谱特征集X中雷达高分辨距离像的功率谱特征个数;π=[π<sub>1</sub>,π<sub>2</sub>,...,π<sub>c</sub>...,π<sub>C</sub>]表示每个聚类的权系数且有<img file="FDA0001114216730000013.GIF" wi="733" he="140" />Mult(·)表示多项分布;w<sub>c</sub>表示第c个聚类中所有M个LVSVM分类器的系数<img file="FDA0001114216730000014.GIF" wi="517" he="70" />M表示目标类别个数;w<sub>cm</sub>表示第c个聚类中第m个LVSVM分类器的系数;λ<sub>cm</sub>表示第c个聚类中的雷达高分辨距离像的功率谱特征对应的第m个LVSVM分类器的隐变量,λ<sub>nm</sub>表示第n个雷达高分辨距离像功率谱特征对应的第m个LVSVM分类器的隐变量且有<img file="FDA0001114216730000015.GIF" wi="363" he="87" />y<sub>m</sub>表示雷达高 分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号;y<sub>nm</sub>表示第n个雷达高分辨距离像的功率谱特征对应于第m个LVSVM分类器的类别标号,且有:若雷达高分辨距离像的功率谱特征x<sub>n</sub>属于第m类目标则y<sub>nm</sub>=+1,否则y<sub>nm</sub>=‑1;<img file="3.GIF" wi="1819" he="173" />表示第n个雷达高分辨距离像的功率谱特征的增广向量;γ表示调和系数;I表示单位矩阵;<img file="FDA0001114216730000029.GIF" wi="95" he="60" />表示高斯分布;(·)<sup>T</sup>表示转置操作;根据dpLVSVM模型推导出雷达高分辨距离像的功率谱特征的概率密度函数,以及dpLVSVM模型各个参数的联合条件后验分布;dpLVSVM模型各个参数的联合条件后验分布即:TSB‑DPM模型中聚类的高斯分布参数<img file="FDA0001114216730000023.GIF" wi="243" he="70" />雷达高分辨距离像的功率谱特征的聚类标号Z、TSB‑DPM模型的截棍参数υ、LVSVM分类器系数<img file="FDA0001114216730000024.GIF" wi="165" he="70" />雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的联合条件后验分布;步骤3,通过步骤2的dpLVSVM模型各个参数的联合条件后验分布,推导各个参数的条件后验分布,即聚类的高斯分布参数<img file="FDA0001114216730000025.GIF" wi="213" he="71" />的条件后验分布、雷达高分辨距离像的功率谱特征的聚类标号Z的条件后验分布、TSB‑DPM模型截棍参数υ的条件后验分布,以及LVSVM分类器系数<img file="FDA0001114216730000026.GIF" wi="147" he="68" />的条件后验分布、雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的条件后验分布;步骤4,设定聚类高斯分布参数<img file="FDA0001114216730000027.GIF" wi="206" he="71" />的初始值、雷达高分辨距离像的功率谱特征的聚类标号Z的初始值、TSB‑DPM模型的截棍参数υ的初始值、LVSVM分类器系数<img file="FDA0001114216730000028.GIF" wi="142" he="71" />的初始值以及雷达高分辨距离像的功率谱特征对应的LVSVM分类器的隐变量λ的初始值;在设定初始值之后,根据设定初始值的参数在步骤3中对应的条件后验分布,按照Gibbs采样技术对设定初始值的参数依次进行采样,总共对设定初始值的参数循环采样I<sub>0</sub>次,I<sub>0</sub>为自然数;步骤5,在对设定初始值的参数循环采样I<sub>0</sub>次之后,从第I<sub>0</sub>+1次开始每间隔S<sub>p</sub>次保存 聚类的高斯分布参数<img file="FDA0001114216730000031.GIF" wi="239" he="63" />雷达高分辨距离像的功率谱特征的聚类标号Z、TSB‑DPM模型的截棍参数υ,以及LVSVM分类器系数<img file="FDA0001114216730000032.GIF" wi="163" he="71" />总共保存T<sub>0</sub>次参数的采样结果;在保存T<sub>0</sub>次参数的采样结果后完成高分辨距离像HRRP的样本训练阶段,同时得到训练的LVSVM分类器和训练的TSB‑DPM模型;步骤6,对测试雷达高分辨距离像进行特征提取得到测试雷达高分辨距离像的功率谱特征<img file="FDA0001114216730000033.GIF" wi="57" he="51" />计算测试雷达高分辨距离像的功率谱特征<img file="FDA0001114216730000034.GIF" wi="36" he="47" />的概率密度函数值,并预先设定拒判门限T<sub>h</sub>,将所述测试雷达高分辨距离像的功率谱特征<img file="FDA0001114216730000035.GIF" wi="38" he="46" />的概率密度函数值与预先设定的拒判门限T<sub>h</sub>比较,根据比较结果判断测试雷达高分辨距离像是否为库外样本;若为库内样本则继续步骤7;步骤7,将保存的T<sub>0</sub>次参数的采样结果中的聚类的高斯分布参数<img file="FDA0001114216730000036.GIF" wi="243" he="63" />截棍参数υ代入测试雷达高分辨距离像的功率谱特征<img file="FDA0001114216730000037.GIF" wi="38" he="47" />的聚类标号<img file="FDA0001114216730000038.GIF" wi="33" he="46" />的条件后验分布,得到测试雷达高分辨距离像的功率谱特征<img file="FDA0001114216730000039.GIF" wi="35" he="54" />的聚类标号<img file="FDA00011142167300000310.GIF" wi="116" he="63" />的条件后验分布,<img file="FDA00011142167300000311.GIF" wi="38" he="47" />为根据保存的第t次采样的参数得到的测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000312.GIF" wi="36" he="54" />的聚类标号,T<sub>0</sub>表示步骤5中设定的保存参数的采样结果的次数;在得到测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000313.GIF" wi="43" he="45" />的聚类标号<img file="FDA00011142167300000314.GIF" wi="119" he="63" />的条件后验分布之后,从测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000315.GIF" wi="34" he="47" />的聚类标号<img file="FDA00011142167300000316.GIF" wi="117" he="63" />的条件后验分布中采样得到测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000317.GIF" wi="40" he="47" />所属的聚类标号<img file="FDA00011142167300000318.GIF" wi="141" he="63" />步骤8,根据测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000319.GIF" wi="35" he="47" />所属的聚类标号<img file="FDA00011142167300000320.GIF" wi="139" he="63" />将测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000321.GIF" wi="38" he="47" />依次输入其所属聚类对应的M个训练的LVSVM分类器中,即,将步骤7得到的测试雷达高分辨距离像的功率谱特征<img file="FDA00011142167300000322.GIF" wi="40" he="47" />所属的聚类标号<img file="FDA00011142167300000323.GIF" wi="118" he="63" />以及步骤5中保存的LVSVM分类器的系数<img file="FDA00011142167300000324.GIF" wi="139" he="63" />代入到训练的LVSVM分类器的判别公式中,输出测试雷达高分辨距离像的目标类别标号<img file="FDA00011142167300000325.GIF" wi="67" he="54" />
地址 710071 陕西省西安市太白南路2号