一种基于多通道比相定位的多人体目标定位跟踪方法与流程

本发明涉及雷达目标定位跟踪技术领域，特别涉及多个人体目标的定位跟踪方法。

背景技术：

在许多探测系统中，运动人体目标的检测和定位跟踪都有着重要的应用。例如，在智能家居领域可以实现房屋监测；在车载雷达领域，可以实现无人驾驶以及防撞等多种功能。

国内外许多研究机构开展了关于运动人体目标定位跟踪的研究。电子科技大学提出了一种基于非相干mimo雷达体制下的联合多目标检测和定位方法(y.ai,w.yi,m.r.morelandeandl.j.kong,“jointmulti-targetdetectionandlocalizationwithanoncoherentstatisticalmimoradar,”internationalconf.oninformationfusion,2014.)，它基于最大似然估计的理论，将多个目标联合最大化的问题分成多个非联合的优化问题，从而消除已经明确的目标所带来的干扰，实现对联合多目标的检测和定位。杭州电子科技大学提出了一种在图像域的目标定位跟踪方法(c.wangandl.zhang,“meanshiftbasedorientationandlocationtrackingoftargets,”2010sixthinternationalconferenceonnaturalcomputation(icnc),2010.)，该方法基于均值漂移算法，利用目标梯度角的概率密度分布为特征，通过多次迭代，实现了目标的方位估计和跟踪。

上述方法有复杂的算法流程和大的计算复杂度，而在实际的目标定位跟踪场合，需要一定的实时性。因此，研究一种具有实时性且定位准确，适用于室内室外场景的多个人体目标定位跟踪方法在雷达目标定位跟踪技术领域具有重要的价值。

技术实现要素：

本发明的发明目的在于：针对上述存在的问题，提供一种适用于多目标的定位方法，基于单发多收(simo)的线性调频连续波(lfmcw)雷达系统，首先利用动目标显示(mti)和去均值的方法对回波数据进行预处理；然后对回波数据的每一周期做二倍插值细化距离维加窗傅里叶变换，根据目标的距离谱，获取目标的距离信息；接着从频谱的目标所在峰值处提取出该目标对应的相位信息；最后，我们将不同接收通道的相位信息融合，从而计算出不同接收通道间的相位差，换算出目标的方位角或俯仰角，结合之前获取的目标距离信息，从而计算出目标在二维或者三维平面的定位结果；得到定位结果后，利用卡尔曼滤波的方法，实现对目标的跟踪。

本发明的一种基于多通道比相定位的多个人体目标定位跟踪方法，基于单发多收的线性调频连续波雷达系统执行以下步骤：

步骤1：回波矩阵预处理：

对雷达接收机各路接收通道的n×l维回波矩阵am分别在慢时间(按行)上进行去均值处理和动目标显示滤波，得到距离-脉冲域矩阵dm，其中n为每路接收通道一次定位时所处理的回波数，l为每个回波信号的采样点数m＝1,…,m，m为雷达接收机的接收通道数，即接收天线数目；

对各距离-脉冲域矩阵dm分别在距离维上(按行)先做加窗处理，再做二倍插值以及傅里叶变换，得到每个扫频周期的回波信号的压缩成sinc时域脉冲信号；

将n个扫频周期的sinc时域脉冲信号组合在一起，得到m路接收通道探测目标对应的距离-脉冲域矩阵zm；

将各距离-脉冲域矩阵zm分别在速度维上(按列)做加窗处理，以及傅里叶变换，得到距离-速度矩阵vm。

步骤2：目标检测：

分别将各距离-速度矩阵vm投影至距离平面，得到投影后的一维行向量pm，即pm的每个元素都是一个潜在目标的距离信息；

分别对m个向量pm进行单元平均恒虚警检测(ca-cfar)，得到检测目标(包含虚警目标)，分别从各向量pm中提取对应各检测目标的元素，得到各检测目标的距离信息其中上标p为初始目标区分符；

基于各距离信息对m路接收通道的检测目标进行融合处理：将距离信息相同的检测目标作为一个目标，并记录对应的距离信息；

对目标进行凝聚处理：两两对比目标的距离信息，若两目标间的距离间距小于人体目标的分布距离范围s，则将比对的两目标作为同一目标，保留幅值较大的目标，删除幅值较小的目标，直至所有目标间的距离间距均大于s；

用表示凝聚处理后的目标数，ri表示目标的距离信息，其中

根据距离信息ri的值，回溯到各距离-速度矩阵v1、v2、….vm中提取出各ri所在距离单元的峰值，得到各目标的速度信息vi；

由于人体目标的运动速度一般小于3m/s，雷达探测信号为lfmcw的时频关系可知，目标运动不会引起速度模糊，因此可以根据公式对目标进行距离解耦和，得到个目标的真实距离ri，其中f0,t,b分别为lfmcw信号的载波频率、信号时宽、信号带宽。

步骤3：目标定位及跟踪：

根据各个目标在距离-速度矩阵v1、v2、….vm对应的频谱峰值，在每个谱峰处提取与所述频谱峰值对应的相位值即第m接收通道回波数据处理后对应于第i个目标的相位值；

对同一目标的m路接收通道的相位值进行两两组合，计算每种组合下的初始角度信息其中的下标c用于区分不同的组合，上标i用于区分不同的检测目标，表示每种组合中同一目标的不同接收通道的相位值，λ表示载波波长，d表示接收通道l、k对应的接收天线的间隔。

优选的，在两两组合时，将接收天线相邻的两个接收通道的相位值组合在一起，再计算对应的初始角度信息

然后，再对同一目标的初始角度信息进行融合处理(例如加权平均)，得到各目标的角度值θi。

根据目标的真实距离ri、角度值θi计算出目标的定位坐标(xi,yi)：

获取同一目标在连续时间序列上的定位坐标(xi,yi)，得到目标的运动轨迹，完成对多人体目标的定位跟踪。

为了实现对目标的三维平面的定位及跟踪，对应的雷达系统的天线阵列为：发射天线位于坐标原点，各接收天线阵列均匀分布于坐标z轴和y轴上。在目标定位及跟踪处理时的步骤如下：

根据各个目标在距离-速度矩阵v1、v2、….vm对应的频谱峰值，在每个谱峰处提取与所述频谱峰值对应的相位值

将对应分布于y轴上的接收天线的接收通道作为第一接收通道，通道数目用m1表示，将对应分布于z轴上的接收天线的接收通道作为第二接收通道，通道数目用m2表示，且m1+m2＝m；则将对应第一接收通道的相位值表示为：对应第二接收通道的相位值表示为：其中m1＝1,…,m1，m2＝1,…,m2；

对同一目标的所有进行两两组合，计算每种组合下的初始俯仰角其中的下标c用于区分不同的组合，上标i用于区分不同的检测目标，表示每种组合中同一目标的不同第一接收通道的相位值，λ表示载波波长，d表示第一接收通道l1、k1对应的接收天线的间隔；并对同一目标的初始俯仰角进行融合处理，得到各目标的俯仰角

对同一目标的所有进行两两组合，计算每种组合下的初始轴偏角其中的下标c用于区分不同的组合，上标i用于区分不同的检测目标，表示每种组合中同一目标的不同第一接收通道的相位值，λ表示载波波长，d表示第一接收通道l2、k2对应的接收天线的间隔；并对同一目标的初始轴偏角进行融合处理，得到各目标的轴偏角

根据目标的真实距离ri、俯仰角轴偏角计算出目标在三维平面的定位坐标(xi,yi,zi)：

获取同一目标在连续时间序列上的定位坐标(xi,yi,zi)，得到目标的运动轨迹。

在分别对进行两两组合时，仅将相邻接收天线的对应通道的进行组合，以便于进一步减少计算量。

为了使目标的运动轨迹更为平滑，以及去除一些异常值，得到目标的运动轨迹后，可以采用卡尔曼滤波的方法对运动轨迹进行滤波处理。

综上所述，由于采用了上述技术方案，本发明的有益效果是：本发明提出了一种单发多收(simo)的雷达系统下，适用于室内和室外环境中多个人体运动目标的检测、定位以及跟踪方法，该方法能够对多目标进行有效的检测定位和跟踪，在提高多目标定位精度的同时，不会造成检测目标的丢失。无论在室内还是室外的环境中，本发明均能够呈现出更好的目标定位效果，同时，结合卡尔曼滤波，本发明能够对目标实现实时的跟踪。本发明能够保证雷达探测系统对多目标定位跟踪的实用效果，在车载雷达等方面，有广泛的应用前景。

附图说明

图1为具体实施方式的处理流程图。

图2为定位场景示意图。

图3为多目标人体走动轨迹示意图

图4为多人体目标定位跟踪结果图。

图5为三维平面目标定位仿真场景示意图。

图6为三维目标定位仿真结果。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面结合实施方式和附图，对本发明作进一步地详细描述。

将本发明用于对多个运动目标的定位和跟踪，其室内实测数据场景示意图如图2所示，坐标系零点位于发射天线位置，4个接收天线沿x轴等间隔放置，相邻两个接收天线之间的水平间隔为半波长d＝λ/2，雷达的发射信号为中心频率24ghz、带宽500mhz的线性调频信号，信号时宽为300us。

论证实验中，三个人体目标分别从不同的起点开始出发，实际的运动轨迹示意图如图3所示。

参见图1，对上述场景采用本发明的多人体目标定位跟踪方法的处理步骤为：

s1：通过一发多收的lfmcw雷达系统回波信号(回波数据)，本实施例中，采用的是1发4收的lfmcw雷达系统，即雷达接收机的接收通道数m＝4，每个接收通道(接收天线)一次定位时所处理的回波数n＝128，每个回波信号的采样点数l＝128；

s2：对4路接收通道所得的n×l维回波矩阵a1、a2、…、a4分别在慢时间上(按行)进行去均值处理和mti滤波，得到去除了静止背景等零频杂波的距离-脉冲域矩阵d1、d2、…d4；

s3：对d1、d2、…d4分别在距离维上(按行)加窗处理且做二倍插值fft(傅里叶变换)，使每个周期的回波信号压缩成sinc时域脉冲信号，sinc带有了目标的距离信息；再将n个扫频周期上进行的上述操作得到的n个sinc信号组合在一起，得到4路接收通道探测目标对应的距离-脉冲域矩阵z1、z2、…z4；

s4：对z1、z2、…z4分别在速度维(按列)做加窗fft，得到对应的距离-速度矩阵v1、v2、….v4；

即，通过上述步骤s1～s4完成对回波矩阵的预处理。

s5：将步骤s3得到的距离-速度矩阵v1、v2、….v4投影至距离平面，得到投影后的一维行向量p1、p2、….p4；

s6：对p1、p2、….p4进行单元平均-恒虚警检测，获取目标(包含虚警目标)及其距离信息；

s7：对目标进行凝聚处理后，凝聚处理后的目标及其距离信息(r1、r2、r3)，即真实目标的距离信息；

s8：根据r1、r2、r3的值，回溯到v1、v2、….v4中提取出r1、r2、r3所在距离单元的峰值，从而获取目标的速度v1、v2、v3；

再根据(其中i＝1,2,3，f0,t,b为lfmcw信号的载波频率、信号时宽和信号带宽)对目标进行距离解耦和，从而得到探测目标对应的真实距离r1、r2、r3；

根据各个目标在距离-速度矩阵v1、v2、….v4对应的频谱峰值，在每个谱峰处提取出该峰值处数值对应的相位值

基于相邻接收天线所对应的接收通道的相位值计算初始角度信息在对同一目标的初始角度信息求平均，可得到目标的角度值θi；

即，通过上述步骤s5～s8完成目标检测处理。

s9：根据步骤s8得到的真实距离r1、r2、r3，可以计算出目标的定位坐标：

s10：重复步骤s2～s9，可以得到各目标的一系列定位坐标，从而得到目标粗略的目标运动轨迹，为了使目标运动轨迹更为平滑，以及去除一些异常值，采用卡尔曼滤波的方法粗略的目标运动轨迹进行滤波，得到最终的目标运动轨迹，如图4所示。图中目标探测到的运动轨迹与实际的目标运动轨迹基本一致。

上述在二维平面的目标定位和跟踪处理过程可以扩展至三维平面，其三维平面目标定位仿真场景示意图如图5所示，其中发射天线t位于坐标原点o，接收天线r阵列均匀分布于坐标z轴和y轴上，见图5中所示的r1,r3,r4,rm-3,rm-2,rm-1,rm，相邻天线的间隔d＝λ/2。基于二维平面的定位原理，沿y轴分布的接收天线阵列可以探测出三维平面中目标pt(坐标为(xt,yt,zt))的俯仰角沿z轴分布的接收天线阵列可以探测出三维平面中目标pt的轴偏角图中的pt′表示目标pt的水平面投影。结合探测到的目标的真实距离r，可以计算出目标在三维平面的定位坐标：

三维目标定位的理论仿真结果如图6所示。

由上述实测数据结果和仿真结果可知，本发明提供的适用于多个人体目标检测、定位和跟踪的方法不仅能够有效地检测到目标的个数，还能对目标实现准确的定位，其定位误差小于0.5米，更能对目标的运动轨迹进行实时的跟踪，验证了本发明的有效性和实用性。