Next Article in Journal
An Optical Sensory System for Assessment of Residual Cancer Burden in Breast Cancer Patients Undergoing Neoadjuvant Chemotherapy
Next Article in Special Issue
Prediction of Ground Wave Propagation Delay for MF R-Mode
Previous Article in Journal
Coordinate Fault Ride-Through Strategy for Connection of Offshore Wind Farms Using Voltage Source-Converter-Based High-Voltage Direct-Current Transmission under Single Polar Fault
Previous Article in Special Issue
Review of Shoreline Extraction Methods from Aerial Laser Scanning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Note on Coarse Alignment of Gyro-Free Inertial Navigation System

1
Advanced Technology Research Team, Hyundai Rotem, Uiwang-si 16082, Republic of Korea
2
Department of Electronics Engineering, Chungnam National University, Daejeon 34134, Republic of Korea
3
Advanced Technology Laboratory, Microinfinity Co., Ltd., Daejeon 34012, Republic of Korea
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(12), 5763; https://doi.org/10.3390/s23125763
Submission received: 1 May 2023 / Revised: 11 June 2023 / Accepted: 15 June 2023 / Published: 20 June 2023
(This article belongs to the Collection Navigation Systems and Sensors)

Abstract

:
In this note, the feasibility of initial alignment of a gyro-free inertial navigation system (GF-INS) is investigated. Initial roll and initial pitch are obtained using leveling of conventional INS since centripetal acceleration is very small. The equation for the initial heading cannot be used since the GF inertial measurement unit (IMU) cannot directly measure the Earth rate. A new equation is derived to obtain the initial heading from GF-IMU accelerometer outputs. Initial heading is expressed in the accelerometer outputs of two configurations, which satisfies a specific condition among 15 GF-IMU configurations presented in the literature. The initial heading error to arrangement and accelerometer error is quantitatively analyzed from the initial heading calculation equation of GF-INS and the initial heading error analysis of the general INS. The initial heading error is investigated when gyroscopes are used with GF-IMU. The results show that the initial heading error depends more on the performance of the gyroscope than that of the accelerometer, and the initial heading cannot be obtained within a practical error level by using only GF-IMU, even when an extremely accurate accelerometer is used. Therefore, aiding sensors have to be used in order to have a practical initial heading.

1. Introduction

The inertial navigation system (INS) continuously provides the attitude, velocity, and position of a vehicle from the specific force and angular velocity measured by an inertial measurement unit (IMU) [1,2,3]. Since the INS has a high output rate and wide bandwidth, it has been widely used as a navigation system of vehicles with high dynamics, such as rockets, guided missiles and aircraft [4].
The initial attitude, initial velocity, and initial position have to be known before the navigation is carried out since the INS calculates a navigation solution by integrating IMU measurements. The initial velocity and initial position have be externally provided, while the initial attitude can be determined in stationary by itself from the gravity and the Earth rate measured by the IMU [2]. This initial attitude determination without external aiding is called self-alignment [5]. In the gimbaled INS (GINS), the self-alignment is a procedure that aligns physically the platform with the navigation frame using gimbals. The self- alignment in the strap-down INS (SDINS) is a procedure to determine the initial attitude represented by the direction cosine matrix (DCM) between the body frame and the navigation frame or Euler angles, roll, pitch, and heading.
The GF-INS is an inertial navigation system that uses only accelerometers. Shuler first proposed in the 1960s the idea that measures the angular motion of a vehicle using the lever-arm effect of the accelerometer [6]. However, subsequent research was not followed since the quality of the accelerometer was not good enough to completely replace the gyroscope. With the development of micro-electro-mechanical system (MEMS) technology, nanotechnology, and cold atomic technology from the late 1990s, researchers have paid much attention to the GF-INS [7,8].
Before performing navigation, attitude should be initialized in the GF-INS as the conventional INS. Even though the self-alignment algorithm for the conventional INS is available, it cannot be used without modification since the measuring method of angular motion for a vehicle in the GF-INS is fundamentally different from that of the conventional INS.
Vaknin and Klein proposed a partial coarse alignment procedure for the roll and pitch only and developed closed-form expressions of the attitude error of the GF-INS coarse alignment [9]. They compared the partial coarse alignment performances of several GF configurations. Vaknin and Klein claimed that a useful initial heading cannot be obtained, regardless of the accelerometer grade, since the GF-INS cannot directly measure the Earth rate [9].
In this note, a coarse alignment of the GF-INS is presented, and the feasibility of the initial heading is evaluated. The initial heading expression in accelerometer outputs of the GF-IMU is derived from the initial heading expression in the coarse alignment of the conventional INS. The configurations of the GF-IMU, in which the initial heading expression can be obtained, are investigated for the configurations of the GF-IMU in the literature [9,10,11,12,13,14,15,16,17]. The initial heading error to arrangement and accelerometer error is quantitatively analyzed from the initial heading calculation equation of GF-INS and the initial heading error analysis of the general INS. The initial heading error is investigated when gyroscopes are used with GF-IMU.
The organization of this note is as follows: in Section 2, the coarse alignment for the conventional INS is introduced. In Section 3, the coarse alignment of the GF-INS is described. In Section 4, the coarse alignment of the GF-INS is described when gyroscopes are added to the GF-IMU. Finally, concluding remarks and further studies are presented in Section 5.

2. Coarse Alignment of INS

2.1. Leveling

Only gravity is measured by the accelerometer when the vehicle is stationary. The measured gravity is represented in Equation (1) between the body frame and the navigation frame [18,19,20].
f b = C n b f n = C 3 C 2 C 1 f n = 1 0 0 0 cos ϕ sin ϕ 0 sin ϕ cos ϕ cos θ 0 sin θ 0 1 0 sin θ 0 cos θ cos ψ sin ψ 0 sin ψ cos ψ 0 0 0 1 0 0 g = cos θ cos ψ cos θ sin ψ sin θ sin ϕ sin θ cos ψ cos ϕ sin ψ sin ϕ sin θ sin ψ + cos ϕ cos ψ sin ϕ cos θ cos ϕ sin θ cos ψ + sin ϕ sin ψ cos ϕ sin θ sin ψ sin ϕ cos ψ cos ϕ cos θ 0 0 g = g sin θ g sin ϕ cos θ g cos ϕ cos θ = f b x f b y f b z
where ϕ , θ , and ψ represent roll, pitch, and heading, respectively. The symbol g denotes gravity. The specific force is represented in the body frame and the navigation frame as f b and f n , respectively. C n b denotes the DCM between the navigation frame and the body frame. C 1 , C 2 , and C 3 correspond to DCMs between two frames when the frame is rotated ψ , θ , and ϕ around the yaw, pitch, and roll axes, respectively.
Roll can be obtained as Equation (2) from Equation (1) [18,19,20].
ϕ = tan 1 f b y f b z
where f b y and f b z are the y -axis and the z -axis component of the specific force, respectively. Pitch θ is obtained from Equation (3) [18,19,20].
θ = tan 1 f b x f b y 2 + f b z 2
where f b x is the x -axis component of the specific force.

2.2. Gyrocompassing

Only the Earth rate is measured by the gyroscope when the vehicle is stationary. The Earth rate is represented as Equation (4) in the body frame and in the navigation frame [18,19,20].
ω i e b = C n b ω i e n = C 3 C 2 C 1 Ω cos L 0 Ω sin L
where Ω and L are the magnitude of the Earth rate and the latitude, respectively. ω i e b and ω i e n are the Earth rate represented in the body frame and navigation frame, respectively. From Equation (4), the Earth rate ω i e r represented in the r frame, which is rotated from the navigation frame by the heading ψ , is given in Equation (5) [18,19,20].
ω i e r = C 1 ω i e n = Ω cos L cos ψ Ω cos L sin ψ Ω sin L = C 3 C 2 1 ω i e b = ω i e b x cos θ + ω i e b y sin ϕ sin θ + ω i e b z cos ϕ sin θ ω i e b y cos ϕ ω i e b z sin ϕ ω i e b x sin θ + ω i e b y sin ϕ cos θ + ω i e b z cos ϕ cos θ = ω i e r x ω i e r y ω i e r z
From Equation (5), the heading is obtained from Equation (6) [18,19,20].
ψ = tan 1 ω i e r y ω i e r x = tan 1 ω i e b y cos ϕ ω i e b z sin ϕ ω i e b x cos θ + ω i e b y sin ϕ sin θ + ω i e b z cos ϕ sin θ

3. Coarse Alignment of GF-INS

3.1. Roll and Pitch

The k -th accelerometer output of the GF-IMU with N accelerometers is given in Equation (7) [11].
y k = d k b T f k b = d k b T a b + ω i b b × ω i b b × r k b + ω ˙ i b b × r k b g b
where d k b denotes the sensing direction of the k -th accelerometer. f k b is the specific force of the k -th accelerometer represented in the body frame. a b is the acceleration at the center of the gravity of the vehicle represented in the body frame. ω i b b is the angular velocity of the body frame with respect to the inertial frame represented in the body frame. ω ˙ i b b is the angular acceleration of the body frame with respect to the inertial frame represented in the body frame. r k b is the position of the k -th accelerometer represented in the body frame. g b is the gravity vector represented in the body frame. If the vehicle is stationary on Earth, then a b = 0 , ω ˙ i b b = 0 , and ω i b b = ω i e b + ω e b b = ω i e b . Therefore, the k -th accelerometer output can be represented in Equation (8).
y k = d k b T f k b = d k b T ω i e b × ω i e b × r k b g b
In Equation (8), the magnitude of the centripetal acceleration is much less than the gravity as Equation (9) since the magnitude of the vector ω i e b , ω i e b , is 7.292115 × 1 0 5   rad / s [9].
ω i e b × ( ω i e b × r k b ) g b
As a result of this, Equation (8) can be represented in Equation (10) [9].
y k [ d k b ] T f k b = [ d k b ] T g b
The calculated specific force f ^ b from N accelerometers outputs can be obtained from Equation (11).
f ^ b = f ^ b x f ^ b y f ^ b z T = g ^ b = D T D 1 D T y ~
where D is the sensing direction matrix of the accelerometers in Equation (12).
D = d 1 b T d 2 b T d N b T T
y ~ is the accelerometer measurement vector in Equation (13).
y ~ = y ~ 1 y ~ 2 y ~ N T
Inserting the specific force in Equation (11) into Equations (2) and (3), roll and pitch can be obtained.

3.2. Heading

In order to obtain the heading, let us change Equation (6) to Equation (14).
ψ = tan 1 ω i e r y ω i e r x = tan 1 ω i e b x ω i e r y ω i e b x ω i e r x = tan 1 ω i e b x ω i e b y cos ϕ ω i e b x ω i e b z sin ϕ ω i e b x 2 cos θ + ω i e b x ω i e b y sin ϕ sin θ + ω i e b x ω i e b z cos ϕ sin θ
Equation (8) can be expressed as Equation (15).
y k + [ d k b ] T g b = [ d k b ] T f k b + [ d k b ] T g b = [ d k b ] T [ ω i e b × ( ω i e b × r k b ) ] = [ d k b ] T Ω i e b Ω i e b r k b = d k b x d k b y d k b z T 0 ω i e b z ω i e b y ω i e b z 0 ω i e b x ω i e b y ω i e b x 0 0 ω i e b z ω i e b y ω i e b z 0 ω i e b x ω i e b y ω i e b x 0 ( r k b ) x ( r k b ) y ( r k b ) z = d k b x d k b y d k b z T ω i e b y 2 ω i e b z 2 ω i e b x ω i e b y ω i e b x ω i e b z ω i e b x ω i e b y ω i e b x 2 ω i e b z 2 ω i e b y ω i e b z ω i e b x ω i e b z ω i e b y ω i e b z ω i e b x 2 ω i e b y 2 ( r k b ) x ( r k b ) y ( r k b ) z
where ω i e b x , ω i e b y , and ω i e b z are the x -, y , and z -axis component of ω i e b , respectively. Ω i e b is the skew symmetric matrix of the vector ω i e b . d k b x , d k b y , and d k b z are the x , y , and z -axis component of the sensing direction of the k -th accelerometer, respectively.
Rearranging Equation (15), Equation (16) can be obtained.
y k + d k b T g b = ( r k b ) y d k b x + ( r k b ) x d k b y ( r k b ) z d k b x + ( r k b ) x d k b z ( r k b ) z d k b y + ( r k b ) y d k b z ( r k b ) y d k b y ( r k b ) z d k b z ( r k b ) x d k b x ( r k b ) z d k b z ( r k b ) x d k b x ( r k b ) y d k b y T ω i e b x ω i e b y ω i e b x ω i e b z ω i e b y ω i e b z ω i e b x 2 ω i e b y 2 ω i e b z 2                                     = h k 1 h k 2 h k 3 h k 4 h k 5 h k 6 x
where x is given in Equation (17).
x = x 1 x 2 x 3 x 4 x 5 x 6 T = ω i e b x ω i e b y ω i e b x ω i e b z ω i e b y ω i e b z ω i e b x 2 ω i e b y 2 ω i e b z 2 T
If all the N accelerometer outputs of the GF-MU are put together, Equation (18) can be obtained from Equation (16).
a = H x = a 1 a 2 a N T = y 1 + d 1 b T g b y N + d N b T g b = h 1 1 h 1 2 h 1 3 h 1 4 h 1 5 h 1 6 h N 1 h N 2 h N 3 h N 4 h N 5 h N 6 x
From Equation (18), x can be obtained as Equation (19).
x = H T H 1 H T a
The heading can be expressed as Equation (20) from Equation (14) and x of Equation (19).
ψ = tan 1 x 1 cos ϕ x 2 sin ϕ x 4 cos θ + x 1 sin ϕ sin θ + x 2 cos ϕ sin θ

3.3. Accelerometer Configuration and Possiblity of Initial Heading Calculation

In this section, it is investigated, using Equation (20) in Section 3.2, whether the initial heading can be calculated for the accelerometer arrangements of GF-IMUs. It is known that at least six accelerometers are required to obtain the navigation information of a rigid body motion using accelerometers [10]. The navigation performance of a GF-INS is known to depend on the arrangement of the accelerometers, including the number, location, and sensing direction of the accelerometers. These results are found in references [9,10,11,12,13,14,15,16,17]. In this note, the possibility of the coarse alignment is investigated for the arrangements with 6 to 12 accelerometers, which have been considered in the existing research results.
Figure 1, Figure 2, Figure 3 and Figure 4 show 6 arrangements with 6 accelerometers, 1 arrangement with 7 accelerometers, 6 arrangements with 9 accelerometers, and 2 arrangements with 12 accelerometers of GF-IMUs, respectively. In Figure 1, Figure 2, Figure 3 and Figure 4, the black arrows indicate the x , y , and z -axes of the body frame. Blue dots indicate the positions of accelerometers and blue arrows the sensing directions of accelerometers.
If x can be obtained from Equation (18), the initial heading can be obtained from Equation (19). In order to obtain x , the rank of the matrix H T H has to be 6. The ranks of the matrix H T H for the arrangements in Figure 1, Figure 2, Figure 3 and Figure 4 are listed in Table 1.
It can be seen from Table 1 that the initial heading can be obtained from Equation (19) when the accelerometer arrangements are given in Figure 3e and Figure 4a. x 1 , x 2 , and x 4 for the initial heading calculation for the accelerometer arrangement in Figure 3e are given in Equations (21)–(23), respectively.
x 1 = 3 22 r 3 a 1 4 a 3 3 a 4 + 3 a 5 + 4 a 6 8 a 7 + 3 a 8
x 2 = 3 22 r 4 a 1 + 2 a 3 + 7 a 4 + 4 a 5 2 a 6 + 4 a 7 7 a 8
x 4 = 3 22 r 5 a 1 3 a 3 + 6 a 4 6 a 5 + 3 a 6 + 5 a 7 + 5 a 8
where a k ( k = 1 , , 8 ) is given in Equation (24). r denotes the distance from the center of gravity (accelerometer 9) to accelerometers 1 to 8.
a k = y k d k b T g b
x 1 , x 2 , and x 4 for the initial heading calculation for the accelerometer arrangement in Figure 4a are given in Equations (25)–(27), respectively.
x 1 = 1 2 r a 5 + a 7
x 2 = 1 2 r a 6 + a 10
x 4 = 1 2 r a 4 a 8 a 12
where a k ( k = 4 , 5 , 6 , 7 , 8 , 10 , 12 ) is given in Equation (24). r denotes the distance from the center of gravity (accelerometers 1 to 3) to accelerometers 4 to 12.

3.4. Initial Heading Error of Generic INS

The coarse alignment error characteristic of generic INS has been studied in many literature publications [1,3,5,21,22]. In this section, it is briefly introduced with the aim of deriving the initial heading error of the GF-INS.
The gyroscope output in the navigation frame can be expressed in Equation (28) when the vehicle is stationary.
ω ~ i e n = C b n ω i e b + δ ω i e b = C b n ω i e b + C b n δ ω i e b
where ω ~ i e n denotes the gyroscope output in the navigation frame and δ ω i e b the gyroscope output error in the body frame. If the roll, pitch, and heading are sufficiently small, C b n can be expressed in Equation (29).
C b n 1 δ ψ δ θ δ ψ 1 δ ϕ δ θ δ ϕ 1
where δ ϕ , δ θ , and δ ψ denote sufficiently small roll, pitch, and heading, respectively. In order to make the description simple, it is assumed that the body frame and the navigation frame are the same, i.e., ω i e b = ω i e n . Then, Equation (28) can be written in Equation (30).
ω ~ i e n = ω ~ i e n N ω ~ i e n E ω ~ i e n D T                         Ω cos L δ θ Ω sin L + δ ω i e b x δ ψ Ω cos L + δ ϕ Ω sin L + δ ω i e b y δ θ Ω cos L Ω sin L + δ ω i e b z
In order to obtain the initial heading, the east component of the gyroscope output ( ω ~ i e n E ) should be zero, as in Equation (31).
ω ~ i e n E = δ ψ Ω cos L + δ ϕ Ω sin L + δ ω i e b y = 0
δ ψ can be written in Equation (32) when δ ϕ = 0 in Equation (31).
δ ψ = δ ω i e b y Ω cos L
It can be seen from Equation (32) that the initial heading error in the generic INS is proportional to the gyroscope error of the east direction (the pitch axis direction) and inversely proportional to Ω cos L when the body frame is the same as the navigation frame. Therefore, the heading error is minimum when the vehicle is located at the equator, and the heading error increases as the vehicle becomes nearer to the pole ( L = ± 90 °).

3.5. Initial Heading Error of GF-INS

Let us now consider the initial heading error of the GF-INS with the previous initial heading error of the general INS.
The accelerometer output y ~ k of the GF-IMU can be represented in Equation (33).
y ~ k = y k + δ y k
where δ y k denotes the error of the k -th accelerometer output. If Equation (18) is used, a ^ can be represented in Equation (34) from the accelerometer output y ~ k in Equation (33) and g ^ b in Equation (11).
a ^ = y ~ 1 + d 1 b T g ^ b y ~ N + d N b T g ^ b = H x ^
where x ^ is the estimate of x from the accelerometer outputs. Error δ x of x ^ can be represented in Equation (35) from Equations (18) and (34).
δ x = x ^ x = H T H 1 H T a ^ a = H T H 1 H T y ˜ 1 + d 1 b T g ^ b y ˜ N + d N b T g ^ b y 1 + d 1 b T g b y N + d N b T g b   = H T H 1 H T δ y 1 + d 1 b T δ g b δ y N + d N b T δ g b = H T H 1 H T δ a 1 δ a N T = H T H 1 H T δ a
where δ a and δ g b denote the errors of a ^ and g ^ b , respectively. The covariance of δ x , Σ δ x can be obtained in Equation (36) from Equation (35).
Σ δ x = H T Σ δ a 1 H 1
where the covariance of δ a , Σ δ a is given in Equation (37) since cov δ a i , δ a j = 0 if i j .
Σ δ a = σ δ a 1 2 0 0 0 0 0 0 σ δ a N 2
where σ δ a k 2 denotes the variance of δ a k . Let us calculate the variances of δ x 1 , δ x 2 , and δ x 4 given in Equations (21)–(27). If accelerometers with the same specification are used, Equation (36) can be represented in Equation (38) since the condition σ δ a 1 2 = = σ δ a N 2 = σ δ a 2 is valid.
Σ δ x = σ δ a 2 H T H 1
The variances of δ x 1 , δ x 2 , and δ x 4 are listed in Table 2 for Figure 3e and Figure 4a from Equation (37).
The errors of x 1 ^ , x 2 ^ , and x 4 ^ in Equation (35) are given in Equation (39) when the body frame is the same as the navigation frame.
δ x 1 δ x 2 δ x 4 = x 1 ^ x 2 ^ x 4 ^ x 1 x 2 x 4 = ω ˜ i e b x ω ˜ i e b y ω ˜ i e b x ω ˜ i e b z ω ˜ i e b x 2 ω i e b x ω i e b y ω i e b x ω i e b z ω i e b x 2 = ( Ω cos L + δ ω i e b x ) δ ω i e b y Ω sin L + δ ω i e b z Ω cos L + δ ω i e b x 0 Ω 2 cos L sin L Ω cos L 2 Ω cos L δ ω i e b y Ω sin L δ ω i e b x + Ω cos L δ ω i e b z 2 Ω cos L δ ω i e b x
Equation (32) can be represented as Equation (40) using Equation (39).
δ ψ = δ ω i e b y Ω cos L = Ω cos L δ ω i e b y Ω cos L 2 = δ x 1 Ω cos L 2
It can be seen from Equation (39) that the heading error is proportional to δ x 1 , which is the error of ω i e b x ω i e b y obtained from Equation (19), and it is impossible to get the heading at the pole as with Equation (32). As shown in Table 2, δ x 1 depends on the grade and arrangement of the accelerometers.
Let us examine the heading error quantitatively by inserting the error specification of actual accelerometers into Equation (40). Table 3 shows the error specification of QA3000, a representative navigation-grade accelerometer manufactured by Honeywell Inc., USA [23] and Absolute Quantum Gravimeter, which is known as the most accurate accelerometer manufactured by Muquans Inc. (Talence, France) [24].
Let us obtain the initial heading error for the accelerometers in Table 3 when L = 0 ° and r = 1 m. In order to check the effect of the accelerometer error on the heading error, δ x can be expressed in the accelerometer error δ y as Equation (41) when g ^ b = g b in Equation (30).
δ x = H T H 1 H T δ y 1 δ y N T = H T H 1 H T δ y
Since σ δ a 2 = σ δ y 2 = σ δ y i 2 ( i = 1 , , N ), 1 σ of δ x 1 can be obtained when the specifications of the accelerometers in Table 3 are inserted into the expression in Table 2. By inserting these values into Equation (40), the 1 σ values of the initial heading errors are given in Table 4.
It can be seen from Table 4 that the practical initial heading of GF-INS cannot be obtained, even when the absolute quantum gravimeter, which is the most accurate accelerometer at the present technology level, and/or the navigation-grade accelerometer are used since the heading error is very large. For the 1 σ of the heading error in the arrangement of Figure 4a to be less than 5°, the accuracy of the accelerometer has to be less than σ δ y = 6.696 × 1 0 11 G. This value is 15 times more accurate than the absolute quantum gravimeter in Table 3. It can be seen from these results that it is difficult to have an initial heading with the required performance in the GF-INS, and these results are similar to those mentioned by Vaknin and Klein, in which the initial heading cannot be obtained [7,9].

4. Coarse Alignment of GF-INS with Gyroscope

According to the results in Section 3, even when the GF-IMU is configured with extremely accurate accelerometers, it is impossible to obtain an initial heading in the GF-INS. In this section, the heading error is examined when gyros are added to the GF-INS.
The gyroscope output g ~ n can be expressed in Equation (42) in the body frame.
g ~ n = g n + δ g n
where δ g n , ( n = x , y , o r z ) denotes the n -axis gyroscope output error in the body frame.
The relationship between the gyroscope output and x in Equation (18) is given in Equation (43).
b = b 1 b 2 b 3 b 4 b 5 b 6 T = g x g y g x g z g y g z g x 2 g y 2 g z 2 T = I 6 × 6 x
where I 6 × 6 denotes a 6 × 6 unit matrix. Inserting Equation (42) into Equation (43), the error of b ~ is given in Equation (44).
δ b = b ~ b = δ b 1 δ b 2 δ b 3 δ b 4 δ b 5 δ b 6 = g ~ x g ~ y g ~ x g ~ z g ~ y g ~ z g ~ x 2 g ~ y 2 g ~ z 2 g x g y g x g z g y g z g x 2 g y 2 g z 2 = g x δ g y + g y δ g x + δ g x δ g y g x δ g z + g z δ g x + δ g x δ g z g y δ g z + g z δ g y + δ g y δ g z 2 g x δ g x + δ g x 2 2 g y δ g y + δ g y 2 2 g z δ g z + δ g z 2
It can be observed from Equation (44) that the error of b ~ is proportional to the gyroscope outputs g x , g y , and g z .
For the simplicity of the description, consider the case that L = 0 ° and the body frame is the same as the navigation frame. In this case, Equation (44) can be represented in Equation (45) since g x g y g z T = Ω 0 0 T .
δ b = b ~ b = δ b 1 δ b 2 δ b 3 δ b 4 δ b 5 δ b 6 = Ω δ g y + δ g x δ g y Ω δ g z + δ g x δ g z δ g y δ g z 2 Ω δ g x + δ g x 2 δ g y 2 δ g z 2 Ω δ g y Ω δ g z 0 2 Ω δ g x 0 0
x can be estimated in Equation (46) using a ^ calculated from the accelerometer outputs and gyroscope outputs b ~ .
x ^ = H a u g T W H a u g 1 H a u g T W a ^ b ~ = H a u g T W H a u g 1 H a u g T W c
where W denotes the weighting matrix given in Equation (47).
W = Σ δ c 1 = diag σ δ a 2 , , σ δ a 2 N , σ δ b 1 2 , , σ δ b 6 2 1
where Σ δ c is the covariance of δ c , and σ δ b k 2 is the variance of δ b k . The measurement matrix H a u g is given in Equation (48).
H a u g = H a H b = h 1 1 h 1 6 h N 1 h N 6 I 6 × 6
where H a and H b are the measurement matrices of the accelerometer output and the gyroscope output, respectively.
The number of arrangements of gyroscopes, in which the gyroscopes are located at one or two axes among the roll, pitch, and yaw axes, is six in the body frame. One gyroscope can be located at the roll, pitch, or yaw axis, and two gyroscopes can be located at the roll and pitch axes, roll and yaw axes, or pitch and yaw axes. b ~ , W , and H b in Equation (46) for the gyroscope arrangements are listed in Table 5. The error covariance of the estimated x in Equation (46) is given in Equation (49).
Σ δ x = H a u g T W H a u g 1
For the accelerometer arrangement in Figure 4a, the heading error can be checked when the initial heading is obtained using the accelerometer outputs and the added gyroscope outputs. The variances of δ x 1 , δ x 2 , and δ x 4 from Equation (49) are listed in Table 6.
It can be seen from Table 6 that the variance of δ x 1 is the same as the case when only accelerometers are used, except for the case when the added gyroscopes are located at the roll and pitch axes. Since the heading error is determined by δ x 1 as shown in Equation (39), it can be seen that the initial heading accuracy is determined by the performance of the gyros located at the roll and pitch axes when the gyroscopes are added to the GF-INS.
Let us obtain the variance of δ x 1 to the grade of the gyroscope in order to have the initial heading error. First, the variance of δ b 1 in Equation (45) is given in Equation (50).
σ δ b 1 2 = var δ b 1 = var Ω δ g y = Ω 2 var δ g y = Ω 2 σ δ g y 2
When the gyroscopes with the same error specifications are located at the roll and pitch axes, the variance of δ x 1 is in Equation (51) since σ δ g x 2 = σ δ g y 2 = σ δ g 2 .
σ δ x 1 2 = σ δ a 2 σ δ b 1 2 σ δ a 2 + 2 r 2 σ δ b 1 2 = Ω 2 σ δ a 2 σ δ g 2 σ δ a 2 + 2 r 2 Ω 2 σ δ g 2
As shown in Equation (50), the variance of δ x 1 depends on the gyroscope performance ( σ δ g 2 ). Table 7 shows the bias error of Honeywell’s tactical-grade IMU, HG1700AG58, and navigation-grade gyro, GG1320AN [25,26]. It is known that the initial heading cannot be obtained when sub-tactical-grade gyroscopes that cannot accurately measure the Earth rate are used. For this reason, the specifications of Honeywell’s tactical- and navigation-grade gyros, which are widely used in guided weapons and navigation systems, are used in the error calculation.
Table 8 shows the 1 σ values of the heading error calculated using Equation (51) for the accelerometer error specifications from Table 3 and the gyroscope error specifications in Table 7 into Equation (51) when L = 0 ° and r = 1 m, as in Section 3.4. This is the same result as that of the generic INS gyro-compassing [1,4]. It can be seen from this result in Table 8 that the initial heading error is more dependent on the performance of the gyro than the performance of the accelerometer. This suggests that the practical initial heading cannot be obtained by using only GF-IMU without additional sensors, such as a gyroscope.

5. Concluding Remarks and Further Studies

In this note, the feasibility of obtaining the initial attitude of the GF-INS has been investigated by deriving the equation for the initial attitude from the accelerometer outputs. The initial heading is expressed in the accelerometer outputs of two configurations, which satisfies a specific condition among 15 GF-IMU configurations presented in the literature. In particular, it has been clarified that it is essential in the initial heading alignment as the general INS to accurately measure the Earth rate by analytically deriving the initial heading error due to the arrangement and error of the accelerometers. The initial heading errors have been checked for two kinds of GF-IMUs, which are composed of representative navigation-grade accelerometers and absolute quantum gravimeters, which are known as the most accurate accelerometers. The errors have been calculated by inserting error specifications into initial heading error equations. The results show that the initial heading cannot be obtained within a practical error level by using only GF-IMU due to limitations of the accelerometer production technology. In addition to this, the initial heading error has been checked when gyroscopes are used with the GF-IMU. It has been checked from this result that the initial heading accuracy depends more on the performance of the gyroscope than that of the accelerometer
It can be expected from the results that it is difficult to determine the initial heading using only the GF-IMU, even though the manufacturing technology of the accelerometer has improved. If the initial heading is needed in the GF-INS, non-inertial adding sensors, such as a magnetic compass or the attitude-determination GNSS receiver with two antennas, can be combined with the GF-IMU.

Author Contributions

Conceptualization, S.-H.O. and D.-H.H.; methodology, S.-H.O. and D.-H.H.; software, H.K.; validation, H.K. and J.-H.S.; writing—original draft preparation, H.K.; writing—review and editing, D.-H.H.; supervision, D.-H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No.2019R1A2C1088134).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Titterton, D.H.; Weston, J.L. Strapdown Inertial Navigation Technology, 2nd ed.; The Institution of Electrical Engineers: Piscataway, NJ, USA, 2004; pp. 17–301. [Google Scholar]
  2. Groves, P.D. Principles of GNSS, Inertial, and Multi Sensor Integrated Navigation Systems, 2nd ed.; Artech House: London, UK, 2013; pp. 15–160. [Google Scholar]
  3. Pitman, G.R. Inertial Guidance; John Wiley & Sons Inc: New York, NY, USA; Hoboken, NJ, USA, 1962; pp. 47–206. [Google Scholar]
  4. Bar-Itzhack, I.Y.; Berman, N. Control Theoretic Approach to Inertial Navigation Systems. AIAA J. Guid. Control. Dyn. 1998, 11, 237–245. [Google Scholar] [CrossRef]
  5. Britting, K.R. Inertial Navigation Systems Analysis; John Wiley & Sons Inc: New York, NY, USA; Hoboken, NJ, USA, 1971; pp. 153–209. [Google Scholar]
  6. Schuler, A.R.; Grammatikos, A.; Fegley, K.A. Measuring Rotational Motion with Linear Accelerometers. IEEE Trans. Aerosp. Electron. Syst. 1967, AES-3, 465–472. [Google Scholar] [CrossRef]
  7. Nusbaum, U.; Klein, I. Control Theoretic Approach to Gyro-Free Inertial Navigation Systems. IEEE Aerosp. Electron. Syst. Mag. 2017, 32, 168–170. [Google Scholar] [CrossRef]
  8. Pachter, M.; Welker, T.C.; Huffman, R.E. Gyro free INS theory. J. Inst. Navig. 2013, 60, 85–96. [Google Scholar] [CrossRef]
  9. Vaknin, E.; Klein, I. Coarse leveling of gyro-free INS. J. Gyroscopy Navig. 2016, 7, 145–151. [Google Scholar] [CrossRef]
  10. Chen, J.-H.; Lee, S.-C.; Debra, D.B. Gyroscope Free Strapdown Inertial Measurement Unit by Six Linear Accelerometers. J. Guid. Control. Dyn. 1994, 17, 286–290. [Google Scholar] [CrossRef]
  11. Padgaonkar, A.J.; Krieger, K.W.; King, A.I. Measurement of Angular Acceleration of a Rigid Body Using Linear Accelerometers. J. Appl. Mech. 1975, 42, 552–556. [Google Scholar] [CrossRef] [Green Version]
  12. Park, S.; Tan, C.W.; Park, S. A scheme for improving the performance of a gyroscope-free inertial measurement unit. Sens. Actuators A: Phys. 2005, 121, 410–420. [Google Scholar] [CrossRef]
  13. Fangjun, Q.; Xu, J.; Sai, J. A New Scheme of Gyroscope Free Inertial Navigation System Using 9 Accelerometers. In Proceedings of the 2009 International Workshop on Intelligent Systems and Applications, Wuhan, China, 23–24 May 2009. [Google Scholar]
  14. Shi, Z.; Cheng, Z.; Yang, J. Study of Independent Initial Alignment Based on Gyroscope Free Strapdown Inertial Navigation System. In Proceedings of the 2010 Chinese Control and Decision Conference, Xuzhou, China, 26–28 May 2010; pp. 2198–2202. [Google Scholar]
  15. Edwan, E.; Knedlik, S.; Loffeld, O. Constrained Angular Motion Estimation in a Gyro-Free IMU. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 596–610. [Google Scholar] [CrossRef]
  16. Wang, X.; Xiao, L. Gyroscope-reduced inertial navigation system for flight vehicle motion estimation. Adv. Space Res. 2017, 59, 413–424. [Google Scholar] [CrossRef]
  17. Santiago, A.A. Extended Kalman Filtering Applied to a Full Accelerometer Strapdown Inertial Measurement Unit. Master’s Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, September 1992. [Google Scholar]
  18. Farrell, J.A. Aided Navigation—GPS with High Rate Sensors; McGrawHill: New York, NY, USA, 2008; pp. 416–421. [Google Scholar]
  19. Noureldin, A.; Karamat, T.B.; Georgy, J. Fundamentals of Inertial Navigation Satellite-Based Positioning and Their Integration; Springer: New York, NY, USA, 2013; pp. 163–166. [Google Scholar]
  20. Kayton, M.; Fried, W.R. Avionics Navigation Systems, 2nd ed.; John Wiley & Sons Inc: New York, NY, USA, 1997; pp. 382–386. [Google Scholar]
  21. Yeon, F.-J. Error Analysis of Analytic Coarse Alignment Methods. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 334–337. [Google Scholar] [CrossRef]
  22. Silva, F.O. Generalized error analysis of analytical coarse alignment formulations for stationary SINS. Aerosp. Sci. Technol. 2018, 79, 500–505. [Google Scholar] [CrossRef]
  23. Honeywell. Q-FLEX QA-3000 ACCELEROMETER. Available online: https://aerospace.honeywell.com/us/en/products-and-services/product/hardware-and-systems/sensors/qa-3000-single-axis-quartz-accelerometer?utm_source=google&utm_medium=cpc&utm_campaign=23-aero-ww-dsa-products&utm_content=dyn-en-lp&gclid=EAIaIQobChMIw_z0z9Om_wIV18SWCh0cnwx8EAAYASAAEgKlDvD_BwE (accessed on 3 June 2023).
  24. Ménoret, V.; Vermeulen, P.; Moigne, N.L.; Bonvalot, S.; Bouyer, P.; Landragin, A.; Desruelle, B. Gravity measurements below 10−9 g with a transportable absolute quantum gravimeter. Sci. Rep. 2018, 8, 12300. [Google Scholar] [CrossRef] [PubMed]
  25. Honeywell. HG1700 Inertial Measurement Unit. Available online: https://aerospace.honeywell.com/us/en/products-and-services/product/hardware-and-systems/sensors/hg1700-inertial-measurement-unit (accessed on 3 June 2023).
  26. Honeywell. GG1320AN Digital Laser Gyro. Available online: https://aerospace.honeywell.com/us/en/products-and-services/product/hardware-and-systems/sensors/gg1320an-digital-ring-laser-gyroscope (accessed on 3 June 2023).
Figure 1. Six accelerometers arrangements for GF-INS. (a) Configuration 1, (b) Configuration 2, (c) Configuration 3, (d) Configuration 4, (e) Configuration 5, (f) Configuration 6.
Figure 1. Six accelerometers arrangements for GF-INS. (a) Configuration 1, (b) Configuration 2, (c) Configuration 3, (d) Configuration 4, (e) Configuration 5, (f) Configuration 6.
Sensors 23 05763 g001
Figure 2. Seven accelerometers arrangement for GF-INS. Configuration 7.
Figure 2. Seven accelerometers arrangement for GF-INS. Configuration 7.
Sensors 23 05763 g002
Figure 3. Nine accelerometers arrangements for GF-INS. (a) Configuration 8, (b) Configuration 9, (c) Configuration 10, (d) Configuration 11, (e) Configuration 12, (f) Configuration 13.
Figure 3. Nine accelerometers arrangements for GF-INS. (a) Configuration 8, (b) Configuration 9, (c) Configuration 10, (d) Configuration 11, (e) Configuration 12, (f) Configuration 13.
Sensors 23 05763 g003aSensors 23 05763 g003b
Figure 4. Twelve accelerometers arrangements for GF-INS. (a) Configuration 14, (b) Configuration 15.
Figure 4. Twelve accelerometers arrangements for GF-INS. (a) Configuration 14, (b) Configuration 15.
Sensors 23 05763 g004
Table 1. Rank of H T H .
Table 1. Rank of H T H .
Accelerometer ConfigurationNumber of Accelerometers Rank H T H
Figure 1a [9,10,12]63
Figure 1b [9]63
Figure 1c [9]65
Figure 1d [11]62
Figure 1e [11]63
Figure 1f [11]63
Figure 2 [16]74
Figure 3a [12]93
Figure 3b [9]93
Figure 3c [13,14]93
Figure 3d [13]93
Figure 3e [9]96
Figure 3f [9]95
Figure 4a [9,15]126
Figure 4b [17]123
Table 2. Variances of δ x 1 , δ x 2 and δ x 4 .
Table 2. Variances of δ x 1 , δ x 2 and δ x 4 .
ConfigurationFigure 3eFigure 4a
Variance of δ x k
σ δ x 1 2 9 11 r 2 σ δ a 2 1 2 r 2 σ δ a 2
σ δ x 2 2 21 22 r 2 σ δ a 2 1 2 r 2 σ δ a 2
σ δ x 4 2 45 44 r 2 σ δ a 2 3 4 r 2 σ δ a 2
Table 3. Specification of accelerometers.
Table 3. Specification of accelerometers.
Model NumberManufacturerError Spec. (1 σ Bias)Type
QA3000Honeywell Inc., USA 25 × 1 0 6 G (G = 9.8 m/s2)Quartz Pendulum
Absolute Quantum
Gravimeter
Muquans, France 1 0 9 GLaser-Cooled Atom
Table 4. The 1 σ value of initial heading error.
Table 4. The 1 σ value of initial heading error.
ConfigurationFigure 3eFigure 4a
Model Number
QA3000 2.39 × 1 0 6 ° 1.87 × 1 0 6 °
Absolute Quantum
Gravimeter
95.51°74.67°
Table 5. Value of b ~ W and Hb for gyro configurations.
Table 5. Value of b ~ W and Hb for gyro configurations.
Gyro Axis b ~ W H b
Roll b 4 diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 4 2 1 0 0 0 1 0 0
Pitch b 5 diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 5 2 1 0 0 0 0 1 0
Yaw b 6 diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 6 2 1 0 0 0 0 0 1
Roll and Pitch b 1 b 4 b 5 T diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 1 2 , σ δ b 4 2 , σ δ b 5 2 1 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0
Roll and Yaw b 2 b 4 b 6 T diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 2 2 , σ δ b 4 2 , σ δ b 6 2 1 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1
Pitch and Yaw b 3 b 5 b 6 T diag σ δ a 1 2 , , σ δ a N 2 , σ δ b 3 2 , σ δ b 5 2 , σ δ b 6 2 1 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1
Table 6. Variances of δ x 1 , δ x 2 and δ x 4 .
Table 6. Variances of δ x 1 , δ x 2 and δ x 4 .
Variance   of   δ x k σ δ x 1 2 σ δ x 2 2 σ δ x 4 2
Gyro Axis
Roll 1 2 r 2 σ δ a 2 1 2 r 2 σ δ a 2 3 σ δ a 2 σ δ b 4 2 3 σ δ a 2 + 4 r 2 σ δ b 4 2
Pitch 1 2 r 2 σ δ a 2 1 2 r 2 σ δ a 2 σ δ a 2 2 σ δ a 2 + 3 r 2 σ δ b 5 2 r 2 3 σ δ a 2 + 4 r 2 σ δ b 5 2
Yaw 1 2 r 2 σ δ a 2 1 2 r 2 σ δ a 2 σ δ a 2 2 σ δ a 2 + 3 r 2 σ δ b 6 2 r 2 3 σ δ a 2 + 4 r 2 σ δ b 6 2
Roll and Pitch σ δ a 2 σ δ b 1 2 σ δ a 2 + 2 r 2 σ δ b 1 2 1 2 r 2 σ δ a 2 σ δ a 2 σ δ b 4 2 2 σ δ a 2 + r 2 3 σ δ b 5 2 3 σ δ a 2 + r 2 3 σ δ a 2 σ δ b 4 2 + 3 σ δ a 2 σ δ b 5 2 + 4 r 2 σ δ b 4 2 σ δ b 5 2
Roll and Yaw 1 2 r 2 σ δ a 2 σ δ a 2 σ δ b 2 2 σ δ a 2 + 2 r 2 σ δ b 2 2 σ δ a 2 σ δ b 4 2 2 σ δ a 2 + 3 r 2 σ δ b 6 2 2 σ δ a 2 + r 2 3 σ δ a 2 σ δ b 4 2 + 3 σ δ a 2 σ δ b 6 2 + 4 r 2 σ δ b 4 2 σ δ b 6 2
Pitch and Yaw 1 2 r 2 σ δ a 2 1 2 r 2 σ δ a 2 σ δ a 2 σ δ a 2 2 + 2 r 2 σ δ a 2 σ δ b 5 2 + 2 r 2 σ δ a 2 σ δ b 6 2 + 3 r 4 σ δ b 5 2 σ δ b 6 2 r 2 2 σ δ a 2 2 + 3 r 2 σ δ a 2 σ δ b 5 2 + 3 r 2 σ δ a 2 σ δ b 6 2 + 4 r 4 σ δ b 5 2 σ δ b 6 2
Table 7. Specification of gyroscope.
Table 7. Specification of gyroscope.
GyroError Spec. (1 σ Bias)Type
RLG in HG1700AG58
(Honeywell Inc., USA)
1°/hTactical-Grade RLG
GG1320AN
(Honeywell Inc., USA)
0.003°/hNavigation-Grade RLG
Table 8. The 1 σ value of heading error.
Table 8. The 1 σ value of heading error.
GyroRLG in HG1700AG58GG1320AN
Accelerometer
QA30003.809°0.011°
Absolute Quantum
Gravimeter
3.804°0.011°
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, H.; Son, J.-H.; Oh, S.-H.; Hwang, D.-H. Note on Coarse Alignment of Gyro-Free Inertial Navigation System. Sensors 2023, 23, 5763. https://doi.org/10.3390/s23125763

AMA Style

Kim H, Son J-H, Oh S-H, Hwang D-H. Note on Coarse Alignment of Gyro-Free Inertial Navigation System. Sensors. 2023; 23(12):5763. https://doi.org/10.3390/s23125763

Chicago/Turabian Style

Kim, Heyone, Jae-Hoon Son, Sang-Heon Oh, and Dong-Hwan Hwang. 2023. "Note on Coarse Alignment of Gyro-Free Inertial Navigation System" Sensors 23, no. 12: 5763. https://doi.org/10.3390/s23125763

APA Style

Kim, H., Son, J. -H., Oh, S. -H., & Hwang, D. -H. (2023). Note on Coarse Alignment of Gyro-Free Inertial Navigation System. Sensors, 23(12), 5763. https://doi.org/10.3390/s23125763

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