Next Article in Journal
G-Quadruplex DNAzyme Molecular Beacon for Amplified Colorimetric Biosensing of Pseudostellaria heterophylla
Next Article in Special Issue
An Aerial-Ground Robotic System for Navigation and Obstacle Mapping in Large Outdoor Areas
Previous Article in Journal
A Flexible Sensor Technology for the Distributed Measurement of Interaction Pressure
Previous Article in Special Issue
About Non-Line-Of-Sight Satellite Detection and Exclusion in a 3D Map-Aided Localization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Scheme for DVL-Aided SINS In-Motion Alignment Using UKF Techniques

1
College of Mechatronics and Automation, National University of Defense Technology, Changsha 410073, China
2
School of Surveying and Geospatial Engineering, University of New South Wales, Sydney, NSW 2052, Australia
*
Author to whom correspondence should be addressed.
Sensors 2013, 13(1), 1046-1063; https://doi.org/10.3390/s130101046
Submission received: 7 November 2012 / Revised: 25 December 2012 / Accepted: 5 January 2013 / Published: 15 January 2013
(This article belongs to the Special Issue New Trends towards Automatic Vehicle Control and Perception Systems)

Abstract

:
In-motion alignment of Strapdown Inertial Navigation Systems (SINS) without any geodetic-frame observations is one of the toughest challenges for Autonomous Underwater Vehicles (AUV). This paper presents a novel scheme for Doppler Velocity Log (DVL) aided SINS alignment using Unscented Kalman Filter (UKF) which allows large initial misalignments. With the proposed mechanism, a nonlinear SINS error model is presented and the measurement model is derived under the assumption that large misalignments may exist. Since a priori knowledge of the measurement noise covariance is of great importance to robustness of the UKF, the covariance-matching methods widely used in the Adaptive KF (AKF) are extended for use in Adaptive UKF (AUKF). Experimental results show that the proposed DVL-aided alignment model is effective with any initial heading errors. The performances of the adaptive filtering methods are evaluated with regards to their parameter estimation stability. Furthermore, it is clearly shown that the measurement noise covariance can be estimated reliably by the adaptive UKF methods and hence improve the performance of the alignment.

1. Introduction

With the development of high-frequency, multi-beam Doppler sonar, which can provide bottom velocity measurements with a precision of 0.3% or less with a update rate of up to 5Hz, a wide variety of Doppler-based navigation techniques have been developed [1]. A typical navigation method for AUVs with minimal sonar use is based on a high quality SINS combined with occasional use of a Doppler Velocity Log [2]. This is an open-loop system so that the initial alignment is of great importance to subsequent navigation operation. In the case of an AUV, the initial alignment is more difficult. That is because the external aiding sensors such as GPS which provide geodetic-frame observations are unavailable for most of the time [3]. Therefore, it is an essential task to achieve an accurate alignment using DVL aiding within a very short period of time.
Due to the random wave motions as well as the dynamics of the vehicle, it is difficult to estimate the attitude to the accuracy of within a few degrees in a short period with the existing coarse alignment methods [4]. Therefore, in-motion SINS alignment with large initial misalignments is always a challenge and thus needs to be considered. The difficulty for the DVL-aided alignment is that DVL provides the velocity measurements in the Doppler instrumental frame and hence could not be used as the measurement for alignment directly. There are mainly two alignment schemes to solve this problem for small misalignments [5,6].The first method is to establish the INS error dynamics in the body frame, so the velocity of the Doppler can be used as the measurements directly [5]. However, the INS error equations will include the unavoidable attitude error. Therefore, whether this method could be used for large misalignments problem needs to be analyzed. The other method is to establish the INS error dynamics in the navigation frame [6] and this scheme is shown in Figure 1. The velocity of Doppler is transformed to the navigation frame by the attitude matrix obtained from INS. Then it can be used as the measurements for the Kalman filter. The problem here is that the attitude errors are usually very large before the alignment is finished. It will cause a large measurement error that may lead to the divergence of the Kalman filter [5]. But the authors in [5] failed to notice that the measurement error could be compensated in the measurement model. By employing this scheme, a new alignment model which allows large misalignments is proposed in this study. From experimental results, the proposed alignment model is shown to be effective with any initial attitude.
The other main effort to deal with the large misalignments problem is based on such nonlinear filtering methods as the so-called extended Kalman filter (EKF), unscented Kalman filter (UKF), and particle filter (PF). Among these nonlinear filtering methods, the UKF is wildly used due to its elimination of the cumbersome derivation and low computational complexity [7]. The UKF is based on the unscented transformation (UT) which is founded on the intuition that an approximation of a probability distribution is easier than that of an arbitrary nonlinear function [8]. It is able to capture the true mean and covariance of the Gaussian Random Variable to the 3rd order accuracy [9] and hence proved to be superior to EKF [10]. The most obvious shortcoming of EKF is that it requires the computation of Jocobian matrices and linear approximations of the EKF can be very inaccurate in some scenarios, leading to filter instability [9,11]. Similar to the KF [1214] and EKF, covariance parameters also play an important role in the filtering performance. Therefore, several adaptive filtering techniques are developed to improve the performance of the UKF [15]. The covariance-matching method which is based on the concept of making the elements of the actual innovation covariance matrix consistent with their theoretical values has been shown to be one of the most promising techniques for KF [16]. This method is extended to UKF for its application in nonlinear systems [17]. However, it could not be guaranteed that the resulting measurement noise covariance is positive definite. Inspired by the work in [18], another covariance-matching adaptive filtering method based on the filter residual sequence is presented. This study first evaluates the impacts of adaptive filtering methods on the parameter estimation stability with different window sizes and different initial measurement noise covariance matrix. Experimental results demonstrate that the measurement noise covariance can be estimated reliably by the adaptive filtering methods. The performance of the adaptive filtering methods are also compared and analyzed.
The rest of this paper is organized as follows: Section 2 is devoted to the presentation of the nonlinear DVL-aided IMU alignment model which can tolerate large misalignments. Section 3 presents the mathematical formulas of the UKF and adaptive UKF techniques. In Section 4, the performance of the proposed algorithms are evaluated and compared with real experimental data. The conclusions are drawn in Section 5.

2. Alignment Model

2.1. INS Error Dynamics Model

The nonlinear SINS error model proposed in [19,20] will be used in this paper. The nominal error state includes the velocity error (δvc), the attitude error (ψ), the accelerometer bias (∇b) and the gyro bias (εb). The formulas of this model are given as follows:
δ v ˙ c = ( I C p c ) C b p f b ( 2 Ω i e c + Ω e c c ) δ v c δ Ω e c c v c + C b c b + δ g c
ψ ˙ = ( I C c p ) ω i c c C b p ε b
˙ b = 0
ε ˙ b = 0
where the superscript c denotes the local level frame at the computed position and its orientation is north-up-east (NUE). p, i and b are the platform frame, the inertial frame and the body frame respectively. C is the direction cosine matrix. f b is the accelerometer measurement. δgc is the gravity error caused by the position error and it can be ignored during the process of the alignment. Ω i e c is the skew matrix of ω i e c which is the angle rate of the Earth frame relative to the inertial frame. ω i e c is given by:
ω i e c = [ ω i e cos l ω i e sin l 0 ]
where: l is the geographic latitude; ωie is the earth's rotation angular velocity. Ω e c c is the skew matrix of ω e c c which is the angle rate of the navigation frame relative to the Earth frame. ω e c c is given by:
ω e c c = [ V E R n V E tan l R n V N R m ]
where Rn and Rm are the transverse and meridian radius of the curvature respectively; the subscripts E and N denote the east and north components in the c frame. ω i c c is the angle rate of the navigation frame relative to the inertial frame. It can be obtained by:
ω i c c = ω i e c + ω e c c
Define:
[ s x c x s y c y s z c z ] = [ sin ( ψ x ) cos ( ψ x ) sin ( ψ y ) cos ( ψ y ) sin ( ψ z ) cos ( ψ z ) ]
C p c is given by:
C p c = [ c z c y s x s y c x s z c y c x s y + s x s z c y s z c x c z s x c z c z s y s x c y + c x s z s y c x c y s x s z s y ]
In the error model presented above, all the three misalignment angles are assumed to be large. For real time applications, it is often the case that there are a large uncertainty in heading angle and low uncertainties in leveling angles [20].

2.2. Measurement Model

The velocity of Doppler in the local level frame v d c can be described as follows:
v d c = C b c C d b v d
where C d b is a constant matrix which translates the velocity of Doppler vd from the Doppler instrumental frame to the vehicle body frame. It needs to be calibrated before a mission is conducted [21]. Supposing the error of v d c is mainly caused by the platform misalignments, the DVL measurement in error v d c can be described as follows:
v ˜ d c = C b c C d b v d + δ C b c C d b v d = v d c + δ C b c C d b v d
where the perturbation of the attitude matrix δ C b c is given by [22]:
δ C b c = ( I C c p ) C b c
Inserting Equation (12) into Equation (11), it yields:
v ˜ d c = v d c + ( I C c p ) C b c C d b v d
Differentiating the velocity of INS and DVL, the measurement model is given below:
v ˜ INS c v ˜ d c = v INS c + δ v INS c v d c ( I C c p ) C b c C d b v d = δ v INS n ( I C c p ) C b c C d b v d

3. UKF Techniques

3.1. UKF in Additive Noise Case

The considered nonlinear discrete-time system with additive noise is presented as follows [23]:
x k = f ( x k 1 ) + w k 1
z k = h ( x k ) + υ k
where xkRn is the state vector; zkRm is the measurement vector; f(·) and h(·) are nonlinear functions; wk and υk are the uncorrelated zero-mean Gaussian white sequences and their conversances are:
E [ w k w j T ] = δ k j Q k E [ υ k υ j T ] = δ k j R k E [ w k υ j T ] = 0
With a view of reducing the computational burden, the non-augmented UKF is widely used in such additive noise cases [7,2426]. The resemblance between the UKF and the KF is that the implementations of the two algorithms all consist of the prediction of the state mean and covariance, and then the update with the measurements [27]. The implementation of the non-augmented UKF algorithm is given as follows [25,26]:
  • Initialization:
    x ˆ 0 = E [ x 0 ] , P 0 = E [ ( x ˆ 0 x 0 ) ( x ˆ 0 x 0 ) T ]
  • Time-updating:
    χ k 1 = [ x ˆ k 1 [ x ˆ k 1 ] L + γ P k 1 [ x ˆ k 1 ] L γ P k 1 ]
    χ i , k k 1 = f ( χ i , k 1 )
    x ˆ k k 1 = i = 0 2 L W i ( m ) χ i , k k 1
    P k k 1 = i = 0 2 L W i ( c ) ( χ i , k k 1 x ˆ k k 1 ) ( χ i , k k 1 x ˆ k k 1 ) T + Q k 1
    χ k k 1 = [ x ˆ k k 1 [ x ˆ k k 1 ] L + γ P k k 1 [ x ˆ k k 1 ] L γ P k k 1 ]
    η i , k k 1 = h ( χ i , k k 1 )
    z ˆ k k 1 = i = 0 2 L W i ( m ) η i , k k 1
  • Measurement-updating:
    P z ˆ k z ˆ k = i = 0 2 L W i ( c ) ( η i , k k 1 z ˆ k k 1 ) ( η i , k k 1 z ˆ k k 1 ) T + R k
    P x ˆ k z ˆ k = i = 0 2 L W i ( c ) ( χ i , k k 1 x ˆ k k 1 ) ( η i , k k 1 z ˆ k k 1 ) T
    K k = P x ˆ k z ˆ k P z ˆ k z ˆ k 1
    x ˆ k = x ˆ k k 1 + K k ( z k z ˆ k k 1 )
    P k = P k k 1 K k P z ˆ k z ˆ k K k T
The parameters for calculating the sigma-points are given as follows:
{ λ = α 2 ( L + κ ) L γ = L + λ W 0 ( m ) = λ / ( L + λ ) W 0 ( c ) = λ / ( L + λ ) + ( 1 α 2 + β ) W i ( m ) = W i ( c ) = 1 / [ 2 ( L + λ ) ] , ( i = 1 , 2 , , 2 L )
where W(m) and W(c) represent the mean weight and covariance weight, respectively [15]; L is the dimension of the state vector; P k k 1 is the i th column of the matrix square root of Pkk−1; α determines the spread of the sigma points around and is usually set to a small positive value (e.g., 1e-3); κ is a secondary scaling parameter which is usually set to 0; β is used to incorporate prior knowledge of the distribution of x (for Gaussian distributions, β = 2 is optimal) [28].

3.2. Innovation-Based Adaptive UKF

The innovation sequence dk which is the difference between the measurement zk and its predicted value kk−1 is given as follows:
d k = z k z ˆ k k 1
where kk−1 is obtained from Equation (25). By matching the covariance matrix of the innovation sequence with its theoretical form, Pkk can be calculated as follows [17]:
P z ˆ k z ˆ k = 1 N j = k N + 1 k d j d j T
N is the length of the sliding sampling window. By replacing the Equation (26) with the Equation (33), the Adaptive UKF (AUKF) algorithm is obtained. This technique is similar to the traditional innovation-based adaptive Kalman filter (AKF) [16,29]. By comparing the Equation (33) with (26), the estimation of the measurement noise covariance can be obtained as follows:
R k = 1 N j = k N + 1 k d j d j T i = 0 2 L W i ( c ) ( η i , k k 1 z ˆ k k 1 ) ( η i , k k 1 z ˆ k k 1 ) T
As can be seen from Equation (34), the estimated measurement noise covariance is the difference between two positive definite matrices, it cannot be guaranteed that the resulting matrix Rk is positive definite. This trouble can easily cause the failure of the filter in real time applications.

3.3. Residual-Based Adaptive UKF

The residual sequence could also be used with the purpose of obtaining a realistic estimator of the measurement noise covariance. The residual sequence εk is defined as follows:
ε k = z k z ˆ k
where k can be obtained by the estimated values k (not the predicted values k|k−1) by an extra UT:
χ k = [ x ˆ k [ x ˆ k ] L + γ P k [ x ˆ k ] L γ P k ]
η i , k = h ( χ i , k )
z ˆ k = i = 0 2 L W i ( m ) η i , k
by extending concept of the residual-based AKF [18] to UKF, a recursive estimator of Rk can be obtained as:
R k = 1 N j = k N + 1 k ε j ε j T + i = 0 2 L W i ( c ) ( η i , k z ˆ k ) ( η i , k z ˆ k ) T
It can be used in the computation of epoch k + 1. Compared to Equation (34), as Equation (39) is the sum of two positive definite matrices, the estimated covariance matrix is always positive definite. A slight disadvantage of this method is that it requires some extra computation for k which is not generated by the standard UKF process. The residual-based method utilizes the measurement of the previous N epoch, whereas the innovation-based method utilizes current and the previous N-1 epoch. Therefore, the innovation sequence and the residual sequence have to be ergodic and stationary over the N steps. Otherwise, the performance of UKF will be degraded.

4. Experimental Results and Discussions

4.1. Test Configuration

The ship-mounted experimental data were collected to evaluate the performance of the in-motion alignment. The experiment was carried out in Yangzi River. The equipped sensors are listed as follows:
  • IMU: Consists of three ring laser gyroscopes with drift rate 0.01° / h(1σ) and three quartz accelerometers with bias 5×10−5g(1σ). Its update rate is 200 Hz.
  • Bottom-lock Doppler: Provides three-axis transformation velocities with accuracy ±5‰ of speed and update rates up to 1 Hz.
  • GPS receiver: Provides velocity with precision of about 0.1m/s, position with precision of about 10 m, and update rates up to 1 Hz.
In the experiment, the IMU and the GPS receiver were set up on a vessel. The DVL module was put beneath 1m underwater. The fixing and level arm parameters of devices are shown in Figure 2. The flowchart of the alignment procedure is shown in Figure 3. As the update rate of the INS is much higher than DVL, the measurement-updating is executed only when the DVL is available. By using the close-loop filtering scheme, the filtering corrections are fed back every measurement-updating. If the filtering epoch is smaller than the sampling window size of innovation or residual, AUKF is working at the mode of UKF.

4.2. Alignment Results by UKF

Figure 4 shows the error curves of heading with initial heading errors from 10 degrees to 180 degrees. As can be seen from the figure, all the error curves converged with time. The proposed alignment model is shown to be effective with any initial attitude.
Figure 5 shows convergence time for different initial heading errors with a converged heading error of less than 0.1 degree. It is clearly shown that larger initial heading error needs more time to converge. In addition, if the initial heading error is larger than 110 degrees. It is hard to guarantee that the final converged heading error would reach 0.1 degree within the 900 s alignment. Figure 6 shows the heading error comparison with extra large initial heading errors. The alignment time was extended to 3,000 s. Though all the estimates converged with time, the heading errors were 0.23 degree, 0.72 degree and 0.73 degree with initial heading error of 160 degrees, 170 degrees and 180 degrees respectively. It is clearly shown that the heading errors converged faster at the beginning and then the speed of the convergence slowed down. This was caused by the fading of the Kalman Filter gain [30]. Therefore, further study is still needed to improve the performance of the UKF. Moreover, providing a relatively accurate initial attitude will be beneficial to the alignment.

4.3. Measurement Noise Covariance Estimation

A test is done by intentionally adding large initial attitude errors (100 degrees for heading, 1 degree for roll and pitch). Figures 7 and 8 show the estimations of measurement noise covariance with different window sizes by the innovation-based and the residual-based AUKF respectively. The trends of the estimations obtained by these two methods are similar. As can be seen from the figures, the estimations of the measurement noise are unstable at the beginning and then converge with time. Finally, all the estimates converge to the value of around 0.01 (m2/s2). It can also be seen from the figures that the estimation oscillation becomes obvious when shorter window sizes are used. It illustrates that a short window size may lead to unstable estimation. This is similar to the conclusions of AKF [16].
A test is done by intentionally adding large initial attitude errors (100 degrees for heading, 1 degree for roll and pitch). Figures 7 and 8 show the estimations of measurement noise covariance with different window sizes by the innovation-based and the residual-based AUKF respectively. The trends of the estimations obtained by these two methods are similar. As can be seen from the figures, the estimations of the measurement noise are unstable at the beginning and then converge with time. Finally, all the estimates converge to the value of around 0.01 (m2/s2). It can also be seen from the figures that the estimation oscillation becomes obvious when shorter window sizes are used. It illustrates that a short window size may lead to unstable estimation. This is similar to the conclusions of AKF [16].
Figures 9 and 10 illustrate the estimates of the measurement noise covariance with different initial R by the AUKF techniques. The window size used for estimation was 100. It is clearly shown that the initial values of R have some impacts only during the filter convergence periods. Finally, all the curves converged and closely matched each other. The requirement of a priori knowledge of the measurement noise covariance is alleviated. In addition, it is clearly shown by Figures 710, the estimated measurement noise covariance obtained by the innovation based method almost approaches that obtained by the residual based case. There is no decided difference between them.

4.4. Performance Evaluation of the Adaptive UKF Techiniques

A test was designed to evaluate the performance of the estimated measurement noise covariance. The initial attitude error was 100 degrees for heading, 1 degree for roll and pitch, respectively. As can be seen from Figures 7 and 8, the values of the estimated measurement noise covariance were around 0.01 (m2/s2). Therefore, this value was applied in the alignment [R = diag(0.01,0.01,0.01) (m2/s2)]. In addition, the alignment was also executed with a larger R [R = diag(0.1,0.1,0.1) (m2/s2)] and a smaller R [R = diag(0.001,0.001,0.001) (m2/s2)]. Figure 11 shows the error curves of the heading with different R. As can be seen from the figure, all the estimates converged with time. But the heading error with R value of 0.01 (m2/s2) converged more rapidly than that with the value of 0.1 (m2/s2) and 0.001 (m2/s2).As shown in Table 1, it took 676 seconds for the heading error to converge within 0.1 degree with R value of 0.01 (m2/s2) while the convergence time for R value of 0.1 (m2/s2) and 0.001 (m2/s2) were 766 s and 800 s respectively. In a sense, the estimated measurement noise covariance is proved to be realistic. Furthermore, it's clearly shown that the measurement noise covariance plays an important part in the performance of the UKF. Once the R is correctly estimated, it can improve the performance of the alignment.
The tests were designed to compare the performance of the UKF and AUKF techniques for their applications in the DVL-aided SINS alignment problem. An example is shown in Figure 12. In the test, the initial attitude error was 100 degrees for heading and 1 degree for leveling. The initial R value applied in UKF and AUKF was 0.01 (m2/s2). The window size for both the innovation-based and residual-based AUKF was 100. As can be seen from Figure 12, the estimated results of the UKF and the AUKF were different after 100 s when the AUKF methods were switched to the adaptive algorithms. Due to the setting of R [0.01 (m2/s2)] given for the UKF was very close to the optimal values obtained from the above tests, the trends of the error curves obtained by the UKF and the AUKF methods are similar. Partial magnification of the heading errors are shown in Figure 13. It is clearly shown that the heading error estimated by the UKF closely matches that predicted by the AUKF methods with the increase of the time. It means that the distinctions between the performances of the UKF and the AUKF are small if a appropriate measurement noise covariance [R = diag(0.01,0.01,0.01) (m2/s2)] is applied. However, in the case of underwater vehicles, it's very difficult to obtain a prior knowledge of the measurement noise covariance due to unstable Doppler measurements. Though a proper R is employed, AUKF methods still performed a little better UKF. Compared with 676 seconds, it took 667 and 664 seconds for the heading errors to converge within 0.1 degree by the innovation based and the residual based AUKF, respectively. As can be seen from Figure 13, their final heading errors are also smaller than UKF. In addition, there is no notable contrast between the AUKF methods. The error curve that predicted by innovation based AUKF is exactly parallel to that obtained by residual based AUKF.
Regarding the attitude obtained from a high precision INS/GPS integration as reference, Figure 14 shows the heading accuracy comparison with different window sizes in the 900 s alignments. It is clearly shown that the estimated results meet the requirement of the alignment of 0.1 degree for the heading error. Choosing an optimal window size is of great importance. However, how to determine the optimal window size is still a challenge issue. It can now only be determined by simulations or experiences. As can be seen from Figure 14, the trends of the curves obtained by the innovation-based and residual-based AUKFs are similar. Accompanying with oscillations, the heading accuracy slowly increased until the window size was around 100. After that, the heading accuracy slowly decreased to 0.03 degree with larger the window size. Then, the heading accuracy remained around 0.03 degree. In addition, the oscillations was obvious when the window size was smaller than 100. This was caused by the unstable estimation of the measurement noise covariance with small window sizes. The result presented in Figure 14 meets our knowledge that the small window size will lead to the biased estimation of the measurement noise covariance while very large window size may cause the adaptive filter loosing the ability of adaptation [29]. Moreover, for this specific case, a slight difference between the innovation-based and residual-based AUKF methods is that the window size of 95 performed best for innovation-based AUKF while window size of 100 performed best for the residual-based AUKF.

5. Conclusions

This paper has presented a new alignment scheme for the DVL-aided SINS in-motion alignment which allows large initial misalignments. From the experimental data, it has been clearly shown that the proposed alignment model can be applied for the DVL-aided SINS in-motion alignment with any initial heading errors. As the measurement noise covariance is of great importance to the performance of the UKF, the covariance-matching methods applied in AKF have been extended for use in the Adaptive UKF. By using innovation-based and residual-based AUKF techniques, the measurement noise covariance can be estimated reliably and hence can improve the performance of the UKF. Its performance has been demonstrated with field experimental data.

Acknowledgments

The first author is sponsored by China Scholarship Council (CSC) for his joint Ph.D. research training at the University of New South Wales, Sydney, Australia.

References

  1. Kinsey, J.C.; Eustice, R.; Whitcomb, L.L. A Survey of Underwater Vehicle Navigation: Recent Advances and New Challenges. Proceedings of IFAC Conference on Maneuvering and Control of Marine Craft, Lisbon, Portugal, September 2006.
  2. Stutters, L.; Liu, H.; Tiltman, C.; Brown, D.J. Navigation technologies for autonomous underwater vehicles. IEEE Trans. Syst. Man Cy. C 2008, 38, 581–589. [Google Scholar]
  3. Li, W.; Wang, J.; Wu, W.; Lu, L. A fast SINS initial alignment scheme for underwater vehicle applications. J. Navig. 2012. [Google Scholar] [CrossRef]
  4. Silson, P.M.G. Coarse alignment of a ship's strapdown inertial attitude reference system using velocity loci. IEEE Trans. Instrum. Meas. 2011, 38, 581–589. [Google Scholar]
  5. Gao, W.; Zhang, X.; Zhao, G.; Ben, Y. A Fine Alignment Method about Doppler-assisted SINS. Proceedings of the 2010 IEEE International Conference on Information and Automation, Harbin, China, 20–23 June 2010; pp. 2333–2337.
  6. Ben, Y.; Zhu, Z.; Li, Q.; Wu, X. DVL Aided Fine Alignment for Marine SINS. Proceedings of the 2011 IEEE International Conference on Mechatronics and Automation, Beijing, China, 7–10 August 2011; pp. 1630–1635.
  7. Liu, Y.; Yu, A.; Zhu, J.; Liang, D. Unscented Kalman filtering in the additive noise case. Sci. China Technol. Sc. 2010, 53, 929–941. [Google Scholar]
  8. Kwangjin, K.; Chan, G. Non-symmetric unscented transformation with application to in-flight alignment. Int. J. Control. Autom. Syst. 2010, 8, 776–781. [Google Scholar]
  9. Wang, Q.; Li, Y.; Rizos, C.; Li, S. The UKF and CDKF for Low-Cost SDINS/GPS in-Motion Alignment. Proceedings of International Symposium on GPS/GNSS, Yokohama, Japan; 2008; pp. 441–448. [Google Scholar]
  10. Seo, W.; Hwang, S.; Park, J.; Lee, J. Precise outdoor localization with a GPS–INS integration system. Robotica 2012. [Google Scholar] [CrossRef]
  11. Shin, E. Estimation Techniques for Low-Cost Inertial Navigation. Ph.D. Thesis, Department of Geomatics Engineering, University of Calgary, Calgary, AB, Canada, May 2005. [Google Scholar]
  12. Petovello, M.G.; Cannon, M.E.; Lachapelle, G. Kalman Filter Reliability Analysis Using Different Update Strategies. Proceedings of the CASI Annual General Meeting, Montreal, QC, Canada, 28–30 April 2003.
  13. Shin, E. Accuracy Improvement of Low Cost INS/GPS for Land Application. M.Sc. Thesis, Department of Geomatics Engineering, University of Calgary, Calgary, AB, Canada, December 2001. [Google Scholar]
  14. Angrisano, A.; Petovello, M.; Pugliano, G. Benefits of combined GPS/GLONASS with low-cost MEMS IMUs for vehicular urban navigation. Sensors 2012, 12, 5134–5158. [Google Scholar]
  15. Zhang, S. An Adaptive Unscented Kalman filter for Dead Reckoning Systems. Proceedings of International conference on Information Engineering and Computer Science, Beijing, China, 19–20 December 2009; pp. 1–4.
  16. Ding, W.; Wang, J.; Rizos, C. Improving adaptive Kalman estimation in GPS/INS integration. J. Navig. 2007, 60, 517–529. [Google Scholar]
  17. Qu, C.; Xu, H.; Tan, Y. SINS/CNS integrated navigation solution using adaptive unscented Kalman filtering. Int. J. Comput. Appl. Technol. 2011, 41, 109–116. [Google Scholar]
  18. Wang, J. Stochastic modeling for RTK GPS/Glonass positioning. J. Inst. Navig. 2000, 46, 297–305. [Google Scholar]
  19. Fang, J.; Yang, S. Study on innovation adaptive EKF for in-flight Alignment of airborne POS. IEEE Trans. Instrum. Meas. 2011, 60, 1378–1388. [Google Scholar]
  20. Kong, X.; Nebot, E.M.; Durrant-Whyte, H. Development of a Non-Linear Psi-Angle Model for Large Misalignment Errors and its Application in INS Alignment and Calibration. Proceedings of the 1999 IEEE International Conference on Robotics and Automation, Detroit, MI, USA, May 1999; pp. 1430–1435.
  21. Kinsey, J.C.; Whitcomb, L.L. Adaptive identification on the group of rigid-body rotations and its application to underwater vehicle navigation. IEEE Trans. Robot. 2007, 23, 124–136. [Google Scholar]
  22. Han, S.; Wang, J. A novel initial alignment scheme for low-cost INS aided by GPS for land vehicle application. J. Navig. 2010, 63, 663–680. [Google Scholar]
  23. Xiong, K.; Zhang, H.; Chan, C. Performance evaluation of UKF-based nonlinear filtering. Automatica. 2006, 42, 261–270. [Google Scholar]
  24. Wu, Y.; Wu, M.; Hu, D.; Hu, X. Unscented Kalman filtering for additive noise case: augmented versus non-augmented. IEEE. Signal Proc. Lett. 2005, 12, 357–360. [Google Scholar]
  25. Van, R. Sigma-Point Kalman Filters for Probabilistic Inference in Dynamic State-Space Models. Ph.D. Thesis, Oregon Health & Science University, Portland, April 2004. [Google Scholar]
  26. Yan, G.; Yan, W.; Xv, D. Application of simplified UKF in SINS initial alignment for large misalignment angles. J. Chin. Inert. Technol. 2008, 16, 253–264. [Google Scholar]
  27. Xu, J.; Jing, Y.; Dimirovski, G.; Ban, Y. Two-Stage Unscented Kalman Filter for Nonlinear Systems in the Presence of Unknown Random Bias. Proceeding of American Control Conference, Seattle, Washington USA, 11–13 June 2008; pp. 3530–3535.
  28. Wan, E.A.; Van, R. The Unscented Kalman Filter for Nonlinear Estimation. Proceeding of the IEEE 2000 Adaptive Systems for Signal Processing, Communications, and Control Symposium, Lake Louise, AB, Canada, 1–4 October 2000; pp. 153–158.
  29. Almagbile, A.; Wang, J.; Ding, W. Evaluating the performance of adaptive Kalman filter methods in GPS/INS integration. CPGPS 2010, 9, 33–40. [Google Scholar]
  30. Yang, Y.; Gao, W. An optimal adaptive Kalman filter. J. Geodesy. 2006, 80, 177–183. [Google Scholar]
Figure 1. DVL-aided IMU alignment scheme.
Figure 1. DVL-aided IMU alignment scheme.
Sensors 13 01046f1
Figure 2. Fixing of experimental devices.
Figure 2. Fixing of experimental devices.
Sensors 13 01046f2
Figure 3. Alignment flowchart.
Figure 3. Alignment flowchart.
Sensors 13 01046f3
Figure 4. Heading error comparison with different initial heading errors.
Figure 4. Heading error comparison with different initial heading errors.
Sensors 13 01046f4
Figure 5. Convergence time comparison with different initial heading errors.
Figure 5. Convergence time comparison with different initial heading errors.
Sensors 13 01046f5
Figure 6. Heading error comparison with extra large initial heading errors.
Figure 6. Heading error comparison with extra large initial heading errors.
Sensors 13 01046f6
Figure 7. Estimation of measurement noise covariance with different window sizes by innovation-based AUKF.
Figure 7. Estimation of measurement noise covariance with different window sizes by innovation-based AUKF.
Sensors 13 01046f7
Figure 8. Estimation of measurement noise covariance with different window sizes by residual-based AUKF.
Figure 8. Estimation of measurement noise covariance with different window sizes by residual-based AUKF.
Sensors 13 01046f8
Figure 9. Estimation of measurement noise covariance with different initial R values by innovation-based AUKF.
Figure 9. Estimation of measurement noise covariance with different initial R values by innovation-based AUKF.
Sensors 13 01046f9
Figure 10. Estimation of measurement noise covariance with different initial R values by residual-based AUKF.
Figure 10. Estimation of measurement noise covariance with different initial R values by residual-based AUKF.
Sensors 13 01046f10
Figure 11. Heading error comparisons with different R.
Figure 11. Heading error comparisons with different R.
Sensors 13 01046f11
Figure 12. Attitude error comparison between UKF and AUKFs with initial attitude error of [1°, 100°, 1°].
Figure 12. Attitude error comparison between UKF and AUKFs with initial attitude error of [1°, 100°, 1°].
Sensors 13 01046f12
Figure 13. Partial magnification of the heading error.
Figure 13. Partial magnification of the heading error.
Sensors 13 01046f13
Figure 14. Heading accuracy comparison with different window sizes.
Figure 14. Heading accuracy comparison with different window sizes.
Sensors 13 01046f14
Table 1. Performance comparisons with different R values.
Table 1. Performance comparisons with different R values.
R Value (m2/s2)Heading Accuracy (°)Convergence Time (s)
1e-10.0629766
1e-20.0282676
1e-30.0367800

Share and Cite

MDPI and ACS Style

Li, W.; Wang, J.; Lu, L.; Wu, W. A Novel Scheme for DVL-Aided SINS In-Motion Alignment Using UKF Techniques. Sensors 2013, 13, 1046-1063. https://doi.org/10.3390/s130101046

AMA Style

Li W, Wang J, Lu L, Wu W. A Novel Scheme for DVL-Aided SINS In-Motion Alignment Using UKF Techniques. Sensors. 2013; 13(1):1046-1063. https://doi.org/10.3390/s130101046

Chicago/Turabian Style

Li, Wanli, Jinling Wang, Liangqing Lu, and Wenqi Wu. 2013. "A Novel Scheme for DVL-Aided SINS In-Motion Alignment Using UKF Techniques" Sensors 13, no. 1: 1046-1063. https://doi.org/10.3390/s130101046

APA Style

Li, W., Wang, J., Lu, L., & Wu, W. (2013). A Novel Scheme for DVL-Aided SINS In-Motion Alignment Using UKF Techniques. Sensors, 13(1), 1046-1063. https://doi.org/10.3390/s130101046

Article Metrics

Back to TopTop