Next Article in Journal
A Hybrid Scheme for Fine-Grained Search and Access Authorization in Fog Computing Environment
Previous Article in Journal
Surface Acoustic Wave Sensors for Hydrogen and Deuterium Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions

College of Automation, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1416; https://doi.org/10.3390/s17061416
Submission received: 11 April 2017 / Revised: 11 June 2017 / Accepted: 14 June 2017 / Published: 16 June 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
An accurate initial alignment must be required for inertial navigation system (INS). The performance of initial alignment directly affects the following navigation accuracy. However, the rapid convergence of meridians and the small horizontal component of rotation of Earth make the traditional alignment methods ineffective in polar regions. In this paper, from the perspective of global inertial navigation, a novel alignment algorithm based on pseudo-Earth frame and backward process is proposed to implement the initial alignment in polar regions. Considering that an accurate coarse alignment of azimuth is difficult to obtain in polar regions, the dynamic error modeling with large azimuth misalignment angle is designed. At the end of alignment phase, the strapdown attitude matrix relative to local geographic frame is obtained without influence of position errors and cumbersome computation. As a result, it would be more convenient to access the following polar navigation system. Then, it is also expected to unify the polar alignment algorithm as much as possible, thereby further unifying the form of external reference information. Finally, semi-physical static simulation and in-motion tests with large azimuth misalignment angle assisted by unscented Kalman filter (UKF) validate the effectiveness of the proposed method.

1. Introduction

With the development of worldwide transportation, the importance of polar navigation is increasing. The polar navigation technology can ensure the safety and reliability of vehicles when they sail in polar regions. However, due to the special geographical and electromagnetic characteristics [1], common navigation methods such as satellite navigation, radio navigation and geomagnetic navigation do not always work efficiently in polar regions [2,3]. As a result, inertial navigation system (INS), which is a highly autonomous and self-contained navigation system, becomes an optimal choice for polar navigation [4,5].
However, in polar regions, the rapid convergence of meridians makes the traditional north-oriented mechanization ineffective because the characteristic that the y-axis always points toward the north would result in both mechanical and mathematical difficulties in a platform instrument or strapdown attitude matrix [6,7]. As a result, the wander navigation system [8], grid navigation system [9], and transverse navigation system [10] have been designed and developed to solve the problem. Nonetheless, any one of these navigation systems does not have the ability of global navigation. Consequently, a combination of them is usually applied to achieve the global navigation [11]. To benefit from the understandability and familiarity, the local geographic frame is being now widely used as the navigation frame in mid-low latitude regions. As a result, the combination of the north-oriented navigation system and a wander navigation system, or a grid navigation system, or a transverse navigation system is currently selected to achieve the global navigation. Then, a need to transform the navigation parameters relative to local geographic frame to other navigation frames has also developed for global navigation. Therefore, if the initial navigation parameters relative to local geographic frame are obtained at the end of polar alignment phase, it would be more convenient to access the following polar navigation phase without any additional workloads and complexities.
Since strapdown INS (SINS) is a dead-reckoning navigation system, an accurate initial alignment is required before navigation [12]. The performance of initial alignment directly affects the following navigation accuracy. For SINS, since initial velocity and position can be easily obtained by external reference navigation systems such as GPS, the purpose of initial alignment is to determine the coordinate transformation matrix from the body frame to computational navigation frame [13,14]. In the last few decades, an amount of work has been done for initial alignment of SINS. However, most of the previous research works are focused on the initial alignment of mid-low latitude. The polar alignment is still a severe challenge for global navigation of SINS.
For initial alignment operated in the polar regions, two major problems need to be addressed. Firstly, the navigation errors with the traditional north-oriented frame are enlarged sharply [8]. In traditional north-oriented mechanization, the propagation of azimuth error is proportional to both linear movement error and tangent function of latitude. As a result, the attitude reference relative to local geographic frame would lose its effectiveness for polar alignment. Therefore, a suitable attitude reference frame needs to be selected for polar alignment. Secondly, the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude [15]. The degree of observability of azimuth would become smaller [16]. As a result, the initial alignment with a horizontal frame would be not achieved by self-alignment in and near the North/South Pole. Moreover, it is also difficult to obtain an accurate coarse alignment of azimuth in polar regions. Therefore, a kind of valid external observation information of filter needs to be selected to achieve polar alignment. In addition, the alignment modeling with large azimuth misalignment angle should also be designed.
For selection of attitude reference frame in polar regions, the wander frame, the grid frame, and the transverse frame have been proposed for polar navigation. However, the wander attitude reference frame just can be applied in most of polar regions because the definition of wander angle is invalid at 90-degree latitude [17]. In the grid frame [18], the grid reference line is set for determining the azimuth angle. The grid attitude reference frame is introduced for polar alignment. In the transverse frame [19], a transverse attitude reference frame is defined there. According to above description, either the grid frame or the transverse frame can be selected as attitude reference frame for polar alignment [18,19,20,21]. However, just the attitude matrix relative to the grid frame or transverse frame is usually acquired when the initial alignment with these two attitude reference frames is completed. As a result, some additional transformation may need to be added to access the following polar navigation phase. Because the wander navigation system, or grid navigation system, or transverse navigation system may be selected as the following navigation system. Moreover, the transformation of attitude matrix with different reference frames would be inaccurate due to the position errors. The influence of position errors on the alignment accuracy is proportional to both the latitude error and secant function of latitude [22]. Consequently, near the pole regions, the influence would become very large. Thus, the influence of position errors on transformation accuracy must be considered in polar regions. Furthermore, the transformation of attitude matrix is also difficult due to the cumbersome derivation and computational complexity [17,18,19,20,21,22].
On the other hand, for polar alignment, the external observation information of filter also needs to be selected validly. Due to the small horizontal component of rotation of Earth, some observation information of filter applied usually in mid-low latitude would be invalid in polar regions, especially in and near pole. For initial alignment, the external position or velocity information is usually selected as the observation information of filter [23,24]. However, more time would be required to achieve the initial alignment when the external position information is selected as the observation information. The accelerated maneuver might also need to be conducted for polar alignment due to the smaller degree of observability of azimuth [25]. Hence, the external velocity reference is preferred selection for the observation information of filter. For velocity reference information, there are usually two forms: velocity information projected in body frame ( V b ) and velocity information projected in the navigation frame ( V n ). However, polar alignment assisted by V b would also be difficult to achieve due to the smaller horizontal component of rotation of Earth. As a result, V n should be selected as the observation information of filter to achieve the polar alignment [19,21]. This means that polar alignment with the grid frame or transverse frame needs different forms of velocity reference information due to the fact that the different attitude reference frames are employed in the grid alignment or transverse alignment. As a result, the grid alignment or transverse alignment would be disadvantageous to achieve a unified polar alignment algorithm and to unify the form of external reference information further.
According to above analysis and considering the features of global navigation of SINS, a polar alignment algorithm, which can directly obtain the attitude matrix relative to local geographic frame when the initial alignment is finished, is expected to be proposed. The transformation of attitude matrix is not influenced by position errors. Moreover, a unified polar alignment algorithm is also expected to be explored as much as possible. Then the form of external observation information of filter can also be unified for polar alignment in future. Furthermore, it would also be advantage to implement the acquirement of the observation information in a form from the external navigation device. On the other hand, since the definition of coordinate system just like that in the traditional Earth frame, the transverse frame is more immediate and convenient [19]. Therefore, an improved SINS mechanization based on pseudo-Earth frame [26], which is a generalized transverse Earth frame, is proposed in this paper. Then, a novel alignment algorithm based on pseudo-Earth frame and backward process is proposed to achieve the initial alignment in polar regions.
The purpose of this research is expected to provide a new idea to solve the problem of polar alignment and explores a unified polar alignment algorithm as much as possible. The rest of this paper is organized as follows: the proposed polar alignment scheme based on pseudo-Earth frame and backward process is presented in Section 2. In Section 3, the dynamic error modeling with large azimuth misalignment angle is designed. The nonlinear model of polar alignment is given in Section 4. The semi-physical static simulation and in-motion tests with the proposed method are presented in Section 5, and Section 6 is the conclusion.

2. The Polar Alignment Scheme Based on Pseudo-Earth Frame

2.1. The Mechanization of SINS with Pseudo-Earth Frame

The pseudo-Earth frame ( p -frame) can be obtained by revolving Earth-Centered Earth-Fixed frame ( e -frame). As shown in Figure 1, the origin of the p -frame ( o p x p y p z p ) is still at the Earth’s mass center; N p denotes the pseudo-North Pole; the o p x p axis overlaps with previous o e z e axis; the o p y p axis is pointing from the center of Earth to the projection point of initial position of vehicle S on the equatorial plane; and the o p z p axis is perpendicular to the o p y p axis and conforms to the right hand rule, while laying in the equatorial plane. The detailed definition of the p -frame can be found in [26]. Supposing the initial longitude, latitude and altitude of vehicle are λ 0 , L 0 and h 0 , respectively. Then, the rotating matrix between the p -frame and e -frame can be written as:
C e p = C x ( λ 0 90 ) C y ( 90 )
where C x ( λ 0 90 ) = [ 1 0 0 0 cos ( λ 0 90 ) sin ( λ 0 90 ) 0 sin ( λ 0 90 ) cos ( λ 0 90 ) ] ; C y ( 90 ) = [ cos ( 90 ) 0 sin ( 90 ) 0 1 0 sin ( 90 ) 0 cos ( 90 ) ] .
For the definitions of pseudo-Earth coordinate system, they are similar to the ones in traditional Earth coordinate system [26]. The position of vehicle is determined by pseudo longitude λ p , pseudo latitude L p , and pseudo altitude h p . The pseudo geographic frame o g p x g p y g p z g p (namely, E p N p U p ) is also similar to the East–North–Up geographic frame (namely, g -frame), which is shown in Figure 1. The initial pseudo position of vehicle can also be given by:
[ L p 0 λ p 0 h p 0 ] T = [ 0 90 L 0 h 0 ] T
According to above definitions, the pseudo-Earth frame is a generalized transverse Earth frame built with the initial position of vehicle in fact [17]. Therefore, the pseudo navigation mechanization of SINS is also similar to the one with transversal Earth coordinate system. In pseudo navigation mechanization, the g p -frame is selected as navigation frame ( n -frame). Only the angular velocity of p -frame relative to inertial frame ( i -frame) differs from that in the local-level coordinate system. As a result, the pseudo navigation equations can be written by [17,26]:
C ˙ b n = C b n ( ω n b b × )
V ˙ p n = C b n f b ( 2 ω i p n + ω p n n ) × V p n + g n
L ˙ p = V N p R + h p
λ ˙ p = V E p sec L p R + h p
h ˙ p = V U p
where ( × ) denotes the skew symmetric matrix of ( ); ω a b c denotes the rotation angular rate of frame b relative to frame a represented in the c frame; ω i p n = [ Ω sin λ p Ω sin L p cos λ p Ω cos L p cos λ p ] T ; Ω denotes the angular rate of Earth rotation; ω p n n = [ V N p / ( R + h p ) V E p / ( R + h p ) ω p n n = V E p tan L p / ( R + h p ) ] T ; R denotes the Earth radius; ω n b b = ω i b b C n b ω i n n , ω i n n = ω i p n + ω p n n , C n b = ( C b n ) T ; ω i b b denotes the outputs of gyro; f b denotes the specific force; and g n denotes the gravity vector in the n-frame, g n [ 0 0 g ] T .
Since the pseudo-Earth frame is a generalized transverse Earth frame and the pseudo latitude at the initial position (namely, L p 0 ) is always zero from Equation (2), the pseudo-Earth frame can solve the problem of attitude reference frame of polar alignment. In addition, the polar alignment algorithm based on pseudo-Earth frame is also easy to implement because it is similar to the polar alignment algorithm with the transverse Earth frame. However, if the strapdown attitude matrix relative to g -frame is obtained by direct transformation with both the pseudo attitude matrix and the relationship of them, the attitude matrix would be still influenced by the position errors. The cumbersome derivation and computational complexity must also be confronted. Nevertheless, by computation, it is easily find that the initial transformation matrix between g p -frame and g -frame at initial position is a constant matrix and never change with the initial position of vehicle. As a result, it would be very advantageous to achieve the transformation of attitude matrix between g -frame and g p -frame without the influence of position errors and the cumbersome computation. The initial transformation matrix C g g p ( 0 ) can also be given by:
C g g p ( 0 ) = C p g p ( 0 ) C e p C g e ( 0 )
where C p g p ( 0 ) = [ sin λ p 0 cos λ p 0 0 sin L p 0 cos λ p 0 sin L p 0 sin λ p 0 cos L p 0 cos L p 0 cos λ p 0 cos L p 0 sin λ p 0 sin L p 0 ] ; C g e ( 0 ) = [ sin λ 0 sin L 0 cos λ 0 cos L 0 cos λ 0 cos λ 0 sin L 0 sin λ 0 cos L 0 sin λ 0 0 cos L 0 sin L 0 ] .
Substituting Equations (1) and (2) into Equation (6), the initial transformation matrix C g g p ( 0 ) can be rewritten by:
C g g p ( 0 ) = [ 0 1 0 1 0 0 0 0 1 ]
From Equation (7), we can obtain that the g p -frame is identical to the South–East–Up geographic frame at initial position. As a result, if the pseudo attitude matrix of vehicle at initial position is obtained when the initial alignment is completed, it would be equal to the attitude matrix relative to local geographic frame is obtained without influence of transformation errors and cumbersome computation. Moreover, the coarse alignment method needs also not to be redesigned, and the traditional coarse alignment algorithm can be applied directly. Therefore, the pseudo navigation mechanization would be simpler to solve the problem of attitude reference frame for polar alignment. Nevertheless, the pseudo attitude matrix of vehicle at initial position is expected to be obtained when the initial alignment based on pseudo navigation mechanization is completed.
Remark. 
If the initial latitude of vehicle relative to e -frame is 90°, the initial longitude of vehicle can choose the arbitrary value with [ 180 , 180 ] . Then, the pseudo-Earth frame is built based on the selected longitude. The initial azimuth angle can still be obtained by the initial pseudo azimuth angle from Equation (7). Nevertheless, we should notice that the azimuth angle relative to g -frame in pole is the angle between the vehicle direction and the selected longitude.

2.2. Backward Process for Pseudo Navigation System

With electronic technology development, the capabilities of the navigation computer become large. The backward process of data has been proposed for initial alignment [27,28,29]. By utilizing the saved data of inertial measurement unit (IMU) with the opposite time sequence, the backward process can accelerate the process of alignment. Nevertheless, here the requirement that the pseudo attitude matrix of vehicle at initial position is obtained at the end of alignment phase is expected to be realized by utilizing the saved IMU data. If the pseudo starpdown attitude matrix and the pseudo position in backward process are identical to the forward ones at the same time, the pseudo position with backward process can return to initial position. Then, the initial attitude matrix relative to g -frame can be obtained from Equation (7). On the other hand, the backward process can still be applied to decrease the alignment time.
From Equations (3)–(5), the discrete form of forward pseudo navigation equations can be given as follows:
C b k n = C b k 1 n ( I + T ω n b k 1 b × )
V p k n = V p k 1 n + T C b k 1 n f k 1 b T ( 2 ω i p k 1 n + ω p n k 1 n ) × V p k 1 n + T g k 1 n
L p k = L p k 1 + T V N p k 1 R + h p k 1
λ p k = λ p k 1 + T V E p k 1 sec L p k 1 R + h p k 1
h p k = h p k 1 + T V U p k 1
where k is the discrete time index for forward computation, and k = 0 , 1 , 2 , m ; A k denotes the computed value of A at time t k , while A k 1 denotes the computed value of A at time t k 1 ; T denotes the sampling period; I denotes the identity matrix.
Then, through further transposition, the backward pseudo navigation equations can also be obtained:
C b k 1 n = C b k n ( I + T ω n b k 1 b × ) 1 C b k n ( I T ω n b k 1 b × )
V p k 1 n = V p k n T C b k 1 n f k 1 b + T ( 2 ω i p k 1 n + ω p n k 1 n ) × V p k 1 n T g k 1 n
L p k 1 = L p k T V N p k 1 R + h p k 1
λ p k 1 = λ p k T V E p k 1 sec L p k 1 R + h p k 1
h p k 1 = h p k T V U p k 1
Moreover, assuming the discrete time index for backward computation is l ( l = 0 , 1 , 2 , m ), the backward pseudo navigation equations can be written as follows:
C b l n = C b l 1 n ( I + T ω n b l 1 b × )
V p l n = V p l 1 n + T C b l 1 n f l 1 b T ( 2 ω i p l 1 n + ω p n l 1 n ) × V p l 1 n + T g l 1 n
L p l = L p l 1 + T V N p l 1 R + h p l 1
λ p l = λ p l 1 + T V E p l 1 sec L p l 1 R + h p l 1
h p l = h p l 1 + T V U p l 1
where the superscript “ ” denotes this term is for the backward process.
Defining ξ denotes the ξ th sample with the normal time sequence ( ξ = 1 , 2 , m ) and η denotes the η th sample with the opposite time sequence ( η = 1 , 2 , m ). According to the forward computed process and backward computed process above, the corresponding relations of time index can be shown in Figure 2.
According to the analysis above, the backward strapdown attitude matrix should equal to the forward one at the same time. As a result, the corresponding terms of Equations (11) and (14) should also be equal. Comparing Equation (11) and Equation (14), we can obtain:
{ C b l n = C b k 1 n C b l 1 n = C b k n
T ω n b l 1 b = T ω n b k 1 b
Considering the corresponding relations of time index from Figure 2, we further have:
l + k 1 = m
ω n b b ( η ) = ω n b b ( ξ )
Obviously, ξ = k , η = l . Combining Equations (19) and (20) can be rewritten by:
ω n b b ( l ) = ω n b b ( m l + 1 )
On the other hand, it is also expected that the backward pseudo position could be identical to the forward ones at the same time. Thus, comparing the corresponding terms of Equations (13) and (16), we can similarly obtain:
{ l + k 1 = m V E p ( η ) = V E p ( ξ ) V N p ( η ) = V N p ( ξ ) V U p ( η ) = V U p ( ξ )
Then, we can further have:
V p n ( l ) = V p n ( m l + 1 )
From Equation (23), we can conclude that the sign of the pseudo velocity in backward process should be opposite to the corresponded one in forward process. Therefore, after taking opposite, Equation (12) should equal to Equation (15), correspondingly. Then we can get:
{ l + k 1 = m f b ( η ) = f b ( ξ ) g n ( η ) = g n ( ξ ) [ 2 ω i p n ( η ) + ω p n n ( η ) ] × V p n ( η ) = [ 2 ω i p n ( ξ ) + ω p n n ( ξ ) ] × V p n ( ξ )
From Equation (24), we have:
f b ( l ) = f b ( m l + 1 )
g n ( l ) = g n ( m l + 1 )
[ 2 ω i p n ( l ) + ω p n n ( l ) ] × V p n ( l ) = [ 2 ω i p n ( m l + 1 ) + ω p n n ( m l + 1 ) ] × V p n ( m l + 1 )
From Equations (25) and (26), the signs of the outputs of accelerometer and the gravity in backward process should be the same as the forward ones. Moreover, substituting Equation (23) into Equation (27), we can obtain:
ω i p n ( l ) = ω i p n ( m l + 1 )
ω p n n ( l ) = ω p n n ( m l + 1 )
From Equation (28), we can conclude that the sign of the angular rate of rotation of Earth in backward process needs to be taken negative. In addition, combining with Equations (21), (28), and (29), the outputs of gyro in backward process can be obtained by:
ω i b b ( l ) = ω i b b ( m l + 1 )
From Equation (30), the sign of outputs of gyro should be opposite to the forward one in the backward process.
According to above derivation, the expectation that the pseudo starpdown attitude matrix and the pseudo position in backward process should be the same as the forward ones at the same time can be realized. While the pseudo velocity of backward process is opposite to the forward one at the same time. Moreover, the signs of the outputs of gyro and the angular rate of rotation of Earth in backward process should be opposite to the corresponded ones with the forward process. The others are all the same as the forward ones. As a result, the expectation that the pseudo attitude matrix of vehicle at initial position is obtained at the end of alignment phase can be realized by backward process with pseudo navigation mechanization. Further, the initial transformation matrix C g g p ( 0 ) can be utilized to obtain the strapdown attitude matrix relative to g -frame with a simple transformation logic and without influence of position errors.

2.3. The Proposed Polar Alignment Scheme with Pseudo-Earth frame

From the analysis above, the alignment algorithm based on pseudo-Earth frame and backward process can be applied to solve the problem of polar alignment. The attitude matrix relative to local geographic frame is obtained when the initial alignment is completed. Moreover, it is also easy to implement without the influence of position errors and the cumbersome computation since the initial transformation matrix is a fixed known constant matrix and never influenced by initial position of vehicle from Equation (7).
Supposing the IMU data are saved in first forward fine alignment and the sampling index of saved data sequence is 1 to m. Moreover, the computed pseudo attitude matrix C b g p ( m ) 1 is also obtained at the end of first forward process. Subsequently, the backward process of the filter estimation is carried out from m to 1, and the pseudo attitude matrix C b g p ( 0 ) 2 is obtained at the end of first backward process. Then, the forward process and backward process based on pseudo-Earth frame is repeatedly conducted until the initial alignment is completed. Comparing with the traditional application of the backward process, one and only difference is that the proposed alignment algorithm should end with the backward process the filter estimation from m to 1. Meanwhile, the computed pseudo attitude matrix of time 0 (namely, C b g p ( 0 ) n ) is obtained. As a result, the pseudo attitude matrix at initial position is obtained by the backward process. The attitude matrix relative to g -frame is also obtained from Equation (7). Then, the navigation parameters relative to g -frame can be transformed to the corresponding reference frame of the following polar navigation system by the existing transformation program without any additional workloads and complexities. In addition, since the reciprocating process spends some time, the IMU data also need to be saved. The sampling index of saved data in the delay duration is represented from m to s. According to the transformational initial navigation parameters, the navigation calculation is conducted from the first sample to the sth sample by the selected polar navigation system, such as the wander navigation system, grid navigation system and transverse navigation system and so on. Then, the SINS accesses the navigation phase. The proposed polar alignment scheme based on pseudo-Earth frame and backward process is shown in Figure 3.
In addition, the implementation of the proposed alignment algorithm is also illustrated in Figure 4. As shown in Figure 4, the traditional coarse alignment algorithm can be still applied to the proposed polar alignment process. Only an initialization with simple transformation is conducted to access the next fine alignment associated with pseudo-Earth frame. Moreover, since the initial pseudo latitude is always zero, the alignment model with the pseudo-Earth frame would be always equivalent to the traditional alignment model of equatorial regions. As a result, the initial alignment with the pseudo-Earth frame can solve the attitude reference problem of polar alignment. Furthermore, since the initial transformation matrix C g g p ( 0 ) is a fixed known constant matrix, the backward process of SINS is proposed here, thereby obtaining the pseudo attitude matrix of vehicle at initial position at the end of initial alignment. As a result, the initial transformation matrix can be cleverly applied to obtain the strapdown attitude matrix relative to g -frame with a simple transformation and without influence of position errors. Thus, there is less cumbersome derivation and low computational complexity in the process of transformation, thereby decreasing the complexity of program and the burden of computer. Moreover, the backward process can still accelerate the filtering estimation and decrease the alignment time. On the other hand, the strapdown attitude matrix relative to g -frame is obtained at the end of alignment phase. It would also be more convenient to access any one of the following polar navigation systems. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, it would make it possible to achieve the unification of polar alignment algorithm. Therefore, the proposed polar alignment algorithm has superior performance.

3. Pseudo Navigation Error Equations with Large Azimuth Misalignment Angle

Considering that an accurate coarse alignment of azimuth is difficult to obtain in polar regions, the dynamic error modeling with large azimuth misalignment angle is designed in this section. Defining the pseudo misalignment angle between computed navigational frame ( n -frame) and n -frame is ϕ p = [ ϕ x p ϕ y p ϕ z p ] T . Here, the pseudo-pitch and pseudo-roll misalignment angles (namely, ϕ x p and ϕ y p ) are small misalignment angles. While the pseudo-azimuth angle should be considered as a large misalignment angle. Then, the rotating matrix between n -frame and n -frame can be written as [30]:
C n n = C ϕ y p C ϕ x p C ϕ z p = [ c ϕ y p c ϕ z p s ϕ y p s ϕ x p s ϕ z p c ϕ y p s ϕ z p + s ϕ y p s ϕ x p c ϕ z p s ϕ y p c ϕ x p c ϕ x p s ϕ z p c ϕ x p c ϕ z p s ϕ x p s ϕ y p c ϕ z p + c ϕ y p s ϕ x p s ϕ z p s ϕ y p s ϕ z p c ϕ y p s ϕ x p c ϕ z p c ϕ y p c ϕ x p ]
where s ϕ j p and c ϕ j p represent the sine function and cosine function, respectively, and j = x , y , z .
Defining the angular velocity of n -frame relative to n -frame is ω n n n , and the relationship between ω n n n and ϕ p can be given by [30]:
ω n n n = C ϕ y p C ϕ x p [ 0 0 ϕ ˙ z p ] + C ϕ y p [ ϕ ˙ x p 0 0 ] + [ 0 ϕ ˙ z p 0 ] = A [ ϕ ˙ x p ϕ ˙ y p ϕ ˙ z p ]
From Equation (32), we have:
ϕ ˙ p = A 1 ω n n n
where
A 1 = 1 c ϕ x p [ c ϕ y p c ϕ x p 0 s ϕ y p c ϕ x p s ϕ y p s ϕ x p c ϕ x p c ϕ y p s ϕ x p s ϕ y p 0 c ϕ y p ]
Only considering that the pseudo-azimuth angle is a large misalignment angle here, Equations (31) and (34) can be rewritten by:
C n n = [ c ϕ z p s ϕ z p ϕ y p s ϕ z p c ϕ z p ϕ x p ϕ y p c ϕ z p + ϕ x p s ϕ z p ϕ y p s ϕ z p ϕ x p c ϕ z p 1 ]
A 1 = [ 1 0 ϕ y p 0 1 ϕ x p ϕ y p 0 1 ]
Due to the existence of misalignment angles and gyro drifts, the differential equation of actual computed pseudo strapdown attitude matrix can be written as:
C ˙ b n = C b n ( ω ˜ n b b × )
where
ω ˜ n b b = ( ω i b b + ε b ) C n b ( ω i n n + δ ω i n n )
where ε b denotes the gyro drifts, δ ω i n n denotes the calculated error of ω i n n .
Since
C b n = C n n C b n
differentiating both sides of Equation (39) gives:
C ˙ b n = C ˙ n n C b n + C n n C ˙ b n
Substituting Equations (3), (37), and C ˙ n n = C n n ( ω n n n × ) into Equation (40), we can obtain:
C n n ( ω n n n × ) C b n = C b n ( ω n b b × ) C n n C b n ( ω ˜ n b b × )
Multiplying the left by C n n and the right by C n b , we further have:
( ω n n n × ) = C b n ( ω n b b × ) C n b C b n ( ω ˜ n b b × ) C n b
Then, Equation (42) can be rewritten by:
ω n n n = C b n ω n b b C b n ω ˜ n b b
Substituting Equation (38) and ω n b b = ω i b b C n b ω i n n into Equation (43), we can obtain:
ω n n n = ( I C n n ) ω ˜ i n n + C n n δ ω i n n C b n ε b
From Equation (33), the pseudo attitude error equation with large azimuth misalignment angle can be obtained by:
ϕ ˙ p = A 1 ω n n n = A 1 [ ( I C n n ) ω ˜ i n n + C n n δ ω i n n C b n ε b ]
where
δ ω i n n = δ ω i p n + δ ω p n n
{ δ ω i p n = [ ω 11 ω 12 ω 13 ] T ω 11 = Ω cos λ p δ λ p ω 12 = Ω cos L p cos λ p δ L p + Ω sin L p sin λ p δ λ p ω 13 Ω sin L p cos λ   p δ L p Ω cos L p sin λ p δ λ p
{ δ ω p n n = [ ω 21 ω 22 ω 23 ] T ω 21 = V N p ( R + h p ) 2 δ h p δ V N p R + h p ω 22 = V E p ( R + h p ) 2 δ h p + δ V E p R + h p ω 23 = V E p R + h p sec 2 L p δ L p V E p tan L p ( R + h p ) 2 δ h p + δ V E p R + h p tan L p
Since the measurement errors and calculation errors exist, the differential equation of actual computed pseudo velocity can be written as:
V ˜ ˙ p n = C b n f ˜ b ( 2 ω ˜ i p n + ω ˜ p n n ) × V ˜ p n + g ˜ n
where V ˜ p n = V p n + δ V p n , f ˜ b = f b + b , ω ˜ i p n = ω i p n + δ ω i p n , ω ˜ p n n = ω p n n + δ ω p n n , g ˜ n = g n + δ g n , and b denotes the zero-bias of accelerometer.
Then the pseudo velocity error equation can be obtained as:
δ V ˙ p n = V ˜ ˙ p n V ˙ p n
Ignoring δ g n and second order small quantities, substituting Equations (4) and (49) into Equation (50), can obtain:
δ V ˙ p n = ( I C n n ) C b n f ˜ b ( 2 δ ω i p n + δ ω p n n ) × V ˜ p n ( 2 ω ˜ i p n + ω ˜ p n n ) × δ V p n + C n n C b n b
Defining L ˜ p = L p + δ L p , λ ˜ p = λ p + δ λ p , h ˜ p = h p + δ h p , we can obtain the pseudo position error equations from Equation (5). They are can be expressed as:
{ δ L ˙ p = L ˜ ˙ p L ˙ p = V ˜ N p R + h ˜ p δ h p + δ V N p R + h ˜ p δ λ ˙ p = λ ˜ ˙ p λ p = V ˜ E p tan L ˜ p sec L ˜ p R + h ˜ p δ L p V ˜ E p sec L ˜ p ( R + h ˜ p ) 2 δ h p + δ V E p R + h ˜ p sec L ˜ p δ h ˙ p = h ˜ ˙ p h ˙ p = δ V U p

4. The Design of Nonlinear Filter Model for Polar Alignment

Aiming at the nonlinear error models derived above, a suitable nonlinear filtering should be selected to replace the standard Kalman filter (KF). Because the extended KF (EKF) has to bear a number of drawbacks such as cumbersome derivation, evaluation of Jacobian matrices, and larger linearized errors and so on [31,32], as a promising substitute for EKF [33], unscented KF (UKF) is employed to accomplish the polar alignment here.
Assuming the measurement errors of accelerometer and gyro are mainly composed of the constant error and zero mean Gaussian white noise. Then we have ε ˙ b = 0 and ˙ b = 0 . Consequently, the filtering state with pseudo frame can be set as:
X ( t ) = [ ( δ P p ) T ( ϕ p ) T ( δ V p n ) T ( b ) T ( ε b ) T ] T
where δ P p denotes the pseudo position errors, and δ P p = [ δ L p δ λ p δ h p ] T .
According to the nonlinear pseudo error equations in Section 3, the state equation of the filtering model of initial alignment can be obtained as:
{ δ P ˙ p = [ δ L ˙ p δ λ ˙ p δ h ˙ p ] T ϕ ˙ p = A 1 [ ( I C n n ) ω ˜ i n n + C n n δ ω i n n C b n ε b ] δ V ˙ p n = [ I ( C n n ) T ] C b n f ˜ b ( 2 δ ω i p n + δ ω p n n ) × V ˜ p n ( 2 ω ˜ i p n + ω ˜ p n n ) × δ V p n + C n n C b n b ˙ b = 0 3 × 1 ε ˙ b = 0 3 × 1
where 0 3 × 1 denotes the zero matrix of indicated dimension.
Further, the general nonlinear system dynamic model can be written by:
X ˙ ( t ) = f ( X , t ) + W ( t )
where the specific expressions of f ( X , t ) can refer to Equation (54), and W ( t ) is the white system process noise with zero mean and covariance Q ( t ) .
In addition, the differences of the pseudo velocity projected in g p -frame between the SINS and external navigation device should be selected as the observations of filter here since the velocity information projected in body frame ( V b ) would be invalid caused by the smaller horizontal component of rotation of Earth. Then, the measurement equation can also be given by:
Z ( t ) = H X ( t ) + V ( t )
where H = [ 0 3 × 6 I 3 × 3 0 3 × 6 ] is the observation matrix, and V ( t ) is the white noise with zero mean and covariance R ( t ) .
On the other hand, the block diagram of polar alignment assisted by UKF is also shown in Figure 5. For application of UKF to initial alignment with large azimuth misalignment angle, the specific derivation and implement can be found in [30,31,33].

5. Simulation and Test

In order to confirm the performance of the proposed polar alignment scheme, both semi-physical static simulation and in-motion tests, which are assisted by UKF under the large azimuth misalignment angle, are conducted in this section.

5.1. Semi-Physical Static Simulation

Due to the restriction of geography, the semi-physical static simulation is conducted. As attributed to the rotation of Earth, the real static angular velocity measured by gyro at different geographical positions can be exactly known. The angular velocity measured by gyro can be expressed as:
ω ˜ i b b = C g b ω i e g + ε b
where ω ˜ i b b denotes the measured angular velocity; C g b = ( C b g ) T and C b g denotes the attitude matrix of IMU relative to g -frame; ω i e g = [ 0 Ω cos L Ω sin L ] T , L is the latitude in e -frame.
Therefore, the static outputs of gyro in polar regions can be simulated by the measured angular velocity in non-polar regions from Equation (57). The static experimental data are collected with fiber optical gyroscope (FOG) SINS, which fixed in a tri-axis high-precision turntable. The attitude reference of IMU is provided by turntable. The gyro drift and the accelerometer bias of test IMU are less 0.01°/h, 100 μg respectively. The local longitude and latitude is 126.6778° and 45.7778°. The pitch, roll and azimuth of IMU relative to the g -frame are 0°, 0° and 30° respectively. The static outputs of IMU for 1200 s are collected, and the update frequency of IMU is 100 Hz. Moreover, according to Equation (57), the original gyro’s outputs and the modified gyro’s outputs with 85° latitude can be shown in Figure 6.
Then, the static alignment with the pseudo-Earth frame is carried out, which is assisted by UKF under the condition of large azimuth misalignment angle. In this simulation, the initial pseudo misalignment angles ϕ x p , ϕ y p and ϕ z p are chosen as 1°, 1° and 40°, respectively. Since the position of the vehicle never changes in the static alignment, the strapdown attitude matrix relative g -frame can be obtained directly from Equation (7) without the influence of position errors and the cumbersome computation. As a result, the application of backward process to polar alignment has not been presented here. Through the pseudo attitude matrix and the Equation (7), the static alignment results relative g -frame is obtained directly and is shown in Figure 7.
As shown in Figure 7, the estimated errors of level misalignment angles and azimuth misalignment angle are convergent with time. The level alignment errors can converge rapidly. However, the convergent time of azimuth error is about 1000 s. Because the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude, the degree of observability of azimuth would become smaller in polar regions. More time is required for azimuth alignment in polar regions. Moreover, it would also be difficult to achieve self-alignment near the pole. Even so, a better alignment result is still obtained with the proposed pseudo navigation mechanization at 85° latitude. The alignment error of azimuth is −0.3607° in 1200 s. Therefore, the proposed pseudo navigation mechanization can solve the problem of traditional attitude reference frame in polar regions caused by SINS mechanization.

5.2. In-Motion Alignment Tests

Due to the restriction of geography, the actual test data of in-motion alignment are difficult to obtain. Consequently, only a vehicle test of low latitude, which is assisted by UKF under the condition of large azimuth misalignment angle, is firstly conducted to verify the feasibility of the proposed alignment scheme. Then, the polar alignment tests with the simulation data are carried out to validate the superior performance of the proposed method. In addition, in the process of polar alignment, the external velocity observation information projected in the navigation frame (namely, g p -frame) should be employed to fulfill in-motion alignment tests since the horizontal component of rotation of Earth is very small in and near pole regions.
Firstly, the vehicle test of low latitude is executed by equipped with a receiver of differential GPS and two different grades fiber optical gyroscope (FOG) SINSs as shown in Figure 8. The gyro constant drift and the accelerometer constant bias of test IMU are less 0.01°/h, 100 μg respectively. The local longitude and latitude are 126.6778° and 45.7778°, respectively. In total, 250 s vehicle data were collected. The initial pseudo misalignment angles ϕ x p , ϕ y p and ϕ z p are chosen as 1°, 1° and 20°, respectively. In addition, since this paper mainly focuses on the design of polar alignment scheme, the application of backward process to accelerate the alignment process more has not been presented. Only a backward process is conducted to implement polar alignment and accelerate the alignment process here. For the application of backward process to decrease the alignment time further, more details can be found in [27,28,29]. Here, the filtering estimation of initial alignment firstly works with test data by the forward process. Subsequently, the filtering estimation is carried out with the opposite time sequence. The alignment errors and the calculated pseudo position are shown in Figure 9 and Figure 10, respectively.
As shown in Figure 9, the estimated errors of three-axis misalignment angles are convergent with time. The pseudo position can also return to the initial position with the backward process when the initial alignment is completed from Figure 10. Although there are still some differences between the initial value and the final value of pseudo position, and the drift errors of pseudo position are unavoidable because the external velocity information is selected as the observation information of filter in the process of initial alignment, the influence of the drift errors of pseudo position on the initial transformation matrix is very small and can be neglected since the initial pseudo latitude is always zero. Therefore, the pseudo strapdown attitude matrix of vehicle at initial position is obtained at the end of alignment phase. Then, the initial attitude matrix relative to g -frame can be obtained directly with Equation (7). There is less cumbersome derivation and low computational complexity in the process of transformation since the initial transformation matrix is a fixed known constant matrix. Consequently, we can conclude that the proposed polar alignment scheme would be feasible for in-motion alignment and has less complexity.
Next, the superior performance of the proposed method for polar alignment is verified by simulation data due to restriction of geography. In polar in-motion alignment simulations, the outputs of gyro and accelerometer are generated by a strapdown INS simulator. The constant and random drifts of gyro are 0.01 °/h and 0.001 °/h, respectively, and the constant and random biases of the accelerometer are 100 μg and 10 μg, respectively. The initial values of attitude relative to g -frame are all set to be zero, and a constant velocity of 10m/s in the along-vehicle direction is applied all the scenarios. According the definitions of pseudo-Earth frame, the vehicle would sail along the pseudo equator with the constant velocity of 10 m/s. Two groups of in-motion data with the different initial latitudes 87° and 90° are obtained. The initial longitudes are all selected as 126°. The length of test data with the initial latitude 87° is 600 s, while the length of test data with the initial latitude 90° is 1200 s. The update frequency of IMU data is 100 Hz. In the process of initial alignment, the initial pseudo position errors L p , λ p , and h p are all set to be 1 m. The initial pseudo velocity errors are all 0.1 m/s. The initial pseudo misalignment angles ϕ x p , ϕ y p , and ϕ z p are chosen as 1°, 1°, and 70°, respectively. The alignment results of 87° latitude and 90° latitude are shown in Figure 11 and Figure 12, respectively. The calculated pseudo positions with the pseudo navigation mechanization are also shown in Figure 13 and Figure 14. In addition, a comparison with the traditional method based on transverse frame at 87° latitude is also carried out. The final alignment result relative to g -frame based transverse frame is obtained by direct transformation with the calculated transverse position from transverse SINS calculation. The mean and variance of 30 alignment errors in 1200 s are shown in Table 1.
As shown in Figure 11 and Figure 12, at the end of the first forward process, the azimuth alignment errors of 87° and 90° are −13.16° in 600 s and −18.5° in 1200 s, respectively. The alignment errors of azimuth are still larger. While the alignment errors of azimuth are −2.885° and −4.736°, respectively, at the end of the backward process. As a result, the proposed alignment scheme with the backward process extends the length of dada equivalently and can utilize the data adequately to decrease the alignment time. Moreover, the alignment performance of the proposed algorithm is also very favorable. The alignment accuracy can meet the requirement of polar alignment. Therefore, the proposed alignment algorithm can solve the problem of polar alignment caused by SINS mechanization.
On the other hand, in Figure 13 and Figure 14, it is clear that the pseudo position returns to the initial position, when the initial alignment is completed. Although there are still some differences, the influence of the drift errors of pseudo position on the initial transformation matrix C g g p ( 0 ) is very small and can be neglected since the initial pseudo latitude is always zero. Therefore, the initial transformation matrix C g g p ( 0 ) can be utilized to obtain the strapdown attitude matrix with a simple transformation logic and without influence of position errors. Moreover, with the simple transformation logic, there have less the cumbersome derivation and low computational complexity in the process of transformation, thereby decreasing the complexity of program and the burden of computer. Here, we should notice that the azimuth angle relative to g -frame is the angle between the vehicle direction and the 126° longitude when the initial latitude is 90°. In addition, in Table 1, we can see that the mean and variance of alignment errors based on the proposed algorithm are also smaller as compared with the traditional method. This is because the proposed algorithm eliminates the influence of position errors on transformational accuracy. Therefore, the proposed algorithm has superior performance. On the other hand, since the strapdown attitude matrix relative to g -frame is obtained at the end of alignment phase, it would also be more convenient to access any one of the following polar navigation systems without any additional workloads and complexities. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, the proposed polar alignment algorithm would be advantageous to unify the forms of polar alignment and the external observation information of filter.
Furthermore, comparing Figure 11 and Figure 12, more time is required for polar alignment with the initial latitude 90°. This is the other key problem, which must be confronted in polar alignment. The degree of observability of azimuth is smaller in polar regions since the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude. As a result, the future efforts will focus on the improvement of the degree of observability of azimuth, thereby improving the polar alignment performance more.

6. Conclusions

In this paper, from the perspective of global inertial navigation, a generalized transverse Earth frame built based on the initial position of vehicle (namely, pseudo-Earth frame) is proposed as the attitude reference frame for polar alignment. Then, a polar alignment algorithm based on the pseudo-Earth frame and backward process is proposed to implement polar alignment. Through the proposed algorithm, the attitude matrix relative to local geographic frame is obtained without influence of position errors when the initial alignment is completed, and there is less cumbersome derivation and low computational complexity. As a result, it is very convenient to access any one of the following polar navigation systems without any additional workloads and complexities. Moreover, it would also be advantageous to unify the forms of polar alignment and the external observation information of filter, thereby further unifying the implemented form of the reference information acquired from the external navigation device. In addition, the backward process can still accelerate the filtering estimation and decrease the alignment time. Finally, the semi-physical static simulation and in-motion tests with large azimuth misalignment angle assisted by UKF are carried out to verify the proposed polar alignment algorithm, and the results exhibit the superior performance of the proposed alignment algorithm in polar regions.

Acknowledgments

This work is supported by Ministry of Science and Technology of China (grant number 2014DFR10010).

Author Contributions

All the authors contributed to this work. This idea is originally from the discussion among a team consisting of Meng liu, Yanbin Gao and Guangchun Li. Guangchun Li and Yanbin Gao provided much experimental support; Meng Liu conducted the simulation and experiments and finished writing this manuscript; Xingxing Guang conducted some work on the modification of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bame, S.J.; Hones, E.W.; Akasofu, S.I. Geomagnetic storm particles in the high-latitude magneto tail. J. Geophys. Res. 1971, 76, 7566–7583. [Google Scholar] [CrossRef]
  2. Ngwira, C.M.; Mckinnell, L.A.; Cilliers, P.J. GPS phase scintillation observed over a high-latitude Antarctic station during solar minimum. J. Atmos. Solar-Terr. Phys. 2010, 72, 718–725. [Google Scholar] [CrossRef]
  3. Andalsvik, Y.L.; Jacobsen, K.S. Observed high-latitude GNSS disturbances during a less-than-minor geomagnetic storm. Radio Sci. 2014, 49, 1277–1288. [Google Scholar] [CrossRef]
  4. Anderson, E.W. I—Navigation in Polar Regions. J. Navig. 1957, 10, 156–161. [Google Scholar]
  5. Titterton, D.; Weston, J. Strapdown Inertial Navigation Technology, 2nd ed.; IET Inc.: Chicago, IL, USA, 2004. [Google Scholar]
  6. Core, G.D.; Nastro, V. A World-wide Mechanization in Inertial Navigation Systems. J. Navig. 1986, 39, 3. [Google Scholar] [CrossRef]
  7. Li, Q.; Ben, Y.; Yu, F. System reset of transversal strapdown INS for ship in polar region. Measurement 2015, 60, 247–257. [Google Scholar] [CrossRef]
  8. Salychev, O.S. Applied Inertial Navigation Problems and Solutions; The BMSTU Press: Moscow, Russia, 2004. [Google Scholar]
  9. Maclure, K.C. Polar Navigation. Arctic 1949, 24, 183–194. [Google Scholar] [CrossRef]
  10. Fox, W.A.W. Transverse Navigation: An Alternative to the Grid System. J. Navig. 1949, 2, 25–35. [Google Scholar] [CrossRef]
  11. Qi, Z.; Yue, Y.Z.; Zhang, X.D. Indirect grid inertial navigation mechanization for transpolar aircraft. J. Chin. Inert. Technol. 2014, 22, 18. [Google Scholar]
  12. Li, W.; Wu, W.; Wang, J.; Lu, L. A Fast SINS Initial Alignment Scheme for Underwater Vehicle Applications. J. Navig. 2013, 66, 181–198. [Google Scholar] [CrossRef]
  13. Xiong, J.; Guo, H.; Yang, Z.H. A Two-Position SINS Initial Alignment Method Based on Gyro Information. Adv. Space Res. 2014, 53, 1657–1663. [Google Scholar] [CrossRef]
  14. Silson, P.M.G. Coarse Alignment of a Ship’s Strapdown Inertial Attitude Reference System Using Velocity Loci. IEEE Trans. Instrum. Meas. 2011, 60, 1930–1941. [Google Scholar] [CrossRef]
  15. Qin, Y.Y. Inertial Navigation; Science Inc.: Beijing, China, 2006. [Google Scholar]
  16. Wu, Y.; Zhang, H.; Wu, M. Observability of Strapdown INS Alignment: A Global Perspective. IEEE Trans. Aerosp. Electron. Syst. 2011, 48, 78–102. [Google Scholar]
  17. Li, Q.; Ben, Y.; Sun, F. Transversal Strapdown INS and Damping Technology for Marine in Polar Region. In Proceedings of the 2014 IEEE/ION Position, Location and Navigation Symposium—PLANS 2014, Monterey, CA, USA, 5–8 May 2014; pp. 1365–1370. [Google Scholar]
  18. Sun, F.; Yang, X.; Ben, Y.; Liu, S.; Li, M. Polar Transfer Alignment Research Based on Inverse Coordinate System. J. Proj. Rocket. Missiles Guid. 2014, 34, 179–182. [Google Scholar]
  19. Yao, Y.Q.; Xu, X.S.; Li, Y. Transverse Navigation under the Ellipsoidal Earth Model and its Performance in both Polar and Non-polar areas. J. Navig. 2015, 1, 1–18. [Google Scholar] [CrossRef]
  20. Cheng, J.; Wang, T.; Guan, D. Polar transfer alignment of shipborne SINS with a large misalignment angle. Meas. Sci. Technol. 2016, 27, 35–101. [Google Scholar] [CrossRef]
  21. Li, Q.; Ben, Y.; Yu, F. Transversal Strapdown INS based on Reference Ellipsoid for Vehicle in Polar Region. IEEE Trans. Veh. Technol. 2016, 65, 7791–7795. [Google Scholar] [CrossRef]
  22. Yan, G.; Yan, W.; Xu, D. On Reverse Navigation Algorithm and Its Application to SINS Gyro-Compass In-Movement Alignment. In Proceedings of the 2008 CCC 27th Chinese Control Conference, Kunming, China, 16–18 July 2008; pp. 724–729. [Google Scholar]
  23. Hu, J.; Cheng, X. A new In-Motion Initial Alignment for Land-Vehicle SINS/OD Integrated System. In Proceedings of the Position, Location and Navigation Symposium-PLANS 2014, Monterey, CA, USA, 5–8 May 2014; pp. 407–412. [Google Scholar]
  24. Alahyari, A.; Rozbahani, S.G.; Habibzadeh, A. INS/DVL Positioning System using Kalman Filter. Aust. J. Basic Appl. Sci. 2011, 5, 1123–1129. [Google Scholar]
  25. Feng, W.; Qin, Y.; Zhang, J. In-Flight Alignment Algorithm for Radar Aided SINS in Polar Regions. J. Northwest. Polytech. Univ. 2014, 32, 131–136. [Google Scholar]
  26. Liu, M.; Gao, Y.; Li, G. An Improved Alignment Method for the Strapdown Inertial Navigation System (SINS). Sensors 2016, 16, 621. [Google Scholar] [CrossRef]
  27. Gao, W.; Ben, Y.; Zhang, X.; Li, Q.; Yu, F. Rapid Fine Strapdown INS Alignment Method under Marine Mooring Condition. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 2887–2896. [Google Scholar] [CrossRef]
  28. Chang, L.; Hu, B.; Li, Y. Backtracking Integration for Fast Attitude Determination-Based Initial Alignment. IEEE Trans. Instrum. Meas. 2015, 64, 795–803. [Google Scholar] [CrossRef]
  29. Estimation Techniques for Low-Cost Inertial Navigation. Available online: http://webee.technion.ac.il/control/info/Projects/Students/2010/Guy%20Rosenthal%20and%20Nir%20Ben%20Zrihem/WWW/unSDINS_dissertation.pdf (accessed on 31 May 2005).
  30. Chang, L.; Hu, B.; Li, A. Strapdown inertial navigation system alignment based on marginalised unscented kalman filter. IET Sci. Meas. Technol. 2013, 7, 128–138. [Google Scholar] [CrossRef]
  31. Fang, J.; Yang, S. Study on Innovation Adaptive EKF for In-Flight Alignment of Airborne POS. IEEE Trans. Instrum. Meas. 2011, 60, 1378–1388. [Google Scholar]
  32. Julier, S.J.; Uhlmann, J.K.; Durrant-Whyte, H.F. A new approach for filtering nonlinear systems. Am. Control Conf. Proc. IEEE 1995, 3, 1628–1632. [Google Scholar]
  33. Li, J.; Song, N.; Yang, G. Fuzzy adaptive strong tracking scaled unscented Kalman filter for initial alignment of large misalignment angles. Rev. Sci. Instrum. 2016, 87, 790–798. [Google Scholar] [CrossRef]
Figure 1. Thepseudo coordinate system.
Figure 1. Thepseudo coordinate system.
Sensors 17 01416 g001
Figure 2. The corresponding relations of time index between forward process and backward process.
Figure 2. The corresponding relations of time index between forward process and backward process.
Sensors 17 01416 g002
Figure 3. The proposed polar alignment scheme.
Figure 3. The proposed polar alignment scheme.
Sensors 17 01416 g003
Figure 4. The implementation of the proposed alignment algorithm.
Figure 4. The implementation of the proposed alignment algorithm.
Sensors 17 01416 g004
Figure 5. The block diagram of polar alignment with unscented Kalman filter (UKF).
Figure 5. The block diagram of polar alignment with unscented Kalman filter (UKF).
Sensors 17 01416 g005
Figure 6. The original outputs of gyro and the modified outputs of gyro with 85° latitude.
Figure 6. The original outputs of gyro and the modified outputs of gyro with 85° latitude.
Sensors 17 01416 g006
Figure 7. The semi-physical static alignment errors with pseudo navigation mechanization at 85° latitude.
Figure 7. The semi-physical static alignment errors with pseudo navigation mechanization at 85° latitude.
Sensors 17 01416 g007
Figure 8. The test vehicle.
Figure 8. The test vehicle.
Sensors 17 01416 g008
Figure 9. The alignment errors of vehicle test in low latitude regions.
Figure 9. The alignment errors of vehicle test in low latitude regions.
Sensors 17 01416 g009
Figure 10. The calculated pseudo position of vehicle test in low latitude regions.
Figure 10. The calculated pseudo position of vehicle test in low latitude regions.
Sensors 17 01416 g010
Figure 11. The in-motion alignment errors based on the proposed polar alignment algorithm with initial latitude 87°.
Figure 11. The in-motion alignment errors based on the proposed polar alignment algorithm with initial latitude 87°.
Sensors 17 01416 g011
Figure 12. The in-motion alignment errors based on the proposed polar alignment algorithm with initial latitude 90°.
Figure 12. The in-motion alignment errors based on the proposed polar alignment algorithm with initial latitude 90°.
Sensors 17 01416 g012
Figure 13. The calculated pseudo-position with initial latitude 87°.
Figure 13. The calculated pseudo-position with initial latitude 87°.
Sensors 17 01416 g013
Figure 14. The calculated pseudo-position with initial latitude 90°.
Figure 14. The calculated pseudo-position with initial latitude 90°.
Sensors 17 01416 g014
Table 1. Statistics of alignment errors based on the two methods at 87° latitude in 1200 s.
Table 1. Statistics of alignment errors based on the two methods at 87° latitude in 1200 s.
ErrorsThe Proposed MethodTraditional Transverse Method
MeanVarianceMeanVariance
Pitch0.003108°5 × 10−70.005417°4.25 × 10−5
Roll−0.04375°7.63 × 10−40.07080°2.6 × 10−4
Azimuth−2.869°3.26 × 10−4−3.127°8.64 × 10−3

Share and Cite

MDPI and ACS Style

Gao, Y.; Liu, M.; Li, G.; Guang, X. Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions. Sensors 2017, 17, 1416. https://doi.org/10.3390/s17061416

AMA Style

Gao Y, Liu M, Li G, Guang X. Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions. Sensors. 2017; 17(6):1416. https://doi.org/10.3390/s17061416

Chicago/Turabian Style

Gao, Yanbin, Meng Liu, Guangchun Li, and Xingxing Guang. 2017. "Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions" Sensors 17, no. 6: 1416. https://doi.org/10.3390/s17061416

APA Style

Gao, Y., Liu, M., Li, G., & Guang, X. (2017). Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions. Sensors, 17(6), 1416. https://doi.org/10.3390/s17061416

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