Open AccessTechnical Note

Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements

Asaf Tal

^1,2,*,

Itzik Klein

² and

Reuven Katz

Faculty of Mechanical Engineering, Technion-Israel Institute of Technology, Haifa 3200003, Israel

Rafael Advanced Defense Systems Ltd, Haifa 3102102, Israel

Author to whom correspondence should be addressed.

Sensors 2017, 17(2), 415; https://doi.org/10.3390/s17020415

Submission received: 29 November 2016 / Revised: 26 January 2017 / Accepted: 14 February 2017 / Published: 22 February 2017

(This article belongs to the Special Issue Inertial Sensors and Systems 2016)

Download

Browse Figures

Figure 1
Scheme of loosely coupled and tightly coupled approaches for INS/DVL fusion. AUV: Autonomous underwater vehicle; DVL: Doppler velocity log; INS: inertial navigation system. "> Figure 2
Technion autonomous underwater vehicle six degrees of freedom (TAUV 6 DOF) simulation layout. IMU: inertial measurement unit. "> Figure 3
Top level diagram of the TAUV navigation system. EKF: extended Kalman filter. GPS: global positioning system. "> Figure 4
Platform frame and tangent frame in respect to the TAUV. "> Figure 5
DVL orientation setup in the TAUV platform frame. "> Figure 6
Extended loosely coupled (ELC) shame for utilizing partial measurements from the DVL. "> Figure 7
Block diagram for implementation of ELC approach in the TAUV navigation system. "> Figure 8
Simulation trajectories used to examine the ELC approaches. "> Figure 9
Root mean square (RMS) errors of velocity for all approaches obtained from trajectory #1. TC: tightly coupled. "> Figure 10
RMS errors of attitude for all approaches obtained from trajectory #1. "> Figure 11
RMS errors of velocity for all approaches obtained from trajectory #2. "> Figure 12
RMS errors of attitude for all approaches obtained from trajectory #2. "> Figure 13
RMS errors of velocity for all approaches obtained from trajectory #3. "> Figure 14
RMS errors of attitude for all approaches obtained from trajectory #3. ">

Versions Notes

Abstract

The Technion autonomous underwater vehicle (TAUV) is an ongoing project aiming to develop and produce a small AUV to carry on research missions, including payload dropping, and to demonstrate acoustic communication. Its navigation system is based on an inertial navigation system (INS) aided by a Doppler velocity log (DVL), magnetometer, and pressure sensor (PS). In many INSs, such as the one used in TAUV, only the velocity vector (provided by the DVL) can be used for aiding the INS, i.e., enabling only a loosely coupled integration approach. In cases of partial DVL measurements, such as failure to maintain bottom lock, the DVL cannot estimate the vehicle velocity. Thus, in partial DVL situations no velocity data can be integrated into the TAUV INS, and as a result its navigation solution will drift in time. To circumvent that problem, we propose a DVL-based vehicle velocity solution using the measured partial raw data of the DVL and additional information, thereby deriving an extended loosely coupled (ELC) approach. The implementation of the ELC approach requires only software modification. In addition, we present the TAUV six degrees of freedom (6DOF) simulation that includes all functional subsystems. Using this simulation, the proposed approach is evaluated and the benefit of using it is shown.

Keywords:

Doppler velocity log/inertial navigation system (DVL/INS) fusion; loosely coupled; tightly coupled; autonomous underwater vehicle six degrees of freedom (AUV 6DOF) simulation; partial measurement

1. Introduction

Underwater operations attract great attention for environmental matters and potential resources, as well as scientific interest. Therefore, the need for underwater robotic systems has become more apparent [1]. With the emergence of inspection-class autonomous underwater vehicles (AUVs), navigation and navigational accuracy are becoming increasingly important in order for the AUV to complete its task. Without an operator in the loop, the vehicle must use sensors to determine its position, velocity and orientation. Most AUVs employ an inertial navigation system (INS) as their main navigation sensor [2,3,4]. This is for many reasons; one of which is that the INS is a standalone system that can provide all of the required navigation data: position, velocity and orientation. However, even with a high-grade INS, the navigation solution drifts in time due to measurement errors of its inertial sensors. Therefore, INSs are usually aided by other external sensors or data such as the Doppler velocity log (DVL) for velocity, magnetometers for heading and depth/pressure sensor for altitude [5,6]. For positioning, global positioning system (GPS) is commonly used in land or air vehicles, yet it cannot be used underwater [7,8] since the electromagnetic signals decay very quickly in the water. A common method for overcoming the localization problem of underwater vehicles is to use acoustic positioning methods such as long baseline (LBL) or ultra-short baseline (USBL) [9]. One drawback of such approaches is the need for a priori deployment of beacons. An alternative method for obtaining localization of underwater vehicles is geophysical navigation which uses physical features of the AUV’s environment to produce an estimate of the location of the AUV [10,11].

A team at the Technion–Israel Institute of Technology is currently developing an AUV named the Technion autonomous underwater vehicle (TAUV). The TAUV project goal is to develop and to produce a small autonomous underwater vehicle, which will serve as a technology demonstrator and a platform for various research programs. The TAUV diameter is 300 mm and its length is 3000 mm. The TAUV has four rudders installed in an “x” configuration, and its maximum operation depth is 200 m. The TAUV navigation system based on an INS aided by a DVL, magnetometer, pressure sensor (PS), and GPS when available [12]. When submerged, the main aiding sensor is the DVL which provides velocity measurements that can help estimate, in addition to the vehicle velocity, some of the orientation and inertial sensor error states depending on the vehicle dynamics [13,14].

There are two main approaches for sensor fusion between INS and other sensors: (1) loosely coupled (LC); and (2) tightly coupled (TC) [15]. We shall refer here to INS/DVL fusion, but the underlying principles are the same for other sensors as well. A top-level block diagram of LC and TC approaches is presented in Figure 1. The raw data of the DVL is the relative velocity in each beam direction. In the LC approach, using parameter estimation the DVL raw data is used to calculate the vehicle velocity, which in turn is compared with its INS counterpart in the navigation filter. The advantage of this method is the simplicity of integration and the ability to combine any off-the-shelf INS with any DVL. However, in order for the DVL to calculate vehicle velocity, it must operate in bottom lock, which refers to the condition when a sufficient number of beam measurements (at least three) are available. In the TC approach, the DVL raw data is directly used in the navigation filter. That is, each beam measurement of the DVL is compared with its calculated INS counterpart and independently integrated into the navigation filter. Therefore, there is no need for a bottom lock stage, and aiding may be applied even with a single beam measurement.

The selected TAUV INS can only receive from the DVL an external velocity vector measurement [16]. Therefore, INS/DVL fusion with the TC approach is not possible, and thus the LC approach must be used. In some situations—such as operation in close proximity to the seafloor, or when experiencing extreme tilt angle or beam malfunction—one or more of the DVL beams may not provide the reflection required for determining velocity [17]. In such cases of partial beam measurements, the DVL fails to maintain bottom lock and is unable to estimate the AUV velocity. Thus, the TAUV INS navigation solution will drift in time.

In this paper, we propose the extended loosely coupled (ELC) approach for calculating vehicle velocity when only partial DVL measurements are available. This is made possible by using the partial measured raw data from the DVL combined with additional information. This calculated vehicle velocity is fed back to the TAUV INS and used for aiding the INS, as in the regular LC approach.

Related approaches in similar situations where partial measurements are available can be found in the literature for INS/GPS fusion or LBL positioning. In [18], a single beacon measurement and two virtual ranging measurements were used to triangulate and solve the vehicle position with the LBL method. Another method is to utilize the knowledge of the vehicle dynamics under a specific scenario. For example, in [19] the authors utilized the vehicle dynamics in a straight trajectory scenario to improve navigation performance with low-cost GPS. In [20] the construction of virtual GPS satellites was suggested in order to facilitate GPS receiver position and velocity solutions in cases of partial GPS availability, thereby enabling the fusion of GPS/INS in the LC approach.

In order to examine and demonstrate the improvement of the navigation performance under the ELC approach, we developed the TAUV six degrees of freedom (6DOF) simulation, which includes all functional sub-systems. This simulation consists of the AUV guidance, navigation and control subsystems and uses a complete hydrodynamic model. Simulation results show great improvement in the navigation performance using the proposed approach compared to the standalone TAUV INS solution.

The rest of the paper is organized as follows: Section 2 describes the TAUV 6DOF simulation and in particular its navigation system. In Section 3 the ELC approach is derived. In Section 4 6DOF simulation results are presented to demonstrate the effectiveness of the partial measurements approach on the TAUV navigation system. Finally, Section 5 provides the conclusions.

2. TAUV 6DOF Simulation

A top-level block diagram of the TAUV 6DOF simulation with its main subsystems is presented in Figure 2. The purpose of each subsystem is discussed next.

Hydrodynamic and Motion block: This model calculates the forces and moments applied on the TAUV for the 6DOF rigid body motion model. The motion model solves the vehicle kinematics to produce the true values of the vehicle’s position, velocity, attitude, and acceleration. The model contains a complete hydrodynamic model derived for the TAUV.

Control system block: The control system is used to stabilize the vehicle attitude by determining the angles of the four TAUV servos. The vehicle velocity is stabilized by controlling the thruster rotation speed. The TAUV control system contains four proportional integrative derivative (PID) controllers: three for the attitude control and one for the velocity control.

Guidance system block: The guidance laws in the TAUV are based on target tracking methods [21]. This subsystem determines the vehicle steering command and the desired velocity in order to follow the AUV’s desired trajectory.

Sensors generator block: The block calculates sensors’ outputs, including their errors, by statistical means. The sensors are: accelerometers, gyros, DVL, PS, magnetometer, and GPS.

Navigation block: The navigation block represents the navigation system of the TAUV, including the navigation equations and navigation filter. The TAUV navigation system is presented in detail in the following section.

2.1. TAUV Navigation System

In Figure 3, a top-level block diagram of the TAUV navigation system is shown. The accelerometer and gyro measurements are integrated into the INS system in order to produce its standalone solution for position, velocity, and attitude. The extended Kalman filter (EKF) [22,23,24] fuses the INS with measurements from the other sensors: the magnetometer, PS, DVL and GPS. By comparing the INS state calculation and aiding the sensor measurements, the EKF estimate corrects the error state vector comprising the INS states and inertial sensor error terms. Notice that each aiding sensor operates at a different sampling rate, yet in each time instance that a measurement arrives, it is fused with the INS.

The navigation equations and the EKF are expressed in a coordinate frame, known as the platform frame. The platform frame is located at the center of buoyancy of the vehicle. The navigation frame is determined by the INS frame. For simplicity of the dynamics equations, we assume that the platform frame and the INS frame are located at the same point at the AUV and with identical orientation. The TAUV platform frame and the tangent frame are represented in Figure 4. The vehicle velocities

(u, v, w)

are expressed in the platform frame, and the vehicle attitudes expressed in the tangent frame by the roll, pitch, and yaw

(ϕ, θ, ψ)

Euler angles.

2.1.1. Navigation Filter

An error state EKF is implemented in the TAUV navigation system. The error state vector

δ x

is defined by the difference between the true state vector and Kalman filter estimate and expressed by [17]:

δ x_{k} = x_{k} - {\hat{x}}_{k}

(1)

where k is the time step. The standard INS error state vector is in accordance with [22]:

δ x = {[\begin{matrix} {(δ r_{t / p}^{t})}^{T} & δ ρ^{T} & {(δ v_{t / p}^{p})}^{T} & δ b_{a}^{T} & δ b_{g}^{T} \end{matrix}]}^{T} \in ℝ^{15 \times 1}

(2)

where

δ r_{t / p}^{t}

is the position error vector expressed in the tangent frame; the vector

δ ρ

is the orientation error also referred to as the frame misalignment error;

δ v_{t / p}^{p}

is the velocity error vector expressed in the platform frame;

δ b_{a}

is the accelerometer biases error expressed in the platform frame; and

δ b_{g}

is the gyro biases error expressed in the platform frame. The INS error state model has the following form:

δ \dot{x} (t) = F_{I N S} (t) δ x + G_{I N S} (t) w_{I N S} (t)

(3)

where

F_{I N S} (t)

is the matrix that relates the error state to the dynamic state [17]:

F_{I N S} = [\begin{matrix} 0_{3 \times 3} & - [{\hat{R}}_{p}^{t} {\hat{v}}_{t / p}^{p} \times] & {\hat{R}}_{p}^{t} & 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & - {\hat{Ω}}_{i / e}^{t} & 0_{3 \times 3} & 0_{3 \times 3} & - {\hat{R}}_{p}^{t} \\ 0_{3 \times 3} & F_{32} & F_{33} & - I_{3} & - [{\hat{v}}_{t / p}^{p} \times] \\ 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} & 0_{3 \times 3} \end{matrix}] \in ℝ^{15 \times 15}

(4)

{\hat{R}}_{p}^{t}

is the estimated transformation matrix from platform frame to tangent frame, and

{\hat{Ω}}_{i / e}^{t}

is a skew symmetric matrix of the earth’s rotation rate expressed in the tangent frame. The expressions for

F_{32}

F_{33}

can be found for example in [17]. The matrix

G_{I N S} (t)

relating the error state to the noise process is defined by:

G_{I N S} (t) = [\begin{matrix} \begin{matrix} 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & - {\hat{R}}_{p}^{t} \\ - I_{3} & - [{\hat{v}}_{t / p}^{p} \times] \\ 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & 0_{3 \times 3} \end{matrix} & \begin{matrix} 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & 0_{3 \times 3} \\ 0_{3 \times 3} & 0_{3 \times 3} \\ I_{3} & 0_{3 \times 3} \\ 0_{3 \times 3} & I_{3} \end{matrix} \end{matrix}] \in ℝ^{15 \times 12}

(5)

The process noises vector is defined by:

w_{I N S} (t) = {[\begin{matrix} \begin{matrix} w_{a}^{T} & w_{g}^{T} \end{matrix} & \begin{matrix} w_{b a}^{T} & w_{b g}^{T} \end{matrix} \end{matrix}]}^{T} \in ℝ^{12 \times 1}

(6)

where

w_{a}

w_{g}

are the accelerometers and gyros noise vectors are modeled as zero-mean white Gaussian processes,

w_{b a}

is the noise vector of the accelerometer biases modeled as random walk and

w_{b g}

is the noise vector of the gyro biases modeled as random walk.

In the linear case, the measurement model is:

δ y = H \times δ x_{k} + n_{w}

(7)

where

n_{w}

is zero-mean Gaussian white noise. The matrix

H

relates the measurement to the error state and is known as the measurement matrix. For each aiding sensor and for each measurement, a measurement rejection algorithm is applied, by examining if the data is within range of three standard deviations compared to the INS estimation.

In the presented navigation system model, we augment the INS state vector from Equation (2) with five additional error states of the DVL: four to model the DVL biases of each beam, and one for common DVL scale factor error. Thus, the augmented error state vector:

δ x = {[\begin{matrix} \begin{matrix} {(δ r_{t / p}^{t})}^{T} & δ ρ^{T} & {(δ v_{t / p}^{p})}^{T} & δ b_{a}^{T} & δ b_{g}^{T} & δ b_{D V L}^{T} \end{matrix} & δ S F_{D V L} \end{matrix}]}^{T} \in ℝ^{20 \times 1}

(8)

where

δ b_{D V L}

is the DVL biases error expressed in the DVL frame and

S F_{D V L}

is the DVL scale factor error. The corresponding dynamics matrix F is given by:

F (t) = [\begin{matrix} F_{I N S} & 0_{15 \times 5} \\ 0_{5 \times 15} & 0_{5 \times 5} \end{matrix}] \in ℝ^{20 \times 20}

(9)

The corresponding noise process matrix G is expressed via:

G (t) = [\begin{matrix} G_{I N S} & 0_{15 \times 5} \\ 0_{5 \times 12} & I_{5} \end{matrix}] \in ℝ^{20 \times 17}

(10)

where the filter process noises vector is modified via:

w (t) = {[\begin{matrix} \begin{matrix} \begin{matrix} w_{a}^{T} & w_{g}^{T} \end{matrix} & \begin{matrix} w_{b a}^{T} & w_{b g}^{T} \end{matrix} \end{matrix} & \begin{matrix} w_{b D V L}^{T} & w_{s f D V L}^{T} \end{matrix} \end{matrix}]}^{T} \in ℝ^{17 \times 1}

(11)

where

w_{b D V L}

is a noise vector modeled as random walk of the DVL biases, and

w_{s f D V L}

is zero-mean Gaussian white noise of the DVL scale factor.

2.1.2. Doppler Solution

The DVL measures the Doppler frequency shift for each depth cell and each beam and then computes the component of relative flow velocity in the direction of each acoustic beam. Based on [25], the relative velocity of each beam can be assumed as:

V_{r e l} = (F_{D} + b_{F_{D}} + n_{F_{D}}) \frac{C (1 + S F_{c})}{2 F_{s}} 1000

(12)

where

F_{D}

is the Doppler frequency shift,

b_{F_{D}}

is the frequency shift bias,

n_{F_{D}}

is the frequency shift noise, C is the velocity of sound in water at transducer face,

F_{s}

is the transmitted acoustic frequency,

S F_{c}

is the scale factor error of the velocity of sound in water, and

n_{F_{D}}

is zero-mean Gaussian white noise. The scale factor is due to the variation of the velocity of sound under different temperatures and salinity level of the water. From Equation (12) we can observe that the DVL measurement model has three main errors: (1) bias error (2) scale factor error and (3) white noise.

The DVL assembly in the TAUV is constructed in an “x” configuration relative to the platform frame (“Janus Doppler configuration”), as shown in Figure 5. Generally, each DVL beam direction in the vehicle frame is defined by:

b_{i} = [\begin{matrix} c {\tilde{ψ}}_{i} s \tilde{θ} \\ s {\tilde{ψ}}_{i} s \tilde{θ} \\ c \tilde{θ} \end{matrix}]

(13)

where

b_{i}

represents the direction of each transducer/beam of the DVL expressed in the platform frame, for i = 1,2,3,4,

\tilde{θ}

is the pitch angle (relative to platform frame) of each DVL beam, and

{\tilde{ψ}}_{i}

is the yaw angle (relative to platform frame) of each DVL beam for i = 1,2,3,4. In the TAUV DVL,

\tilde{θ}

is equal to 20°.

The expression for yaw angle of each beam direction in the TAUV can be written as:

{\tilde{ψ}}_{i} = (i - 1) 90 ° + 45 °, i = 1, 2, 3, 4

(14)

The DVL beam velocity is modeled as:

{\tilde{y}}_{v_{i}} = [v_{t / p}^{p} (1 + S F_{D V L}) + ω_{t / p}^{p} \times l_{D_{i}}] \times b_{i} + n_{v_{i}} + b_{D V L, i}

(15)

where

l_{D_{i}}

is the position vector of each transducer expressed in the platform frame,

ω_{t / p}^{p}

is the angular rate vector of the platform expressed in the platform frame,

b_{D V L}

is a 4 × 1 vector representing the bias of each beam,

S F_{D V L}

is the scale factor error, and

n_{v_{i}}

is measurement noise. Let A be a matrix contacting all beam direction vectors:

A = [\begin{matrix} b_{1}^{T} \\ b_{2}^{T} \\ b_{3}^{T} \\ b_{4}^{T} \end{matrix}]

(16)

Thus, by using Equation (15), the DVL-based velocity measurement of the vehicle can be expressed by:

{\tilde{v}}_{t / p}^{p} = {\tilde{y}}_{v_{i}} {(A^{T} A)}^{- 1} A^{T}

(17)

Note that the velocity solution in Equation (17) is based on four beams, of which one beam is for redundancy due to the beam configuration. Therefore, in practice only three beams are enough to calculate the vehicle velocity.

2.1.3. DVL Fusion

In the LC fusion technique, the DVL output to the navigation filter is the platform velocity (Equation (17)). The corresponding measurement matrix is:

H = [\begin{matrix} 0 & 0 & I & 0 & 0 & 0 & 0 \end{matrix}] \in ℝ^{20 \times 3}

(18)

and the velocity measurement noise covariance matrix is:

R_{v e l} = {(A^{T} W A)}^{- 1} \in ℝ^{3 \times 3}

(19)

where

W = R_{D V L}^{- 1}

and

R_{D V L}

is the DVL measurement noise covariance matrix [17], defined as:

R_{D V L} = b_{i}^{T} [l_{D_{i}} \times] σ_{g}^{2} I {[l_{D_{i}} \times]}^{T} b_{i} + σ_{v}^{2} \in ℝ^{4 \times 4}

(20)

In the case of small-level arm

l_{D_{i}}

, Equation (20) can be well approximated as:

R_{D V L} \approx d i a g (σ_{v 1}^{2}, σ_{v 2}^{2}, σ_{v 3}^{2}, σ_{v 4}^{2}) \in ℝ^{4 \times 4}

(21)

where

σ_{v_{i}}

is the velocity variance for each beam measurement for i = 1,2,3,4. The variance value of each beam is taken from the DVL spec [26] or calculated in a lab test.

In the TC fusion technique, the DVL raw data is used to aid the INS. That is, each transducer measurement is fused with its INS counterpart (no need to calculate the DVL vehicle velocity, Equation (17)) and the measurement matrix [17] is defined as:

H_{T C} = [\begin{matrix} \begin{matrix} 0 & 0 & A (1 + S \hat{F}) & 0 & A [l_{D_{i}} \times] \end{matrix} & I & A \times {\hat{v}}_{t / p}^{p} \end{matrix}] \in ℝ^{20 \times 4}

(22)

The measurement noise covariance matrix is taken from Equation (21) for each beam.

3. INS/DVL Fusion with the ELC Approach

As mentioned in Section 1, in some situations several or all of the DVL beams do not provide reflection. With no reflection, the beam cannot provide the relative velocity of the AUV Equation (15). In cases where only one beam is not available, the DVL can still calculate the vehicle velocity Equation (17) but without any redundant data. If two beams are not available, the DVL cannot calculate the vehicle velocity. In this case, no velocity aiding is provided to the TAUV INS.

In this section, we present the ELC approach, with four different methods for utilizing the partial measurements from the DVL (instead of not using them at all) together with external information to construct the DVL based velocity estimate. This velocity is used in the TAUV navigation filter, as illustrated in Figure 6. The four ELC methods differ in the external information they employ as elaborated in the following subsections.

3.1. Virtual Beam

This method utilizes the raw data from the DVL and the filter prediction of the velocity error in order to solve for the vehicle velocity. Without the loss of generality, we assume that beams #3 and #4 are not available. Let the matrix

A_{s}

relate the velocity in the platform frame to the velocity in each beam direction. So, for small value of lever arm, following Equation (17) the velocities in beam directions (#1, #2, and #3) are:

y_{v} = A_{s} \times v_{t / p}^{p} A_{s} = {[\begin{matrix} b_{1}^{T} & b_{2}^{T} & b_{3}^{T} \end{matrix}]}^{T}

(23)

From Equation (23), the value of the third beam may be found as:

y_{v 3} = b_{3}^{T} [v_{t / p}^{p}] y_{v 3} \approx b_{3}^{T} {[\begin{matrix} {\hat{u}}_{k - 1} & {\hat{v}}_{k - 1} & {\hat{w}}_{k - 1} \end{matrix}]}^{T}

(24)

where

v_{t / p}^{p}

is taken from the last INS step. In this method we estimate only the value of the third beam. Generally, the value of the fourth beam can be estimated instead. Plugging Equation (24) into Equation (23) we solve the unknown vehicle velocity vector at time k via:

[\begin{matrix} u_{k} \\ v_{k} \\ w_{k} \end{matrix}] = {[A]}^{- 1} [\begin{matrix} y_{v 1} \\ y_{v 2} \\ y_{v 3} \end{matrix}]

(25)

It is important to note that now, theoretically, a correlation between the measurement noise and process noise exists, since the DVL velocity measurement Equation (25) depends on the DVL partial raw data (with measurement noise) but also on the INS velocity solution Equation (24) (with process noise). Therefore:

E 〈 n_{w}, {\bar{w}}^{T} 〉 \neq 0

(26)

Thus, the covariance term Equation (26), which is usually nullified in the EKF formulation, needs to be included. However, in our current analysis we neglect this term.

In order to calculate the corresponding measurement covariance matrix, we need to derive an expression for the constructed third beam measurement standard deviation (STD)

σ_{v 3}

. To that end we use the error covariance matrix

p_{k - 1}

(last step of the fusing process), to define the velocity STD by:

σ_{u} = \sqrt{P_{7, 7}}; σ_{v} = \sqrt{P_{8, 8}}; σ_{w} = \sqrt{P_{9, 9}}

(27)

Using Equation (24) we derive the expression for

σ_{v 3}

σ_{v 3} = n_{f} \sqrt{{(b_{3, 1} σ_{u})}^{2} + {(b_{3, 2} σ_{v})}^{2} + {(b_{3, 3} σ_{w})}^{2}}

(28)

where

b_{i, j} = b_{i} (j)

and

n_{f}

is a factor to compensate for our assumption of neglecting Equation (26). Using the known values of for

σ_{v 1}

and

σ_{v 2}

and Equation (28), the DVL measurement noise covariance matrix is:

R_{V B 3 D} = [\begin{matrix} σ_{v 1}^{2} & 0 & 0 \\ 0 & σ_{v 2}^{2} & 0 \\ 0 & 0 & σ_{v 3}^{2} \end{matrix}] W_{s} = R_{V B 3 D}^{- 1}

(29)

From Equation (29) the velocity measurement noise covariance matrix is expressed via:

R_{V B} = {({A_{s}}^{T} W_{s} A_{s})}^{- 1}

(30)

In summary, in the proposed approach we calculate one velocity beam of the DVL via Equation (24). Combining the two beam measurement from the DVL, we solve the velocity vector of the vehicle using Equation (25).

3.2. Nullfying Sway Velocity

In this method we set the sway velocity to be zero, that is

v_{k} = 0

. This assumption is reasonable for scenarios such as straight line trajectories, which in practice are the AUV trajectories for most of the operating time (although the AUV may be influenced by some disturbances that alter the straight line). Since we assume v = 0, only the other two velocity components (u, w) are to be calculated using the partial DVL data. Similar to the previous method outlined in Section 3.1, we assume that beams #3 and #4 are not available.

Let:

A_{0} = [\begin{matrix} b_{1, 1} & b_{1, 3} \\ b_{2, 1} & b_{2, 3} \end{matrix}]

(31)

This matrix relates the beam velocities #1 and #2 (since beams #3 and #4 are not available) to the vehicle velocity (u and w). Thus, the relation between DVL beam velocities

y_{v 1}

and

y_{v 2}

to vehicle velocities is:

[\begin{matrix} y_{v 1} \\ y_{v 2} \end{matrix}] = A_{0} [\begin{matrix} u \\ w \end{matrix}]

(32)

From Equation (32) the surge and heave velocities are:

[\begin{matrix} u_{k} \\ w_{k} \end{matrix}] = [A_{0}^{- 1}] [\begin{matrix} y_{v 1} \\ y_{v 2} \end{matrix}]

(33)

Note that the matrix

A_{0}

is not singular, because of the independence of the beams direction vectors. The DVL measurement noise covariance matrix is then:

R_{N S V 2 D} = [\begin{matrix} σ_{v 1}^{2} & 0 \\ 0 & σ_{v 2}^{2} \end{matrix}] W_{0} = R_{N S V 2 D}^{- 1}

(34)

and

R_{0} = {(A_{z e r o}^{T} W_{0} A_{z e r o})}^{- 1}

(35)

From Equation (35) we extract the velocity variances of the calculated vehicle velocity components (u and w). Therefore, the velocity measurement noise covariance matrix that is plugged into the navigation filter is:

R_{N S V} = [\begin{matrix} R_{0} [1, 1] & 0 & 0 \\ 0 & ε & 0 \\ 0 & 0 & R_{0} [2, 2] \end{matrix}]

(36)

where

ε

is a parameter that reflects the assumption of v = 0 and has a value close to zero. Equation (36) is used for practical considerations since the TAUV INS requires a complete velocity vector with its corresponding measurement covariance.

In summary, in the proposed approach we assume the sway velocity component is zero. Combining this assumption with the two beams measurements from the DVL, we solve the velocity vector of the vehicle using Equation (33).

3.3. Partial LC Fusing

This method utilizes the DVL setup configuration in order to calculate one component of the AUV velocity, u or v, depending on the active transducer order. Therefore, in this approach only one velocity component measurement is fused into the filter. As in the previous approaches, we assume that beams #3 and #4 are not available. When the DVL setup is in an “x” configuration, as in the TAUV, we can derive the following relations between the components of the DVL beam direction vectors based on Equation (13):

\begin{array}{l} b_{1, 1} = - b_{2, 1} \\ b_{1, 2} = b_{2, 2} \\ b_{1, 3} = b_{2, 3} \end{array}

(37)

If the DVL setup is in a “+” configuration and the available beams are pointing in the same direction, the two components of the velocity vector—u and w or v and w—can be calculated. If the active beams are not pointing in the same direction (e.g., beam #1 is aimed toward the surge direction and beam #2 toward the sway direction) we cannot utilize any data from the DVL using this method. Here, (under the assumption that only beams #1 and #2 are available), we intend to solve the surge velocity u, utilizing the partial measurements from the DVL. Using Equation (37) we can express the velocity of beams #1 and #2 as:

\begin{array}{l} b_{1, 1} u + b_{1, 2} v + b_{1, 3} w = y_{v 1} \\ - b_{1, 1} u + b_{1, 2} v + b_{1, 3} w = y_{v 2} \end{array}

(38)

Subtraction gives:

2 b_{1, 1} u = y_{v 1} - y_{v 2}

(39)

The sway velocity is found using Equation (39):

u_{k} = \frac{y_{v 1} - y_{v 2}}{2 b_{1, 1}}

(40)

where the corresponding measurement noise covariance is found using Equation (40):

R_{u, P L C F} = \frac{σ_{v 1}^{2} + σ_{v 2}^{2}}{4 b_{1, 1}^{2}}

(41)

For practical considerations in the TAUV INS, the velocity measurement noise covariance matrix must have a structure of

R_{P L C F} \in ℝ^{3 \times 3}

. So, the matrix that is passed to the filter is defined by:

R_{P L C F} = [\begin{matrix} R_{u, P L C F} & 0 & 0 \\ 0 & R_{v, P L C F} & 0 \\ 0 & 0 & R_{w, P L C F} \end{matrix}]

(42)

Yet, since only the surge velocity, u, is measured, variances

R_{v, P L C F}

and

R_{w, P L C F}

, corresponding to sway and heave velocities, are:

R_{v, P L C F} \to \infty, R_{w, P L C F} \to \infty

(43)

To summarize, in the proposed approach we calculate the surge velocity component using the two beam measurements Equation (40). The filter ignores the two other velocity components (sway and heave) by setting their measurement variance to infinity. That is, in practice only one velocity component is used to aid the INS.

3.4. Virtual Heave Velocity

This method is an elaboration of the previous method derived in Section 3.3, where we extract the surge velocity u from the DVL raw data, while the DVL allows us access to two measured velocities. In order to utilize the two DVL measurements, we use the velocity prediction from the filter. Therefore, in this approach two vehicle velocity components are used to aid the INS: u is taken from Section 3.3, and v is calculated here. To that end, we sum the two equations from Equation (38) to derive the following expression:

2 b_{1, 2} v + 2 b_{1, 3} w = y_{v 1} + y_{v 2}

(44)

Thus, we have one equation, Equation (44), with two unknown velocity components,

v

and

w

. To calculate them we use the estimated velocity from the navigation filter. Assuming the AUV is travelling in a straight line, the heave velocity (

w

) estimation is probably more accurate due to the gravitation vector in heave (and even more so if the pressure sensor is active to measure the vehicle depth). Using Equation (44) and

w_{k - 1}

from the filter, the sway velocity is:

v_{k} = \frac{y_{v 1} + y_{v 2}}{2 b_{1, 2}} - \frac{b_{1, 3}}{b_{1, 2}} {\hat{w}}_{k - 1}

(45)

The corresponding velocity variance component is:

R_{v, V H V} = \frac{σ_{v 1}^{2} + σ_{v 2}^{2}}{4 b_{1, 2}^{2}} - \frac{b_{1, 3}^{2}}{b_{1, 2}^{2}} {σ_{w}}^{2}

(46)

From Equations (41) and (46) the velocity measurement noise covariance matrix that is plugged into the navigation filter is given in Equation (47). As before, for practical considerations in the TAUV INS, the velocity measurement noise covariance matrix must have a structure of

R_{V H V} \in ℝ^{3 \times 3}

; thus we set

R_{w, V H V} \to \infty

, and the measurement noise covariance matrix is:

R_{V H V} = [\begin{matrix} R_{u, P L C F} & 0 & 0 \\ 0 & R_{v, V H V} & 0 \\ 0 & 0 & R_{w, V H V} \end{matrix}]

(47)

where

R_{u, P L C F}

is defined in Equation (41) and

R_{v, V H V}

is defined in Equation (46).

In summary, in the proposed approach we use the surge velocity component calculated as in Section 3.3. Using the estimated heave velocity component, the sway velocity component is solved via Equation (45). The filter ignores the other velocity component (heave) by setting its measurement variance to infinity. That is, in practice only two velocity components are used to aid the INS.

3.5. Summary of ELC Methods

Table 1 presents a summary of the four ELC methods. VB: virtual beam; NSV: nullifying sway velocity; PLCF: partial loosely coupled fusing; VHV: virtual heave velocity

3.6. ELC Implemenation

The DVL enables access to its raw data and its calculated velocity. In the case of partial measurements, the DVL cannot calculate the vehicle velocity, and only the partial raw data is available. From the raw data, we can calculate the vehicle velocity and the variance from all four approaches as presented in the previous sections. Instead of choosing which approach to use in the TAUV, in Equation (48) we employ a simple selector to pick the best vehicle velocity components out of the four approaches and use each derived corresponding velocity to aid the INS.

\begin{array}{l} R_{u} = \min {R_{u, VB}, R_{u, N S V}, R_{u, P L C F}, R_{u, V H V}} \\ R_{v} = \min {R_{v, VB}, R_{v, N S V}, R_{v, P L C F}, R_{v, V H V}} \\ R_{w} = \min {R_{w, VB}, R_{w, N S V}, R_{w, P L C F}, R_{w, V H V}} \end{array}

(48)

The velocity measurement noise covariance matrix (with the corresponding velocity measurement) that goes to the TAUV INS is:

R_{v e l} = [\begin{matrix} R_{u} & 0 & 0 \\ 0 & R_{v} & 0 \\ 0 & 0 & R_{w} \end{matrix}]

(49)

Figure 7 shows the implementation of the ELC approach and LC approach in the TAUV. The red lines are relevant only in partial measurement scenarios.

4. Analysis and Results

4.1. Simulation Scenarios and Parameters

In this section, we evaluate the ELC approach under three commonly used trajectories of the TAUV and other AUVs. The objective of using several trajectories is to evaluate the navigation performance of each method under different dynamic conditions. Figure 8 shows the chosen trajectories for the analysis. The first trajectory is a straight-line trajectory along the north direction while the AUV travels with constant speed and depth. This trajectory is chosen since in steady state the TAUV will, theoretically, travel in straight lines. The second trajectory follows a figure eight pattern. This trajectory contains permanent maneuvers and also a dive with constant velocity. The third trajectory is a classic seeking trajectory (lawn mower pattern) used for bottom survey. This trajectory begins with a dive at constant rate, and then in the searching phase constant depth is kept.

The simulation and navigation parameters used to produce the results in the paper are summarized in Table 2. The navigation parameters were chosen according to a literature survey [27] and the actual TAUV sensors specifications [16,26].

4.2. Simulation Results

The navigation performance obtained using Monte Carlo runs with the 6DOF simulation is presented in terms of the root mean square (RMS) errors of the velocity vector and the attitude error under the scenarios of trajectories #1, #2 and #3. The accelerometer and gyro biases are not presented since they behave as expected when fusing any other velocity measurements. That is, when the vehicle travels with constant velocity (trajectory #1) only the z-axis accelerometer bias and x and y gyro bias are observable. When maneuvering (trajectories #2 and #3) the estimation performance improves in terms of more observable error states.

We compare the performance of the ELC to the performance of the standalone INS (which is the case in the TAUV with partial DVL measurements scenario) and TC fusion, in order to test the effectiveness of the approach. We note that the analysis was made only for DVL/INS fusion. We expect that when using other available sensors on the TAUV (such as PS or magnetometer), the performance of the proposed approach will improve but this is not the subject of the paper.

In the next figures, the following notations are used: TC (tightly coupled; Section 2.1.3), VB (virtual beam; Section 3.1), NSV (nullifying sway velocity; Section 3.2), PLCF (partial loosely coupled fusing; Section 3.3), VHV (virtual heave velocity Section 3.4).

4.2.1. Trajectory #1 Simulation Results

The velocity RMS errors for all the methods (including TC) are presented in Figure 9. Not presented in the figure is the standalone INS performance, which is 34 m/s of error after 250 s.

All the proposed approaches improved the standalone INS performance; in particular, the NSV reduced the RMS velocity error to 0.05 m/s. In addition, the NSV approach and the VHV approach also performed better than the TC approach. The reason for this result is the assumption that the AUV does not have velocity in the sway direction, and in the case of a straight trajectory this assumption is reasonable. The PLCF obtained the worst performance but still improved the standalone INS solution by 131%.

Figure 10 shows the RMS errors of the orientation misalignment. The standalone INS performance is 1.4° after 250 s. All proposed approaches improved the standalone INS; in particular the NSV and VB approaches achieved 33% improvement. The PLCF attitude error seemed to converge in time, due to the observation of the roll angle (in this method we use surge velocity aiding).

4.2.2. Trajectory #2 Simulation Results

The velocity RMS errors for all the methods (including TC) are presented in Figure 11. Not presented in the figure is the standalone INS performance, which is 31.8 m/s of error after 250 s. All the proposed approaches improved the standalone INS; for example, the PLCF approach reduced the RMS velocity error to 1.2 m/s. In addition, the NSV approach also managed to perform better than the TC approach, but only with minor differences. Notice that all approaches, including TC, achieved better performance than in trajectory #1 due to the maneuvering of the AUV.

Figure 12 shows the RMS error of the attitude errors. The standalone INS performance is 1.26° after 250 s. As in the velocity RMS errors, all proposed approaches improved the standalone INS, in particularly the PLCF, which achieved a 26% improvement and the same performance as TC.

4.2.3. Trajectory #3 Simulation Results

Figure 13 presents the velocity RMS errors of all ELC approaches and the TC approach. The standalone INS velocity error in this trajectory is 36.5 m/s after 250 s. As seen in Figure 10, all the proposed approaches improved the standalone INS, and the VB approach reduced the RMS velocity error to 0.5 m/s. The NSV approach managed to obtain similar performance to that of the TC, although its velocity error oscillates due to the change in the maneuvering direction. In the first 30 s, the velocity errors of methods VB and PLCF drift rapidly in time. The reason for this behavior is that the vehicle performs diving in that time period. After the TAUV completes the diving and reaches its desired depth, the velocity error of the VB approach does not drift in time.

Figure 14 shows the RMS errors of the attitude. The standalone INS performance is 1.29° after 250 s. All proposed approaches improved the standalone INS; in particular, the NSV achieved 38% improvement and even achieved an improvement of 12% relative to the TC approach.

5. Conclusions

In the presented study the ELC approach for INS/DVL fusion with partial DVL measurements was investigated, applying four different methods. The TAUV simulation was developed in order to evaluate the ELC methods.

Results show that the presented approaches significantly improved the navigation performance of the TAUV in cases where only partial measurements from the DVL are available. In order to implement the suggested approach in the TAUV, only software modifications are needed.

Acknowledgments

We would like to acknowledge the Technion Autonomous System and Robotics Program for funding the TAUV project. We would like to thank Gil Iosilevskii for providing us with the TAUV hydrodynamic model. We also gratefully acknowledge the financial support of Rafael Advanced Defense Systems Ltd.

Author Contributions

Asaf Tal, Itzik Klein and Reuven Katz conceived and designed the ELC approach. Asaf Tal, Itzik Klein and Reuven Katz contributed to the writing of the paper. The 6DOF simulation was designed by Asaf Tal.

Conflicts of Interest

The authors declare no conflict of interest.

References

Antonelli, G. Underwater Robots, Motion and Force Control of Vehicle Manipulator Systems, 2nd ed.; Springer: Cassino, Italy, 2005. [Google Scholar]
Panish, R.; Taylor, M. Achieving high navigation accuracy using inertial navigation systems in autonomous underwater vehicles. In Proceedings of the 2011 IEEE OCEANS, Santander, Spain, 6–9 June 2011.
Jalving, B.; Gade, K.; Bovio, E. Integrated inertial navigation systems for auvs for rea applications. Norwegian Defence Research Establishment, NATO Undersea Research Centre. Available online: http://www.navlab.net/Publications/Integrated_Inertial_Navigation_Systems_for_AUVs_for_REA_Applications.pdf (accessed on 16 February 2017).
Kongsberg, Autonomous Underwater Vehicle—AUV the Hugin Family. Available online: https://www.km.kongsberg.com/ks/web/nokbg0397.nsf/AllWeb/A6A2CC361D3B9653C1256D71003E97D5/$file/HUGIN_Family_brochure_r2_lr.pdf (accessed on 5 October 2016).
Alahyari, A.; Rozbahani, S.G.; Habibzadeh, A.; Alahyari, R.; Dousti, M. INS/DVL positioning system using kalman filter. Aust. J. Basic Appl. Sci. 2011, 5, 1123–1129. [Google Scholar]
Lee, C.-M.; Lee, P.-M.; Hong, S.-W.; Kim, S.-M. Underwater navigation system based on inertial sensor and doppler velocity log using indirect feedback kalman filter. Offshore Polar Eng. 2005, 15, 88–95. [Google Scholar]
Kondo, H.; Maki, T.; Ura, T.; Sakamaki, T. AUV navigation based on multi-sensor fusion for breakwater observation. In Proceedings of the ISARC2006 International Symposium on Automation and Robotics in Construction, Tokyo, Japan, 3–5 October 2006.
He, B.; Zhang, H.; Li, C.; Zhang, S.; Liang, Y.; Yan, T. Autonomous navigation for autonomous underwater vehicles based on information filters and active sensing. Sensors 2011, 11, 10958–10980. [Google Scholar] [CrossRef] [PubMed]
Leonard, J.J.; Bennett, A.A.; Smith, C.M.; Feder, H.J.S. Autonomous Underwater Vehicle Navigation; Departure of Ocean Engineering, Massachuesttes Institute of technology: Boston, MA, USA, 1998. [Google Scholar]
Sturrers, L.; Liu, H.; Titman, C.; Brown, D.J. Navigation technologies for autonomous underwater vehicles. IEEE Trans. Syst. 2008, 38, 581–589. [Google Scholar]
Teixeira, F.C.; Pascoal, A.O.M. Geophysical navigation of autonomous underwater vehicles using geomagnetic information. IFAC Proc. Vol. 2016, 41, 178–183. [Google Scholar] [CrossRef]
Livshitz, I. Development of a Navigation System for a Small Autonomous Underwater Vehicle; Technion—Israel Institute of Technology: Haifa, Israel, 2016. [Google Scholar]
Jalving, B.; Gade, K.; Svartveit, K.; Willumsen, A.; Sorhagen, R. Dvl velocity aiding in the hugin 1000 integrated inertial navigation system. Model. Identif. Control 2004, 25, 223–235. [Google Scholar] [CrossRef]
Klein, I.; Diamant, R. Observability analysis of DVL/PS aided ins for a maneuvering auv. Sensors 2015, 15, 26818–26837. [Google Scholar] [CrossRef] [PubMed]
Groves, P.D. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems; Artech House: Boston, MA, USA, 2008. [Google Scholar]
Sbg Systems. Available online: https://www.sbg-systems.com/products/ekinox-ins-high-performance-mems-inertial-system (accessed on 1 October 2016).
Miller, P.A.; Farrell, J.A.; Zhao, Y.; Djapic, V. Autonomous underwater vehicle navigation. IEEE J. Oceanic Eng. 2010, 35, 663–678. [Google Scholar] [CrossRef]
LaPointe, C.E.G. Virtual Long Baseline (Vlbl) Autonomous Underwater Vehicle Navigation Using a Single Transponder; Massachuesttes Institute of Technology: Boston, MA, USA, 2006. [Google Scholar]
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]
Klein, I.; Filin, S.; Toledo, T. A modified loosely-coupled approach. J. Appl. Geodesy 2011, 5, 87–98. [Google Scholar] [CrossRef]
Breivik, M.; Fossen, T.I. Guidance Laws for Autonomous Underwater Vehicles. Available online: http://www.diva-portal.org/smash/record.jsf?pid=diva2%3A346219&dswid=-1553 (accessed on 16 February 2017).
Farrell, J.A. Aided Navigation GPS with High Rate Sensors; McGraw-Hill: New York, NY, USA, 2008. [Google Scholar]
Ribeiro, M.I. Kalman and Extended Kalman Filters: Concept, Derivation and Properties. Institute for Systems and Robotics Instituto Superior Tecnico, 2004. Available online: http://users.isr.ist.utl.pt/~mir/pub/kalman.pdf (accessed on 16 February 2017).
Terejanu, G.A. Extended Kalman Filter Tutorial. Department of Computer Science and Engineering, University at Buffalo. Available online: https://homes.cs.washington.edu/~todorov/courses/cseP590/readings/tutorialEKF.pdf (accessed on 16 February 2017).
Teledyne RD Instruments, Adcp Coordinate Transformation Formulas and Calculations. 2008. Available online: https://wiki.oceannetworks.ca/download/attachments/44236966/adcp+coordinate+transformation_Jan08.pdf?version=1&modificationDate=1423090844000 (accessed on 16 February 2017).
Seapilot Doppler Velocity log 300 khz/600 khz/1200 khz. Available online: http://www.m-b-t.com/fileadmin/redakteur/Ozeanographie/ROWE/LO-RES-SeaPILOT-300-600-1200-DVL-Brochure-27-Dec-2013.pdf (accessed on 10 September 2016).
Kinsey, J.C.; Eustice, R.M.; Whitcomb, L.L. A Survey of Underwater Vehicle Navigation: Recent Advances and New Challenges. IFAC Conference of Manoeuvering and Control of Marine Craft. 2006. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.134.5601&rep=rep1&type=pdf (accessed on 16 February 2017).

Figure 1. Scheme of loosely coupled and tightly coupled approaches for INS/DVL fusion. AUV: Autonomous underwater vehicle; DVL: Doppler velocity log; INS: inertial navigation system.

Figure 2. Technion autonomous underwater vehicle six degrees of freedom (TAUV 6 DOF) simulation layout. IMU: inertial measurement unit.

Figure 3. Top level diagram of the TAUV navigation system. EKF: extended Kalman filter. GPS: global positioning system.

Figure 4. Platform frame and tangent frame in respect to the TAUV.

Figure 5. DVL orientation setup in the TAUV platform frame.

Figure 6. Extended loosely coupled (ELC) shame for utilizing partial measurements from the DVL.

Figure 7. Block diagram for implementation of ELC approach in the TAUV navigation system.

Figure 8. Simulation trajectories used to examine the ELC approaches.

Figure 9. Root mean square (RMS) errors of velocity for all approaches obtained from trajectory #1. TC: tightly coupled.

Figure 10. RMS errors of attitude for all approaches obtained from trajectory #1.

Figure 11. RMS errors of velocity for all approaches obtained from trajectory #2.

Figure 12. RMS errors of attitude for all approaches obtained from trajectory #2.

Figure 13. RMS errors of velocity for all approaches obtained from trajectory #3.

Figure 14. RMS errors of attitude for all approaches obtained from trajectory #3.

Table 1. ELC Methods.

**Table 1.** ELC Methods.
Method	Assumtion	External Information	Calculations	Setting
VB (Section 3.1)	None	Last estimated velocity vector	Velocity vector Equation (24)	None
NSV (Section 3.2)	Zero sway velocity	None	Two velocity components Equation (33)	None
PLCF (Section 3.3)	None	None	One velocity component Equation (40)	Two velocity components with infinity variance
VHV (Section 3.4)	None	Last estimated heave velocity	Two velocity components Equations (40) and (45)	One velocity component with infinity variance

Table 2. TAUV 6DOF simulation and navigation parameters.

**Table 2.** TAUV 6DOF simulation and navigation parameters.
Parameter	Value
Time duration	250 (s)
AUV velocity	2 (m/s)
Accelerometer bias	0.5 (mg/h)
Gyro bias	3 (°/h)
Accelerometers noise	0.072 (m/s/√h)
Gyro noise	0.34 (°/√h)
Accelerometer bias random walk	$1 \times 10^{- 5} (m / s^{2} / \sqrt{s})$
Gyro bias random walk	$2.8 \times 10^{- 5}$ (°/s $/ \sqrt{s}$ )
Position initial error	north: 2 (m), east: 2 (m), height: 2 (m)
Velocity initial error	u: 0.05 (m/s) × v: 0.05 (m/s) × w: 0.05 (m/s)
Attitude initial error	Yaw: 1.14 (°) roll/pitch: 0.57 (°)
IMU rate	150 (Hz)
DVL rate	1 (Hz)
DVL noise	0.042 (m/s)
DVL bias	0.005 (m/s)
DVL bias random walk	$5 \times 10^{- 5}$ (m/s/√s)
DVL scale factor	0.7 (%)
DVL scale factor random walk	$5 \times 10^{- 3}$ (%/√s)
Magnometer noise	Yaw: 5.72 (°) × roll/pitch: 1.15 (°)
Pressure sensor noise	194 (Mpa)
Magnometer rate	0.5 (Hz)
Pressure sensor rate	0.25 (Hz)

© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Tal, A.; Klein, I.; Katz, R. Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements. Sensors 2017, 17, 415. https://doi.org/10.3390/s17020415

AMA Style

Tal A, Klein I, Katz R. Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements. Sensors. 2017; 17(2):415. https://doi.org/10.3390/s17020415

Chicago/Turabian Style

Tal, Asaf, Itzik Klein, and Reuven Katz. 2017. "Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements" Sensors 17, no. 2: 415. https://doi.org/10.3390/s17020415

APA Style

Tal, A., Klein, I., & Katz, R. (2017). Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements. Sensors, 17(2), 415. https://doi.org/10.3390/s17020415

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

Article Menu

Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements

Abstract

1. Introduction

2. TAUV 6DOF Simulation

2.1. TAUV Navigation System

2.1.1. Navigation Filter

2.1.2. Doppler Solution

2.1.3. DVL Fusion

3. INS/DVL Fusion with the ELC Approach

3.1. Virtual Beam

3.2. Nullfying Sway Velocity

3.3. Partial LC Fusing

3.4. Virtual Heave Velocity

3.5. Summary of ELC Methods

3.6. ELC Implemenation

4. Analysis and Results

4.1. Simulation Scenarios and Parameters

4.2. Simulation Results

4.2.1. Trajectory #1 Simulation Results

4.2.2. Trajectory #2 Simulation Results

4.2.3. Trajectory #3 Simulation Results

5. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI