Next Article in Journal
Frame Registration for Motion Compensation in Imaging Photoplethysmography
Next Article in Special Issue
Optimization of Sparse Planar Arrays with Minimum Spacing and Geographic Constraints in Smart Ocean Applications
Previous Article in Journal
Background Light Rejection in SPAD-Based LiDAR Sensors by Adaptive Photon Coincidence Detection
Previous Article in Special Issue
Self-Organized Fast Routing Protocol for Radial Underwater Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity

Science and Technology on Underwater Vehicle Laboratory, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(12), 4339; https://doi.org/10.3390/s18124339
Submission received: 9 November 2018 / Revised: 3 December 2018 / Accepted: 4 December 2018 / Published: 8 December 2018

Abstract

:
In the single-beacon underwater tracking system, vehicles rely on slant range measurements from an acoustic beacon to bound errors accumulated by dead reckoning. Ranges are usually obtained based on a presumed known effective sound velocity (ESV). Since the ESV is difficult to determine accurately, traditional methods suffer from large positioning error. By treating the unknown ESV as a state variable, a novel single-beacon tracking model (the so called “5-sv” model) and an extended Kalman filter (EKF)-based solution method have been discussed to solve the problem of ESV estimation. However, due to the uncertainty of underwater acoustic propagation, the probabilistic characteristics of the ESV uncertainty and acoustic measurement noise are unknown and varying both with time and location. EKF, which runs with presupposed noise parameters, cannot describe the practical noise specifications. To overcome the divergence issue of EKF-based single-beacon tracking methods, this paper proposes an adaptive Kalman filter-based single-beacon tracking algorithm which employs the “5-sv” model as the baseline model. Through numerical examples using simulated and field data, both the filter and smoother results show that while implementing the proposed algorithm, the tracking accuracy can be significantly improved, and the estimated noise parameter agrees well with its true value.

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.

2. The “5-sv” Single-Beacon Tracking Model

A review of the “5-sv” single-beacon underwater tracking model, which is served as the baseline model employed in this paper, is offered here as a refresher. This review is also for acquainting the reader with the notation to be used throughout the article.

2.1. Kinematic Model

By treating the unknown ESV as a state variable, a single-beacon underwater tracking model—referred to as the “5-sv” (5 state variables) model—was proposed in [20,21,22]:
x ˙ ( t ) y ˙ ( t ) v ˙ c x ( t ) v ˙ c y ( t ) v ˙ e ( t ) = v w ( t ) cos φ ( t ) + v c x ( t ) v w ( t ) sin φ ( t ) + v c y ( t ) 0 0 0
where the state variables x ( t ) and y ( t ) represent the horizontal position of the vehicle; v c x ( t ) and v c y ( t ) are the two unknown ocean current components in the x and y directions, respectively; and v e ( t ) is the ESV between the beacon and the vehicle. When the ESV is multiplied by the acoustic signal transit time between two underwater point yields the geometric or slant range between them. In Equation (1), v w ( t ) and φ ( t ) are the vehicle speed and heading through water, which are measured from a speed sensor and an electronic compass, respectively.
While treating v w ( t ) and φ ( t ) as control parameters, and including the process noise w ( t ) , Equation (1) can be rewritten as
x ˙ ( t ) = F ( t ) x ( t ) + L ( t ) u ( t ) + w ( t )
where
F ( t ) = 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , L ( t ) = 1 0 0 1 0 0 0 0 0 0
and
x ( t ) = x ( t ) y ( t ) v c x ( t ) v c y ( t ) v e ( t ) , u ( t ) = v w x ( t ) v w y ( t )
in which v w x = v w cos φ and v w y = v w sin φ are the in-water speed components of the vehicle in the directions of x and y, respectively. The process noise vector w ( t ) R 5 × 1 accounts for the uncertainty of the kinematic model.
Assuming that the control term u ( t ) is constant within the discrete sampling interval Δ t , we obtain the corresponding discrete kinematic equation
x k = A k 1 x k 1 + B k 1 u k 1 + w k 1
where
A k = 1 0 Δ t 0 0 0 1 0 Δ t 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 , B k = Δ t 0 0 Δ t 0 0 0 0 0 0
and
x k = x k y k v c x , k v c y , k v e , k , u k = v w x , k v w y , k , w k = w x , k w y , k w c x , k w c y , k w e , k
in which the subscript k denotes the kth time step. The uncertainty w k associated with the discrete kinematic model is modeled as Gaussian white noise processes, consisting of five components in three groups: (1) w x , k and w y , k , position uncertainty in the x and y directions; (2) w c x , k and w c y , k , ocean current uncertainty in the x and y directions; and (3) w e , k , uncertainty associated with the ESV. The corresponding process noise covariance matrix is:
Q k = Δ t 2 ( σ w 2 cos 2 φ k + σ c 2 ) Δ t 2 σ w 2 cos φ k sin φ k 0 0 0 Δ t 2 σ w 2 cos φ k sin φ k Δ t 2 ( σ w 2 sin 2 φ k + σ c 2 ) 0 0 0 0 0 σ c 2 0 0 0 0 0 σ c 2 0 0 0 0 0 σ e , k 2
where σ w is the standard deviation of v w , k uncertainty, σ c is the standard deviation of w c x , k and w c y , k , and σ e , k is that of w e , k .

2.2. Measurement Model

In the application of the Kalman filter for single-beacon underwater tracking, the measurement equation is a nonlinear function of the state variables, with the general form
m k = h k ( x k ) + v k
in which h k ( · ) is a nonlinear function. It is assumed that the process and measurement noise w k and v k are mutually independent, zero-mean, Gaussian random processes, with covariance matrices Q k and R k , respectively.
In the present application, if the time of emission (TOE) T k e is known and the TOA T k a is measured, then the transit time T k t is given by
T k t = T k a T k e
By treating the transit time T k t as the measured quantity for m k , we have the nonlinear measurement equation
m k = h k ( x k ) + υ t , k
where h k can be expressed in terms of the state variables x k , y k and v e , k to be
h k = ( x k x b ) 2 + ( y k y b ) 2 + ( z k z b ) 2 v e , k
in which ( x k , y k , z k ) and ( x b , y b , z b ) are the positions of the vehicle and the beacon, respectively. Throughout this paper, the vertical position of an underwater vehicle z k assumed to be a known value obtained from a depth sensor. The transit time measurement noise υ t , k is modeled as a white Gaussian noise with variance R t , k = σ t , k 2 .
It is also assumed that the vehicle velocity relative to the ground v g is measured from a device, such as a Doppler velocity log (DVL). From the measured v g and the vehicle velocity through water v w , using v c = v g v w , one has an indirect measurement for the ocean current:
v c , k = v c x , k v c y , k
The corresponding measurement matrix H k associated with the ocean current measurements is simply
H k = 0 0 1 0 0 0 0 0 1 0

2.3. Comparison with Traditional Model

Traditional single-beacon tracking methods always treated the ESV as a known quantity. In reference to Equations (7) and (8), the corresponding discrete state vector and process uncertainty model are:
x k = x k y k v c x , k v c y , k
and
Q k = Δ t 2 ( σ w 2 cos 2 φ k + σ c 2 ) Δ t 2 σ w 2 cos φ k sin φ k 0 0 Δ t 2 σ w 2 cos φ k sin φ k Δ t 2 ( σ w 2 sin 2 φ k + σ c 2 ) 0 0 0 0 σ c 2 0 0 0 0 σ c 2
From Equations (4), (8), (15) and (16), it is clear that an additional parameter σ e is introduced and needed to be properly tuned while implementing EKF based on the “5-sv” model.
Furthermore, when the slant range r k computed from a transit time measurement is used in traditional models, the counterpart of Equation (12) is
h k = ( x k x b ) 2 + ( y k y b ) 2 + ( z k z b ) 2
Interested readers are referred to [20,21,22] for detailed comparison between the “5-sv” model and traditional models.
Compared from Equations (12) and (17), it can be seen that the degree of nonlinearity of the measurement equation in the “5-sv” model is much severer than that of traditional models. In essence, the “5-sv” model has converted the problem of estimating two unknowns ( x k and y k ) from one measurement ( r k ) into the problem of estimating three unknowns ( x k , y k and v e , k ) from one measurement ( T k t ). Difficulties of the Kalman tuning process while using the “5-sv” model also appears in that the setting of ESV uncertainty σ e interacts with the initial setting of position uncertainties. EKF based on the “5-sv” model is much harder to tune, and more likely to diverge than traditional models if the initial settings are improper. Furthermore, the physics of acoustic wave propagation and the underwater acoustic environment creates some unique differences, including longer transit times and severe refraction through the stratified ocean. Thus, uncertainty associated with the acoustic propagation velocity are very difficult to deal with because they cannot be measured and are generally varying with both space and time. Even if an approximate σ e is obtained from experience or data post-processing, divergence issue of EKF based on the “5-sv” model would still occur frequently due to the violation of known noise statistics assumptions. In Section 4.2, a sensitivity study of σ e based on filed data will be done to indicate the effect of σ e on the estimator performance.

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 x ^ k and the corresponding covariance matrix P k are
x ^ k = A k 1 x ^ k 1 + + B k 1 u k 1
and
P k = A k 1 P k 1 + A k 1 T + Q ^ k 1
where Q ^ k 1 is the estimated process noise covariance matrix at the ( k 1 ) th epoch.
Throughout this paper, a variable with a hat, such as x ^ , 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 x ^ k + and the corresponding covariance matrix P k + are
x ^ k + = x ^ k + K k ( m k h ( x ^ k ) )
and
P k + = ( I K k H k ) P k
in which
K k = P k H k T ( H k P k H k T + R ^ k ) 1
is the Kalman gain matrix and R ^ k is the estimated measurement noise covariance matrix. For the single-beacon underwater tracking algorithm discussed herein, H k is the Jaccbian of nonlinear measurement equation h ( ) at the a priori estimate x ^ k
H k ( x ^ k ) = x ^ k x b v ^ e , k r ^ k y ^ k y b v ^ e , k r ^ k 0 0 r ^ k ( v ^ e , k ) 2
in which r ^ k is computed by
r ^ k = ( x ^ k x b ) 2 + ( y ^ k y b ) 2 + ( z ^ k z b ) 2
Noting that z ^ k is available from a depth sensor. The measurement innovation ( m k h ( x ^ k ) ) in Equation (20) is denoted by e ¯ k in this paper.

3.3. Estimation of Unknown R k

From literature [23], the main purpose of the Windows-based adaptive estimation is to estimate the R k 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:
J k = j = k W + 1 k p ( z j | C ¯ v j )
where C ¯ v j = E ( e ¯ j e ¯ j T ) is the covariance matrix of innovation vector at time t j .
Assuming e ¯ j satisfies zero-mean Gauss distribution, Equation (25) can be written as
J k = j = k W + 1 k N ( z j : h ( x ^ k ) , C ¯ v j )
where N ( x : μ , Q ) means the vector x satisfies Gaussian distribution with mean μ and covariance matrix Q . The probability density function is
N ( x : μ , Q ) = 1 | 2 π Q | exp ( 1 2 ( x μ ) T Q 1 ( x μ ) )
The logarithmic form of Equation (26) is taken as
ln J k = j = k W + 1 k 1 2 ( ln | 2 π × C ¯ v j | + e ¯ j T C ¯ v j 1 e ¯ j )
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
C ¯ v j = C ¯ v k , ( k W + 1 ) j k
After multiplying Equation (28) by 2 and neglecting the constant term, the ML criterion of maximizing J within the windows width W becomes the minimization problem which is described as
C ¯ ^ v k = arg min J 1 , k J 1 , k = j = k W + 1 k ln | C ¯ v k | + j = k W + 1 k e ¯ j T C ¯ v k 1 e ¯ j
The minimum value is calculated by the partial of J 1 , k with respect to C ¯ v k
J 1 , k C ¯ v k = W tr { C ¯ v k 1 } j = k W + 1 k e ¯ j T C ¯ v k 1 C ¯ v k 1 e ¯ j
To obtain the above formula, the following two relations from matrix differential calculus have been used
ln | A | x = 1 | A | | A | x = tr { A 1 | A | x }
A 1 x = A 1 | A | x A 1
Based on the relations from vector inner product x T x = tr ( x x T ) and the trace property of matrix, Equation (31) can be written as
J 1 , k C ¯ v k = tr { C ¯ v k 1 ( W C ¯ v k j = k W + 1 k e ¯ j e ¯ j T ) C ¯ v k 1 }
Let
J 1 , k C ¯ v k = 0
then the estimated innovation covariance matrix is obtained as
C ¯ ^ v k = 1 W j = k W + 1 k e ¯ j e ¯ j T
On the other hand, the theoretical innovation covariance matrix can be calculated through the definition of innovation, and have
C ¯ v k = H k P k H k T + R k
Let the estimated value of the innovation covariance matrix equals to the theoretical value, the estimated R k can then be obtained
R ^ k = 1 W j = k W + 1 k e ¯ j e ¯ j T H k P k H k T
This estimated R k will be used for calculating the Kalman gain at epoch k.

3.4. Estimation of Unknown Q k

To estimate the unknown Q k , rewritten Equation (22) as
K k = P k H k T C ¯ ^ v k 1
where C ¯ ^ v k is assumed to predetermined from Equation (36). Since C ¯ ^ v k and P k are symmetric matrices, Equation (39) can be written as
C ¯ ^ v k K k T = H k P k
Plugging the left side of Equation (40) into Equation (21) yields
P k = K k C ¯ ^ v k K k T + P k +
From Equations (41) and (19), the estimated Q ^ k 1 can be get
Q ^ k 1 = K k C ¯ ^ v k K k T A k 1 P k 1 + A k 1 T + P k +
To calculate Q ^ k , the following hypothesis is made:
Hypothesis 2.
The varying period of Q k is much larger than the sampling period Δ t , thus we have
Q k Q k 1
Based on the Hypothesis 2, Equations (36) and (42), the estimated Q k at time t k can be obtained
Q ^ k = 1 W j = k W + 1 k K k e ¯ j e ¯ j T K k T A k 1 P k 1 + A k 1 T + P k +
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: x ^ k 1 + , P k 1 + , { e ¯ k W + 1 , e ¯ k W + 2 , e ¯ k 1 }
Prediction:
1:
x ^ k = A k 1 x ^ k 1 + + B k 1 u k 1
2:
P k = A k 1 P k 1 + A k 1 T + Q ^ k 1

Correction:
1:
r ^ k = ( x ^ k x b ) 2 + ( y ^ k y b ) 2 + ( z ^ k z b ) 2
2:
H k = [ x ^ k x b v ^ e , k r ^ k y ^ k y b v ^ e , k r ^ k 0 0 r ^ k ( v ^ e , k ) 2 ]
3:
e ¯ k = m k h ( x ^ k )
4:
R ^ k = 1 W j = k W + 1 k e ¯ j e ¯ j T H k P k H k T
5:
K k = P k H k T ( H k P k H k T + R ^ k ) 1
6:
x ^ k + = x ^ k + K k e ¯ k
7:
P k + = ( I K k H k ) P k
8:
Q ^ k = 1 W j = k W + 1 k K k e ¯ j e ¯ j T K k T A k 1 P k 1 + A k 1 T + P k +

Return: x ^ k + , P k + , e ¯ k

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 0.3 m/s. The process noise parameter σ e is chosen as 0.1 m/s. σ c and σ w are both set to 0.01 m/s. For clarity, the estimation performance of the unknown R k and Q k are studied in two distinct scenarios, respectively. The Q k is assumed to be precisely available during the estimation of unknown R k , 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) v e , 0 = 1540 m/s for the initial ESV; (3) 0.35 m/s for the initial ocean current in both x and y directions; and (4) σ w = 0.01 m/s; (5) σ c = 0.01 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 Δ H and the estimation error of the ESV Δ v e are chosen as the evaluation indexes of the Monte Carlo simulations, which are defined respectively as:
RMS Δ H 1 M s = 1 M ( x k s x ^ k s ) 2 + ( y k s y ^ k s ) 2 RMS Δ v e 1 M s = 1 M ( v e , k s v ^ e , k s ) 2
in which ( x k s , y k s ) and ( x ^ k s , y ^ k s ) represent respectively the target and the estimated positions at the sth Monte Carlo run, and v e , k s , v ^ e , k s 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 R k

For the first simulation, the measurement noise parameter is unknown to both the ASB and ESB algorithms, while the parameters of Q k are chosen exactly their target value. Specifically, the initial σ e is chosen the same as its target value as 0.1 m/s, while the initial σ t is set as 0.05 s with the same initial offset for both the ASB and ESB algorithms.
The comparison of the RMS Δ H 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 e ¯ t for the ESB and the ASB algorithm is compared in Figure 3. Since e ¯ t is a signed variable, it is compared from one randomly selected Monte Carlo simulation. From Figure 3, it is obviously that e ¯ t 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 σ t .
From the comparison of RMS Δ v e 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 σ t while implementing the ASB algorithm from running 50 Monte Carlo simulations is shown in Figure 5. The estimated σ t agrees well with its target value, while that in the ESB algorithm remains invariant as respected.

4.1.2. Estimate the Unknown Q k

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 σ e is set as 0.5 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 σ t is chosen as zero for consistency.
For comparing the performance of the ASB algorithm with those of the ESB algorithm, the RMS Δ H 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 R k scenario can be observed.
The estimated RMS Δ v e 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 σ e from 50 Monte Carlo simulations is shown in Figure 9. Although the estimated σ e 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: v e , 0 = 1560 m/s, σ c = 0.01 m/s, σ w = 0.01 m/s and σ t = 0.0001 s. First, to demonstrate the Kalman tuning issue of traditional ESB algorithms mentioned in Section 2.3, a sensitivity study of σ e is carried out based on the field data. Shown in Figure 10a is the averaged RMS of the horizontal distance error ARMS Δ H = i = 1 N ( ( x i x ^ i ) 2 + ( y i y ^ i ) 2 ) with σ e varying from 0.01 to 1 m/s. Shown in Figure 10b is the averaged RMS of the estimated ESV error ARMS Δ v e = i = 1 K ( v e , i v ^ e , i ) 2 . N and K are the total number of fixed sampling interval and transit time measurements, respectively. From Figure 10a,b, it is obviously that σ e has a great influence on the tracking performance. Improper σ e 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 Δ H for the filter and the RTS smoother, respectively. The component errors Δ x and Δ y 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.

Author Contributions

Z.-C.D. and X.Y. conceived of and designed the method and performed the numerical study. H.-D.Q. and Z.-B.Z. participated in the experiment and analyzed the data. X.Y. wrote the paper. H.-D.Q. and Z.-B.Z. guided and supervised the overall process.

Funding

This work was funded in part by the National Natural Science Foundation of China under Grant U1713205 and 51809062, and in part by the Stable Supporting Fund of Science and Technology on Underwater Vehicle Laboratory under Grant SXJQR2018WDKT07.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Q.; Zhang, J.; Chemori, A.; Xiang, X. Virtual Submerged Floating Operational System for Robotic Manipulation. Complexity 2018, 2018, 9528313. [Google Scholar] [CrossRef]
  2. Xiang, X.; Yu, C.; Lapierre, L.; Zhang, J.; Zhang, Q. Survey on FuzzyLogicBased Guidance and Control of Marine Surface Vehicles and Underwater Vehicles. Int. J. Fuzzy Syst. 2018, 20, 572–586. [Google Scholar] [CrossRef]
  3. Xiang, X.; Liu, C.; Su, H.; Zhang, Q. On decentralized adaptive full-order sliding mode control of multiple UAVs. ISA Trans. 2017, 71, 196–205. [Google Scholar] [CrossRef] [PubMed]
  4. Xiang, X.; Yu, C.; Zhang, Q. On intelligent risk analysis and critical decision of underwater robotic vehicle. Ocean Eng. 2017, 140, 453–465. [Google Scholar] [CrossRef]
  5. Roper, D.T.; Phillips, A.B.; Harris, C.A.; Salavasidis, G.; Pebody, M.; Templeton, R.; Amma, S.V.S.; Smart, M.; McPhail, S. Autosub long range 1500: An ultra-endurance AUV with 6000 Km range. In Proceedings of the OCEANS 2017, Aberdeen, UK, 19–22 June 2017; pp. 1–5. [Google Scholar] [CrossRef]
  6. Carreras, M.; Hernandez, J.D.; Vidal, E.; Palomeras, N.; Ribas, D.; Ridao, P. Sparus II AUV—A Hovering Vehicle for Seabed Inspection. IEEE J. Ocean. Eng. 2018, 43, 344–355. [Google Scholar] [CrossRef]
  7. Paull, L.; Saeedi, S.; Seto, M.; Li, H. AUV Navigation and Localization: A Review. IEEE J. Ocean. Eng. 2014, 39, 131–149. [Google Scholar] [CrossRef]
  8. Ferreira, B.; Matos, A.; Cruz, N. Single beacon navigation: Localization and control of the MARES AUV. In Proceedings of the OCEANS 2010 MTS/IEEE SEATTLE, Seattle, WA, USA, 20–23 September 2010; IEEE: New York, NY, USA, 2010; pp. 1–9. [Google Scholar]
  9. Wang, S.; Chen, L.; Hu, H.; Gu, D. Single beacon based localization of AUVs using moving horizon estimation. In Proceedings of the 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3–7 November 2013; IEEE: New York, NY, USA, 2013; pp. 885–890. [Google Scholar]
  10. Harris, Z.J.; Whitcomb, L.L. Preliminary study of cooperative navigation of underwater vehicles without a DVL utilizing range and range-rate observations. In Proceedings of the IEEE International Conference on Robotics and Automation, Stockholm, Sweden, 16–21 May 2016; pp. 2618–2624. [Google Scholar]
  11. Claus, B.; Kepper, J.H.; Suman, S.; Kinsey, J.C. Closed loop one way travel time navigation using low grade odometry for autonomous underwater vehicles. J. Field Robot. 2017, 35, 421–434. [Google Scholar] [CrossRef]
  12. Batista, P. GES Long Baseline Navigation With Unknown Sound Velocity and Discrete-Time Range Measurements. IEEE Trans. Control Syst. Technol. 2015, 23, 219–230. [Google Scholar] [CrossRef]
  13. Thomson, D.J.M.; Dosso, S.E.; Barclay, D.R. Modeling AUV Localization Error in a Long Baseline Acoustic Positioning System. IEEE J. Ocean. Eng. 2018, 43, 955–968. [Google Scholar] [CrossRef]
  14. Chen, Y.; Zheng, D.; Miller, P.A.; Farrell, J.A. Underwater Inertial Navigation With Long Baseline Transceivers: A Near-Real-Time Approach. IEEE Trans. Control Syst. Technol. 2016, 24, 240–251. [Google Scholar] [CrossRef]
  15. Miller, P.A.; Farrell, J.A.; Zhao, Y.; Djapic, V. Autonomous Underwater Vehicle Navigation. IEEE J. Ocean. Eng. 2010, 35, 663–678. [Google Scholar] [CrossRef] [Green Version]
  16. Larsen, M.B. Synthetic long baseline navigation of underwater vehicles. In Proceedings of the Oceans, Providence, RI, USA, 11–14 September 2000; Volume 3, pp. 2043–2050. [Google Scholar]
  17. Casey, T.; Guimond, B.; Hu, J. Underwater vehicle positioning based on time of arrival measurements from a single beacon. In Proceedings of the OCEANS 2007, Vancouver, BC, Canada, 29 September–4 October 2007; IEEE: New York, NY, USA, 2007; pp. 1–8. [Google Scholar]
  18. Fallon, M.F.; Papadopoulos, G.; Leonard, J.J.; Patrikalakis, N.M. Cooperative AUV navigation using a single maneuvering surface craft. Int. J. Robot. Res. 2010, 29, 1461–1474. [Google Scholar] [CrossRef]
  19. Walls, J.M.; Eustice, R.M. An origin state method for communication constrained cooperative localization with robustness to packet loss. Int. J. Robot. Res. 2014, 33, 1191–1208. [Google Scholar] [CrossRef] [Green Version]
  20. Zhu, Z.; Hu, S.L.J. Model and Algorithm Improvement on Single Beacon Underwater Tracking. IEEE J. Ocean. Eng. 2017, 43, 1143–1160. [Google Scholar] [CrossRef]
  21. Zhu, Z.; Hu, S.L.J.; Li, H. Kalman-based underwater tracking with unknown effective sound velocity. In Proceedings of the Oceans, Monterey, CA, USA, 19–23 September 2016; pp. 1–9. [Google Scholar]
  22. Zhu, Z.; Hu, S.L.J.; Li, H. Effect on Kalman based underwater tracking due to ocean current uncertainty. In Proceedings of the Autonomous Underwater Vehicles, Tokyo, Japan, 6–9 November 2016; pp. 131–137. [Google Scholar]
  23. Mohamed, A.H.; Schwarz, K.P. Adaptive Kalman Filtering for INS/GPS. J. Geod. 1999, 73, 193–203. [Google Scholar] [CrossRef]
  24. Qin, Y.; Peng, H.; Xi, Y.; Xie, W.; Sun, Y.; Chen, X. An adaptive modeling and asset allocation approach to financial markets based on discrete microstructure model. Appl. Soft Comput. 2016, 43, 390–405. [Google Scholar] [CrossRef]
  25. Zhang, W.; Shi, W.; Ma, Z. Adaptive unscented Kalman filter based state of energy and power capability estimation approach for lithium-ion battery. J. Power Sources 2015, 289, 50–62. [Google Scholar] [CrossRef]
  26. Simon, D. Optimal State Estimation: Kalman, H Infinity, and Nonlinear Approaches; John Wiley & Sons: Hoboken, NJ, USA, 2006; pp. 123–129. [Google Scholar]
  27. Sarkka, S. Bayesian Filtering and Smoothing; Institute of Mathematical Statistics Textbooks, Cambridge University Press: New York, NY, USA, 2013. [Google Scholar]
  28. Gelb, A. Applied Optimal Estimation; MIT Press: Cambridge, MA, USA; London, UK, 1974; pp. 107–118. [Google Scholar]
Figure 1. Simulated Trajectory.
Figure 1. Simulated Trajectory.
Sensors 18 04339 g001
Figure 2. RMS of Horizontal distance error comparison between the ESB and the ASB algorithms using simulation data with unknown R k : (a) Filter; (b) RTS smoother.
Figure 2. RMS of Horizontal distance error comparison between the ESB and the ASB algorithms using simulation data with unknown R k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g002
Figure 3. Comparison of transit time measurement innovation between the ESB and the ASB algorithms using simulation data with unknown R k .
Figure 3. Comparison of transit time measurement innovation between the ESB and the ASB algorithms using simulation data with unknown R k .
Sensors 18 04339 g003
Figure 4. RMS of ESV estimated error by the ESB and the ASB algorithms using simulation data with unknown R k : (a) Filter; (b) RTS smoother.
Figure 4. RMS of ESV estimated error by the ESB and the ASB algorithms using simulation data with unknown R k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g004
Figure 5. Estimated measurement noise of the ESB and the ASB algorithms using simulation data with unknown R k .
Figure 5. Estimated measurement noise of the ESB and the ASB algorithms using simulation data with unknown R k .
Sensors 18 04339 g005
Figure 6. RMS of horizontal distance error comparison between the ESB and the ASB algorithms using simulation data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 6. RMS of horizontal distance error comparison between the ESB and the ASB algorithms using simulation data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g006
Figure 7. Comparison of transit time measurement innovation between the ESB and the ASB algorithms using simulation data with unknown Q k .
Figure 7. Comparison of transit time measurement innovation between the ESB and the ASB algorithms using simulation data with unknown Q k .
Sensors 18 04339 g007
Figure 8. RMS of the ESV estimation error comparison between the ESB and the ASB algorithms using simulation data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 8. RMS of the ESV estimation error comparison between the ESB and the ASB algorithms using simulation data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g008
Figure 9. Estimated process noise comparison between the ESB and the ASB algorithms using simulation data with unknown Q k .
Figure 9. Estimated process noise comparison between the ESB and the ASB algorithms using simulation data with unknown Q k .
Sensors 18 04339 g009
Figure 10. Estimation performance versus σ e from EKF using field data: (a) RMS Δ r ; (b) RMS Δ v e .
Figure 10. Estimation performance versus σ e from EKF using field data: (a) RMS Δ r ; (b) RMS Δ v e .
Sensors 18 04339 g010
Figure 11. Planar position estimates comparison between the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 11. Planar position estimates comparison between the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g011
Figure 12. Horizontal distance error comparison between the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 12. Horizontal distance error comparison between the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g012
Figure 13. Position component errors of the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 13. Position component errors of the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g013
Figure 14. Estimated ESV by the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Figure 14. Estimated ESV by the ESB and the ASB algorithms using field data with unknown Q k : (a) Filter; (b) RTS smoother.
Sensors 18 04339 g014
Table 1. Sampling frequency and noise level of simulated measurements.
Table 1. Sampling frequency and noise level of simulated measurements.
SensorMeasurementSimulated Sampling
Frequency (Hz)
Noise ( σ )
DVL v g x , v g y 1 0.002 m/s
electronic compass φ 10.2
speed sensor v w 1 0.001 m/s
hydrophoneTOA 0.001 s

Share and Cite

MDPI and ACS Style

Deng, Z.-C.; Yu, X.; Qin, H.-D.; Zhu, Z.-B. Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity. Sensors 2018, 18, 4339. https://doi.org/10.3390/s18124339

AMA Style

Deng Z-C, Yu X, Qin H-D, Zhu Z-B. Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity. Sensors. 2018; 18(12):4339. https://doi.org/10.3390/s18124339

Chicago/Turabian Style

Deng, Zhong-Chao, Xiang Yu, Hong-De Qin, and Zhong-Ben Zhu. 2018. "Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity" Sensors 18, no. 12: 4339. https://doi.org/10.3390/s18124339

APA Style

Deng, Z. -C., Yu, X., Qin, H. -D., & Zhu, Z. -B. (2018). Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity. Sensors, 18(12), 4339. https://doi.org/10.3390/s18124339

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