1. Introduction
Autonomous underwater vehicles (AUVs) have been used for a variety of offshore commercial and scientific applications [
1,
2,
3,
4,
5,
6]. For these applications, accurate tracking of AUVs is essential to ensure the accuracy of the gathered data [
7]. Many underwater tracking approaches have been implemented and tested; however, traditional underwater tracking techniques, such as dead reckoning, suffer from unbounded tracking errors. The widely used acoustic approaches, such as long-base line (LBL) system, are restricted by the costly setting up since the location of each deployed beacon in the array must be precisely surveyed before conducting navigation operations. Therefore, tracking AUVs using a single beacon would provide dramatic time and cost savings for underwater vehicle operations, which has been proposed and studied in [
7,
8,
9,
10,
11,
12,
13,
14,
15].
In the single-beacon underwater tracking system, AUVs rely on range measurements from an acoustic beacon with known position to bound errors accumulated by dead reckoning. Early researches on single-beacon tracking can be traced back to Larsen, in which the positioning and velocity error of AUV were treated as state variables, and the extended Kalman filter (EKF) was applied to estimate the vehicle position over time [
16]. For positioning an AUV based on time of arrival (TOA) measurements from a single beacon, Casey et al. established a single-beacon tracking model with unknown ocean current using EKF as the state estimator [
17]. Later, Fallon developed an AUV- USV (unmanned Surface Vehicle) cooperative tracking model using EKF to update the AUV position [
18]. Considering the packet loss problem in underwater communication, Walls proposed a robust origin state method for range-based cooperative tracking based upon a single-beacon tracking scheme, and used the extended information filter (EIF) as the state estimator [
19].
In the single-beacon underwater tracking methods mentioned above, the slant range was commonly adopted as the primary measurement, which was obtained from the measured acoustic transit time in conjunction with a presumed effective sound velocity (ESV) for the acoustic propagation between the beacon and the receiver. A known constant ESV was routinely assumed in [
8,
9,
10,
11,
16,
17,
18,
19]. Since the ESV, in practice, is location dependent, and is very difficult to determine accurately, these methods suffer from large ESV measurement error and consequently has large positioning error. By treating the unknown ESV as a state variable and the transit time as a measurement, Zhu et al. proposed a novel single-beacon localization model (namely the “5-sv” model), and used the EKF as the state estimator [
20,
21,
22]. It was shown that by properly tuning the process and measurement noise parameters, the EKF based on the “5-sv” model could estimate the unknown ESV well, and the tracking performance of the “5-sv” model significantly outperform the traditional model.
By augmenting the unknown ESV as a state variable, the “5-sv” model not only introduced an additional parameter related to the ESV uncertainty needed to be tuned, but also increases the nonlinearity of the measurement model. It had been shown that the EKF performance depended largely on the initial setting of the ESV uncertainty parameter which would stay invariant during the whole estimation process [
20]. Furthermore, while implementing the “5-sv” model, the parameters of the process noise covariance matrix and the measurement noise covariance matrix were difficult to tune, and these were usually predetermined based on experience or data post-processing. However, the huge variation of underwater environment leads to a rapid change in the probabilistic properties of the acoustic measurement noise. The problem of EKF estimates degradation occurs frequently while implementing the “5-sv” model due to the violation of known noise statistics assumptions.
To overcome the divergence issue of EKF-based single-beacon tracking methods, this paper will propose an adaptive Kalman filter (AKF)-based single-beacon underwater tracking algorithm. Although the AKF has been intensively investigated to reduce the influence of process noise covariance matrix and measurement noise covariance matrix errors in fields like integrated navigation, economic projection, and chemistry, to the best knowledge of the authors, little literature has been published implementing AKF for single-beacon underwater tracking, especially based on the model with unknown ESV [
23,
24,
25]. Through simulation and field data, the estimation results between using the proposed adaptive algorithm and the traditional EKF will be compared. In addition to the real-time filter tracking, this paper also implements the Rauch-Tung-Striebel (RTS) smoother for post-processing [
26,
27]. Both simulation and field data will be used to study the possible improvement of the proposed adaptive single-beacon underwater tracking algorithm.
3. Windows-Based Adaptive Single-Beacon Tracking Algorithm
Unlike the EKF, the AKF can adjust the noise parameters online based upon the immediate measurements and had proven to be a good solution to improve the robustness of Kalman filters. In this section, a Windows-based adaptive single-beacon tracking algorithm based on the “5-sv” model will be proposed to solve the divergence issue of traditional EKF-based methods.
As presented in many textbooks [
26,
28], the discrete Kalman filter equations include the prediction and correction equations. Compared with EKFs, the Windows-based AKF introduced an additional step for estimating the process and measurement noise covariance matrices based upon the innovation sequence.
3.1. State Prediction
Since the kinematic equation of the “5-sv” model is linear, the prediction equation of the adaptive single-beacon tracking algorithm is the same as the standard Kalman filter, so the estimate of the state vector
and the corresponding covariance matrix
are
and
where
is the estimated process noise covariance matrix at the
th epoch.
Throughout this paper, a variable with a hat, such as , represents the estimate of the variable, superscripts “−” and “+” denote quantities associated with a priori(prediction) and a posteriori(correction) estimates, respectively.
3.2. State Correction
For the correction equations of adaptive single-beacon tracking algorithm, the estimate of the state vector
and the corresponding covariance matrix
are
and
in which
is the Kalman gain matrix and
is the estimated measurement noise covariance matrix. For the single-beacon underwater tracking algorithm discussed herein,
is the Jaccbian of nonlinear measurement equation
at the
a priori estimate
in which
is computed by
Noting that
is available from a depth sensor. The measurement innovation
in Equation (
20) is denoted by
in this paper.
3.3. Estimation of Unknown
From literature [
23], the main purpose of the Windows-based adaptive estimation is to estimate the
through the innovation sequence, in which the innovation is considered to be an intermediate variable. The fundamental principle is making the elements in the actual innovation covariance matrix consistent with their theoretical value. Their actual values are intractable; however, they can be estimated based on maximum likelihood (ML) criterion. The likelihood function is chosen to be the joint conditional probability of measurement given estimated variable within the fixed windows
W:
where
is the covariance matrix of innovation vector at time
.
Assuming
satisfies zero-mean Gauss distribution, Equation (
25) can be written as
where
means the vector
satisfies Gaussian distribution with mean
and covariance matrix
. The probability density function is
The logarithmic form of Equation (
26) is taken as
To simplify the expression, the following hypothesis is made:
Hypothesis 1. For a Windows-based AKF, the innovation covariance matrix is constant within the windows width W, that is After multiplying Equation (
28) by
and neglecting the constant term, the ML criterion of maximizing
J within the windows width
W becomes the minimization problem which is described as
The minimum value is calculated by the partial of
with respect to
To obtain the above formula, the following two relations from matrix differential calculus have been used
Based on the relations from vector inner product
and the trace property of matrix, Equation (
31) can be written as
Let
then the estimated innovation covariance matrix is obtained as
On the other hand, the theoretical innovation covariance matrix can be calculated through the definition of innovation, and have
Let the estimated value of the innovation covariance matrix equals to the theoretical value, the estimated
can then be obtained
This estimated will be used for calculating the Kalman gain at epoch k.
3.4. Estimation of Unknown
To estimate the unknown
, rewritten Equation (
22) as
where
is assumed to predetermined from Equation (
36). Since
and
are symmetric matrices, Equation (
39) can be written as
Plugging the left side of Equation (
40) into Equation (
21) yields
From Equations (
41) and (
19), the estimated
can be get
To calculate , the following hypothesis is made:
Hypothesis 2. The varying period of is much larger than the sampling period , thus we have Based on the Hypothesis 2, Equations (
36) and (
42), the estimated
at time
can be obtained
Algorithm 1 summarizes the detailed procedure of the Windows-based adaptive single-beacon underwater tracking algorithm.
Algorithm 1 Windows-based adaptive single-beacon underwater tracking algorithm |
Input: Prediction:- 1:
- 2:
Correction:- 1:
- 2:
- 3:
- 4:
- 5:
- 6:
- 7:
- 8:
Return: |
4. Numerical Studies
Both simulation and field data are used to evaluate the performance of the Windows-based adaptive single-beacon underwater tracking algorithm (to be referred to as the “ASB” algorithm) against that of the EKF-based algorithm (to be referred to as the “ESB” algorithm). In addition to the real-time filter tracking, this paper also implements the RTS smoother for post-processing. The RTS smoother is the most commonly used fixed-interval smoother, which operates backward in time steps after the complete filtering solution has been obtained [
26,
27].
4.1. Simulation Data
The reader is referred to [
20] regarding the simulation technique for generating TOA data. For consistency, the trajectory in simulation was chosen to be the same as the trajectory of field data (shown in
Figure 1). For the simulation of TOAs, the ESV has been set equal to a constant 1530 m/s during the whole trajectory, with TOEs every 10 seconds. Both ocean current components in the
x and
y directions are
m/s. The process noise parameter
is chosen as
m/s.
and
are both set to
m/s. For clarity, the estimation performance of the unknown
and
are studied in two distinct scenarios, respectively. The
is assumed to be precisely available during the estimation of unknown
, and vice versa.
The sampling frequency of each simulated measurement and the corresponding noise level characterized by its standard deviation are summarized in
Table 1. Notice that the sampling interval of the hydrophone is irregular due to varying acoustic propagation times because of vehicle motion.
In the application of the ASB and the ESB algorithms, the following initial settings have been chosen: (1) a total of 10 meters in both the x and y directions for the initial position offset; (2) m/s for the initial ESV; (3) m/s for the initial ocean current in both x and y directions; and (4) m/s; (5) m/s. Specifically, the innovation window width W is chosen as 10 while implementing the ASB algorithm.
Fifty Monte Carlo simulations have been done to compare the performance of the ESB and the ASB algorithm. The root mean square (RMS) of the horizontal distance error
and the estimation error of the ESV
are chosen as the evaluation indexes of the Monte Carlo simulations, which are defined respectively as:
in which
and
represent respectively the target and the estimated positions at the
sth Monte Carlo run, and
,
represent respectively the target and the estimated ESV at the
sth Monte Carlo run.
M = 50 represents the total number of Monte Carlo runs.
4.1.1. Estimate the Unknown
For the first simulation, the measurement noise parameter is unknown to both the ASB and ESB algorithms, while the parameters of are chosen exactly their target value. Specifically, the initial is chosen the same as its target value as m/s, while the initial is set as s with the same initial offset for both the ASB and ESB algorithms.
The comparison of the
between the ASB and the ESB algorithm is shown in
Figure 2a,b from running the filter and the RTS smoother, respectively. It can be concluded that the tracking error associated with the ASB algorithm is noticeably smaller than that of the ESB algorithm.
Furthermore, the innovation of the transit time measurement
for the ESB and the ASB algorithm is compared in
Figure 3. Since
is a signed variable, it is compared from one randomly selected Monte Carlo simulation. From
Figure 3, it is obviously that
of the ASB algorithm converges rapidly then fluctuates around zero indicating a consistent filter, while that of the ESB algorithm has a large variation due to the improper setting of
.
From the comparison of
shown in
Figure 4, it can be seen that considering both the converging rate and steady state estimation error, the estimation performance of the ESV while implementing the ASB algorithm outperforms that of the ESB algorithm significantly.
Finally, the mean value of the estimated
while implementing the ASB algorithm from running 50 Monte Carlo simulations is shown in
Figure 5. The estimated
agrees well with its target value, while that in the ESB algorithm remains invariant as respected.
4.1.2. Estimate the Unknown
In the second scenario, the process noise parameter is unknown to both the ASB and ESB algorithms, while the parameter of measurement noise is chosen exactly their target value. Specifically, the initial is set as m/s for both the ASB and ESB algorithms with the same initial offset. To eliminate the possible impact of measurement noise, the hydrophone is kept clean without random noise contamination while generating the simulation data, and the initial is chosen as zero for consistency.
For comparing the performance of the ASB algorithm with those of the ESB algorithm, the
for both two algorithms are shown in
Figure 6a,b from running the filter and the RTS smoother, respectively. Clearly, using the ASB algorithm can significantly enhance the tracking accuracy over the ESB algorithm.
The comparison of the measurement innovations is shown in
Figure 7, in which similar features with the unknown
scenario can be observed.
The estimated
is shown in
Figure 8a,b from running the filter and the RTS smoother, respectively. Once again, the estimated ESV of ASB algorithm has a faster converging rate and smaller steady state error than that of ESB algorithm.
For completeness, the mean value of the estimated unknown
from 50 Monte Carlo simulations is shown in
Figure 9. Although the estimated
exist large variation at the beginning, the ASB algorithm can estimate its target value after a certain period of time.
4.2. Field Data
The field data used to demonstrate the efficiency of the ASB algorithm was collected from a surface boat equipped with a hydrophone. The boat also had access to the GPS, so the ground truth trajectory of the vehicle was available. A beacon was mounted at the sea floor with a surveyed location. Since the slant range between the vehicle and beacon could be computed from their known locations, the ground truth of every instant ESV could be computed from dividing the slant range by the corresponding transit time, obtained from the TOA measurement minus a known TOE.
While implementing the filter and the RTS smoother on the field data, the following initial quantities were chosen:
m/s,
m/s,
m/s and
s. First, to demonstrate the Kalman tuning issue of traditional ESB algorithms mentioned in
Section 2.3, a sensitivity study of
is carried out based on the field data. Shown in
Figure 10a is the averaged RMS of the horizontal distance error
with
varying from
to 1 m/s. Shown in
Figure 10b is the averaged RMS of the estimated ESV error
.
N and
K are the total number of fixed sampling interval and transit time measurements, respectively. From
Figure 10a,b, it is obviously that
has a great influence on the tracking performance. Improper
in the ESB algorithm will degenerate the estimation accuracy drastically.
The comparison of the estimated trajectories between the ASB and ESB is shown in
Figure 11a,b from running the filter and the RTS smoother, respectively. Similarly,
Figure 12a,b are the corresponding comparison of the horizontal distance error
for the filter and the RTS smoother, respectively. The component errors
and
for the ASB and ESB are shown in
Figure 13. From
Figure 11,
Figure 12 and
Figure 13, it can be concluded that using the ASB could greatly improve the tracking accuracy.
Shown in
Figure 14 is the estimated ESV from the filter and RTS smoother for these two algorithms, together with the target ESV. Similarly, the estimated ESV while implementing the ASB algorithm agrees better with its target value than the ESB algorithm.
5. Conclusions
To eliminate the range measurement error induced by imprecise knowledge of ESV, the so called “5-sv” single-beacon underwater tracking model was recently proposed. However, for applying the EKF based on the specific “5-sv” model, the tracking performance depended largely on the initial knowledge of process and measurement noise statistics, and filter divergence occurred frequently due to the huge variation of underwater acoustic propagation.
To overcome the limitation of EKF-based methods, this paper proposed an adaptive Kalman filter-based single-beacon underwater tracking algorithm basing upon the “5-sv” model. Through numerical examples of using simulated and field data, both the filter and RTS smoother results suggested that the tracking accuracy was significantly improved while implementing the proposed adaptive algorithm. The estimated process and measurement noise parameters also agreed well with their true value.