CN109444840B - 一种基于机器学习的雷达杂波抑制方法 - Google Patents
一种基于机器学习的雷达杂波抑制方法 Download PDFInfo
- Publication number
- CN109444840B CN109444840B CN201811471743.1A CN201811471743A CN109444840B CN 109444840 B CN109444840 B CN 109444840B CN 201811471743 A CN201811471743 A CN 201811471743A CN 109444840 B CN109444840 B CN 109444840B
- Authority
- CN
- China
- Prior art keywords
- cluster
- data
- clutter
- distance
- clustering
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000010801 machine learning Methods 0.000 title claims abstract description 12
- 230000001629 suppression Effects 0.000 title claims abstract description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 125000001475 halogen functional group Chemical group 0.000 claims description 5
- 238000011478 gradient descent method Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 10
- 238000007781 pre-processing Methods 0.000 abstract description 3
- 238000007796 conventional method Methods 0.000 abstract 1
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
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
- 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
- G01S7/414—Discriminating targets with respect to background clutter
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)
Abstract
本发明公开了一种基于机器学习的雷达杂波抑制方法,实测数据产生包含杂波和目标的数据集作为输入数据,对输入数据进行预处理,初步处理后选取合适特征和初始权重进行加权,然后对其聚类,将聚为不同类别结果进行评价,再反馈到特征加权模块调整权值,直到找到评价指标最优的一组权重和聚类数。最佳聚类数中,一类为目标,其余类为不同类型杂波,可将实时雷达回波数据与之进行比较,实现杂波抑制,从而得到目标信号。相比于以往利用杂波图抑制杂波的方法,实现结构简单,实现成本低,效果更好,且具有很高的简捷性、实时性和通用性。
Description
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种基于机器学习的雷达杂波抑制方法。
背景技术
随着低空空域的逐步开放以及人们低空活动的日益增加,低空空域的交通状况越来越复杂,为了维持低空飞行秩序,保证低空飞行器安全,实施低空空域飞行器的监视,已经成为空管发展的必然趋势。而低慢小目标RCS小,速度慢,且信杂比低,回波易淹没在杂波中,具有“管控难、侦测难、处置难”等特点。另外,低慢小探测雷达不同于一般雷达,需要对超低空杂波进行处理,要在成片的杂波中很快的分辨出目标回波也是不容易的。由于出现大量的杂波,给雷达终端处理***带来很大的难度。故将杂波抑制是一个亟待解决的问题,对低慢小目标探测雷达来说至关重要。
目前,基于机器学***均值),然后迭代的进行分配点和更新类簇中心点的步骤,直至类簇中心点的变化很小,或者达到指定的迭代次数。这种算法思想比较简单,但是初始的K值和初始类簇中心点的选组对于聚类效果影响较大,且该方法仅适用于球形分布数据。
发明内容
发明目的:本发明提供一种基于机器学习的雷达杂波抑制方法,实现结构简单,聚类效果好,实现成本低,具有很高的通用性和简捷性。
技术方案:本发明所述的一种基于机器学习的雷达杂波抑制方法,包括以下步骤:
(1)将通过生成对抗网络产生的包含杂波和目标的数据集作为输入数据,并对输入数据作归一化处理;
(2)选取合适特征,构建特征向量,并对其进行特征加权;
(3)对步骤(2)的特征向量进行快速聚类,将聚为不同类别结果进行评价;
(4)将步骤(3)的评价结果,反馈到步骤(2),利用梯度下降法直到找到评价指标戴维森保丁指数最优的一组权重和聚类数。
步骤(2)所述特征主要包括距离、方位、俯仰、RCS、速度、差通道能量、保护通道。
步骤(2)所述的对初始权重进行加权通过以下公式实现:
其中,xi是第i个数据的m维特征向量,xi=(xi1,xi2,···xit,···xim),i=1,2,···,n。
所述步骤(3)包括以下步骤:
(31)给定用于确定截断距离dc的参数t∈(0,1),算距离dij,并令
dji=dij,i<j,i,j∈IS
(33)对非聚类中心数据点进行归类,按照ρ的值从大到小遍历;
(34)若nc>1,则将每个cluster中的数据点进一步分为cluster core和clusterhalo。
所述步骤(4)所述的戴维森保丁指数计算如下:
其中,||xiWi-xjWj||2对应于簇Ci与簇Cj中心点的距离,avg(C)对应于簇C内样本间的平均距离,k为类别数。
有益效果:与现有技术相比,本发明的有益效果:1、本发明采用的快速聚类算法,相比于常用的K-means算法,只需要考虑点与点之间的距离,基于密度计算,且可以得到很好的非球形聚类效果,通过降低分辨率后的特征加权,也降低了计算复杂度;2、本发明实现结构简单,聚类效果更好,实现成本低,具有更好的通用性和简捷性。
附图说明
图1为本发明的流程图
图2基于密度的新型聚类算法流程图;
图3数据点的γ值分布情况;
图4DBI随聚类数的变化情况。
具体实施方式
下面结合附图对本发明作进一步详细描述。
本发明提供一种基于机器学习的雷达杂波抑制方法,如图1所示。
1、将通过生成对抗网络产生的包含杂波和目标的数据集作为输入数据,并对输入数据作归一化处理。
实测数据既包含杂波也包含目标,将其作为真实的数据集,利用生成网络造出“假”数据集,并通过判别网络来判别是真实数据集还是假样本集,最终的输出效果是判别网络给出结果是一个接近0.5的值,说明“假”数据集已经十分贴合真实数据集了。故可以用其产生大量“假”数据集,作为本***输入数据。对输入数据作归一化处理。
2、选取合适特征,构建特征向量,并对其进行特征加权
特征在于选取距离、方位、俯仰、RCS、速度、差通道能量、保护通道作为特征向量。在特征加权模块中,对于给定的数据集S={x1,x2,···,xi,···,xn},其中,xi是描述第i个数据的m维特征向量,xi=(xi1,xi2,···xit,···xim),i=1,2,···,n。数据集中数据点本身等于xi的个数与数据点总数的比值反映了该数据点的统计分布情况,所以基于特征加权的计算方法为:
3、对其进行快速聚类,将聚为不同类别结果进行评价
聚类的核心在于该算法是基于密度的,对于每个数据点xi,提出两个变量—数据点i所在位置的局部密度ρi和到更高密度点的最短距离δi。这两个量都只取决于两个数据点xi和xj之间的距离dij,在本方法中,需要计算数据点的两两距离,对于大样本的数据集来说,计算成本较高,但特征加权使得数据量减少,降低了计算复杂度,基于特征加权的数据点距离计算表示如下:
对于任一数据点xi,选取高斯核的局部密度ρi计算方式如下:
dc为截断距离,需事先指定,局部密度表示所有样本点中(不含i点)和i点距离在dc之内的样本点的数量,与xi距离小于dc的数据点越多,ρi的值越大。最短距离δi表示比i点密度高的所有点中最近的距离,对应的点也被称为最近邻点,定义如下:
这样可以保证在ρ值相同的情况下,经降序排列后,只将排在最前面的数据点作为聚类中心。聚类中心同时具有较大的ρ值和δ值,原聚类算法中画出ρ和δ的二维决策图。在本方法中,考量这两个值的乘积即γi=ρiδi,i∈IS,取较大的数据点作为聚类中心。故先对γ降序排序并以下标为横轴,值为纵轴画在二维坐标平面内。显然,值越大越有可能是聚类中心,故从非聚类中心到聚类中心有一个明显的跳跃,用数值检测可以判断出来,避免人为的选取中心。具体步骤,如图2所示,主要分为聚类预处理即选取聚类中心和聚类算法实施。为描述方便,先引入一些记号作为算法实施步骤的基础,待聚类的数据集包含nc个cluster。
(1)初始化及预处理:给定用于确定截断距离dc的参数t∈(0,1);计算距离dij,并令dji=dij,i<j,i,j∈IS;确定截断距离dc;计算并生成其降序排列计算和并对γ值降序排列,以序号-值为横纵坐标画在二维平面内,如图3所示。
(3)对非聚类中心数据点进行归类,按照ρ的值从大到小遍历。
(4)若nc>1,则将每个cluster中的数据点进一步分为cluster core和clusterhalo。初始化标记hi=0,i∈IS;为每一个cluster生成一个平均局部密度上限标识cluster halo。
4、将步骤3的评价结果,反馈到步骤2,利用梯度下降法直到找到评价指标DBI最优的一组权重和聚类数。
实施聚类算法后用戴维森保丁指数(Davies-Bouldin Index,DBI)指标需对其结果进行评价,将其作为目标优化函数,并反馈到特征加权模块,利用梯度下降法迭代的更新权重值,直至找出最优的权重和类别数。DBI计算如下:
其中,||xiWi-xjWj||2对应于簇Ci与簇Cj中心点的距离,avg(C)对应于簇C内样本间的平均距离,k为类别数。最佳聚类数中,一类为目标,其余类为不同类型杂波,实时雷达回波数据与之进行比较,实现杂波抑制,从而得到目标信号。DBI的值越小,聚类效果越优。如图4所示,DBI随着聚类数的变化情况,可以看出,聚类数为3时,DBI最小,效果最佳。
因此本发明提出的杂波抑制方法无需采集杂波图实现对消,具有实现结构简单,聚类效果更好,实现成本低的优势,具有很高的通用性和简捷性。
Claims (4)
1.一种基于机器学习的雷达杂波抑制方法,其特征在于,包括以下步骤:
(1)将通过生成对抗网络产生的包含杂波和目标的数据集作为输入数据,并对输入数据作归一化处理;
(2)选取合适特征,构建特征向量,并对其进行特征加权;
(3)对步骤(2)的特征向量进行快速聚类,将聚为不同类别结果进行评价;
(4)将步骤(3)的评价结果,反馈到步骤(2),利用梯度下降法直到找到评价指标戴维森保丁指数最优的一组权重和聚类数;
所述步骤(3)包括以下步骤:
(31)基于特征加权的数据点距离计算两个数据点xi和xj之间的距离dij,并令dji=dij,i<j,i,j∈IS:
(33)对于任一数据点xi,选取高斯核的局部密度ρi,对非聚类中心数据点进行归类,按照局部密度ρ的值从大到小遍历,ρi计算方式如下:
其中,dc为截断距离;
(34)nc为cluster的个数,若nc>1,则将每个cluster中的数据点进一步分为clustercore和cluster halo。
2.根据权利要求1所述的一种基于机器学习的雷达杂波抑制方法,其特征在于,步骤(2)所述特征主要包括距离、方位、俯仰、RCS、速度、差通道能量、保护通道。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811471743.1A CN109444840B (zh) | 2018-12-04 | 2018-12-04 | 一种基于机器学习的雷达杂波抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811471743.1A CN109444840B (zh) | 2018-12-04 | 2018-12-04 | 一种基于机器学习的雷达杂波抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444840A CN109444840A (zh) | 2019-03-08 |
CN109444840B true CN109444840B (zh) | 2020-12-25 |
Family
ID=65555248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811471743.1A Active CN109444840B (zh) | 2018-12-04 | 2018-12-04 | 一种基于机器学习的雷达杂波抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444840B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2734071C1 (ru) * | 2020-04-16 | 2020-10-12 | Акционерное общество «Всероссийский научно-исследовательский институт радиотехники» (АО «ВНИИРТ») | Межобзорное устройство картографирования пассивных помех для импульсно-доплеровских радиолокационных станций и способ его осуществления |
CN112698281A (zh) * | 2020-12-09 | 2021-04-23 | 北京百度网讯科技有限公司 | 雷达数据处理方法、***、装置、电子设备和存储介质 |
CN112597820A (zh) * | 2020-12-10 | 2021-04-02 | 南京长峰航天电子科技有限公司 | 一种基于雷达信号分选的目标聚类方法 |
CN113702923B (zh) * | 2021-08-27 | 2023-07-18 | 电子科技大学 | 一种雷达杂波模拟方法及电子设备 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4001820A (en) * | 1971-08-05 | 1977-01-04 | The United States Of America As Represented By The Secretary Of The Navy | Multiple hyperplane recognizer |
CN102254020A (zh) * | 2011-07-22 | 2011-11-23 | 西安电子科技大学 | 基于特征权重的全局k-均值聚类方法 |
CN104766100A (zh) * | 2014-10-22 | 2015-07-08 | 中国人民解放军电子工程学院 | 基于机器学习的红外小目标图像背景预测方法与装置 |
CN106093946A (zh) * | 2016-05-27 | 2016-11-09 | 四川九洲空管科技有限责任公司 | 一种适用于场面监视雷达的目标凝聚方法及装置 |
CN107064894A (zh) * | 2017-01-11 | 2017-08-18 | 南京御达电信息技术有限公司 | 一种基于深度学习的杂波抑制方法 |
EP3239656A1 (fr) * | 2016-04-29 | 2017-11-01 | Thales | Procédé d'optimisation de la détection de cibles marines et radar aéroporté mettant en oeuvre un tel procédé |
CN107390187A (zh) * | 2017-07-28 | 2017-11-24 | 西安电子科技大学 | 基于混合高斯分布的杂波抑制方法 |
CN108805028A (zh) * | 2018-05-05 | 2018-11-13 | 南京理工大学 | 基于电磁强散射点的sar图像地面目标检测与定位方法 |
CN108872948A (zh) * | 2018-08-01 | 2018-11-23 | 哈尔滨工业大学 | 一种高频地波雷达电离层杂波抑制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5933245B2 (ja) * | 2011-12-15 | 2016-06-08 | 三菱電機株式会社 | クラッタ抑圧装置 |
CN104091038A (zh) * | 2013-04-01 | 2014-10-08 | 太原理工大学 | 基于大间隔分类准则的多示例学习特征加权方法 |
CN104239895B (zh) * | 2014-09-03 | 2017-05-24 | 西安电子科技大学 | 基于特征降维的sar目标鉴别方法 |
-
2018
- 2018-12-04 CN CN201811471743.1A patent/CN109444840B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4001820A (en) * | 1971-08-05 | 1977-01-04 | The United States Of America As Represented By The Secretary Of The Navy | Multiple hyperplane recognizer |
CN102254020A (zh) * | 2011-07-22 | 2011-11-23 | 西安电子科技大学 | 基于特征权重的全局k-均值聚类方法 |
CN104766100A (zh) * | 2014-10-22 | 2015-07-08 | 中国人民解放军电子工程学院 | 基于机器学习的红外小目标图像背景预测方法与装置 |
EP3239656A1 (fr) * | 2016-04-29 | 2017-11-01 | Thales | Procédé d'optimisation de la détection de cibles marines et radar aéroporté mettant en oeuvre un tel procédé |
CN106093946A (zh) * | 2016-05-27 | 2016-11-09 | 四川九洲空管科技有限责任公司 | 一种适用于场面监视雷达的目标凝聚方法及装置 |
CN107064894A (zh) * | 2017-01-11 | 2017-08-18 | 南京御达电信息技术有限公司 | 一种基于深度学习的杂波抑制方法 |
CN107390187A (zh) * | 2017-07-28 | 2017-11-24 | 西安电子科技大学 | 基于混合高斯分布的杂波抑制方法 |
CN108805028A (zh) * | 2018-05-05 | 2018-11-13 | 南京理工大学 | 基于电磁强散射点的sar图像地面目标检测与定位方法 |
CN108872948A (zh) * | 2018-08-01 | 2018-11-23 | 哈尔滨工业大学 | 一种高频地波雷达电离层杂波抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109444840A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444840B (zh) | 一种基于机器学习的雷达杂波抑制方法 | |
CN108564192B (zh) | 一种基于气象因子权重相似日的短期光伏功率预测方法 | |
CN110443281B (zh) | 基于hdbscan聚类的文本分类自适应过采样方法 | |
CN105005789B (zh) | 一种基于视觉词汇的遥感图像地物分类方法 | |
CN107576949B (zh) | 基于密度权重和混合核函数的svdd雷达目标一维距离像识别方法 | |
CN108363810A (zh) | 一种文本分类方法及装置 | |
CN110751121B (zh) | 基于聚类与sofm的无监督雷达信号分选方法 | |
CN112149715B (zh) | 一种基于深度学习的无人机空战威胁评估方法 | |
CN113095442A (zh) | 基于半监督学习在多维度雷达数据下的冰雹识别方法 | |
CN109150830B (zh) | 一种基于支持向量机和概率神经网络的层次入侵检测方法 | |
CN111860359B (zh) | 一种基于改进随机森林算法的点云分类方法 | |
CN103235954A (zh) | 一种基于改进AdaBoost算法的地基云图识别方法 | |
CN112613450A (zh) | 一种增强在困难样本上表现的3d目标检测方法 | |
CN115220007A (zh) | 一种针对姿态识别的雷达点云数据增强方法 | |
CN115908930A (zh) | 一种基于改进cfwpso-svm的前视声呐图像识别分类方法 | |
CN117725448A (zh) | 气象导航信号特征的聚类分析方法 | |
CN114219228A (zh) | 一种基于em聚类算法的体育场疏散评价方法 | |
Jiang et al. | A method of forecasting thunderstorms and gale weather based on multisource convolution neural network | |
CN115797814B (zh) | 基于无人机的输电线路树障隐患实时评估计算方法 | |
Su et al. | Unsupervised k-means combined with SOFM structure adaptive radar signal sorting algorithm | |
CN109934344B (zh) | 一种改进的基于规则模型的多目标分布估计方法 | |
CN115034690B (zh) | 一种基于改进模糊c-均值聚类的战场态势分析方法 | |
CN117079120A (zh) | 一种基于改进ga算法的目标识别模型优化方法 | |
CN110163437A (zh) | 基于DPK-means的日前光伏发电功率预测方法 | |
CN115146742A (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210906 Address after: 210000 108-25, science and technology entrepreneurship R & D incubation complex building (Building 5), Nanjing Baixia high tech Industrial Park, No. 5, Yongzhi Road, Qinhuai District, Nanjing, Jiangsu Province Patentee after: Nanjing Liuji Photoelectric Technology Research Institute Co.,Ltd. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee before: Nanjing University of Aeronautics and Astronautics |
|
TR01 | Transfer of patent right |