Mathematical Problems in Engineering

Volume 2017 (2017), Article ID 5752870, 9 pages

https://doi.org/10.1155/2017/5752870

## Tracking Maneuvering Group Target with Extension Predicted and Best Model Augmentation Method Adapted

Air and Missile Defense College, Air Force Engineering University, Shaanxi, China

Correspondence should be addressed to Linhai Gan

Received 13 March 2017; Revised 3 June 2017; Accepted 14 June 2017; Published 24 September 2017

Academic Editor: Vladimir Turetsky

Copyright © 2017 Linhai Gan and Gang Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The random matrix (RM) method is widely applied for group target tracking. The assumption that the group extension keeps invariant in conventional RM method is not yet valid, as the orientation of the group varies rapidly while it is maneuvering; thus, a new approach with group extension predicted is derived here. To match the group maneuvering, a best model augmentation (BMA) method is introduced. The existing BMA method uses a fixed basic model set, which may lead to a poor performance when it could not ensure basic coverage of true motion modes. Here, a maneuvering group target tracking algorithm is proposed, where the group extension prediction and the BMA adaption are exploited. The performance of the proposed algorithm will be illustrated by simulation.

#### 1. Introduction

Groups are structured objects and formations of entities moving in a coordinated manner [1]. Limited by poor sensor resolution and little requirement in application, the objects to be tracked are considered as point sources in conventional sense [2]. With ever-increasing sensor resolution capability and requirements for tracking a crowd, a herd, or an air fleet, group target tracking has been a hot topic in recent years. An increasing capability of sensor technology makes it possible to learn more feature information about a group. However, there will be a fluctuation number of detections for a group of closely spaced targets, accounting for limited sensor resolution, which will prevent a successful tracking of the individual targets [2]. Thus, there are some new challenges in tracking a group target compared with tracking point one.

The methods for group target tracking mainly include tracking via Poisson likelihoods [3, 4], group’s evolution modeled by a dynamic Bayesian network [5], the approach with random matrices [2, 6–11], group tracking using random finite sets [12, 13], and so on. A random matrix method was proposed by Koch [8] in 2008, which described the extension of the group by a symmetrical positive definite (SPD) random matrix, and characterized by simple filter equations, a small number in preset parameters, strong robustness, and so on. On the basis of [8], an interacting multiple model (IMM) structure was introduced by Feldmann [2, 7], and the performance for maneuvering extended objects tracking was improved; moreover, with the measurement noise being considered in the innovation covariance, the application area of the algorithm was expanded. A suitable measurement model to address the radar’s measurement noise and its conversion from polar coordinates to Cartesian coordinates was proposed in [14]. A Bayesian framework was developed in [15], within which the probability density function (PDF) of the object state and extension and the probability mass function of the object class were jointly obtained; a joint tracking and classification algorithm for an extended object was proposed using random matrix consequently. Within the random matrix framework, a nonellipsoidal extended object (NEO) was approximated by multiple ellipsoidal subobjects, with each subobject described by a random matrix in [11, 16]. In [17], the orientation of the extension in the next time was predicted by the extension estimation of the last two moments; the assumption that the extension of the adjacent two moments keeps invariant was no longer valid, but there was no analysis in theory about this. A multiple group target tracking algorithm under clutter environment was proposed based on random matrix method in [18].

IMM approach is one of the most wildly used tracking algorithms for maneuvering target; with the interacting of models, fusion estimation for maneuvering target was obtained. With a fixed structure, many more models are needed for an efficient tracking performance, as the target maneuverings in a more complex mode, which will lead to an increase in computation burden and model competition consequently. One of the variable structure multiple model (VSMM) methods called best model augmentation (BMA) was proposed in [19], in which the state estimation model set was adjusted adaptively at each step, which will alleviate the problem of the IMM approach referred above.

The paper is organized as follows. Section 2 gives a derivation of a RM method for group target tracking under Bayesian framework, where the group extension prediction is available. Section 3 describes the BMA approach and its modified form. Section 4 outlines the procedure of the proposed algorithm. The performance of the algorithm is presented via simulations in Section 5. Finally, conclusions are given in Section 6.

#### 2. Bayesian Group Target Tracking

##### 2.1. Background

Considering a target group with an ellipsoidal shape, we use a SPD random matrix to describe the extension of scan . Consider a 2-dimensional extension, , can be expressed as , with and denoting the main axis and minor axis of the ellipse, respectively, and is a rotation matrix.

It is assumed that in each scan there is a random number of independent position measurements:where and denote the state to be estimated and the th sensor measurement of scan , respectively. Moreover, denotes the measurement matrix, is the th measurement noise and is assumed to be a zero mean normally distributed random vector with covariance , where is a scaling factor, is the predicted extended state, and is the covariance matrix of the sensor measurement error.

In a Bayesian view, a group target tracking algorithm is an iterative updating scheme for conditional probability densities at each time [8]. Its product representation can be written as

or

with its Bayes’ formula aswhere and are used to denote the set of the measurements in a particular scan and for the accumulated measurements set among the time. The conditional probability densities and describe the kinematical object properties and the extended properties, respectively. And denotes the likelihood function and is a predicted joint probability density function.

With the assumption that the extension does not tend to change over time in [8], that is, , moreover , where represents the estimation of the extended state. The joint prediction of kinematic state and extension can be derived from (2):

The prediction of kinematic state and extension can be expressed, respectively, as follows.

(A) Kinematic state is as follows: where and denote prediction and estimation of kinematic state, respectively; is a -dimension identity matrix; and denote the predicted covariance and estimated covariance, respectively; denotes the state transform matrix, that is, for a constant velocity model, , where is the sample period; denotes the covariance of process noise; stands for Kronecker products.

(B) Extension is as follows: where denotes a Wishart distribution function; denotes freedom degree, with and denoting the evolution parameters of extension; denotes the sensor sample period; denotes “Generalized Beta Type II” density.

In reality, the assumption that extension does not tend to change over time is not so valid in many situations, for example, when the group target is moving in a turning mode, the orientation of the group extension will change timely; a proper prediction to the extension may improve the tracking performance.

##### 2.2. Innovation Prediction

It can be learned heuristically from [2] that the joint prediction can be derived as

While the group target has a relatively large extension and a dense scattering, the kinematic state may have little influence on the extension. So, posterior probability density of the extension is nearly independent from the kinematic state, but mainly dependent on the measurement. Thus we make the following assumption.

*Assumption 1. *The posterior probability density of extension can be approximated by marginal probabilityThuswhere denotes the inverse Wishart density of* d*-dimension SPD random matrix , denotes freedom degree; , and is the estimation of the freedom degree. Thus, (8) can be approximated as

###### 2.2.1. Kinematic State

It can be learned from (3) and (11) that the prediction and estimation of extension and kinematic state can be calculated, respectively. Assume that the process noise is a zero mean normally distributed random vector with covariance . Thus, the prediction of kinematic state is the same as the standard Kalman filter’s; and can be given by

###### 2.2.2. Extension

Naval vessels, Submarines, or ground moving convoys show a clear orientation [9]. As for airplane formation, the angle between its moving orientation and the extension appearance will keep invariant (or a relatively small change) in general. Thus, the group’s extension can be predicted according to its moving orientation in many applications. Assume thatwhere denotes a predicted rotation matrix, while in a 2-dimensional situation, it satisfies: where denotes the prediction of the angle that the group target rotates in two adjacent moments. and denote the velocity in - and -axis, respectively, and are parts of kinematic state . According to Assumption 1, and are approximately independent. Thus, can be expressed as (15), according to the property of Wishart distribution: Thus the distribution of can be derived as [8]

##### 2.3. State Update

With the regularities of distribution being invariant, the state update procedure is similar to [2].wherewith and denoting the estimation of the kinematic state and its covariance, respectively.

*Extension Update*where denotes the estimation of the extension. And

#### 3. Improved BMA Approach

##### 3.1. Rationale of BMA

Models with fixed parameter are used to form the candidate model set of BMA, whose structure is different from that of the basic model set. The model who can best match the true motion mode of the target in the candidate model set will be activated. Moreover, the activated model will be selected to form a state estimation model set (SEMS) of a particular moment combination with the basic model set. All the estimations of the models in SEMS are fused to output the final estimation [19].

Assume that denotes the measurement sequence; the SEMS iswhere denotes a small set of necessary models and it is adopted to ensure a basic coverage of true modes; denotes the activated model set at , which is activated from to improve the estimation performance; denotes the candidate model set of .

In order to select the model to be activated**,** a Kullback-Leibler (KL) criterion will be adopted to quantify the difference between a candidate model and the true motion mode; the one who has minimum divergence in KL with true mode will be selected to replace the activated model set of the previous. The SEMS is updated consequently.

The discrepancy of model from true mode can be provided by a metric :where with denoting the matrix trace; denotes the sequence of model set through , and denotes the dimension of , with acting as a common variable of and . Particularly, denotes the measurement vector; denotes the measurement sequence. As the true model is unknown in practice, is approximated by where denotes the model transition probability.

If where OtherwiseThe model to be activated from is

##### 3.2. BMA Adaptation

As shown in Section 3.1, the SEMS of BMA are mainly composed of the basic model set, which keeps invariant all over the tracking period. This will bring the following two problems:①With the basic model set selected offline, more prior knowledge is needed to ensure a basic coverage of true modes, which is a difficult task in practice.②To cover a complicate maneuvering mode, the basic model set should be augmented, which will make only a part of model efficiently match the maneuvering at each moment; the other models just do little help in improving the tracking performance, moreover, will bring a heavy burden in computation and model competition.

To overcome the drawbacks of BMA method mentioned above, an adaptive structure is considered for a better performance in maneuvering tracking. Suppose thatwhere denotes the model set eliminated from , which is then added to form at ; denotes the optimal model set selected from , which is activated at . Equation (29) shows the update of the estimation model set and the candidate model set; moreover, the basic model set of step can be updated as .

The metric represents the discrepancy between model and true motion mode , which can be calculated just as (22). Moreover, will be selected from either or .

Ifthen

Otherwise, the activated model is

The model eliminated from SEMS is

#### 4. Fusion Procedure

*(**1) Reinitialization*. ① Mixed probability iswhere denotes the model probability of step; , with denoting the transition probability and denoting the number of model in the state estimation model set.

② Mixed estimation is

*(**2) Update*. ① The prediction , , can be calculated by (12)~(13).

② Likelihood function iswhere

Here is the abbreviation of .

③ state and covariance update is as follows.

The estimations , , , , and can be calculated according to Section 2.3.

With .

④ Model probability update is as follows:

*(**3) Fusion Estimation*

#### 5. Simulation Results

Assume that the initial kinematic state of the group is ; initial extended state . The parameter denotes measurement ratio. The measure sustains 40 s, with sampling interval and measurement error covariance . The group take a constant turning (CT) with , , and being angular speeds of time periods 9~16 s, 17~24 s, and 25~32 s, respectively, and move with a constant velocity (CV) mode in other time periods. The track of the group is as shown in Figure 1, where the black spots denote the real measurements of the group target; the ellipse shows the confidence region of the group extension with confidence level being 0.9.