CN109782216B - 一种二维空间谱峰值搜索的简易方法 - Google Patents

一种二维空间谱峰值搜索的简易方法 Download PDF

Info

Publication number
CN109782216B
CN109782216B CN201811603867.0A CN201811603867A CN109782216B CN 109782216 B CN109782216 B CN 109782216B CN 201811603867 A CN201811603867 A CN 201811603867A CN 109782216 B CN109782216 B CN 109782216B
Authority
CN
China
Prior art keywords
dimensional
spectrum
data
sequence
dimensional space
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
CN201811603867.0A
Other languages
English (en)
Other versions
CN109782216A (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.)
CETC 20 Research Institute
Original Assignee
CETC 20 Research Institute
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 CETC 20 Research Institute filed Critical CETC 20 Research Institute
Priority to CN201811603867.0A priority Critical patent/CN109782216B/zh
Publication of CN109782216A publication Critical patent/CN109782216A/zh
Application granted granted Critical
Publication of CN109782216B publication Critical patent/CN109782216B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明提供了一种二维空间谱峰值搜索的简易方法,获取电磁波信号的二维空间谱并离散化,计算空间谱的每一行中所有列数据的最大值,利用相邻值比较法确认每一个位置数据对是否对应二维谱空间谱的峰值;若为数据峰值,则保留该组数据对;否则,剔除该组数据对;空间谱峰值的位置为保留下来的数据对所对应的角度。本发明根据系统的对响应速度以及方向搜索精确度的具体要求来对Δ进行灵活配置,将复杂且耗时的二维空间的扫描搜索,通过两个一维空间搜索的方式对其进行了简化操作,能够降低信号处理的复杂度,节省搜索时间,提高电磁波探测系统的响应的实时性。

Description

一种二维空间谱峰值搜索的简易方法
技术领域
本发明属于电磁波空间谱估计领域,涉及一种谱峰值搜索方法。
背景技术
电磁波探测系统通过扫描来获取的电磁信号的空间谱强度。扫描的方式可以仅针对方位或者俯仰进行一维扫描,也可以同时进行方位和俯仰的二维扫描。
对于一维扫描获取的电磁空间谱,用一维相邻值比较法可以确定某点是否为空间谱峰值,具体方法是判断该点的空间谱强度是否同时大于左右相邻点的空间谱强度。对于二维扫描获取的电磁空间谱,用二维相邻值比较法可以确定某点是否为空间谱峰值,具体方法是判断该点的空间谱强度是否同时大于方位上相邻点和俯仰上的相邻点的空间谱强度。
虽然利用二维相邻值比较法进行二维空间谱搜索,可以获取电磁波的来向。但是,二维搜索方法的搜索过程需要遍历所有的空间方向,搜索过程十分耗时,故而会降低电磁波探测系统响应和预警速度。
发明内容
为了克服现有技术的不足,本发明提供一种二维空间谱峰值搜索的简易方法,可以有效的提升搜索速度,提高系统处理的实时性。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1,获取电磁波信号的二维空间谱
Figure GDA0003699158590000011
其中,θ为俯仰角,
Figure GDA0003699158590000012
为方位角,0°≤θ≤90°,
Figure GDA0003699158590000013
步骤2,对
Figure GDA0003699158590000014
进行离散化,离散化后二维谱空间谱Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°;m和n分别为θ/Δ和
Figure GDA0003699158590000015
数据区间内的整数值,在二维谱空间谱Q(m,n)中,m表示行序号,m的最大值用M表示,n表示列序号,n的最大值用N表示;
步骤3,计算Q(m,n)的每一行中所有列数据的最大值,组成一维序列Ri,0≤i≤M;
步骤4,利用相邻值比较法计算一维序列Ri的极大值点的位置序号,位置序号组成序列X;
步骤5,计算Q(m,n)的每一列中所有行数据的最大值,组成一维序列Lj,0≤j≤N;
步骤6,利用相邻值比较法计算一维序列Lj的极大值点的位置序号,位置序号组成序列Y;
步骤7,将序列X每个元素与序列Y的每一个元素组成一组位置数据对;
步骤8,根据二维相邻值比较法,确认每一个位置数据对是否对应二维谱空间谱为Q(m,n)的峰值;若为数据峰值,则保留该组数据对;否则,剔除该组数据对;
步骤9,空间谱峰值的位置为步骤8中保留下来的数据对所对应的角度。
本发明的有益效果是:
1.本发明可以通过调整Δ值,控制系统的搜索响应速度。Δ值越小则峰值搜索越精确,Δ值越大则搜索速度越快,可以根据系统的对响应速度以及方向搜索精确度的具体要求来对Δ进行灵活配置。
2.本发明将复杂且耗时的二维空间的扫描搜索,通过两个一维空间搜索的方式对其进行了简化操作,能够降低信号处理的复杂度,节省搜索时间,提高电磁波探测系统的响应的实时性。
附图说明
图1是二维空间谱
Figure GDA0003699158590000021
示意图,其中(a)是空间谱的3D显示,(b)是局部空间谱等高线显示;
图2是序列Ri的曲线图;
图3是序列Li的曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明通过空间谱估计方法,例如Capon、MUSIC等方法,获得电磁波在不同来向上的空间谱响应;通过扫描计算分析二维空间谱的峰值的位置,获得空间电磁波的来向信息。在电磁探测系统使用过程中,通过空间谱估计方法,对空间电磁波信息的侦察和预警,并为电磁反制提供保障信息。
具体实现如下:
步骤1:获取电磁波信号的二维空间谱
Figure GDA0003699158590000022
其中,θ为俯仰角,
Figure GDA0003699158590000023
为方位角,0°≤θ≤90°,
Figure GDA0003699158590000024
步骤2:对二维空间谱
Figure GDA0003699158590000025
进行离散化,离散化后二维谱空间谱为Q(m,n),Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°。m和n分别为θ/Δ和
Figure GDA0003699158590000031
数据区间内的整数值。在二维谱空间谱Q(m,n)中,m表示行序号,m的最大值用M表示,n表示列序号,n的最大值用N表示。
步骤3:计算Q(m,n)的每一行中所有列数据的最大值,组成一维序列Ri,0≤i≤M。
步骤4:利用相邻值比较法计算一维序列Ri的极大值点的位置序号,位置序号组成序列X。
步骤5:计算Q(m,n)的每一列中所有行数据的最大值,组成一维序列Lj,0≤j≤N。
步骤6:利用相邻值比较法计算一维序列Lj的极大值点的位置序号,位置序号组成序列Y。
步骤7:将序列X每个元素与序列Y的每一个元素组成一组位置数据对。
步骤8:根据二维相邻值比较法,确认每一个位置数据对是否对应二维谱空间谱为Q(m,n)的峰值。若为数据峰值,则保留该组数据对;否则,剔除该组数据对。
步骤9:空间谱峰值的位置为步骤8中保留下来的数据对所对应的角度。
本发明的实施例包括以下步骤:
步骤1:获取电磁波信号的二维空间谱
Figure GDA0003699158590000032
其中,θ为俯仰角,
Figure GDA0003699158590000033
为方位角,0°≤θ≤90°,
Figure GDA0003699158590000034
见图1所示。
步骤2:对二维空间谱
Figure GDA0003699158590000035
进行离散化,离散化后二维谱空间谱为Q(m,n),Q(m,n)=P(m·Δ,n·Δ),其中,角度离散间隔Δ为1°。m和n为整数,m∈[0,90],n∈[0,359]。
步骤3:计算Q(m,n)的每一行中所有列数据的最大值,组成一维序列Ri,0≤i≤90,如图2所示。
步骤4:利用相邻值比较法计算一维序列Ri的极大值点的位置序号,位置序号组成序列X,X={40,50,70}。
步骤5:计算Q(m,n)的每一列中所有行数据的最大值,组成一维序列Lj,0≤j≤N,如图3所示。
步骤6:利用相邻值比较法计算一维序列Lj的极大值点的位置序号,位置序号组成序列Y,Y={40,120,180}。
步骤7:将序列X每个元素与序列Y的每一个元素组成一组位置数据对,分别为[40,40],[40,120],[40,180],[50,40],[50,120],[50,180],[70,40],[70,120],[70,180]。
步骤8:根据二维相邻值比较法,确认每一个位置数据对是否对应二维谱空间谱为Q(m,n)的峰值。若为数据峰值,则保留该组数据对;否则,剔除该组数据对。经过筛选剔除,保留的数据对分别为[40,120],[50,120],[50,180],[70,40];
步骤9:空间谱峰值的位置分别为(俯仰40°方位120°)、(俯仰50°方位120°)、(俯仰50°方位180°)、(俯仰70°方位40°)。
至此,一种二维空间谱峰值搜索的简易方法的设计完成。

Claims (1)

1.一种二维空间谱峰值搜索的简易方法,其特征在于包括下述步骤:
步骤1,获取电磁波信号的二维空间谱
Figure FDA0003699158580000011
其中,θ为俯仰角,
Figure FDA0003699158580000012
为方位角,0°≤θ≤90°,
Figure FDA0003699158580000013
步骤2,对
Figure FDA0003699158580000014
进行离散化,离散化后二维谱空间谱Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°;m和n分别为θ/Δ和
Figure FDA0003699158580000015
数据区间内的整数值,在二维谱空间谱Q(m,n)中,m表示行序号,m的最大值用M表示,n表示列序号,n的最大值用N表示;
步骤3,计算Q(m,n)的每一行中所有列数据的最大值,组成一维序列Ri,0≤i≤M;
步骤4,利用相邻值比较法计算一维序列Ri的极大值点的位置序号,位置序号组成序列X;
步骤5,计算Q(m,n)的每一列中所有行数据的最大值,组成一维序列Lj,0≤j≤N;
步骤6,利用相邻值比较法计算一维序列Lj的极大值点的位置序号,位置序号组成序列Y;
步骤7,将序列X每个元素与序列Y的每一个元素组成一组位置数据对;
步骤8,根据二维相邻值比较法,确认每一个位置数据对是否对应二维谱空间谱为Q(m,n)的峰值;若为数据峰值,则保留该组数据对;否则,剔除该组数据对;
步骤9,空间谱峰值的位置为步骤8中保留下来的数据对所对应的角度。
CN201811603867.0A 2018-12-26 2018-12-26 一种二维空间谱峰值搜索的简易方法 Active CN109782216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811603867.0A CN109782216B (zh) 2018-12-26 2018-12-26 一种二维空间谱峰值搜索的简易方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811603867.0A CN109782216B (zh) 2018-12-26 2018-12-26 一种二维空间谱峰值搜索的简易方法

Publications (2)

Publication Number Publication Date
CN109782216A CN109782216A (zh) 2019-05-21
CN109782216B true CN109782216B (zh) 2022-08-23

Family

ID=66498432

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811603867.0A Active CN109782216B (zh) 2018-12-26 2018-12-26 一种二维空间谱峰值搜索的简易方法

Country Status (1)

Country Link
CN (1) CN109782216B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111722185B (zh) * 2020-05-15 2023-10-13 深圳市微纳感知计算技术有限公司 一种特征声音的定位方法、装置及设备
CN117805721A (zh) * 2024-02-29 2024-04-02 成都大公博创信息技术有限公司 车载复杂信号监测测向方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6493380B1 (en) * 1999-05-28 2002-12-10 Nortel Networks Limited System and method for estimating signal time of arrival
JP2004309250A (ja) * 2003-04-04 2004-11-04 Shimadzu Corp クロマトグラフ用データ処理装置
CN103116162A (zh) * 2012-12-14 2013-05-22 西安电子科技大学 基于目标空间稀疏性的高分辨声呐定位方法
CN103412287A (zh) * 2013-09-01 2013-11-27 西安电子科技大学 基于lvd的线性调频信号参数估计方法
CN107783129A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 一种旋翼无人机防撞毫米波雷达信号处理方法
CN107783123A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 无人驾驶汽车复杂环境防撞毫米波雷达信号处理系统及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6493380B1 (en) * 1999-05-28 2002-12-10 Nortel Networks Limited System and method for estimating signal time of arrival
JP2004309250A (ja) * 2003-04-04 2004-11-04 Shimadzu Corp クロマトグラフ用データ処理装置
CN103116162A (zh) * 2012-12-14 2013-05-22 西安电子科技大学 基于目标空间稀疏性的高分辨声呐定位方法
CN103412287A (zh) * 2013-09-01 2013-11-27 西安电子科技大学 基于lvd的线性调频信号参数估计方法
CN107783129A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 一种旋翼无人机防撞毫米波雷达信号处理方法
CN107783123A (zh) * 2016-08-25 2018-03-09 大连楼兰科技股份有限公司 无人驾驶汽车复杂环境防撞毫米波雷达信号处理系统及方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An accurate real-time algorithm for spectrum peaks search in 2D MUSIC;Jin-feng Chen 等;《 2011 International Conference on Multimedia Technology》;20110830;全文 *
一种新的DOA估计的高分辨率算法;徐豫西 等;《计算机工程与应用》;20151230;全文 *
一种采用稀疏表示的快速空时自适应方法;解虎 等;《西安电子科技大学学报(自然科学版)》;20151030;第42卷(第5期);全文 *
准实时二维MUSIC谱估计方法研究;吴世龙;《舰船电子工程》;20060820(第04期);全文 *

Also Published As

Publication number Publication date
CN109782216A (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
CN110363158B (zh) 一种基于神经网络的毫米波雷达与视觉协同目标检测与识别方法
CN102032875B (zh) 一种基于图像处理的电缆护套厚度测量方法
US11443437B2 (en) Vibe-based three-dimensional sonar point cloud image segmentation method
CN109993800B (zh) 一种工件尺寸的检测方法、装置及存储介质
CN109655805B (zh) 一种基于扫描线段重合长度估计的激光雷达定位方法
CN104200461B (zh) 基于互信息图像选块和sift特征的遥感图像配准方法
CN105335973B (zh) 运用于带钢加工生产线的视觉处理方法
WO2016055031A1 (zh) 直线检测、图像处理的方法及相关装置
CN103292701A (zh) 基于机器视觉的精密器件在线尺寸测量方法
CN111781608B (zh) 一种基于fmcw激光雷达的运动目标检测方法及系统
JP7015723B2 (ja) 物体検出装置、物体検出システム、及び物体検出方法
CN105865344A (zh) 一种基于机器视觉的工件尺寸测量方法和装置
CN109782216B (zh) 一种二维空间谱峰值搜索的简易方法
JP2012189445A (ja) 物体検出装置及び物体検出方法
CN103727930A (zh) 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法
CN110009680B (zh) 基于圆特征及异面特征点的单目图像位置、姿态测量方法
CN104050660A (zh) 一种测量工件圆形边缘的方法
CN112130142A (zh) 一种复杂运动目标微多普勒特征提取方法及系统
CN106771329B (zh) 一种无人机减速过程中运行速度的检测方法
CN115267718A (zh) 基于点云拼接的环视雷达360°探测实现方法
CN104463896B (zh) 基于核相似区分布特性的图像角点检测方法及系统
CN104050659A (zh) 一种测量工件线性边缘的方法
CN104700401B (zh) 一种基于K‑Means聚类法的图像仿射变换控制点选取方法
CN104142495A (zh) 一种基于频谱的斜视sar点目标插值与剖面截取方法
CN110287514B (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