Next Article in Journal
Exploiting Big Data for Experiment Reporting: The Hi-Drive Collaborative Research Project Case
Previous Article in Journal
Smarter Evolution: Enhancing Evolutionary Black Box Fuzzing with Adaptive Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Cooperative Target State Estimation and Tracking Optimization Method of Multi-UUV

1
College of Intelligent Systems Science and Engineering, Harbin Engineering University, Harbin 150001, China
2
Qingdao Innovation and Development Base, Harbin Engineering University, Qingdao 266000, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(18), 7865; https://doi.org/10.3390/s23187865
Submission received: 27 April 2023 / Revised: 28 May 2023 / Accepted: 1 September 2023 / Published: 13 September 2023
(This article belongs to the Section Navigation and Positioning)

Abstract

:
This work studied two sub-problems of the cooperative state estimation and cooperative optimization of tracking paths in multiple unmanned underwater vehicle (multi-UUV) cooperative target tracking. The mathematical model of each component of the multi-UUV cooperative target tracking system was established. According to the target bearing-only information obtained by each unmanned underwater vehicle’s (UUV) detection, the extended Kalman filter algorithm based on interacting with multiple model bearing-only data was used to estimate the target state in a distributed way, and the federal fusion algorithm was used to fuse the estimated results of each UUV. The fused target state was predicted, and, based on the predicted target state, to achieve the persistent tracking of the target, the particle swarm optimization algorithm was used for the online collaborative optimization of the UUV tracking path. The simulation results showed that the multi-UUV distributed fusion filtering algorithm could obtain a better target state estimation effect, and the online path collaborative optimization method based on the prediction of the target state could achieve persistent target tracking.

1. Introduction

Target tracking is a hot research direction in the field of UUV applications. In the underwater environment, the UUV mainly depends on its own passive sonar to judge the relative orientation information of a target according to the noise radiation of the target [1]. Under the background of a complex underwater environment, less detection information, and uncertain target motion, it is often very difficult to persistently track a target by only relying on a single UUV. However, multi-UUV cooperative target tracking makes up for these shortcomings. Multi-UUV cooperatively estimates the target state and fuses the estimated results to achieve a high-precision estimation of the target’s motion state, which greatly reduces the target loss probability in the tracking process. While multi-UUV cooperative target tracking brings these advantages, it also brings two major challenges: one is how to fuse the estimation results of multi-UUVs; and the second is how to optimize the tracking path online according to the estimation results to achieve better detection [2].
In reference [3], a path-planning strategy combining the Lyapunov principle and a collision avoidance potential function was proposed to solve the problem of unmanned aerial vehicle (UAV) ground target tracking in the presence of obstacles and motion threats, which achieves collision avoidance while target tracking. Reference [4] used an interactive multi-model algorithm based on a second-order Markov chain to track maneuvering targets and verified the effectiveness of this algorithm compared to other multi-model algorithms. Reference [5] proposed a multi-UAV path-planning algorithm based on the spatial refined voting mechanism (SRVM) and particle swarm optimization (PSO). Reference [6] mainly studied a multiple autonomous underwater vehicle (AUV) system for cooperative tracking. Reference [7] analyzed the research status of cooperative estimation in the navigation field and presented a novel methodology of a fading cubature Kalman filter (CKF) with an augmented mechanism to address the problems involved in tightly coupled inertial navigation system/celestial navigation system (INS/CNS) integration for hypersonic vehicle (HV) navigation. Finally, it was verified that this method has a good performance in HV navigation with tightly coupled INS/CNS integration.
In this paper, a distributed target state estimation algorithm based on IMM-EKF is proposed to solve the problem of multi-UUV cooperative target tracking. Each UUV performs target state estimation in a distributed manner and the estimation results of each UUV are fused through a federal fusion structure. A collaborative path optimization method based on particle swarm optimization (PSO) is proposed to solve the problem of UUV tracking path optimization. One-step prediction is performed based on the fusion filtering results, the predicted target state is used as the path optimization input to optimize the UUV tracking path online in real time, and the tracking performance is improved.

2. Problem Description and Modeling

2.1. Problem Description

Firstly, the definition of the multi-UUV cooperative target tracking problem is given as follows:
Definition 1.
(Multi-UUV cooperative target tracking problem) assumes that there are multiple UUVs and moving targets with limited prior knowledge in the space. Each UUV carries sensors and can detect target information, which requires each UUV to cooperatively complete the target state estimation according to the detection information. According to the target state estimation result and each UUV space constraint, the tracking path is adjusted to complete the target tracking task.
Definition 1 can be modeled by mathematical formulation. The first is the mathematical model of fusion estimation:
max ( λ ( P ( X t [ k ] | Z i [ k ] , X t [ k 1 ] ) ) ) s . t . X t [ k ] = f ( X t [ k 1 ] ) + Q [ k ]     Z i [ k ] = h i ( X t [ k ] , X u , i [ k ] ) + ω [ k ] , i = 1 N ,
where P ( ) is the conditional probability description of multi-UUV cooperative target state fusion estimation, λ ( ) is the corresponding scalar function. X t [ k ] is the target state variable at time k , and X t [ k 1 ] is the target state variable at time k 1 . X u , i [ k ] is the state variable of the i -th UUV at time k . f ( ) is the target state transition function. Z i [ k ] represents the detection information of the i -th UUV at time k . h i ( ) is the nonlinear detection function of the i -th UUV. Q [ k ] represents the process noise of the target model and ω [ k ] is the detection noise of the sonar. N is the total number of UUVs executing tracking tasks.
The result of the target state fusion estimation will be used as the basis for multi-UUV tracking path optimization. The tracking optimization problem can be described as:
X c o m [ k ] = max ( i = 1 N J s e n ( X t [ k ] , X u , i [ k ] ) ) ,
where X c o m [ k ] is the instruction sequence of the UUV at time k and J s e n ( ) is the efficiency function of sonar detection, which is described as follows:
J s e n ( X t [ k ] , X u , i [ k ] ) = 1 , 0 , T h e   t a r g e t   i s   n o t   l o s t T h e   t a r g e t   h a s   b e e n   l o s t .
The above mathematical model is understood as follows. The purpose of multi-UUV cooperative target state fusion estimation is, under the condition that the target state X t [ k 1 ] at time k 1 and multi-UUV detection information Z i [ k ] at time k are known, that the estimation result X t [ k ] of the target state at time k is optimal. The control instruction sequence of a multi-UUV solved by the optimization algorithm should strive to ensure that more UUVs meet the requirements of persistent tracking. In other words, on the premise that the target is detected by as many UUVs as possible, the optimal control sequence for a multi-UUV is solved.
According to Definition 1, the multi-UUV cooperative target tracking problem in this paper can be described as follows:
First of all, the tracking plane is two-dimensional and the target is mobile. In the process of motion, the target may change its motion mode at any time. The detection method is bearing-only detection. The UUV carries passive sonar, which can detect the relative bearing information of the target. The tracking mode is active tracking, and the UUV will adjust the tracking path online according to the target state fusion estimation results during the tracking process. The requirement of the tracking task is persistent tracking, that is, to keep the target within the detection range of each UUV at all times. The schematic diagram of multi-UUV cooperative target tracking is shown in Figure 1.

2.2. Dynamic and Kinematic Models of UUV

Assuming that the UUV only moves in a two-dimensional plane, the following dynamic equation can be established for the UUV [8]:
u = F m 11 d 11 m 11 u + m 22 m 11 v r v = m 11 u r m 22 d 22 m 22 v r = T m 33 d 33 m 33 r + m 11 m 22 m 33 u v ,
where m 11 = m X u , m 22 = m X v , m 33 = I z N r , d 11 = X u , d 22 = Y v , d 33 = N r , u , and v are the longitudinal and transverse velocities of the UUV, respectively; r is the bow angular velocity of the UUV; m is the mass of the UUV; I z is the fixed torque of the UUV on the shaft; Z is the longitudinal thrust of the UUV propeller; T is the UUV bow torque; and X , Y , and N are the viscous hydrodynamic coefficients.
The UUV dynamic model describes the relationship between force, torque, and UUV velocity. In order to describe the relationship between the UUV’s velocity and position, and between the bow’s angular velocity and heading, a kinematic model of the UUV is established as follows:
x ˙ = u cos ψ v sin ψ y ˙ = u sin ψ + v cos ψ ψ ˙ = r ,
where x and y are the positions of the UUV in the geodetic coordinate system and ψ is the heading of the UUV.

2.3. Target Motion Model

Assume that the target only moves in the two-dimensional plane, and its state variable is X = x t x t y t y t T . The mathematical description of the target motion is:
X k + 1 = F X k + G Q k ,
where F is the state transition matrix, G is the control matrix, and Q k is the Gaussian white noise with a mean of 0 and covariance of σ 2 . The state transition matrix includes the uniform motion model and uniform turning model. Its matrix representation is as follows:
F C V = 1 0 T 0 0 1 0 T 0 0 1 0 0 0 0 1 ,   G C V = T 2 2 0 0 T 2 2 T 0 0 T , F C T = 1 sin ω T ω 0 1 cos ω T ω 0 cos ω T 0 sin ω T 0 1 cos ω T ω 1 sin ω T ω 0 sin ω T 0 cos ω T ,   G C T = T 2 2 0 0 T 2 2 T 0 0 T
where T is the sampling time and ω is the turning angular velocity.

2.4. Bearing-Only Detection Model

It is assumed that the UUV carries passive sonar, and the detection mode is bearing-only detection. The detection information does not consider other detection information except the relative azimuth, and the detection plane is in the two-dimensional plane, without considering the depth factor. The detection angle of the UUV’s passive sonar at time k can be described as follows:
Z [ k ] = arctan x t x u y t x u + ω [ k ] ,
Among them, x t and y t are the target positions at k , x u and y u are the UUV positions at k , and ω [ k ] is the detection noise of the passive sonar at k .

3. Fusion Estimation and Tracking Optimization Solving Framework

3.1. Tracking Optimization Approach Based on Target Prediction States

According to the different target states used for UUVs, the tracking optimization approach includes two aspects: an approach based on current target states and an approach based on target prediction states. The tracking optimization approach based on target prediction states is adopted in this paper.
In the process of target tracking, if only the target’s current states are considered, there is a great risk of target loss. As shown in Figure 2, assuming that path optimization is performed every 3 s, the target may be out of the detection range of the UUV after two seconds, and the UUV still executes the instructions obtained based on X t [ k ] .
According to the above problem, some studies have addressed this predicament through a diverse range of techniques, such as random weighting, Sage windowing, and adaptive factor adjustments [9,10]. In this paper, another solution is reached on the basis of the current target state fusion estimation results to predict the state of the target at the future time based on the target prediction state to optimize the tracking path of the UUV, which is shown in Figure 3. X t [ k + 1 ] is the target prediction state on the k + 1 moment based on the estimation results X t [ k ] . In addition, the tracking structure based on the target predicted state is shown in Figure 4, which greatly reduces the probability of target loss in the tracking process.

3.2. Finite Centralized Distributed Solving Framework

There are three commonly used control structures in multi-UUV cooperative control systems: centralized, distributed, and hybrid [11]. In the centralized architecture, there is a central node, and all UUV nodes communicate with the central node directly. In the distributed control structure, there is no central node, and UUVs communicate with each other. In the hybrid structure, there is a manual control center, and all UUVs communicate with the control center and can also communicate with each other.
At present, the multi-UUV system is not advanced enough to perform tasks completely autonomously without central intervention, so the existence of a central control node is inevitable. At the same time, it is hoped that each UUV has a certain autonomous ability, and can communicate and cooperate equally. Therefore, the hybrid control structure is an ideal control structure applied to multiple UUV systems. Combined with the above closed-loop tracking task process based on predicting the target state, a finite-set distributed solution framework for solving the multi-UUV cooperative target tracking problem is given, as shown in Figure 5.
By solving the framework shown in Figure 5, it can be seen that each is a closed loop. Multiple UUVs communicate in the task execution process through collaborative target state estimation and tracking path optimization, and each UUV in the task process is linked, in both the artificial control of monitoring and management. In this solving process, there are two main points of collaboration between UUVs:
(1)
The target state information is exchanged between the UUVs to complete the fusion estimation of the target state;
(2)
According to the respective position information and target estimation information fusion between the UUVs, the tracking path optimization is completed.
The two collaborations between UUVs also reflect the two keyword problems of solving the multi-UUV cooperative target tracking problem, the multi-UUV cooperative target state fusion estimation problem, and the multi-UUV tracking path cooperative optimization problem.

4. Distributed Cooperative Target State Fusion Estimation Method

4.1. Structure of Distributed Fusion Estimation Based on IMM-EKF and Federated Fusion

In the process of tracking a target, the sonar detects the orientation information of the target. According to the bearing-only detection model of UUV sonar, the observation equation is nonlinear, so the nonlinear filtering method should be used to estimate the target state. Common nonlinear filtering methods are: extended Kalman filter, unscented Kalman filter, and particle filter [12]. Considering the practical application of UUV target tracking and the advantages and disadvantages of the various filtering methods, the extended Kalman filter algorithm with a simple algorithm, strong engineering implementation ability, and better real-time performance is selected in this paper [13,14].
In the process of target tracking, due to the variable motion of targets, it is often difficult for a single model to match the real motion of a target. The interacting multiple model (IMM) is a good and computationally efficient method for dealing with the model-matching problem of maneuvering target tracking. Reference [15] applied the IMM algorithm to the vehicle state estimation problem and achieved good results. Combined with the extended Kalman filter algorithm, this paper uses the extended Kalman filter algorithm based on the interacting multiple model (IMM-EKF) to estimate the target state.
The choice of fusion structure also has an important impact on the problem of multi-UUV cooperative estimation. There are three commonly used fusion structures: a distributed structure, a centralized structure, and a hierarchical structure [16]. According to the previous problem description, the task requirement of the multi-UUV cooperative target tracking studied in this paper is persistent tracking, that is, ensuring that the target is within the detection range of each UUV at all times. Assuming that the detection range of the UUVs is equal to the communication range, communication between the UUVs must be possible on the premise of satisfying the persistent tracking. Because all UUVs can communicate with each other, each UUV can make use of the detection information of all UUVs when performing fusion estimation, and as long as each UUV uses the same fusion rule, the consistency of the estimation results is not considered.
In summary, the structure of distributed fusion estimation based on IMM-EKF and federated fusion is shown in Figure 6.

4.2. IMM-EKF Estimation Method

When the target state estimation is conducted and the distributed target state estimation method is used, each UUV according to itself detects the target information independently of the target state estimation based on the IMM-EFK algorithm, which is passed through mutual communication between the estimated results and estimated covariance, finally according to their own estimation information and the other UUV estimation information fusion estimation. The specific steps of the IMM-EKF algorithm are as follows.
(1)
Input interaction
We first calculate the mixing probability from model i to model j as follows:
u i j ( k 1 | k 1 ) = i = 1 r u i ( k 1 ) / c j _ ,
where the normalization constant c j _ is calculated by the following equation:
c j _ = i = 1 r p i j u i ( k 1 ) ,
where r represents the number of models in the target motion model set, p i j represents the transition probability from model i to model j , and u i ( k 1 ) represents the probability of model i at time k 1 .
The mixed-state estimates and mixed covariance estimates for model j can be further obtained:
X 0 j ^ ( k 1 | k 1 ) = i = 1 r X i ^ ( k 1 | k 1 ) u i ( k 1 | k 1 ) ,
P 0 j ( k 1 | k 1 ) = i = 1 r u i j ( k 1 | k 1 ) { P i ( k 1 | k 1 ) + [ X i ^ ( k 1 | k 1 ) X 0 j ^ ( k 1 | k 1 ) ] [ X i ^ ( k 1 | k 1 ) X 0 j ^ ( k 1 | k 1 ) ] T } ,
where j = 1 , 2 , , r , X i ^ ( k 1 | k 1 ) , and P i ( k 1 | k 1 ) are the state and covariance estimates of model j at time k 1 , respectively.
(2)
Extended Kalman filter
Taking X i ^ ( k 1 | k 1 ) , P i ( k 1 | k 1 ) , and probe information Z as the inputs, the extended Kalman filter is performed.
State prediction:
X j ^ ( k | k 1 ) = F j ( k 1 ) X 0 j ^ ( k 1 | k 1 ) ,
Error covariance prediction:
P j ( k | k 1 ) = F j P 0 j ( k 1 | k 1 ) F j T + G j Q j G j T ,
The Jacobian matrix is calculated. Under the bearing-only detection model, the target observation equation is nonlinear, and the Jacobian matrix needs to be solved to linearize the approximation. The Jacobian matrix is solved as follows:
H = Z x t , Z x t , Z y t , Z y t , Z x t , Z y t ,
calculate the Kalman gain:
K j ( k ) = P j ( k | k 1 ) H T H P j ( k | k 1 ) H T + R 1 ,
status updates:
X j ^ ( k | k ) = X j ^ ( k | k 1 ) + K j ( k ) Z ( k ) a r c t a n x t x u y t y u ,
covariance update:
P j ( k | k ) = I K j ( k ) H ( k ) P j ( k | k 1 ) ,
where F j ( k 1 ) is the state transition matrix in the target motion model set, G j is the control matrix, Q j is the process noise, R is the detection noise, x t , x t , y t , y t , x t , and y t are the state variables of the target, x t and y t are the target predicted positions obtained according to the target model, and x u and y u are the UUV positions.
(3)
Model probability update
Calculate the likelihood function for model j :
Λ j ( k ) = 1 ( 2 π ) n / 2 | S j ( k ) | 1 / 2 exp 1 2 v j T S j 1 ( k ) v j ,
among them:
v j ( k ) = Z ( k ) H ( k ) X j ^ ( k | k 1 ) ,
S j ( k ) = H ( k ) P j ( k | k 1 ) H ( k ) T + R ( k ) ,
the probability of model j is:
u j ( k ) = Λ j ( k ) c j _ / c ,
where the normalization constant is:
c = j = 1 r Λ j ( k ) c j _ ,
(4)
Output interaction
Based on the model probabilities, the final state estimate and covariance are obtained using the weighted summation of the estimated results for each filter as follows:
X ^ ( k | k ) = j = 1 r X j ^ ( k | k ) u j ( k ) ,
P ( k | k ) = j = 1 r u j ( k ) P j ( k | k ) + X j ^ ( k | k ) X ^ ( k | k ) X j ^ ( k | k ) X ^ ( k | k ) T .

4.3. Federated Fusion Estimation Method

In order to improve the accuracy of the target state estimation with multi-UUV, this paper adopts the collaborative estimation method of federal fusion. After each UUV completes the state estimation, it transmits the obtained target state and covariance information to all the other UUVs by broadcasting. Each UUV performs the weighted fusion of the received target state information and the target state information obtained via its own estimation according to the corresponding covariance. The specific fusion rules are as follows:
X g ^ = P g i = 1 N P i i 1 X i ^ ,
P g = i = 1 N P i i 1 1 ,
where N denotes the number of sub-filters, x i ^ denotes the local state estimate of the i -th sub-filter, and p i i denotes the corresponding covariance. x g ^ is the final target state estimation result after fusion.
Finally, the distributed cooperative target state fusion estimation algorithm based on IMM-EKF can be described as follows:
Finally, the distributed cooperative target state fusion estimation algorithm based on IMM-EKF can be described in Table 1:

5. Tracking Path Cooperative Optimization Method Based on PSO

5.1. Optimization Function Design for Persistent Target Tracking Task

The goal of the target tracking task is not to lose the target, that is, to keep the target in the detection range of the UUV sensor, and on this premise, to optimize the tracking path of the UUV to achieve the better detection of the target and reduce the target state estimation error. Firstly, according to the task requirements and practical engineering applications, the objective function of the optimization problem, namely, the fitness function of the particle swarm optimization algorithm, can be designed as follows:
f = J p e r s + J c r o s + J s a f e ,
The fitness function is composed of three parts, and each part is a penalty function. In the process of seeking the optimal solution in line with the optimization objective, the smaller the fitness function value of the particle position, the smaller the penalty value of the position, indicating that the position is better. The detailed description of the three-part penalty function is as follows:
(1)
Persistent tracking penalty function
J p e r s [ k ] = i = 0 N ω s e n i ,
J p e r s [ k ] represents the persistent tracking penalty function value at time k . When the penalty function value is 0, it means that all UUVs can observe the target. N represents the number of UUVs of the tracker and ω s e n i is the persistent tracking penalty factor. The expression is as follows:
ω s e n i = 0 , X t [ k ] Ω s e n i 1 , X t [ k ] Ω s e n i ,
where i = 1 N , N is the number of tracker UUVs, X t [ k ] is the position of the target at time k , Ω s e n i is the sensor detection range of the i -th UUV, X t [ k ] Ω s e n i means that the target at time k is within the sensor detection range of the i -th UUV, and X t [ k ] Ω s e n i means that the target at time k is not within the sensor detection range of the i -th UUV.
(2)
Path-crossing penalty function
J c r o s ( X c o m [ k ] , X u u v [ k ] ) = i = 0 , j = 0 N ω c r o s , ( i j ) ,
In the formula, X u u v [ k ] represents the position of the UUV at time k , X c o m [ k ] represents the instruction position of the future UUV optimized by the particle swarm optimization algorithm at time k , N represents the number of tracker UUVs, and ω c r o s is the path-crossing penalty factor, specifically expressed as:
ω c r o s = 0 ,   R o u t e i R o u t e j = 1 ,   R o u t e i R o u t e j ,
where i , j = 1 N and i j , R o u t e i , and R o u t e j are the paths of the i -th UUV and j -th UUV, respectively. R o u t e i R o u t e j = means that the two UUV paths have no crossing and R o u t e i R o u t e j means that the paths of the two UUVs have crossings.
As for the specific judgment method for path crossing, the current position of the UUVs, and the corresponding instruction position of the particle swarm optimization income as a path segment, there are several trackers with several paths according to the quick rejection in mathematical experiments, judging the line segment intersection method to determine whether the paths intersect, and then judging whether the instruction of the particle swarm optimization algorithm is reasonable. The specific path-crossing judgment method process is described as follows:
Firstly, the current position of the UUV is taken as the starting point of the path segment and the corresponding instruction position of the UUV is taken as the end point of the path segment. Each path segment is represented by the starting point and the end point. Then, the path segment of the i -th UUV can be represented as R o u t e i = { X i [ k ] , X c o m i [ k ] } , so that there are several UUVs and path segments. Observe the cross-product of four points.
(3)
Instruction safety penalty function
J s a f e ( X t [ k ] , X c o m [ k ] ) = ω s a f e ,
In the formula, X t [ k ] is the position of the target at time k , X u u v [ k ] represents the position of the UUV at time k , and ω s a f e is the instruction safety penalty factor, as detailed below:
ω s a f e = 0 , 1 , X c o m [ k ] Ω s a f e X c o m [ k ] Ω s a f e ,
where Ω s a f e represents the safe area. The specific meaning of the safe area is that the safe distance should be kept between the UUV instructions, the safe distance should be kept between the UUV instructions and the target position, and the instruction should not be within the prohibited navigation area. X c o m [ k ] Ω s a f e indicates that the instruction position of the optimized UUV meets the safety requirements and X c o m [ k ] Ω s a f e indicates that the optimized instruction position does not meet the safety requirements.

5.2. Cooperative Optimization Algorithm for Tracking Path Based on PSO

Particle swarm optimization (PSO) is a swarm intelligence optimization algorithm based on birds searching for food [17]. The core of the particle swarm optimization algorithm is two formulas, namely, the velocity update formula and the position update formula:
V i D k + 1 = ω V i D k + c 1 r 1 ( P i X i D k ) + c 2 r 2 ( P g X i D k ) ,
X i D k + 1 = X i D k + V i D k + 1 ,
where V i D k + 1 is the D -dimensional velocity of the particle at time k + 1 , X i D k + 1 is the D -dimensional position of the particle at time k + 1 , ω is the inertia weight, c 1 and c 2 are the learning factors, and r 1 and r 2 are random numbers. In practice, the following Equation (33) is often used to update ω [18]:
ω = ( ω m a x ω m i n ) i t e r m a x i t e r i t e r + ω m i n ,
In the formula, i t e r m a x is the maximum iteration number of the algorithm, i t e r is the current iteration number, ω m a x is the initial inertia weight, ω m i n is the minimum inertia weight, and the value is generally ω m a x = 0.9 , ω m i n = 0.4 .
In this paper, combined with the tracking path optimization problem of three UUVs’ cooperative target tracking, the three UUVs are regarded as a particle, and each particle contains the two-dimensional coordinates of the three UUVs, so the particle dimension is two-dimensional. The optimal solution of the optimization problem is the two-dimensional coordinates of the three UUVs that meet the task requirements. In each iteration, each particle is evaluated by the fitness function designed in the previous section, and the particles are updated by Equations (34)–(36). Finally, the optimal solution is obtained or the minimum accuracy requirement is achieved through layer-after-layer iteration.
The detailed steps of the multi-UUV cooperative target tracking path optimization algorithm based on particle swarm optimization are shown in Table 2:
In this table, X u [ k ] is the current position of the UUV, N is for tracking the target number of UUVs, X t [ k ] is for the prediction of the target position, f i t ( i ) is the i -th particle fitness value of the current position, p a r t i l e ( i ) . p o s i t i o n is the current position of the i -th particle, p a r t i l e ( i ) . p o s i t i o n B e s t is the best position of the i -th particle, p a r t i l e ( i ) . f i t B e s t is the fitness value of the best position of the i -th particle, p o p u l a t i o n . p o s i t i o n B e s t denotes the best position of the population, p o p u l a t i o n . f i t B e s t denotes the fitness value of the best position of the population, p a r t i l e ( i ) . V is the current velocity of the particle, i t e r denotes the current iteration number of the algorithm, i t e r m a x is the maximum number of iterations of the algorithm, Δ is the accuracy of the current solution, and Δ m i n is the particle swarm optimization accuracy threshold. When the accuracy requirements are met, it is considered that the optimal solution of the problem has been found and the iteration can be ended. In this paper, according to the characteristics of the fitness function designed in Section 5.1, when the optimal fitness value of the population is unchanged for N consecutive beats, it is considered that the requirements have been met.

5.3. Target State Prediction Methods

In Section 3.1 of the paper, the closed-loop target tracking task process with two different path optimization bases was discussed. In the process of tracking path optimization, if the algorithm uses a long time or the UUV communication takes a long time, the path optimization based on the current target fusion estimation results may be at risk of target loss. It is more reasonable to predict the state of the target in the future based on the estimated state of the target fusion and use the prediction results as the basis for path optimization, so the method can greatly reduce the probability of target loss. The specific prediction method is as follows.
Firstly, given the number of models j in each IMM-EKF filter and the initial probability μ i [ 0 ] of the model, after each filtering, the model probability μ i [ k ] is updated. According to μ i [ k ] , the motion model F j [ k ] that most conforms to the real motion state of the target at the current time k is judged, and the target prediction state is:
X t [ k + T f o ] = F j [ k ] X t [ k ] ,
F j [ k ] is the target model with the maximum corresponding probability μ i j in the IMM-EKF algorithm, X t [ k ] is the fusion estimation state of the target at time k , and T f o is the prediction duration.

6. Simulations

In order to verify the feasibility of the method, this section combines the algorithms in Table 1 and Table 2 to simulate and verify the whole process of multi-UUV cooperative target tracking, and the specific simulation design is described below.
(1)
The total simulation time was 1500 beats and the sampling time was T = 0.5   s .
(2)
Target setting. The target state variable was set to X t = x t y t x t y t T and the CV model and CT model established in Section 2.3 were used for the target motion model. In the simulation process, the target performed a turning motion with an angular velocity π / 600 ( rad/s ) from 150 s to 225 s, a turning motion with an angular velocity π / 500 ( rad/s ) from 375 s to 400 s, a turning motion with an angular velocity π / 800 ( rad/s ) from 525 s to 600 s, a turning motion with an angular velocity π / 900 ( rad/s ) from 600 s to 675 s, and a straight-line motion with uniform speed in the rest of the time. The initial state variable of the target was set to X t [ 0 ] = [ 0 m , 2.5 m/s , 0 m , 2.3 m/s ] , and the motion plane was two-dimensional.
(3)
UUV setting. In the simulation, using three UUVs for target tracking, a PID controller was used to achieve UUV movement control, using a dynamic model for the UUV heading and X , Y directional thrust. Each UUV carries passive sonar, which can only detect the target’s bearing; the detection radius was set to 300 m, the detection azimuth was Z [ k ] , and the detection noise was Gaussian white noise with a mean of 0 and variance of ( 0.01 rad ) 2 . The UUVs can communicate with each other by underwater sound with a communication radius of 300 m.
(4)
The settings of the distributed cooperative target state fusion estimation algorithm based on IMM-EKF. Each UUV used the IMM-EKF algorithm to estimate the target state. Additionally, UUVs can be federally fused locally, according to the target information sent by other UUVs and the local target information. Among them, IMM-EKF contained CV and CT models, the initial model probability was μ [ 0 ] = [ 0.5 , 0.5 ] , and the Markov probability transition matrix of the system was set as follows:
P i j = 0.99 0.01 0.01 0.99 ,
The system noise for all three filters was set as follows:
Q [ k ] = 0.01 2 d i a g ( [ 1 m , 1 m , 1 m/s , 1 m/s ] ) .
(5)
The settings of the collaborative optimization algorithm for tracking a path based on particle swarm optimization. The particle swarm size was set to 200 particles, the maximum number of iterations was 1000, the learning factors were c 1 = 2.8 and c 2 = 1.3 , and the inertia weight was updated using Equation (36); the initial inertia weight was 0.9, and the minimum inertia weight was 0.4. The iteration was terminated when the optimal fitness of the population was unchanged for 100 consecutive beats.
Firstly, the current target fusion estimation state was used as the basis for path optimization, and the simulation results were as follows:
In the zoom-in position of Figure 7, it can be seen that UUV1 lost the target, which did not meet the task requirement of persistent tracking.
Then, the above parameters remained the same. Taking the target prediction state as the input basis for path optimization, the prediction time was set to half of the particle swarm optimization time, and the simulation result is shown as follows.
According to the analysis of the multi-UUV cooperative target tracking simulation diagram in Figure 8, the overall navigation path of the UUV was consistent with the target navigation path. It can be seen from the local enlarged image that the target was within the detection range of the three UUVs in the linear motion and the turning motion, which met the task requirements of persistent tracking, and the UUV navigated towards the command position optimized by the particle swarm at each time.
Figure 9 and Figure 10 show the estimated results and estimated root mean square error (RMSE) of the target position. It can be seen that the IMM-EKF algorithm could also obtain good estimation results for a single UUV when the target moved in a uniform straight line. However, when the target moved in a turn, the RMSE of the single UUV estimation results fluctuated greatly, and the estimation results were not ideal. However, the RMSE of the three UUV estimation results was significantly reduced after using the federal fusion filtering algorithm. This shows that the IMM-EKF algorithm could obtain better motion state estimation results when the target moved in a uniform straight line, but when the target turned, the single UUV estimation results fluctuated greatly. The federal fusion could make up for this defect. It could also obtain good estimation results when the target turned and moved.
Figure 11 shows the real-time distance between each UUV and the target. Figure 12 shows the real-time distance between the UUVs. As can be seen in Figure 12, there would be no collisions between the UUVs in the tracking process (due to the large range of distance variation between the UUVs, the coordinates of the points with small values are marked in the figure to show that there was no collision, even when the distance between the UUVs was small), which indicates that the fitness function designed in Section 5.1 played the effect of collision avoidance between the UUVs. Figure 11 shows that the distance between each UUV and the target during the tracking process was always less than 300 m, which realizes the purpose of persistent tracking.

7. Conclusions

The distributed cooperative target state fusion estimation algorithm based on IMM-EKF proposed in this paper can estimate target motion state according to bearing-only information well, obtain a smooth estimation result in a target’s straight motion and turning motion, and has a high estimation accuracy. The cooperative path optimization method based on the particle swarm optimization algorithm can obtain the path points that meet the requirements of persistent tracking according to the target fusion estimation results. The optimization function designed in this study not only completed the task requirements of persistent tracking, but also had a beneficial effect on UUV collision avoidance, such as the path-crossing penalty factor and instruction safety penalty factor. The simulation results show that the distributed cooperative target state fusion estimation algorithm based on IMM-EKF and cooperative path optimization method based on particle swarm optimization can ensure that multiple UUVs can track a target without losing it and avoid collisions when the target moves in a uniform straight line or a turning direction.

Author Contributions

Conceptualization, T.C. and Q.Q.; methodology, T.C. and Q.Q.; writing—original draft preparation, T.C. and Q.Q.; writing—review and editing, T.C. and Q.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United States Department of Defense. Unmanned Systems Integrated Roadmap FY2011-2036; United States Department of Defense: Washington, DC, USA, 2011.
  2. Wang, L.; Su, F.; Zhu, H.; Shen, L. Active sensing based cooperative target tracking using UAVs in an urban area. In Proceedings of the 2010 2nd International Conference on Advanced Computer Control, Shenyang, China, 27–29 March 2010; pp. 486–491. [Google Scholar]
  3. Peng, Z.; Chen, Z. Ground target tracking and collision avoidance for UAV based guidance vector field. J. Adv. Comput. Intell. 2015, 19, 277–283. [Google Scholar] [CrossRef]
  4. Lan, J.; Li, X.R.; Jilkov, V.P.; Mu, C. Second-order Markov chain based multiple-model algorithm for maneuvering target tracking. IEEE Trans. Aero. Electron. Syst. 2013, 49, 3–19. [Google Scholar] [CrossRef]
  5. Liu, Y.; Zhang, X.; Zhang, Y.; Guan, X. Collision free 4D path planning for multiple UAVs based on spatial refined voting mechanism and PSO approach. Chin. J. Aeronaut. 2019, 32, 1504–1519. [Google Scholar]
  6. Shinzaki, D.; Gage, C.; Tang, S.; Moline, M.; Wolfe, B.; Lowe, C.G.; Clark, C. A multi-AUV system for cooperative tracking and following of leopard sharks. In Proceedings of the 2013 IEEE International Conference on Robotics and Automation, Karlsruhe, Germany, 6–10 May 2013; pp. 4153–4158. [Google Scholar]
  7. Gao, B.; Li, W.; Hu, G.; Zhong, Y.; Zhu, X. Mahalanobis distance-based fading cubature Kalman filter with augmented mechanism for hypersonic vehicle INS/CNS autonomous integration. Chin. J. Aeronaut. 2022, 35, 114–128. [Google Scholar]
  8. Do, K.D.; Pan, J. Control of Ships and Underwater Vehicles: Design for Underactuated and Nonlinear Marine Systems; Springer: London, UK, 2009. [Google Scholar]
  9. Gao, S.; Wei, W.; Zhong, Y.; Subic, A. Sage windowing and random weighting adaptive filtering method for kinematic model error. IEEE Trans. Aero. Electron. Syst. 2015, 51, 1488–1500. [Google Scholar] [CrossRef]
  10. Gao, B.; Gao, S.; Zhong, Y.; Hu, G.; Gu, C. Interacting multiple model estimation-based adaptive robust unscented Kalman filter. Int. J. Control Autom. Syst. 2017, 15, 2013–2025. [Google Scholar] [CrossRef]
  11. Fu, X.; Pan, J.; Wang, H.; Gao, X. A formation maintenance and reconstruction method of UAV swarm based on distributed control. Aerosp. Sci. Technol. 2020, 104, 105981. [Google Scholar] [CrossRef]
  12. Saha, M.; Ghosh, R.; Goswami, B. Robustness and sensitivity metrics for tuning the extended Kalman filter. IEEE Trans. Instrum. Meas. 2014, 63, 964–971. [Google Scholar] [CrossRef]
  13. Zhu, X.; Gao, B.; Zhong, Y.; Gu, C.; Choi, K.S. Extended Kalman filter based on stochastic epidemiological model for COVID-19 modelling. Comput. Biol. Med. 2021, 137, 104810. [Google Scholar] [CrossRef] [PubMed]
  14. Zhu, X.; Gao, B.; Zhong, Y.; Gu, C.; Choi, K.S. Extended Kalman filter for online soft tissue characterization based on Hunt-Crossley contact model. J. Mech. Behav. Biomed. 2021, 123, 104667. [Google Scholar] [CrossRef]
  15. Xu, Y.; Zhang, W.; Tang, W.; Liu, C.; Yang, R.; He, L.; Wang, Y. Estimation of Vehicle State Based on IMM-AUKF. Symmetry 2022, 14, 222. [Google Scholar] [CrossRef]
  16. Lee, D.J. Unscented information filtering for distributed estimation and multiple sensor fusion. In Proceedings of the AIAA Guidance, Navigation and Control Conference and Exhibit, Honolulu, Hawaii, 18–21 August 2008. [Google Scholar]
  17. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the IEEE International Conference on Neural Networks, Piscataway, NJ, USA, 27 November–1 December 1995; pp. 1942–1948. [Google Scholar]
  18. Mojtahedi, A.; Baybordi, S.; Fathi, A.; Yaghubzadeh, A. A hybrid particle swarm optimization and genetic algorithm for model updating of a pier-type structure using experimental modal analysis. China Ocean Eng. 2020, 34, 697–707. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of multi-UUV cooperative target tracking.
Figure 1. Schematic diagram of multi-UUV cooperative target tracking.
Sensors 23 07865 g001
Figure 2. Tracking optimization approach based on target current states.
Figure 2. Tracking optimization approach based on target current states.
Sensors 23 07865 g002
Figure 3. Tracking optimization approach based on target prediction states.
Figure 3. Tracking optimization approach based on target prediction states.
Sensors 23 07865 g003
Figure 4. Tracking structure based on target predicted state.
Figure 4. Tracking structure based on target predicted state.
Sensors 23 07865 g004
Figure 5. Finite centralized distributed solving framework.
Figure 5. Finite centralized distributed solving framework.
Sensors 23 07865 g005
Figure 6. Structure of distributed fusion estimation based on IMM-EKF and federated fusion.
Figure 6. Structure of distributed fusion estimation based on IMM-EKF and federated fusion.
Sensors 23 07865 g006
Figure 7. Cooperative target tracking result of three UUVs (optimized based on current state).
Figure 7. Cooperative target tracking result of three UUVs (optimized based on current state).
Sensors 23 07865 g007
Figure 8. Cooperative target tracking result of three UUVs.
Figure 8. Cooperative target tracking result of three UUVs.
Sensors 23 07865 g008
Figure 9. Estimated result of target position.
Figure 9. Estimated result of target position.
Sensors 23 07865 g009
Figure 10. Estimated root mean square errors of target position.
Figure 10. Estimated root mean square errors of target position.
Sensors 23 07865 g010
Figure 11. Real-time distance between each UUV and target.
Figure 11. Real-time distance between each UUV and target.
Sensors 23 07865 g011
Figure 12. Real-time distance among UUVs.
Figure 12. Real-time distance among UUVs.
Sensors 23 07865 g012
Table 1. Distributed cooperative target state fusion estimation algorithm based on IMM-EKF.
Table 1. Distributed cooperative target state fusion estimation algorithm based on IMM-EKF.
Step 1 Initialization. Set each UUV X i ^ [ 0 ] , P i [ 0 ] , p i j , μ i [ 0 ] i = 1 , , N .
Step 2 IMM-EKF. Each AUV is executed:
  Step 2.1 Input interaction. Use Equations (8)–(11) to obtain the mixed input under
     each target model.
  Step 2.2 EKF. Use the mixed input obtained in the previous step to obtain the
     estimation results under each model according to Equations (12)–(17).
  Step 2.3 Model probability update. Obtain the model probability at the current
     moment μ i [ k ] .
  Step 2.4 Mixed output. The local estimation results of each UUV are obtained using
     Equations (23) and (24) X ^ i ( k | k ) and P i ( k | k ) .
Step 3 UUV communication. Each UUV sends X ^ i ( k | k ) and P i ( k | k ) to the other UUVs
  through underwater acoustic communication.
Step 4 Fusion estimation. Each UUV obtains the final target state estimation results X g ^
  and P g through Equations (25) and (26).
Step 5 Repeat Step 2–Step 4 until the task is over.
Table 2. Cooperative optimization algorithm of tracking path based on PSO.
Table 2. Cooperative optimization algorithm of tracking path based on PSO.
Step 1: Input X u [ k ] , N , X t [ k ] .
Step 2: Initialization population.
   Set c 1 , c 2 , ω , V m a x ;
   Set particle & population;
Step 3: f i t ( i ) = J p e r s + J c r o s + J s a f e
Step 4: if ( f i t ( i ) p a r t i l e ( i ) . f i t B e s t )
       p a r t i l e ( i ) . p o s i t i o n B e s t = p a r t i l e ( i ) . p o s i t i o n
      if ( f i t ( i ) p o p u l a t i o n . f i t B e s t )
         p o p u l a t i o n . p o s i t i o n B e s t = p a r t i l e ( i ) . p o s i t i o n
      end
     end
Step 5: Update p a r t i l e ( i ) . V & p a r t i l e ( i ) . p o s i t i o n according to Equations (34) and (35);
Step 6: if ( i t e r = i t e r m a x or Δ Δ m i n )
     Step 7;
     else
     Step 3;
     end
Step 7: Output X c o m .
Step 8: UUV execute X c o m .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, T.; Qi, Q. Research on the Cooperative Target State Estimation and Tracking Optimization Method of Multi-UUV. Sensors 2023, 23, 7865. https://doi.org/10.3390/s23187865

AMA Style

Chen T, Qi Q. Research on the Cooperative Target State Estimation and Tracking Optimization Method of Multi-UUV. Sensors. 2023; 23(18):7865. https://doi.org/10.3390/s23187865

Chicago/Turabian Style

Chen, Tao, and Qi Qi. 2023. "Research on the Cooperative Target State Estimation and Tracking Optimization Method of Multi-UUV" Sensors 23, no. 18: 7865. https://doi.org/10.3390/s23187865

APA Style

Chen, T., & Qi, Q. (2023). Research on the Cooperative Target State Estimation and Tracking Optimization Method of Multi-UUV. Sensors, 23(18), 7865. https://doi.org/10.3390/s23187865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop