一种车载毫米波FMCW MIMO雷达快速成像方法

黄以兰， 晋良念，2， 刘庆华

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

0 引言

车载毫米波雷达作为汽车主动安全领域的关键传感器部件，可有效穿透雾、烟、灰尘，实现全天时、全天候工作负荷要求。这种雷达采用的体制现有合成孔径雷达(SAR)和多发多收(MIMO)雷达等。相比于SAR，MIMO雷达是利用发射和接收天线间的位置关系形成更大孔径的虚拟天线阵列。考虑到车载雷达的实现复杂度及半导体器件成本限制，本文采用时分复用MIMO(TDM-MIMO)雷达工作机制，该方式对硬件性能需求较低、实现简单，减少MIMO雷达的数量和复杂性，降低了成本。

传统的成像方法包括后向投影(Back Projection，BP)、延时求和(Delay-And-Sum，DAS)和数字波束形成(DBF)等。文献[4]基于SAR成像几何模型，结合车载平台推导了BP算法，通过仿真实验验证了其在车载前视阵列雷达的适用性。文献[5]提出了在距离向进行距离压缩与方位向进行DAS，最后结合距离-方位信息得到二维成像结果。文献[6]提出的DBF通过加权求和使多路信号变为一路信号，使得目标方向的信号得到增强，同时对非目标方向的信号进行抑制。

尽管这些算法能够快速而简单地实现目标的成像，但是都存在分辨率低且旁瓣较高的问题。为此，采用高分辨成像方法是必要的。文献[7]和[8]分别介绍了在获得距离向信息后采用ESPRIT、MUSIC超分辨算法，虽然相比于传统算法有着更好的角度分辨率，但这些算法需要对接收信号的协方差矩阵作特征分解，需要的数据帧数较多，对均匀线阵中天线的数量也有着较高的要求。文献[9]提出了IAA在FMCW MIMO雷达系统中的应用，该方法在帧数较少的情况下可使角分辨率从DAS波束形成的10°提高到5°，但存在算法复杂度较高，无法实现工程实时性的问题。

针对上述问题，本文提出了一种车载毫米波FMCW MIMO雷达IAA高分辨成像的快速实现方法，该方法首先基于MIMO机制采用FFT来获取目标的距离信息，然后对存在目标的距离单元内利用傅里叶算子特性和GS因子分解进行数据协方差矩阵及其逆矩阵的快速计算，利用快速Toeplitz矩阵向量乘法有效计算目标反射系数，迭代次数少(10次左右)、计算量适中。

1 信号模型

雷达系统模型如图1所示，采用时分复用FMCW MIMO雷达(时分复用带来的相位误差补偿方法具体见文献[10])，根据MIMO雷达原理，以第一个接收阵元为参考点，有个发射阵元，各阵元距离参考点的间距为t,=1,2,…,；接收阵元有个，各阵元距离参考点的间距为r,=2,3,…,-1。假设第个阵元发射信号为

()=exp(j2π+jπ)

(1)

式中，为信号时间，为载波频率，=为调频斜率，为信号带宽，为调频周期。

图1 系统模型

将成像区域划分为个距离单元,,…,，个角度单元,,…,，则对于第个接收阵元来说接收到该区域内所有目标的回波信号可以表示为

(2)

式中:(,)表示加权指示函数，若在(,)位置上有目标，则(,)≠0，反之(,)=0；(,)表示目标与第个发射阵元和第个接收阵元之间的传播时延，其为

(3)

)+j2π(,)]

(4)

对中频信号进行采样得到离散采样序列：

)+j2π(,)]

(5)

式中，为采样间隔，=1,2,…,为采样点数。对()计算点FFT并取单边谱，得到该通道的一维距离像为

exp[j2π(,)],=1,2,…,

(6)

从式(6)可以看出，谱峰位于处，即=(,)。将式(3)代入式(6)并化简，则

(7)

(8)

式中，()为对应第个距离单元角度目标的散射系数。

最后，将所有发射阵元和接收阵元组合的回波数据堆叠，并将距离的产生的双程时延补偿，得到第个距离单元的测量向量，即

(9)

=+

(10)

2 算法描述

根据文献[11]所述，IAA通过最小化下面的加权最小二乘(WLS)代价函数来求解式(10)，即

(-()())(-()())

(11)

式中，为干扰数据协方差矩阵，且

(12)

这里的为数据协方差矩阵，即

(13)

将式(11)最小化并代入式(12)，再根据矩阵求逆引理，得到目标散射系数()为

(14)

根据式(13)可知，依赖于()的值。首先，使用DAS估计值初始化IAA，即

(15)

如式(13)、式(14)所示， IAA的计算负担主要来自每一次迭代的矩阵及其逆矩阵的计算和()中分子、分母的计算，因此需要的存储空间大，计算复杂度也高。接下来，根据()的傅里叶矩阵特性快速计算，然后根据快速Toeplitz矩阵向量乘法分别计算()中的分子、分母，以减少每次迭代的计算量。令

(16)

(17)

2.1 R的快速计算

根据公式(13)，令()=|()|，=diag([()() …()])，则可以表示为

(18)

将代入式(18)，可得

(19)

(20)

由式(20)可以看出，是Toeplitz矩阵，令

1,…,-1

(21)

则式(19)的可以简化为

(22)

2.2 φD(θk)的快速计算

(23)

定义线性预测和相关的预测误差如下：

(24)

(25)

由式(24)给出的线性方程组可以通过Levinson-Durbin算法。根据文献[13]，的GS因子分解表达式可表示为

=(,)(,)-(,)(,)

(26)

式中，是下三角Toeplitz矩阵，其中的，和×维移位矩阵分别为

(27)

(28)

(29)

2.3 φN(θk)的快速计算

将式(26)的代入式(17)，可得

(,)(,)]()=

(30)

(31)

定义=(-+1,-+2,…,)，则其可以表示为

(32)

由式(32)可知,c是两个下三角Toeplitz矩阵与向量乘积的相减。因此,c也可以通过快速Toeplitz矩阵向量乘法来快速计算。令ΦN=[φN(θ1),…,φN(θk),…,φN(θK)]T,根据式(31),ΦN也可以类似于φDθk 那样,通过对c进行FFT运算来实现其快速计算。2.4 算法流程及复杂度分析 算法流程如表1所示。首先采用FFT计算估计初始值,迭代计算式(13)、式(14)的过程中通过FFT运算快速计算R,根据GS因子分解求解R-1,采用Toeplitz矩阵向量乘法的快速算法及am(θk)的傅里叶算子特性快速计算φD(θk)、φN(θk),直至算法收敛,结束迭代。

表1 FIAA算法流程

从前面的分析可以知道，文献[9]提出的IAA每次迭代过程都需要对()和进行更新，总共需要的计算复杂度为(2++)，其中对()的求解过程需要进行，其计算复杂度为()，而本文提出的FIAA在计算的过程中的复杂度仅为()。此外，通过Levinson-Durbin算法求解式(24)需要计算()次，根据GS因子分解求解需要5(2)次。求解()、()过程中两次使用的快速Toeplitz矩阵向量乘法需要7(2)次，加上FFT运算需要的3()次，FIAA共需要计算复杂度为(+12(2)+3())。所以，整体看来本文给出的FIAA处理计算复杂度都有较大程度的降低。

3 仿真与实验结果分析

3.1 仿真结果分析

为保证距离和角度的解算范围以及对多目标的分辨能力，要对车载毫米波雷达系统的工作参数做出限定和选取。系统的距离及角度分辨率分别为

(33)

为了评估雷达性能并比较DAS波束形成、文献[9]提出的IAA和本文提出的FIAA算法，在MATLAB中对场景目标进行仿真，参数如表2所示。

表2 系统参数

假设目标场景的距离单元数为=256，角度范围是-60°～60°，角度单元间隔为1°，即角度单元数=121。设置3个点目标，参数如表3所示，图2给出了不同算法功率谱对比，其中图2(a)～(c)分别为DAS波束形成、IAA、FIAA算法的功率谱图。目标的成像结果如图3所示，其中图3(a)～(c)分别为DAS波束形成、IAA、FIAA算法的仿真数据的距离-角度二维成像图。从这些图中可以看出，DAS波束形成算法只能分辨出目标1和目标3，分辨不出相距5°的两个目标，由此可见分辨率大约为10°，而且旁瓣较高。而IAA与FIAA算法可以很清晰地分辨出3个目标，其分辨率可达到5°，旁瓣较低。因此，相比于DAS波束形成算法具有很高的旁瓣和很低的分辨率，很难分辨具体的方位角， IAA与FIAA算法都获得了较好的角度估计以及有效地降低了旁瓣。

表3 目标参数

(a) DAS波束形成

(a) DAS波束形成

下面对算法的运行时间进行对比，以单个目标为例，设置单目标距离=12 m，角度=10°；以多个目标为例则设置多目标参数如表3所示。表4给出了DAS波束形成、文献[9]提出的IAA与本文提出的FIAA分别在单目标场景及多目标场景中运行时间的对比结果。从表中可以看出，无论是在单目标还是多目标场景，DAS波束形成算法所需时间最短，实时性最高；IAA所需时间最长，不符合车载雷达实时性的需求；而FIAA的运行速度虽比不上DAS波束形成，但相比于IAA，在单目标场景中FIAA的运算时间减少7倍左右，在多目标场景中运算时间减少了4倍左右。由此可以说明所提算法性能较好，有效地减少了运行时间，更适合应用于要求实时性的车载场景中。

表4 DAS波束形成、IAA、FIAA的对比

综上所述，DAS波束形成运行时间短，但存在分辨率低和高旁瓣的问题，文献[9]提出的IAA虽然提高了分辨率但算法复杂度高且运算时间较长，本文提出的FIAA在保证有效提高分辨率、降低旁瓣的基础上，能够减少算法的运行时间，提高性能。

3.2 实测数据处理结果分析

实际探测场景如图4所示，雷达系统主要由 TI(德州仪器)公司的高性能毫米波雷达前端IWR6843 评估板卡和DCA1000数据采集卡组成。数据采集卡连续接收雷达前端输出的数字差频信号并通过 USB 接口将原始回波差频数据发送给计算机存储，为FMCW MIMO 雷达信号处理提供实测数据。

图4 实验场景

实验过程中，首先对暗室无目标场景进行测量，为后续数据预处理作准备。然后测量的是相同距离3 m处不同角度(=0°，=20°)的双目标，经过去除耦合波处理后结果如图5所示，其中图5(a)～(c)分别为DAS波束形成、IAA、FIAA算法的实测数据的距离-角度二维成像图。从图中可以看出，相比于DAS波束形成算法，FIAA得到的方位图像具有更清晰的亮点以及更好的角度分辨率。在运行时间上，采用DAS波束形成算法、文献[9]提出的IAA、本文提出的FIAA分别为0.002，0.106和0.022 s，FIAA比IAA的运行时间缩短了5倍左右，可以发现本文方法在提高分辨率的基础上有效地降低了运行时间，满足了车载场景对实时性的需求。

(a) DAS波束形成

4 结束语

本文提出的车载毫米波FMCW MIMO雷达快速成像方法，解决了DAS波束形成算法测角分辨率较低、旁瓣多和子空间类算法需要的数据帧数较多以及运算复杂度较高的问题。仿真数据和实验数据的处理结果表明：相比于波束形成与其他高分辨率算法，FIAA得到的方位图像具有更清晰的亮点以及更好的角度分辨率。在今后的工作中，继续探索车载毫米波雷达高分辨测角方法的快速算法，并将其扩展到 3D成像中，致力于将该方法更完善地应用于更复杂真实的环境中。