毫米波FMCW MIMO雷达三维点云成像方法

晋良念， 王 燃

(1. 桂林电子科技大学信息与通信学院， 广西桂林 541004；2. 广西无线宽带通信与信号处理重点实验室, 广西桂林 541004)

0 引言

近年来，自动驾驶、智慧交通、区域监控以及安防等领域的需求飙升，这些领域的相关应用场景要求传感器能够探测环境中目标的空间三维信息。相比于其他感知环境的传感器如相机、激光雷达来说，商业的毫米波雷达有着独特的优势。不仅因其大小适合安装在汽车的各个位置，能很好地适应汽车和工业领域对雷达设备体积的要求，而且探测距离远，对各种恶劣天气环境的鲁棒性强，在复杂工作环境下能正常运行。而毫米波FMCW MIMO雷达高分辨的3D点云成像技术将有助于提高上述应用中目标检测、分类、预警等功能的效率与准确性。

现有毫米波FMCW MIMO雷达3D点云成像方法包括三维快速傅里叶变换(3D FFT)、快速傅里叶变换-二维多重信号分类(FFT-2D MUSIC)和三维多重信号分类(3D MUSIC)等算法。其中，3D FFT是一种基于快速傅里叶变换算法的三维参数估计方法。FFT-2D MUSIC方法利用FFT估计距离、2D MUSIC估计方位角和俯仰角，与3D FFT相比，提高了空间角分辨能力。3D MUSIC算法在此基础上将2D MUSIC扩展到3D MUSIC，提出一种距离-多普勒-方位角联合估计方法，它通过三维张量构造新的信号模型，并利用前向后向空间平滑技术对数据进行叠加和平滑，获得了三维参数的高分辨。但是，汽车和工业毫米波雷达系统往往方位角的视场(FOV)很大，距离分辨率高，导致划分的方位角网格数和距离单元数多。诸如上述结合MUSIC的点云成像方法，既需要已知源数量和大量的快照，又要对协方差矩阵的奇异值求解，以及联合估计的多层循环搜索，都极大增加了算法复杂度，算法运行的时间往往成倍增加。为了满足实时性的要求，TI公司提出了一种 FFT测距加Capon测角的三维点云成像方法，一方面利用Capon测角没有特征值分解带来的复杂运算，另一方面通过加入CFAR检测算法，检测存在目标的单元，来避免多维循环搜索，这种信号处理框架大大降低了方法运行的时间。

尽管如此，Capon算法的角分辨不足，因而为了改善到达角(DOA)估计的性能，文献[10]提出一种基于渐近最小方差的稀疏迭代方法(SAMV)。该方法在1D角度上划分网格，利用采样数据协方差矩阵与真实协方差矩阵的差异，根据AMV准则，推导出网格功率的迭代公式，从而得到更好的角度谱，有着比Capon和MUSIC更佳的分辨能力。本文在SAMV的基础上，将其扩展至2D，提出一种适用于毫米波FMCW MIMO雷达的三维点云成像方法。首先对差拍信号做FFT，得到1D距离像；然后在每个距离单元的水平方向上划分方位角网格，构造过完备的方向矢量矩阵，从AMV准则的角度出发，得出网格功率的迭代求解公式，由此得到距离-方位角热图；接着对热图的两个维度分别做1D CASO-CFAR，得到目标的距离和方位角单元；最后根据检测出来的单元位置信息，在俯仰角上划分网格，采用SAMV算法估计网格的信号功率，经过检测得到俯仰角。相较于现有方法，该方法在天线数量受限和有限快拍的情况下，能明显提高方位角和俯仰角上的分辨率。

1 信号模型

FMCW MIMO雷达系统模型如图1所示，有3个发射阵元～和4个接收阵元～，其中和的间距为2(为雷达信号的波长)，与和的间距均为，并向轴正方向平移2；到的阵元间距为2。为了方便后续的建模，按照文献[11]，将图1的MIMO阵列结构等效为图2的12个虚拟阵元收发共置来进行处理，对应的坐标位置矩阵=[,…,,…,]∈，其中=⊗+⊗，∈和∈分别是发射阵元与接收阵元的3D坐标矩阵。

图1 MIMO阵列模型

图2 等效虚拟阵列模型

MIMO阵列采用时分发射模式。以为例，假设在一个调频周期内发射的信号为

()=exp[j(2π+π)],0≤≤

(1)

式中，为信号时间，为每个Chirp的调频周期，为调频斜率，为中心频率。根据图2，接收的回波信号可以表示为

(2)

式中，为目标的反射系数，为发射天线到接收天线的信号传播时延，包括距离和角度时延，即

(3)

式中，为目标和原点的径向距离，=[sincos,coscos,sin]为目标的单位方向矢量。中的和为目标的方位角和俯仰角，取值范围为方位角和俯仰角的有效视场角和内，即∈(-2,2),∈(-2,2)。

将回波信号与发射信号进行混频，得到差拍信号为

exp[j2π(+++

(4)

式中，平方项和二次项相对于其他项可以忽略，又因的值太小，exp[j2π()]项也可以忽略，由此改写为

()=exp[j2π(·)]×

exp[j2π(2+2)]+

()

(5)

对()进行采样得到离散采样序列：

()

(6)

式中，为采样速率，=1,2,…,为采样点数。若将式(6)推广到个目标，设目标距离分别为,,…,，所在的方位角和俯仰角分别为(,),(,),…,(,)，则式(6)的信号模型可以进一步写为

()

(7)

()

(8)

2 算法描述

为了从式(8)得到3D高分辨点云像，需经过一维距离像估计、距离-方位角热图生成、CFAR目标检测和俯仰角估计等一系列的信号处理，具体流程如图3所示。

图3 三维点云成像方法流程框图

2.1 一维距离像估计

如图2所示，就第个虚拟阵元，先对第

=1,2,…,

(9)

2.2 距离-方位角热图生成

(10)

(11)

=1,…,,=1,…,

(12)

()=()+()

(13)

上述式(13)为用2D SAMV联合方位角和俯仰角估计时的信号模型，将其用于方位角网格功率估计时，需令=0，而且等于1和2，等于1到4。将()重新组合，得到

(14)

()=()+()

(15)

(16)

式中，()=diag((),…,(),…,())，对应的()为方位角上第个网格上的信号功率，()为噪声功率，()和()为待求的参数， 对协方差矩阵()矢量化，即

(())=(())=()

(17)

式中，()=[(),…,(),…,(),()]，为系数矩阵，即=[(,0)⊗(,0),…,(,0)⊗(,0),…,(,0)⊗(,0),()]，这里的⊗代表Kronecker积。对矩阵()的矢量表示可以看作未知的参数矢量()中元素的线性组合。由于Chirp数有限，可以得到采样协方差矩阵为

(18)

(19)

(20)

(21)

(22)

式中，为迭代次数。根据Capon波束形成算法对网格功率的估计，有=1/，则式(21)、式(22)可改写为

(23)

(24)

在上述迭代过程中，初始值可以通过周期图法来求解，即

(25)

(26)

=[(1),…,(),…,()]∈×

(27)

2.3 CASO-CFAR检测

将式(27)以距离维度和方位角维度分别送入CASO-CFAR检测器，其中=[(),…,(),…,()]，=[(1),…,(),…,()]，对每个方位角网格上的所有距离单元做1D CASO-CFAR，检测出存在目标的距离单元，再对距离单元处的所有方位角网格做1D CASO-CFAR，检测出存在目标的距离-方位角网格单元(,d)，得到对应的(,)。由于在应用中我们往往需要将距离较近的值能够检测出来，文献[15]中表明，单元平均选小恒虚警检测(CASO-CFAR)相对于单元平均恒虚警检测(CA-CFAR)，使得距离较近的目标检测出来而不出现目标遮蔽的现象。具体处理结构如图4所示。

图4 CASO-CFAR处理结构图

2.4 俯仰角估计

接下来，对检测的距离-方位角网格单元(,d)进行俯仰角估计时，根据式(13)，需令=，=d，且等于1到3，等于1到4，将()重新组合，得到

(28)

()=()+()

(29)

(30)

(31)

3 仿真与实验结果分析

3.1 雷达参数设计

为了对近距离区域内人员等目标的3D位置参数进行探测，需要进行合理的毫米波雷达参数设计。雷达参数主要包括MIMO雷达天线数量及阵列设计、调频起始频率、调频斜率、ADC采样点数、ADC采样率、工作波长和每帧的调频周期数，这些参数由应用场景的最大测量距离、距离分辨率等指标来确定，以下指标分别满足公式：

=×

(32)

(33)

(34)

(35)

雷达硬件支持的参数有上限，将其作为雷达参数设计的先决条件，包括发射功率决定的最大测量距离、最大的ADC采样率与最小的ADC采样率、最大的扫频带宽、最大的调频斜率与最小的调频斜率、水平及垂直方向最大视场和等。具体的参数设置如表1所示。

表1 雷达系统和目标参数

3.2 仿真结果分析

设置4个点目标，参数设置如表2所示。

表2 仿真目标三维参数

图5和图6分别为文献[9]方法和本文方法生成的距离-方位角热图以及经过CASO-CFAR检测的结果。从图中可以看出，文献[9]的方法已经无法在分辨出方位角上相差5°的目标1、2和目标3，只在目标1、2的方位角处有亮点，而方位角相差更大的目标4处也有明显的亮点。同时该方法的旁瓣较高，导致经过CFAR检测之后，真实距离-方位角单元的周围许多单元也被检测出来。相比之下，本文方法不仅在4个目标的真实方位角位置处均有明显的亮点，而且旁瓣较低，真实位置周围只有极少的单元被检测出来。图7和图8分别为两种方法在检测单元处估计的俯仰角网格功率谱图。从图中可以看出，文献[9]方法在所有检测单元处都只有一个谱峰，无法分辨出不同俯仰角的目标1和2，而本文方法的俯仰角功率谱在不同的检测单元上有两个不同位置的亮点。图9和图10分别为两种方法在空间直角坐标系中生成的3D点云像。从图中可以看出，文献[9]方法生成的点云像只能分辨出两个目标，而本文方法4个目标的点云都很清晰分辨出来，成像效果很好，只有目标4的俯仰角估计有5°的偏移。由此可见，相比较文献[9]的方法，本文提出的方法不仅在方位角上获得了更好的分辨率，而且在只有两排天线的条件下，俯仰角上依然可以得到一定的分辨率。

(a) 距离-方位角热图

(a) 距离-方位角热图

图7 文献[9]方法生成的俯仰网格功率谱图

图8 本文方法生成的俯仰网格功率谱图

图9 文献[9]方法生成的3D点云像

图10 本文方法生成的3D点云像

3.3 实测数据处理结果

实验所采用雷达系统由TI(德州仪器)公司的毫米波雷达IWR6843ISK评估板卡、毫米波传感器承载卡平台MMWAVEICBOOST和数据采集卡DCA1000组成。评估板卡可以设置雷达系统参数，天线系统如图1所示，将差频信号在设置的采样率下进行I/Q采样，之后将采样好的信号通过工业载板传输给DCA1000，DCA1000将处理好的差频采样复信号打包成二进制文件通过网口传输给电脑。

实验时的场景如图11所示，在距离雷达5 m处的两个人体目标，人体目标的高度大约为1.7 m，人体的宽度大约为0.5 m，厚度大约为0.2 m，两人相距为1 m，雷达的高度为1 m，图12和图13分别为两种方法生成的距离方位热图及CFAR检测结果。从图中可以看出，文献[9]中的方法，已无法分辨在方位角上相差5°的两个人体目标，而本文方法依然能够清晰地分辨出两个目标。图14和图15分别为两种方法在检测单元处估计出来的俯仰角网格功率谱图。从图中可以看出，文献[9]方法由于Capon测角算法的分辨率有限，且雷达在俯仰维上只有两排天线，无法有效地估计目标所在的俯仰角范围，而本文方法依然在目标的俯仰角范围有明显的亮斑。图16和图17分别为两种方法的三维点云成像结果。从图中可以看出，文献[9]方法只能分辨出一个目标的点云像，且因俯仰角上分辨率很低，只能取峰值位置，导致点云成像结果和目标轮廓相差很大。而本文方法能清晰地生成两个目标的点云，同时由于SAMV方法估计出来的是网格的真实功率，经过检测算法后，三维点云结果和目标的轮廓十分相似。

图11 实验场景

(a) 距离-方位角热图

(a) 距离-方位角热图

图14 文献[9]方法生成的俯仰网格功率谱图

图15 本文方法生成的俯仰网格功率谱图

图16 文献[9]方法生成的3D点云像

图17 本文方法生成的3D点云像

4 结束语

本文提出一种毫米波FMCW MIMO雷达高分辨三维点云成像方法，通过毫米波雷达的大带宽，得到高分辨的距离信息，引入SAMV方法，提升雷达的测角能力，结合CASO-CFAR算法，减少点云成像的时间。通过仿真和实验结果表明，相对于现有方法，该方法即使在天线数目有限的商业毫米波雷达上，依然有着更高的分辨率，兼顾了实时性的要求。本方法生成的近距离人体目标三维点云像，点云轮廓更加清晰、准确，便于后续基于点云的目标分类、姿势识别和计数等雷达应用的开发。