Next Article in Journal
Centrifugal Microfluidic System for Nucleic Acid Amplification and Detection
Next Article in Special Issue
Inertial Sensor Error Reduction through Calibration and Sensor Fusion
Previous Article in Journal
Surface Plasmon Resonance Sensor Based on Ethylene Tetra-Fluoro-Ethylene Hollow Fiber
Previous Article in Special Issue
Research on the Rapid and Accurate Positioning and Orientation Approach for Land Missile-Launching Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Analytic Alignment Method for a SINS

1
School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
2
School of Computer Science and Electronic Engineering, University of Essex, Essex CO4 3SQ, UK
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(11), 27930-27953; https://doi.org/10.3390/s151127930
Submission received: 25 July 2015 / Revised: 3 October 2015 / Accepted: 28 October 2015 / Published: 4 November 2015
(This article belongs to the Special Issue Inertial Sensors and Systems)

Abstract

:
Analytic alignment is a type of self-alignment for a Strapdown inertial navigation system (SINS) that is based solely on two non-collinear vectors, which are the gravity and rotational velocity vectors of the Earth at a stationary base on the ground. The attitude of the SINS with respect to the Earth can be obtained directly using the TRIAD algorithm given two vector measurements. For a traditional analytic coarse alignment, all six outputs from the inertial measurement unit (IMU) are used to compute the attitude. In this study, a novel analytic alignment method called selective alignment is presented. This method uses only three outputs of the IMU and a few properties from the remaining outputs such as the sign and the approximate value to calculate the attitude. Simulations and experimental results demonstrate the validity of this method, and the precision of yaw is improved using the selective alignment method compared to the traditional analytic coarse alignment method in the vehicle experiment. The selective alignment principle provides an accurate relationship between the outputs and the attitude of the SINS relative to the Earth for a stationary base, and it is an extension of the TRIAD algorithm. The selective alignment approach has potential uses in applications such as self-alignment, fault detection, and self-calibration.

1. Introduction

If the Strapdown inertial navigation system (SINS) remains static on the ground with limited vibrations, the attitude of the SINS with respect to the Earth can be obtained directly from the gravity and rotational velocity vectors of the Earth if they are not collinear [1,2,3,4]. This approach is called analytic alignment (AA). The problem is about finding the transformation matrix from one coordinate system to the other when the components of two abstract vectors are given in two different coordinate systems. TRIAD [5,6] is an algorithm that does just that. A minimum of two non-coplanar vector pairs are required for a solution using TRIAD algorithm. The TRIAD algorithm is applied to the AA for a SINS, and then it is called the analytic coarse alignment method [1,2,3,4]. If the vibrations are so severe that the outputs of the gyroscopes are much larger than the rotational velocity of the Earth, AA cannot be performed. There are many approaches to the alignment problem including the inertial frame alignment method [7,8], Kalman fine alignment [9,10,11], and compass alignment [12,13]. All of these methods require part or all of the inertial navigation process, where the velocity and attitude of the SINS update with the time. Because the average velocity of the SINS is at zero for the base in the presence of vibrations, the attitude could then be computed. These alignment methods for the base in the presence of vibrations can be classified as inertial navigation computational alignment (INCA). However, INCA requires more alignment time and involves massive calculation compared to the AA method. Self-alignment methods can be divided into AA and INCA based on the principles used.
Because vibration is common in actual systems, INCA has been widely studied for several decades [7,8,9,10,12]. There have been few studies on AA, mainly because it has difficulty achieving alignment in the presence of vibrations. Reference [1] presents a basic AA method, analytic coarse alignment. One additional vector can be generated with the vector cross product between the gravity and rotational velocity vectors of the Earth. Then, a 3 × 3 matrix can be formed given the three vectors. In a SINS, the attitude information is stored either as a coordinate transformation matrix or as attitude angles (roll, pitch and yaw). The coordinate transformation matrix can be obtained directly through a one-time matrix calculation. Jiang [2] noted that there are three other vectors that can be generated from the cross-product of the gravity and rotational velocity vectors. Among these five vectors, four lie in the plane of the meridian and the other is perpendicular to this plane. The analytic coarse alignment problem is then to form a 3 × 3 matrix that consists of three linearly independent vectors. There are six possible sets that can be used to build a 3 × 3 matrix for the analytic coarse alignment. Two simple and significant sets are analyzed in [2] in detail. Although both methods are derived from the same measurements of the local gravity and rotational velocity vectors, their error formulations are not completely identical. The six analytic coarse alignment methods and a direct method for a SINS on stationary base were compared using error analyses in [3]. The error characteristics of the direct method can be evaluated because the direct method is equivalent to one of the six analytic coarse alignment methods.
It has been the standard practice for many years to use all six outputs from the IMU for self-alignment. A new AA method called selective alignment in ground platforms is proposed which requires only three outputs from the IMU and a few properties from the remaining outputs such as the sign and the approximate value. The selective alignment method has potential uses in problems such as optimal AA, fault detection, and self-calibration.
The selective alignment method shows that two non-coplanar vector pairs are redundant for determining the transformation matrix, however, a minimum of two non-coplanar vector pairs are required for a solution using TRIAD algorithm. The selective alignment may be applied in other attitude determination problems such as spacecraft attitude estimation where TRIAD algorithm is a basic method [6,14]. The formulas for the selective alignment are given through geometrical analysis in this paper. The principle could also be obtained through algebra analysis and we give a brief description about this in a recent conference paper [15]. After all, the geometrical analysis for the selective alignment shows us a much more intuitive and figurative relationship between the IMU outputs and its attitude relative to the Earth on the stationary base.
The remainder of the paper is organized as follows: We review the analytic coarse alignment method in Section 2. The basic principle of selective alignment is presented in Section 3. To prevent a zero value in the denominator or a negative value in the radicand at certain extreme attitudes, certain adjustments are required, the details of which are presented in Section 4. The results of simulations and experiments to test the validity and the use of the selective alignment method are presented in Section 5. Section 6 provides concluding remarks.

2. Analytic Coarse Alignment

The analytic coarse alignment is an application of the TRIAD algorithm to the AA for a SINS [2,6].
The process of finding the matrix using TRIAD is as follows [2,6,16]. Let u 1 and u 2 denote the column vectors whose elements are, respectively, the components of the two abstract vectors when resolved in one coordinate system (typically a body frame), and let v 1 and v 2 denote the column vectors whose elements are, respectively, the components of the abstract vectors when resolved in the other coordinate system (typically a reference frame). The reference-to-body coordinate transformation matrix A satisfies:
u 1 = A v 1 , u 2 = A v 2
The algorithm calls for the computation of the following column vectors in the body frame:
r 1 = u 1 / | u 1 | , r 2 = u 2 / | u 2 | , r 3 = ( u 1 × u 2 ) / | u 1 × u 2 | , r 4 = r 3 × r 1 , r 5 = r 3 × r 2
and the following corresponding column vectors in the reference frame:
s 1 = v 1 / | v 1 | , s 2 = v 2 / | v 2 | , s 3 = ( v 1 × v 2 ) / | v 1 × v 2 | , s 4 = s 3 × s 1 , s 5 = s 3 × s 2
of which four vectors are at the same plane except s 3 which is perpendicular to the plane. The coordinate matrix A can be obtained by
A = [ r 3 , r i , r j ] [ s 3 , s i , s j ] 1
where i , j = 1 , 2 , 4 , 5 , and i < j . There are six different sets to yield A .
When the SINS is static with respect to the Earth, there are only two inputs for the IMU: the gravity and rotational velocity vectors of the Earth, which are denoted by g and w i e , respectively, and are shown in Figure 1. A local ENU (East-North-Up) frame is used as the reference frame, or navigation frame, denoted by O x n y n z n . The body frame of the SINS is denoted by O x b y b z b .
Figure 1. The gravity and rotational velocity vectors in the navigation frame.
Figure 1. The gravity and rotational velocity vectors in the navigation frame.
Sensors 15 27930 g001
Neglecting sensor errors and base vibrations, the following measurements are available:
{ f b = C n b g n w b = C n b w i e n
where f b contains the outputs of a three-axis accelerometer with f b = [ f x , f y , f z ] T , w b contains the outputs of a three-axis gyroscope with w b = [ w x , w y , w z ] T , g n = [ 0 , 0 , g ] T and g is the gravitational acceleration, w i e n = [ 0 , w i e cos L , w i e sin L ] T with w i e and L denoting the rotation rate of the Earth and the latitude, respectively, and C n b is the navigation-to-body coordinate transformation matrix. From Equations (4) and (5), the following two equations can be obtained, and these are typically used in analytic coarse alignment methods [2,3].
C n b = [ f b , w b , f b × w b ] [ g n , w i e n , g n × w i e n ] 1
C n b = [ f b , f b × w b , ( f b × w b ) × f b ] [ g n , g n × w i e n , ( g n × w i e n ) × g n ] 1

3. Selective Alignment

If one of the IMU outputs, f y for example, is known, the angle between the y b -axis and g can be computed. If w y is also given, the angle between the y b -axis and w i e can be determined. As shown in Figure 2a, there are two possible solutions for the y b -axis because usually there will be two lines of intersection between the two conical surfaces.
After the y b -axis is obtained, the O x b y b z b frame can be obtained if one more axis is computed. For example, if w z is also known, the z b -axis will be on the conical surface with its center at w i e . In addition, because the z b -axis is normal to the y b -axis, the z b -axis must be in the plane that is normal to the y b -axis. There are usually two lines of intersection of a plane and a conical surface that represent the z b -axis, as shown in Figure 2b.
Given these facts, there are usually four possible results for the O x b y b z b frame given f y , w y , and w z .
Figure 2. Existence of the selective alignment solution: (a) Two lines that represent the y b -axis will be formed by the intersections of the two conical surfaces; (b) Two lines that represent the z b -axis will be formed by the intersection of a plane and a conical surface.
Figure 2. Existence of the selective alignment solution: (a) Two lines that represent the y b -axis will be formed by the intersections of the two conical surfaces; (b) Two lines that represent the z b -axis will be formed by the intersection of a plane and a conical surface.
Sensors 15 27930 g002

3.1. Derivation

Now let us derive mathematical expressions for the four solutions of O x b y b z b . As shown in Figure 3a, O x n y n z n denotes the navigation coordinates, as indicated by the subscript n , and O x b y b z b denotes the body coordinates, as indicated by the subscript b . Let O A = ( x x , y x , z x ) , O B = ( x y , y y , z y ) , and O C = ( x z , y z , z z ) be coordinates in the navigation frame, and let O A = ( 1 , 0 , 0 ) , O B = ( 0 , 1 , 0 ) , and O C = ( 0 , 0 , 1 ) be coordinates in the body frame. After ascertaining O x b y b z b , C n b can be obtained by
C n b = [ x x y x z x x y y y z y x z y z z z ]
Figure 3. (a) The definitions of O A , O B , and O C ; (b) O B , namely, the y b -axis, can be determined given w y and f y . Here, the y b -axis is not usually in the O y z plane.
Figure 3. (a) The definitions of O A , O B , and O C ; (b) O B , namely, the y b -axis, can be determined given w y and f y . Here, the y b -axis is not usually in the O y z plane.
Sensors 15 27930 g003
Assume three outputs of the IMU are given, which are f y , w y , and w z . If w y is given, as shown in Figure 3b, we have
cos B O D = w y w i e
Referring to Figure 3b, in the navigation frame the following equation holds:
O B O D = O B O D cos B O D
Because O B and O D are unit vectors and O D = w i e n w i e , substituting Equation (8) into Equation (9) gives
y y cos L + z y sin L = w y w i e
If f y is given, the angle between O B and the z -axis is known, which implies
z y = f y g
Combining Equations (10) and (11), and O B = 1 yields
{ y y cos L + z y sin L = w y w i e z y = f y g x y 2 + y y 2 + z y 2 = 1
Solving Equation (12) yields two sets of solutions:
{ x y = ± 1 y y 2 z y 2 y y = w y w i e cos L z y tan L z y = f y g
Similarly, if w z is given, we have
y z cos L + z z sin L = w z w i e
As shown in Figure 3a, because O C is the unit vector of the positive z b -axis and O C O B , the following equations hold
{ x z 2 + y z 2 + z z 2 = 1 x y x z + y y y z + z y z z = 0 y z cos L + z z sin L = w z w i e
Solving Equation (15) yields two sets of solutions
{ x z = y y y z + z y z z x y y z = w z w i e cos L z z tan L z z = b ± b 2 4 a c 2 a
where { a = ( x y 2 + y y 2 ) tan 2 L 2 y y z y tan L + x y 2 + z y 2 x y 2 b = 2 w z ( x y 2 + y y 2 ) tan L + 2 w z y y z y w i e x y 2 cos L c = y y 2 w z 2 x y 2 w i e 2 cos 2 L + w z 2 w i e 2 cos 2 L 1 .
Then, O A can be obtained as shown in Equation (17)
O A = O B × O C
Hence, from Equation (7), C n b can be obtained. There are four sets of solutions for C n b given f y , w y and w z .
Next, we consider the case where f y , w y and f x are given. Likewise, given f x rather than w z , we have
z x = f x g
Because O A is the unit vector of the positive x b -axis and O A O B , the following equations hold
{ x x 2 + y x 2 + z x 2 = 1 x y x x + y y y x + z y z x = 0 z x = f x g
Solving Equation (19) yields two sets of solutions
{ x x = b ± b 2 4 a c 2 a y x = x y x x + z y z x y y z x = f x g
where { a = x y 2 + y y 2 b = 2 x y z y f x g c = ( y y 2 + z y 2 ) f x 2 g 2 y y 2 .
Therefore, O C is given by
O C = O A × O B
Hence, C n b can be obtained from Equation (7). There are four sets of solutions for C n b given f y , w y , and f x .
The coordinate transformation matrix C n b has been derived in two typical cases. For both cases, one output of the three-axis accelerator and one output of the three-axis gyroscope, both in the same axis, should be selected. The case where the remaining measurement is another output of the three-axis gyroscope is denoted by FWW, e.g., f y , w y , and w z , and the case where the remaining measurement is another output of the three-axis accelerator is denoted by FWF, e.g., f y , w y , and f x .
There are six possible sets for FWW and six possible sets for FWF, as shown in Table 1.
Table 1. The six possible sets for FWW and FWF.
Table 1. The six possible sets for FWW and FWF.
123456
Output 1 f x f x f y f y f z f z
Output 2 w x w x w y w y w z w z
FWWOutput 3 w y w z w z w x w x w y
FWFOutput 3 f y f z f z f x f x f y
For FWW, the three selected outputs are denoted as f α , w α and w β , where the subscripts α , β and γ represent x , y or z and should be different. The body frame is given by
{ x α = ± 1 y α 2 z α 2 y α = w α w i e cos L z α tan L z α = f α g
{ x β = y α y β + z α z β x α y β = w β w i e cos L z β tan L z β = b ± b 2 4 a c 2 a
[ x γ y γ z γ ] = { [ x α , y α , z α ] T × [ x β , y β , z β ] T , if M [ x β , y β , z β ] T × [ x α , y α , z α ] T , if N
where a, b and c in Equation (23) are { a = ( x α 2 + y α 2 ) tan 2 L 2 y α z α tan L + x α 2 + z α 2 x α 2 b = 2 w β ( x α 2 + y α 2 ) tan L + 2 w β y α z α w i e x α 2 cos L c = y α 2 w β 2 x α 2 w i e 2 cos 2 L + w β 2 w i e 2 cos 2 L 1 , M in Equation (24) denotes { α = x β = y or { α = y β = z or { α = z β = x and N in Equation (24) denotes { α = y β = x or { α = z β = y or { α = x β = z .
For FWF, the three selected outputs are denoted as f α , w α and f β . The subscripts α , β and γ denote x , y or z and should be different. The body frame is given by Equations (22), (24), and (25):
{ x β = b ± b 2 4 a c 2 a y β = x α x β + z α z β y α z β = f β g
where a , b , and c in Equation (25) are { a = x α 2 + y α 2 b = 2 x α z α f β g c = ( y α 2 + z α 2 ) f β 2 g 2 y α 2 .

3.2. Choice of Solutions

After ascertaining O x b y b z b , C n b can be obtained from Equation (7). There are four possible solutions, but only one is appropriate. Given the four possible solutions of C n b , the vectors f b and w b can be computed from Equation (5). Comparing the calculated values of f b and w b with the original six outputs of the IMU, the appropriate solution can be easily ascertained amongst the four possible solutions.
To demonstrate the method for choosing the correct solution among the four possible solutions, one group of simulated outputs are used with L = 32.026372 ° , g = 9.8  m/s 2 , and w i e = 7.2921158 × 10 5  rad/s . f y , w y , and w z are selected to compute the attitude through the selective alignment method. Given the four possible solutions of the attitude (namely, C n b ), the four possible sets of the six outputs can be computed from Equation (5), which are shown in Table 2.
Table 2. The four possible sets of six inertial measurement unit (IMU) outputs calculated with f y , w y , and w z used for selective alignment according to Equations (22)–(24).
Table 2. The four possible sets of six inertial measurement unit (IMU) outputs calculated with f y , w y , and w z used for selective alignment according to Equations (22)–(24).
Angular Velocity ( 10 4 rad / s )Specific Force ( m / s 2 )
w x w y w z f x f y f z
1−0.25500.67820.0824−3.14973.35188.6536
20.25500.67820.08247.61923.3518−5.1724
30.25500.67820.08243.14973.35188.6536
4−0.25500.67820.0824−7.61923.3518−5.1724
Original outputs−0.25500.67820.0824−3.14973.35188.6536
It is easy to discern that the attitude for Set 1 in Table 2 is the correct solution by comparing the calculated and measured values of the outputs using one of the following sets of criteria:
  • The signs of f x and f z
  • The signs of w x and f z
  • The approximate value of f x .
The values compared will depend on the outputs selected. Through the selective alignment method, the attitude can be determined from any three outputs of the IMU in Table 1 in addition to at least one force component and one angular velocity component in the same axis. Some limited information from the remaining three outputs such as the signs and approximate value is also essential. Furthermore, if, for example, f x , f y , and w z are known, f z can be calculated from Equation (26). Then, the attitude by the selective alignment method can be obtained.
{ f x 2 + f y 2 + f z 2 = g 2 w x 2 + w y 2 + w z 2 = w i e 2

4. Exception in Selective Alignment

In practice, there are errors in the IMU outputs such as bias, and vibrations in the base will also affect the IMU outputs and thus the alignment. So, it is possible that x α in Equation (22) will have an imaginary value, and thus there is no solution. Furthermore, there are singularities in Equations (23) and (25) when x α and y α equal zero, respectively. Thus, the problems are divided into two classes, those in which the denominator should be nonzero and those in which the radicand should be nonnegative. We consider the denominator problem first.

4.1. Nonzero Denominator

For FWW, if x α equals zero, from Equation (15) we have
{ x β 2 + y β 2 + z β 2 = 1 y α y β + z α z β = 0 y β cos L + z β sin L = w z w i e
Solving Equation (27) yields
{ x β = ± 1 y β 2 z β 2 y β = w z w i e cos L z β tan L z β = y α w β w i e ( y α sin L z α cos L )
In Equation (28), ( y α sin L z α cos L ) must be nonzero. If ( y α sin L z α cos L ) = 0 , then because x α = 0 we have [ x α , y α , z α ] T = [ 0 , cos L , sin L ] T or [ x α , y α , z α ] T = [ 0 , cos L , sin L ] T , which indicates that the α b -axis coincides with w i e . Thus, the β -axis should be in the plane that is perpendicular to w i e . Because w β was also selected and must be zero, the β -axis must also be in the plane that is perpendicular to w i e . Two coincident planes have an infinite number of lines of intersection, so only the α b -axis can be obtained in this particular case.
For FWF, if y α is equal to zero, from Equation (19) we have
{ x β 2 + y β 2 + z β 2 = 1 x α x β + z α z β = 0 z β = f β g
Solving Equation (29) yields
{ x β = z α z β x α y β = ± 1 x β 2 z β 2 z β = f β g
In Equation (30), x α should be non-zero. If x α = 0 , then we have [ x α , y α , z α ] T = [ 0 , 0 , 1 ] T or [ x α , y α , z α ] T = [ 0 , 0 , 1 ] T , which indicates that the α -axis coincides with g . Thus, the β -axis should be in the plane that is perpendicular to g . Because f β was selected and must be zero, the β -axis must also be in the plane that is perpendicular to g . Two coincident planes have an infinite number of lines of intersection, so only the a b -axis can be obtained in this particular case.

4.2. Nonnegative Radicand

To ensure that the value of x α in Equation (22) is not imaginary, ( y α 2 + z α 2 ) should be less than 1. However, it is possible that f α can be slightly greater than g because of sensor errors or base vibrations. And, errors in w α may also result in a negative radicand. Constructing algebraic rules to prevent a negative value of the radicand is complicated and abstract. Graphically, a negative value in the radicand occurs when the two conical surfaces in Figure 2a do not intersect or the plane and the conical surface in Figure 2b do not intersect. Therefore, negative values in the radicand can be avoided if these surfaces intersect.
If the angle between the α b -axis and g is η , then cos η = f α g . Given w α and denoting the angle between the α b -axis and w i e as θ , we have cos θ = w α w i e . Assuming η is accurate, θ must lie in a certain range that depends on η , or the two conical surfaces shown in Figure 2a cannot intersect. Referring to Figure 4, given η , let A O P i = η . Then, θ must satisfy B O P i θ B O P i , i = 1 , 2 , 3 . The range of θ based on η is given by
{ π 2 L η θ π 2 L + η 0 η < π 2 L π 2 + L + η θ π 2 L + η π 2 L η < π 2 + L π 2 + L + η θ 3 π 2 + L η π 2 + L η π
Figure 4. The critical point for the intersection of the two conical surfaces, where P i is the reflection of P i about the z n -axis and O C is the reflection of O B about the y n -axis.
Figure 4. The critical point for the intersection of the two conical surfaces, where P i is the reflection of P i about the z n -axis and O C is the reflection of O B about the y n -axis.
Sensors 15 27930 g004
Figure 5a shows the range of θ for various values of η from 0 to π according to Equation (31).
Figure 5. (a) The range of θ as a function of η , where the point ( η , θ ) must be within the rectangle; (b) The range of w α as a function of f α , where the point ( f α , w α ) must be within the ellipse.
Figure 5. (a) The range of θ as a function of η , where the point ( η , θ ) must be within the rectangle; (b) The range of w α as a function of f α , where the point ( f α , w α ) must be within the ellipse.
Sensors 15 27930 g005
f α and w α can be obtained through f α = g cos η and w α = w i e cos θ . In this case, g and w i e can be considered constants. Given η and θ , we can plot the range of w α as a function of f α . It follows that ( f α , w α ) must be within the ellipse shown in Figure 5b to ensure that the two conical surfaces in Figure 2a intersect and the radicand in Equation (22) will not be negative.
After ascertaining the α b -axis, given another w β or f β , the β b -axis can be obtained. As shown in Figure 2b, the β b -axis is formed by the intersection of the conical surface and the plane. As shown in Figure 6, M N is a plane normal to the α b -axis. The conical surfaces A O A and C O C do not intersect the plane M N , but the conical surface B O B does. It follows that only if B is in the area between O M and O N will the conical surface B O B intersect the M N plane. In other words, the conical surface B O B is defined by w β or f β , which should satisfy
{ FWF : f α 2 + f β 2 g 2 FWW : w α 2 + w β 2 w i e 2
If Equation (32) holds, the conical surface resulting from f β or w β intersects the plane normal to the α b -axis and the radicand in Equations (23) and (25) will not be negative.
Figure 6. Intersection with the plane normal to the α b -axis and g ( w i e ) . The conical surfaces A O A , B O B , C O C are centered about g ( w i e ) and the points M , N are reflections of M , N , respectively, about g ( w i e ) .
Figure 6. Intersection with the plane normal to the α b -axis and g ( w i e ) . The conical surfaces A O A , B O B , C O C are centered about g ( w i e ) and the points M , N are reflections of M , N , respectively, about g ( w i e ) .
Sensors 15 27930 g006

4.3. Exception Handling in the Algorithm

To ensure that the selective alignment method performs robustly, various exceptions must be addressed. In the first stage, where the α -axis is determined, if | f α | is slightly greater than g because of bias in the accelerometer or base vibrations when the α -axis is nearly parallel to g , f α should be substituted with f α | f α | g . Similarly, if | w α | is slightly greater than w i e because of drift in the gyroscope or base vibrations when the α -axis is nearly parallel to w i e , w α should be substituted with w α | w α | w i e . In addition, ( f α , w α ) should be within the area shown in Figure 5b; otherwise, f α is substituted with f α + Δ f α or w α is substituted with w α + Δ w α to make that ( f α , w α ) lie in the boundary of the ellipse shown in Figure 5b. In the second stage, where β -axis is obtained, f β or w β must satisfy Equation (32); if not, the algorithm substitutes for f β or w β so that the equality holds.
In practice, the algebraic operations may not result in values that are identically zero because of the limits of floating-point operations on the computer, so values within a certain tolerance around zero should be considered zero in the algorithm. We use 10 12 < x < 10 12 to judge if x is equal to zero.
Usually, four results can be obtained from the three selective outputs based on the selective alignment. There might be only two results or one result when there is only one intersection in Figure 2. If no intersection occurs because of bias errors or vibration, the value of the relevant selected output is changed to ensure an intersection. Because all six original outputs are given, the correct value of the attitude is picked out by minimizing the following cost function for ease of computer programming, although two signs or one approximate value from the remaining outputs could be enough to obtain the right result, as explained in Section 3.2.
Figure 7. Flowchart for the selective alignment algorithm.
Figure 7. Flowchart for the selective alignment algorithm.
Sensors 15 27930 g007
δ = | w b C n b w i e n w i e | + | f b C n b g n g |
where w b and f b are the original outputs of the IMU, C n b is computed from Equation (7), which may have four different values, and | | is the norm operator.
To handle exceptions, a binary number with five bits, b 5 b 4 b 3 b 2 b 1 , is used. If there are no exceptions, b 5 b 4 b 3 b 2 b 1 = 00000 ; otherwise,
  • b 1 = 1 if f α > g
  • b 2 = 1 if w α > w i e
  • b 3 = 1 if there is no intersection with the α b -axis
  • b 4 = 1 if in FWW, w α 2 + w β 2 > w i e 2 , or if in FWF, f α 2 + f β 2 > g 2
  • b 5 = 1 if in FWW, α b -axis w i e , or if in FWF, α b -axis g (The symbol “ ” denotes “approximate to”)
The exceptions rarely happen unless the SINS is level (one of the body axes is along the direction of the gravity) or one of the body axes is parallel to the axis of rotation of the Earth. If b 5 = 1 , the attitude cannot be obtained. The other four exceptions arise because of measurement errors or base vibrations, which do not significantly affect the selective alignment after corresponding exception handling. Figure 7 shows a flowchart for the selective alignment algorithm.

5. Simulations and Experiments

Simulations and experiments were conducted to validate the performance of the proposed selective alignment method.

5.1. Static Alignment Example

Similar to the coordinate transformation matrix C n b , the Euler angles, i.e., roll, pitch, and yaw, can also be used to express the attitude [4].
Simulations were conducted to test the selective alignment method. Because the alignment is usually performed in a short period of time, only the bias and noise in the IMU were considered. With the SINS remaining static on the ground, the AA methods including the selective alignment method and the analytic coarse alignment method were used to obtain the attitude angles. The values of the various parameters were chosen as follows:
Simulation conditions for the static alignment:
L = 32 ° , g = 9.8 m / s 2 , w i e = 7.2921158 × 10 5 rad / s , [ roll , pitch , yaw ] T = [ 10 °, 10 °, 10 °] T ,
[ x , y , z ] T = [ 10 4 , 10 4 , 10 4 ] T × g , [ ε x , ε y , ε z ] T = [ 0.01 , 0.01 , 0.01 ] T   ° /h ,
[ σ f x , σ f y , σ f y ] T = [ 5 × 10 5 , 5 × 10 5 , 5 × 10 5 ] T × g , [ σ w x , σ w y , σ w z ] T = [ 0.005 , 0.005 ,0.005] T   ° /h
where i is the bias for the i -axis accelerometer, ε i is the drift (assumed constant) for the i -axis gyroscope, σ f i is the standard deviation of the white noise for the i -axis accelerometer, and σ w i is the standard deviation of the white noise for the i -axis gyroscope. The sampling rate was 100 Hz, and the IMU outputs were sampled for one minute. The average values were computed and used to test the AA method. There were twelve sets as a result of the various possibilities in the selective alignment method. The alignment results are shown in Table 3.
Table 3. The attitude angles obtained using analytic alignment (AA) methods in the simulation.
Table 3. The attitude angles obtained using analytic alignment (AA) methods in the simulation.
f x w x f y w y f z w z Roll ( ° )Pitch ( ° )Yaw ( ° ) b 5 b 4 b 3 b 2 b 1
Selective Alignment9.994310.00589.957200000
9.997010.09109.972500000
9.99339.97329.951400000
9.99319.96649.950200000
9.994310.00589.748500000
10.494710.00589.748500000
9.960610.00589.748500000
9.938610.00589.748500000
9.99339.97329.898900000
10.00159.96499.946400000
9.960610.00589.710600000
9.965110.00139.736500000
Analytic Coarse AlignmentEquation (6a)9.993310.00589.9451-
Equation (6b)9.993310.00589.9445-
True attitude101010-
Where the contents from left to right in the six grids denote f x , w x , f y , w y , f z , and w z , respectively. The corresponding contents are selected when being painted.
Next, assume [ roll , pitch , yaw ] T = [ 0 ° , 0 ° , 0 ° ] T , which is more likely to occur in practice, and the values of the other parameters were left unchanged. The results are shown in Table 4.
Table 4. The attitude angles obtained using AA methods in the simulation.
Table 4. The attitude angles obtained using AA methods in the simulation.
f x w x f y w y f z w z Roll ( ° )Pitch ( ° )Yaw ( ° ) b 5 b 4 b 3 b 2 b 1
Selective Alignment−0.00580.0058359.958800000
−0.00580.0715359.958800000
−0.00580359.958801000
−0.0058−0.0451359.958800000
−0.00580.0058000100
−0.07170.0058000100
00.0058001100
00.0058001100
−0.0058−0.04477.339710101
0−0.0451359.955200101
0.04470.0058262.599810101
0−0.0451001101
Analytic Coarse AlignmentEquation (6a)−0.00580.0058359.9588-
Equation (6b)−0.00580.0058359.9588-
True attitude000-
From Table 4, it can be observed that intolerant errors of the yaw angle occur due to the b 5 = 1 exception in ( f z , w z , f x ) and ( f z , w z , f y ) for the selective alignment because they are quite close to the extreme state α b -axis g in the FWF mode where the β axis is uncertain when the SINS is level, as explained in Section 4.1. The other exceptions except the b 5 = 1 exception do not significantly affect the selective alignment. From Table 3, no exception occurs, which indicates that the exceptions rarely happen unless the SINS is level or one of the body axes is parallel to the axis of rotation of the Earth. In general, the selective alignment and the analytic coarse alignment could achieve consistent alignment accuracy, and the alignment accuracy would be decreased when two gyroscope outputs are selected for the selective alignment.

5.2. Fault Detection Example

The selective alignment method may be potentially an effective approach for sensor fault detection such as that prior to take-off in aircraft. That is, by determining the alignment with independent sensors, individual sensor failures could be identified on the ground. A simulation was performed with self-alignment and assuming a 0.1 % -level IMU. The accuracy of the z -axis gyroscope was assumed to be degraded.
Simulation conditions for fault detection:
[ x , y , z ] T = [ 10 4 , 10 4 , 10 4 ] T × g , [ ε x , ε y , ε z ] T = [ 0.01 , 0.01 , 1 ] T   ° /h ,
[ σ f x , σ f y , σ f y ] T = [ 5 × 10 5 , 5 × 10 5 , 5 × 10 5 ] T × g , [ σ w x , σ w y , σ w z ] T = [ 0.005 , 0.005 , 0.005 ] T   ° /h ,
[ roll , pitch , yaw ] T = [ 0 ° , 0 ° , 30 ° ] T
Figure 8. The alignment results using selective alignment for fault detection.
Figure 8. The alignment results using selective alignment for fault detection.
Sensors 15 27930 g008
The IMU outputs were sampled for ten seconds. The average values were computed, and these are denoted as w b and f b . The difference between the norm of w b and w i e was 0.556   ° /h , and the difference between the norm of f b and g was 1.0121 × 10 4 g . It is obvious that several of the outputs of the three-axis gyroscope are faulty. The selective alignment method can be further used to identify the erroneous output. The alignment results are shown in Figure 8. Because the platform was level, the two sets of outputs ( f z , w z , f x ) and ( f z , w z , f y ) were not selected for use in the selective alignment method in case of the b 5 = 1 exception.
Because the outputs of the three-axis accelerometer are correct, the roll and pitch angles obtained through the analytic coarse alignment method are rational. From Figure 8, once w z is selected for the selective alignment method, the roll angle or the pitch angle will have obvious errors. It can be observed that the z -axis output of the three-axis gyroscope is severely affected.

5.3. Self-Calibration Example

An error in the z -axis gyroscope will cause errors mainly in the yaw angle in the inertial navigation process. To reach a balance between cost and accuracy, an IMU may include one high-accuracy z -axis gyroscope and two other less accurate gyroscopes. Using the selective alignment method with the outputs of the higher-accuracy z -axis gyroscope and the accelerometers, the attitude can be obtained. Then, the errors in the other, lower-accuracy gyroscopes can be compensated.
A simulation was performed to test the validity of the self-calibration procedure. The SINS was assumed to be static on the ground, and the SINS was assumed not to be level to avoid the b 5 = 1 exception.
Simulation conditions for the self-calibration:
[ x , y , z ] T = [ 10 4 , 10 4 , 10 4 ] T × g , [ ε x , ε y , ε z ] T = [ 0.01 , 0.01 , 1 ] T   ° /h ,
[ σ f x , σ f y , σ f y ] T = [ 5 × 10 5 , 5 × 10 5 , 5 × 10 5 ] T × g , [ σ w x , σ w y , σ w z ] T = [ 0.05 , 0.05 ,0.005] T   ° /h ,
[ roll , pitch , yaw ] T = [ 20 °, 20 °,20°] T
Table 5. The AA results for the self-calibration test.
Table 5. The AA results for the self-calibration test.
f x w x f y w y f z w z Roll ( ° )Pitch ( ° )Yaw ( ° )
Selective Alignment19.994320.006119.5203
20.355722.532620.5163
19.992219.990019.5140
20.000320.051619.5379
19.994320.006118.5666
27.982320.006118.5666
19.976020.006118.5666
20.665020.006118.5666
19.992219.990019.8515
19.894820.086819.5674
19.976020.006119.8044
19.645720.330018.8464
Analytic Coarse AlignmentEquation (6a)19.992220.006119.4363
Equation (6b)19.992220.006119.4388
True attitude202020
The IMU outputs were sampled for one minute. The average values were computed and then used for the compensation. The alignment results are given in Table 5. It can be observed that the attitude for the selective alignment method using ( f z , w z , f x ) produced the best result. Having computed C n b through the ( f z , w z , f x ) option of selective alignment, the constant drifts of the two lower-accuracy gyroscopes were obtained from Equation (34). The drift rates were ε ^ x = 0.0694   ° /h and ε ^ y = 0.0894   ° /h , and the drift errors of the two lower-accuracy gyroscopes are well compensated without multi-position.
ε ^ = w b C n b w i e n

5.4. Alignment in the Static Vehicle

The experiment was performed at the National University of Defense Technology (NUDT). The SINS used in the tests, shown in Figure 9, was built by NUDT and includes three 90-type mechanically dithered ring laser gyroscopes. The accuracy of the gyroscope is 0.01   ° /h , and the accuracy of the accelerometer is 100 ug.
Figure 9. The vehicular experiment scene and the related equipment.
Figure 9. The vehicular experiment scene and the related equipment.
Sensors 15 27930 g009
First, the vehicle remained static and the engine was not working. The outputs of the IMU were sampled for ten minutes. The average of the last two minutes’ data was computed and then used for the AA. The alignment results using fine alignment method given the data of ten minutes were considered as the true attitude. The alignment results are shown in Table 6. The errors of the attitude for the 14 AA methods are shown in Figure 10, where δ γ , δ θ and δ ψ denote the errors of the roll angle, pitch angle and yaw angle, respectively.
From Table 6 and Figure 10, when ( f x , f z , w z ) or ( f y , f z , w z ) are selected to conduct the selective alignment. The error of the yaw angle is worse. That is because they are close to the b 5 = 1 exception as described in Section 4.3. Besides, if f y is selected for the selective alignment, the pitch angles are the same. If the three selections from the six IMU outputs for the selective alignment include two gyroscope outputs or exclude f x or f y , the accuracy of the roll or pitch is worse, as shown in Table 6.
Table 6. The attitude angles obtained using AA methods.
Table 6. The attitude angles obtained using AA methods.
f x w x f y w y f z w z Roll ( ° )Pitch ( ° )Yaw ( ° )
Selective Alignment1−0.07441.6172292.8975
2−0.07441.6537292.8974
3−0.07441.7247292.8973
4−0.07441.6760292.8974
5−0.07441.6172292.9183
6−0.08991.6172292.9183
7−0.60411.6172292.9183
8−0.09971.6172292.9183
9 −0.07441.7247292.1757
10−0.05361.7255292.8692
11−0.60411.6172274.1607
12−0.05391.7255292.8565
Analytic Coarse Alignment13Equation (6a)−0.07441.6172292.9162
14Equation (6b)−0.07441.6172292.9145
True attitude−0.07451.6173292.9206
Figure 10. The errors of the attitude obtained using AA methods.
Figure 10. The errors of the attitude obtained using AA methods.
Sensors 15 27930 g010

5.5. Alignment in the Vehicle with Vibration

Because vibration is common in actual systems, the alignment in the base with vibration is studied. The vehicle remained static for ten minutes, then fire the engine and let the engine idle. The outputs of the three-axis gyroscope were greatly affected because of vibration caused by the idle engine. And the yaw angle obtained through AA could be greatly influenced. The outputs from the IMU for two minutes during the time when the engine was idling were used for the application of the AA. And the alignment results using fine alignment method given the outputs of the whole time were considered as the true attitude.
Since the six outputs from the IMU have a redundancy for the AA. The accuracy of the alignment could be increased compared to other AA methods by using two or more sets of selections for the selective alignment and grouping the alignment results in a certain way. The roll and pitch angles can be determined due to the outputs of the three-axis accelerometer, while outputs from the three-axis gyroscope are required to obtain the yaw angle. The outputs of the three-axis gyroscope are much more easily influenced than the three-axis accelerometer in the base with vibration. The following AA method based on the selective alignment is presented. The yaw angle obtained through the selective alignment with the selection of ( f x , f y , w x ) is denoted by ψ x = h x ( w x ) . The yaw angle obtained through the selective alignment with the selection of ( f x , f y , w y ) is denoted by ψ y = h y ( w y ) . Approximatively, w x ~ N ( u x , σ x 2 ) , w y ~ N ( u y , σ y 2 ) , where N ( u t , σ t 2 ) denotes the normal distribution with the mean u t and the standard deviation σ t , t = x , y . u t was estimated as the average value of the data for two minutes from the output in t -axis gyroscope, and σ t 2 was estimated as the computed variance value from the same data divided by the number of the samples. The variance of the yaw angle could be decreased using
ψ x y = ( 1 k ) ψ x + k ψ y
where k = D ( ψ x ) D ( ψ x ) + D ( ψ y ) , and it results in D ( ψ x y ) = D ( ψ x ) D ( ψ y ) D ( ψ x ) + D ( ψ y ) . D ( ) is the variance operator.
D ( ψ i ) was computed using a simple method by constructing three-element array. The principle of the method is described as follows. If w ~ N ( u , σ 2 ) , then
E ( w ) = u , D ( w ) = σ 2 , E ( w 2 ) = u 2 + σ 2 , E ( w 3 ) = u 3 + 3 u σ 2
where E ( ) is the expect operator. Construct w = [ u 3 2 σ , u , u + 3 2 σ ] T , which satisfies
E ( w ) = u , D ( w ) = σ 2 , E ( w 2 ) = u 2 + σ 2 , E ( w 3 ) = u 3 + 3 u σ 2
It can be observed that the estimation accuracy using the three-element array can reach the third order.
Similarly, in order to compute the variance of the yaw angle due to the influenced gyroscope output, construct w t = [ u t 3 2 σ t , u t , u t + 3 2 σ t ] T as the outputs from one-axis gyroscope, the variance of the yaw angle using the selective alignment method can be estimated by
D ( ψ t ) = 1 3 i = 1 3 { h t [ w t ( i ) ] ψ t ¯ }
where ψ t ¯ = 1 3 i = 1 3 h t [ w t ( i ) ] , t = x , y .
The experiment was conducted for fifty times. Assume the yaw angle using the i -th method shown in Figure 11 for the j -th experiment is ψ ^ i j , and the corresponding true value using fine alignment method is ψ j . | δ i | and σ i are used to test the average alignment error and the stability, respectively, for the i -th method shown in Figure 11. They are given by
| δ i | = 1 50 j = 1 50 | ψ ˜ i j | , σ i = 1 50 1 j = 1 50 ( ψ ˜ i j ψ ˜ ¯ i ) 2
where ψ ˜ i j = ψ ^ i j ψ j , ψ ˜ ¯ i = 1 50 j = 1 50 ψ ˜ i j .
The results are shown in Figure 11. It can be observed that the yaw angle using Equation (35) based on the selective alignment achieves the best accuracy and stability among the AA methods. However, it has to be noted that AA could usually not achieve better alignment result in the base with obvious vibration than INCA presented in the introduction section. Nevertheless, AA has its vitality for the self-alignment of the SINS on the ground, especially in the base with limited vibration, where both AA and INCA could achieve consistent alignment precision and AA requires less alignment time and has much lower computational complexity. Moreover, the principle of the selective alignment shows that two vectors are redundant to determine the coordinate transformation matrix, while at least two vectors are required for the TRIAD algorithm. Given the idea, some more improvement may be achieved in some other applications where TRIAD algorithm is used such as the spacecraft attitude determination problem.
Figure 11. The errors of the yaw angle obtained using AA methods.
Figure 11. The errors of the yaw angle obtained using AA methods.
Sensors 15 27930 g011

6. Conclusions

A new AA method called selective alignment is proposed, which uses only three outputs of the IMU and a few properties from the remaining outputs such as the sign and the approximate value to obtain the attitude of the SINS relative to the Earth in the static base on the ground. Simulations and experiments were conducted, and the results demonstrate that the selective alignment can be used for self-alignment, fault detection and self-calibration, and the selective alignment could achieve better alignment result than the analytic coarse alignment which is based on the TRIAD algorithm. Selective alignment algorithm, as an extension of the TRIAD algorithm, implies that two non-coplanar vectors are redundant for determining the transformation matrix, which could be studied further in other attitude determination problems where TRIAD algorithm is applied such as spacecraft attitude estimation.

Acknowledgments

This study was supported by the Funding Program for Graduate Student Research and Innovation in Jiangsu Province, China (Grant No. CXLX11_0239). The authors would like to thank Guo Wei of the National University of Defense Technology for his invaluable assistance with the experiments.

Author Contributions

C. Tan conceived and designed the study. Y. Wang, X. Zhu and Z. Wu gave suggestions for improving the system scheme. Y. Su offered the need for the experiment. C. Tan wrote the paper. D. Gu reviewed the manuscript. Y. Su and X. Zhu oversaw the project. The authors approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Britting, K.R. Inertial Navigation Systems Analysis; John Wiley & Sons: New York, NY, USA, 1971. [Google Scholar]
  2. Jiang, Y.F. Error analysis of analytic coarse alignment methods. Aerosp. Electron. Syst. Soc. 1998, 34, 334–337. [Google Scholar] [CrossRef]
  3. Zhao, H.; Shang, H.; Wang, Z.; Jiang, M. Comparison of initial alignment methods for sins. In Proceedings of the 9th World Congress on Intelligent Control and Automation (WCICA, 2011), Taipei, Taiwan, 21–25 June 2011; pp. 42–47.
  4. Li, J.; Fang, J.; Du, M. Error analysis and gyro-bias calibration of analytic coarse alignment for airborne POS. IEEE Instrum. Meas. Soc. 2012, 61, 3058–3064. [Google Scholar] [CrossRef]
  5. Harold, D.K.; Harold, D.K. A passive system for determining the attitude of a satellite. Aiaa J. 1963, 2, 1350–1351. [Google Scholar]
  6. Shuster, M.D.; OH, S.D. Three-axis attitude determination from vector observations. J. Guid. Control Dyn. 1981, 4, 70–77. [Google Scholar] [CrossRef]
  7. Silson, P.M. Coarse alignment of a ship’s strapdown inertial attitude reference system using velocity loci. IEEE Instrum. Meas. Soc. 2011, 60, 1930–1941. [Google Scholar] [CrossRef]
  8. 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 on Position, Location and Navigation Symposium, Monterey, CA, USA, 5–8 May 2008; pp. 961–965.
  9. Ali, J.; Ullah Baig Mirza, M.R. Initial orientation of inertial navigation system realized through nonlinear modeling and filtering. Measurement 2011, 44, 793–801. [Google Scholar] [CrossRef]
  10. Chang, L.; Hu, B.; Li, A.; Qin, F. Strapdown inertial navigation system alignment based on marginalised unscented kalman filter. IET Sci. Meas. Technol. 2013, 7, 128–138. [Google Scholar] [CrossRef]
  11. 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]
  12. Ben, Y.; Li, Q.; Zhang, Y.; Huo, L. Time-varying gyrocompass alignment for fiber-optic-gyro inertial navigation system with large misalignment angle. Opt. Eng. 2014, 53. [Google Scholar] [CrossRef]
  13. Liu, X.; Xu, X.; Wang, L.; Liu, Y. A fast compass alignment method for sins based on saved data and repeated navigation solution. Measurement 2013, 46, 3836–3846. [Google Scholar] [CrossRef]
  14. Shuster, M.D. The optimization of triad. J. Astronaut. Sci. 2007, 55, 245–257. [Google Scholar] [CrossRef]
  15. Tan, C.; Wang, Y.; Zhu, X.; Su, Y.; Wei, G. Improved alignment method for a sins using two vector measurements. In Proceedings of the Fifth International Conference on Instrumentation, Measurement, Computer, Communication and Control (IMCCC), Qinhuangdao, China, 18–20 September 2015. (accepted).
  16. Bar-Itzhack, I.Y.; Harman, R.R. Optimized triad algorithm for attitude determination. J. Guid. Control Dyn. 1997, 20, 208–211. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Tan, C.; Zhu, X.; Su, Y.; Wang, Y.; Wu, Z.; Gu, D. A New Analytic Alignment Method for a SINS. Sensors 2015, 15, 27930-27953. https://doi.org/10.3390/s151127930

AMA Style

Tan C, Zhu X, Su Y, Wang Y, Wu Z, Gu D. A New Analytic Alignment Method for a SINS. Sensors. 2015; 15(11):27930-27953. https://doi.org/10.3390/s151127930

Chicago/Turabian Style

Tan, Caiming, Xinhua Zhu, Yan Su, Yu Wang, Zhiqiang Wu, and Dongbing Gu. 2015. "A New Analytic Alignment Method for a SINS" Sensors 15, no. 11: 27930-27953. https://doi.org/10.3390/s151127930

APA Style

Tan, C., Zhu, X., Su, Y., Wang, Y., Wu, Z., & Gu, D. (2015). A New Analytic Alignment Method for a SINS. Sensors, 15(11), 27930-27953. https://doi.org/10.3390/s151127930

Article Metrics

Back to TopTop