
divergence if the initial state estimate is not close to the true
initial state.
The most common way this problem arises is through
initialization of the heading angle. In order to provide
an accurate initial heading angle to the EKF, a two-step
approach consisting of coarse alignment followed by fine
alignment is commonly used [4]. Other approaches to the
handling of large heading error with an EKF include [6],
where an EKF model that can handle large initial alignment
errors is developed. This approach still uses an EKF, which
requires linearization, and must also switch to a different
inertial navigation model when alignment error is reduced
below a certain level, much like the coarse and fine alignment
approach. In this work, use of the UKF-M allows for filter
convergence in the presence of large initial heading error
without any special modifications to the filter or system
model.
Other variations of the Kalman filter have been proposed
in the literature for nonlinear filtering, such as the unscented
Kalman filter (UKF) [7]. Instead of linearizing the state
propagation and measurement equations, the UKF uses a set
of points, named sigma points, distributed such that their
empirical mean and covariance estimate the state mean and
covariance. These sigma points are propagated through the
nonlinear equations in the prediction and update steps. In
this way, the UKF can be accurate to a higher order than
the EKF and does not require linearization of the inertial
navigation equations. Approaches to using the UKF for aided
inertial navigation system have been explored in works such
as [8], [9], and [10]. However, implementation challenges
exist when using the UKF for inertial navigation. If vehicle
roll, pitch, and heading are represented in the state vector as
Euler angles, the algorithm must compensate for the fact that
0 degrees and 360 degrees are equivalent. When propagating
sigma points through the inertial navigation equations, one or
more of the sigma points may wrap, leading to an inaccurate
representation of the state mean and covariance.
The unscented Kalman filter on manifolds (UKF-M),
proposed in [11], is a modification of the UKF that is able
to directly handle states that evolve on the manifold SO(3),
such as vehicle attitude. With vehicle attitude represented by
an rotation matrix in SO(3), the effect of wrapping of Euler
angles on sigma points can be avoided. In [2], the authors
present an inertial navigation algorithm for AUVs based on
the UKF-M. The authors demonstrate the ability of the filter
to accurately estimate the state covariance through Monte-
Carlo simulations and show that the filter can run in real
time on an AUV. In this work, we build upon the results in
[2] by using the UKF-M for inertial navigation using sensor
models with lever arm compensation and validating a the
filter in real time with data from an AUV survey mission.
III. INERTIAL NAVIGATION EQUATIONS
We adopt the notation and navigation frame inertial navi-
gation formulations presented in [4]. The inertial navigation
equations involve three reference frames: the Earth-centered
Earth-fixed (ECEF) frame, the body frame, and the naviga-
tion frame.
The origin of the Earth-centered Earth-fixed (ECEF) ref-
erence frame is attached to the geometric center of the Earth.
The positive z-axis of the ECEF reference frame is aligned
with the Earth’s axis of rotation and passes through the north
pole. The positive x-axis passes through the intersection of
the equator and the prime-meridian. The y-axis completes
a right-handed coordinate system [1]. Note that, since the
ECEF frame is fixed to the earth, it is rotating with respect
to an inertial frame [3]. Values expressed in this frame are
indicated with an e.
The origin of the body frame is fixed to a point on the
AUV with the positive x-axis passing through the nose of
the AUV, the positive z-axis pointing down, and the y-axis
completing a right-handed coordinate system. In this work,
the body frame is taken to be the IMU measurement frame.
Values expressed in the body frame are denoted with b.
The navigation frame is a coordinate frame that is tangent
to the earth’s surface with an origin co-located with the origin
of the body frame. The navigation frame’s north direction
points toward true north, the down direction points toward
the earth’s center perpendicular to the tangent plane, and the
east direction points east orthogonal to the north and down
directions. Values in the navigation frame are notated with
an n.
Fig. 2. Illustration of the inertial navigation reference frames at a moment
in time. The navigation frame origin is fixed to the body frame origin with
position described by its geodetic latitude (Lb) and longitude (λb).
The nonlinear kinematic equations that relate attitude,
velocity, and position of the body frame to angular velocity
and linear acceleration of the body frame are
˙
Cn
b=Cn
bΩb
ib −(Ωn
ie +Ωn
en)Cn
b
˙
vn
eb =fn
ib +gn
b−(Ωn
en + 2Ωn
ie)vn
eb
˙
pb=Tp
r(n)vn
eb
(1)
where Cn
bis the rotation matrix from the body frame to
the navigation frame, vn
eb is the velocity of the body frame
relative to the ECEF frame, expressed in the navigation