Next Article in Journal
Robust Adaptive Finite-Time Synergetic Tracking Control of Delta Robot Based on Radial Basis Function Neural Networks
Previous Article in Journal
Physical Interpretation of Nanofluid (Copper Oxide and Silver) with Slip and Mixed Convection Effects: Applications of Fractional Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method

1
Navigation Research Center, Nanjing University of Aeronautics & Astronautics, Nanjing 211106, China
2
School of Air Transportation and Engineering, Nanhang Jincheng College, Nanjing 211156, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(21), 10862; https://doi.org/10.3390/app122110862
Submission received: 4 October 2022 / Revised: 22 October 2022 / Accepted: 24 October 2022 / Published: 26 October 2022

Abstract

:
In order to improve the continuity and smoothness of transpolar flight and optimize autonomous navigation performance, a SINS/SRS/CNS (strapdown inertial navigation system/spectral redshift navigation system/celestial navigation system) multi-information fusion global autonomous navigation method based on parameter conversion was studied for this article. The global autonomous navigation scheme based on multi-information fusion was designed. The principle of spectral redshift navigation was studied. On this basis, the system equations of the SINS/SRS/CNS multi-information fusion global autonomous navigation system were established in the middle–low latitudes and high latitudes. Furthermore, the navigation and filter parameter conversion relationships between the geographic navigation coordinate frame and the grid navigation coordinate frame were derived. The simulation and experiment verified that the SINS/SRS/CNS multi-information fusion global autonomous navigation method with parameter conversion can effectively improve the accuracy and smoothness and realize non-oscillation switching in transpolar navigation. In the vehicle experiment, the proposed algorithm improved the horizontal position accuracy by more than 29% compared with the multi-information fusion global autonomous navigation method without filter parameter conversion.

1. Introduction

High-altitude and long-endurance flight requires UAVs (unmanned aerial vehicles) to have global autonomous navigation capability [1]. The unique geographical environment of the polar regions puts forward higher requirements for the continuity and smoothness of navigation. SINS (the strapdown inertial navigation system) is applicable to global navigation because of its concealment, real-time operation, and independence from complex environments. However, with the increase in latitude, the east–north–up navigation coordinate frame extensively used in the middle–low latitudes is ineffective because of meridian convergence [2,3,4]. Consequently, the transverse geographic coordinate frame that rotates the traditional longitude and latitude lines by 90° and the grid coordinate frame that is based on the Greenwich meridian plane are presently used [5,6,7]. In addition, the SINS errors gradually accumulate, making it difficult to accomplish a UAV high-precision and long-endurance global flight mission. Consequently, it is essential to introduce other navigation information for SINS correction [8].
In recent years, many studies have been carried out on polar navigation. A fault-tolerant grid SINS/DVL (Doppler velocity log)/USBL (ultra-short baseline) integrated algorithm was proposed to improve the reliability and accuracy of polar navigation for ocean space applications [9]. The integrated navigation scheme of transverse SINS/DVL based on the virtual sphere model was proposed, and the navigation accuracy could meet the requirements of an autonomous underwater vehicle [10]. The rotation modulation inertial navigation system algorithm under a grid coordinate frame was studied, and the applied basis of polar navigation was established [11]. An interacting dual-model-based adaptive filter was proposed in INS/GNSS/DVL integrated navigation to overcome nonstationary noise and certain non-Gaussian factors, and ship sailing simulations demonstrated that the proposed architecture and algorithms enhanced navigational performance in the Arctic [12]. A SINS/CNS/GPS integrated navigation system that combines the azimuth navigation algorithm and the grid navigation algorithm was obtained, and the feasibility in the high latitudes was verified [13]. Related studies have achieved remarkable results. However, autonomous navigation technology of multi-information fusion for an airborne application background and transregional switching between middle and high latitudes need further research.
To meet the demand for fully autonomous, high-accuracy, and long-endurance navigation, correcting SINS errors with CNS (the celestial navigation system) has already become one research focus of aerospace engineering [14]. The U.S. Air Force has developed various typical airborne SINS/CNS integrated navigation products to meet the requirements of strategic reconnaissance missions [15]. SRS (the spectral redshift navigation system) is a novel autonomous navigation technology that uses the Doppler shift effect of astronomical spectra to realize aircraft navigation with preferable real-time operation and autonomy [16]. The fusion of SINS with SRS and CNS could guarantee the reliability of the navigation system while maintaining a high level of autonomy.
In this article, a SINS/SRS/CNS multi-information fusion global autonomous navigation method based on parameter conversion is proposed. Our main contribution is to establish the SINS/SRS/CNS multi-information fusion system equations in different navigation coordinate frames and to derive the navigation and filter parameter conversion relationships between different coordinate frames, ensuring the high autonomy and precision of the UAV global navigation system. The sections of this article are arranged as follows: The design of the scheme of multi-information fusion global autonomous navigation is detailed in Section 2. Section 3 analyzes the principle of spectral redshift navigation. The system equations of the middle–low-latitude geographic coordinate frame and the high-latitude grid coordinate frame are also established in Section 3. Based on Section 3, Section 4 is devoted to deriving the navigation and filter parameter conversion relationships between different navigation coordinate frames. Section 5 describes the simulation and experiment that were performed to verify the proposed algorithm. Section 6 discusses the main conclusion of this article.

2. Design of the Multi-Information Fusion Global Autonomous Navigation Scheme

For the SINS/SRS/CNS multi-information fusion global autonomous navigation, SINS was selected as the major system, while SRS and CNS were selected as auxiliary systems. The velocity provided by SRS and the latitude and longitude provided by CNS are used to correct the navigation errors of SINS. Moreover, a BA (barometric altimeter) is used to compensate for the inability of CNS to measure altitude.
A high-altitude and long-endurance UAV needs to traverse between polar and nonpolar regions during global flights. Consequently, it is necessary to consider the navigation and filter parameter conversion relationships between different navigation coordinate frames.
In summary, the multi-information fusion global autonomous navigation scheme for UAVs is shown in Figure 1.

3. The Multi-Information Fusion Autonomous Navigation System Based on SINS/SRS/CNS

First, the frames used in this section are defined:
i frame—inertial coordinate frame. It is fixed with the inertial space and can be used as a reference to describe the Earth and carrier angular motion.
g frame—east–north–up geographic coordinate frame. It is the navigation coordinate frame in the nonpolar regions.
e frame—Earth coordinate frame. It is fixed with the Earth and can be used as a reference to describe the carrier position.
b frame—body coordinate frame. The direction relationship between the b frame and the navigation system is the carrier attitude.
G frame—grid coordinate frame. It is the navigation coordinate frame in the polar regions. The angle between the grid north and the geographic north is the grid angle, σ .
c frame—calculation coordinate frame. It is the navigation coordinate frame described by SINS, which contains the platform error angle.
The relation diagram for the above frames is shown in Figure 2.

3.1. The Principle of Spectral Redshift Navigation

When the photon frequency emitted by a celestial body produces a color change in the spectrum, the spectrum emitted by the celestial body differs from the spectrum the carrier receives. The redshift is an essential parameter of the spectrum, including the relative velocity. The redshift, z , is defined as:
z = λ λ 0 λ 0 = f 0 f m f m
where λ 0 is the initial wavelength of the spectrum lines, λ is the measured wavelength; f 0 is the initial frequency of the spectral lines; and f is the measured frequency.
The redshift equation can be obtained from the equation:
1 + z = f 0 f m = 1 + | v | cos θ / c 1 | v | 2 / c 2
where v is the velocity of the carrier relative to the celestial body in the i frame; θ is the angle between the astro-carrier vector and v ; and c is the light velocity.
According to the ephemeris of the celestial bodies and the carrier attitude, combined with the measurement of spectral redshift, the flight velocity of the carrier can be obtained [17,18]. The basic principle of spectral redshift autonomous navigation is as follows. The schematic diagram of spectral redshift navigation is shown in Figure 3:
During flight, the carrier can detect the light signal emitted by the celestial body. From the Doppler principle, it can be known that the spectral frequency received by the carrier is not equal to the frequency the celestial body emitted, so the relative velocity between the carrier and the celestial body can be calculated by measuring the spectral redshift. When the number of non-collinear celestial bodies detected by the carrier is not less than three, the carrier velocity can be obtained by geometrical relations. The solution process is as follows [16,19,20].
According to Equation (2), when the carrier moves relative to the celestial body, it can be obtained that:
| v | cos θ = ( 1 + z ) c 2 | v | 2 c
The above equation is also the radial velocity of the carrier relative to the celestial body:
v r = | v | cos θ = ( 1 + z ) c 2 | v | 2 c
Equation (4) is applied to the first reference celestial body:
v r 1 = ( 1 + z 1 ) c 2 | v p v 1 | 2 c
where v r 1 is the radial velocity of the carrier relative to the first reference celestial body; z 1 is the redshift of the carrier relative to the first reference celestial body; v p is the velocity of the carrier in the i frame; and v 1 is the velocity of the first reference celestial body in the i frame.
The equation group can be listed by selecting three reference celestial bodies:
{ v r 1 = ( 1 + z 1 ) c 2 | v p v 1 | 2 c v r 2 = ( 1 + z 2 ) c 2 | v p v 2 | 2 c v r 3 = ( 1 + z 3 ) c 2 | v p v 3 | 2 c
According to the geometrical relation, v p and v r 1 , v r 2 , v r 3 satisfy:
{ v r 1 = ( v p v 1 ) u 1 v r 2 = ( v p v 2 ) u 2 v r 3 = ( v p v 3 ) u 3
where v 1 , v 2 , v 3 are the velocities of reference celestial bodies in the i frame, obtained from the ephemeris, and u 1 , u 2 , u 3 are the unit vectors between the celestial bodies and the carrier, observed by the star sensor.
By substituting Equation (7) into Equation(6), we obtain:
{ ( v p v 1 ) u 1 ( 1 + z 1 ) c 2 | v p v 1 | 2 + c = 0 ( v p v 2 ) u 2 ( 1 + z 2 ) c 2 | v p v 2 | 2 + c = 0 ( v p v 3 ) u 3 ( 1 + z 3 ) c 2 | v p v 3 | 2 + c = 0
Then, the carrier velocity in the i frame can be obtained.
The core equipment of SRS is the optical spectrum instrumentation. The spectral parameters can be measured by the optical spectrum instrumentation. Then, the redshift and the carrier velocity can be calculated. The calculated principle shows that the carrier velocity is obtained by geometrical relations, not by circulative iteration. This means that its error would not accumulate over time, and the correlation time is short. Therefore, the error of SRS can be regarded as white noise, and related studies have been conducted in this way [16,19,20].

3.2. The System Equations of the SINS/SRS/CNS Multi-Information Fusion Autonomous Navigation System in the Middle–Low Latitudes

3.2.1. State Equations

In the middle–low latitudes, the g frame was chosen as the navigation coordinate frame. The SINS error equations in the g frame can be obtained as follows [21].
The mathematical platform error angle equation can be denoted as:
ϕ ˙ g = ϕ g × ( ω i e g + ω e g g ) + ( δ ω i e g + δ ω e g g ) C b g ε b
where ϕ g is the mathematical platform error angle in the g frame; ω i e g is the projection of the Earth’s rotation angular velocity, ω i e , in the g frame; ω e g g is the angular velocity of the g frame to the e frame, projected in the g frame; δ ω i e g and δ ω e g g are the errors of ω i e g and ω e g g ; C b g is the direction cosine matrix from the b frame to the g frame; and ε b is the gyroscope drift in the b frame.
The velocity error equation can be denoted as:
δ v ˙ g = f g × ϕ g ( 2 δ ω i e g + δ ω e g g ) × v g ( 2 ω i e g + ω e g g ) × δ v g + C b g b
where δ v g is the velocity error in the g frame; f g represents the specific force in the g frame; and b is the accelerator drift in the b frame.
The position error, δ p g , comprises the latitude error, δ L ; the longitude error, δ λ ; and the altitude error, δ h . The equations of position error are as follows:
{ δ L ˙ = δ v N R M + h δ λ ˙ = δ v E R N + h sec L + v E R N + h sec L t g L δ L δ h ˙ = δ v U
where R M and R N are the radii of curvatures in meridian and prime vertical, respectively, and V E and V N are the east velocity and north velocity, respectively.
The gyroscope and accelerator drift error equations can be denoted as
{ ε ˙ b = 0 ˙ b = 0
According to Equations (9)–(12), select ϕ g , δ v g , δ p g , ε b , and b as the state vector:
X g ( t ) = [ ϕ g δ v g δ p g ε b b ] T
Then, the state equation of the SINS/SRS/CNS multi-information fusion autonomous navigation system in the middle–low latitudes can be established:
X ˙ g ( t ) = A g ( t ) X g ( t ) + G g ( t ) W g ( t )
where X g ( t ) is the state vector; A g ( t ) is the state transition matrix; G g ( t ) is the system noise coefficient matrix; and W g ( t ) is the system noise vector.

3.2.2. Measurement Equations

There are two sets of measurements: one is the SINS/SRS velocity measurement, and the other is the SINS/CNS position measurement.
(1)
SINS/SRS velocity measurement
Given that the velocity provided by SRS is in the i frame, converting it to the g frame in the middle–low latitudes is essential:
v S R S g = C i g v p = C e g C i e v p
where C i g is the direction cosine matrix from the i frame to the g frame; C i e is the direction cosine matrix from the i frame to the e frame; C e g is the direction cosine matrix from the e frame to the g frame; and v S R S g is the velocity in the g frame calculated by SRS.
The velocity differences between SINS and SRS were chosen as one measurement to construct the velocity measurement equation:
Z v g ( t ) = v S I N S g v S R S g = H v g ( t ) X g ( t ) + V v g ( t )
where v S I N S g is the velocity in the g frame calculated by SINS; H v g ( t ) = [ 0 3 × 3 diag [ 1 1 1 ] 0 3 × 9 ] is the measurement matrix of the velocity; and V v g ( t ) is the velocity noise.
(2)
SINS/CNS position measurement
The position differences between SINS, CNS, and BA were selected as the other measurement to construct the position measurement equation:
Z p g ( t ) = [ L S I N S L C N S λ S I N S λ C N S h S I N S h B A ] = H p g ( t ) X g ( t ) + V p g ( t )
where L S I N S , λ S I N S , h S I N S are the latitude, longitude, and altitude obtained by SINS; L C N S , λ C N S are the latitude and longitude obtained by CNS; h BA is the altitude obtained by BA; H p g ( t ) = [ 0 3 × 6 diag [ 1 1 1 ] 0 3 × 6 ] is the measurement matrix of the position; and V p g ( t ) is the position noise.

3.3. The System Equations of the SINS/SRS/CNS Multi-Information Fusion Autonomous Navigation System in the High Latitudes

3.3.1. State Equations

In the high latitudes, the G frame was selected as the navigation coordinate frame. The coordinates in the e frame are used to represent the position of the UAV. Therefore, the error equations are as follows [9]:
{ ϕ ˙ G = ϕ G × ( ω i e G + ω e G G ) + C v δ v G C b G ε b δ v ˙ G = f G × ϕ G + ( v G × C v ( 2 ω i e G + ω e G G ) × ) δ v G + C b G b δ P ˙ e = C G e δ v G C G e ( v G × ) C P δ P e
where ϕ G is the mathematical platform error angle in the G frame; v G is the velocity in the G frame, and δ v G is the error; δ P e is the position error in the e frame; ω i e G is the projection of the Earth’s rotation angular velocity, ω i e , in the G frame; ω e G G is the angular velocity of the G frame to the e frame, projected in the G frame; f G represents the specific force in the G frame; C b G is the direction cosine matrix from the b frame to the G frame; C G e is the direction cosine matrix from the G frame to the e frame;
C v = 1 ( R M + h ) ( R N + h ) [ 0 1 0 1 0 0 0 cot L sin σ 0 ] ; σ is the angle between the grid north and the geographic north; C P = ξ ( R M + h ) ( R N + h ) [ sin L 0 cos λ cos L sin λ sin L cos 2 L 1 / ξ 2 sin λ sin L cos L cos L sin L sin λ 0 cos 2 L sin λ cos λ ] , and ξ = 1 1 cos 2 L sin 2 λ .
The gyroscope and accelerator drift error equations are the same as Equation (12).
According to Equations (18) and (12), select ϕ G , δ v G , δ P e , ε b , and b as the state vector:
X G ( t ) = [ ϕ G δ v G δ p e ε b b ] T
Then, the state equation of the SINS/SRS/CNS multi-information fusion autonomous navigation system in the high latitudes can be established:
X ˙ G ( t ) = A G ( t ) X G ( t ) + G G ( t ) W G ( t )
where X G ( t ) is the state vector; A G ( t ) is the state transition matrix; G G ( t ) is the system noise coefficient matrix; and W G ( t ) is the system noise vector.

3.3.2. Measurement Equations

(1)
SINS/SRS velocity measurement
In the high latitudes, it is required to convert the velocity in the i frame calculated by SRS to the G frame:
v S R S G = C i G v p = C e G C i e v p
where C i G is the direction cosine matrix from the i frame to the G frame; C e G is the direction cosine matrix from the e frame to the G frame; and v S R S G is the velocity in the G frame calculated by SRS.
The velocity differences between SINS and SRS were chosen as one measurement to construct the velocity measurement equation:
Z v G ( t ) = v S I N S G v S R S G = H v G ( t ) X G ( t ) + V v G ( t )
where v S I N S G is the velocity in the G frame calculated by SINS; H v G ( t ) = [ 0 3 × 3 diag [ 1 1 1 ] 0 3 × 9 ] is the measurement matrix of the velocity; and V v G ( t ) is the velocity noise.
(2)
SINS/CNS position measurement
In the state equations of the high latitudes, the UAV position is represented by the coordinates in the e frame. Therefore, it is necessary to convert the latitude and longitude from CNS and the altitude from BA to the e frame:
{ x C N S = ( R N + h B A ) cos L C N S cos λ C N S y C N S = ( R N + h B A ) cos L C N S sin λ C N S z C N S = [ R N ( 1 f 2 ) + h B A ] sin L C N S
where x C N S ,   y C N S ,   z C N S are the position coordinates in the e frame obtained by CNS and BA.
The position differences between SINS, CNS, and BA in the e frame were selected as the other measurement:
Z p G ( t ) = [ x S I N S x C N S y S I N S y C N S z S I N S z C N S ] = H p G ( t ) X G ( t ) + V p G ( t )
where x S I N S ,   y S I N S ,   z S I N S are the position coordinates in the e frame obtained by SINS, H p G ( t ) = [ 0 3 × 6 diag [ 1 1 1 ] 0 3 × 6 ] is the measurement matrix of the position; and V p G ( t ) is the position noise.

4. Parameter Conversion of the SINS/SRS/CNS Multi-Information Fusion Autonomous Navigation System

4.1. Navigation Parameter Conversion

Different navigation coordinate frames are adopted for polar regions and nonpolar regions. Therefore, the navigation parameters need to be switched between the g frame and the G frame during long-duration and transpolar flights.
C g G is the direction cosine matrix from the g frame to the G frame, and its transpose matrix C G g is the direction cosine matrix from the G frame to the g frame:
C g G = [ cos σ sin σ 0 sin σ cos σ 0 0 0 1 ]
When the UAV enters the polar regions, the navigation coordinate frame switches from the g frame to the G frame. The polar attitude matrix, C b G , and the polar velocity, v G , can be obtained by C g G :
{ C b G = C g G C b g v G = C g G v g
On the other hand, when the UAV leaves the polar regions, the navigation coordinate frame switches from the G frame to the g frame. The nonpolar attitude matrix, C b g , and the nonpolar velocity, v g , can be obtained by C G g :
{ C b g = C G g C b G v g = C G g v G
The conversion relationship between the position coordinates in the e frame and latitude, longitude, and altitude can be obtained from Equation (23).

4.2. Filter Parameter Conversion

The navigation coordinate frame switches while the UAV enters/leaves the polar regions, and the filter structure and parameters also change. The system equations in the polar/nonpolar regions were established (Section 3), and the coefficient matrices in different coordinate frames were also determined. The system noise and measurement noise are only related to the sensor accuracy, not to the navigation coordinate frame. Beyond that, to switch between the different system equations, it is necessary to convert the filter state vector and the covariance matrix.
By comparing Equations (13) and (19), it can be known that the inertial sensor errors, ε b and b , are the same in the g frame and the G frame, so no conversion is required.
The conversion relationship between ϕ g and ϕ G can be derived as follows.
Due to the calculation error, the navigation coordinate frame described by SINS is the c frame, containing the platform error angle. The platform error angle, ϕ G , relative to the G frame is usually a slight angle, so the first-order approximation is feasible. It can be deduced that:
C G c = I ( ϕ G × )
Then,
δ C b G = C b c C b G = C G c C b G C b G = ( ϕ G × ) C b G
Derive both sides of C b G = C g G C b g
δ C b G = δ C g G C b g + C g G δ C b g = ( ϕ g G G × ) C g G C b g C g G ( ϕ g × ) C b g = ( ϕ g G G × ) C b G C g G ( ϕ g × ) C G g C b G
where ϕ g G G = [ 0 0 δ σ ] T ; δ σ = sin σ cos σ cos L sin L δ L + 1 cos 2 σ cos 2 L sin L δ λ .
By comparing Equations (29) and (30), we can obtain:
( ϕ G × ) = C g G ( ϕ g × ) C G g + ( ϕ g G G × )
This can be reorganized as:
ϕ G = C g G ϕ g + ϕ g G G = C g G ϕ g + [ 0 0 δ σ ] = C g G ϕ g [ 0 0 sin σ cos σ cos L sin L δ L + 1 cos 2 σ cos 2 L sin L δ λ ] = C g G ϕ g [ 0 0 0 0 0 0 sin σ cos σ cos L sin L 1 cos 2 σ cos 2 L sin L 0 ] [ δ L δ λ δ h ] = C g G ϕ g [ 0 0 0 0 0 0 sin σ cos σ cos L sin L 1 cos 2 σ cos 2 L sin L 0 ] δ p g C g G ϕ g T 1 δ p g
The conversion relationship between δ v g and δ v G is as follows:
δ v G = C g G δ v g + δ C g G v g = C g G δ v g ( ϕ g G G × ) C g G v g = C g G δ v g [ ( 0 0 δ σ ) × ] C g G v g = C g G δ v g [ 0 δ σ 0 δ σ 0 0 0 0 0 ] [ cos σ sin σ 0 sin σ cos σ 0 0 0 1 ] [ v E v N v U ] = C g G δ v g + [ sin σ v E cos σ v N cos σ v E sin σ v N 0 ] δ σ = C g G δ v g + [ sin σ v E cos σ v N cos σ v E sin σ v N 0 ] [ sin σ cos σ cos L sin L 1 cos 2 σ cos 2 L sin L 0 ] [ δ L δ λ δ h ] = C g G δ v g + [ sin σ v E cos σ v N cos σ v E sin σ v N 0 ] [ sin σ cos σ cos L sin L 1 cos 2 σ cos 2 L sin L 0 ] δ p g C g G δ v g + T 2 δ p g
The position in the middle–low latitudes is expressed by latitude, longitude, and altitude, with the error of δ p g = [ δ L δ λ δ h ] T . Meanwhile, the position in the high latitudes is expressed by the coordinates in the e frame, with the error of δ p e = [ δ x δ y δ z ] T .
Depending on the relation between latitude, longitude, altitude and the coordinates in the e frame described in Equation (23), it can be obtained that:
{ δ x = δ [ ( R N + h ) cos L cos λ ] = ( R N + h ) ( sin L cos λ δ L cos L sin λ δ λ ) + cos L cos λ δ h δ y = δ [ ( R N + h ) cos L sin λ ] = ( R N + h ) ( sin L sin λ δ L + cos L cos λ δ λ ) + cos L sin λ δ h δ z = δ [ [ R N ( 1 f 2 ) + h ] sin L ] = [ R N ( 1 f 2 ) + h ] cos L δ L + sin L δ h
Equation (34) can be written as:
δ p e = [ ( R N + h ) sin L cos λ ( R N + h ) cos L sin λ cos L cos λ ( R N + h ) sin L sin λ ( R N + h ) cos L cos λ cos L sin λ [ R N ( 1 f 2 ) + h ] cos L 0 sin L ] δ p = T 3 δ p g
When the UAV enters the polar regions, the navigation coordinate frame switches from the g frame to the G frame. Combining Equations (32), (33), and (35), the polar state vector can be obtained as:
X G ( t ) = [ C g G 0 T 1 0 0 0 C g G T 2 0 0 0 0 T 3 0 0 0 0 0 I 3 × 3 0 0 0 0 0 I 3 × 3 ] X g ( t ) = Φ X g ( t )
According to the definition of the variance matrix and Equation (36):
P G ( t ) = E { ( X G ( t ) X ^ G ( t ) ) ( X G ( t ) X ^ G ( t ) ) T } = E { ( Φ X g ( t ) Φ X ^ g ( t ) ) ( Φ X g ( t ) Φ X ^ g ( t ) ) T } = E { Φ ( X g ( t ) X ^ g ( t ) ) ( X g ( t ) X ^ g ( t ) ) T Φ T } = Φ E { ( X g ( t ) X ^ g ( t ) ) ( X g ( t ) X ^ g ( t ) ) T } Φ T = Φ P g ( t ) Φ T
where P g ( t ) is the covariance matrix in the g frame and P G ( t ) is the covariance matrix in the G frame.
Similarly, when the UAV leaves the polar regions, the state vector and covariance matrix in the g frame can be obtained as:
{ X g ( t ) = Φ - 1 X G ( t ) P g ( t ) = Φ - 1 P G ( t ) ( Φ T ) - 1
In summary, it is only necessary to switch the navigation coordinate frame in the software and algorithm according to the calculated latitude. The navigation hardware devices do not need any conversion.

5. Experiment and Discussion

5.1. Simulation and Analysis

To validate the performance of the SINS/SRS/CNS multi-information fusion global autonomous navigation method, a simulation and an analysis were performed, combined with the system equations established in Section 3. The parameters of the navigation sensors for the simulation are shown in Table 1.
The latitude of the original location was 87.5°. The UAV flew at uniform velocity and arrived at latitude 88° at 212 s. The navigation coordinate frame was switched from the g frame to the G frame. The proposed SINS/SRS/CNS multi-information fusion autonomous navigation method with filter parameter conversion was compared with that without filter parameter conversion. Figure 4a–c demonstrate the position error, velocity error, and attitude error generated by these two methods, respectively.
It can be seen from Figure 4 that the errors of multi-information fusion autonomous navigation without filter parameter conversion oscillate when entering the polar regions. In the attitude errors, the yaw error without the filter parameter conversion was the largest, and the oscillation amplitude reached up to 3’. Due to the low observability of the yaw during a uniform flight, the change in filter structure was prone to oscillation. The oscillation amplitude of the velocity error was 0.3 m/s. The oscillations in the velocity error and position error were smaller than that of the attitude error and could converge quickly. That is because the multi-information fusion autonomous navigation filter took the velocity provided by SRS and the position provided by CNS as the measurements. Therefore, the velocity and position could be corrected directly. During the conversion of the navigation coordinate frame, the errors of the multi-information fusion autonomous navigation method with the filter parameter conversion were lower, and there were no significant oscillations.
The RMSEs (root-mean-square errors) of the SINS/SRS/CNS multi-information fusion autonomous navigation with and without filter parameter conversion are shown in Figure 5.
Based on the results demonstrated in Figure 4 and Figure 5, the navigation performance of the multi-information fusion autonomous navigation method with filter parameter conversion was better than that without filter parameter conversion. It could realize smooth switching between different coordinate frames and improve the navigation accuracy.

5.2. Experiment and Analysis

In view of the restriction of experimental conditions, actual high-altitude and long-endurance flight experiments cannot be conducted in our research center. Considering that RTK (real-time kinematics) can provide high-accuracy velocity and position, the following experimental scheme was designed. The high-accuracy velocity and position offered by RTK were taken as the true values of the SRS velocity, CNS position, and BA altitude. The CNS position error, BA altitude error, and SRS velocity error simulated in Section 5.1 could be superimposed on the true values to obtain the experimental data. On this basis, a SINS/RTK integrated navigation vehicle experiment in a natural system was conducted to verify the parameter conversion method proposed in Section 4. The experimental vehicle is shown in Figure 6.
The parameters of the vehicle experiment are listed in Table 2.
The experiment trajectory is demonstrated by the red line in Figure 7. The approximate location was east longitude 118.790° and north latitude 31.939°. The navigation coordinate frame was the g frame in the beginning, and it switched to the G frame after 300 s. The transpolar navigation was imitated by coordinate system switching.
Figure 8 illustrates the positions generated by the different methods.
Figure 8 indicates that the filter parameter conversion method proposed in this article can achieve smooth switching between different coordinate systems without oscillation.
The horizontal position RMSEs of these two methods are shown in Figure 9.
Figure 9 also shows that the proposed multi-information fusion autonomous navigation method with filter parameter conversion can effectively reduce the navigation errors caused by coordinate frame switching. Compared with the method without filter parameter conversion, the horizontal position accuracy of the proposed method was enhanced by about 29%.

6. Conclusions

To improve the continuity and smoothness of high-altitude and long-endurance UAV transpolar flight and to optimize autonomous navigation performance, a SINS/SRS/CNS multi-information fusion global autonomous navigation method was proposed. The system equations in the middle–low latitudes and high latitudes were established. The navigation and filter parameter conversion relationships between different coordinate frames were derived. The simulation and experiment verified that the SINS/SRS/CNS multi-information fusion global autonomous navigation method with parameter conversion can effectively improve the accuracy and smoothness and realize non-oscillation switching in transpolar navigation. In the vehicle experiment, the proposed algorithm improved the horizontal position accuracy by more than 29% compared with the multi-information fusion global autonomous navigation method without filter parameter conversion.

Author Contributions

Conceptualization, B.Z. and Q.Z.; methodology, B.Z.; software, B.Z.; validation, B.Z., Q.Z. and C.G.; formal analysis, X.Z. and W.Q.; investigation, C.G. and X.Z.; resources, J.L.; data curation, C.G. and W.Q.; writing—original draft preparation, B.Z.; writing—review and editing, Q.Z.; visualization, B.Z.; supervision, Q.Z.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grants 61533008 and 61603181), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (No. KYCX18_0268), the Qing Lan Project of the Jiangsu Higher Education Institutions (No. 2022), and the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (No. 19KJD590001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank all the editors and anonymous reviewers for their helpful comments and valuable remarks.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sheridan, I. Drones and Global Navigation Satellite Systems: Current Evidence from Polar Scientists. R. Soc. Open Sci. 2020, 7, 191494. [Google Scholar] [CrossRef] [Green Version]
  2. Huang, L.; Xu, X.; Ge, H.R. Single-Axis Rotation Modulation Transverse SINS Based on Virtual Sphere Model in Polar Region. IEEE Sens. J. 2022, 22, 13442–13450. [Google Scholar] [CrossRef]
  3. Zhao, W.Y.; Zhang, G.T.; Liu, W. Working Performance Analysis on Polar Navigation of Semi-analytical Inertial Navigation System. In Proceedings of the 2021 IEEE Conference on Telecommunications, Optics and Computer Science (TOCS), Shenyang, China, 10–11 December 2021. [Google Scholar]
  4. Zhao, B.; Zeng, Q.H.; Liu, J.Y. A New Polar Alignment Algorithm Based on the Huber Estimation Filter with the Aid of BeiDou Navigation Satellite System. Int. J. Distrib. Sens. Netw. 2021, 17, 15501477211004115. [Google Scholar] [CrossRef]
  5. Zhao, W.Y.; Zhang, G.T.; Liu, W. Research on High-latitude Transfer Alignment Technology. In Proceedings of the 2021 International Conference on Intelligent Computing, Automation and Systems (ICICAS), Chongqing, China, 29–31 December 2021. [Google Scholar]
  6. Zhang, F.B.; Gao, X.H.; Song, W.B. A Vision Aided Initial Alignment Method of Strapdown Inertial Navigation Systems in Polar Regions. Sensors 2022, 22, 4691. [Google Scholar] [CrossRef] [PubMed]
  7. Cui, W.T.; Ben, Y.Y.; Zhang, H.X. A Review of Polar Marine Navigation Schemes. In Proceedings of the 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), Portland, OR, USA, 20–23 April 2020. [Google Scholar]
  8. Song, L.J.; Zhao, W.L.; Cheng, Y.X. Based on Grid Reference Frame for SINS/CNS Integrated Navigation System in the Polar Regions. Complexity 2019, 2019, 2164053. [Google Scholar]
  9. Zhao, L.; Kang, Y.Y.; Cheng, J.H. A Fault-Tolerant Polar Grid SINS/DVL/USBL Integrated Navigation Algorithm Based on the Centralized Filter and Relative Position Measurement. Sensors 2019, 19, 3899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Huang, L.; Xu, X.; Zhao, H.M. Transverse SINS/DVL Integrated Polar Navigation Algorithm Based on Virtual Sphere Model. Math. Probl. Eng. 2020, 2020, 8892750. [Google Scholar] [CrossRef]
  11. Zhao, C.L.; Wu, W.Q.; Lian, J.X. Research on Rotating Modulation Inertial Navigation System Error Characteristics Simulation Method in Polar Area. In Proceedings of the Guidance, Navigation & Control Conference (CGNCC), Yantai, China, 8–10 August 2014. [Google Scholar]
  12. Sun, J.L.; Wu, Z.L.; Yin, Z.D. Adaptive Filtering and Temporal Alignment Based Fusion Algorithm for Navigation Systems in the Arctic Region. IEEE Syst. J. 2018, 13, 2022–2033. [Google Scholar] [CrossRef]
  13. Song, L.J.; Yang, G.Q.; Zhao, W.L. The Inertial Integrated Navigation Algorithms in the Polar Region. Math. Probl. Eng. 2020, 2020, 5895847. [Google Scholar] [CrossRef]
  14. Yang, S.J.; Feng, W.W.; Wang, S. A SINS/CNS Integrated Navigation Scheme with Improved mathematical horizon reference. Measurement 2022, 195, 111028. [Google Scholar] [CrossRef]
  15. Tian, M. Review of Polar Integrated Navigation Algorithm. J. Phys. Conf. Ser. 2019, 1213, 032017. [Google Scholar]
  16. Wei, W.H.; Gao, Z.H.; Gao, S.S. A SINS/SRS/GNS Autonomous Integrated Navigation System Based on Spectral Redshift Velocity Measurements. Sensors 2018, 18, 1145. [Google Scholar] [CrossRef]
  17. Collins, J.; Conger, R. MANS-Autonomous Navigation and Orbit Control for Communications Satellites. In Proceedings of the 15th International Communicatons Satellite Systems Conference and Exhibit, San Diego, CA, USA, 28 February–3 March 1994. [Google Scholar]
  18. Tai, F.; Noerdlinger, P.D. A Low Cost Autonomous Navigation System. In Proceedings of the Guidance and Control Conference, Keystone, CO, USA, 4–8 February 1989. [Google Scholar]
  19. Gao, G.L.; Gao, S.S.; Hong, G.Y. A Robust INS/SRS/CNS Integrated Navigation System with the Chi-Square Test-Based Robust Kalman Filter. Sensors 2020, 20, 5909. [Google Scholar] [CrossRef]
  20. Gao, G.L.; Gao, S.S.; Hu, G.L. Spectral Redshift Observation-based SINS/SRS/CNS Integration with an Adaptive Fault-tolerant Cubature Kalman Filter. Meas. Sci. Technol. 2021, 32, 095103. [Google Scholar] [CrossRef]
  21. Chen, W.N.; Zeng, Q.H.; Liu, J.Y. Seamless Autonomous Navigation Based on the Motion Constraint of the Mobile Robot. Ind. Robot. 2017, 44, 178–188. [Google Scholar] [CrossRef]
Figure 1. The multi-information fusion global autonomous navigation scheme diagram for UAVs.
Figure 1. The multi-information fusion global autonomous navigation scheme diagram for UAVs.
Applsci 12 10862 g001
Figure 2. The relation between the frames.
Figure 2. The relation between the frames.
Applsci 12 10862 g002
Figure 3. The schematic diagram of spectral redshift navigation, where r 1 , r 2 , r 3 , r p represent the position vectors of celestial body 1, celestial body 2, celestial body 3, and carrier P in the inertial coordinate system, respectively; r 1 p , r 2 p , r 3 p represent the position vectors from three celestial bodies to the carrier; v r 1 , v r 2 , v r 3 are the radial velocities of the carrier relative to three reference celestial bodies; and v p is the velocity vector of the carrier in the i frame.
Figure 3. The schematic diagram of spectral redshift navigation, where r 1 , r 2 , r 3 , r p represent the position vectors of celestial body 1, celestial body 2, celestial body 3, and carrier P in the inertial coordinate system, respectively; r 1 p , r 2 p , r 3 p represent the position vectors from three celestial bodies to the carrier; v r 1 , v r 2 , v r 3 are the radial velocities of the carrier relative to three reference celestial bodies; and v p is the velocity vector of the carrier in the i frame.
Applsci 12 10862 g003
Figure 4. Error curves of the two methods: (a) position error curves; (b) velocity error curves; and (c) attitude error curves.
Figure 4. Error curves of the two methods: (a) position error curves; (b) velocity error curves; and (c) attitude error curves.
Applsci 12 10862 g004aApplsci 12 10862 g004b
Figure 5. The RMSEs of multi-information fusion autonomous navigation: (a) the RMSEs of the position; (b) the RMSEs of the velocity; and (c) the RMSEs of the attitude.
Figure 5. The RMSEs of multi-information fusion autonomous navigation: (a) the RMSEs of the position; (b) the RMSEs of the velocity; and (c) the RMSEs of the attitude.
Applsci 12 10862 g005aApplsci 12 10862 g005b
Figure 6. The experimental vehicle.
Figure 6. The experimental vehicle.
Applsci 12 10862 g006
Figure 7. The experimental trajectory.
Figure 7. The experimental trajectory.
Applsci 12 10862 g007
Figure 8. The positions of the two methods with and without filter parameter conversion.
Figure 8. The positions of the two methods with and without filter parameter conversion.
Applsci 12 10862 g008
Figure 9. The horizontal position RMSEs of multi-information fusion autonomous navigation.
Figure 9. The horizontal position RMSEs of multi-information fusion autonomous navigation.
Applsci 12 10862 g009
Table 1. Sensor rate settings.
Table 1. Sensor rate settings.
ParametersValue
Gyro constant bias (deg/h)0.01
Accelerometer constant bias (ug)100
Gyro angular random walk (deg/sqrt(h))0.001
Accelerometer velocity random walk (ug/sqrt(Hz))10
CNS position error (m)20
BA altitude error (m)10
SRS velocity error (m/s)0.5
Table 2. The parameters of the vehicle experiment.
Table 2. The parameters of the vehicle experiment.
ParametersValue
Gyro constant bias (deg/h)0.015
Accelerometer constant bias (ug)85
Gyro angular random walk (deg/sqrt(h))0.0007
Accelerometer velocity random walk (ug/sqrt(Hz)) 10
Horizontal positioning error (m)20
Altitude error (m)10
Velocity error (m/s)0.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, B.; Zeng, Q.; Liu, J.; Gao, C.; Zhu, X.; Qiao, W. A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method. Appl. Sci. 2022, 12, 10862. https://doi.org/10.3390/app122110862

AMA Style

Zhao B, Zeng Q, Liu J, Gao C, Zhu X, Qiao W. A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method. Applied Sciences. 2022; 12(21):10862. https://doi.org/10.3390/app122110862

Chicago/Turabian Style

Zhao, Bin, Qinghua Zeng, Jianye Liu, Chunlei Gao, Xiaoling Zhu, and Wei Qiao. 2022. "A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method" Applied Sciences 12, no. 21: 10862. https://doi.org/10.3390/app122110862

APA Style

Zhao, B., Zeng, Q., Liu, J., Gao, C., Zhu, X., & Qiao, W. (2022). A Novel SINS/SRS/CNS Multi-Information Fusion Global Autonomous Navigation Method. Applied Sciences, 12(21), 10862. https://doi.org/10.3390/app122110862

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