CN114879139B - 一种车载4d毫米波雷达的联合测角方法、装置及相关设备 - Google Patents
一种车载4d毫米波雷达的联合测角方法、装置及相关设备 Download PDFInfo
- Publication number
- CN114879139B CN114879139B CN202210818690.6A CN202210818690A CN114879139B CN 114879139 B CN114879139 B CN 114879139B CN 202210818690 A CN202210818690 A CN 202210818690A CN 114879139 B CN114879139 B CN 114879139B
- Authority
- CN
- China
- Prior art keywords
- target
- dimensional
- array
- angle
- row
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
-
- 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/415—Identification of targets based on measurements of movement associated with the target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请公开了一种车载4D毫米波雷达的联合测角方法、装置及相关设备,该方法包括:基于目标二维收发阵列中的回波数据矩阵,确定二维测角矩阵,并对其各行进行FFT或DBF处理,得到水平维度压缩后的二维测角矩阵;从该水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于该目标行确定各目标在水平方向的角度信息和索引信息;最后基于索引信息确定目标列,并利用超分辨率算法对目标列进行竖直方向的测角,得到竖直方向的角度信息。本申请无需对二维测角矩阵中的所有行以及所有列进行FFT或DBF,只需要二维测角矩阵中存在至少一行数目足够多的阵元,则可以通过合适的收发天线阵元排布,实现较为有效的角度测量。
Description
技术领域
本申请涉及毫米波雷达技术领域,更具体地说,是涉及一种车载4D毫米波雷达的联合测角方法、装置及相关设备。
背景技术
汽车行业正朝着更智能的方向发展,从高级辅助驾驶系统(Advanced DriverAssistance System,ADAS)到面向更高级别的自动驾驶,这些功能的实现离不开为决策系统提供周边环境信息感知的各类传感器。毫米波雷达有着探测距离远、可以测速、可以全天时全天候工作、成本低等优势,已经成为了汽车最为核心的传感器之一。
传统的毫米波雷达传感器只具备水平维度的测角功能,自动驾驶等级的提高对传感器的要求也更高,于是催生了近几年开始走向市场的4D毫米波雷达。4D雷达可以感知的目标物包括距离、速度、水平角度、竖直角度这四个维度的信息。相较于传统的车载4D毫米波雷达,4D雷达采用了二维面阵,具备了竖直向的测角能力,同时两个维度的测角分辨率也更高。
由于车载4D毫米波雷达是安装于汽车上的,其尺寸不能过大,这与4D雷达要求的极高的测角分辨率(一般要求水平向和竖直向的测角分辨率在1°左右)是相矛盾的:更高角分辨率的获得需要更大的天线面板尺寸。因此,如何利用有限面板尺寸的天线来对目标物进行角度测量,并使其满足角分辨率要求是一个亟待解决的问题。
发明内容
有鉴于此,本申请提供了一种车载4D毫米波雷达的联合测角方法、装置及相关设备,以实现在面板尺寸受限的前提下进行有效测角。
为实现上述目的,本申请第一方面提供了一种车载4D毫米波雷达的联合测角方法,包括:
基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向上的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
优选地,所述基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵的过程,包括:
对各回波数据矩阵进行二维快速傅里叶变换2D-FFT,得到各距离-多普勒谱;
对所述距离-多普勒谱进行非相关累积处理,得到累积后的距离-多普勒谱;
对各距离-多普勒谱进行恒虚警率CFAR检测,得到各目标的距离和速度索引;
基于所述各距离-多普勒谱以及所述各目标的距离和速度索引,获取各目标对应的二维测角矩阵。
优选地,所述利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息的过程,包括:
利用迭代自适应算法IAA或正交匹配追踪算法OMP对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
优选地,基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向的索引信息的过程,包括:
基于所述目标行,利用找极大值法确定各目标在水平方向的角度信息以及在水平方向上的索引信息。
优选地,所述目标二维收发阵列为非均匀天线阵列,包括若干作为发射天线的发射阵元以及若干作为接收天线的接收阵元;
优选地,所述目标二维收发阵列中,一半以上的发射阵元设置在同一行,各接收阵元排列成两行,且每行中相邻的接收阵元之间的间隔一致。
优选地,所述目标二维收发阵列包括12个发射阵元和16个接收阵元;所述12个发射阵元排列成4行,其中6个发射阵元设置在第2行或第3行,另外6个发射阵元均匀排列在另外3行;所述两行接收阵元分别设置在发射阵元的上、下两侧。
本申请第二方面提供了一种车载4D毫米波雷达的联合测角装置,包括:
数据矩阵确定单元,用于基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
数据矩阵压缩单元,用于对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
水平角度确定单元,用于从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向的角度信息以及在水平方向的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
竖直角度确定单元,基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
本申请第三方面提供了一种车载4D毫米波雷达的联合测角设备,包括:存储器和处理器;
所述存储器,用于存储程序;
所述处理器,用于执行所述程序,实现上述的车载4D毫米波雷达的联合测角方法的各个步骤。
本申请第四方面提供了一种存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时,实现如上述的车载4D毫米波雷达的联合测角方法的各个步骤。
经由上述的技术方案可知,本申请基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵。其中,所述二维测角矩阵包含各目标的位置和速度信息。接着,从所述二维测角矩阵中确定包含实际阵元的数量最多的目标行。其中,每一所述阵元对应有一发射天线和一接收天线。由于所述目标行包含阵元的数量最多,意味着所述目标行包含足够的水平方向上的角度信息。从而,通过对所述目标行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,可以得到各目标在水平方向的角度信息以及在水平方向的索引信息。然后,基于每一目标的所述索引信息,从所述二维测角矩阵中确定所述目标的目标列,最后利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。本申请无需对二维测角矩阵中的所有行以及所有列进行FFT或DBF,且对每一行、每一列的元素稀疏度要求不高,只需要二维测角矩阵中存在至少一行数目足够多的阵元,则可以通过合适的收发天线阵元排布,实现较为有效的角度测量,并且满足角度分辨率的要求。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1示例了本申请实施例公开的非均匀排布的二维收发阵列的示意图;
图2示例了本申请实施例公开的二维虚拟阵列的示意图;
图3为本申请实施例公开的车载4D毫米波雷达的联合测角方法的示意图;
图4示例了本申请实施例公开的水平方向上测角的效果示意图;
图5示例了本申请实施例公开的竖直方向上测角的效果示意图;
图6示例了采用传统的二维FFT进行测角的效果示意图;
图7示例了采用本申请的方法进行联合测角的效果示意图;
图8为本申请实施例公开的车载4D毫米波雷达的联合测角装置的示意图;
图9为本申请实施例公开的车载4D毫米波雷达的联合测角设备的示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
为便于理解,下面首先介绍本申请提供的车载4D毫米波雷达的联合测角方法所采用的二维收发阵列。请参阅图1,本申请实施例提供的二维收发阵列采用了非均匀的面阵设计,其包含了12个发射天线以及16个接收天线,该二维收发阵列以多发多收(Multiple-Input Multiple-Output,MIMO)的工作模式应用于车载4D毫米波雷达系统中。其中,多发多收是指雷达的天线面板使用按照一定的设计排布的多个发射天线和多个接收天线来完成电磁波信号的发射和接收,这种方法能以产生虚拟孔径的方式实现更高的分辨率,以实现天线孔径的拓展。
如图1所示,该16个接收阵元平均分成两组,每组接收阵元为间隔的均匀排布,分别设置于面阵的第1行与第6行,两行之间间距为。该12个发射天线设置于面阵的第2至第5行,这4行之间的竖直方向上的位置关系满足最小冗余阵列设计原则。具体地,这4行之间的竖直方向上的位置关系应该是满足(0,1,4,6)的比例关系,对应到以竖直向的最小值来说,这4行之间的间隔即为:,以及。因此,在不考虑天线尺寸的前提下,如图1所示的二维收发阵列所占尺寸为:水平方向上,竖直方向上。以车载毫米波雷达行业规定的77GHz工作频率计算,该二维收发阵列对应的面板尺寸约为:99.43mm*67.79mm,符合对车载毫米波雷达尺寸的要求。
该二维收发阵列可以产生的二维虚拟阵列如图2所示。具体地,MIMO模式下,12发16收可以产生12*16=192个虚拟阵元,但是在该二维收发阵列中,第3行和第7行的最中间阵元会发射重合(对应图2中的两个使用黑色填充的阵元),于是实际产生的非重合阵元只有190个。
下面介绍本申请实施例提供的车载4D毫米波雷达的联合测角方法。请参阅图3,本申请实施例提供的车载4D毫米波雷达的联合测角方法可以包括如下步骤:
步骤S101,基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵。
其中,发射天线使用时分的模式或其他模式发射线性调频信号(chirp信号),单帧下每个天线发射个chirp,接收天线对每个chirp的采样点数为个,因此,每一虚拟阵元对应有1个大小的回波数据矩阵,每一接收天线对应有个大小的回波数据矩阵,其中,为发射天线的个数。
示例性地,对应于图1中的二维收发阵列,每一接收阵元对应有12个,对应的192个虚拟阵元总共可得到192个大小的回波数据矩阵,其中舍弃前述对应位置重复的两个回波数据矩阵,可以得到有效的190个回波数据矩阵。
对各回波数据矩阵分别基于距离、速度的维度进行FFT,接着基于天线维度进行累积,最后通过CFAR检测等处理可以得到二维测角矩阵。可以理解的是,该二维测角矩阵包含各目标的位置和速度信息。
步骤S102,对该二维测角矩阵的各行进行FFT或DBF处理,得到水平维度压缩后的二维测角矩阵。
其中,通过该二维测角矩阵的各行进行FFT(fast Fourier transform,快速傅里叶变换)或DBF(Digital Beam Forming,数字波束成形)处理,实现了从行的维度对该二维测角矩阵进行压缩,经过压缩后的该水平维度压缩后的二维测角矩阵包含了各目标在水平方向上的角度信息以及各目标在水平方向的索引信息。
步骤S103,从水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于该目标行确定各目标在水平方向的角度信息以及在水平方向的索引信息。
可以理解的是,每一阵元对应有一发射天线和一接收天线,在接收阵元均匀排布的前提下,包含阵元数量最多的行,意味着该行所对应的发射阵元数量最多,包含水平方向上的更多的角度信息,以该行作为目标行,即可满足水平方向上的角度计算。对应于如图2所示的二维虚拟阵列,可以将第3行或第7行确定为目标行。
基于该目标行,通过找极大值等方法可以得到各目标在水平方向的角度信息以及在水平方向的索引信息。
步骤S104,基于每一目标的索引信息,从水平维度压缩后的二维测角矩阵中确定该目标的目标列,并利用超分辨率算法对各目标的目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
在车载雷达的实际应用中,一般不使用超分辨率算法做角度测量,主要有两个原因:一是包括经典的MUSIC(Multiple Signal Classification,多重信号分类算法)、ESPRIT(Estimating Signal Parameter via Rotational Invariance Techniques,基于旋转不变技术的信号参数估计)等算法在内的超分辨算法需要多快拍的数据(对应到车载毫米波雷达中就是需要多帧的数据),二是基本所有的超分辨算法都需要较大的计算量。考虑到实际应用场景中雷达和目标都在运动,多快拍下目标与雷达的相对位置关系早已发生改变,利用多快拍的数据进行角度测量会导致较大的偏差。此外车载毫米波雷达对测量结果的更新率有要求(一般为20Hz),所以超分辨算法很少用于车载毫米波雷达中的角度测量。
而本申请实施例在竖直向使用超分辨算法进行测角的依据包括:
1、选用的超分辨算法只基于单快拍(单帧)的数据。
2、不需要对全部的列进行超分辨运算,而只是对水平向存在目标的角度位置进行超分辨率运算。
一般而言,因为上述数据是从单个距离和速度门下得到的,可以通过设置带宽等参数使距离分辨率和速度分辨率足够小,可以认为同一个距离和速度门下最多只允许有4个左右目标,也即需要在竖直方向做超分辨测角的列不会超过4列。
此外,从阵列设计的角度来说,在竖直向使用超分辨测角,在达到同样的角度分辨率的前提下,超分辨算法相较于FFT或波束成形的算法,其对阵列稀疏度的要求更低,也即允许使用较少的阵元而能达到更高的角分辨率,更加有利于满足车载雷达的尺寸要求。
本申请基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵。其中,所述二维测角矩阵包含各目标的位置和速度信息。接着,从所述二维测角矩阵中确定包含实际阵元的数量最多的目标行。其中,每一所述阵元对应有一发射天线和一接收天线。由于所述目标行包含阵元的数量最多,意味着所述目标行包含足够的水平方向上的角度信息。从而,通过对所述目标行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,可以得到各目标在水平方向的角度信息以及在水平方向的索引信息。然后,基于每一目标的所述索引信息,从所述二维测角矩阵中确定所述目标的目标列,最后利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。本申请无需对二维测角矩阵中的所有行以及所有列进行FFT或DBF,且对每一行、每一列的元素稀疏度要求不高,只需要二维测角矩阵中存在至少一行数目足够多的阵元,则可以通过合适的收发天线阵元排布,实现较为有效的角度测量,并且满足角度分辨率的要求。
可以理解的是,本申请上述实施例公开的车载4D毫米波雷达的联合测角方法适用于所有的二维收发阵列,图1所示的二维收发阵列是较优的一种。二维收发阵列中接收阵元、发射阵元之间的具体间距,主要取决于想要达到的最大不模糊测角范围和角度分辨率,当最大不模糊测角范围这一系统参数确定后,相邻阵元之间的最小间隔就会被确定:比如给出的阵列中,水平方向的最大不模糊范围设置为,此时依据前述的雷达的最大不模糊测角范围公式,得到水平方向的阵元的最小间隔应该为。
基于此,在本申请的一些实施例中,目标二维收发阵列可以为非均匀天线阵列,包括若干作为发射天线的发射阵元以及若干作为接收天线的接收阵元。在该目标二维收发阵列中,水平方向上的发射阵元或接收阵元的间距最小值为,竖直方向上的发射阵元或接收阵元的间距最小值为,其中,为所采用的毫米波的波长。
为贴合本申请实施例提供的车载4D毫米波雷达的联合测角方法,在水平方向保证使用FFT或波束成形进行测角的准确性,在二维收发阵列尺寸受限的情况下,可以将一半的发射天线设置在同一行,使其产生如图2中第3行和第7行那样足够密的虚拟阵元。
基于此,在本申请的一些实施例中,在目标二维收发阵列中,一半以上的发射阵元设置在同一行,各接收阵元排列成两行,且每行中相邻的接收阵元之间的间隔一致。
而竖直方向由于可以使用超分辨算法,对阵列稀疏度的要求低一些,可以对发射阵元按照最小冗余阵的间隔进行排布,实际的排布方式可以有多种,但对应的测角分辨率可能不一样。
基于此,在本申请的一些实施例中,目标二维收发阵列可以包括12个发射阵元和16个接收阵元。其中,该12个发射阵元排列成4行,其中6个发射阵元设置在第2行或第3行,另外6个发射阵元均匀排列在另外3行。两行接收阵元分别设置在发射阵元的上、下两侧。
在本申请的一些实施例中,上述步骤S101基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵的过程,可以包括:
S1,对各回波数据矩阵进行2D-FFT,得到各距离-多普勒谱。
其中,通过对各回波数据矩阵进行2D-FFT(2-Dimention Fast FourierTransform,二维快速傅里叶变换),可以完成回波数据在距离维度、速度维度的压缩。
S2,对距离-多普勒谱进行非相关累积处理,得到累积后的距离-多普勒谱。
其中,该非相关累积处理实际上是对S1中得到的数据基于天线维度进行累积,从而得到每个通道下的数据。
S3,对各距离-多普勒谱进行CFAR(Constant False Alarm Rate,恒虚警率)检测,得到各目标的距离和速度索引。
具体地,根据非相干积累结果进行恒虚警率检测,可以得到这一帧内所有可能的目标。根据CFAR检测的可能目标在距离-多普勒谱上的对应位置,可以计算出目标的距离和速度信息。
S4,基于各距离-多普勒谱以及各目标的距离和速度索引,获取各目标对应的二维测角矩阵。
示例性地,对应于图1中的二维收发阵列,可以得到大小为18*87的二维测角矩阵。其中,由于图2中所示的二维虚拟阵列实际上只有190个有效阵元,因此,对于二维测角矩阵中阵元空缺的元素可以采用0来填充。
在本申请的一些实施例中,上述步骤S103基于该目标行确定各目标在水平方向上的角度信息以及在水平方向的索引信息的过程,可以包括:
基于该目标行,利用找极大值法确定各目标在水平方向的角度信息以及在水平方向上的索引信息。
在本申请的一些实施例中,上述步骤S104利用超分辨率算法对各目标的目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息的过程,可以包括:
利用IAA(Iterative Adaptive Approach,迭代自适应算法)或OMP(OrthogonalMatching Pursuit,正交匹配追踪算法)对各目标的目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
为验证本申请上述实施例提供的车载4D毫米波雷达的联合测角方法的有效性,在水平方向上设置了3个目标,目标分别位于0°、1.15°和58°。其中,设置两个目标间距为1.15°是为了验证二维收发阵列在水平方向上有能力分辨出来角度间距为1.15°的两个目标,从而验证水平方向上的角分辨率;设置另外一个目标在58°的位置是为了验证二维收发阵列在水平方向可以检测到58°的目标,从而验证水平方向上的最大不模糊测角范围。
其中,D为天线孔径,此处等于,为载频对应的波长,此处取77GHz对应的波长,为方位角,不同方位角下角分辨率不一样,一般取=0,得到阵列的角分辨率,于是由公式可知所设计阵列的水平向角分辨率约为1.15°,如图4所示,仿真的结果显示设置的位于0°和1.15°处的目标可以被分开,验证了阵列设计的正确性。
竖直方向上设置了3个目标,目标分别位于0°、1.2°和18°,如图1所示的二维收发阵列中,竖直向的阵列孔径大小为,不使用超分辨算法时,结合前面的公式,可得其理论的角分辨率约为2.32°,而从图5所示的仿真结果来看,位于0°和1.2°两处的目标可以被区分开,也即实现了竖直方向的超分辨,经多次仿真实验验证,在图1所示的二维收发阵列下,竖直方向可以达到约1.2°的角分辨率。
此外,为体现本申请实施例提供的车载4D毫米波雷达的联合测角方法与传统的2D-FFT的测角方法效果对比,对随机选取的同一组4个目标进行测角,其中,这些目标分别位于:(-25°,12°)、(6°,8°)、(0°,0°)、(20°,-5°),两种方法的二维处理结果分别如图6和图7所示。
可以看到两种方法都可以完成对设置目标的探测:对应二维平面中四个最大的极值点。但是传统的方法由于阵列的稀疏度较高,导致其它位置处的旁瓣很高,实际应用中一些发射能量较弱的目标极易被较高的旁瓣盖住,在找目标的过程中如果方法不当很容易导致目标找错。而本申请实施例提供的车载4D毫米波雷达的联合测角方法可以将旁瓣压制到极低的水平,二维的测角效果得到了极大改善。
下面对本申请实施例提供的车载4D毫米波雷达的联合测角装置进行描述,下文描述的车载4D毫米波雷达的联合测角装置与上文描述的车载4D毫米波雷达的联合测角方法可相互对应参照。
请参见图8,本申请实施例提供的车载4D毫米波雷达的联合测角装置,可以包括:
数据矩阵确定单元21,用于基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
数据矩阵压缩单元22,用于对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
水平角度确定单元23,用于从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向的角度信息以及在水平方向的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
竖直角度确定单元24,基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
在本申请的一些实施例中,数据矩阵确定单元21基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵的过程,可以包括:
对各回波数据矩阵进行二维快速傅里叶变换2D-FFT,得到各距离-多普勒谱;
对所述距离-多普勒谱进行非相关累积处理,得到累积后的距离-多普勒谱;
对各距离-多普勒谱进行恒虚警率CFAR检测,得到各目标的距离和速度索引;
基于所述各距离-多普勒谱以及所述各目标的距离和速度索引,获取各目标对应的二维测角矩阵。
在本申请的一些实施例中,水平角度确定单元23基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向的索引信息的过程,可以包括:
基于所述目标行,利用找极大值法确定各目标在水平方向的角度信息以及在水平方向上的索引信息。
在本申请的一些实施例中,竖直角度确定单元24利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息的过程,可以包括:
利用迭代自适应算法IAA或正交匹配追踪算法OMP对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
在本申请的一些实施例中,所述目标二维收发阵列为非均匀天线阵列,包括若干作为发射天线的发射阵元以及若干作为接收天线的接收阵元;
在本申请的一些实施例中,所述目标二维收发阵列中,一半以上的发射阵元设置在同一行,各接收阵元排列成两行,且每行中相邻的接收阵元之间的间隔一致。
在本申请的一些实施例中,所述目标二维收发阵列包括12个发射阵元和16个接收阵元;所述12个发射阵元排列成4行,其中6个发射阵元设置在第2行或第3行,另外6个发射阵元均匀排列在另外3行;所述两行接收阵元分别设置在发射阵元的上、下两侧。
本申请实施例提供的车载4D毫米波雷达的联合测角装置可应用于车载4D毫米波雷达的联合测角设备,如计算机等。可选的,图9示出了车载4D毫米波雷达的联合测角设备的硬件结构框图,参照图9,车载4D毫米波雷达的联合测角设备的硬件结构可以包括:至少一个处理器31,至少一个通信接口32,至少一个存储器33和至少一个通信总线34。
在本申请实施例中,处理器31、通信接口32、存储器33、通信总线34的数量为至少一个,且处理器31、通信接口32、存储器33通过通信总线34完成相互间的通信;
处理器31可能是一个中央处理器CPU,或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit),或者是被配置成实施本申请实施例的一个或多个集成电路等;
存储器32可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatilememory)等,例如至少一个磁盘存储器;
其中,存储器33存储有程序,处理器31可调用存储器33存储的程序,所述程序用于:
基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向上的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
可选的,所述程序的细化功能和扩展功能可参照上文描述。
本申请实施例还提供一种存储介质,该存储介质可存储有适于处理器执行的程序,所述程序用于:
基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向上的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
可选的,所述程序的细化功能和扩展功能可参照上文描述。
综上所述:
本申请基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵。其中,所述二维测角矩阵包含各目标的位置和速度信息。接着,从所述二维测角矩阵中确定包含实际阵元的数量最多的目标行。其中,每一所述阵元对应有一发射天线和一接收天线。由于所述目标行包含阵元的数量最多,意味着所述目标行包含足够的水平方向上的角度信息。从而,通过对所述目标行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,可以得到各目标在水平方向的角度信息以及在水平方向的索引信息。然后,基于每一目标的所述索引信息,从所述二维测角矩阵中确定所述目标的目标列,最后利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。本申请无需对二维测角矩阵中的所有行以及所有列进行FFT或DBF,且对每一行、每一列的元素稀疏度要求不高,只需要二维测角矩阵中存在至少一行数目足够多的阵元,则可以通过合适的收发天线阵元排布,实现较为有效的角度测量,并且满足角度分辨率的要求。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间可以根据需要进行组合,且相同相似部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (7)
1.一种车载4D毫米波雷达的联合测角方法,其特征在于,包括:
基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向上的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息;
所述目标二维收发阵列为非均匀天线阵列,包括若干作为发射天线的发射阵元以及若干作为接收天线的接收阵元;
所述目标二维收发阵列中,一半以上的发射阵元设置在同一行,各接收阵元排列成两行,且每行中相邻的接收阵元之间的间隔一致;
所述目标二维收发阵列包括12个发射阵元和16个接收阵元;所述12个发射阵元排列成4行,其中6个发射阵元设置在第2行或第3行,另外6个发射阵元均匀排列在另外3行;两行接收阵元分别设置在发射阵元的上、下两侧。
2.根据权利要求1所述的方法,其特征在于,所述基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵的过程,包括:
对各回波数据矩阵进行二维快速傅里叶变换2D-FFT,得到各距离-多普勒谱;
对所述距离-多普勒谱进行非相关累积处理,得到累积后的距离-多普勒谱;
对各距离-多普勒谱进行恒虚警率CFAR检测,得到各目标的距离和速度索引;
基于所述各距离-多普勒谱以及所述各目标的距离和速度索引,获取各目标对应的二维测角矩阵。
3.根据权利要求1所述的方法,其特征在于,所述利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息的过程,包括:
利用迭代自适应算法IAA或正交匹配追踪算法OMP对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息。
4.根据权利要求1所述的方法,其特征在于,基于所述目标行确定各目标在水平方向上的角度信息以及在水平方向的索引信息的过程,包括:
基于所述目标行,利用找极大值法确定各目标在水平方向的角度信息以及在水平方向上的索引信息。
5.一种车载4D毫米波雷达的联合测角装置,其特征在于,包括:
数据矩阵确定单元,用于基于目标二维收发阵列中各接收天线接收到的回波数据矩阵,确定各目标对应的二维测角矩阵,所述二维测角矩阵包含各目标的位置和速度信息;
数据矩阵压缩单元,用于对所述二维测角矩阵的各行进行快速傅里叶变换FFT或进行数字波束成形DBF处理,得到水平维度压缩后的二维测角矩阵;
水平角度确定单元,用于从所述水平维度压缩后的二维测角矩阵中确定包含实际阵元的数量最多的目标行,并基于所述目标行确定各目标在水平方向的角度信息以及在水平方向的索引信息,其中,每一所述阵元对应有一发射天线和一接收天线;
竖直角度确定单元,基于每一目标的所述索引信息,从所述水平维度压缩后的二维测角矩阵中确定所述目标的目标列,并利用超分辨率算法对各目标的所述目标列进行竖直方向上的测角,得到各目标在竖直方向的角度信息;
所述目标二维收发阵列为非均匀天线阵列,包括若干作为发射天线的发射阵元以及若干作为接收天线的接收阵元;
所述目标二维收发阵列中,一半以上的发射阵元设置在同一行,各接收阵元排列成两行,且每行中相邻的接收阵元之间的间隔一致;
所述目标二维收发阵列包括12个发射阵元和16个接收阵元;所述12个发射阵元排列成4行,其中6个发射阵元设置在第2行或第3行,另外6个发射阵元均匀排列在另外3行;两行接收阵元分别设置在发射阵元的上、下两侧。
6.一种车载4D毫米波雷达的联合测角设备,其特征在于,包括:存储器和处理器;
所述存储器,用于存储程序;
所述处理器,用于执行所述程序,实现如权利要求1~4中任一项所述的车载4D毫米波雷达的联合测角方法的各个步骤。
7.一种存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如权利要求1~4中任一项所述的车载4D毫米波雷达的联合测角方法的各个步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210818690.6A CN114879139B (zh) | 2022-07-13 | 2022-07-13 | 一种车载4d毫米波雷达的联合测角方法、装置及相关设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210818690.6A CN114879139B (zh) | 2022-07-13 | 2022-07-13 | 一种车载4d毫米波雷达的联合测角方法、装置及相关设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114879139A CN114879139A (zh) | 2022-08-09 |
CN114879139B true CN114879139B (zh) | 2022-09-23 |
Family
ID=82683228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210818690.6A Active CN114879139B (zh) | 2022-07-13 | 2022-07-13 | 一种车载4d毫米波雷达的联合测角方法、装置及相关设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114879139B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115327473B (zh) * | 2022-10-17 | 2022-12-16 | 广东大湾区空天信息研究院 | 等效多快拍的4d毫米波雷达测角方法、装置及相关设备 |
CN116106847B (zh) * | 2023-01-18 | 2023-11-07 | 珠海微度芯创科技有限责任公司 | 毫米波雷达二维联合超分辨测角方法、装置及存储介质 |
CN116299301B (zh) * | 2023-05-17 | 2023-07-21 | 广东大湾区空天信息研究院 | 一种毫米波雷达rcs实时估测方法、装置及相关设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN111693978A (zh) * | 2020-05-13 | 2020-09-22 | 复旦大学 | 基于mimo毫米波雷达的散点检测方法 |
CN111896947A (zh) * | 2019-10-14 | 2020-11-06 | 镇江盛益系统科技有限公司 | 一种用于汽车毫米波雷达的快速超分辨跟踪系统及方法 |
CN112946582A (zh) * | 2021-04-02 | 2021-06-11 | 华域汽车系统股份有限公司 | 一种基于车载mimo雷达的天线阵列及其使用方法 |
CN113314832A (zh) * | 2021-06-15 | 2021-08-27 | 东南大学 | 一种毫米波车载mimo雷达天线阵列装置及设计方法 |
CN113504522A (zh) * | 2021-07-06 | 2021-10-15 | 北京航空航天大学 | 一种基于发射天线随机切换的时空解耦与超分辨测角方法 |
CN114252881A (zh) * | 2021-12-08 | 2022-03-29 | 广州极飞科技股份有限公司 | 一种目标检测方法、装置、存储介质及终端设备 |
CN114265058A (zh) * | 2021-12-13 | 2022-04-01 | 中汽创智科技有限公司 | Mimo雷达目标测角方法、装置、电子设备及存储介质 |
CN114447631A (zh) * | 2022-01-14 | 2022-05-06 | 深圳市道通智能汽车有限公司 | 一种发射天线阵列、发射天线系统及毫米波雷达 |
CN114609626A (zh) * | 2022-03-16 | 2022-06-10 | 电子科技大学 | 一种车载毫米波雷达多目标检测方法 |
-
2022
- 2022-07-13 CN CN202210818690.6A patent/CN114879139B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111896947A (zh) * | 2019-10-14 | 2020-11-06 | 镇江盛益系统科技有限公司 | 一种用于汽车毫米波雷达的快速超分辨跟踪系统及方法 |
CN110673086A (zh) * | 2019-10-31 | 2020-01-10 | 上海无线电设备研究所 | 一种基于数字阵列雷达的二维角度超分辨方法 |
CN111693978A (zh) * | 2020-05-13 | 2020-09-22 | 复旦大学 | 基于mimo毫米波雷达的散点检测方法 |
CN112946582A (zh) * | 2021-04-02 | 2021-06-11 | 华域汽车系统股份有限公司 | 一种基于车载mimo雷达的天线阵列及其使用方法 |
CN113314832A (zh) * | 2021-06-15 | 2021-08-27 | 东南大学 | 一种毫米波车载mimo雷达天线阵列装置及设计方法 |
CN113504522A (zh) * | 2021-07-06 | 2021-10-15 | 北京航空航天大学 | 一种基于发射天线随机切换的时空解耦与超分辨测角方法 |
CN114252881A (zh) * | 2021-12-08 | 2022-03-29 | 广州极飞科技股份有限公司 | 一种目标检测方法、装置、存储介质及终端设备 |
CN114265058A (zh) * | 2021-12-13 | 2022-04-01 | 中汽创智科技有限公司 | Mimo雷达目标测角方法、装置、电子设备及存储介质 |
CN114447631A (zh) * | 2022-01-14 | 2022-05-06 | 深圳市道通智能汽车有限公司 | 一种发射天线阵列、发射天线系统及毫米波雷达 |
CN114609626A (zh) * | 2022-03-16 | 2022-06-10 | 电子科技大学 | 一种车载毫米波雷达多目标检测方法 |
Non-Patent Citations (3)
Title |
---|
An Efficient Algorithm Based on Frequency Scaling for THz Stepped-Frequency SAR Imaging;Chuang Wang等;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20220321;第60卷;1-15 * |
基于HRRP序列的钝头倒角锥目标微动特性分析及参数估计;苏楠等;《电子与信息学报》;20190731;第41卷(第7期);1751-1757 * |
基于压缩感知的弹药自动装填系统数据采集方法;刘曦等;《南京理工大学学报》;20161031;第40卷(第5期);544-548 * |
Also Published As
Publication number | Publication date |
---|---|
CN114879139A (zh) | 2022-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114879139B (zh) | 一种车载4d毫米波雷达的联合测角方法、装置及相关设备 | |
CN109946665B (zh) | 基于阵列雷达的获取真实目标的方法 | |
KR102445130B1 (ko) | 레이더 표적의 횡방향 상대 속도 성분을 확인하기 위한 방법 및 장치 | |
CN114280593B (zh) | 一种雷达系统及车辆 | |
CN115327473B (zh) | 等效多快拍的4d毫米波雷达测角方法、装置及相关设备 | |
CN111679266B (zh) | 汽车毫米波雷达稀疏阵列栅瓣虚假目标识别方法及系统 | |
JP7027579B2 (ja) | 自動車用のmimoレーダセンサ | |
US20240053469A1 (en) | Method and device to process radar signal | |
CN110865364B (zh) | 雷达的目标解算方法及终端设备 | |
US11327170B2 (en) | Azimuth and elevation radar imaging with single-dimension antenna arrays of radar system | |
CN113109781B (zh) | 波达方向估计方法、雷达和可移动设备 | |
US20200284901A1 (en) | Millimeter-wave radar for unmanned aerial vehicle swarming, tracking, and collision avoidance | |
WO2021129581A1 (zh) | 一种信号处理方法及装置 | |
Gao et al. | Perception through 2d-mimo fmcw automotive radar under adverse weather | |
US20210323560A1 (en) | Vehicle speed calculation method, system, device, and storage medium | |
US8194250B2 (en) | Wideband RF detection | |
CN116359908A (zh) | 点云数据增强方法、装置、计算机设备、系统及存储介质 | |
US11181614B2 (en) | Antenna array tilt and processing to eliminate false detections in a radar system | |
CN116500620A (zh) | 毫米波雷达的数据处理方法、装置、存储介质及无人车 | |
CN111175753A (zh) | 一种车载防撞雷达分波目标探测方法 | |
JP7160561B2 (ja) | 方位演算装置及び方位演算方法 | |
US20230375690A1 (en) | Efficient Direction Of Arrival Estimation Using Low Rank Approximation | |
CN114265058A (zh) | Mimo雷达目标测角方法、装置、电子设备及存储介质 | |
CN114879138B (zh) | 一种毫米波雷达二维角度计算方法、装置及相关设备 | |
CN116106847B (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 |