基于置信规则推理方法的雷达目标跟踪

刘威, 陈先桥，初秀民

(1. 武汉理工大学 计算机科学与技术学院，湖北 武汉 430063； 2. 武汉理工大学 智能交通系统研究中心，湖北 武汉 430063； 3.水路公路交通安全控制与装备教育部工程研究中心，湖北 武汉 430063)



基于置信规则推理方法的雷达目标跟踪

刘威1,2,3, 陈先桥1,2,3，初秀民2,3

(1. 武汉理工大学 计算机科学与技术学院，湖北 武汉 430063； 2. 武汉理工大学 智能交通系统研究中心，湖北 武汉 430063； 3.水路公路交通安全控制与装备教育部工程研究中心，湖北 武汉 430063)

摘要：针对雷达目标跟踪中数据相关联的问题，提出一种基于雷达目标航向和速度对雷达数据进行数据关联的方法，解决目标密度较大时所出现的目标误跟或丢失问题，并采用基于置信规则库(belief rule-based，BRB)的方法建立了目标速度与航向的置信规则库，通过证据推理(evidential reasoning，ER)算法解决了不同量纲数据之间的冲突问题，从而验证了此方法的可靠性。实验结果表明，其真实目标置信度可达到95%以上。最后，运用此方法将得到的真实目标点进行Kalman滤波，实现对雷达目标的跟踪，通过Monte Carlo仿真验证了此方法具有很好的目标跟踪效果。

关键词：雷达目标跟踪；置信规则库；证据推理；航向；航速

网络出版地址：http://www.cnki.net/kcms/detail/23.1390.u.20160421.1040.028.html

随着内河水运信息化的飞速发展，为了保证航运业安全有效地运行，雷达被广泛运用于追踪船舶运行状态方面，从而避免船舶发生碰撞造成不必要的经济损失。特别是长江干线，现在已经成为了世界上运输量最大并且最繁忙的通航河流。因此为了防止航行事故的发生，对船舶进行交通管制显得尤为重要。普通船用雷达一般通过人工标绘来获取船舶的航速、航向等航行数据。随着技术的不断发展，现在的船用雷达普遍具有自动标绘功能(automatic radar plotting aid，ARPA)，可以自动获取上述数据，对运动的目标进行提取与跟踪。然而，近年来，随着跨临河建筑物的快速增加、船舶交通流量的不断增加，雷达所采集到的ARPA目标数据含有大量噪声，给目标跟踪带来了很大干扰。如何去除这些干扰，以达到更好的跟踪效果，成为了一个亟需解决的问题。

目标跟踪的含义是利用远程测量数据，通过使用数据关联算法和Kalman滤波器来估计运动目标的状态[1]。根据跟踪算法所采用的结构，可以采用不同类型关联程序[2-5]，一旦数据关联完成，线性的或扩展的卡尔曼滤波器便对目标的位置和速度给出进一步的估计值。雷达数据处理过程中，数据关联是目标跟踪的核心部分。文献[2]通过采用Kalman滤波器与最近邻算法相结合实现了对车辆轨迹的追踪。文献[6]通过将聚类分析方法运用在概率数据关联算法中，以降低目标在复杂环境下的关联错误率，提高了目标量测的稳定性。最近邻算法[7]原理简单，计算量小；概率关联算法[8]与标准卡尔曼滤波在最大存储量上几乎相等，甚至基本不变，因此其算法容易实现。但在目标密度较大时，这两种方法都容易出现丢失目标或误跟目标的问题。

本文根据获取的APRP目标，提出一种基于目标的速度和航向对目标数据进行关联，去除真实目标中的噪声干扰，并将速度和航向作为BRB(Belief rule-based)系统的前置属性，建立置信规则库，通过ER算法来演绎推理，以此验证该方法的可靠性，最后采用Kalman滤波实现对目标的跟踪算法。

1基于速度和航向的目标航迹检测

本文所采集的目标数据均为ARPA雷达返回的雷达数据。虽然ARPA设备具备一定的自动识别目标能力，但是在内河水域中它会将礁石、桥梁等目标都识别为船只。为此，在滤波之前，必须将这些错误的目标点剔除。本文提出将选取的目标点先进行速度和航向的检测，如果满足条件，则将此目标点作为录取点来进行Kalman滤波，从而形成一个航迹上的跟踪点。

1) 速度检测。

常见的船舶速度值为0～20 m/s，其转弯速度或加速度为2 (°)/s。因此，在选取目标点时，首先对其速度进行检测，如果目标点的速度太小或超过20 m/s，则认为此目标点为伪目标点，并将其舍去，只有速度在最大速度和最小速度之间的目标点，才给予录取。但是仅仅依靠速度上的检测是不够的。如图1所示，满足速度检测条件的点迹有A、B。

图1　速度检测Fig.1　Speed detection

根据最近邻域准则A点离预测目标点最近，因此A点应该是所要关联的点迹。但因为船舶不可能在很短的时间内做出接近90°的转弯到达A点，因而A点不是真实目标点迹，而相对较远的B点才是所要选取的目标点。因此还需要对目标点迹做出符合航向上实际转向的关联，也就是航向检测。

2) 航向检测。

当观测点符合速度检测后，再对其进行航向检测，如果观测点也符合航向检测，最后将符合两者条件的点作为最终可以相关联的目标点，进行目标点迹关联。航向检测图如图2所示。

图2　航向检测Fig.2　Course detection

图2中，A、B、C为真实目标点迹，Z1、Z2为目标观测点迹；CF为航迹BC的延长线，以CF为中线，两边分别有角度为θ的航向限制角，即CD、CE与CF的夹角。设目标观测点与CF的夹角为α，当α在θ范围内时，则满足航向检测，否则不满足。图2中可以看到，观测点Z1与CF的夹角在θ范围内，故满足航向检测条件，可以进行目标点关联；而观测点Z2与CF的夹角大于θ，故不满足航向检测条件，不予关联，将其舍弃。

对于目标航迹关联的整体流程图如图3所示。

图3　目标航迹关联流程图Fig.3　Target track association flowchart

2BRB推理验证航迹检测法

表1　速度与航向样本表

在BRB系统中，其推理过程包括：多种形式信息的输入转换、规则激活权重的计算以及使用ER聚集BRB规则得到最终置信度。

2.1BRB输入转换

一个前提属性值的输入转换指的是将此值转换成不同的置信度，并将这些置信度分配给此前提属性的不同参考值。这就相当于将一个输入值转换成一个对应于前提属性参考值的置信度分布。特别地，一个前提属性Pi的输入值(它的置信度假设为xi)可以通过以下置信分布来转

High=(1-Medium),Low=0.00

(2)

Medium=(1-Low),High=0.00

(3)式中：H代表最大值，xi代表输入值，M代表中间值，L代表最小值。

在进行输入转换之前，需要设定前提属性的参考值，假设两个前提属性的参考值分布如下：速度={(Fast，10.4)，(Normal，5.3)，(Low，0.4)}；航向={(Large，10.8)，(Normal，4.7)，(Small，1.2)}。然后结合式(8)、(9)完成对置信度αin的计算。比如，选择一个样本，速度=6 m/s，航向=4°，则速度隶属于{(Fast，10.4)，(Normal，5.3)，(Low，0.4)}的程度为(0.1473，0.8527，0)，即Medium=(10.4-6)/(10.4-5.3)=0.862，High=1-0.8627，Low=0；同理，航向隶属于(Large，10.8)，(Normal，4.7)，(Small，1.2)的程度为(0，0.800，0.200)。

2.2计算激活权重

一般地，在置信规则中，连接符号“∧”通常被用来代表前置属性的逻辑关系。这就意味着只有当规则中所有前因被激活时，得到的结果才具有可信度。根据上述置信分布的计算，第k条规则下激活权重ωk可以由式(4)来计

2.3ER推理输出

ER方法是用来集合所有L规则下的前提属性数据分组，从而通过给定的前提属性Pi输入值获得结果属性中每一个参考值的信任程度。本文采用分析的ER算法[11]，其产生的输出结果O(Y)由结果属性的参考值组成，如式(5)所示。

(5)式中：βj代表的是对应于一个结果参考值Cj的置信度。βj的计算由ER算法的分析方程式(6)来求得。

最终的合并结果或通过ER推理所得到的输出表达式为：{(C1,β1),(C2,β2),…,(CN,βN)}，其中βj是隶属于结果属性中第j个参考值Cj的最终置信度。本文中选取的输出结果分布分别是真实目标、虚假目标，根据表1的样本数据，结合式(1)～(6)，其训练后最终的置信规则库如表2所示。

(6)

表2中的置信度规则反映了速度、航向与真实目标之间的因果关系，以规则5为例，其含义是：在当前时刻，速度的取值在正常范围内，而航向也在其规定的正常范围内时，其确定的真实目标的置信度为93.49%，虚假目标的置信度为6.51%，且该规则成立的置信度为1，也就是说完全相信此规则下所得到的结果是成立的。其他的规则可以用同样的方法来解释。因此，文中所提出的以速度和航向来检测目标航迹的方法具有较高可靠性。

3目标航迹Kalman滤波

Kalman滤波器[14]输入的是笛卡尔坐标系下的目标位置，在目标跟踪过程中，通过时间tn-1下的系统状态，结合时间tn下的测量参数，利用Kalman滤波进行当前时刻tn下系统状态估计。该滤波器以递归的方式来计算它所得到的解[5]。

当前时间标记下的系统状态包含的是笛卡尔坐标系下目标点的位置和速度[1]：

(7)

假设目标是匀速运动，则基于先前状态下的时刻t的状态预测方程可以由(8)式来表述：

(8)

式(8)等价于：

(9)

式中：Wtn-1是期望值为零和协方差矩阵为Ex的高斯噪声状态向量，A是状态转移矩阵，T代表的是测量周期。

本文中，对实验的结果只考虑目标位置的估计，因此测量矩阵只包含了目标位置的信息：

(10)

因此观测方程可以表示为

(11)

式中n为高斯零均值和协方差矩阵Ey的测量噪声。

Kalman方程描述当前的状态和状态估计更新的预测为

(12)

协方差为测量噪声Ey和状态的误差估计Ex具有以下含义：当测量值受噪声的影响较大时，可以给测量噪声Ey赋予较大的权值。这样处理即默认更多地依赖于基于模型的滤波器的预测值，而不是噪声Ey的预测值，同时也假设所依靠的模型有较高的精确度，因此协方差矩阵具有更小的值。

4仿真结果及分析

使用真实的测量数据在MATLAB环境下进行目标点航迹关联以及Kalman滤波跟踪的仿真。航迹起始位置为(150m,100m)，初始速度为(40m/s,30m/s)，采样间隔T=2s，MonteCarlo仿真次数为150次。仿真流程图如图4所示。

图4　目标航迹仿真流程图Fig.4　Target tracks simulation flowchart

图4中p代表观测点的位置，vp表示径向速度，θ表示航向。在每一次时间戳t，对上述这些数据进行检测关联，速度大于20m/s或小于0的点舍去；不在航向θ以内的点舍去，在这里选取的航向范围为：0°<θ<60°；最后结合最小邻域法选取最合适的目标点，并将其送入Kalman滤波器中进行目标跟踪，最终显示目标的航迹。通过获取的ARPA雷达数据得到船舶的位置信息，其位置绘制如图5所示。

图5　目标的真实位置Fig.5　Real position of targets

图5中有900个ARPA目标点，可以明显看到目标的轨迹位于笛卡尔坐标的平面内，但是在这些目标点中，并不全是真实的目标点，还有些点属于杂波干扰所导致的虚假目标点，这些目标点需要被剔除。接下来的步骤是数据关联，在一定的范围内，不符合速度和航向判断条件的目标被去除掉，只有有效的目标点被显示出来。其数据关联后的目标点迹如图6所示。

图6　数据关联后的目标点迹Fig.6　Target trajectory after data association

图6中可以看到在进行数据的关联后，之前的干扰点以及不符合判断条件的目标点都被很好地剔除掉，为后续的目标跟踪及航迹显示提供了更加真实可靠的目标点。在进行Kalman滤波跟踪后所得到的目标轨迹平滑度取决于测量以及状态噪声协方差矩阵的选择。通过实验可知，要想得到平滑度较好的轨迹，需要对测量噪声矩阵分配较大的权值。Kalman滤波后的目标轨迹如图7所示。

图7　Kalman滤波跟踪后的目标轨迹Fig.7　Target trajectory after Kalman filtering

从以上实验结果看出，利用一个简单的线性Kalman滤波器，结合基于目标点径向速度、航向的检测(结合最近邻关联算法)的算法简单可行，且能得到令人满意的跟踪效果。图8给出三种算法对目标跟踪的RMSE均方根误差，其中(a)为概率数据互联算法(PDAF)，它主要用于解决杂波环境下的目标跟踪问题，并且误跟和丢失目标的概率较小，计算量相对较小。可以看到本文所提出的算法具有更高的跟踪精度。此外，表3给出了150次MonteCarlo仿真实验中，三种算法所用的耗时，进一步论证相对于其他两种算法，本文算法的计算复杂度更小，跟踪实时性高。

表3　目标跟踪中各算法的用时

图8　三种算法的RMSE比较Fig.8　Three algorithms RMSE for target tracking

5结论

1)建立了基于BRB的雷达目标识别系统，该系统与传统的IF-THEN规则库相比较，其在知识的描述上更加丰富，应用上更具有现实意义。结果表明通过速度和航向来检测目标航迹具有较高的置信度。

2)为了实现真实目标的轨迹跟踪，构建了两个重要的信号处理结构单元：数据关联单元和Kalman滤波跟踪单元。数据关联阶段主要完成更新被测目标现有的轨迹，通过对目标点径向速度和航向的判断，来实现虚假目标点的删除，其目的是尽可能地减少错误目标轨迹的生成。Kalman滤波器在本文中被用来消除测量噪声，使得真实目标轨迹更加平滑。

3)将本文算法与其他传统算法进行了对比分析，实验结果表明，所提算法不仅在跟踪精度和鲁棒性上更好，而且用时也相对更少。

总之，通过仿真实验分析，本文所提算法无论是在雷达目标跟踪精度上，还是在可操作性上，均符合实际工程的要求，可以为雷达目标的跟踪与实现提供一种新的途径。

参考文献：

[1]BAR-SHALOM Y, LI X R, KIRUBARAJAN T. Estimation with applications to tracking and navigation: theory, algorithms and software[M]. New York: John Wiley & Sons, 2001.

[2]MACAVEIU A, CMPEANU A. Automotive radar target tracking by Kalman filtering[C]//Proceedings of the 2013 11th International Conference on Telecommunication in Modern Satellite, Cable and Broadcasting Services. Nis, Serbia: IEEE, 2013, 2: 553-556.

[3]陈志敏, 薄煜明, 吴盘龙, 等. 拟蒙特卡罗粒子滤波改进算法及其在雷达目标跟踪中的应用[J]. 应用科学学报, 2012, 30(6): 607-612.

CHEN Zhimin, BO Yuming, WU Panlong, et al. Improved quasi-monte-carlo particle filtering and its application to radar target tracking[J]. Journal of applied sciences, 2012, 30(6): 607-612.

[4]KUHN H W. The Hungarian method for the assignment problem[J]. Naval research logistics, 2005, 52(1): 7-21.

[5]SIMON D. Kalman filtering with state constraints: a survey of linear and nonlinear algorithms[J]. IET control theory and application, 2010, 4(8): 1303-1318.

[6]李秀良. 基于聚类分析的雷达航迹跟踪数据关联算法[J]. 指挥信息系统与技术, 2010, 1(3): 62-65.

LI Xiuliang. Data association algorithm in target tracking process of radar based on clustering analysis[J]. Command information system & technology, 2010, 1(3): 62-65.

[7]ABBAS M A, Shoukry A A. CMUNE: a clustering using mutual nearest neighbors algorithm[C]//Proceedings of the 2012 11th international conference on information science, signal processing and their applications (ISSPA). Montreal, QC: IEEE, 2012: 1192-1197.

[8]叶晓雪, 王新民, 李俨, 等. 修正概率关联算法在航迹关联中的应用[J]. 计算机仿真, 2011, 28(3): 19-21, 30.

YE Xiaoxue, WANG Xinmin, LI Yan, et al. Application of modified PDA algorithm in dense clutter[J]. Computer simulation, 2011, 28(3): 19-21, 30.

[9]李彬, 王红卫, 杨剑波, 等. 基于置信规则推理的库存控制方法[J]. 华中科技大学学报:自然科学版, 2011(7): 76-79.

LI Bin, WANG Hongwei, YANG Jianbo, et al. Belief rule-based inference method for inventory control[J]. Journal of Huazhong University of Science and Technology :natural science edition, 2011(7): 76-79.

[10]YANG Jianbo, LIU Jun, WANG Jin, et al. Belief rule-base inference methodology using the evidential reasoning approach-RIMER[J]. IEEE transactions on systems, man, and cybernetics-part a: systems and humans, 2006, 36(2): 266-285.

[11]YANG Jianbo, XU Dongling. Evidential reasoning rule for evidence combination[J]. Artificial intelligence, 2013, 205: 1-29.

本文引用格式：

刘威, 陈先桥，初秀民. 基于置信规则推理方法的雷达目标跟踪[J]. 哈尔滨工程大学学报， 2016, 37(6): 826-831.

LIU Wei，Chen Xianqiao，CHU Xiumin. Radar target tracking via belief rule-based methodology [J]. Journal of Harbin Engineering University, 2016, 37(6): 826-831.

收稿日期：2015-03-11.

基金项目：国家自然科学基金面上项目(51179146).

作者简介：刘威(1987-),男, 博士研究生； 通信作者：陈先桥，E-mail: chenxq1121 @163.com.

DOI:10.11990/jheu.201503029

中图分类号：TP391

文献标志码：A

文章编号：1006-7043(2016)06-0826-06

Radar target tracking via belief rule-based methodology

LIU Wei1,2,3，CHEN Xianqiao1,2,3，CHU Xiumin2,3

(1. School of Computer Science and Technology, Wuhan University of Technology, Wuhan 430063, China; 2. Intelligent Transportation System Research Center, Wuhan University of Technology, Wuhan 430063, China; 3. Engineering Research Center for Transportation Safety, Ministry of Education, Wuhan 430063, China)

Abstract：To address the data association problems in radar target tracking, a method of associating radar data according to radar target course and speed is proposed. It aims to troubleshoot erroneous targeting or target loss in the event of a large density of targets. The method adopts the belief rule-based (BRB) method to establish the rule base of the navigation velocity and course of radar targets. Through evidential reasoning it can mitigate the conflict between different dimensional data and thus proves reliable. The experimental results show that it can reach a real objective confidence of more than 95%. Finally, the real target points obtained by this method are filtered using a Kalman filter to track the radar targets. Monte Carlo simulation verifies that this method shows excellent radar target tracking.

Keywords：radar target tracking; belief rule-based; evidential reasoning; azimuth angle; radial velocity

网络出版日期：2015-04-21.

陈先桥(1961-), 男,教授,博士生导师.