CN107238817B - 一种参数自适应设置及自动调整的雷达辐射源信号分选方法 - Google Patents
一种参数自适应设置及自动调整的雷达辐射源信号分选方法 Download PDFInfo
- Publication number
- CN107238817B CN107238817B CN201710539573.5A CN201710539573A CN107238817B CN 107238817 B CN107238817 B CN 107238817B CN 201710539573 A CN201710539573 A CN 201710539573A CN 107238817 B CN107238817 B CN 107238817B
- Authority
- CN
- China
- Prior art keywords
- algorithm
- class
- classes
- sample data
- signal
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23211—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with adaptive number of clusters
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种参数自适应设置及自动调整的雷达辐射源信号分选方法。该方法首先在特征向量的选取中结合脉冲到达角的同时引入脉内特征——对称Holder系数进行分选,然后方法基于迭代自组织数据分析算法,算法简称为ISODATA,根据特征向量自身的特征,利用搜索直方图峰值统计的方法设置预期类数,同时对ISODATA算法类内各数据分布的标准差的上限和各类中心间最小距离的下限进行初步设定,最后在ISODATA算法运行时对类内各数据分布的标准差的上限和各类中心间最小距离的下限控制设定自适应准则达到准确最终分选的目的,通过这种方式解决信号分选问题。
Description
技术领域
本发明涉及一种雷达辐射源信号分选方法。
背景技术
随着现代雷达技术的快速发展,尤其是各式新体制雷达的广泛应用,战场电磁环境日益复杂多变,其主要表现为同时同地的辐射源信号严重交错,信号调制形式复杂多样,这使得由基本参数所描述的特征空间边界严重交叠,基于传统五大基本参数的辐射源信号分选基本不能满足实际的需要。
脉内特征分析方法为减小特征空间的交叠开辟了新的研究思路,它使得提取新的信号参数以进一步充分表征雷达辐射源信号特征的想法成为现实。
ISODATA算法的优点在于通过逐步分化的过程可以认识到大部分待分类对象在分类伊始中不太可能显示的本质属性,它是根据样本的差异来自行决定最终的聚类数目,这使得分类更加科学。但是ISODATA算法也有其明显的缺点:传统ISODATA聚类算法参数的预设往往是根据经验设置,而现代战场电磁环境下,雷达辐射源信号种类数量繁多,变幻莫测,根本无法获得先验信息对算法的参数进行准确预设,从而导致聚类的结果差强人意。
因此,本发明把脉内特征参数——对称Holder系数补充到分选参数当中,通过自适应调整ISODATA算法的方式,合理地实现雷达辐射源信号的分选。
发明内容
为了解决雷达辐射源信号的准确分选问题,本发明公开了一种参数自适应设置及自动调整的雷达辐射源信号分选方法。该方法采用的特征向量是脉冲到达角和对称Holder系数这两个参数,基于迭代自组织数据分析算法,算法简称为ISODATA,首先根据特征向量自身的特征,利用搜索直方图峰值统计的方法设置预期类数,同时对类内各数据分布的标准差的上限和各类中心间最小距离的下限进行初步设定,在算法运行时对类内各数据分布的标准差的上限和各类中心间最小距离的下限控制设定自适应准则,通过这种方式解决信号分选问题。
本发明的详细技术措施步骤如下:
第一步:特征向量的选取
选用雷达辐射源信号脉冲到达角和对称Holder系数构成特征向量。
(一)雷达辐射源脉冲到达角的提取
提取雷达信号脉冲流的脉冲,获取雷达信号脉冲流脉冲描述字中的脉冲到达角。
(二)对称Holder系数的提取
提取雷达信号脉冲流的脉冲,并对提取的信号进行预处理,主要包括对信号频域的转换,带宽和能量的归一化。
根据下式,计算脉冲的频谱和矩形信号以及三角形信号的对称Holder系数HU和HT。
(1)与矩形信号U(k)的对称Holder系数
(2)与三角形信号T(k)的对称Holder系数
式中,F(i)为信号的频谱,同时选定p=5,q=5。
构造特征向量B=(DOA,HT)和D=(DOA,HT,HU),其中B用于搜索直方图峰值统计,D用于最终的辐射源信号分选。
第二步:搜索直方图峰值统计的方法来预估聚类的大致数目
得到特征向量B后,通过搜索直方图峰值统计的方法设置ISODATA算法的预期类数。
记特征向量B={Bi}={B1,B2,...,BN}(i=1,2,...,N),特征向量的维数为2维。
划分区间由样本数量决定。设每一维划分区间数为q,
记DOAmax,DOAmin,HTmax,HTmin分别为特征向量中DOA和HT的最大值和最小值。
则向量DOA维区间宽度DOAΔ和向量HT维区间宽度HTΔ分别为
至此将特征向量B所在的平面区域划分为Q=q2个小区间,每个小区间可以表示为
{bij=(DOAi,HTj)|DOAi∈[DOAm,DOAm+1],HTj∈[HTn,HTn+1],m,n∈[1,q]} (8)
式中,DOAm=DOAmin+(m-1)DOAΔ,HTn=HTmin+(m-1)HTΔ。
然后分别统计落于不同小区间中向量元素的个数,到区间统计值的直条高度。将每一个划分好区间中统计值与相邻区间进行比较,若该区间的统计值大于其所有相邻区间的统计值,则记该区间为一个峰值区间。最后统计峰值区间的数目即为峰值数目S。
预设类数设置为直方图峰值的数目即
第三步:类内各数据分布的标准差的上限和各类中心间最小距离的下限的预先设置
得到特征向量D,预先设置类内各数据分布的标准差的上限θS和各类中心间最小距离的下限θD。
首先对D中每一维进行属性归一化:
式中:uj为属性归一化样本数据xi'j第j维的中心。
继续求取样本数据第j维的分布稀疏性:
λj是样本数据中第j维规范后的标准化偏差。
由此设置初始值:
式中,S为预期类数;
然后直接求样本数据自身第j维的方差
由此设置初始值
第四步:根据自适应准则运行ISODATA算法进行分选
对于整个算法的运行,采用自适应原则,即通过设置类内各数据分布的标准差的上限θS和各类中心间最小距离的下限θD的变换准则,达到自动逼近最佳聚类效果的目的。
设置参数后运行ISODATA算法。
每次聚类的算法结束后,分别计算提取结果的类内距离ΘS'和类间距离ΘD'。
式中:c'为算法运行产生的分类数,Cp为归类数据的集合,vi为第i类的中心,xn为属于Cp类的数据,vj为与vi不同的类的中心。
第一次聚类后仅改变θS,设定自适应原则如下
再次运行ISODATA算法
并设定自适应的准则如下:
而退出自适应迭代的条件为:
Δs=|θs'-θs|=0 (21)
ΔD=|θD'-θD|=0 (22)
式中:ΔS、ΔD为算法前后连续两次运行结果的类内距离以及类间距离的差值。
退出自适应迭代后即完成分选,进而统计分选结果。
至此,一种参数自适应设置及自动调整的雷达辐射源信号分选方法的整个过程结束。
本发明具有如下有益效果:
①分选准确率高。本发明中,预期类数设定的相对准确以及自适应准则的设定能够提高准确率,另外,引入对称Holder系数这一脉内特征,大大增加了特征向量的可分性,所以分选效果相对理想。
附图说明
图1是基于自适应调整的雷达辐射源信号分选方法的流程图;
图2是雷达样本数据的分布图;
图3是搜索直方图峰值统计效果图;
图4为理想分类图;
图5为算法分选效果图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细描述。
采用图1所示的方法,选取5种不同调制类型的且参数较为接近的雷达辐射源信号进行仿真实验,工作参数变化类型及取值范围如表1所示。
表1雷达辐射源信号特征参数表
(1)特征向量的选取
对雷达辐射源信号进行特征向量提取,每种类型的雷达产生500个样本数据。
(2)搜索直方图峰值统计的方法来预估聚类的大致数目
样本数据分布于平面如图2所示,然后经由搜索直方图峰值统计,得到直方图峰值图如图3所示,分析后峰值数目为
S=5
则预设分类数目
(3)类内各数据分布的标准差的上限和各类中心间最小距离的下限的预先设置
继续求取样本数据第j维的分布稀疏性:
得到
λ1=0.2330,λ2=0.2248,λ3=0.2596
由此设置初始值:
然后直接求样本数据自身第j维的方差
得到
σ1=1.7908,σ2=1.9398,σ3=2.6530
(4)根据自适应原则运行ISODATA算法进行分类
2)设置基本参数
每个分类中允许的最少数据个数:θn=4;
迭代时最多可合并的类的对数:L=2
最多允许的迭代运算的次数:I=70;
预期类数:C=8;
类内各数据分布的标准差的上限:θs=0.35816;
不同类中心间最小距离的下限:θD=1.124。
2)运行ISODATA算法
得到第一次结果后计算类内距离和类间距离
Θs=1067.1
ΘD=2.4872
改变类内各数据分布的标准差的上限后再次运行算法
Θs=1067.1
ΘD=2.4872
改变参数后再次运行算法
Θs=1067.1
ΘD=2.4872
此时
Δs=|θs'-θs|=0
ΔD=|θD'-θD|=0
退出算法,得到分选结果如表2所示。样本数据实际分布如图4所示,方法分类如图5所示,图中不同的标签代表不同的类。
表2仿真结果
Claims (4)
1.一种参数自适应设置及自动调整的雷达辐射源信号分选方法,其特征在于该方法的步骤如下:
①选用雷达辐射源信号脉冲到达角和对称Holder系数构成特征向量,其中脉冲到达角以及信号频谱和三角形信号的对称Holder系数构成搜索直方图峰值统计的特征向量;脉冲到达角以及信号频谱和矩形信号、三角形信号的对称Holder系数构成分选特征向量;
②利用搜索直方图峰值统计的方法确定的预期类数,预期类数为峰值数目的1.5倍;
③预先设置ISODATA算法中的类内各数据分布的标准差的上限θS和各类中心间最小距离的下限θD,具体为
首先对样本数据中每一维进行属性归一化:
式中:uj为归一化样本数据xij'第j维的中心;
继续求取样本数据第j维的分布稀疏性:
λj是样本数据中第j维规范后的标准化偏差;
由此设置各类中心间最小距离的下限θD的初始值:
式中,S为预期类数;
然后直接求样本数据自身第j维的方差
由此设置类内各数据分布的标准差的上限θS的初始值:
④根据自适应准则运行ISODATA算法进行分选;
所述的自适应准则,即通过设置类内各数据分布的标准差的上限和各类中心间最小距离的下限的变换准则,达到自动逼近最佳聚类效果的目的,具体为
设置参数后运行ISODATA算法,
每次聚类的算法结束后,分别计算提取算法结果的类内距离ΘS'和类间距离ΘD',
式中:c'为算法运行产生的分类数,Cp为归类数据的集合,vi为第i类的中心,xn为属于Cp类的数据,vj为与vi不同的类的中心
第一次聚类后仅改变θS,设定自适应原则如下
更改参数后再次运行ISODATA算法
并设定自适应的准则如下:
而退出自适应迭代的条件为:
ΔS=|θs'-θs|=0
ΔD=|θD'-θD|=0
式中:ΔS、ΔD为算法前后连续两次运行结果的类内距离以及类间距离的差值。
2.根据权利要求1所述的雷达辐射源信号分选方法,其特征在于,所述步骤①中的信号频谱具体为:接收到的信号的进行快速傅里叶变换后生成的谱线图。
3.根据权利要求1所述的雷达辐射源信号分选方法,其特征在于,所述步骤①中的三角形信号和矩形信号具体为:选取三角形信号和矩形信号的长度应与信号频谱的长度匹配。
4.根据权利要求1所述的雷达辐射源信号分选方法,其特征在于,所述步骤②中的搜索直方图峰值统计具体为:将一个向量矩阵XN×M中的各个元素在图中用直条的形式描述出来,若向量矩阵的列数为1,得到的直方图就是一个二维图;列数为2,得到的直方图就是一个三维图,进而将直方图推广到向量矩阵列数为M(M≥3)的情况,得到的直方图就是一个M+1的图,进而可以得到区间统计值的直条高度,将每一个划分好区间中统计值与相邻区间进行比较,若该区间的统计值大于其所有相邻区间的统计值,则记该区间为一个峰值区间;最后统计峰值区间的数目即为峰值数目。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710539573.5A CN107238817B (zh) | 2017-07-04 | 2017-07-04 | 一种参数自适应设置及自动调整的雷达辐射源信号分选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710539573.5A CN107238817B (zh) | 2017-07-04 | 2017-07-04 | 一种参数自适应设置及自动调整的雷达辐射源信号分选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107238817A CN107238817A (zh) | 2017-10-10 |
CN107238817B true CN107238817B (zh) | 2020-03-03 |
Family
ID=59991398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710539573.5A Active CN107238817B (zh) | 2017-07-04 | 2017-07-04 | 一种参数自适应设置及自动调整的雷达辐射源信号分选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107238817B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107967490A (zh) * | 2017-12-08 | 2018-04-27 | 中国船舶重工集团公司第七二四研究所 | 一种基于直方图极小值点的脉冲描述字自适应聚类方法 |
CN108535704B (zh) * | 2018-04-10 | 2022-01-21 | 贵州理工学院 | 一种基于自适应二维聚类的信号预分选方法 |
CN109270497B (zh) * | 2018-10-28 | 2023-03-31 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 雷达脉冲信号的多维参数预分选方法 |
CN109670255B (zh) * | 2018-12-26 | 2020-04-07 | 北京华大九天软件有限公司 | 一种时序参数聚类的典型仿真条件推荐方法 |
CN109683143A (zh) * | 2019-03-07 | 2019-04-26 | 西安电子科技大学 | 雷达信号分选方法、装置、计算机设备及可存储介质 |
CN109839619B (zh) * | 2019-03-15 | 2020-10-16 | 北京应用物理与计算数学研究所 | 基于自适应分桶的雷达信号粗分选方法、***及存储介质 |
CN110806563B (zh) * | 2019-11-19 | 2021-07-06 | 西南交通大学 | 基于雷达脉冲混叠程度判定的辐射源信号聚类分选方法 |
CN111832632B (zh) * | 2020-06-28 | 2024-01-23 | 中国人民解放军战略支援部队信息工程大学 | 基于高次频谱对称Holder系数的雷达信号分选方法及*** |
CN112507855A (zh) * | 2020-12-04 | 2021-03-16 | 国网浙江省电力有限公司武义县供电公司 | 一种基于瞬时包络等势星球图的通信辐射源个体识别方法 |
CN112733613A (zh) * | 2020-12-18 | 2021-04-30 | 国网浙江省电力有限公司武义县供电公司 | 基于希尔伯特变换与赫尔德系数特征的辐射源识别方法 |
CN112986925B (zh) * | 2021-02-01 | 2024-03-12 | 东南大学 | 一种基于图像特征的雷达脉冲序列分选方法 |
CN113033629B (zh) * | 2021-03-09 | 2022-08-05 | 中南大学 | 基于改进布谷鸟算法的雷达信号分选方法及装置 |
CN117423338B (zh) * | 2023-12-18 | 2024-03-08 | 卓世未来(天津)科技有限公司 | 一种数字人交互对话方法和*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0732668A1 (fr) * | 1995-03-17 | 1996-09-18 | Alcatel | Méthode d'extraction de contours utilisant une analyse multi-fractale |
CN101902648A (zh) * | 2002-04-19 | 2010-12-01 | 德罗普莱特科技公司 | 小波变换***,方法和计算机程序产品 |
CN103489005A (zh) * | 2013-09-30 | 2014-01-01 | 河海大学 | 一种基于多分类器融合的高分辨率遥感影像分类方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993377B2 (en) * | 2002-02-22 | 2006-01-31 | The Board Of Trustees Of The University Of Arkansas | Method for diagnosing heart disease, predicting sudden death, and analyzing treatment response using multifractal analysis |
-
2017
- 2017-07-04 CN CN201710539573.5A patent/CN107238817B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0732668A1 (fr) * | 1995-03-17 | 1996-09-18 | Alcatel | Méthode d'extraction de contours utilisant une analyse multi-fractale |
CN101902648A (zh) * | 2002-04-19 | 2010-12-01 | 德罗普莱特科技公司 | 小波变换***,方法和计算机程序产品 |
CN103489005A (zh) * | 2013-09-30 | 2014-01-01 | 河海大学 | 一种基于多分类器融合的高分辨率遥感影像分类方法 |
Non-Patent Citations (3)
Title |
---|
Recognition of Radar Emitter Signals Based on SVD and AF Main Ridge Slice;Qiang Guo et al.;《JOURNAL OF COMMUNICATIONS AND NETWORKS》;20151031;第491-498页 * |
基于ISODATA聚类的复杂体制雷达信号分选;张洪亮 杨承志;《电子信息对抗技术》;20101130;第34-38页 * |
基于对称Holder系数的雷达辐射源信号特征分析;王海华 沈晓峰;《空间电子技术》;20101231;第55-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107238817A (zh) | 2017-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107238817B (zh) | 一种参数自适应设置及自动调整的雷达辐射源信号分选方法 | |
CN106371610B (zh) | 一种基于脑电信号的驾驶疲劳的检测方法 | |
CN107515895B (zh) | 一种基于目标检测的视觉目标检索方法与*** | |
CN108596362B (zh) | 基于自适应分段聚合近似的电力负荷曲线形态聚类方法 | |
CN109583436A (zh) | 一种基于毫米波雷达的手势识别*** | |
CN105913081B (zh) | 基于改进的PCAnet的SAR图像分类方法 | |
CN102236901B (zh) | 基于图论聚类和色彩不变空间的目标跟踪方法 | |
CN107220977B (zh) | 基于模糊聚类的有效性指标的图像分割方法 | |
CN110992341A (zh) | 一种基于分割的机载LiDAR点云建筑物提取方法 | |
CN104200114B (zh) | 流式细胞仪数据快速分析方法 | |
CN113537316B (zh) | 一种基于4d毫米波雷达点云的车辆检测方法 | |
CN102663447B (zh) | 基于判别相关分析的跨媒体检索方法 | |
CN102289673A (zh) | 基于偏最小二乘法的高光谱遥感图像波段选择方法 | |
CN108846346B (zh) | 基于密度峰值聚类算法的神经元峰电位在线分类方法 | |
CN105203446B (zh) | 基于概率分布细胞分类统计方法 | |
CN112462355B (zh) | 一种基于时频三特征提取的对海目标智能检测方法 | |
CN113050797A (zh) | 一种通过毫米波雷达实现手势识别的方法 | |
CN111401168A (zh) | 一种无人机的多层雷达特征提取与选择方法 | |
CN105930859B (zh) | 基于线性流形聚类的雷达信号分选方法 | |
CN109001702B (zh) | 无载波超宽带雷达人体动作识别方法 | |
CN109102022B (zh) | 一种基于支持向量机的光变曲线分类方法 | |
CN103942415B (zh) | 一种流式细胞仪数据自动分析方法 | |
CN110969121A (zh) | 一种基于深度学习的高分辨率雷达目标识别算法 | |
CN103218623B (zh) | 基于自适应近邻保持鉴别投影的雷达目标特征提取方法 | |
Dou et al. | Aircraft recognition in high resolution SAR images using saliency map and scattering structure features |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20200207 Address after: 264001 Research and Academic Department, 188 Erma Road, Zhifu District, Yantai City, Shandong Province Applicant after: Naval Aviation University of PLA Address before: 264001 Yantai City, Zhifu Province, No. two road, No. 188, Department of research, Applicant before: Naval Aeronautical Engineering Institute PLA |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |