Abstract

The major advantage of the passive multiple-target tracking is that the sonars do not emit signals and thus they can remain covert, which will reduce the risk of being attacked. However, the nonlinearity of the passive Doppler and bearing measurements, the range unobservability problem, and the measurement to target data association uncertainty make the passive multiple-target tracking problem challenging. To deal with the target to measurement data association uncertainty problem from multiple sensors, this paper proposed a batch recursive extended Rauch-Tung-Striebel smoother- (RTSS-) based probabilistic multiple hypothesis tracker (PMHT) algorithm, which can effectively handle a large number of passive measurements including clutters. The recursive extended RTSS which consists of a forward filter and a backward smoothing is used to deal with the nonlinear Doppler and bearing measurements. The target range unobservability problem is avoided due to using multiple passive sensors. The simulation results show that the proposed algorithm works well in a passive multiple-target tracking system under dense clutter environment, and its computing cost is low.

1. Introduction

Passive multiple-target tracking has gained more and more attention in the fields of military and civilian, such as navigation, monitoring and early warning, and salvage [13]. How to discover targets timely and to track targets accurately becomes one of the hot topics. The aim of the multiple-target tracking is to estimate the expected targets’ states, such as position, velocity, and acceleration, from the linear or nonlinear measurements [46]. The advantage of the passive target tracking is that the sonar can remain covert, which will reduce the risk of being attacked. The challenges for underwater passive multiple-target tracking include that the measurements are usually nonlinear, the target range may be unobservable, and the measurement to target data association is complex [5, 6].

In general, the passive measurements include bearing, Doppler, and bearing rate. So the passive multiple-target tracking is a typically nonlinear tracking problem [715]. One method to handle the nonlinear measurements is the pseudo linearization estimation (PLE) algorithm [1618]. The other approach is the recursive nonlinear Bayesian filter and smoother, such as the extended Kalman filter (EKF) [1922], the unscented Kalman filter (UKF) [23, 24], the cubature Kalman filter (CKF) [11, 12], and the particle filter (PF) [16]. EKF locally linearizes the targets’ state and measurement equations using the first-order Taylor series expansion of the nonlinear transformations around the predicted target state. The merit of the EKF is that it is very easy to understand and apply to target tracking problem duo to its linear approximations to the nonlinear system or measurement function, and its computational cost is lower than that of other nonlinear filter, such as UKF and PF. The disadvantage of the EKF is that the tracking performance may be not good if the system model is seriously nonlinear or non-Gaussian.

The UKF uses a deterministic sampling method to capture targets’ posterior distribution of mean and covariance based on the unscented transform. So the UKF has better tracking capability than the EKF to some extent, but its computational cost is larger than that of the EKF. The PF algorithm forms the Monte Carlo approximation to the solution of the Bayesian filter. It uses a set of particle samples to approximate the target distribution. The PF is used in the cases that the tracking system is highly nonlinear or non-Gaussian. The Bayesian filters only estimate the state of the targets from the history of measurements. On the contrary, the Bayesian smoother can use the forthcoming measurements to estimate the current state of targets. Corresponding to the Bayesian filter, the Bayesian smoother includes Rauch-Tung-Striebel smoother (RTSS), extended RTSS, unscented RTSS, and particle smoother [25].

The other challenge for passive sonar multiple-target tracking is the complexity and particularity of the underwater environment, which is characterized by many of false alarms caused by reverberation and multipath effect [3, 4]. Furthermore, the low target detection will cause tracking uncertainty in target location. All of those will cause a data association uncertainty problem. To handle the data association uncertainty problem, some algorithms are proposed, such as the multiple hypothesis trackers (MHT) [26], the joint PDA filter (JPDA) [27], the probabilistic MHT (PMHT) [28, 29], the random finite set framework-based probability hypothesis density (PHD) [30], the cardinalized PHD (CPHD) [3133], and multi-Bernoulli filter (MBF). The MHT makes all the data association hypotheses probability optimal. The JPDA makes multiple hypotheses into a single hypothesis and performs the Kalman update with composite measurements. The PMHT is based on the expectation maximization (EM), which optimizes the multiple-target states’ maximum a posteriori (MAP) estimation [3436]. Different from the data association algorithms, the PHD, CPHD, and MBF are based on the random finite set theory, which makes all the measurements a measurement set and all the targets a target set.

The information entropy theories are also used to estimate the target states. The fuzzy -means clustering method based on maximum information entropy and the probabilistic data association filter (PDA) is proposed in [37], which uses a value optimized by the maximum information entropy to represent the measurement to target association probability. The multiple-target tracking problem is also solved by the maximum entropy intuitionistic fuzzy data association [38], cross entropy [39], maximum fuzzy entropy-based Gaussian clustering algorithm [40], entropy distribution and game theory based on the probability hypothesis density (PHD) method [41], maximum entropy fuzzy based on the fire-fly and PF [42], and the distributed cross entropy-based -generalized labelled multi-Bernoulli filter [43].

As for the target range unobservability problem, the target range is observable from bearing and Doppler measurements if and only if the bearings are not constant [44, 45]. So in order to avoid the range unobservability problem, this paper uses multiple passive sensors to track targets. In addition, the target tracking accuracy of using multiple sensors is better than that of the single sensor generally.

The most commonly used passive measurement is bearing. In this paper, in order to improve the multiple targets’ range observability and tracking performance, we introduce the nonlinear Doppler measurement and use multiple sensors. The extended RTSS method is used to deal with the nonlinear Doppler and bearing measurements. The batch recursive multiple-sensor PMHT algorithm is used to handle the measurement to target data association complexity problem.

The remainder of this paper is as follows. The passive multiple-target tracking system model and measurement model are given in Section 2. Section 3 develops the multiple-sensor PMHT algorithm which is suitable for multiple-target tracking under dense clutter environment. The simulation result is given in Section 4. At last, a summary is given in Section 5.

2. System Model and Measurement Model

We consider the passive multiple-target tracking problem in a two-dimensional space.

2.1. Tracking System Model

Assume there are targets in the tracking space. The target’s state is , where and are the location of target in the and coordinate, respectively, and and are the velocity of target . Assume the tracking time is from to . The overview of the multiple-target tracking using Doppler and bearing measurements is shown in Figure 1.

Assume that all the targets move according to the following constant velocity (CV) model or constant acceleration (CA) model [46]. where is the system process noise which is assumed as Gaussian white noise. is the system state transition matrix. For the CV and CA models, is denoted as and , respectively, and where is the sampling interval.

The process noise covariance matrices for the CV and CA models are denoted as and , respectively, which are given by where is the process noise intensity.

2.2. Measurement Model for Multiple Sensors

Assume there are passive sensors in the target tracking space, and the sensor’s states is , in which and are the location of the sensor and and are the velocity of passive sensor .

The Doppler and bearing measurements are nonlinear with respect to targets’ states and sensors’ states, which are given by in which is the measurement noise with a covariance matrix .

The Doppler and bearing measurement function is given by where is the sound speed in water, is targets’ radiation frequency, and is bearing from target and sensor .

3. Multiple-Sensor PMHT

3.1. PMHT for Multiple Sensors

The following notations are used in this section.

is the total number of passive observers.

is the total number of targets.

is the number of measurements.

is the total number of tracking time.

is the set of target states at time .

is the set of measurements at time .

is the set of measurement-target assignments at time .

is the set of the target state for time .

is the set of measurements for time .

is the set of measurement-target assignments for time .

The PMHT is a Bayesian framework-based batch recursive algorithm, which obtains the maximum a posteriori estimation of the target states based on the expectation maximization (EM) method [47].

Let indicate that measurement from sensor at time is associated with target . Assume the measurement to target assignment is independent from each other.

Define the prior probability of the measurement from the sensor and target as

The MAP estimate of is

In order to calculate the MAP, define the following auxiliary function as where is the number of EM iteration.

The goal is to maximize the auxiliary function over using an initialized target state . In each EM iteration step, the goal of the PMHT is to solve the following equation.

We have

According to the conditional probabilistic theory,

From (12), (13), and (14), the posterior association probability is in which

According to the conditional probabilistic theory,

Substituting equation (17) into equation (10), the defined auxiliary function is calculated as

The has no effect with regard to maximizing the auxiliary function in (18). So it can be removed from (18).

So the auxiliary function is given by

In order to maximize the auxiliary function, we compute its derivative and set the derivative to zero.

It turns out that the is given by applying the extended RTSS method, in which the synthetic measurement and are calculated by

The implementation of the batch recursive multiple-sensor extended RTSS-based PMHT is depicted in Table 1.

3.2. Recursive Extended RTSS

The recursive extended RTSS is a linearized RTS smoother, which is based on analogous approximation to the EKF. It is a Gaussian approximation to the Bayesian smoother for the nonlinear target state and measurement model. The extended RTSS firstly performs the extended Kalman filter; then, a Kalman smoother is applied. That is, the extended RTSS consists of a forward filter and a backward smoother, which is summarized as follows:

Forward filter:

For ,

Time update

Measurement update where the filter gain is

Backward smoothing:

For , where

For the Doppler and bearing measurement, the measurement Jacobian matrix is given by where .

3.3. Data Fusion for Multiple Sensors

For the tracking problem using multiple sensors, we need to fuse the target and measurement information from multiple sensors. The most used data fusion methods are parallel fusion and sequential fusion [48].

In this paper, we use the sequential data fusion for multiple sensors which updates the target state and covariance by each sensor sequentially. The multiple sensors deal with the measurements one after another with an intermediate target state and covariance at one time.

Let and denote the target state estimation and corresponding covariance processed by the sensor . The Doppler and bearing measurements from the first passive sensor are used to calculate the first intermediate target state and covariance for each of the targets using a single sensor filter. Then, the measurements from the next passive sensor are used to update the intermediate target state and covariance as follows: where is the filter gain, and

The sequential fusion implementation of multiple sensors is shown in Figure 2.

4. Simulation

Assume the targets move in the two-dimensional space, and there are two static passive sensors. The two passive sensors deal with Doppler and measurements from targets and clutter. The two passive sensors located at (0, -3500) m and (0, -5500) m, respectively.

Assume that the clutter number is a Poisson distribution, and they are uniformly distributed in the Doppler and bearing measurement space. The average number of clutter in each sampling scan is 20 in the Doppler and bearing measurements domain. The target detection probability is the same for all targets; here, we set it to 0.8. The tracking time is 160 s with a sampling interval 1 s. The Monte Carlo run is 200. The process noise intensity is 1 m. The Doppler measurement noise variance is 1 Hz, and the bearing noise variance is .

4.1. Case of the CV Model

The five targets’ initial position and velocity for the CV model are given in Table 2.

For the CV model, the simulated initial position and velocity are generated from the position and velocity of ground truth with an estimation error such that and , where  m and  m/s.

The tracking scenario of true target trajectories and PMHT-estimated trajectories for the CV model is shown in Figure 3. The true Doppler and measurements for two static passive sensors without clutter and the PMHT synthetic Doppler and bearing measurements under dense clutter are given in Figures 4 and 5, respectively.

In Figure 3, the extended RTSS-based PMHT algorithm can track all the five targets effectively; even if the target 1 and target 4 tracks are away from the true trajectory, the proposed algorithm can get the true target track after several time scans in the middle tracking time scan.

As seen in Figures 4 and 5, the density of Doppler and bearing measurement from clutter is high, and the synthetic Doppler and bearing measurement of the proposed algorithm broadly consists of the true Doppler and bearing measurements without clutter, which means that the proposed algorithm has good declutter ability.

The position and velocity RMSE of the five targets are given in Figures 6 and 7. Accordingly, the average position and velocity RMSE are shown in Tables 3 and 4. As shown in Figures 6 and 7, the position RMSE and velocity RMSE for target 1 are increasing in the middle sampling scans, and as time goes on, the RMSE is decreased to a low lever, which is similar to other targets. This also can be seen from Tables 3 and 4. The average position and velocity RMSE is small which can meet the tracking accuracy requirement.

This paper uses the average normalized estimation error squared (ANEES) to evaluate the consistency of the proposed algorithm. For one target, the ANEES is defined as follows: where is Monte Carlo runs, is the true target state, is the estimated target state, and is the target state covariance.

The ANEES for fives targets is shown in Figure 8. As can be seen in Figure 8, the consistency of the proposed algorithm is good.

4.2. Case of the CA Model

The four targets’ initial position, velocity, and acceleration for the CA model are given in Table 5.

For the CA model, the simulated initial position , velocity , and acceleration are generated from the position , velocity , and acceleration of ground truth with an estimation error such that , , and , where  m,  m/s, and  m/s2.

The tracking scenario of the CA model for target true trajectories and the estimated tracks is shown in Figure 9. Similar to the CV model, the estimated tracks are consistent with the true targets’ trajectories.

The position RMSE and velocity RMSE for the CA model are shown in Figures 10 and 11. As shown in Figures 10 and 11, the position RMSE and velocity RMSE for target 3 are increasing in the middle sampling scans, and as time goes on, the RMSE is decreased to a low lever. The average position and velocity RMSE for the CA model is small which can meet the tracking accuracy requirement.

5. Conclusion

The major advantage of the passive sonars multiple-target tracking is that the sonars do not emit signals, and thus they can remain covert, which will reduce the risk of being attacked. But there are also challenges. Firstly, the Doppler and bearing measurements are nonlinear which makes the multiple-target tracking difficult. Secondly, the target states may be unobservable. Thirdly, the underwater environment is with dense clutter which will cause the measurement to target data association uncertainty problem. To deal with those problems, this paper proposed the extended RTSS-based batch PMHT method for multiple sensors and applied it to the passive multiple-sensor tracking system under dense clutter environment. This paper uses the extended RTSS algorithm to handle the nonlinear Doppler and bearing measurements. Multiple passive sonars are used to avoid the target state range unobservable problem. The multiple-sensor batch PMHT is used to deal with the data association uncertainty problem under dense clutter. The experiment results demonstrated that the proposed extended RTSS-based multiple-sensor PMHT algorithm can track multiple targets efficiently in the dense clutter environment, and the computing time is low.

Data Availability

The simulation condition data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (61703333, 62076201, and U1934222), the Natural Science Basic Research Plan in Shaanxi Province of China (2019JQ-746, 2019-740), the Key Laboratory of Shaanxi Provincial Department of Education (20JS088), and the Science and Technology Project of Beilin District (GX2017).