Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Study of Impact Damage in PVA-ECC Beam under Low-Velocity Impact Loading Using Piezoceramic Transducers and PVDF Thin-Film Transducers
Previous Article in Journal
The Design of a Single-Bit CMOS Image Sensor for Iris Recognition Applications
Previous Article in Special Issue
A Novel Grid SINS/DVL Integrated Navigation Algorithm for Marine Application
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 18
Issue 2
10.3390/s18020670
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Error Analysis of the K-Rb-21Ne Comagnetometer Space-Stable Inertial Navigation System

by
Qingzhong Cai
1,*,
Gongliu Yang
1,
Wei Quan
1,
Ningfang Song
1,
Yongqiang Tu
1 and
Yiliang Liu
2
1
School of Instrument Science and Opto-electronics Engineering, Beihang University, Beijing 100191, China
2
National Remote Sensing Center of China, Beijing 100036, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 670; https://doi.org/10.3390/s18020670
Submission received: 7 November 2017 / Revised: 7 February 2018 / Accepted: 21 February 2018 / Published: 24 February 2018
(This article belongs to the Special Issue Inertial Sensors for Positioning and Navigation)
Browse Figures
Figure 1
<p>Structure of the space-stable platform.</p> "> Figure 2
<p>Mechanization of space-stable platform Inertial Navigation System (INS).</p> "> Figure 3
<p>Single error test result in a numerical simulation of the mathematical model.</p> "> Figure 4
<p>Numerical simulation based on semi-physical data.</p> "> Figure 5
<p>Raw data of the comagnetometer.</p> "> Figure 6
<p>Position error in numerical simulation based on semi-physical data.</p> ">
Versions Notes

Abstract

:
According to the application characteristics of the K-Rb-21Ne comagnetometer, a space-stable navigation mechanization is designed and the requirements of the comagnetometer prototype are presented. By analysing the error propagation rule of the space-stable Inertial Navigation System (INS), the three biases, the scale factor of the z-axis, and the misalignment of the x- and y-axis non-orthogonal with the z-axis, are confirmed to be the main error source. A numerical simulation of the mathematical model for each single error verified the theoretical analysis result of the system’s error propagation rule. Thus, numerical simulation based on the semi-physical data result proves the feasibility of the navigation scheme proposed in this paper.

1. Introduction

With the development of quantum physics, the atomic gyroscope, based on quantum sensing technology, has gained widespread attention. The Inertial Navigation System (INS) with atomic gyroscopes has great potential in long-term independent navigation applications, such as underwater deep-sea robots [1,2]. There are two types of atomic gyroscope in theory [3]: Atomic Interferometry Gyroscopes (AIGs) [4,5,6,7], and Atomic Spin Gyroscopes (ASGs) [8,9,10,11,12,13,14,15]. The ASG does not need much laser power and linewidth, has no need for high vacuum, and is insensitive to gravity. As a result, it has been widely researched in recent years. The research on ASG has transitioned from the theoretical research phase to the prototype development and experimenting phase [16,17,18,19]. Fang et al. deduced the dynamics equation of the ASG [20]. The random drift of the atomic spin gyroscope was modelled and optimised by Quan et al. [21], Duan et al. [22], and Zou et al. [23]. The scale factor of the ASG was calibrated by Zou et al. [24]. The cross-talk effect of a dual-axis ASG was suppressed by Jiang et al. [25].
Based on the atomic Spin-Exchange Relaxation Free (SERF), the hybrid optical pumping can measure with high sensitivity the magnetic field or angular velocity [26,27,28]. Significant breakthroughs in the research of the K-Rb-21Ne comagnetometer prototype have provided a base for studying comagnetometer navigation methods. Because the angular velocity measurement is based on the spin atomic’s inertia and precession, it has a small measuring range just like the mechanical gyroscope does, the signal detection is performed using linear polarised light, which causes the comagnetometer to not be torqued on as an optical gyroscope. Thus, the comagnetometer cannot be used in either a traditional strap-down INS or a north-slaved local-level platform INS. Thus, it is necessary to carry out theoretical research on the comagnetometer INS, verify its principle and navigation accuracy to promote comagnetometer optimal design, and analyse the early stage requirements of the comagnetometer.
The comagnetometer INS is aimed at long-term, high-accuracy independent navigation for the submarine. The comagnetometer is similar to the electrostatic gyroscope (ESG), and thus can adopt a space-stable platform system as its navigation scheme [29,30,31,32]. The gyroscope works in a space-stable inertial frame when using a space-stable platform that is suitable for the comagnetometer INS, because it is conducive to the gyroscope’s optimal accuracy without a large measuring range or torqueing to the gyroscope. Existing research results on space-stable INS mainly apply to missile navigation systems [33,34,35,36], which does not apply to the long-term navigation of the comagnetometer INS. Thus, this paper aims to design a special space-stable navigation scheme for the comagnetometer INS, present the requirements of the comagnetometer prototype, and analyse the theoretical accuracy by determining the error propagation function.
This paper is organised as follows. In section two, the structure of a space-stable platform and its requirements of the comagnetometer are discussed. In section three, the mechanization equation of the comagnetometer space-stable INS is introduced. In section four, the theoretical accuracy of the comagnetometer INS is analysed. In section five, the semi-physical simulation results are presented as well as the accuracy verification. The last section concludes and discusses the findings.

2. Characteristic and Error Model of the Comagnetometer

2.1. Characteristic

The dynamics of the comagnetometer can be described by Bloch equations [3,16,20]. In steady state, keeping the leading order terms of each magnetic field and light shift, if the pump beam is in the z-direction and the probe beams are in the x and y direction, the Bloch equations are [16]:
P x e = P z e γ e R t o t e · Ω y γ n P y e = P z e γ e R t o t e · Ω x γ n
where, P x e , P y e and P z e are the polarization of electron spin of the alkali metal atom in x-, y- and z-direction; γ e is the electron gyromagnetic ratio and γ n is the nuclear gyromagnetic ratio; R t o t e is the total relaxation rate of electron spin. Ω x and Ω y are the measurement rotation in x- and y-direction.
Since the probe beam in x-direction is used to detect P x e , it can sense the rotation in y-direction. Similarly, the probe beam in y-direction is used to detect P y e and senses the rotation in x-direction. If both probe beams are used in ASG, the ASG is actually a dual-axis gyroscope [17].
The transfer function of the comagnetometer deduced by the Bloch equations can be simplified as [20]:
G ( s ) = k ω n 2 s 2 − 2 ζ ω n s + ω n 2
where, k , ζ and ω n are comprised of ASG-related parameters.
It can be seen from Equations (1) and (2) that the transfer function of the comagnetometer is similar to a dual-axis rate gyroscope. As previously mentioned, the angular velocity measurement is based on the spin atomic inertia and precession, similar to the mechanical gyroscope.
According to the measurement principle of the comagnetometer, the nuclear spin angular momentum (NSAM) helps the electron spin angular momentum (ESAM) to avoid being affected by the magnetic field, but the strength of the compensation field generated by the NSAM is limited to a quick automatic compensation of the magnetic field [1]. The measuring range cannot be very large—with an experimental value of ±5°/s on the existing prototype—thus, the comagnetometer does not apply to a strap-down INS like the traditional optical gyroscope does.
The signal detection of the comagnetometer is based on linear polarised light detection, which is very similar to the optical gyroscope. Regarding the torqueing condition, because the comagnetometer’s rotor consists of innumerable atomics, which is different from a traditional mechanic rotor, high accuracy torqueing cannot be exerted on spin atomic like a mechanical gyroscope does. Therefore, it is similar to an optical gyroscope in that the comagnetometer does not apply to a traditional north-slaved local-level platform INS, which requires torqueing on the gyroscope.

2.2. Error Model

Taking the interference factors of magnetic interference and light into account, the leading terms of transverse electron spin polarization can be simply expressed as [25]:
P x e = − γ e P z e R t o t e D s · [ Ω y γ + γ e Ω x γ R t o t e ( δ B z + L z ) − ( δ B z + L z ) γ e L x − R t o t e L y ] P y e = γ e P z e R t o t e D s · [ Ω x γ − γ e Ω y γ R t o t e ( δ B z + L z ) − R t o t e L x + ( δ B z + L z ) γ e L y ] ,
where, D s = R t o t e 2 + γ e ( δ B z + L z ) 2 , γ = γ e γ n / ( γ e − Q γ n ) , Q is the electron slowing-down factor due to hopping between hyperfine levels, δ B z is the magnetic interference in z-axis direction, L x , L y and L z are the light shifts in the x-, y,- and z-axis direction.
The error model can be transformed as
P x e = S y Ω y + α y x Ω x + ε y P y e = S x Ω x + α x y Ω y + ε x
where, S x = γ e P z e R t o t e γ D s and S y = − γ e P z e R t o t e γ D s are the scale factor, α x y = − γ e 2 P z e γ D s ( δ B z + L z ) and α y x = − γ e 2 P z e R t o t e γ D s ( δ B z + L z ) are the misalignment, ε x = − γ e P z e R t o t e 2 D s L x + γ e 2 P z e R t o t e ( δ B z + L z ) D s L y and ε y = γ e 2 P z e R t o t e ( δ B z + L z ) D s L x + γ e P z e R t o t e 2 D s L y are the bias of comagnetometer.

3. Comagnetometer Space-Stable Platform

3.1. Structure of the Space-Stable Platform

As explained above, the application characteristic of the comagnetometer does not match the scheme of the widely researched north-slaved, local-level platform INS, or the strap-down INS. The space-stable INS adopts a space-stable platform instead of exerting precessional torque on the gyroscope. On the space-stable platform, the gyroscope can operate freely to help the comagnetometer reach the optimal accuracy. Thus, the space-stable platform navigation scheme is considered.
The structure of the space-stable platform is shown in Figure 1. The platform is composed of the stable element, inner gimbal, middle gimbal, and outer gimbal from the inside out. Two dual-axis comagnetometers and an accelerometer triad are installed on the stable element. One comagnetometer has an atomic spin direction parallel to the polar-axis (called a polar-axis gyroscope), and the other one has an atomic spin direction parallel to the equator (called the equator gyroscope). The polar-axis gyroscope produces a two-degree-of-freedom angle, or angular velocity signals, to control the inner and middle gimbal axis. The equator gyroscope produces a one-degree-of-freedom angle, or angular velocity signals, to control the stable element axis. The outer gimbal axis is taken as a redundant axis controlled by the angle sensor output signal on the inner gimbal axis, or stable element axis; it is used to avoid a gimbal lock when the three axes are located in the same plane. The comagnetometers are not torqued on and the stable element functions in a space-stable state.

3.2. Discussion of the Requirements for the Comagnetometer Prototype

At present, the comagnetometer prototype operates in an angular velocity measurement mode based on gyroscopic precession. There are two problems that apply to the servo control of the space-stable platform of the comagnetometer that must be solved. One problem is dynamic performance, including the damping factor and the bandwidth. Although the sampling rate is set higher, the bandwidth of the comagnetometer prototype cannot be higher than 10 Hz, which cannot ensure the stability of the platform; yet according to the design experience of the platform INS, it requires a bandwidth higher than 100 Hz. The other problem is measurement error accumulation. Because the platform is controlled by the angular velocity signal, the attitude error of the platform accumulates during the whole navigation process. It requires the comagnetometer to have a much higher stability.
In the future, the comagnetometer prototype can be improved and will operate in an angle measurement mode based on gyroscopic inertia. As the signal can be detected by linearly polarised light, the measuring frequency will be high enough (1 kHz) to ensure the platform’s stability.

4. Error Analysis of the Comagnetometer Space-Stable INS

4.1. Mechanization Equation

Before introducing the mechanization of the space-stable platform system, the main coordinate frames are given first:
(a) Inertial frame (i-frame)
The inertial system (i-frame) is a right-hand orthogonal coordinate system with the origin at the geocentre. The Xi-axis is parallel to the intersecting line of the equatorial plane, and the local meridian plane at initial time (t = 0), pointing outside from the geocentre; the Zi-axis coincides with the Earth’s spin axis.
(b) Earth frame (e-frame)
The Earth system (e-frame) is also a right-hand orthogonal coordinate system with the origin at the geocentre. The three axes are fixed with the Earth; the Xe-axis is parallel to the intersecting line of the local meridian plane at initial time (t = 0) and the equatorial plane, pointing outside from the centre of the Earth; the Ze-axis coincides with the Earth’s spin axis.
(c) Geographic frame (n-frame)
In this paper, the geographic frame is defined with the north-, east-, down-pointing axes: the origin is the centroid of the vehicle; the Xn-axis, also known as the north (N) axis, is perpendicular to the meridian plane which contains the vehicle, directly toward the north; the Yn-axis, also known as the east (E) axis, points to geodetic east; and the Zn-axis, also known as the down (D) axis, is directly downward along the geocentric vertical through the vehicle.
Since the space-stable platform does not torque on the gyroscope, the platform frame (p-frame) coincides with the i-frame, and the mechanization can take either the i-frame or the e-frame as the navigation calculation frame. In this paper, the e-frame is adopted, and the mechanization diagram of the space-stable platform INS is shown in Figure 2 [30,31,32,33,34,35,36].
In the e-frame, the navigation calculation equation of the INS is:
r ˙ e = v e v ˙ e = C p e f p − 2 [ ω i e e ] v e + g e C ˙ p e = − [ ω i e e ] C p e + C p e [ ω d p ]
where, [ • ] stands for the skew symmetric matrix of a small angle vector, for example, if α = [ α 1 α 2 α 3 ] T , then [ α ] = [ 0 − α 3 α 2 α 3 0 − α 1 − α 2 α 1 0 ] ; r e is the position vector from the center of earth to the vehicle in e-frame; v e is the velocity vector in e-frame; C p e is the transition matrix from p-frame to e-frame; f p is the specific force measured by the accelerometer triad in p-frame; ω i e e is the Earth's rotation angular velocity vector in e-frame; g e is the gravitational acceleration vector in e-frame; ω d p is the angular velocity of comagnetometer drift in p-frame.
The position (longitude λ, latitude L and altitude h) in n-frame can be iterated using the position vector r e = [ r x e r y e r z e ] T in e-frame:
[ L λ h ] = [ arctan ( r z e ( r x e ) 2 + ( r y e ) 2 ( 1 + e 2 R N sin L r z e ) ) λ 0 + arctan ( r y e r x e ) ( r x e ) 2 + ( r y e ) 2 cos L − R N ] ,
where e is the eccentricity of the Earth; RN is the radius of curvature of the prime vertical plane, λ0 is the initial longitude.
The velocity vector in the n-frame can be solved using the velocity vector v e in the e-frame:
v n = C e n v e ,
where C e n = [ − cos λ sin L − sin λ sin L cos L − sin λ cos λ 0 − cos λ cos L − sin λ cos L − sin L ] .
Then the attitude matrix C b n is calculated using the matrix C b p measured using platform frame angle and the matrix C p e from p-frame to e-frame:
C b n = C e n C p e C b p .
Finally, the attitude can be calculated using C b n .

4.2. Error Model of the System

According to the space-stable mechanization equation, the error equation of the comagnetometer INS can be derived by the Ψ -angle error model as:
δ r ˙ e = δ v e δ v ˙ e = − 2 [ ω i e e ] δ v e − [ ω i e e ] [ ω i e e ] δ r e + [ f e ] Ψ e + C p e δ f p Ψ ˙ e = − [ ω i e e ] Ψ e + C p e δ ω d p ,
where, δ r e is the position error vector of the vehicle in e-frame; δ v e is the velocity error vector relative to the Earth in e-frame; Ψ e = [ φ x e φ y e φ z e ] T presents the drift error angle of the p-frame relative to the calculation frame (e-frame); δ f p stands for measurement error vector of the accelerometer triad in p-frame; and δ ω d p is the angular velocity error of comagnetometer drift.
Multiplying C e n by both sides of (9), and considering the displacement angular velocity ρ n from n-frame to e-frame, the Ψ -angle error equation can be transformed to n-frame:
δ r ˙ n = δ v n − [ ρ n ] δ r n δ v ˙ n = − ( 2 [ ω i e n ] + [ ρ n ] ) δ v n − [ ω i e n ] [ ω i e n ] δ r n + [ f n ] Ψ n + δ f n Ψ ˙ n = − [ ω i e n + ω e n n ] Ψ n + C p n δ ω d p ,
where, ρ n = [ v E R N + h − v N R M + h − v E R N + h tan L ] T , R M is radius of curvature in meridian; R N is radius of curvature in prime vertical; h is the altitude; v E and v N are the velocity components in east- and north- direction.

4.3. Error Propagation Rule

Simplifying (10), and omitting the cross-coupling items of the vertical and horizontal loops, the position error equation of horizontal loop is [31]:
δ r ¨ N = − 2 ω i e sin L δ r ˙ E − ( ω s 2 − ω i e 2 sin 2 L ) δ r N + g φ y e + δ f N δ r ¨ E = 2 ω i e sin L δ r ˙ N − ( ω s 2 − ω i e 2 ) δ r E − g ( − φ x e sin L + φ z e cos L ) + δ f E ,
where ω s = g / R is the Schuler oscillating angular frequency.
It can be seen that the space-stable platform INS has similar Schuler, Foucault and Earth oscillation with the north-slaved local-level INS.
Because the Schuler and Foucault oscillation errors are all restrained by horizontal damping in long-term navigation, they are omitted from the equation solution for simplification. Assuming that the initial attitude error is Ψ 0 e = [ φ x 0 e φ y 0 e φ z 0 e ] T , the comagnetometer measurement error is δ ω d p = [ δ ω x δ ω y δ ω z ] T , the accelerometer measurement error is δ f p = [ δ f x δ f y δ f z ] T , the time domain solution of the space-stable platform INS’s position error is:
δ L = − φ x 0 e sin ( ω i e t ) + φ y 0 e cos ( ω i e t ) − δ ω x t sin ( ω i e t ) + δ ω y t cos ( ω i e t )   + δ f x sin L cos ( ω i e t ) / g + δ f y sin L sin ( ω i e t ) / g − δ f z cos L / g δ λ = φ x 0 e cos ( ω i e t ) sin L + φ y 0 e sin ( ω i e t ) sin L − φ z 0 e + δ ω x t cos ( ω i e t ) sin L   + δ ω y t sin ( ω i e t ) sin L − δ ω z t cos L + δ f x sin ( ω i e t ) / g − δ f y cos ( ω i e t ) / g
Then the main error source of the space-stable platform INS in long-term navigation is analysed from (12). Position error-related items caused by initial attitude error and accelerometer bias are both bounded errors, while the position error caused by bias infinitely increases with the navigation time increasing. Specifically, the z-axis bias causes a longitude error linearly increasing with time, which is similar to the one caused by up- and north-axis biases in the north-slaved, local-level INS. The x- and y-axis biases cause longitude and latitude errors oscillating with the Earth period, and their amplitudes increase with the navigation time linearly; this is an error characteristic that has never been identified in the north-slaved, local-level INS. It is because the inertial instruments in the two types of INSs have different orientations relative to the Earth’s gravitational field.
According to (4), the measurement error of the comagnetometer in the INS is:
δ ω d p = ( δ S + M ) ω i e e + ε .
where, δ S = [ δ S x 0 0 0 δ S y 0 0 0 δ S z ] is the scale factor error matrix, M = [ 0 α x y α x z α y x 0 α y z α z x α z y 0 ] is the misalignment matrix and ε = [ ε x ε y ε z ] T is the bias matrix of comagnetometer.
In e-frame, the Earth's rotation angular velocity is ω i e e = [ 0 0 ω i e ] T . With (12) and (13), we can extract the divergent errors as:
δ L = − ( α x z ω i e + ε x ) t sin ( ω i e t ) + ( α y z ω i e + ε y ) t cos ( ω i e t ) δ λ = ( α x z ω i e + ε x ) t cos ( ω i e t ) sin L + ( α y z ω i e + ε y ) t sin ( ω i e t ) sin L − ( δ S z ω i e + ε z ) t cos L
It can be seen from (14) that the main error source of the comagnetometer space-stable platform INS in long-term navigation includes the biases of the three gyroscopes, the scale factor error of the z-axis gyroscope, and the misalignment of the x- and y-axis gyroscope in the z-direction. The difference from a north-slaved, local-level INS is that the latitude error of the space-stable INS can also increase unboundedly with time in long-term navigation, as does the amplitude of the Earth periodical position error.
The divergent position errors caused by comagnetometer errors are summarised in Table 1.

5. Simulation Tests

5.1. Numerical Simulation of the Mathematical Model

To verify the results of error propagation analysis, a numerical simulation of the mathematical model for each error term is carried out. By setting x-, y- and z-axis biases as 0.0001°/h, the accelerometer biases as 1 μg, the scale factor error as 1 ppm, and the misalignment as 2″ respectively, a 168 h (7 days) long-term navigation simulation is performed. The single error test results are presented in Figure 3a–h. It can be seen from the figures that the x- and y-axis biases cause longitude and latitude errors oscillating with the Earth period, and that the oscillations amplify with the navigation time linearly; the z-axis bias causes a longitude error linearly increasing with time, which is not affected by the Earth period. The accelerometer biases do not cause any unbounded error increasing with time. The scale factor error of the z-axis causes a longitude error similar to that of the z-axis bias; the misalignment due to the x- and y-axis non-orthogonal with the z-axis causes longitude and latitude errors, similar to that of the x- and y-axis bias. The simulation result verified the error propagation and its period analysis result found in Section 4.3.
The max position errors of the numerical simulation are shown in Table 2. It can be seen that the simulation results are consistent with the calculation result by taking simulation conditions into the error expressions in Table 1, which proves the analysis result about the error propagation rule. Thus, it can be inferred that the three-axis biases, the scale factor of the z-axis and the misalignment due to the x- and y-axis non-orthogonal with the z-axis are the main error source in long-term navigation of the comagnetometer INS.

5.2. Numerical Simulation Based on Semi-Physical Data

The design of the semi-physical simulation test is shown in Figure 4. The random error is extracted from the raw data of the comagnetometer prototype. A group of 7-h comagnetometer raw data with a sampling frequency of 200 Hz is shown in Figure 5. The simulation data is generated by superimposing the numerical simulation data of the mathematical model and the real random error. Taking this semi-physical data as input, the off-line simulation of the system is performed.
First, add the comagnetometer random error with the mean square error of 0.01°/h into the simulation data. The position error result of the 7-h navigation result is shown in Figure 6a. The position accuracy can reach 1 nmile/4 h, which preliminarily proves the feasibility of using the comagnetometer in a space-stable platform INS. Then, decrease the comagnetometer random error to the mean square error of 0.0001°/h proportionally and recycle the raw data. The position error result of 168 h (7 days) navigation is shown in Figure 6b. The position accuracy can reach 1 nmile/7 days, which preliminarily proves the potential of comagnetometer use in long-term navigation.

6. Conclusions and Discussion

In this paper, a space-stable navigation scheme, as well as its mechanization, is proposed according to the comagnetometer’s application characteristic. By establishing the error model of the space-stable system and analysing its error propagation rule, some useful conclusions are derived about theoretical positioning accuracy in long-term independent navigation using the comagnetometer INS: The x- and y-axis biases and the misalignment due to the x- and y-axis non-orthogonal with the z-axis cause longitude and latitude errors oscillating with the Earth period, and amplify linearly with navigation time; the bias and the scale factor error of the z-axis causes longitude error increasing with time. A numerical simulation of the mathematical model results of each single error verified the theoretical analysis result of the system’s error propagation. The three-axis biases, the scale factor of the z-axis and the misalignment of the x- and y-axis non-orthogonal with the z-axis are the main error source in long-term navigation of the comagnetometer INS. For an INS with the position accuracy of 1 nmile/7 days, the bias, the scale factor and the misalignment of the INS should meet the performance specification requirement with an accuracy better than 0.0001°/h, 1 ppm and 2″, respectively. A numerical simulation of the semi-physical data proves that the navigation scheme proposed in this paper is feasible.
This paper provides a theoretical basis for comagnetometer INS design and also presents the requirements of the comagnetometer prototype. The bandwidth of the angular velocity measurement mode and the development of an angle measurement mode represent key points for future research.

Acknowledgments

This work was supported by the National Key R&D Program of China (Grant No. 2016YFB0501600), the National Equipment Pre-research Foundation (Grant No. 6140517020101) and the National Natural Science Foundation of China (Grants No. 61227902, No. 61673041, and No. 61374210).

Author Contributions

Q. Cai, G. Yang and N. song conceived and desined the experiments; W. Quan contributed data; Y. Tu designed the structure; Q. Cai performed the experiments and analyzed the data; Y. Liu wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fang, J.; Qin, J. Advances in atomic gyroscopes: A view from inertial navigation applications. Sensors 2012, 12, 6331–6346. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, H.G.; Williams, T.C. Strategic inertial navigation systems-High-accuracy inertially stabilized platforms for hostile environments. IEEE Control Syst. 2008, 28, 65–85. [Google Scholar] [CrossRef]
  3. Dong, H.; Fang, J.; Qin, J.; Chen, Y. Analysis of the electrons-nuclei coupled atomic gyroscope. Optics Communications. Opt. Commun. 2011, 284, 2886–2889. [Google Scholar] [CrossRef]
  4. Durfee, D.S.; Shaham, Y.K.; Kasevich, M.A. Long-term stability of an area-reversible atom-interferometer sagnac gyroscope. Phys. Rev. Lett. 2005, 97, 387–392. [Google Scholar] [CrossRef] [PubMed]
  5. Dubetsky, B.; Kasevich, M.A. Atom interferometer as a selective sensor of rotation or gravity. Phys. Rev. A 2006, 74, 154. [Google Scholar] [CrossRef]
  6. Gustavson, T.L.; Landragin, A.; Kasevich, M.A. Rotation sensing with a dual atom-interferometer Sagnac gyroscope. Class. Quantum Gravity 2000, 17, 2385–2398. [Google Scholar] [CrossRef]
  7. Dutta, I.; Savoie, D.; Fang, B.; Venon, B.; Alzar, C.; Geiger, R.; Landragin, A. Continuous cold-atom inertial sensor with 1 nrad/sec rotation stability. Phys. Rev. Lett. 2016, 116, 183003. [Google Scholar] [CrossRef] [PubMed]
  8. Du, J.; Rong, X.; Zhao, N.; Wang, Y.; Yang, J.; Liu, B. Preserving electron spin coherence in solids by optimal dynamical decoupling. Nature 2010, 461, 1265–1268. [Google Scholar] [CrossRef] [PubMed]
  9. Shi, F.; Zhang, Q.; Wang, P.; Sun, H.; Wang, J.; Rong, X.; Chen, M.; Ju, C.; Reinhard, F.; Chen, H.; et al. Protein imaging. Single-protein spin resonance spectroscopy under ambient conditions. Science 2015, 347, 1135–1138. [Google Scholar] [CrossRef] [PubMed]
  10. Kornack, T.W.; Ghosh, R.K.; Romalis, M.V. Nuclear spin gyroscope based on an atomic comagnetometer. Phys. Rev. Lett. 2005, 95, 230801. [Google Scholar] [CrossRef] [PubMed]
  11. Allred, J.C.; Lyman, R.N.; Kornack, T.W.; Romalis, M.V. High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation. Phys. Rev. Lett. 2002, 89, 130801. [Google Scholar] [CrossRef] [PubMed]
  12. Kitching, J.; Donley, E.A.; Hodby, E.; Shkel, A.; Eklund, E.J. Compact Atomic Magnetometer and Gyroscope Based on a Diverging Laser Beam. U.S. Patent 7872473, 18 January 2011. [Google Scholar]
  13. Fang, J.; Chen, Y.; Zou, S.; Liu, X.; Hu, Z.; Quan, W.; Yuan, H.; Ding, M. Low frequency magnetic field suppression in an atomic spin co-magnetometer with a large electron magnetic field. J. Phys. B Atomic Mol. Opt. Phys. 2016, 49, 065006. [Google Scholar] [CrossRef]
  14. Li, R.; Fan, W.; Jiang, L.; Duan, L.; Quan, W.; Fang, J. Rotation sensing using a K-Rb-21Ne comagnetometer. Phys. Rev. A 2016, 94, 032109. [Google Scholar] [CrossRef]
  15. Dong, H.; Gao, Y. Comparison of Compensation Mechanism Between an NMR Gyroscope and an SERF Gyroscope. IEEE Sens. J. 2017, 17, 4052–4055. [Google Scholar] [CrossRef]
  16. Qin, J.; Fang, J.; Wan, S. Experimental design of a dual axis Atomic Spin Gyroscope. In Proceedings of the 2012 8th IEEE International Symposium on Instrumentation and Control Technology (ISICT), London, UK, 11–13 July 2012; pp. 24–27. [Google Scholar]
  17. Ochoa, J.E.; Ochoa, J.E.; Correa, M.M.; Valencia, A.M.; Mcewen, J.G.; Gallo, J.A.; Bilo, G. Stable three-axis nuclear spin gyroscope in diamond. Phys. Rev. A 2012, 86, 12184–12188. [Google Scholar]
  18. Fang, J.; Qin, J.; Wan, S.; Chen, Y.; Li, R. Atomic spin gyroscope based on 129Xe-Cs comagnetometer. Chin. Sci. Bull. 2013, 58, 1512–1515. [Google Scholar] [CrossRef]
  19. Zhang, H.; Zou, S.; Chen, X. Optimal modeling on magnetic shielding cylinder of atomic spin gyroscope considering transverse shielding factor. In Proceedings of the 2015 IEEE Metrology for Aerospace (MetroAeroSpace), Benevento, Italy, 4–5 June 2015; pp. 6–8. [Google Scholar]
  20. Fang, J.; Wan, S.; Yuan, H. Dynamics of an all-optical atomic spin gyroscope. Appl. Opt. 2013, 52, 7220–7227. [Google Scholar] [CrossRef] [PubMed]
  21. Quan, W.; Lv, L.; Liu, B. Modeling and optimizing of the random atomic spin gyroscope drift based on the atomic spin gyroscope. Rev. Sci. Instrum. 2014, 85, 113104. [Google Scholar] [CrossRef] [PubMed]
  22. Duan, L.; Quan, W.; Jiang, L.; Fan, W.; Ding, M.; Hu, Z.; Fang, J. Common-mode noise reduction in an atomic spin gyroscope using optical differential detection. Appl. Opt. 2017, 56, 7734–7740. [Google Scholar] [CrossRef] [PubMed]
  23. Zou, S.; Zhang, H.; Chen, X. Modeling and filter algorithm analysis of all-optical atomic spin gyroscope’s random drift. In Proceedings of the 2015 IEEE Metrology for Aerospace (MetroAeroSpace), Benevento, Italy, 4–5 June 2015; pp. 207–219. [Google Scholar]
  24. Zou, S.; Zhang, H.; Chen, X.; Chen, Y.; Fang, J. A novel calibration method research of the scale factor for the all-optical atomic spin inertial measurement device. J. Opt. Soc. Korea 2015, 19, 415–420. [Google Scholar] [CrossRef]
  25. Jiang, L.; Quan, W.; Li, R.; Duan, L.; Fan, W.; Wang, Z.; Liu, F.; Xing, L.; Fang, J. Suppression of the cross-talk effect in a dual-axis K-Rb-21Ne comagnetometer. Phys. Rev. A 2017, 95, 062103. [Google Scholar] [CrossRef]
  26. Quan, W.; Li, Y.; Liu, B. Simultaneous measurement of magnetic field and inertia based on hybrid optical pumping. EPL 2015, 110, 60002. [Google Scholar] [CrossRef]
  27. Fang, J.; Wang, T.; Zhang, H.; Li, Y.; Cai, H.-W. In-situ measurement of magnetic field gradient in a magnetic shield by a spin-exchange relaxation-free magnetometer. Chin. Phys. B 2015, 24, 060702. [Google Scholar] [CrossRef]
  28. Karaulanov, T.; Savukov, I.; Kim, Y.J. Spin-exchange relaxation-free magnetometer with nearly parallel pump and probe beams. Meas. Sci. Technol. 2016, 27, 055002. [Google Scholar] [CrossRef]
  29. Elton, O.L.; Moore, J.P. Marine ESG Navigation as a Capability for the Present. Navigation 1973, 20, 126–136. [Google Scholar] [CrossRef]
  30. Dushman, A.; Sandberg, H. Accuracy of ESG monitor/sins inertial navigation system. In Proceedings of the Ocean 73 IEEE International Conference on Engineering in the Ocean Environment, Seattle, WA, USA, 25–28 September 1973; pp. 602–610. [Google Scholar]
  31. Gao, Z. Error Propagation Property of Inertial Navigation System. In Inertial Navigation System Technology, 1st ed.; Tsinghua University Press: Beijing, China, 2012; pp. 279–289. [Google Scholar]
  32. Wu, Q.; Han, F. New optimal approach to space-stable inertial navigation system. In Proceedings of the 2011 10th International Conference on Electronic Measurement & Instruments (ICEMI), Chengdu, China, 16–19 August 2011; pp. 296–299. [Google Scholar]
  33. Nash, R.A.; Levine, S.A.; Roy, K.J. Error analysis of space-stable inertial navigation systems. IEEE Trans. Aerosp. Electron. Syst. 1971, 7, 617–629. [Google Scholar] [CrossRef]
  34. Hutchinson, C.E.; Nash, R.A. Comparison of Error Propagation in Local-Level and Space-Stable Inertial Systems. IEEE Trans. Aerosp. Electron. Syst. 1971, 7, 1138–1142. [Google Scholar] [CrossRef]
  35. Nash, R.A.; Hutchinson, C.E. Altitude damping of space-stable inertial navigation systems. IEEE Trans. Aerosp. Electron. Syst. 1973, 9, 18–27. [Google Scholar] [CrossRef]
  36. Hutchinson, C.E.; Fagan, J.H. Kalman filter design considerations for space-stable inertial navigation systems. IEEE Trans. Aerosp. Electron. Syst. 1973, 9, 306–319. [Google Scholar] [CrossRef]
Sensors 18 00670 g001 550
Figure 1. Structure of the space-stable platform.
Figure 1. Structure of the space-stable platform.
Sensors 18 00670 g001
Sensors 18 00670 g002 550
Figure 2. Mechanization of space-stable platform Inertial Navigation System (INS).
Figure 2. Mechanization of space-stable platform Inertial Navigation System (INS).
Sensors 18 00670 g002
Sensors 18 00670 g003 550
Figure 3. Single error test result in a numerical simulation of the mathematical model.
Figure 3. Single error test result in a numerical simulation of the mathematical model.
Sensors 18 00670 g003
Sensors 18 00670 g004 550
Figure 4. Numerical simulation based on semi-physical data.
Figure 4. Numerical simulation based on semi-physical data.
Sensors 18 00670 g004
Sensors 18 00670 g005 550
Figure 5. Raw data of the comagnetometer.
Figure 5. Raw data of the comagnetometer.
Sensors 18 00670 g005
Sensors 18 00670 g006 550
Figure 6. Position error in numerical simulation based on semi-physical data.
Figure 6. Position error in numerical simulation based on semi-physical data.
Sensors 18 00670 g006
Table
Table 1. Position errors caused by comagnetometer errors.
Table 1. Position errors caused by comagnetometer errors.
Error TermLatitude ErrorLongitude Error
ε x − t sin ( ω i e t ) t cos ( ω i e t ) sin L
ε y t cos ( ω i e t ) t sin ( ω i e t ) sin L
ε z - − t cos L
α x z − ω i e t sin ( ω i e t ) ω i e t cos ( ω i e t ) sin L
α y z ω i e t cos ( ω i e t ) ω i e t sin ( ω i e t ) sin L
δ S z - − ω i e t cos L
Table
Table 2. Max position error in numerical simulation of the mathematical model.
Table 2. Max position error in numerical simulation of the mathematical model.
Error TermMax Latitude ErrorMax Longitude Error
ε x = 0.0001 ∘ / h 0.98 nmile0.65 nmile
ε y = 0.0001 ∘ / h 1.01 nmile0.65 nmile
ε z = 0.0001 ∘ / h 0.02 nmile0.8 nmile
α x z = 2 ″ 1.41 nmile0.95 nmile
α y z = 2 ″ 1.45 nmile0.92 nmile
δ S z = 1 p p m 0.02 nmile1.18 nmile

Share and Cite

MDPI and ACS Style

Cai, Q.; Yang, G.; Quan, W.; Song, N.; Tu, Y.; Liu, Y. Error Analysis of the K-Rb-21Ne Comagnetometer Space-Stable Inertial Navigation System. Sensors 2018, 18, 670. https://doi.org/10.3390/s18020670

AMA Style

Cai Q, Yang G, Quan W, Song N, Tu Y, Liu Y. Error Analysis of the K-Rb-21Ne Comagnetometer Space-Stable Inertial Navigation System. Sensors. 2018; 18(2):670. https://doi.org/10.3390/s18020670

Chicago/Turabian Style

Cai, Qingzhong, Gongliu Yang, Wei Quan, Ningfang Song, Yongqiang Tu, and Yiliang Liu. 2018. "Error Analysis of the K-Rb-21Ne Comagnetometer Space-Stable Inertial Navigation System" Sensors 18, no. 2: 670. https://doi.org/10.3390/s18020670

APA Style

Cai, Q., Yang, G., Quan, W., Song, N., Tu, Y., & Liu, Y. (2018). Error Analysis of the K-Rb-21Ne Comagnetometer Space-Stable Inertial Navigation System. Sensors, 18(2), 670. https://doi.org/10.3390/s18020670

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