CN107576949A - 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法 - Google Patents

基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法 Download PDF

Info

Publication number
CN107576949A
CN107576949A CN201710728138.7A CN201710728138A CN107576949A CN 107576949 A CN107576949 A CN 107576949A CN 201710728138 A CN201710728138 A CN 201710728138A CN 107576949 A CN107576949 A CN 107576949A
Authority
CN
China
Prior art keywords
kernel function
ker
svdd
training sample
distance
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
CN201710728138.7A
Other languages
English (en)
Other versions
CN107576949B (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 CN201710728138.7A priority Critical patent/CN107576949B/zh
Publication of CN107576949A publication Critical patent/CN107576949A/zh
Application granted granted Critical
Publication of CN107576949B publication Critical patent/CN107576949B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法,属于雷达目标识别领域。由于K型核函数具有泛化能力强,有利于提取训练数据的全局特性,避免了径向基核函数复杂的指数运算,同时具有多项式核函数计算量小和径向基核函数逼近精度高的优点。而径向基核函数具有很好的局部特性,因此本发明将K型核函数与径向基核函数进行组合替代传统的SVDD算法中的核函数;同时,本发明采用基于截断距离的局部密度算法在高维核特征空间中计算支持向量与训练样本数据之间的局部密度,根据密度的分布,调整超闭球的形状,有效的提高了雷达一维距离像单类目标的识别性能。

Description

基于密度权重和混合核函数的SVDD雷达目标一维距离像识别 方法
技术领域
本发明适用于雷达目标识别领域,具体涉及一种基于密度权重和混合核函数的SVDD雷达目标一维距离像识别。
背景技术
雷达高分辨距离像(High Resolution Range Profile,HRRP)是待识别目标散射点子回波沿雷达视线方向上投影向量和,反映了目标散射点在雷达视线上的分布情况。相对于低分辨率雷达所获得的目标雷达散射截面积(Radar cross-section,RCS)而言,能够获得更多关于目标的结构和形状等信息;相对于合成孔径雷达(Synthetic ApertureRadar,SAR)图像和逆合成孔径雷达(Inverse Synthetic Aperture Radar,ISAR)图像而言,具有易于获取和存储量小等特点,因而在雷达自动目标识别领域(Radar AutomaticTarget Recognition,RATR)受到了广泛应用。
支持向量数据描述(SVDD)是一种基于贝叶斯最优决策理论的单类分类器,其基本思想是通过非线形变换将数据点映射到高维特征空间,寻求一个包含大部分高维特征空间数据点的最小超闭球(超球体)。由于SVDD算法具有复杂度低、扩充性强,以及对训练样本数据规模上要求不高等优点,已经广泛应用于故障诊断、异常检测、目标识别等多种领域。随后提出了很多SVDD改进算法,如基于聚类的SVDD目标识别算法、基于密度权重的SVDD目标识别算法、基于多核的SVDD目标识别算法等等。
基于密度权重的SVDD目标识别类型的算法都是在数据域根据训练数据的相对分布位置,计算局部密度作为SVDD惩罚因子的权重来实现的。但是SVDD算法是将训练样本数据依据核函数映射到高维特征空间,原始数据域的密度分布并不完全等价于在高维特征空间中的密度分布,影响了SVDD算法支持向量的选择和超闭球的形状,导致识别性能的下降。同时核函数的选择是SVDD算法的核心,一旦核函数选定,映射的高维特征空间就固定不变,SVDD算法的计算,支持向量的选择都是在该空间中完成。传统SVDD算法通常是以径向基核函数为核函数,但径向基核函数存在计算量大、泛化性能较弱等缺点,从而导致整个识别算法性能的下降。
发明内容
本发明的发明目的在于:针对上述存在的问题,提出了一种基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法,以提高雷达一维距离像单类目标的识别性能。
本发明的基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法包括下列步骤:
步骤1:设置支持向量数据描述SVDD算法的核函数:
设置K型核函数Ker_k(p,q)、径向基核函数Ker_rbf(p,q),由Ker_k(p,q)、Ker_rbf(p,q)的加权求和得到组合核函数Ker(p,q)=W·Ker_rbf(p,q)+(1-W)·Ker_k(p,q),其中W∈(0,1)为组合核函数系数,p、q为对应的自变量,即雷达目标一维距离像数据;
如待处理的雷达目标一维距离像数据为g维列向量,则其中,pm、qm表示p、q的任一维,λ>0反应了K型核函数的宽度。
径向基核函数Ker_rbf(p,q)表达式通常为:σ反应了径向基型核函数的宽度。
步骤2:对训练样本(雷达目标一维距离像数据)进行高维核特征空间映射处理:
遍历所有训练样本,计算任意训练样本间的K型核函数映射结果Ker_k(xi,xj)、径向基核函数映射结果Ker_rbf(xi,xj),其中xi、xj表示训练样本集的任意训练样本,i,j=1,2,…,N,N表示训练样本数;
由Ker(xi,xj)=W·Ker_rbf(xi,xj)+(1-W)·Ker_k(xi,xj)得到训练样本的高维核特征空间映射结果Ker(xi,xj);
步骤3:计算支持向量数据描述SVDD算法的边界:
在满足约束条件A的前提下,对优化目标进行求解得到SVDD算法的边界;
所述约束条件A为:0≤βij≤C且其中β12,...,βN为每个训练样本对应的拉格朗日系数,C为预设的惩罚参数,且C>0;
步骤4:计算每个训练样本的核距离系数:
将拉格朗日系数满足0<βi<C的训练样本作为支持向量;
基于对应的组合核函数Ker(xi,xj)计算各支持向量与所有训练样本的核距离D(SVn,xi)=Ker(SVn,SVn)-2Ker(SVn,xi)+Ker(xi,xi),其中下标n为支持向量标识符;
对各支持向量的核距离进行概率分布的归一化处理,将归一化处理结果作为支持向量的训练样本的核距离系数,并将非支持向量的训练样本核距离系数置为0;
其中,对任意支持向量的核距离D(SVn,xi)进行概率分布的归一化处理具体为:
根据公式得到归一化处理结果其中 NSV表示支持向量的数目;在计算时,为了避免核距离D(SVn,xi)为0的情况(即训练样本xi就是支持向量SVp),当xi=SVn时,令核距离D(SVn,xi)为当前支持向量SVn到其它训练样本(除SVn外的所有训练样本)的核距离的最小值;
步骤5:对待识别对象z(雷达目标一维距离像数据)进行识别处理:
对待识别对象z进行高维核特征空间映射处理:
计算待识别对象z的K型核函数Ker_k(z,z)、径向基核函数Ker_rbf(z,z);
由Ker(z,z)=W·Ker_rbf(z,z)+(1-W)·Ker_k(z,z)得到待识别对象z的高维核特征空间映射结果Ker(z,z);
计算第一距离d(z)(待识别对象z与超球体球心的距离):
计算第二距离dker_density(z)(基于密度分布的待识别对象z到超球体球心的距离):
对第一、二距离进行加权求和,得到融合距离测度result(z),基于所述融合距离测度result(z)判断待识别对象z是否位于SVDD算法的边界外,若是,则待识别对象z为非真实目标数据;否则为真实目标数据。优选的方式为:将对应各支持向量的超球体半径(支持向量与超球体球心的距离得到对应超球半径)的均值作为真实目标数据与非真实目标数据的判决阈值,若result(z)大于该判决阈值,则z为非真实目标数据;否则为真实目标数据。
综上所述,由于采用了上述技术方案,本发明的有益效果是:由于K型核函数具有泛化能力强,有利于提取训练数据的全局特性,避免了径向基核函数复杂的指数运算,同时具有多项式核函数计算量小和径向基核函数逼近精度高的优点。而径向基核函数具有很好的局部特性,因此本发明将K型核函数与径向基核函数进行组合替代传统的SVDD算法中的核函数;同时,本发明采用基于截断距离的局部密度算法在高维核特征空间中计算支持向量与训练样本数据之间的局部密度,根据密度的分布,调整超闭球的形状,有效的提高了雷达一维距离像单类目标的识别性能。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式对本发明作进一步地详细描述。
实施例
为了验证本发明的有效性,进行如下仿真实验。
对5种仿真飞机目标的一维距离像数据进行了识别,包括AH64、AN26、B52、B1B、F15五种型号的飞机。雷达带宽为400MHZ,工作频率为6GHz。在目标姿态角0°~30°范围内每隔0.1°的一维距离像中,取AH64一维距离像的0°、0.2°、0.4°、...、30°作为训练样本数据;将AN26、B52、B1B、F15一维距离像的0°~15°范围的一维距离像分别与AH64中0°~30°范围内余下的样本数据组合成测试数据。
对5种目标的识别采用现有的SVDD和本发明提出的基于密度权重和混合核函数的SVDD雷达目标一维距离像识别算法进行对比,其中参数设定:径向基核函数的参数固定σ=10,惩罚因子参数固定C=0.1;K型核函数的参数λ=0.001,组合核函数系数W=0.01,融合距离测度result(z)的加权系数δ设置为0.7,即result(z)=δ·d(z)+(1-δ)·dker_density(z),分别在(0,1)、(0,1)、[0.7,1)内取λ,W,δ,步长分别为0.001、0.01、0.01,使用网格搜索法找到算法的最佳识别率,识别结果如表1所示。
表1两种方法对5种仿真飞机目标一维距离像的识别结果(%)
可以看出,AH64作为训练数据,剩余的AH64样本数据和AN26、B52、B1B、F15组成测试数据时,本发明基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法在同样σ、C的情况下,识别率都高于传统的SVDD算法。因此,可以得出结论,本发明的基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法比单一的SVDD算法具有更好的识别效果。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。

Claims (3)

1.基于密度权重和混合核函数的SVDD雷达目标一维距离像识别方法,其特征在于,包括下列步骤:
步骤1:设置支持向量数据描述SVDD算法的核函数:
设置K型核函数Ker_k(p,q)、径向基核函数Ker_rbf(p,q),由Ker_k(p,q)、Ker_rbf(p,q)的加权求和得到组合核函数Ker(p,q)=W·Ker_rbf(p,q)+(1-W)·Ker_k(p,q),其中W∈(0,1)为组合核函数系数,p、q为对应的自变量;
步骤2:对训练样本进行高维核特征空间映射处理:
遍历所有训练样本,计算任意训练样本间的K型核函数映射结果Ker_k(xi,xj)、径向基核函数映射结果Ker_rbf(xi,xj),其中xi、xj表示训练样本集的任意训练样本,i,j=1,2,…,N,N表示训练样本数;
由Ker(xi,xj)=W·Ker_rbf(xi,xj)+(1-W)·Ker_k(xi,xj)得到训练样本的高维核特征空间映射结果Ker(xi,xj);
步骤3:计算支持向量数据描述SVDD算法的边界:
在满足约束条件A的前提下,对优化目标进行求解得到SVDD算法的边界;
所述约束条件A为:0≤βij≤C且其中β12,...,βN为每个训练样本对应的拉格朗日系数,C为预设的惩罚参数,且C>0;
步骤4:计算每个训练样本的核距离系数:
将拉格朗日系数满足0<βi<C的训练样本作为支持向量;
基于对应的Ker(xi,xj)计算各支持向量与所有训练样本的核距离D(SVn,xi)=Ker(SVn,SVn)-2Ker(SVn,xi)+Ker(xi,xi),其中下标n为支持向量标识符;
对各支持向量的核距离进行概率分布的归一化处理,将归一化处理结果作为支持向量的训练样本的核距离系数,并将非支持向量的训练样本核距离系数置为0;
步骤5:对待识别对象z进行识别处理:
对待识别对象z进行高维核特征空间映射处理:
计算待识别对象z的K型核函数Ker_k(z,z)、径向基核函数Ker_rbf(z,z);
由Ker(z,z)=W·Ker_rbf(z,z)+(1-W)·Ker_k(z,z)得到待识别对象z的高维核特征空间映射结果Ker(z,z);
计算第一距离
计算第二距离
对第一、二距离进行加权求和,得到融合距离测度result(z),基于所述融合距离测度result(z)判断待识别对象z是否位于SVDD算法的边界外,若是,则待识别对象z为非真实目标数据;否则为真实目标数据。
2.如权利要求1所述的方法,其特征在于,对任意支持向量的核距离D(SVn,xi)进行概率分布的归一化处理具体为:
根据公式得到归一化处理结果
其中NSV表示支持向量的数目;在计算step1(n)时,若xi=SVn,则令核距离D(SVn,xi)为支持向量SVn到除SVn外的所有训练样本的核距离的最小值。
3.如权利要求1或2所述的方法,其特征在于,步骤5中,将对应各支持向量的超球体半径的均值作为真实目标数据与非真实目标数据的判决阈值,若result(z)大于所述判决阈值,则z为非真实目标数据;否则为真实目标数据。
CN201710728138.7A 2017-08-23 2017-08-23 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法 Active CN107576949B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710728138.7A CN107576949B (zh) 2017-08-23 2017-08-23 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710728138.7A CN107576949B (zh) 2017-08-23 2017-08-23 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法

Publications (2)

Publication Number Publication Date
CN107576949A true CN107576949A (zh) 2018-01-12
CN107576949B CN107576949B (zh) 2020-03-27

Family

ID=61034451

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710728138.7A Active CN107576949B (zh) 2017-08-23 2017-08-23 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法

Country Status (1)

Country Link
CN (1) CN107576949B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416393A (zh) * 2018-03-16 2018-08-17 西安电子科技大学 一种基于高分辨一维距离像的目标识别优化方法
CN108985268A (zh) * 2018-08-16 2018-12-11 厦门大学 基于深度迁移学习的归纳式雷达高分辨距离像识别方法
CN109117793A (zh) * 2018-08-16 2019-01-01 厦门大学 基于深度迁移学习的直推式雷达高分辨距离像识别方法
CN109901557A (zh) * 2019-03-28 2019-06-18 山东科技大学 一种基于数据域描述的可变控制限故障检测方法
CN109993224A (zh) * 2019-03-27 2019-07-09 中国人民解放军战略支援部队航天工程大学 基于深度学习与多核学习的geo卫星形状与姿态识别方法
CN110135444A (zh) * 2019-06-17 2019-08-16 电子科技大学 一种目标库属性一维像判别最优本征子空间特征提取方法
CN110320894A (zh) * 2019-08-01 2019-10-11 陕西工业职业技术学院 一种准确划分混叠区域数据类别的火电厂制粉系统故障检测方法
CN113156416A (zh) * 2021-05-17 2021-07-23 电子科技大学 一种基于多重核字典学习的未知目标判别方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810704A (zh) * 2014-01-23 2014-05-21 西安电子科技大学 基于支持向量机和判别随机场的sar图像变化检测方法
CN104008403A (zh) * 2014-05-16 2014-08-27 中国人民解放军空军装备研究院雷达与电子对抗研究所 一种svm(矢量机)模式的多目标识别判定方法
CN104112143A (zh) * 2014-07-23 2014-10-22 大连民族学院 基于加权超球支持向量机算法的图像分类方法
CN106203481A (zh) * 2016-06-30 2016-12-07 中国人民解放军军械工程学院 基于混合核rvm的电子装备状态预测方法
CN106886213A (zh) * 2017-03-13 2017-06-23 北京化工大学 一种基于核相似度支持向量数据描述的间歇过程故障检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103810704A (zh) * 2014-01-23 2014-05-21 西安电子科技大学 基于支持向量机和判别随机场的sar图像变化检测方法
CN104008403A (zh) * 2014-05-16 2014-08-27 中国人民解放军空军装备研究院雷达与电子对抗研究所 一种svm(矢量机)模式的多目标识别判定方法
CN104112143A (zh) * 2014-07-23 2014-10-22 大连民族学院 基于加权超球支持向量机算法的图像分类方法
CN106203481A (zh) * 2016-06-30 2016-12-07 中国人民解放军军械工程学院 基于混合核rvm的电子装备状态预测方法
CN106886213A (zh) * 2017-03-13 2017-06-23 北京化工大学 一种基于核相似度支持向量数据描述的间歇过程故障检测方法

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416393A (zh) * 2018-03-16 2018-08-17 西安电子科技大学 一种基于高分辨一维距离像的目标识别优化方法
CN108985268B (zh) * 2018-08-16 2021-10-29 厦门大学 基于深度迁移学习的归纳式雷达高分辨距离像识别方法
CN108985268A (zh) * 2018-08-16 2018-12-11 厦门大学 基于深度迁移学习的归纳式雷达高分辨距离像识别方法
CN109117793A (zh) * 2018-08-16 2019-01-01 厦门大学 基于深度迁移学习的直推式雷达高分辨距离像识别方法
CN109117793B (zh) * 2018-08-16 2021-10-29 厦门大学 基于深度迁移学习的直推式雷达高分辨距离像识别方法
CN109993224A (zh) * 2019-03-27 2019-07-09 中国人民解放军战略支援部队航天工程大学 基于深度学习与多核学习的geo卫星形状与姿态识别方法
CN109993224B (zh) * 2019-03-27 2021-02-02 中国人民解放军战略支援部队航天工程大学 基于深度学习与多核学习的geo卫星形状与姿态识别方法
CN109901557A (zh) * 2019-03-28 2019-06-18 山东科技大学 一种基于数据域描述的可变控制限故障检测方法
CN109901557B (zh) * 2019-03-28 2021-06-22 山东科技大学 一种基于数据域描述的可变控制限故障检测方法
CN110135444A (zh) * 2019-06-17 2019-08-16 电子科技大学 一种目标库属性一维像判别最优本征子空间特征提取方法
CN110135444B (zh) * 2019-06-17 2022-02-08 电子科技大学 一种目标库属性一维像判别最优本征子空间特征提取方法
CN110320894A (zh) * 2019-08-01 2019-10-11 陕西工业职业技术学院 一种准确划分混叠区域数据类别的火电厂制粉系统故障检测方法
CN113156416A (zh) * 2021-05-17 2021-07-23 电子科技大学 一种基于多重核字典学习的未知目标判别方法
CN113156416B (zh) * 2021-05-17 2022-05-17 电子科技大学 一种基于多重核字典学习的未知目标判别方法

Also Published As

Publication number Publication date
CN107576949B (zh) 2020-03-27

Similar Documents

Publication Publication Date Title
CN107576949A (zh) 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法
FARSHCHIAn et al. The Pareto distribution for low grazing angle and high resolution X-band sea clutter
Du et al. A two-distribution compounded statistical model for radar HRRP target recognition
CN105607045B (zh) 一种欺骗式干扰下组网雷达的优化布站方法
CN105137418B (zh) 基于全邻模糊聚类的多目标跟踪与数据互联方法
CN109752710B (zh) 一种基于稀疏贝叶斯学习的快速目标角度估计方法
CN108256436A (zh) 一种基于联合分类的雷达hrrp目标识别方法
CN109840542B (zh) 基于极化特征的自适应维度决策树分类方法
CN107843875A (zh) 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法
CN107229032A (zh) 一种构建四阵元立体阵列的方法和装置
CN106021697A (zh) 一种快速相控阵雷达时间-能量资源联合管理方法
Tivive et al. An improved SVD-based wall clutter mitigation method for through-the-wall radar imaging
CN111711432B (zh) 一种基于ukf和pf混合滤波的目标跟踪算法
CN110083915A (zh) 一种高频区雷达目标散射中心正向自动化建模方法
CN113900103A (zh) 突发性强对流灾害天气自适应快速识别预警改进算法
CN107133648A (zh) 基于自适应多尺度融合稀疏保持投影的一维距离像识别方法
CN104268574A (zh) 一种基于遗传核模糊聚类的sar图像变化检测方法
CN109597048A (zh) 基于二维卷积神经网络的米波雷达doa估计方法
CN106408571B (zh) 一种基于最优模糊因子选择的可变类遥感图像分割方法
CN110852019A (zh) 一种雷达目标的散射中心建模方法
CN105184829A (zh) 一种紧密空间目标检测和高精度质心定位方法
CN113255603A (zh) 基于黎曼流形监督降维的增强矩阵恒虚警率检测方法
CN104268565A (zh) 基于回归学习的景象匹配区选取方法
CN116030252A (zh) 融合边缘信息的模糊迭代sar图像超像素分割方法
CN103699811A (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