Next Article in Journal
Distance Calibration between Reference Plane and Screen in Direct Phase Measuring Deflectometry
Previous Article in Journal
A Weld Position Recognition Method Based on Directional and Structured Light Information Fusion in Multi-Layer/Multi-Pass Welding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fast Alignment of SINS for Marching Vehicles Based on Multi-Vectors of Velocity Aided by GPS and Odometer

School of Instrumentation Science and Opto-electronics Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(1), 137; https://doi.org/10.3390/s18010137
Submission received: 8 November 2017 / Revised: 26 December 2017 / Accepted: 2 January 2018 / Published: 5 January 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In the strap-down inertial navigation system (SINS), the initial attitude matrix is acquired through alignment. Though there were multiple valid methods, alignment time and accuracy are still core issues, especially regarding the condition of the motion carrier. Inspired by the idea of constructing nonlinear vectors by velocity in a different coordinate frame, this paper proposes an innovative alignment method for a vehicle-mounted SINS in motion. In this method, the core issue of acquiring the attitude matrix is to calculate the matrix between the inertial frame and the initial body frame, which can be constructed through the nonlinear velocity vectors’ information from the GPS and the odometer at different moments, which denominate the multi-vector attitude determination. The possibility of collinearity can easily be avoided by a turning movement. The characteristic of propagation of error is analyzed in detail, based on which an improved method is put forward to depress the effect of random noise. Compared with the existing alignment methods, this method does not use the measurement information of accelerometers. In order to demonstrate its performance, the method is compared with the two-position alignment method and the traditional two-stage alignment method. Simulation and vehicle-based experiment results show that the proposed alignment method can establish an attitude reference in 100 s with an azimuth error of less than 0.06°, and that the accuracy does not have a strong correlation with the accelerometer.

1. Introduction

Alongside the significant advances in the inertial sensors and navigation technologies, strap-down inertial navigation systems (SINS) have reduced in structure complexity and become a key member of the navigation system [1]. The strap-down inertial navigation system, simplified as SINS, is fixedly installed in the carrier, and the gyroscopes and accelerometers immediately sense and detect the dynamic characteristics of the carrier. The navigation solution is achieved through a mathematical platform, which is different from the gimbaled inertial navigation system. The SINS is an independent navigation system that has been used in a wide range of applications, covering the navigation of marines, vehicles, missiles, and aircraft [2].
It is widely known that the state update of the SINS is achieved based on numerical integration [3]; therefore, it is necessary to know the initial navigation parameters, including velocity, position, and attitude for navigation calculation. The initial alignment of the SINS is mainly about the determination of the initial attitude matrix between the body frame and the navigation frame, since the velocity and position can be easily obtained by external reference [4]. Initial alignment is one of the core technologies for the SINS; the requirements of initial alignment are high accuracy and short time [1]. It is urgently expected to achieve an accurate alignment within a short time; however, it is difficult and contradictory to achieve such an endpoint, especially for the SINS on a moving base [5,6,7].
The existing alignment methods typically comprise two steps: coarse alignment and fine alignment. The goal of a coarse alignment is to provide a roughly known initial attitude within a short time, while the fine alignment mainly drives the misalignment angles to zero [8]. An analytic method using the known gravity and Earth rate signals in the local level frame and the measurements obtained by using accelerometers and gyros to obtain the system attitude directly is often used for coarse alignment when the SINS is in stationary conditions [9]. Meanwhile, optimal estimation methods based on filter are used to realize the fine alignment [10].
To improve the self-alignment performance of the SINS in stationary bases, the multi-position alignment method was put forward, as it can improve the observability of the SINS by changing the position or equivalently by rotating vehicle body [2,11]. In contrast, the two-position method is mostly defined in the static base. It can evaluate some bias and factor errors of the inertial sensors, and achieve higher alignment accuracy. However, two parameters must be considered in order to employ the method effectively: the minimum number of positions required to obtain the best alignment, and the optimum configuration for a given multi-position alignment [11].
Besides alignment in stationary conditions, a substantial amount of research work about initial alignment on a moving base has been launched. Transfer alignment is one of the most important methods to align inertial navigation systems on moving bases [12,13,14,15]. The alignment of a slave system with respect to a master reference is called transfer alignment, which is widely used for the SINS in ships and planes for its speed and high accuracy compared with other alignment methods [16]. Choosing a reasonable matching method is one of the key points of transfer alignment. The transfer alignment method has entered a rapid development and widely used period after the “velocity plus attitude” matching method, which was first put forward by Kain and Clouiter in 1989 [16]. Transfer alignment models include linear and nonlinear models, which respectively correspond with the small and large misalignment angles [17].
In addition to transfer alignment, some other research studies about alignment methods on an in-motion base are also done. Gravity and Earth rational rate cannot be accurately measured due to the high frequency noise in the shacking base; as a result, some methods to eliminate the noise from the outputs of inertial sensors are proposed to achieve higher accuracy [18,19,20,21]. Since the outputs of inertial sensors installed on the board directly are vulnerable to external interference, many researchers devoted their energy to improving the performance of the initial alignment with external disturbance.
Since 2000, the articles from IXBlue Company claimed that their newest product—Octans—can fulfill alignment within five minutes under any swinging condition by observing the gravitational apparent motion. As a result, various alignment methods have been proposed based on the idea of tracking gravitational apparent motion [17,21,22,23,24,25,26]. Similar to the analytic method, the basic idea of both methods for obtaining initial attitude is to construct noncollinear vectors by tracking gravitational apparent motion in an inertial frame [26,27,28]. There are two problems caused by the alignment mechanism: how to avoid collinear possibility between two vectors, and how to recognize the gravitational apparent motion.
Based on the idea of constructing noncollinear vectors, many research studies have attempted to achieve the alignment in marching bases with much more simple methods. With the aid of velocity from the GPS, Dr. Peter M. G. Silson and Dr. Simon Jordan introduced a method to accomplish the alignment with the velocity vectors obtained from the GPS and the integrated accelerations [27]. For the low-cost internal navigation system (INS), Han, Songlai introduced an alignment method with two-stage Kalman filtering; with the coarse alignment and the fine alignment, the proposed alignment approach can complete the initial alignment more quickly and more accurately [29]. For the application in landed carriers, Chang, Lubin introduced an attitude estimation-based alignment method, which could attenuate the disturbance in the odometer to a certain extent, and with a low-pass finite impulse response digital filter, the novel method does have a better performance compared with the attitude determination-based method [30]. However, there is still much space to shorten the alignment time and improve the accuracy. The information from the GPS and the odometer could be combined organically to achieve much better performance in the initial alignment of the SINS.
In this paper, a novel alignment method for the marching vehicle is proposed based on the idea of matrix decomposition. In this method, the constant matrix from the initial body frame to the inertial frame is calculated through vectors constructed by the velocity from the GPS in the navigation frame, and the velocity from the odometer in the body frame. Using the outputs of gyros tracks the changes of body frame, and using the position information (latitude and longitude) from GPS and the alignment time calculates the matrix from the inertial frame to the navigation frame. Compared with the alignment method based on the tracking gravitational apparent motion, this method can easily avoid the collinear possibility by a turning movement. Compared with the current alignment for a marching vehicle [31], this method can fulfill alignment in a short time, and does not need the coarse alignment stage before marching. Compared with the existing alignment methods, this method does not use the measurement information of accelerometers. The outstanding characteristic of this method is that the accuracy of the azimuth is superior, especially with low-precision inertial sensors, and the initial alignment is completed totally in motion. Through the simulation, vehicle experiments, and the comparison with the two-position alignment method, the new alignment mechanism can achieve higher azimuth accuracy in much shorter time.

2. New Alignment Mechanism Based on the GPS and the Odometer

2.1. Coordinate Frame Definitions

The coordinate frames used in this paper are defined as follows:
  • The e frame is the Earth-fixed coordinate frame, whose origin is the center of the Earth. The x e axis is in the equatorial plane and points to the prime meridian, and the z e axis is parallel to the rotation axis of the Earth. The y e axis completes the right-handed coordinate system.
  • The e 0 frame is also the Earth-fixed coordinate frame, whose origin is the center of the Earth. The x e 0 axis is in the equatorial plane, and points to the meridian of the vehicle’s position at the beginning of initial alignment, and the z e 0 axis is parallel to the rotation axis of the Earth. The y e 0 axis completes the right-handed coordinate system.
  • The i 0 frame is formed by fixing the e 0 frame at the beginning of alignment in inertial space.
  • The n frame is the instantaneous navigation coordinate frame, which is the local level coordinate frame. The x n axis points to the east, the y n axis points to the north, and the z n axis points upwards.
  • The b frame is the instantaneous body coordinate frame, which is defined as the common “Right–Forward–Up” frame.
  • The i b 0 frame is formed by fixing the b frame at the beginning of alignment in the inertial space.

2.2. Constructing Nonlinear Vectors by Turning Movement

The alignment method of the SINS in this paper, which is aided by the GPS and the odometer, is designed for a marching vehicle on land. The crucial point, whose first step is to construct nonlinear vectors, is how to figure out the matrix between the i 0 frame and the i b 0 frame.
The i 0 frame and the n frame on the Earth are shown as Figure 1. When the vehicle moves along a straight line for a short period of time, the n frame can be recognized as unchangeable but for its origin, thus the velocity of vehicle in the n frame is parallel in this period of time, whose projections in the i 0 frame are collinear consequently. Then, how can the nonlinear vectors be obtained in a short period of time?
Obviously, if the vehicle moves with a turning at a certain moment, the vehicle’s velocity in the n frame before turning is not parallel to that after turning. It is assumed that: the vehicle at position A moves with a ground speed of V n ( t 1 ) at t 1 moment, and the direction of the velocity is along the y n axis, as shown in Figure 2; then, the vehicle moves along the smooth curve from position A to position B at t 2 moment, where its ground speed is recorded as V n ( t 2 ) , and the motion directs to the reverse of the axis x n . It is obvious that V n ( t 1 ) and V n ( t 2 ) are not parallel, and the projections of them in the i 0 frame are nonlinear vectors as a consequence.
It is worth noting that the nonlinear vectors are obtained by a turning movement with the angle of 90°, as shown in Figure 2. It is quite understandable that V n ( t 1 ) and V n ( t 2 ) will be parallel again if the turning angle becomes 180°, which should be avoided. Theoretically, if only the turning angle is larger than 0° and smaller than 180°, the nonlinear vectors can be obtained. In fact, in order to cut down the calculation error and depress the interference from measurement errors of velocity, the turning angle should have a large discrepancy relative to 0° and 180°.

2.3. Alignment Algorithm

The attitude matrix C b n ( t ) changes with time when the vehicle is in motion. In order to obtain it, the matrix C b n ( t ) can be decomposed into two parts, as follows [32]:
C b n ( t ) = C i 0 n ( t ) C b i 0 ( t )
where C N M is the attitude matrix between N and M .
The matrix C i 0 n ( t ) in Equation (1) can be calculated as follows [33]:
C i 0 n ( t ) = C e n ( t ) ( C e n 0 ) T C e 0 n 0 C i 0 e 0 ( t )
The matrix C e n ( t ) C e n 0 C e 0 n 0 can be acquired through the latitude and longitude of the vehicle, while the matrix C i 0 e 0 ( t ) can be determined constantly through the alignment time. Therefore, the specific expression of matrix C i 0 n is as follows:
C i 0 n ( t ) = [ sin [ Δ λ + ω i e t ] cos [ Δ λ + ω i e t ] 0 sin L t cos [ Δ λ + ω i e t ] sin L t sin [ Δ λ + ω i e t ] cos L t cos L t cos [ Δ λ + ω i e t ] cos L t sin [ Δ λ + ω i e t ] sin L t ]
where Δ λ can be determined by Δ λ = λ t λ 0 , and λ 0 is the meridian of the vehicle’s position at the beginning of initial alignment; L t and λ t are the latitude and longitude at the current time t , respectively; λ 0 , L t and λ t can be obtained from the GPS, or the L t and λ t can be computed through λ 0 and the velocity from the GPS; and ω i e refers to the Earth rotation angular velocity.
The matrix C b i 0 ( t ) indicates the transformation relation between the body frame and the initial inertial coordinate frame. Furthermore, C b i 0 ( t ) can be decomposed as follows [32]:
C b i 0 ( t ) = C i b 0 i 0 ( t ) C b i b 0 ( t )
With the information from the gyroscope, the matrix C b i b 0 ( t ) in Equation (4) can be determined constantly through the attitude-updating algorithm [24]:
C ˙ b i b 0 ( t ) = C b i b 0 ( t ) [ ω i b 0 b b ( t ) × ]
where the ω i b 0 b b is the gyro measurement value, and [•×] represents the skew-symmetric matrix of vector •.
According to Equations (1) and (5), the problem of solving C b n ( t ) in the SINS alignment can be attributed to the computing of C i b 0 i 0 . Inspired by the idea of tracking gravitational apparent motion to construct noncollinear vectors [22,25], C i b 0 i 0 can be obtained if the projection of the identical physical vector is known in the i 0 frame and the i b 0 frame synchronously. Then, the specific expression of C i b 0 i 0 will be deduced in detail.
In the SINS aided by the GPS and the odometer, the velocity of a vehicle in the n frame can be obtained from the GPS, while the odometer can provide the velocity of a vehicle in the b frame through misalignment compensation. The misalignment angles between the odometer and the b-frame, computed as α = [ α θ α γ α ψ ] T , are calibrated in advance. The velocity of a vehicle in the n frame and the b frame are defined respectively as follows:
{ V n ( t ) = [ V E n V N n V U n ] T V b ( t ) = [ V x b V y b V z b ] T = C b m [ 0 V m 0 ] T
where V E n , V N n , and V U n denote the vehicle’s velocity in the n frame pointing east, north, and upwards, respectively; V x b , V y b , and V z b denote the vehicle’s velocity in the b frame, and V m is the velocity from the odometer pointing forwards; and C b m denotes the transform matrix between the odometer and the b-frame. In most cases, the C b m could be simplified as follows, considering that α θ , α γ and α ψ are always small quantities.
C b m = [ 1 α ψ α γ α ψ 1 α θ α γ α θ 1 ]
It is obvious that the velocity’s projection of a vehicle in the i 0 frame ( V i 0 ( t ) ) can be obtained by V n ( t ) , while that in the i b 0 frame, ( V i b 0 ( t ) ) can be calculated through V b ( t ) as follows:
{ V i 0 ( t ) = C n i 0 ( t ) V n ( t ) V i b 0 ( t ) = C b i b 0 ( t ) V b ( t )
where C n i 0 ( t ) = ( C i 0 n ( t ) ) T can be obtained by Equation (3); and C b i b 0 ( t ) can be obtained by Equation (5).
The relationship between V i 0 ( t ) and V i b 0 ( t ) is explicit, as follows:
V i b 0 ( t ) = C i 0 i b 0 V i 0 ( t )
where C i 0 i b 0 = ( C i b 0 i 0 ) T . The matrix can be computed through a series of noncollinear vectors, which can be acquired by choosing a different time t j in the alignment stage. Therefore, the formula is as follows:
V j i b 0 ( t ) = C i 0 i b 0 V j i 0 ( t )    j = 1 , 2 , m
In order to compute the optimal solution of C i 0 i b 0 , construct the indicator function as follows [24]:
J * ( C i 0 i b 0 ) = 1 2 j = 1 m w j | V j i b 0 C i 0 i b 0 V j i 0 | 2 = min
where w j is the weighting coefficient, and the sum of squares | V j i b 0 C i 0 i b 0 V j i 0 | can be calculated as follows:
| V j i b 0 C i 0 i b 0 V j i 0 | 2 = ( V j i b 0 C i 0 i b 0 V j i 0 ) T ( V j i b 0 C i 0 i b 0 V j i 0 ) = [ ( V j i b 0 ) T ( V j i 0 ) T ( C i 0 i b 0 ) T ] ( V j i b 0 C i 0 i b 0 V j i 0 ) = | V j i b 0 | 2 ( V j i b 0 ) T C i 0 i b 0 V j i 0 ( V j i 0 ) T ( C i 0 i b 0 ) T V ˜ i r + ( V j i 0 ) T ( C i 0 i b 0 ) T C i 0 i b 0 V j i 0 = | V j i b 0 | 2 + | V j i 0 | 2 2 ( V j i b 0 ) T C i 0 i b 0 V j i 0
Therefore, Equation (10) is equivalent to the formula below:
J * ( C i 0 i b 0 ) = 1 2 j = 1 m w j | V j i b 0 C i 0 i b 0 V j i 0 | 2 = 1 2 j = 1 m w j ( | V j i b 0 | 2 + | V j i 0 | 2 ) j = 1 m w j ( V j i b 0 ) T C i 0 i b 0 V j i 0 = min
The indicator function can be reconstructed as follows:
J ( C i 0 i b 0 ) = j = 1 m w j ( V j i b 0 ) T C i 0 i b 0 V j i 0 = max
The simplified computing of the indicator function is stated as Equation (14):
J ( C i 0 i b 0 ) = tr ( [ w 1 ( V 1 i b 0 ) T w 2 ( V 2 i b 0 ) T w m ( V m i b 0 ) T ] C i 0 i b 0 [ V 1 i 0 V 2 i 0 V m i 0 ] ) = tr ( C i 0 i b 0 [ V 1 i 0 V 2 i 0 V m i 0 ] [ w 1 ( V 1 i b 0 ) T w 2 ( V 2 i b 0 ) T w m ( V m i b 0 ) T ] ) J ( C i 0 i b 0 ) = tr ( C i 0 i b 0 j = 1 m w j V j i 0 ( V j i b 0 ) T ) = max
With the reconstructed indicator function as Equation (14), the matrix [ j = 1 m w j V j i 0 ( V j i b 0 ) T ] can be decomposed to U D V T , which is known as singular value decomposition, supposing det ( j = 1 m w j V j i 0 ( V j i b 0 ) T ) > 0 . The optimal attitude matrix can be computed as follows:
C i 0 i b 0 = U V T
Generally speaking, C i b 0 i 0 in Equation (15) cannot be orthogonal for the calculated errors caused by measurement errors of the GPS and the odometer, and this C i b 0 i 0 should be further orthogonalized as follows [22]:
( C i b 0 i 0 ) o = C i b 0 i 0 [ ( C i b 0 i 0 ) T C i b 0 i 0 ] 1 2
where (•)o represents the result of matrix • through orthogonalization.
Ultimately, substituting Equations (3), (4), and (16) into Equation (1), the initial alignment can be fulfilled.
It deserves much attention that V i 0 ( t j ) and V i 0 ( t j ) should be noncollinear, which can be guaranteed by the turning movement of the vehicle, as analyzed in Section 3.1. It is obvious that accelerometer measurements are not used in this alignment method, which is the different from the existing alignment method.

3. Error Analysis and the Improved Alignment Method

3.1. Error Analysis

The novel alignment method is detailed as presented above, and it is essential to figure out the influence factor of the initial errors. The real-time propagation error of the alignment procedure based on a multi-vector may be somewhat complicated, and it may be more important to get the elements relating to the misalignment. To obtain the error equations, similar to the coarse inertial alignment method, two different moments are selected arbitrarily to execute the alignment errors.
Without loss of generality, two specified moments, T l and T m , are selected (assuming T l < T m ), and the simplified calculation of C i 0 i b 0 can be determined as follows, consulting the coarse alignment method based on the inertial frame [28].
C i 0 i b 0 = ( C i b 0 i 0 ) T = [ [ V i b 0 ( T l ) ] T [ V i b 0 ( T l ) × V i b 0 ( T m ) ] T [ V i b 0 ( T l ) × V i b 0 ( T m ) × V i b 0 ( T l ) ] T ] 1 [ [ V i 0 ( T l ) ] T [ V i 0 ( T l ) × V i 0 ( T m ) ] T [ V i 0 ( T l ) × V i 0 ( T m ) × V i 0 ( T l ) ] T ]
According to the character of the orthorhombic matrix, the new expression of C i b 0 i 0 is as follows:
C i b 0 i 0 = ( C i 0 i b 0 ) T = C i 0 C i b 0
where the C i 0 and C i b 0 are defined in Equation (19):
{ C i 0 = [ V i 0 ( T l ) , V i 0 ( T l ) × V i 0 ( T m ) , V i 0 ( T l ) × V i 0 ( T m ) × V i 0 ( T l ) ] C i b 0 = [ V i b 0 ( T l ) , V i b 0 ( T l ) × V i b 0 ( T m ) , V i b 0 ( T l ) × V i b 0 ( T m ) × V i b 0 ( T l ) ] 1
Therefore, the attitude matrix at T m can be rewritten as Equation (20):
C b n ( T m ) = C i 0 n ( T m ) C b i 0 ( T m ) = C i 0 n ( T m ) C i b 0 i 0 C b i b 0 ( T m ) = C i 0 n ( T m ) C i 0 C i b 0 C b i b 0 ( T m )
Taking account of the short interval between T l and T m , as the initial alignment is typically fulfilled in minutes or shorter, it is rational to assume that the longitude and latitude remain constant throughout the alignment procedure. Thus, there exists an approximate equality of C i 0 n ( T m ) C i 0 n ( T l ) . Therefore, the attitude matrix at T m is finally determined as Equation (21):
C b n ( T m ) = C i 0 n ( T m ) C i 0 C i b 0 C b i b 0 ( T m ) = C n C i b 0 C b i b 0 ( T m )
where C n = [ V n ( T l ) , V n ( T l ) × V n ( T m ) , V n ( T l ) × V n ( T m ) × V n ( T l ) ] .
Assuming that the variation of the Euler angles of the body frame at T l and T m are Δ θ , Δ γ and Δ ψ , the following equation is acceptable.
Δ C = C b ( T m ) b ( T l ) = C i b 0 b ( T l ) C b i b 0 ( T m ) = [ cos ( Δ γ ) cos ( Δ ψ ) + sin ( Δ γ ) sin ( Δ ψ ) sin ( Δ θ ) cos ( Δ γ ) sin ( Δ ψ ) + sin ( Δ γ ) cos ( Δ ψ ) sin ( Δ θ ) sin ( Δ γ ) cos ( Δ θ ) sin ( Δ ψ ) cos ( Δ θ ) cos ( Δ ψ ) cos ( Δ θ ) sin ( Δ θ ) sin ( Δ γ ) cos ( Δ ψ ) cos ( Δ γ ) sin ( Δ ψ ) sin ( Δ θ ) sin ( Δ γ ) sin ( Δ ψ ) cos ( Δ γ ) cos ( Δ ψ ) sin ( Δ θ ) cos ( Δ γ ) cos ( Δ θ ) ]
In order to simplify the calculation process, the angles Δ θ and Δ γ can be assumed to be zero, taking account of the motion forms in most cases. Thus, the transformation matrix Δ C can be rewritten as follows:
Δ C = [ cos ( Δ ψ ) sin ( Δ ψ ) 0 sin ( Δ ψ ) cos ( Δ ψ ) 0 0 0 1 ]
Based on Equations (21) and (23), it is understandable to yield Equation (24):
C b = C i b 0 C b i b 0 ( T m ) = [ ( Δ C ) T V b ( T l ) , ( ( Δ C ) T V b ( T l ) ) × V b ( T m ) , ( ( Δ C ) T V b ( T l ) ) × V b ( T m ) × ( ( Δ C ) T V b ( T l ) ) ] 1
Therefore, the explicit expression of C b n ( t 2 ) can be calculated through the following equation:
C b n ( T m ) = C i 0 n ( T m ) C i 0 C i b 0 C b i b 0 ( T m ) = C n C b
The attitude angle between the n frame and the b frame is defined as ϕ = [ θ γ φ ] T , where θ , γ , and φ denote pitch, roll, and yaw respectively. It is widely known that ϕ can be expressed as follows:
{ sin θ = C 32 tan γ = C 31 C 33 tan φ = C 12 C 22
where C M N is the element of C b n ( T m ) at row M and column N .
According to Equations (21), (23), (24), and (26), the attitude angle can be computed as follows:
{ sin θ = V U n ( T m ) V y b ( T m ) tan γ = V y b ( T m ) V U n ( T l ) V y b ( T l ) V U n ( T m ) cos ( Δ ψ ) V E n ( T l ) V N n ( T m ) V E n ( T m ) V N n ( T l ) tan φ = V E n ( T m ) V N n ( T m )
Taking the time derivative on both sides of Equation (27), it gives:
{ d θ = V y b ( T l ) d V U n ( T m ) V U n ( T m ) d V y b ( T m ) ( V y b ( T m ) ) 2 cos θ d γ = P d Q + Q d P [ ( V E n ( T m ) V N n ( T l ) V E n ( T l ) V N n ( T m ) ) 2 ] ( sec γ ) 2 d φ = V N n ( T m ) d V E n ( T m ) V E n ( T m ) d V N n ( T m ) ( V N n ( T m ) ) 2 ( sec φ ) 2
where P = ( V E n ( T m ) V N n ( T l ) V E n ( T l ) V N n ( T m ) ) , Q = V U n ( T m ) V y b ( T l ) cos ( Δ ψ ) V U n ( T l ) V y b ( T m ) , and d ( ) denotes differential coefficient of .
The following equations are understandable, considering that the value of pitch and roll are small.
{ sin θ 0 , cos θ 1 sin γ 0 , cos γ 1 V U n ( t ) 0
The relationship among V n ( t ) , V b ( t ) , and ϕ is obviously as follows:
{ V E n ( t ) = V y b ( t ) cos θ sin φ V N n ( t ) = V y b ( t ) cos θ cos φ V U n ( t ) = V y b ( t ) sin θ
Substituting Equations (29) and (30) into Equation (28) yields:
{ d θ = d V U n ( T m ) V y b ( T m ) d γ = d V U n ( T m ) V y b ( T m ) tan ( Δ ψ ) d V U n ( T l ) V y b ( T l ) sin ( Δ ψ ) d φ = cos φ d V E n ( T m ) V y b ( T m ) sin φ d V N n ( T m ) V y b ( T m )
Alignment errors are mainly affected by measurement errors of velocity from the GPS, the turning angle, and the vehicle’s velocity participating in calculation. Alignment errors are also affected by the odometer’s measurement errors, which provide the forward velocity as V y b ( t ) . It is understandable that Δ ψ should be far from 0° and 180°, and high speed leads to small errors.

3.2. Improved Method

According to the analysis in Section 3.1, the main error resource is the measurement errors of speed from the GPS, mostly with the character of stochastic error. It is universally known that integration can depress the effect of random noise; based on this, an improved method to calculate C i b 0 i 0 is put forward.
Taking the time integral on both sides of Equation (7) yields:
{ u = V i 0 ( t ) d t = C n i 0 ( t ) V n ( t ) d t υ = V i b 0 ( t ) d t = C b i b 0 ( t ) V b ( t ) d t
The indicator function of Equation (14) can be rewritten as follows:
J ( C i 0 i b 0 ) = tr ( C i 0 i b 0 j = 1 m w j u j ( υ j ) T ) = max
Through singular value decomposition, the optimal attitude matrix C i 0 i b 0 can be computed in the same way as Equation (15). This C i 0 i b 0 should also be further orthogonalized, as in Equation (16).

4. Simulation and Vehicle Test

4.1. Simulation

4.1.1. Simulation Settings

In order to demonstrate the alignment method proposed in this paper, an experiment is carried out and the motive trajectory could be generated through simulating the actual motion of a vehicle. The vehicle is assumed to have translational movement with a turning at an arbitrary moment among the motion. Firstly, the vehicle is assumed to be heading forward at the speed of 50 m/s at beginning; in the next 25 s, there are kinematics e.g., acceleration and deceleration, pitching and rolling motion (approximately within one degree); then the vehicle starts turning at the speed of 10 m/s; a turning with 50° will be finished in 10 s, i.e., the yaw angle changes from 70° to 120° ( Δ ψ = 50 ° ); the vehicle begins accelerating with the unchangeable acceleration for 10 s; then, with several motions, such as pitching and rolling; all of the kinematics are terminated at the 50th s, when the speed of the vehicle is 60 m/s; then, the vehicle moves along a straight line at this velocity, which lasts for the rest of the time. The whole simulation time can be changed considering the requirements of different methods.
The traditional two-position alignment method can evaluate sensor errors, and in most cases, it can represent a relatively high alignment accuracy of the SINS with specified inertial sensor precision. The process of the two-position method is different from the new method, while it can be the benchmark for the novel alignment method. Therefore, we conducted the simulation as a comparison reference for the new method. The process of two-position alignment is described as: firstly, keep still for 30 s to fulfill the coarse alignment, which can obtain the initial attitude with large errors. Then, turning 90° changes the position only by angular movement, among which a liner Kalman filter is used to drive the misalignment angles to zero.
The initial position of alignment is set to λ 0 = 116.34° and L 0 = 39.98°, and the initial attitude is assumed as ϕ 0 = [ 1 ° 1 ° 70 ° ] T . The sample period of gyro and accelerometer is 0.001 s or 1 millisecond, while the update rate of the velocity from the GPS and the odometer is 10 Hz. The sensor errors are defined in Table 1. The random errors of the GPS in different directions are recognized as irrelevant.
With the above setting, the theoretical outputs of the gyros, the accelerometer, the GPS, and the odometer can be generated by back-stepping the navigation solution of the SINS. When the errors in Table 1 are added to these theoretical outputs, the data sensors output can be produced. In order to evaluate the alignment accuracy, the theoretical values of attitude are used as references. The difference between alignment attitude and the references are regarded as alignment errors.

4.1.2. Simulation Results and Analysis

Error Analysis Validation

The alignment accuracy is mainly related to the measurement errors of the velocity information obtained through the GPS and the odometer, as stated in Section 3.1. To testify to the analysis results, a comparison could be conducted between the “alignment error” and the “calculation error”, which represent the error of the alignment process and the error computed in Equation (31), respectively. The raw alignment method stated in Section 2.3, without the integration of the velocity information, is conducted with the situation in Section 4.1.1.
The simulation results are shown in Figure 3, in which (a)–(c) denote the alignment errors of pitch, roll, and yaw, respectively, and the red squares denote the alignment errors from 75 to 120 s, while the green triangles are the calculation errors in Equation (31). The abscissa represents time, whose unit is second, while the vertical coordinates represents errors with respect to the unit of degree. The statistics of the alignment errors and calculation errors are shown in Table 2. There is little difference between the “alignment error” and the “calculation error”. The simulation results illustrate that the error analysis, deduced in Section 3.1, is reasonable.

Alignment Mechanism Validation

Based on the simulation condition set above, a simulation experiment is carried out. The results of the novel method are shown in Figure 4, in which (a)–(c) denote the alignment errors of pitch, roll, and yaw, respectively, and the accompanying diagram shows the detailed convergence process of the corresponding attitude angle. The abscissa represents time in seconds, while the vertical coordinate represents attitude errors in degrees. The output attitude is not available until the 30th second of the alignment duration. Before the output attitude is available at the 30th s, the pitch, roll, and yaw errors just represent the actual attitude of the carrier. The spikes in (a) and (b) reflect the changing of pitch and roll in the marching of the carrier, as stated in Section 4.1.1. The simulation results illustrate that the novel method can fulfill the initial alignment for motion vehicles in 100 s, with horizontal attitude error about 0.015°, and azimuth error about 0.03°. Comparing with the simulation results of the raw method, as shown in Figure 4, the improved method does have a higher alignment accuracy.
The frequently-used two-position alignment method is carried out to compare with the novel method, and the variation of attitude errors are presented in Figure 5. The convergence time of the horizontal attitude angles is about 100 s, while it requires approximately 250 s for the azimuth to accomplish its constriction. However, the two-position alignment method does have a higher horizontal angle accuracy, and practically the same azimuth precision as the novel method.
Without loss of generality, the simulation of the novel alignment method and the two-position method is promoted 25 times to illustrate the validity and the repeatability of the alignment results. The simulation conditions vary accordingly through changing the maneuver manners of the vehicles, e.g., the pitching or rolling amplitude, the turning angle, the forward speed, etc. The attitude errors are at the 300th s are regarded as alignment errors for the two-position alignment method, while the attitude errors are at the 100th s for the novel method. Therefore, the statistic characteristics of the simulation results are presented in Figure 6 and Figure 7 for the novel method and the two-position alignment method, respectively.
For the two-position alignment method, the attitude errors are definitely controlled within 0.01° for the horizontal angle and 0.4° for the azimuth (peak to peak). Taking account of the inertial sensors, including gyros and accelerometers with biases of 0.1°/h and 500 μg, respectively, the alignment results practically achieve the theoretical precision. Meanwhile, as the novel alignment method does not have strong correlation with the accuracy of gyros and accelerometers, the attitude errors converge to 0.03° for the horizontal angle and 0.04° for the azimuth, and the alignment results are mostly affected by the turning angle, the velocity of the GPS, and the odometer.
Comparing the alignment time and accuracy of these two methods, the novel procedure, as presented at Section 3, can fulfill the entire alignment process in 100 s, which is much shorter than the traditional two-position method. The horizontal results are slightly worse, while the azimuth does have more excellent performance than the two-position method, supposing the sensor errors exhibited in Table 1. In consideration of the simulation condition, the improved method can fulfill the initial alignment totally in motion, and do have superior azimuth accuracy and a much shorter time, which markedly improves the maneuverability of the vehicle.

4.2. Vehicle Test

The vehicle test was conducted to validate the actual performance of the novel alignment mechanism. The navigation system in the experiment is manufactured by our own research group, and the parameters of inertial sensors are similar with Table 1. The GNSS receiver is purchased from NovAtel, including the GPS-702-GG and the OEM615 as receiver. The velocity accuracy of the receiver can achieve 0.03 m/s (RMS), and the horizonal position accuracy is 1.5 m (RMS). The odometer is purchased from the state enterprise and it was utilized in many fields. The output of the odometer is pulse signal, and the scale factor parameter is under 0.2%.
As attitude reference, the high-precision SINS was chosen to conduct the experiment, and the main characteristics are shown in Table 3. Considering the superior pose-maintaining performance in the navigation phase, the attitude output of the high-precision navigation system could be the reference of the novel method.
The experiment program can be described as follows:
(a)
The testing navigation system (simplified as System I) and the reference system (simplified as System II) were mounted in the testing car, as shown in Figure 8. The receiver of the GPS is integrated in the navigation system, the information from the GPS can be shared by System I and System II. The velocity information from the odometer is pulse signal, which can be collected by System I immediately.
(b)
The misalignment of the two systems could be obtained through optical sighting, as these two systems had installed the optical prism. The misalignment must be compensated when processing the output data.
(c)
System I and System II were launched simultaneously and with a strick on the mounting platform. System II would start its alignment and navigation procedure in normal mode. The testing carrier could moving out after the alignment of System II finished. All of the data from System I and System II, including the inertial sensor, the GPS, and the odometer, should be saved in the computer for further processing.
With the assistance of synchronous processing, the alignment errors of System I compared with System II are presented in Figure 9, and the installation error of the two systems was already taken into account. The experiment result indicates that the horizontal errors are about 0.04°, and the azimuth error is about 0.06°. Compared with the simulation results, it demonstrates that the performance of the horizontal attitude is worse, while the azimuth error is in accordance with the simulation. Further, the detailed graph shows that there is somewhat of an angle drift in the post-alignment period. The poor performance of pitch and roll, along with the drift of the attitude, could attribute to the misalignment between the odometer and the SINS. The wheel slip could also impact the velocity calculation, and therefore affect the alignment results. However, the azimuth result is far superior to the theoretical accuracy of the SINS with the applied gyros and the accelerometers. Furthermore, the alignment time is confined to 100 s, which is shorter than most marching alignment procedures, and especially important in motion operation for the vehicles.
The comparison with the traditional alignment method could be carried out with the offline sensors data. The most frequently used alignment method for marching carriers is the fine alignment based on Kalman filter along with the coarse inertial alignment. This method could be described as having two stages: coarse inertial alignment period, consuming about 3 min; fine alignment based on Kalman filter, continuing calculation until the convergence of the attitude finished. With the same sensor data, the alignment processing was conducted, and the alignment errors relating to System II are presented in Figure 10. The coarse alignment did have a relatively large misalignment, and the fine alignment could converge the attitude to a smaller value. The convergence time for the attitude is a long period. There must be ~300–400 s for the attitude to approach the final value. The alignment accuracy is about 0.08° for horizontal angles, and 0.8° for the azimuth, which is worse than the theoretical value because of the dynamics.
Another comparison is also carried out with the alignment method proposed in [27], where the alignment is accomplished through the velocity vector from the GPS and integrated accelerations. As stated in [27], the method is more suitable when the carriers are marching with an approximate circumference motion. While for linear moving with turning in the heading of the carrier, the performance of the method is not ideal: the alignment time is long with low precision. The experiment result for the azimuth is shown in Figure 11. The alignment result does vary in the alignment period, and the final error of the azimuth is about 0.5°.
With the testing result of the novel alignment method and the comparison with typical methods, it can be concluded that: the novel alignment method could achieve an initial attitude establishment; it could accomplish the alignment within 100 s, which is much shorter than the other methods; the horizontal accuracy could achieve 0.03°, and the azimuth precision could achieve 0.06°. As a result, the novel method is suitable for situations when the carriers are in motion and it requires that the odometer should be calibrated in advance. However, the changing in the carrier environment could impact the performance of the method, which is necessary for further researching.

5. Conclusions

Tackling the requirement of fast initial alignment for a marching vehicle, a novel alignment procedure is proposed for the SINS aided by the GPS and the odometer. The matrix C i 0 i b 0 , the core problem of the computing of the attitude matrix C b n , can be achieved through the nonlinear velocity vectors information from the GPS and the odometer at different moments, denoted as the multi-vector attitude determination. Comparing with the existing alignment methods, this approach does not use the measurement information of accelerometers, while the velocity of the vehicle in the n frame and the b frame is needed synchronously. Compared with the existing alignment methods for a marching vehicle, the initial alignment can be fulfilled in shorter time, and does have superior performance, especially with relatively low-precision inertial sensors. Besides, the only requirement of this method for a marching vehicle is a turning movement in the alignment stage, which can be easily realized in the real conditions. The qualitative relation between measurement errors and alignment errors is given, and the alignment results are mostly impacted by the velocity from the GPS and the odometer. To cope with the random noise of velocity from the GPS, an improved method is put forward based on the error equation. A comparison with the two-position alignment method is developed to demonstrate its performance. According to the theoretical analysis, the simulation, and the vehicle test, the initial alignment can be fulfilled in 100 s, and the alignment azimuth error is smaller than 0.06° with the condition that the bias and random errors of gyros are 0.1° and 0.01°, respectively. Compared with the two-position alignment method, the accuracy of the horizontal angle is slightly worse, while the azimuth is much better and the alignment time is shorter with the same condition of gyros and accelerometers, particularly low-precision inertial sensors. The novel alignment method proposed could be applied in many fields where the carriers need to fulfill the alignment in a marching process. The application could include chariots, warships, and submarines, which are always in moving condition as a form of protection.

Author Contributions

Longjun Ran, Lailiang Song and Chunxi Zhang conceived the idea. Longjun Ran and Lailiang Song conducted the simulation and the vehicle test. Longjun Ran wrote the paper. Lailiang Song reviewed and edited the manuscript. All authors read and approved this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gaho, A.A.; Qureshi, M.S. The generic investigation of initial alignment of SINS. J. Space Technol. 2012, 1, 11–16. [Google Scholar]
  2. Yu, H.; Wu, W.; Wu, M.; Yu, M.; Hao, M. Stochastic observability-based analytic optimization of SINS multiposition alignment. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2181–2192. [Google Scholar] [CrossRef]
  3. Wang, Y.; Sun, F.; Zhang, Y.; Liu, H.; Min, H. Central difference particle filter applied to transfer alignment for SINS on missiles. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 375–387. [Google Scholar] [CrossRef]
  4. Titterton, D.; Weston, J.L. Strapdown Inertial Navigation Technology; IET: Stevenage, UK, 2004. [Google Scholar]
  5. Jiang, Y.F. Error analysis of analytic coarse alignment methods. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 334–337. [Google Scholar] [CrossRef]
  6. Pecht, E.; Mintchev, M.P. Observability analysis for INS alignment in horizontal drilling. IEEE Trans. Instrum. Meas. 2007, 56, 1935–1945. [Google Scholar] [CrossRef]
  7. Jiancheng, F.; Sheng, Y. Study on innovation adaptive EKF for in-flight alignment of airborne POS. IEEE Trans. Instrum. Meas. 2011, 60, 1378–1388. [Google Scholar] [CrossRef]
  8. Wu, M.; Wu, Y.; Hu, X.; Hu, D. Optimization-based alignment for inertial navigation systems: Theory and algorithm. Aerosp. Sci. Technol. 2011, 15, 1–17. [Google Scholar] [CrossRef]
  9. Sun, F.; Sun, W.; Wu, L. Coarse alignment based on IMU rotational motion for surface ship. In Proceedings of the 2010 IEEE/ION Position Location and Navigation Symposium (PLANS), Indian Wells, CA, USA, 4–6 May 2010; pp. 993–998. [Google Scholar]
  10. Yu, H.; Zhu, H.; Gao, D.; Yu, M.; Wu, W. A stationary north-finding scheme for an azimuth rotational IMU utilizing a linear state equality constraint. Sensors 2015, 15, 4368–4387. [Google Scholar] [CrossRef] [PubMed]
  11. Lee, J.G.; Park, C.; Park, H.W. Multiposition alignment of strapdown inertial navigation system. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 1323–1328. [Google Scholar] [CrossRef]
  12. Dai, H.; Dai, S.; Cong, Y.; Zhao, G.; Wu, G. Rapid transfer alignment of laser SINS using quaternion based angular measurement. Optik-Int. J. Light Electron Opt. 2013, 124, 4364–4368. [Google Scholar] [CrossRef]
  13. Wendel, J.; Metzger, J.; Trommer, G.F. Rapid transfer alignment in the presence of time correlated measurement and system noise. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Providence, RI, USA, 16–19 August 2004; pp. 1–12. [Google Scholar]
  14. Hao, Y.; Xiong, Z.; Wang, W.; Sun, F. Rapid transfer alignment based on unscented Kalman filter. In Proceedings of the American Control Conference, Minneapolis, MN, USA, 14–16 June 2006. [Google Scholar] [CrossRef]
  15. Groves, P.D. Optimising the transfer alignment of weapon INS. J. Navig. 2003, 56, 323–335. [Google Scholar] [CrossRef]
  16. Liu, X.; Xu, X.; Liu, Y.; Wang, L. A fast and high-accuracy transfer alignment method between M/S INS for ship based on iterative calculation. Measurement 2014, 51, 297–309. [Google Scholar] [CrossRef]
  17. Chen, Y.; Zhao, Y. New rapid transfer alignment method for SINS of airborne weapon systems. J. Syst. Eng. Electron. 2014, 25, 281–287. [Google Scholar] [CrossRef]
  18. Sun, F.; Sun, W. Mooring alignment for marine SINS using the digital filter. Measurement 2010, 43, 1489–1494. [Google Scholar] [CrossRef]
  19. Lv, S.; Xie, L.; Chen, J. New techniques for initial alignment of strapdown inertial navigation system. J. Frankl. Inst. 2009, 346, 1021–1037. [Google Scholar]
  20. El-Sheimy, N.; Nassar, S.; Noureldin, A. Wavelet de-noising for IMU alignment. Aerosp. Electron. Syst. Mag. 2004, 19, 32–39. [Google Scholar] [CrossRef]
  21. Li, Z.; Wang, J.; Gao, J.; Li, B.; Zhou, F. A vondrak low pass filter for IMU sensor initial alignment on a disturbed base. Sensors 2014, 14, 23803–23821. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, X.; Zhao, Y.; Liu, X.; Yang, Y.; Song, Q.; Liu, Z. An improved self-alignment method for strapdown inertial navigation system based on gravitational apparent motion and dual-vector. Rev. Sci. Instrum. 2014, 85, 125108. [Google Scholar] [CrossRef] [PubMed]
  23. Li, Q.; Ben, Y.; Yang, J. Coarse alignment for Fiber Optic Gyro SINS with external velocity aid. Optik-Int. J. Light Electron Opt. 2014, 125, 4241–4245. [Google Scholar] [CrossRef]
  24. Liu, X.; Liu, X.; Song, Q.; Yang, Y.; Liu, Y.; Wang, L. A novel self-alignment method for SINS based on three vectors of gravitational apparent motion in inertial frame. Measurement 2015, 62, 47–62. [Google Scholar] [CrossRef]
  25. Ma, L.; Wang, K.; Shao, M. Initial alignment on moving base using GPS measurements to construct new vectors. Measurement 2013, 46, 2405–2410. [Google Scholar] [CrossRef]
  26. Gu, D.; El-Sheimy, N.; Hassan, T.; Syed, Z. Coarse alignment for marine SINS using gravity in the inertial frame as a reference. In Proceedings of the 2008 IEEE/ION Position, Location and Navigation Symposium, Monterey, CA, USA, 5–8 May 2008; pp. 961–965. [Google Scholar]
  27. Silson, P.M.G.; Jordan, S. A novel inertial coarse alignment method. In Proceedings of the AIAA Guidance, Navigation and Control Conference, Toronto, ON, Canada, 2–5 August 2010. [Google Scholar]
  28. Feng, S.; Tong, C. Accuracy analysis of coarse alignment based on gravity in inertial frame. Chin. J. Sci. Instrument. 2011, 32, 2409–2415. [Google Scholar]
  29. Han, S.; Wang, J. A novel initial alignment scheme for low-cost INS aided by GPS for land vehicle applications. J. Navig. 2010, 63, 663–680. [Google Scholar] [CrossRef]
  30. Chang, L.B.; He, H.Y.; Qin, F.J. In-Motion Initial Alignment for Odometer-Aided Strapdown Inertial Navigation System Based on Attitude Estimation. IEEE Sens. J. 2017, 17, 766–773. [Google Scholar] [CrossRef]
  31. Zhao, C.; Qin, Y.; Yan, G. On-the-move alignment for strap-down inertial navigation system. In Proceedings of the International Conference on Information and Automation, Changsha, China, 20–23 June 2008; pp. 1428–1432. [Google Scholar]
  32. Sun, F.; Lan, H.; Yu, C.; El-Sheimy, N.; Zhou, G.; Cao, T.; Liu, H. A robust self-alignment method for ship’s strapdown INS under mooring conditions. Sensors 2013, 13, 8103–8139. [Google Scholar] [CrossRef] [PubMed]
  33. Jianping, M.A. A new GPS aided initial alignment algorithm for SINS on moving base. Chin. J. Sens. Actuators 2010, 23, 1656–1661. [Google Scholar]
Figure 1. The i 0 frame and the n frame on the Earth.
Figure 1. The i 0 frame and the n frame on the Earth.
Sensors 18 00137 g001
Figure 2. Turning movement in the n frame.
Figure 2. Turning movement in the n frame.
Sensors 18 00137 g002
Figure 3. Attitude error for raw method: the green triangles are the calculation errors by Equation (31); the red squares denote the alignment errors.
Figure 3. Attitude error for raw method: the green triangles are the calculation errors by Equation (31); the red squares denote the alignment errors.
Sensors 18 00137 g003
Figure 4. Variation of attitude error with novel alignment method: (a) the error variation of pitch; (b) the error variation of roll; and (c) the error variation of yaw.
Figure 4. Variation of attitude error with novel alignment method: (a) the error variation of pitch; (b) the error variation of roll; and (c) the error variation of yaw.
Sensors 18 00137 g004
Figure 5. Variation of attitude error among two-position alignment stage: (a) the error variation of pitch; (b) the error variation of roll; and (c) the error variation of yaw.
Figure 5. Variation of attitude error among two-position alignment stage: (a) the error variation of pitch; (b) the error variation of roll; and (c) the error variation of yaw.
Sensors 18 00137 g005
Figure 6. Repeatability of the novel alignment method: (a) the error discreteness of pitch; (b) the error discreteness of roll; (c) the error discreteness of yaw.
Figure 6. Repeatability of the novel alignment method: (a) the error discreteness of pitch; (b) the error discreteness of roll; (c) the error discreteness of yaw.
Sensors 18 00137 g006
Figure 7. Repeatability of the two-position alignment method: (a) The error discreteness of pitch; (b) the error discreteness of roll; and (c) the error discreteness of yaw.
Figure 7. Repeatability of the two-position alignment method: (a) The error discreteness of pitch; (b) the error discreteness of roll; and (c) the error discreteness of yaw.
Sensors 18 00137 g007
Figure 8. Devices mounted in the testing car for alignment experiment. UPort 1650 (MOXA, Taiwan): USB-to-serial converter.
Figure 8. Devices mounted in the testing car for alignment experiment. UPort 1650 (MOXA, Taiwan): USB-to-serial converter.
Sensors 18 00137 g008
Figure 9. Attitude errors of the novel alignment method in a moving vehicle test: (a) the error curves of pitch; (b) the error curves of roll; and (c) the error curves of yaw.
Figure 9. Attitude errors of the novel alignment method in a moving vehicle test: (a) the error curves of pitch; (b) the error curves of roll; and (c) the error curves of yaw.
Sensors 18 00137 g009
Figure 10. Attitude errors of the traditional alignment method in moving vehicle test: (a) the error curves of pitch; (b) the error curves of roll; (c) the error curves of yaw.
Figure 10. Attitude errors of the traditional alignment method in moving vehicle test: (a) the error curves of pitch; (b) the error curves of roll; (c) the error curves of yaw.
Sensors 18 00137 g010
Figure 11. Yaw errors of the alignment method in [27].
Figure 11. Yaw errors of the alignment method in [27].
Sensors 18 00137 g011
Table 1. Sensor errors settings.
Table 1. Sensor errors settings.
GyroBias0.1 °/h
Random Walk 0.01   ° / h
Updating Frequency1 kHz
AccelerometerBias500 μg
Updating Frequency1 kHz
GPSVelocity accuracy0.03 m/s (RMS)
Horizontal position accuracy2 m (RMS)
Updating Frequency10 Hz
Odometer biasScale Factor0.2%
Updating Frequency10 Hz
Table 2. Statistics of attitude error for alignment and calculation.
Table 2. Statistics of attitude error for alignment and calculation.
StatisticsPitch (°)Roll (°)Yaw (°)
MeanAlignment error−0.0036−0.1516−0.0969
Calculation error−0.0036−0.1497−0.0065
StdAlignment error0.28230.00190.1574
Calculation error0.28199.8 × 10−40.3137
Table 3. Sensor errors of the high-precision strap-down inertial navigation system (SINS).
Table 3. Sensor errors of the high-precision strap-down inertial navigation system (SINS).
GyroBias0.01 °/h
Random Walk 0.001   ° / h
Updating Frequency1 kHz
AccelerometerBias50 μg
Updating Frequency1 kHz

Share and Cite

MDPI and ACS Style

Zhang, C.; Ran, L.; Song, L. Fast Alignment of SINS for Marching Vehicles Based on Multi-Vectors of Velocity Aided by GPS and Odometer. Sensors 2018, 18, 137. https://doi.org/10.3390/s18010137

AMA Style

Zhang C, Ran L, Song L. Fast Alignment of SINS for Marching Vehicles Based on Multi-Vectors of Velocity Aided by GPS and Odometer. Sensors. 2018; 18(1):137. https://doi.org/10.3390/s18010137

Chicago/Turabian Style

Zhang, Chunxi, Longjun Ran, and Lailiang Song. 2018. "Fast Alignment of SINS for Marching Vehicles Based on Multi-Vectors of Velocity Aided by GPS and Odometer" Sensors 18, no. 1: 137. https://doi.org/10.3390/s18010137

APA Style

Zhang, C., Ran, L., & Song, L. (2018). Fast Alignment of SINS for Marching Vehicles Based on Multi-Vectors of Velocity Aided by GPS and Odometer. Sensors, 18(1), 137. https://doi.org/10.3390/s18010137

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