Next Article in Journal
High Sensitivity Refractometer Based on TiO2-Coated Adiabatic Tapered Optical Fiber via ALD Technology
Next Article in Special Issue
A Matrix-Based Proactive Data Relay Algorithm for Large Distributed Sensor Networks
Previous Article in Journal / Special Issue
Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements

1
Faculty of Navigation and Naval Weapons, Polish Naval Academy, Smidowicza 69, Gdynia 81-103, Poland
2
Faculty of Civil and Environmental Engineering, Department of Geodesy, Gdansk University of Technology, Narutowicza 11/12, Gdansk 80-233, Poland
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(8), 1279; https://doi.org/10.3390/s16081279
Submission received: 16 May 2016 / Revised: 1 August 2016 / Accepted: 2 August 2016 / Published: 15 August 2016
(This article belongs to the Special Issue Advances in Multi-Sensor Information Fusion: Theory and Applications)

Abstract

:
The article presents a study of the accuracy of estimating the position coordinates of BAUV (Biomimetic Autonomous Underwater Vehicle) by the extended Kalman filter (EKF) method. The fusion of movement parameters measurements and position coordinates fixes was applied. The movement parameters measurements are carried out by on-board navigation devices, while the position coordinates fixes are done by the USBL (Ultra Short Base Line) system. The problem of underwater positioning and the conceptual design of the BAUV navigation system constructed at the Naval Academy (Polish Naval Academy—PNA) are presented in the first part of the paper. The second part consists of description of the evaluation results of positioning accuracy, the genesis of the problem of selecting method for underwater positioning, and the mathematical description of the method of estimating the position coordinates using the EKF method by the fusion of measurements with on-board navigation and measurements obtained with the USBL system. The main part contains a description of experimental research. It consists of a simulation program of navigational parameter measurements carried out during the BAUV passage along the test section. Next, the article covers the determination of position coordinates on the basis of simulated parameters, using EKF and DR methods and the USBL system, which are then subjected to a comparative analysis of accuracy. The final part contains systemic conclusions justifying the desirability of applying the proposed fusion method of navigation parameters for the BAUV positioning.

1. Introduction

Determining the position coordinates with high accuracy during the process of conducting marine underwater navigation is becoming a significant problem. This is mainly due to the fact that the Global Navigational Satellite System (GNSS) cannot be used. Sometimes using satellite navigation during emersion could be also difficult, especially in port areas, where infrastructure facilities and other vessels can block and reflect navigation signals [1,2]. Due to these difficulties, it becomes necessary to use dedicated navigation systems (NS) in such cases [3,4,5,6,7,8,9,10]. These dedicated navigation systems generally:
  • use devices to carry out navigational parameter measurements in water, such as log, hydrostatic pressure sensor, or echo-sounder,
  • estimate the position coordinates as a result of the fusion of various navigational parameters, e.g., speed with course or submergence [11,12,13,14,15,16,17,18,19,20].
In the case of BAUV, a selection of devices and methods for determining the position can be made dependent on the parameters of its operation, shape, and size of the hull (construction), and type of propulsion used.
A BAUV, imitating fish in terms of the design and movement, is currently being prepared at the Polish Naval Academy (PNA) (Figure 1).
Starting from the bow to the tail fin, the BAUV consists of the following modules [9]:
  • module with camera and “looking” forward echo sounder (“wet” compartment),
  • module of the sensors (in the upper part: USBL, hydro modem and “looking” up echo sounder; in the bottom part: sonar and “looking” down echo sounder—“wet” compartment),
  • module of the lateral fins (“dry” compartment),
  • module of electronics and batteries (the batteries located in the bottom part, which can move along the longitudinal axis of the vehicle giving the possibility of trimming; two computers PC-104 and power management system mounted in the upper part—the whole compartment is “dry”),
  • module of the caudal fin consisting of two segments (the tail fin driven by an electric motor with nominal power 250 W, rotary motion converted into an oscillating motion).
The range of swimming is designed for about 2 nautical miles and assumed optimum cruise speed is around 1 m/s (2 knots). The swimming range is an important operating parameter and it is a determining factor in the construction of the designed NS.
Several NS variants have been designed for the conceptual stage of the BAUV. They differ from each other, mainly in applied navigation devices and in algorithms of the position coordinates estimation. The conceptual NS described in this publication consists of components installed on the BAUV:
  • GPS Aided Inertial Navigation System (GPS/INS) “VN-200” [21]—due to the operating environment, it is not possible to use GPS, therefore this subsystem is referred to as INS “VN-200” in subsequent parts of this article,
  • “ALIZE” electro-magnetic log [22],
  • ”WIKA S-20” hydrostatic pressure sensor [23],
  • USBL “MicronNav System” transponder system [24],
  • “MICRON DATA MODEM” slave hydromodem [25],
and components located on the seabed:
  • USBL “MicronNav System” transceiver system,
  • “MICRON DATA MODEM” master.
The USBL system will be dedicated to measure the approximate position coordinates (sending to the BAUV by hydromodem). The INS “VN-200” will measure the course, while the log will measure speed, and the pressure sensor will be used to submergence determination. These parameters will then be used to estimate submerged BAUV position by the extended Kalman filter. The BAUV moves using wave motion.
The NS constructed in this way will be subject to a comparative evaluation of the dead reckoning NS (that uses only the course and speed measurement), and the NS based only on the USBL system. The basic criterion for this assessment is the accuracy of the estimated position coordinates.

2. Evaluation of the Accuracy of Determining the Position Coordinates Using the USBL System

Determining the horizontal coordinates using USBL system is based on the measurement of two parameters: direction α (yaw) and the distance between the Transceiver and the Transponder r (Figure 2).
The mean error M x y and the mean error ellipse (i.e., the length of the a and b axes) of determined horizontal coordinates ( x ,   y ) using the USBL system can be calculated by applying the law of mean errors propagation [26,27,28]. Knowing the functions of the individual test result:
α = arc tg x y
r =   x 2 + y 2  
the equations of position lines mean error can be represented by:
σ l ( α ) = σ α / ( α x ) 2 + ( α y ) 2
σ l ( r ) = σ r / ( r x ) 2 + ( r y ) 2
that enables determination of:
M x y = 1 sin θ σ l ( α ) 2 + σ l ( r ) 2 = ( σ α · r ) 2 + σ r 2
a = σ α · r
b = σ r
where θ —the intersection angle of the position lines. In the case of USBL it always equals 90°; σ α —α (yaw) mean measurement error; σ r r (distance) mean measurement error; r —distance between the Transceiver and Transponder; a —the length of the major axis of the mean error ellipse; b —the length of the minor axis of the mean error ellipse.
Based on the dependence Equation (5), information provided by the manufacturer, and the values of the mean errors of measurement, σ α = 3 ° σ r = 0.2 m, an area map has been compiled on the accuracy of position coordinates, determined by the Tritech USBL “MicronNav System” (Figure 3) [24].
Figure 3 shows that the mean errors of position coordinate increases with the measured distance. It reaches 26.2 m at a distance of 500 m. This is the distance limit at which the USBL “MicronNav System” can work [24].
Figure 4 shows the mean error ellipse of position coordinates determined in the same direction, in one hundred meter intervals.
It shows clearly that the coordinate error increases substantially in a perpendicular direction to the measurement direction with increasing distance. However, it is constant in direction of the measurement being carried out, and is only 0.2 m.

3. Evaluating the Accuracy of Determining the Position Coordinates Using Dead Reckoning

If the impact of current, drift, and waving of sea on the BAUV movement can be avoided, the only sources of errors are the devices used to dead reckoning. In the case of the BAUV, these are the Inertial Navigation System, INS “VN-200”, which determines course over ground (COG) [21] and the electromagnetic log “ALIZE”, which determines speed through water (STW). In this case, the STW is equal to the speed over ground (SOG) [22].
By measuring using these devices COG ( k ) , SOG ( k ) at the time of ( k ) , vector x h of horizontal coordinates and their accuracy can be calculated for the time ( k + 1 ) . Functions f ( x h ( k ) , u h ( k ) , w h ( k ) ) can be used for this purpose describing the non-linear BAUV movement model:
x h ( k + 1 )   = [ f x ( k + 1 ) f y ( k + 1 ) ] = [ x ( k ) y ( k ) ] + [ Δ t ( k ) · SOG ( k ) · sin COG ( k ) Δ t ( k ) · SOG ( k ) · cos COG ( k ) ] + [ w x ( k ) w y ( k ) ]
where x ( k ) , y ( k ) —dead reckoned horizontal coordinates of BAUV position at the time of k ; COG ( k ) —the BAUV course over ground at the time of k ; SOG ( k ) —the BAUV speed over ground at the time of k ; Δ t ( k ) —duration between the time of k and k + 1; w x ( k ) , w y ( k ) —the so-called intentional interference in determining the coordinates at the time of k (expressed as a zero mean normal distribution N[0, 1]); and calculated on the basis of Equation (8) the covariance matrix
P ( k + 1 ) = F x ( k + 1 ) P ( k ) F x ( k + 1 ) T + Q ( k )
here: F x ( k + 1 ) = [ f x x f x y f y x f y y ] = [ 1 0 0 1 ] —the so-called matrix system is calculated as Jacobian matrix from the function f ( x ( k ) , u ( k ) , 0 ) ; P ( k ) —the covariance matrix, determined at the time of k .
Q ( k ) = [ f x C O G f x SOG f y COG f y SOG ] [ σ C O G 2 0 0 σ S O G 2 ] [ f x C O G f x SOG f y COG f y SOG ] T
σ C O G —the course over ground mean measurement error; σ S O G —the speed over ground mean measurement error.
Knowing the covariance matrix of initial (the previous) position coordinates P( k ) and the covariance matrix of the vector of growth coordinates P( k + 1 / 1 ), taking into account the errors resulting from the operation of dead reckoning navigation, the covariance matrix of the coordinates of the current BAUV position can be determined
P ( k + 1 ) = P ( k ) + P ( k + 1 / 1 )
where
P ( k + 1 ) = [ σ x 2 σ x y σ y x σ y 2 ]
σ x 2 = Δ t 2 · ( S O G ( k ) · σ C O G · cos COG ( k ) ) 2 + Δ t 2 · ( σ S O G · sin COG ( k ) ) 2
σ y 2 = Δ t 2 · ( S O G ( k ) · σ C O G · sin COG ( k ) ) 2 + Δ t 2 · ( σ S O G · cos COG ( k ) ) 2
σ x y = σ y x = Δ t 2 · sin ( 2 · COG ( k ) ) · ( σ S O G 2 ( S O G · σ C O G ) 2 ) / 2
Using this matrix, in turn, the mean error of coordinates can easily be calculated
M x y = σ x 2 + σ y 2
and the parameters of the mean error ellipse:
a = 1 2 · ( σ x 2 + σ y 2 + p )
b = 1 2 · ( σ x 2 + σ y 2 p )
τ = 1 2   arc tg 2 σ x y   σ x 2 σ y 2
where
p = ( σ x 2 σ y 2 ) 2 + 4 · σ x y 2
τ —direction angle of ellipse.
Relying on the information provided by the manufacturer of log “ALIZE” that σ S O G = 0.5 kn and determined arbitrarily σ C O G = 3 ° for INS “VN-200”, using Equations (10) and (11), a mean error graph of position coordinates determined using dead reckoning was prepared (Figure 5).
Figure 5 shows that the mean error of position coordinates increases linearly as a function of distance. Every hundred meters it increases constantly at 25.54 m, reaching a value of 127.7 m after traveling 500 m.
Figure 6 shows the mean error ellipses of position coordinates determined using the dead reckoning method after traveling 200 and 500 m (their dependencies No. 12–13 have been calculated).
Figure 6 shows that the largest coordinate error occurs in the direction according to COG. The length of the axis b is equal to the length of the mean error ellipse of coordinates determined using USBL system (Figure 4). This is due to the fact that in both methods of determining the position coordinates the positioning line is a result of the direction measurement that is burdened with the same error value σ C O G   = σ α .

4. The Genesis of the Problem

The choice between the USBL system and the method of dead reckoning based on the accuracy criterion is difficult. In some cases, the USBL system can be identified as being better (e.g., when measurements are carried out in small distances from the transceiver), while in others, the method of dead reckoning (e.g., immediately after updating coordinates by GNSS fixes). If a USBL system is selected, you give up additional COG and SOG. If you choose the DR method measurements r , α , β (range, yaw, pitch) are not taken into account. It is certain that skilful use of a greater number of measurements can improve their accuracy [29].
In the case of the NS BAUV, the different availability of measurements is assumed, i.e.,
  • at regular intervals of time and with high frequency—the course and speed determined by navigation devices installed on the BAUV (not less than once per second),
  • at irregular intervals and with low frequency—the direction and distance to the Transceiver of the appointed USBL system and transmitted by hydro-modem (less than once per second—because of the difficulties in carrying out USBL system measurements and sending the results of these measurements by hydro-modem caused interference in the propagation of acoustic waves in the water).
These parameters will then be used to determine the position coordinates using the extended Kalman filter (EKF).
In its action, the filter may combine the various measurements (e.g., the course with the distance to a navigation mark), carried out at different time intervals, taking into account their errors—adjusted dynamically by the matrix of weights. Thanks to this mechanism, measurements which are obviously erroneous (e.g., carried out in the USBL system, on the border of its operation) can be firmly suppressed, and accurate measurements amplified (e.g., SOG and COG determined in a short period of time from the moment of updating the coordinates by GNSS fixes).
However, will the coordinates obtained by the EKF method, a combination of parameter measurements of BAUV motion and the direction and distance of the USBL system, be more accurate than either the coordinate USBL system or the DR method in all cases?
The rest of this article presents a study on the positioning accuracy of BAUV by the EKF method compared to the USBL and DR positioning system, aimed at finding answers to the above question.

5. Description of Data Fusion Using the Extended Kalman Filter

Let us carry out a fusion of data using the extended Kalman filter combining COG ( k ) and SOG ( k ) measurements, changes in submergence Δ z ( k ) carried out every one second with 3D coordinate system measurements of the USBL system, carried out at different time intervals (in some moments k + 1 ).
Let us assume that the results of the USBL measurement system will create the so-called vector of observation, z ( k + 1 ) = [ r ( k + 1 )   α ( k + 1 )   β ( k + 1 ) ] T , described by the function [30]:
h ( x ( k + 1 ) , ν ( k + 1 ) ) = [ f r ( k + 1 ) f α ( k + 1 ) f β ( k + 1 ) ] = [ r ( k + 1 ) α ( k + 1 ) β ( k + 1 ) ] + [ ν r ( k + 1 ) ν α ( k + 1 ) ν β ( k + 1 ) ]
where
r ( k + 1 ) = ( x 1 ( k + 1 ) x ( k + 1 ) ) 2 + ( y 1 ( k + 1 ) y ( k + 1 ) ) 2 + ( z 1 ( k + 1 ) z ( k + 1 ) ) 2
α ( k + 1 ) = arc tg x 1 ( k + 1 ) x ( k + 1 ) y 1 ( k + 1 ) y ( k + 1 )
β ( k + 1 ) =  arc tg z 1 ( k + 1 ) z ( k + 1 ) ( x 1 ( k + 1 ) x ( k + 1 ) ) 2 + ( y 1 ( k + 1 ) y ( k + 1 ) ) 2
ν r ( k + 1 ) , ν α ( k + 1 ) , ν β ( k + 1 ) —values of measurement errors (with zero mean normal distribution); x ( k + 1 ) , y ( k + 1 ) , z ( k + 1 ) —fixed coordinates of transceiver position; x 1 ( k + 1 ) , y 1 ( k + 1 ) , z 1 ( k + 1 ) —variable coordinates of the transponder position (BAUV). And on the basis of COG ( k ) , SOG ( k ) measurements and the submergence changes Δ z ( k ) , vector x of BAUV’s 3D position coordinates, using the function f ( x ( k ) , u ( k ) , w ( k ) ) :
x ( k + 1 )   = [ f x ( k + 1 ) f y ( k + 1 ) f z ( k + 1 ) ] = [ x ( k ) y ( k ) z ( k ) ] + [ Δ t ( k ) · SOG ( k ) · sin COG ( k ) Δ t ( k ) · SOG ( k ) · cos COG ( k ) Δ z ( k ) ] + [ w x ( k ) w x ( k ) w z ( k ) ]
In pursuing data fusion, based on functions h ( x ( k + 1 ) , 0 ) i f ( x ( k ) , u ( k ) , 0 ) we calculate the estimated vector of coordinates x ^ ( k + 1 ) and the covariance matrix P ( k + 1 ) , using the following dependencies:
x ^ ( k + 1 ) = f ( x ( k ) , u ( k ) , 0 )
P ( k + 1 ) = F ( k + 1 ) P ( k ) F ( k   +   1 ) T + Q ( k )
y ( k + 1 ) = z ( k + 1 ) h ( x ^ ( k + 1 ) , 0 )
S ( k + 1 ) ( H ( k + 1 ) P ( k + 1 ) H ( k + 1 ) T + R ( k + 1 )
K ( k + 1 ) = P ( k + 1 ) H ( k + 1 ) T   S ( k + 1 ) 1
x ^ ( k + 1 ) = x ^ ( k + 1 ) + K ( k + 1 ) y ( k + 1 )
P ( k + 1 ) + ( I K ( k + 1 ) H ( k + 1 ) ) P ( k + 1 )
where x ^ ( k + 1 )   i   P ( k + 1 ) —estimated vector of coordinates of the BAUV position and its covariance matrix, determined a priori for the time of k + 1; x ^ ( k + 1 )   i   P ( k + 1 ) —estimated vector of coordinates of the BAUV position and its covariance matrix, determined a posteriori for the time of k + 1, F ( k + 1 ) = [ f x x f x y f x z f y x f y y f y z f z x f z y f z z ] = [ 1 0 0 0 1 0 0 0 1 ] —the so-called matrix system is calculated as a function of the Jacobian matrix f ( x ( k ) , u ( k ) , 0 ) , Q ( k ) = [ σ x 2 σ x y σ x z σ y x σ y 2 σ y z σ z x σ z y σ z 2 ] = [ f x C O G f x SOG f x Δ z f y COG f y SOG f y Δ z f z COG f z SOG f z Δ z ] [ σ C O G 2 0 0 0 σ S O G 2 0 0 0 σ Δ z 2 ] [ f x C O G f x SOG f x Δ z f y COG f y SOG f y Δ z f z COG f z SOG f z Δ z ] T —matrix of noises of the state vector at the time of k (with adopted mean errors σ C O G , σ S O G , σ Δ z of COG and SOG measurements and changes in submergence),
σ x z = σ y z = σ z x = σ z y = 0
σ z 2 = σ Δ z 2
P ( k )   —updated covariance matrix of the state vector used in subsequent time k + 1, H ( k + 1 ) = [ f r x f r y f r z f α x f α y f α z f β x f β y f β z ] = [ x 1 ( k + 1 ) x ( k + 1 ) l 3 ( k + 1 ) y 1 ( k + 1 ) y ( k + 1 ) l 3 ( k + 1 ) z 1 ( k + 1 ) z ( k + 1 ) l 3 ( k + 1 ) y 1 ( k+1 )y( k+1 ) l 2 ( k+1 ) x 1 ( k+1 )x( k+1 ) l 2 ( k+1 ) 0 ( x 1 ( k+1 )x( k+1 ) )( z 1 ( k+1 )z( k+1 ) ) l 2 ( k+1 ) l 3 ( k+1 ) ( y 1 ( k+1 )y( k+1 ) )( z 1 ( k+1 )z( k+1 ) ) l 2 ( k+1 )    l 3 ( k+1 ) l 2 ( k+1 ) l 3 ( k+1 ) ] —Jacobian matrix with function h ( x ( k ) , 0 ) ,
l 3 ( k + 1 ) = ( x 1 ( k + 1 ) x ( k + 1 ) ) 2 + ( y 1 ( k + 1 ) y ( k + 1 ) ) 2 + ( z 1 ( k + 1 ) z ( k + 1 ) ) 2
l 2 ( k + 1 ) = ( x 1 ( k + 1 ) x ( k + 1 ) ) 2 + ( y 1 ( k + 1 ) y ( k + 1 ) ) 2
R ( k + 1 ) =   [ σ r 2 0 0 0 σ α 2 0 0 0 σ β 2 ] —matrix of noises of the observation vector in time of k + 1 (adopted mean measurement errors: σ r —range, σ α —yaw, σ β —pitch); I —identity matrix.

6. Research and Analysis of the Obtained Results

Let us now do an experiment to software simulation of the BAUV passage along the straight section with simultaneous:
  • generating COG —measurements of INS carried out by “VN-200”, SOG —carried out with log “ALIZE” and Δ z —carried out with a hydrostatic pressure sensor “S-20 WIKA”;
  • generating by r , α , β measurements—carried out by the USBL “MicronNav System”;
  • calculating the BAUV estimated coordinates of positions, parallel to the USBL system, with the method of dead reckoning and using the extended Kalman filter.
As a measuring testing ground, we assume a sea area around the USBL transponder system (Figure 7).
Let the BAUV swims near the USBL transceiver system (located at a depth of 6 m) overcoming the test section with COG =   220 ° , SOG =   1  m / s , submerged at z   =   3  m . Meanwhile, τ = 0   s ,   600   s simulated measurements will be carried out by the on-board BAUV and USBL system, i.e.,
  • at fixed moments k + 1 (every one second) COG ( k + 1 ) =   COG w ( k + 1 ) + Δ C O G , SOG ( k + 1 ) = SOG w ( k + 1 ) + Δ S O G and Δ z ( k + 1 ) =   Δ z w ( k + 1 ) + Δ Δ z ;
  • at some moments   k + 1   (every 5, 20 and 100 s) r ( k + 1 ) =   r w ( k + 1 ) + Δ r , α ( k + 1 ) = α w ( k + 1 ) + Δ α , β ( k + 1 ) = β w ( k + 1 ) + Δ β .
Each simulated value of the measurement is obtained by adding the reference measurement and the error (e.g., simulated COG ( k + 1 ) =   COG w ( k + 1 ) + Δ C O G ). Measurement error Δ will be treated as a random variable with uniform distribution, the value of which will be contained within three mean errors (3σ) of measurements made in the device (i.e., the probability of its occurrence will be at 99%).
Reference measurements of the USBL system r w ( k + 1 ) , α w ( k + 1 ) , β w ( k + 1 ) will be determined relative to the reference positions consecutively occupied by the BAUV, which moves with undisturbed movement. On the other hand, the BAUV reference position will be dead reckoned (DR) on the basis of reference COG w ( k + 1 ) , SOG w ( k + 1 ) and Δ z w ( k + 1 ) , which will be located exactly at the measuring section (see Figure 7).
For realigning the experiment for calculating the accepted mean error values close to the actual values (provided by the manufacturer and determined on the basis of their own—selected arbitrarily):
  • σ C O G = 8 ° —for INS “VN-200”,
  • σ S O G = 0.5 kn—for the electromagnetic log “ALIZE”,
  • σ Δ z = 0.25  m —for the hydrostatic pressure sensor “WIKA S-20”,
  • σ r = 0.2  m , σ α = σ β = 3 ° —for the USBL “MicronNav System” system.
The experiment will commence at the moment of k = 0 and will be completed when the BAUV swims about 600 m. The extended Kalman filter will be “fine-tuned” at the start of the experiment and the coordinates of the initial BAUV position will be determined with a mean error of about 3 m (it is assumed that the coordinates were determined before submergence using GNSS) [31]. In the calculations, changes of COG , SOG , and Δ z caused by the BAUV wave motion (swinging, on a sine curve) will be omitted.

6.1. Test No. 1

The test consisted of making three and then one hundred passages along the test section.
COG , SOG and Δ z were measured every one second.
r , α , β were measured every 5 s.

First Passage

In Figure 8 determined routs using different methods of positioning are presented. Figure 9, Figure 10, Figure 11 and Figure 12 present graphs depicting simulated measurement errors ( Δ C O G ,   Δ S O G ,   Δ r ,   Δ α , Δ β ) during the first passage of the test section (these graphs will not be displayed during the presentation of the results of further testing). Figure 9, Figure 10, Figure 11 and Figure 12 clearly show that the dispersion of the measurement errors are random and limited. The error value does not exceed 3 σ fluctuating:
  • Δ C O G in the range of 15 ° ,   15 ° ,
  • Δ S O G in the range of 0.8  m / s ,   0.8  m / s ,
  • Δ α , Δ β in the range of 9 ° ,   9 ° ,
  • Δ r in the range of 0.6   m ,   0.6  m .
In Figure 13, for comparison, three graphs of the reference position distance from the position estimated by the extended Kalman filter and dead reckoning methods and determined by the USBL system are presented.
They show that the accuracy of the estimated position coordinates of EKF method is the same during the entire test duration. The accuracy of the estimated position coordinates using the DR method decreases as a function of time. On the other hand, the accuracy of the USBL position system grows in the time interval <400 s, 500 s>—it is then that the distance between the transceiver and the transponder does not exceed 50 m.
Table 1 shows the statistical parameters characterizing the accuracy of coordinates of all positions set at the test section.
The results in Table 1 clearly show that the best values of the statistical parameters were obtained by the EKF method.
On the basis of simulation measurements during two successive passages along the test section, the following route graphs were prepared and the reference position distance from the estimated position (Figure 14, Figure 15, Figure 16 and Figure 17) as well as the statistical parameters (Table 2 and Table 3).

Second Passage

Figure 14. Determined routes (single passage).
Figure 14. Determined routes (single passage).
Sensors 16 01279 g014
Figure 15. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 15. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g015
Table 2. Parameters of statistical positioning errors by EKF, DR, and USBL system (single passage).
Table 2. Parameters of statistical positioning errors by EKF, DR, and USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.16.42.3
DR0.122.611.8
USBL0.281.623.6

Third Passage

Figure 16. Determined route (single passage).
Figure 16. Determined route (single passage).
Sensors 16 01279 g016
Figure 17. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 17. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g017
Table 3. Parameters of statistical positioning errors by EKF, DR and USBL system (single passage).
Table 3. Parameters of statistical positioning errors by EKF, DR and USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.15.31.8
DR0.214.04.7
USBL0.786.121.6
On the basis of graphs analysis of Figure 14, Figure 15, Figure 16 and Figure 17, as well as statistical parameters from Table 2 and Table 3, it can be stated that the accuracy of estimated position coordinates using the EKF method is higher than those estimated using the DR method and determined by the USBL system. The reference position distance from the estimated one by EKF method does not exceed 7 m. Unfortunately, on the graphs it can be noticed that time intervals have occurred in which the reference position is further from the position estimated by EKF method, in comparison to the estimated position by DR method or determined using the USBL system. One can ask oneself a question here: how often do such time intervals occur? Generalized results of statistical analysis can answer that question. It is based on two additional indicators.
The former describes the average ξ distance to the reference position with estimated calculation on the basis of n distance ξ obtained for n passages along the test section at the same moment of   k , in accordance with the formula:
ξ ( k ) = i = 0 n ζ i ( k )
Figure 18 and Figure 19 presented the graph of the determined route, average distance to the reference position with the estimated, established on the basis of a hundred distance measurements carried out at the same moment, on each of the hundred BAUV passages of the test section.

One Hundred Passages

Generalizing the results of test No. 1, based on Figure 19, it can be stated that the average value of the reference position distance from the estimated position is constant in the case of EKF, it increases linearly for DR and decreases and then increases linearly in the case of the USBL system.
The second indicator to the generalized analysis is based on the histogram showing the frequency of the same distance to the reference position from the estimated position. The histogram is presented in Figure 20.
It results from the histogram that a two meter distance from the estimated position by EKF occurred most often. In addition, as a function of the distance to 5 m, EKF histogram columns occur most often. DR and USBL histograms are flattened as a function of distance. This proves that the incidence of very scattered distance values (up to 15 m) is highly probable. The statistical parameters of the positioning errors using EKF, DR methods and the USBL system are presented in Table 4.
The results in Table 4 clearly show that the best values of the statistical parameters were obtained by the EKF method.

6.2. Test No. 2

The test depended on carrying out three, followed by a hundred passages in the test section. COG , SOG and Δ z were measured every one second. r , α , β were measured every 20 s.

First Passage

The statistical parameters of the positioning errors using EKF, DR methods and the USBL system for single passage are presented in Table 5 (for first passage) in Table 6 (for second passage) and in Table 7 (for third passage).

Second Passage

Table 6. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Table 6. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.17.52.5
DR0.126.79.7
USBL1.464.321.0

Third Passage

The analysis of Figure 21, Figure 22 and Figure 23 shows that the accuracy of the estimated position coordinates using EKF method improves almost every 20 s, when measurements from the USBL system are included in the calculations. Even if the position determined by the USBL system is significantly distant from the reference position, the accuracy of the estimated position using EKF method increases. This may be due to the fact that EKF suppresses the significant and obviously erroneous measurement direction of α ( k + 1 ) , β ( k + 1 ) and reinforces the accurate measurement of the distance of r ( k + 1 ) .
On the basis of analysis of the graphs in Figure 21, Figure 22 and Figure 23, and the statistical parameters from Table 6 and Table 7, it can be concluded that the accuracy of the estimated position coordinates using EKF method is greater than the estimated ones using DR method and determined by USBL system. The maximum reference position distance from the estimated ones using EKF method does not exceed 8 m. However, as in test No. 1, it can be seen on the graphs that there were intervals in which the reference position is more distant from the estimated position by the EKF method, compared to the estimated position by the DR method or that determined by the USBL system.
Table 7. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Table 7. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.36.12.2
DR0.315.36.2
USBL1.372.621.1

One Hundred Passages

Figure 24 presents determined routs during one hundred passages using different methods of positioning. Graph of the reference position distance is presented in Figure 25.
Generalizing the results of Test No. 2, they are similar to Test No. 1 based on the indicators presented in Figure 26 and Table 8.

6.3. Test No. 3

The test consisted of carrying out three, and then a hundred passages along the test section. COG , SOG and Δ z were measured every one second. r , α , β were measured every 100 s. The statistical parameters of the positioning errors using EKF, DR methods and the USBL system for single passage are presented in Table 9 (for first passage) in Table 10 (for second passage) and in Table 11 (for third passage).

First Passage

Table 9. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Table 9. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance from the Reference Position (m)Maximum Distance from the Reference Position (m)Average Distance from the Reference Position (m)
EKF0.111.64.0
DR0.333.822.0
USBL3.437.518.7

Second Passage

Table 10. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Table 10. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance from the Reference Position (m)Maximum Distance from the Reference Position (m)Average Distance from the Reference Position (m)
EKF0.511.44.9
DR0.528.712.0
USBL7.620.014.0

Third Passage

Table 11. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Table 11. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance from the Reference Position (m)Maximum Distance from the Reference Position (m)Average Distance from the Reference Position (m)
EKF0.29.14.0
DR0.220.38.9
USBL3.522.011.7
Graphs of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system for single passage are presented in Figure 27, Figure 28 and Figure 29. The analysis of these Figures shows that the accuracy of the estimated position coordinates using EKF method improves almost every 100 s, when measurements from the USBL system are included in the calculations. The accuracy of the estimated position coordinates using the EKF method is greater than the estimated ones using the DR method and determined by the USBL system. The reference position distance from the estimated ones using EKF method is slightly larger than in previous tests and reach about 11 m. This is due to the long period of estimating the position using EKF without taking measurements of the USBL system.

One Hundred Passages

Figure 30 presents determined routs during one hundred passages using different methods of positioning. Statistical parameters of positioning errors using EKF, DR methods and the USBL system for these passages are presented in Table 12.
The generalized results of the Test No. 3, clearly show that the estimation of position coordinates using the EKF method gives the best results. In Figure 31, a marked improvement in the accuracy of the position can be seen when calculating additional measurements from the USBL system (at 100 s). The histogram in Figure 32 confirms the exact estimated position coordinates by the EKF method can be most frequently expected.

7. Conclusions

NS of the BAUV, which uses only the dead reckoning method based on the course and speed measurements, determines the position coordinates which error has quickly built up over time. In a study, after 500 s (corresponding to travelling 500 m) the average error of the position coordinates reached a value of 127.7 m (Figure 5).
In turn, the USBL system ensures high accuracy of determining the position coordinates, though only at short distances from the transceiver. In studies, at a distance of 500 m, the average error of determining the position coordinates will reach a value of 26.2 m (Figure 3).
Combining the EKF measurement method in the calculations carried out by devices installed on the BAUV (log, INS, hydrostatic pressure sensor) with distance and direction measurements carried out in the USBL system allows for obtaining the position coordinates significantly more accurate. This happens mainly because the position coordinates obtained are based on low accurate COG and SOG measurements and may be adjusted every now and then by more precise distance measurements of r carried out by the USBL system. This is confirmed by the full results of the three tests presented in the article.
The presented research concerned only method of BAUV positioning during undisturbed motion along the straight line. The obtained results are based only on simulated measurements. Further research will focus on tests in real environment. Sea currents and waving of sea will be taken into consideration. Furthermore, positioning accuracy of BAUV will be assessed during movement along composed trajectories including straight sections as well as curves lines.

Author Contributions

The work presented in the paper was carried out in equal contribution of all authors. Krzysztof Naus conceived the concept and idea. Authors together implemented software and carried out the experiments. The manuscript was also written as a common work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nowak, A. The Proposal to “Snapshot” RAIM method for GNSS vessel receivers working in poor space segment geometry. Pol. Marit. Res. 2015, 22, 3–8. [Google Scholar] [CrossRef]
  2. Janowski, A.; Nowak, A.; Przyborski, M.; Szulwic, J. Mobile Indicators in GIS and GPS positioning accuracy in cities. In Rough Sets and Intelligent Systems Paradigms, Lecture Nots in Computer Science; Springer International Publishing: Berlin, Germany, 2014; pp. 309–318. [Google Scholar]
  3. Naus, K. Electronic navigational chart as an equivalent to image produced by hypercatadioptric camera system. Pol. Marit. Res. 2015, 22, 3–9. [Google Scholar] [CrossRef]
  4. Naus, K.; Waz, M. Precision in determining ship position using the omnidirectional map to visual shoreline image comparative method. J. Navig. 2016, 69, 391–413. [Google Scholar] [CrossRef]
  5. Naus, K.; Waz, M. A simplified navigational chart pyramid dedicated to an autonomous navigational system. Pol. Hyperb. Res. 2012, 3, 139–163. [Google Scholar]
  6. Naus, K.; Szymak, P. Accuracy of positioning autonomous biomimetic underwater vehicle using additional measurement of distances. Annu. Navig. 2015, 22, 65–79. [Google Scholar]
  7. Ribas, D.; Ridao, P.J.; Neira, J.; Tardos, J.D. SLAM using an imaging sonar for partially structured underwater environments. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Beijing, China, 9–15 October 2006; pp. 5040–5045.
  8. Syder, F.D.; Morris, D.D.; Haley, P.H.; Collins, R.T.; Okerholm, A.M. Northrop Grumman. Autonomous River Navigation. Proc. SPIE Mob. Robots XVII 2004, II, 221–232. [Google Scholar]
  9. Szymak, P.; Praczyk, T.; Naus, K.; Szturomski, B.; Malec, M.; Morawski, M. Autonomous Underwater Vehicle with Silent Undulating Propulsion for Underwater ISR; Technical Report on Stage No. 2 from Development Project No. DOBR-BIO4/033/13015/2013 Entitled; Polish Naval Academy: Gdynia, Poland, 2015. [Google Scholar]
  10. Wang, C. Simultaneous Localization, Mapping and Moving Object Tracking. Ph.D. Thesis, Robotics Institute, Carnegie Mellon University, Pittsburgh, PA, USA, 2004. [Google Scholar]
  11. Bovio, E.; Cecchi, D.; Baralli, F. Autonomous underwater vehicles for scientific and naval operations. Annu. Rev. Contr. 2006, 30, 117–130. [Google Scholar] [CrossRef]
  12. Daniel, M. SSS-SLAM: An Object Oriented Matlab Framework for Underwater SLAM Using Side Scan Sonar. In Proceedings of the Days Automatic, Automatic Spanish Committee of the IFAC (CEA-IFAC), Valencia, Spain, 3–5 September 2014; pp. 1–8.
  13. Gomez-Gil, J.; Ruiz-Gonzalez, R.; Alonso-Garcia, S.; Gomez-Gil, F.J. A Kalman Filter implementation for precision improvement in low-cost GPS positioning of tractors. Sensors 2013, 13, 15307–15323. [Google Scholar] [CrossRef] [PubMed]
  14. Hasan, A.H.; Grachev, A.N. Target tracking by adaptive EKF using fast genetic algorithm. Int. J. Inf. Eng. 2014, 4, 50–58. [Google Scholar]
  15. Hegrenæs, Ø.; Hallingstad, O.; Gade, K. Towards model-aided navigation of underwater vehicles. Model. Identif. Control 2007, 28, 113–123. [Google Scholar] [CrossRef]
  16. McColgan, J.; McGookin, W.E. Coordination of multiple biomimetic autonomous underwater vehicles using strategies based on the schooling behaviour of fish. Robotics 2016, 5, 1–23. [Google Scholar] [CrossRef]
  17. Muhaiyuddin, M.D.; Hassan, Z.; Jia, C.; Zhou, H.; Pan, R.; Anvar, A.M. Modelling and simulation of the Autonomous Underwater Vehicle (AUV) robot. In Proceedings of the 21st International Congress on Modelling and Simulation, Gold Coast, Australia, 29 November–4 December 2015; pp. 739–745.
  18. Tian, S.; He, B.; Chen, S.; Liu, G. Application of Forward-Looking Scanning Sonar on Simultaneous Localization and Mapping for an Autonomous Underwater Vehicle; IACSIT Press: Singapore, 2012. [Google Scholar]
  19. Wang, S.; Chen, L.; Hu, H.; Gu, D. Single beacon based localization of AUVs using moving horizon. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Tokyo, Japan, 3–7 November 2013; pp. 885–900.
  20. Wang, S.; Chen, L.; Hu, H.; Gu, D. An optimization based moving horizon estimation with application to localization of autonomous underwater vehicles. Robot. Auton. Syst. 2014, 62, 1581–1596. [Google Scholar] [CrossRef]
  21. INS/GPS VectorNav VN-200 User Manual. Available online: http://www.vectornav.com (accessed on 1 April 2016).
  22. Electromagnetic Speed Log Alize Specification. Available online: http://www.benmarine.fr/files/ALIZE_en.pdf (accessed on 1 April 2016).
  23. Pressure Transmitter WIKA S-20. Available online: http://www.przetworniki24.pl/pl/katalog/przetworniki_ciNSienia/S_20 (accessed on 1 April 2016).
  24. USBL MicronNav System. Available online: http://www.tritech.co.uk/media/products/usbl-tracking-system-micronnav.pdf (accessed on 1 April 2016).
  25. Hydromodem Micron Data Modem. Available online: http://www.tritech.co.uk/product/micron-data-modem (accessed on 1 April 2016).
  26. Naus, K. Evaluation of accuracy the position of the vessel designated stereoscopic cameras system. Rep. Geod. 2010, 88, 79–88. [Google Scholar]
  27. Naus, K. Accuracy in fixing ship’s positions by CCD camera survey of horizontal angles. Geomat. Environ. Eng. 2011, 5, 47–61. [Google Scholar]
  28. Naus, K.; Waz, M. Accuracy in fixing ship’s positions by camera survey of bearings. Geod. Cartogr. 2011, 60, 61–73. [Google Scholar] [CrossRef]
  29. Waz, M.; Naus, K. Ripping Radar Image with the Sea Map Using the Alignment Method; Logistyka, Institute of Logistics and Warehousing: Radom, Poland, 2010; pp. 2009–2019. [Google Scholar]
  30. Kiriy, E.; Buehler, M. Three-State Extended Kalman Filter for Mobile Robot Localization; Technical Report TR-CIM 05.07; McGill University: Montreal, QC, Canada, 2002; pp. 1–17. [Google Scholar]
  31. What is SBAS? Available online: http://egnos-portal.gsa.europa.eu/discover-egnos/about-egnos/what-sbas (accessed on 1 April 2016).
Figure 1. Side view of 3D model of the BAUV built in Polish Naval Academy (3D design by Bogdan Szturomski) [9].
Figure 1. Side view of 3D model of the BAUV built in Polish Naval Academy (3D design by Bogdan Szturomski) [9].
Sensors 16 01279 g001
Figure 2. The principle of determining the horizontal coordinates ( x ,   y ) using USBL system.
Figure 2. The principle of determining the horizontal coordinates ( x ,   y ) using USBL system.
Sensors 16 01279 g002
Figure 3. Mapping the accuracy of the coordinate position of determined USBL system (describes mean error M x y in meters).
Figure 3. Mapping the accuracy of the coordinate position of determined USBL system (describes mean error M x y in meters).
Sensors 16 01279 g003
Figure 4. Changing the parameters of the mean error ellipse as a function of measured distance.
Figure 4. Changing the parameters of the mean error ellipse as a function of measured distance.
Sensors 16 01279 g004
Figure 5. Mean error of coordinates as a function of distance covered (adopted P ( k = 0 ) = [ 0 0 0 0 ] , SOG = 1 m/s, COG = 90°).
Figure 5. Mean error of coordinates as a function of distance covered (adopted P ( k = 0 ) = [ 0 0 0 0 ] , SOG = 1 m/s, COG = 90°).
Sensors 16 01279 g005
Figure 6. Change of parameters of mean error ellipse after covering a distance of 200 and 500 m.
Figure 6. Change of parameters of mean error ellipse after covering a distance of 200 and 500 m.
Sensors 16 01279 g006
Figure 7. Measuring testing ground.
Figure 7. Measuring testing ground.
Sensors 16 01279 g007
Figure 8. Determined routes (single passage).
Figure 8. Determined routes (single passage).
Sensors 16 01279 g008
Figure 9. Graph of simulated measurement error Δ C O G .
Figure 9. Graph of simulated measurement error Δ C O G .
Sensors 16 01279 g009
Figure 10. Graph of simulated measurement error Δ S O G .
Figure 10. Graph of simulated measurement error Δ S O G .
Sensors 16 01279 g010
Figure 11. Graph of simulated measurement error Δ α , Δ β .
Figure 11. Graph of simulated measurement error Δ α , Δ β .
Sensors 16 01279 g011
Figure 12. Graph of simulated measurement error Δ r .
Figure 12. Graph of simulated measurement error Δ r .
Sensors 16 01279 g012
Figure 13. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 13. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g013
Figure 18. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Figure 18. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Sensors 16 01279 g018
Figure 19. Determined routes (one hundred passages).
Figure 19. Determined routes (one hundred passages).
Sensors 16 01279 g019
Figure 20. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Figure 20. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Sensors 16 01279 g020
Figure 21. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 21. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g021
Figure 22. Graph distancing the reference position from the position determined by EKF and DR methods of the USBL system (single passage).
Figure 22. Graph distancing the reference position from the position determined by EKF and DR methods of the USBL system (single passage).
Sensors 16 01279 g022
Figure 23. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 23. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g023
Figure 24. Determined route (one hundred passages).
Figure 24. Determined route (one hundred passages).
Sensors 16 01279 g024
Figure 25. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (hundred passages).
Figure 25. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (hundred passages).
Sensors 16 01279 g025
Figure 26. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Figure 26. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Sensors 16 01279 g026
Figure 27. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 27. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g027
Figure 28. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 28. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g028
Figure 29. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Figure 29. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (single passage).
Sensors 16 01279 g029
Figure 30. Determining the route (one hundred passages).
Figure 30. Determining the route (one hundred passages).
Sensors 16 01279 g030
Figure 31. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Figure 31. Graph of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Sensors 16 01279 g031
Figure 32. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Figure 32. Histogram of the reference position distance from the position estimated by EKF and DR methods and determined by the USBL system (one hundred passages).
Sensors 16 01279 g032
Table 1. Statistical parameters of positioning errors by EKF and DR methods and the USBL system (single passage).
Table 1. Statistical parameters of positioning errors by EKF and DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.16.62.2
DR0.725.815.9
USBL0.375.222.6
Table 4. The statistical parameters of the positioning errors using EKF, DR methods and the USBL system (one hundred passages).
Table 4. The statistical parameters of the positioning errors using EKF, DR methods and the USBL system (one hundred passages).
Positioning MethodMinimum Distance from the Reference Position (m)Maximum Distance from the Reference Position (m)Average Distance from the Reference Position (m)
EKF0.63.32.2
DR0.620.110.6
USBL1.357.023.4
Table 5. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Table 5. Statistical parameters of positioning errors using EKF and DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.27.33.4
DR0.734.910.9
USBL0.468.024.3
Table 8. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Table 8. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (single passage).
Positioning MethodMinimum Distance to the Reference Position (m)Maximum Distance to the Reference Position (m)Average Distance to the Reference Position (m)
EKF0.74.42.9
DR0.719.110.8
USBL1.355.922.7
Table 12. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (one hundred passages).
Table 12. Statistical parameters of positioning errors using EKF, DR methods and the USBL system (one hundred passages).
Positioning MethodMinimum Distance from the Reference Position (m)Maximum Distance from the Reference Position (m)Average Distance from the Reference Position (m)
EKF0.66.44.4
DR0.618.910.7
USBL4.843.020.5

Share and Cite

MDPI and ACS Style

Krzysztof, N.; Aleksander, N. The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements. Sensors 2016, 16, 1279. https://doi.org/10.3390/s16081279

AMA Style

Krzysztof N, Aleksander N. The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements. Sensors. 2016; 16(8):1279. https://doi.org/10.3390/s16081279

Chicago/Turabian Style

Krzysztof, Naus, and Nowak Aleksander. 2016. "The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements" Sensors 16, no. 8: 1279. https://doi.org/10.3390/s16081279

APA Style

Krzysztof, N., & Aleksander, N. (2016). The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements. Sensors, 16(8), 1279. https://doi.org/10.3390/s16081279

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