CN109782216B - 一种二维空间谱峰值搜索的简易方法 - Google Patents
一种二维空间谱峰值搜索的简易方法 Download PDFInfo
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明提供了一种二维空间谱峰值搜索的简易方法,获取电磁波信号的二维空间谱并离散化,计算空间谱的每一行中所有列数据的最大值,利用相邻值比较法确认每一个位置数据对是否对应二维谱空间谱的峰值;若为数据峰值,则保留该组数据对;否则,剔除该组数据对;空间谱峰值的位置为保留下来的数据对所对应的角度。本发明根据系统的对响应速度以及方向搜索精确度的具体要求来对Δ进行灵活配置,将复杂且耗时的二维空间的扫描搜索,通过两个一维空间搜索的方式对其进行了简化操作,能够降低信号处理的复杂度,节省搜索时间,提高电磁波探测系统的响应的实时性。
Description
技术领域
本发明属于电磁波空间谱估计领域,涉及一种谱峰值搜索方法。
背景技术
电磁波探测系统通过扫描来获取的电磁信号的空间谱强度。扫描的方式可以仅针对方位或者俯仰进行一维扫描,也可以同时进行方位和俯仰的二维扫描。
对于一维扫描获取的电磁空间谱,用一维相邻值比较法可以确定某点是否为空间谱峰值,具体方法是判断该点的空间谱强度是否同时大于左右相邻点的空间谱强度。对于二维扫描获取的电磁空间谱,用二维相邻值比较法可以确定某点是否为空间谱峰值,具体方法是判断该点的空间谱强度是否同时大于方位上相邻点和俯仰上的相邻点的空间谱强度。
虽然利用二维相邻值比较法进行二维空间谱搜索,可以获取电磁波的来向。但是,二维搜索方法的搜索过程需要遍历所有的空间方向,搜索过程十分耗时,故而会降低电磁波探测系统响应和预警速度。
发明内容
为了克服现有技术的不足,本发明提供一种二维空间谱峰值搜索的简易方法,可以有效的提升搜索速度,提高系统处理的实时性。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤2,对进行离散化,离散化后二维谱空间谱Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°;m和n分别为θ/Δ和数据区间内的整数值,在二维谱空间谱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.本发明将复杂且耗时的二维空间的扫描搜索,通过两个一维空间搜索的方式对其进行了简化操作,能够降低信号处理的复杂度,节省搜索时间,提高电磁波探测系统的响应的实时性。
附图说明
图2是序列Ri的曲线图;
图3是序列Li的曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明通过空间谱估计方法,例如Capon、MUSIC等方法,获得电磁波在不同来向上的空间谱响应;通过扫描计算分析二维空间谱的峰值的位置,获得空间电磁波的来向信息。在电磁探测系统使用过程中,通过空间谱估计方法,对空间电磁波信息的侦察和预警,并为电磁反制提供保障信息。
具体实现如下:
步骤2:对二维空间谱进行离散化,离散化后二维谱空间谱为Q(m,n),Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°。m和n分别为θ/Δ和数据区间内的整数值。在二维谱空间谱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中保留下来的数据对所对应的角度。
本发明的实施例包括以下步骤:
步骤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.一种二维空间谱峰值搜索的简易方法,其特征在于包括下述步骤:
步骤2,对进行离散化,离散化后二维谱空间谱Q(m,n)=P(m·Δ,n·Δ),其中,Δ为角度离散间隔,0°<Δ≤90°;m和n分别为θ/Δ和数据区间内的整数值,在二维谱空间谱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中保留下来的数据对所对应的角度。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111722185B (zh) * | 2020-05-15 | 2023-10-13 | 深圳市微纳感知计算技术有限公司 | 一种特征声音的定位方法、装置及设备 |
CN117805721B (zh) * | 2024-02-29 | 2024-05-07 | 成都大公博创信息技术有限公司 | 车载复杂信号监测测向方法及系统 |
Citations (6)
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的线性调频信号参数估计方法 |
CN107783123A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 无人驾驶汽车复杂环境防撞毫米波雷达信号处理系统及方法 |
CN107783129A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 一种旋翼无人机防撞毫米波雷达信号处理方法 |
-
2018
- 2018-12-26 CN CN201811603867.0A patent/CN109782216B/zh active Active
Patent Citations (6)
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的线性调频信号参数估计方法 |
CN107783123A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 无人驾驶汽车复杂环境防撞毫米波雷达信号处理系统及方法 |
CN107783129A (zh) * | 2016-08-25 | 2018-03-09 | 大连楼兰科技股份有限公司 | 一种旋翼无人机防撞毫米波雷达信号处理方法 |
Non-Patent Citations (4)
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) | 一种基于神经网络的毫米波雷达与视觉协同目标检测与识别方法 | |
CN109993800B (zh) | 一种工件尺寸的检测方法、装置及存储介质 | |
US11443437B2 (en) | Vibe-based three-dimensional sonar point cloud image segmentation method | |
CN101950419B (zh) | 同时存在平移和旋转情况下的快速图像配准方法 | |
CN111781608B (zh) | 一种基于fmcw激光雷达的运动目标检测方法及系统 | |
CN109839624A (zh) | 一种多激光雷达位置标定方法及装置 | |
CN104200461B (zh) | 基于互信息图像选块和sift特征的遥感图像配准方法 | |
WO2016055031A1 (zh) | 直线检测、图像处理的方法及相关装置 | |
CN103292701A (zh) | 基于机器视觉的精密器件在线尺寸测量方法 | |
CN105865344A (zh) | 一种基于机器视觉的工件尺寸测量方法和装置 | |
CN106157219B (zh) | 基于车载扫描系统的道路中心线提取方法和装置 | |
JP7015723B2 (ja) | 物体検出装置、物体検出システム、及び物体検出方法 | |
CN109782216B (zh) | 一种二维空间谱峰值搜索的简易方法 | |
CN105447867B (zh) | 基于isar图像的空间目标姿态估计方法 | |
CN110580705B (zh) | 一种基于双域图信号滤波检测建筑物边缘点的方法 | |
CN106771329B (zh) | 一种无人机减速过程中运行速度的检测方法 | |
CN104050660A (zh) | 一种测量工件圆形边缘的方法 | |
CN108229560B (zh) | 基于轮廓曲线匹配算法实现数控系统工件定位匹配的方法 | |
CN115267718A (zh) | 基于点云拼接的环视雷达360°探测实现方法 | |
CN104463896B (zh) | 基于核相似区分布特性的图像角点检测方法及系统 | |
CN103337080A (zh) | 基于梯度方向Hausdorff 距离的红外图像和可见光图像配准技术 | |
CN111626236A (zh) | 一种快速的椭圆目标检测方法 | |
CN115436966A (zh) | 一种激光雷达参考水深控制点批量提取方法 | |
CN104050659A (zh) | 一种测量工件线性边缘的方法 | |
CN104700401B (zh) | 一种基于K‑Means聚类法的图像仿射变换控制点选取方法 |
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 |