Next Article in Journal
Wearable Sensor Based Stooped Posture Estimation in Simulated Parkinson’s Disease Gaits
Next Article in Special Issue
An Improved Adaptive Compensation H∞ Filtering Method for the SINS’ Transfer Alignment Under a Complex Dynamic Environment
Previous Article in Journal
Social-Aware Driver Assistance Systems for City Traffic in Shared Spaces
Previous Article in Special Issue
The Elastic Contact and Stability Analysis of an Inertial Micro-Switch with a Spring Stationary Electrode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Monitoring Navigation Method for Cold Atom Interference Gyroscope

1
Institute of Navigation Instrument, Harbin Institute of Technology, Harbin 150001, China
2
College of Automation, Harbin Engineering University, Harbin 150001, China
3
State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(2), 222; https://doi.org/10.3390/s19020222
Submission received: 29 November 2018 / Revised: 3 January 2019 / Accepted: 4 January 2019 / Published: 9 January 2019
(This article belongs to the Special Issue Inertial Sensors and Systems 2018)

Abstract

:
The implementation principle of a typical three-pulse cold atom interference gyroscope is introduced in this paper. Based on its configuration and current research status, the problems of cold atom interference gyro are pointed out. The data-rate is insufficient, and it is difficult to achieve high dynamic measurement. Then, based on these two limitations, a novel design of the monitoring navigation system of the cold atom interference gyroscope (CAIG) and an intermediate-grade inertial measurement unit (IMU) was proposed to obtain the long-term position result without GPS signals, such as the Inertial Navigation System (INS) in underwater vehicles. While the CAIG was used as the external gyro, the bias of IMU and the misalignment angle between the CAIG-frame and the IMU-frame are obtained through filtering technique. The simulation test and field test demonstrated the improvements of the long-term positioning accuracy of the INS.

1. Introduction

Due to its advantages of complete autonomy, undisturbed, large-amount of output information, and robust real-time performance, the Inertial Navigation System (INS) has irreplaceable benefits in both military and civilian fields, and the application requirements continue to promote the development of inertial sensors. Since the accuracy of gyroscopes and accelerometers directly affects the accuracy of positioning and attitude output of Strapdown Inertial Navigation Systems (SINSs) [1,2,3,4], the development of inertial navigation technology can approximately rely on the advance of two inertial sensor technologies, especially for the gyroscope.
The first-generation gyroscope is represented by electrostatic rate gyro and mechanical gyro [5,6,7]. At present, the bias stability of mechanical gyros is better than 10 4 ( /h). The second-generation gyroscope is represented by the ring-laser gyroscope and the fiber-optic gyroscope (FOG), and the bias stability is better than 10 3 ( / h) [8,9]. As a representative of the newest generation of gyro, the cold atom interference gyroscope is based on the Sagnac effect of the atomic de Broglie wave, and its theoretical precision may reach 10 12 ( / h) [10]. Even though the current cold atomic gyro has many problems to be resolved from practical engineering applications, it has significant precision merits and considerable application prospects.
Since the first implementation of the atom-based interferometer at Stanford University in 1991 [11], atomic cooling and steering technology has developed rapidly. At present, with the atomic excitation, capture, manipulation, and the continuous reduction of the limit temperature of the atom cloud, atomic material wave related technologies have received extensive attention, and successively realized new prototypes such as atomic fountain clocks [12,13,14,15,16] and atomic gravimeters [17,18,19,20,21], have both surpassed the precision of traditional ones. It is worth noting that, due to the technology limitations, cold atom interferometers are mainly used for constant measurements, such as gravitation constants, local gravity, and measurement of basic physical constants [22,23,24,25,26]. Based on the measurement principle of CAIG, which is explained in Section 2, its dynamic measurement range is insufficient to perform a dynamic navigation test. In particular, the data-rate of CAIG cannot realize an independent inertial navigation system, and most of them are used to measure the earth’s rotation constant [10,27,28,29,30].
As the newest-generation sensor for rotation measurement, the CAIG with ultra-high precision potential leads to a breakthrough of IMUs, at the same time, it brings opportunities to solve long-term positioning problem of underwater vehicles without GPS signals. The CAIG‘s advantage lies in long-term stability, repeatability, and extremely high expected sensitivity. Unfortunately, its data-rate is too low to compose an INS till now. This insufficient data-rate will lead to algorithmic solution error. At the same time, the limited dynamic range of the CAIG cannot meet the requirement of accurately capturing the dynamic information of the carrier [31,32]. On the other hand, FOG is an advanced inertial sensor for rotation measurement of carriers. It obtains the advantages of high data-rate and large dynamic range. However, due to the limited improvement, the FOG cannot satisfy the long-term accuracy requirements in terms of bias stability and repeatability. To the best of our knowledge, there are no studies about combining CAIG and FOG. In this paper, a novel monitoring navigation method for CAIG is proposed, and dedicated to combining the benefits of the CAIG and FOG to improve the positioning accuracy of the SINS. When the underwater vehicle is in a lower dynamic environment, the CAIG is used to estimate the bias of the FOG. Otherwise, a mature SINS of FOG is used for navigation in a highly dynamic environment. This monitoring navigation method with lower requirements of CAIG’s data rate, and the CAIG’s measurement accuracy is fully utilized under the limited dynamic range.
In Section 2, the configuration of the CAIG is introduced, and the reasons for the lack of dynamic measurement and the data rate are explained. In Section 3, the design of the monitoring navigation system is introduced, and simulation analysis in Section 4 shows the conditions to realize a fast and accurate monitoring system. Section 5 shows the results of land vehicle experiments and ship experiments, with a significant improvement of the positioning accuracy. In order to reduce the amount of calculation and visually show the effect of this monitoring system [33,34], the classic Kalman filter is chosen to obtain the optimal solution of the time-invariant linear system.

2. Three-Pulse Cold Atom Interferometer for Rotational Measurement

The configuration of the rotation measurement based on the cold atom interferometer mainly includes three-pulse design and four-pulse design. In order to easily realize the triaxial measurement, the three-pulse opposite projection configuration is selected as an example in this paper. As shown in Figure 1, the implementation of CAIG is mainly divided into atom cloud preparation, cooling, state selection, and velocity selection process, interference process, and final atomic state detection. Among them, the imprisonment and manipulation of the atom cloud are all realized by laser control technique, which requires a precise optical system for frequency, phase and light emphasizing. In order to achieve better measurement accuracy and reduce the influence of the external environment, we hope that the external environment of the atomic interference gyroscope achieves vacuum, low temperature, magnetic shielding, and vibration isolation.

2.1. Atomic Cooling and Trapping Progress

In order to reduce the interference errors caused by the atom position error, velocity distribution and other factors to the interference result, the atom in the initial state of steam needs to be cooled to the extremely low temperature of the μ K -level and trapped as an easy-to-manipulate atom cloud [35,36,37]. The existence time is prolonged through this progress.
When a precise laser illuminates an atom, as long as the oscillation frequency of the atom is equal to the frequency of the laser, the atom absorbs photons and undergoes a transition to reach an excited state. Depending on the principle of energy conservation, the momentum of the atom itself is greatly reduced. On the other hand, the atomic state in the excited state is extremely difficult to maintain, then it will also spontaneously radiate out photons. These photons are emitted to all directions. Accompanied by photon emission, the momentum of the atoms by further losses. In short, atomic cooling is the process by which the laser’s constant action on the atom slows it down. In the same way, we can achieve precise control of the velocity and direction of atoms through precise control of the laser, so-called the atomic manipulation [38].
In order to obtain a better signal-to-noise ratio, we should prepare the atom cloud to reach the magnitude of 10 8 to 10 9 atoms. As shown in Figure 1, 87 R b vapor becomes an atomic cloud through the action of 2D-MOT and 3D-MOT. After the capture zone, the number of atom clouds is continuously reduced due to factors such as the diffusion of cold atom clouds and the collision of various cavity walls. Furthermore, when there is background atom gas in the vapor chamber, the probe light will emit a large amount of stray fluorescence, which will increase the background noise and decrease the signal-to-noise ratio.
In short, considering the importance of atom cloud preparation and its requirements, the atom cloud preparation process with high reliability is an important guarantee for achieving high precision of the atomic gyro. This is the reason why we use more time to get more atoms into the interference area, eventually led to a lower data rate of CAIG.

2.2. Atomic State-Selection Progress

In the process of atomic interference, the interference is usually selected to the D2 line of 87 R b , so that the 87 R b atom enters the specific cycle between the ground state and the excited state, and gives the rotation information according to the population of atoms in each energy state. To avoid the effects of stray atoms and to improve the signal-to-noise ratio, the progress of atomic state-selection is necessary [39].
Depending on the complex energy state structure of 87 R b , atoms at each state can only absorb the specific frequency and then transit to a excited state. The selected interference progress is between 5 2 S 1 / 2 , F = 1 5 2 S 1 / 2 , F = 2 . Before the interference progress, the pump laser is used to transfer the available atoms of other energy levels 5 2 S 1 / 2 , F = 1 , as is shown in Figure 2.
In addition, there are other m F 0 atoms that will reduce the signal-to-noise ratio of the interference result due to changes in their sensitive magnetic fields. They are generally blown away by repulsion laser, and finally a pure atom cloud is obtained and trapped in a three-dimensional magneto-optical trap, and then they are ready to enter the interference region.

2.3. Atomic Interference Progress

In this section, a three-pulse interferometer is taken as an example to introduce the interference progress of CAIG. When the action laser is detuned to δ 12 δ A C = 0 , the most simplified transition probability with the Raman frequency Ω r = Ω e f f can be obtained by
P g e = sin 2 Ω e f f t 2 = 1 2 ( 1 cos Ω e f f t )
It can be seen that the transition probability fluctuates with the laser duration t. When the duration of the pulse is t = π / 2 Ω r , then half of the atoms will undergo state transition, and this pulse is called ‘ π / 2 -pulse’. The ‘ π / 2 -pulse’ is used as a splitter mirror in the optical interferometry. On the other hand, When the duration of the pulse is t = π / Ω r , then all of the atoms will undergo state transition, and this pulse is called ‘ π -pulse’. The ‘ π -pulse’ is used as a reflecting mirror in the optical interferometry. In the atomic interference region, three pairs of Raman laser are used as a splitter mirror, reflect mirror, and splitter mirror in the optical interferometry. In summary, we can control the atom cloud by controlling the Raman pulses to realize the atomic interferometer shown in Figure 1.
Interferometric gyroscopes [40], including optical interference gyros and cold atomic interference gyros, are based on the Sagnac effect, as shown in (2). The core issue of interferometric gyroscopes is to measure and separate the phase shift caused by the rotation.
Φ S a g n a c = 4 π E h c 2 A · Ω
As shown in Figure 3, the phase shift caused by the rotation caused by the Doppler effect is
Φ D o p p l e r = ( ϕ 1 ϕ 2 A ) ( ϕ 2 B ϕ 3 )
where ϕ i is the phase from the ith pulse, ϕ i = k e f f ( i ) · z ( t ) + ϕ a r b ( i ) , related to the wave vector k e f f , the atomic position z ( t ) , and the random phase of the pulse ϕ a r b . Due to the Doppler effect, a slight deflection θ of the pulse occurs when the pulse acts on the atom, θ = ± Ω T 1 . When we ignore random phases, there are
ϕ 1 k e f f · θ 1 L = k e f f · Ω T · L ϕ 2 A = ϕ 2 B ϕ 3 k e f f · θ 3 · L = k e f f · Ω T · L
Substituting (4) into (3) results in a phase shift Φ D o p p l e r ,
Φ D o p p l e r = ( ϕ 1 ϕ 2 A ) ( ϕ 2 B ϕ 3 ) = 2 k e f f · Ω · T L = 2 k e f f · Ω T · v · T
When the actual rotation axis is not perpendicular to the interference plane, (5) is corrected into (6). Defining equivalent interference area A = ( k e f f / M ) T × v T = ( k e f f / M ) T × L , (5) is converted into Φ D o p p l e r = 2 M A · Ω / . Therefore, the Doppler effect phase shift caused by rotation is equal to the Sagnac phase expression of (2), where is the reduced Planck constant. On the other hand, it can be seen that the interference area A is proportional to the square of the pulse interval time T, that is how the time domain atomic gyro increases its sensitivity.
Φ D o p p l e r = 2 k e f f · ( Ω × v ) T 2
Based on the above principle of rotation measurement, the output phase of the atomic gyro in the laboratory consists of several parts, as shown in (7).
Δ Φ t o t a l = k e f f · g · T 2 + 2 M A Ω 2 k e f f · ( Ω × g ) T 3 + Δ Φ 0
where k e f f · g · T 2 is the phase part caused by gravity, 2 M A Ω / is the phase part caused by rotation, 2 k e f f · ( Ω × g ) T 3 is the phase correction when gravity g is not parallel to the actual rotation axis, and Φ a r b is the Random phase introduced by three pulses, Φ a r b = ϕ a r b ( 1 ) 2 ϕ a r b ( 2 ) + ϕ a r b ( 3 ) .
Typical dual interference loops with opposite-projectile atom clouds are shown in Figure 4, and we can get a phase shift from each loop. As shown in (8), we can process the two-loop data output to obtain the common mode part and the differential mode part.
Φ t o t 1 = k e f f · g · T 2 + 2 m Ω · A 1 2 k e f f · ( Ω × g ) T 3 + Δ Φ 0 Φ t o t 2 = k e f f · g · T 2 + 2 m Ω · A 2 2 k e f f · ( Ω × g ) T 3 + Δ Φ 0
In two interference loops, the equivalent loop area vector A 1 = A 2 . Then the phase shift related to the rotation is successfully extracted by
Φ t o t 1 Φ t o t 2 = 2 m Ω · ( A 1 A 2 ) = 4 m Ω · A 1

2.4. Fluorescence Detection Progress

After the interference process, the population of the atomic state is obtained by fluorescence detection progress [41]. The ratio of atoms in two ground state 5 2 S 1 / 2 , F = 1 and 5 2 S 1 / 2 , F = 2 is calculated and converted into the transition probability result of the interference process.
P a b = | C b ( t ) | 2 = 1 2 ( 1 cos Δ Φ )
The absorption imaging technique is adopted here, and the fluorescence multiplier is received by the photomultiplier tube, and the total number of atoms at different energy states is obtained by time domain integration, in another word, the ratio of the population is obtained. In the specific implementation process, it is necessary to detect the laser light intensity and frequency to obtain precise control, and also needs to isolate stray light isolation, external temperature, and electromagnetic field isolation treatment.
It is worth noting that during the probing process, we need a unique and determined 2 π period from the conversion of the ratio of the population P a b to the phase result Δ Φ . In order to achieve dynamic measurements, it is common to reduce the scale factor to obtain a larger range of interference phase result which is from larger angular velocity variations. In order to maintain the system accuracy of the cold atom interference gyroscope, only a small dynamic measurement range can be obtained, and it still cannot be independently used for the dynamic measurement.

3. The Novel Monitoring Navigation Method of Triaxial CAIG & FOG

The principled shortcomings of CAIG have fully explained above. In this section, a new scheme is proposed to improve the positioning accuracy of the conventional SINS through additional CAIG. Here we take FOG as an example, which is most commonly used in underwater vehicles.

3.1. The Monitoring Navigation Scheme of Triaxial CAIG & FOG

As is known that bias drift of FOG is the main error source affecting the accuracy of SINSs. In the monitoring navigation system, the triaxial cold-atom interference gyroscope is used as an additional gyro to monitor an IMU of FOGs, as is shown in Figure 5, the triaxial FOG and CAIG is fixed to the underwater vehicle with the same body-frame. It is worth noting that we expect the measurement axes of FOG and CAIG to be approximately parallel to ensure the linearity of the system. The difference between the triaxial rotation rate of the CAIG and the FOG is the measurement observations. Through the Kalman filtering, the triaxial bias of the FOG is calculated and corrected, so that the monitoring navigation system works in the closed loop state. Due to the limitation of CAIG in dynamic measurement, the monitoring system achieves the estimation of the gyro through CAIG in a relatively stable process. In the case of higher dynamics, the navigation output of position and attitude are performed by the fiber-optic SINS. In order to reduce the influence of data out-of-synchronization on the performance of the monitoring system, the rotation rate output of the cold atom interference gyroscope and the fiber-optic gyroscope needs to be synchronized by data processing to realize the estimation of the gyro drift. Finally, the positioning accuracy of the long-term navigation based on the fiber-optic SINS can be improved.
In practical engineering, an inevitable misalignment angle between the cold atom interference gyroscope frame (represented by A-frame) and the fiber-optic gyroscope frame (represented by F-frame) exists. The existence of the misalignment angle will directly affect the estimation of the triaxial gyro drift. Considering the constant misalignment angle, the direction cosine matrix C A F between A-frame and F-frame
C A F = I [ φ a × ]
where, [ φ a × ] is the anti-symmetric matrix of misalignment angle φ a = [ φ a x , φ a y , φ a z ] T .
The rotation rate ω i b A b A measured by CAIG is transformed into the F-frame from A-frame,
ω i b F b F = C A F ω i b A b A
As is known, the measured values of gyros are
ω ^ i b F b F = ω i b F b F + ε F ω ^ i b A b A = ω i b A b A + ε A
where, ω ^ i b F b F and ω ^ i b A b A are the measured value of rotation rate of FOG and CAIG, ω i b F b F and ω i b A b A are the actual value of rotation rate in F-frame and in A-frame. ε F is the bias of FOG, ε F = [ ε F x , ε F y , ε F z ] T , and ε A is the bias of CAIG, ε A = [ ε A x , ε A y , ε A z ] T .
In the triaxial CAIG & FOG monitoring navigation system, the misalignment angle φ a , and the bias of FOG ε F is constant, then the system state equation is established as (14),
X k + 1 = Φ k + 1 , k X k + Γ k W k
where, the state variables X = [ φ a x , φ a y , φ a z , ε F x , ε F y , ε F z ] T .
And the observation is the rotation rate difference δ ω between the CAIG & FOG,
δ ω = ω ^ i b F b F ω ^ i b A b A = ε F ε A + ω i b A b A [ φ a × ]
where, the bias of high-precision cold atom interference gyroscope ε A is much smaller than ε F , and can be ignored here. Then (15) can be simplified to
δ ω = ω i b A b A [ φ a × ] + ε F
Organize (16) into a matrix form, the observation equation is established as
Z k = H k X k + V k
where V k is a zero mean valued Gaussian white noises with variance of R k , and the observing matrix is
H k = 0 ω i b A z b A ω i b A y b A 1 0 0 ω i b A z b A 0 ω i b A x b A 0 1 0 ω i b A y b A ω i b A x b A 0 0 0 1

3.2. The Observability Analysis of Monitoring Navigation System

We gave the model of the monitoring system in the previous section, and then the next step is to perform an observability analysis. One of the reasons for observability analysis is the need to determine the efficiency of the Kalman filter designed to estimate the state of the system. On the other hand, if the observability of the system is poor, we cannot obtain an accurate estimate of the state even if the noise level is negligible. In other words, the observability analysis determines the conditions under which the system can obtain an accurate estimate of the state.
In observability analysis theory, the system of interest is the homogeneous system. Consider the observation equation of a continuous system
z j ( t ) = H j x ( t )
where, z j ( t ) is the observation in j-th period, and H j is the observing matrix.
Differentiate the measurement vector n-1 times inside each segment to obtain
z j ( t ) = H j x ( t ) z ˙ j ( t ) = H j x ˙ ( t ) = H j Φ j x ( t ) z ¨ j ( t ) = H j x ¨ ( t ) = H j Φ j 2 x ( t ) z j ( n 1 ) ( t ) = H j x ( n 1 ) ( t ) = H j Φ j n 1 x ( t )
which can be written as
Z j ( t ) = Q j X ( t )
where, Q j is the observability matrix in the j-th period,
Q j T = H j T H j Φ j T H j Φ j 2 T H j Φ j n 1 T
According to the principle of piece-wise constant systems (PWCS) observability analysis theory [42], when the rank of the Total Observability Matrix Q r is equal to the number of observable state variables, the system is fully observable. Furthermore, as the system model shown in Section 3.1, the observability matrix for the j-th period can be written as
Q r = Q 1 Q 2 Φ n 1 Q r 1 r 1 Φ n 1 = Q 1 Q 2 Q r
In our monitoring system, there is
Q r T = H j T , H j Φ T , , H j Φ n 1 T = H j T , H j T , , H j T
H j = 0 ω i b A z b A ( j ) ω i b A y b A ( j ) 1 0 0 ω i b A z b A ( j ) 0 ω i b A x b A ( j ) 0 1 0 ω i b A y b A ( j ) ω i b A x b A ( j ) 0 0 0 1
where, ω i b A x b A ( j ) , ω i b A y b A ( j ) , ω i b A z b A ( j ) corresponds to the triaxial measurements of CAIG.
To simplify the analysis, the respective stripped observability matrix (SOM) U s , 0 ( j ) is
U s , 0 ( j ) = 0 ω i b A z b A ( i ) ω i b A y b A ( i ) 1 0 0 ω i b A z b A ( i ) 0 ω i b A x b A ( i ) 0 1 0 ω i b A y b A ( i ) ω i b A x b A ( i ) 0 0 0 1 0 ω i b A z b A ( j ) ω i b A y b A ( j ) 1 0 0 ω i b A z b A ( j ) 0 ω i b A x b A ( j ) 0 1 0 ω i b A y b A ( j ) ω i b A x b A ( j ) 0 0 0 1
Therefore, if there are any two time periods i, j, (27) is established, r a n k [ U s , 0 ( j ) ] = 6 , and the system’s six state variables are completely observable.
ω i b A x b A ( i ) ω i b A x b A ( j ) 0 ω i b A y b A ( i ) ω i b A y b A ( j ) 0 ω i b A z b A ( i ) ω i b A z b A ( j ) 0
Based on the above observability analysis, common movement mode of underwater vehicles, such as static mode, triaxial sway mode, constant speed mode, and constant acceleration mode is shown here,
  • In the static mode and constant speed mode, r a n k [ U s , 0 ( j ) ] = 3 , which means the system isn’t completely observable. It is worth noting that in this monitoring navigation system, the accurate estimation of the triaxial bias and the misalignment angle is indispensable, otherwise, the monitoring system loses its accuracy.
  • In constant acceleration mode of the underwater vehicle, the latitude and the height change and affect the measurement of earth’s rotation, but these changes will be submerged in the gyro noise and is difficult to be measured.
  • In the case of triaxial sway mode, (27) can be realized in several swing periods. The bias and misalignment angle are completely observable. In fact, underwater surge sway is a necessary condition for our monitoring system to achieve identification. On the other hand, when in land vehicles or other carriers, there may be bumps or other motions similar like swaying, the experiment results will be shown in the field test in Section 5.

4. Simulation Experiment of Monitoring Navigation System

The simulation verification for the monitoring navigation system described in the previous section is shown here. The bias of the cold atom interference gyroscope is constant and negligible. The bias of the FOG is 0.1 / h , and the triaxial misalignment angles between A-frame and F-frame are 1 , 2 , 3 , the data-rate of CAIG is set to 5 Hz , and the FOG‘s data-rate is 100 Hz . The initial value of the Kalman filtering is selected as follows,
X 0 = 0 , 0 , 0 , 0 , 0 , 0 T P 0 = d i a g ( 1 ) 2 , ( 2 ) 2 , ( 3 ) 2 , ( 0.1 / h ) 2 , ( 0.1 / h ) 2 , ( 0.1 / h ) 2
The simulation tests were carried out in Matlab on a computer with an Intel(R) Core(TM) i7-8500U CPU 3.79 GHz and 20 GB of RAM.

4.1. Constant Speed Mode

In order to verify the monitoring effect under constant speed motion, the initial heading angle is set to be 30 , both pitch angle and roll angle is set to be 0 , and the underwater vehicle is at the constant speed of 15 knot . As a result of the simulation, the misalignment angle and the gyro bias are shown in Figure 6.
It can be seen that the simulation results are consistent with the conclusions of the observability analysis. The movement of underwater vehicles at the speed of 15 knot cannot achieve effective excitation, and the monitoring system can not estimate the bias of FOG or misalignment angle in a short time.

4.2. Constant Accelerate Mode

In order to verify the monitoring effect under constant accelerate motion, the initial heading angle is set to be 30 , both pitch angle and roll angle is set to be 0 , the initial speed is set to be 0 m / s , and the acceleration is 0.05 m 2 / s . As a result of the simulation, the misalignment angle and the gyro bias are shown in Figure 7.
It can be seen that the stable bias estimation cannot be achieved within ten minutes, which is consistent with the previous analysis.

4.3. Triaxial Swing Mode

In order to verify the monitoring effect under triaxial swing mode, the initial heading angle, pitch angle, and roll angle are set to be 0 , the initial speed is set to be 0 m / s , and the acceleration is set to be 0 m 2 / s . As it‘s known that the sea condition is divided into 10 grades referred to Douglas sea scale [43,44], several sway parameters of our ship under typical sea conditions are represented in Table 1. In the sea condition of level-2, the wave is small, the wavelength is short, but the waveform is significant. In the sea condition of level-4, the waves have a distinct shape and form white waves everywhere. In the sea condition of level-6, the wind-cut waves on the crests begin to stretch along the wavy slope into bands. In typical sea conditions, the monitoring navigation system based on CAIG/FOG has been verified. Table 2 is the simulation result and the filter convergence time. At the same time, Figure 8 shows that the misalignment angle can be recognized extremely quickly, and even worse sea condition with more obvious three-axis swing can reduce the convergence time.
The above simulation result shows that the misalignment angle between two frames and the FOG‘s bias are not well estimated in the constant speed mode or in the constant acceleration mode. On the other hand, the triaxial swing mode such as ships in typical sea conditions, the estimation algorithm for the misalignment angle between two frames and the FOG‘s bias is with high precision, especially the convergence time of the estimated misalignment angle is extremely short.
Due to the excellent performance of the monitoring navigation system in the case of three-axis sway, we demonstrate the impact of CAIG‘s data rate in Figure 9. It can be seen that the increase in data rate directly corresponds to the decrease of convergence time of the monitoring system, but has no effect on the estimation result. However, in the actual carrier experiment, the difference between the data rates of two gyros may introduce a difference in time label, and then introduce an obvious error into the monitoring navigation system.

5. The Field Test of the Monitoring Navigation Method

In this section, we conducted a land vehicle experiment to initially validate the monitoring navigation algorithm. In the field test, there are bumps, direction-shifting, and vehicle steering during the vehicle’s movement. It is expected that gyro’s bias and misalignment angle estimation can be achieved. In 4 May 2018, the vehicle experiment was carried out in the Science Park of Harbin Institute of Technology, as shown in Figure 10.
Figure 11 shows our experimental system, consisting of two INSs, an uninterruptible power supply (UPS), a stabilized voltage supply (SVS), and a GPS receiver. The GPS signal collected here is used to obtain the real position, and to obtain more real speed information as a reference. The results of the following inertial navigation algorithm are not involved in GPS. Limited by the test conditions, two sets of FOGs with different precisions are used for the field experiment. The reference information in the experiment is from a high-precision fiber gyro of the order of 10 /h with a data-rate of 100 Hz. The fiber-optic gyroscope being monitored is on the order of 10 /h with a data rate of 100 Hz, both gyros are developed by our laboratory. During the experiment, we circled the park six times in one hour, the travel distance was more than 20 km, and the maximum speed of the vehicle exceeded 15 m/s. The monitoring result is shown in Figure 12, it is worth noting that during the experiment, the estimated effect is limited by the dynamic measurement errors of two gyros during the Turns and bumps progress. Referring to the accuracy of the experimental gyro, the estimated value is considered to be stable when the fluctuation of the estimated bias is less than 0.01 /h.
In Figure 13a, the blue dotted line is the speed error by the normal inertial navigation algorithm without the help of GPS signal, the red solid line is the speed error under the monitoring algorithm. In Figure 13b, the blue dotted line is the positioning error of the original inertial navigation algorithm; the red solid line is the positioning error under the monitoring algorithm. Due to the limitation of the accuracy of the FOG ( 10 1 / h ) in this part, the maximum error of the original inertial navigation system is in line with theoretical expectations. As can be seen from Figure 13, the monitoring system achieves a better speed output and a better positioning result.
Furthermore, we conducted an experiment in August 2018 in the sea near Zhanjiang, China, and the ship experiment route is shown in Figure 14, while the experimental device is shown in Figure 15, During the experiment, the sea state is about level-2. The reference information in the experiment is from a high-precision gyro of the order of 10 3 / h with a data-rate of 100 Hz . The fiber-optic gyroscope being monitored is on the order of 10 1 / h with a data rate of 100 Hz. The monitoring result is shown in Figure 16. It can be seen that a stable estimate of the misalignment angle and the triaxial gyro bias can be achieved at 2 min.
In Figure 17a, the blue dotted line is the speed error by the original inertial navigation algorithm, the red solid line is the speed error under the monitoring algorithm. In Figure 17b, the blue dotted line is the positioning error of the original inertial navigation algorithm; the red solid line is the positioning error under the monitoring algorithm. Due to the limitation of the accuracy of the FOG ( 10 1 / h ) in this part, the maximum error of the original inertial navigation system is in line with theoretical expectations. As can be seen from Figure 17, the monitoring system achieves a better speed output result and a significantly improved positioning result.

6. Conclusions

The monitoring navigation system proposed in this paper can achieve an approximate stable estimation of gyro’s bias drift under general dynamic conditions, and compensated to the SINS to obtain a better positioning result. However, due to the defects of CAIG’s dynamic measurement range, it is still necessary to use FOG for navigation in a highly dynamic environment. With the further development of CAIG, the increase in data rate can further accelerate the monitoring estimation progress. Furthermore, a closer combination of FOG and CAIG can be realized, and a rough range of phases can be exactly realized by FOG to achieve high precision under high dynamics measuring.
Considering the current research status of cold atom interference gyroscopes, there are still some limitations in the sea test under this monitoring model. For example, the sampling rate will limit the timing and estimation effect of the position calibration after the monitoring process, and the data delay between the two gyros will directly affect the estimation error. Another important point is that the dynamic performance of the cold atom interference gyroscope is insufficient. Due to CAIG‘s accuracy is proportional to the scale factor, the higher measurement accuracy means a greater possibility that the phase results will appear on other interference stripes, and greatly affecting the measurement performance.

Author Contributions

Conceptualization, Q.L.; Data curation, Q.L.; Formal analysis, L.Z.; Methodology, L.Z.; Project administration, W.G. and R.L.; Supervision, W.G.; Validation, R.L., Z.Y. and S.L.

Funding

This research was funded by National Natural Science Foundation of China (Grant No. 51509063, 61503090, 51709068, 61573117,11674362), Natural Science Foundation of Heilongjiang Province (Grant No. F2015031), the China Postdoctoral Science Foundation (No. 2018M631920), the Fundamental Research Funds for the Central Universities (Grant No. HIT.MKSTISP.2016 03), Postdoctoral Foundation of Heilongjiang Province Government (No. LBH-Z17094, LBH-Z17091), the Outstanding Youth Foundation of Hubei Province of China under Grant No. 2018CFA082, the National Defense Science and Technology Innovation Project of China and the Youth Innovation Promotion Association of Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CAIGCold Atom Interference Gyroscope
FOGFiber-optic Gyroscope
INSInertial Navigation System
SINSStrapdown Inertial Navigation System
IMUInertial Measurement Unit
PWCSPiece-wise Constant Systems

References

  1. Xia, X.W.; Sun, Q. Initial Alignment Algorithm Based on the DMCS Method in Single-Axis RSINS with Large Azimuth Misalignment Angles for Submarines. Sensors 2018, 18, 2123. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, Y.; Yu, F.; Wang, Y.; Wang, K. A Robust SINS/VO Integrated Navigation Algorithm Based on RHCKF for Unmanned Ground Vehicles. IEEE Access 2018, 6, 56828–56838. [Google Scholar] [CrossRef]
  3. Zhang, Y.; Yu, F.; Gao, W.; Wang, Y. An Improved Strapdown Inertial Navigation System Initial Alignment Algorithm for Unmanned Vehicles. Sensors 2018, 18, 3297. [Google Scholar] [CrossRef] [PubMed]
  4. Sun, Q.; Tian, Y.; Diao, M. Cooperative Localization Algorithm based on Hybrid Topology Architecture for Multiple Mobile Robot System. IEEE Internet Things J. 2018. [Google Scholar] [CrossRef]
  5. Hashimoto, M.; Cabuz, C.; Minami, K.; Esashi, M. Silicon resonant angular rate sensor using electromagnetic excitation and capacitive detection. J. Micromech. Microeng. 1995, 5, 219. [Google Scholar] [CrossRef]
  6. Söderkvist, J. Micromachined gyroscopes. Sens. Actuators A Phys. 1994, 43, 65–71. [Google Scholar] [CrossRef]
  7. Lawrence, A. Modern Inertial Technology: Navigation, Guidance, and Control; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  8. Lefevre, H.C. The Fiber-Optic Gyroscope; Artech House: Norwood, MA, USA, 2014. [Google Scholar]
  9. Sun, Q.; Diao, M.; Zhang, Y.; Li, Y. Cooperative Localization Algorithm for Multiple Mobile Robot System in Indoor Environment Based on Variance Component Estimation. Symmetry 2017, 9, 94. [Google Scholar] [CrossRef]
  10. Gustavson, T.; Bouyer, P.; Kasevich, M. Precision rotation measurements with an atom interferometer gyroscope. Phys. Rev. Lett. 1997, 78, 2046. [Google Scholar] [CrossRef]
  11. Kasevich, M.; Chu, S. Atomic interferometry using stimulated Raman transitions. Phys. Rev. Lett. 1991, 67, 181. [Google Scholar] [CrossRef]
  12. Gerginov, V.; Nemitz, N.; Weyers, S.; Schröder, R.; Griebsch, D.; Wynands, R. Uncertainty evaluation of the caesium fountain clock PTB-CSF2. Metrologia 2009, 47, 65. [Google Scholar] [CrossRef]
  13. Levi, F.; Calonico, D.; Calosso, C.E.; Godone, A.; Micalizio, S.; Costanzo, G.A. Accuracy evaluation of ITCsF2: A nitrogen cooled caesium fountain. Metrologia 2014, 51, 270. [Google Scholar] [CrossRef]
  14. Heavner, T.P.; Donley, E.A.; Levi, F.; Costanzo, G.; Parker, T.E.; Shirley, J.H.; Ashby, N.; Barlow, S.; Jefferts, S. First accuracy evaluation of NIST-F2. Metrologia 2014, 51, 174. [Google Scholar] [CrossRef]
  15. Fang, F.; Chen, W.; Liu, K.; Liu, N.; Suo, R.; Li, T. Design of the new NIM6 fountain with collecting atoms from a 3D MOT loading optical molasses. In Proceedings of the 2015 Joint Conference of the IEEE International Frequency Control Symposium & the European Frequency and Time Forum (FCS), Denver, CO, USA, 12–16 April 2015; pp. 492–494. [Google Scholar]
  16. Fang, F.; Li, M.; Lin, P.; Chen, W.; Liu, N.; Lin, Y.; Wang, P.; Liu, K.; Suo, R.; Li, T. NIM5 Cs fountain clock and its evaluation. Metrologia 2015, 52, 454. [Google Scholar] [CrossRef]
  17. Kasevich, M.; Chu, S. Measurement of the gravitational acceleration of an atom with a light-pulse atom interferometer. Appl. Phys. B 1992, 54, 321–332. [Google Scholar] [CrossRef]
  18. Malossi, N.; Bodart, Q.; Merlet, S.; Lévèque, T.; Landragin, A.; Dos Santos, F.P. Double diffraction in an atomic gravimeter. Phys. Rev. A 2010, 81, 013617. [Google Scholar] [CrossRef]
  19. Merlet, S. Détermination aBsolue De G Dans Le Cadre De L’expérience De La Balance Du Watt. Ph.D. Thesis, Observatoire de Paris, Paris, France, 2010. [Google Scholar]
  20. Merlet, S.; Bodart, Q.; Malossi, N.; Landragin, A.; Dos Santos, F.P.; Gitlein, O.; Timmen, L. Comparison between two mobile absolute gravimeters: Optical versus atomic interferometers. Metrologia 2010, 47, L9. [Google Scholar] [CrossRef]
  21. Schmidt, M.; Senger, A.; Hauth, M.; Freier, C.; Schkolnik, V.; Peters, A. A mobile high-precision absolute gravimeter based on atom interferometry. Gyroscopy Navig. 2011, 2, 170. [Google Scholar] [CrossRef]
  22. Fixler, J.B.; Foster, G.; McGuirk, J.; Kasevich, M. Atom interferometer measurement of the Newtonian constant of gravity. Science 2007, 315, 74–77. [Google Scholar] [CrossRef]
  23. Cadoret, M.; De Mirandes, E.; Cladé, P.; Guellati-Khélifa, S.; Schwob, C.; Nez, F.; Julien, L.; Biraben, F. Combination of Bloch oscillations with a Ramsey-Bordé interferometer: New determination of the fine structure constant. Phys. Rev. Lett. 2008, 101, 230801. [Google Scholar] [CrossRef]
  24. Müller, H.; Peters, A.; Chu, S. A precision measurement of the gravitational redshift by the interference of matter waves. Nature 2010, 463, 926. [Google Scholar] [CrossRef]
  25. Rosi, G.; Sorrentino, F.; Cacciapuoti, L.; Prevedelli, M.; Tino, G. Precision measurement of the Newtonian gravitational constant using cold atoms. Nature 2014, 510, 518. [Google Scholar] [CrossRef] [PubMed]
  26. Jin, W. Precision measurement with atom interferometry. Chin. Phys. B 2015, 24, 053702. [Google Scholar]
  27. Gustavson, T.; Landragin, A.; Kasevich, M. Rotation sensing with a dual atom-interferometer Sagnac gyroscope. Class. Quantum Gravity 2000, 17, 2385. [Google Scholar] [CrossRef]
  28. Stockton, J.; Takase, K.; Kasevich, M. Absolute geodetic rotation measurement using atom interferometry. Phys. Rev. Lett. 2011, 107, 133001. [Google Scholar] [CrossRef] [PubMed]
  29. Riou, I.; Mielec, N.; Lefèvre, G.; Prevedelli, M.; Landragin, A.; Bouyer, P.; Bertoldi, A.; Geiger, R.; Canuel, B. A marginally stable optical resonator for enhanced atom interferometry. J. Phys. B Atomic Mol. Opt. Phys. 2017, 50, 155002. [Google Scholar] [CrossRef] [Green Version]
  30. Yao, Z.W.; Lu, S.B.; Li, R.B.; Luo, J.; Wang, J.; Zhan, M.S. Calibration of atomic trajectories in a large-area dual-atom-interferometer gyroscope. Phys. Rev. A 2018, 97, 013620. [Google Scholar] [CrossRef] [Green Version]
  31. Battelier, B.; Barrett, B.; Fouché, L.; Chichet, L.; Antoni-Micollier, L.; Porte, H.; Napolitano, F.; Lautier, J.; Landragin, A.; Bouyer, P. Development of compact cold-atom sensors for inertial navigation. Quantum Opt. Int. Soc. Opt. Photonics 2016, 9900, 990004. [Google Scholar] [Green Version]
  32. Bochkati, M.; Schön, S.; Schlippert, D.; Schubert, C.; Rasel, E. Could cold atom interferometry sensors be the future inertial sensors?—First simulation results. In Proceedings of the Inertial Sensors and Systems (ISS), Karlsruhe, Germany, 19–20 September 2017; pp. 1–20. [Google Scholar]
  33. Li, Y.; Niu, X.; Zhang, Q.; Zhang, H.; Shi, C. An in situ hand calibration method using a pseudo-observation scheme for low-end inertial measurement units. Meas. Sci. Technol. 2012, 23, 105104. [Google Scholar] [CrossRef]
  34. Lacambre, J.B.; Narozny, M.; Louge, J.M. Limitations of the unscented Kalman filter for the attitude determination on an inertial navigation system. In Proceedings of the 2013 IEEE Digital Signal Processing and Signal Processing Education Meeting (DSP/SPE), Napa, CA, USA, 11–14 August 2013; pp. 187–192. [Google Scholar]
  35. Lawall, J.; Bardou, F.; Bouchaud, J.P.; Saubamea, B.; Bigelow, N.; Leduc, M.; Aspect, A.; Cohen-Tannoudji, C. Recent Advances in Subrecoil Laser Cooling. AIP Conf. Proc. 1994, 323, 193–210. [Google Scholar]
  36. Allred, J.; Lyman, R.; Kornack, T.; Romalis, M. High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation. Phys. Rev. Lett. 2002, 89, 130801. [Google Scholar] [CrossRef]
  37. Neuman, K.C.; Block, S.M. Optical trapping. Rev. Sci. Instrum. 2004, 75, 2787–2809. [Google Scholar] [CrossRef] [PubMed]
  38. Ashkin, A. Optical trapping and manipulation of neutral particles using lasers. Proc. Natl. Acad. Sci. USA 1997, 94, 4853–4860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Gillot, J.; Gauguet, A.; Büchner, M.; Vigué, J. Optical pumping of a lithium atomic beam for atom interferometry. Eur. Phys. J. D 2013, 67, 263. [Google Scholar] [CrossRef]
  40. Bouyer, P. The centenary of Sagnac effect and its applications: From electromagnetic to matter waves. Gyroscopy Navig. 2014, 5, 20–26. [Google Scholar] [CrossRef]
  41. Rocco, E.; Palmer, R.; Valenzuela, T.; Boyer, V.; Freise, A.; Bongs, K. Fluorescence detection at the atom shot noise limit for atom interferometry. New J. Phys. 2014, 16, 093046. [Google Scholar] [CrossRef] [Green Version]
  42. Goshen-Meskin, D.; Bar-Itzhack, I. Observability analysis of piece-wise constant systems. I. Theory. IEEE Trans. Aerosp. Electron. Syst. 1992, 28, 1056–1067. [Google Scholar] [CrossRef]
  43. Salvesen, N.; Tuck, E.; Faltinsen, O. Ship motions and sea loads. Trans. SNAME 1970, 78, 250–287. [Google Scholar]
  44. Faltinsen, O. Sea Loads on Ships and Offshore Structures; Cambridge University Press: Cambridge, UK, 1993; Volume 1. [Google Scholar]
Figure 1. System implementation of cold atom interference gyroscope.
Figure 1. System implementation of cold atom interference gyroscope.
Sensors 19 00222 g001
Figure 2. Atomic state-selection Technique.
Figure 2. Atomic state-selection Technique.
Sensors 19 00222 g002
Figure 3. Schematic diagram of the Doppler effect.
Figure 3. Schematic diagram of the Doppler effect.
Sensors 19 00222 g003
Figure 4. Dual interference loops with opposite-projectile atom clouds.
Figure 4. Dual interference loops with opposite-projectile atom clouds.
Sensors 19 00222 g004
Figure 5. The monitoring navigation configuration scheme of triaxial cold atom interference gyroscope (CAIG) & fiber-optic gyroscope (FOG).
Figure 5. The monitoring navigation configuration scheme of triaxial cold atom interference gyroscope (CAIG) & fiber-optic gyroscope (FOG).
Sensors 19 00222 g005
Figure 6. The monitoring system simulation results under constant speed mode. The left figure is the trixial misalignment angle curve, and the right figure is the simulated curve for bias of FOG.
Figure 6. The monitoring system simulation results under constant speed mode. The left figure is the trixial misalignment angle curve, and the right figure is the simulated curve for bias of FOG.
Sensors 19 00222 g006
Figure 7. The monitoring system simulation results under constant accelerate mode. The left figure is the trixial misalignment angle curve, and the right figure is the simulated curve for bias of FOG.
Figure 7. The monitoring system simulation results under constant accelerate mode. The left figure is the trixial misalignment angle curve, and the right figure is the simulated curve for bias of FOG.
Sensors 19 00222 g007
Figure 8. The monitoring system simulation results under triaxial sway mode. The sea condition is level-2, level-4, and level-6, respectively corresponding to the red solid line, blue dashed line, and green dotted line. The left figure is the misalignment angle estimation curve, and the right figure is the simulation estimation curve of the gyro bias.
Figure 8. The monitoring system simulation results under triaxial sway mode. The sea condition is level-2, level-4, and level-6, respectively corresponding to the red solid line, blue dashed line, and green dotted line. The left figure is the misalignment angle estimation curve, and the right figure is the simulation estimation curve of the gyro bias.
Sensors 19 00222 g008
Figure 9. The monitoring system simulation results o uf different dander triaxiata ratel sway mode. The data rate of CAIG is 5 Hz, 20 Hz, and 100 Hz, respectively corresponding to the red solid line, blue dashed line, and green dotted line. The left figure is the misalignment angle estimation curve, and the right figure is the simulation estimation curve of the gyro bias.
Figure 9. The monitoring system simulation results o uf different dander triaxiata ratel sway mode. The data rate of CAIG is 5 Hz, 20 Hz, and 100 Hz, respectively corresponding to the red solid line, blue dashed line, and green dotted line. The left figure is the misalignment angle estimation curve, and the right figure is the simulation estimation curve of the gyro bias.
Sensors 19 00222 g009
Figure 10. The vehicle experiment route in Harbin Institute of Technology in 4 May 2018.
Figure 10. The vehicle experiment route in Harbin Institute of Technology in 4 May 2018.
Sensors 19 00222 g010
Figure 11. The vehicle experimental device.
Figure 11. The vehicle experimental device.
Sensors 19 00222 g011
Figure 12. The parameter estimation result of vehicle monitoring experiment. (a) The estimation of the misalignment angle after starting the monitoring system. (b) The stable estimation details of the misalignment angle. (c) The estimation of the triaxial bias after starting the monitoring system. (d) The stable estimation details of the triaxial bias.
Figure 12. The parameter estimation result of vehicle monitoring experiment. (a) The estimation of the misalignment angle after starting the monitoring system. (b) The stable estimation details of the misalignment angle. (c) The estimation of the triaxial bias after starting the monitoring system. (d) The stable estimation details of the triaxial bias.
Sensors 19 00222 g012
Figure 13. The output result of vehicle monitoring experiment. (a) Velocity errors in the east direction and in the north direction. (b) Positioning errors of INS.
Figure 13. The output result of vehicle monitoring experiment. (a) Velocity errors in the east direction and in the north direction. (b) Positioning errors of INS.
Sensors 19 00222 g013
Figure 14. The ship experiment route near Zhanjiang in August 2018.
Figure 14. The ship experiment route near Zhanjiang in August 2018.
Sensors 19 00222 g014
Figure 15. The experimental device of ship carrier.
Figure 15. The experimental device of ship carrier.
Sensors 19 00222 g015
Figure 16. The parameter estimation result of monitoring experiment on level-2 sea condition. (a) The estimation of the misalignment angle after starting the monitoring system. (b) The stable estimation details of the misalignment angle. (c) The estimation of the triaxial bias after starting the monitoring system. (d) The stable estimation details of the triaxial bias.
Figure 16. The parameter estimation result of monitoring experiment on level-2 sea condition. (a) The estimation of the misalignment angle after starting the monitoring system. (b) The stable estimation details of the misalignment angle. (c) The estimation of the triaxial bias after starting the monitoring system. (d) The stable estimation details of the triaxial bias.
Sensors 19 00222 g016
Figure 17. The positioning result of underwater vehicle monitoring experiment. (a) Velocity errors in the east direction and in the north direction. (b) Positioning errors of INS.
Figure 17. The positioning result of underwater vehicle monitoring experiment. (a) Velocity errors in the east direction and in the north direction. (b) Positioning errors of INS.
Sensors 19 00222 g017aSensors 19 00222 g017b
Table 1. The sway parameters of our ship under different sea conditions.
Table 1. The sway parameters of our ship under different sea conditions.
LevelWave Height (m)RollPitchHeading
AmplitudePeriodAmplitudePeriodAmplitudePeriod
level-20.1∼0.50.50 20.00 s 0.20 30.00 s 0.30 30.00 s
level-41.25∼2.54.75 21.00 s 0.64 10.00 s 0.70 18.61 s
level-64.0∼6.012.57 17.42 s 1.87 10.02 s 1.63 17.03 s
Table 2. The simulation result of monitoring system under traxial sway mode.
Table 2. The simulation result of monitoring system under traxial sway mode.
Sea ConditionsBias ( /h)Misalignment Angle ( )Convergence Time (s)
ε Fx ε Fy ε Fz ϕ ax ϕ ay ϕ az
level-20.09544480.09968460.10144491.0002.0003.00074.60
level-40.09963740.09940610.10054861.0002.0003.00064.40
level-60.09823970.09680950.09899471.0002.0003.00044.80

Share and Cite

MDPI and ACS Style

Zhang, L.; Gao, W.; Li, Q.; Li, R.; Yao, Z.; Lu, S. A Novel Monitoring Navigation Method for Cold Atom Interference Gyroscope. Sensors 2019, 19, 222. https://doi.org/10.3390/s19020222

AMA Style

Zhang L, Gao W, Li Q, Li R, Yao Z, Lu S. A Novel Monitoring Navigation Method for Cold Atom Interference Gyroscope. Sensors. 2019; 19(2):222. https://doi.org/10.3390/s19020222

Chicago/Turabian Style

Zhang, Lin, Wei Gao, Qian Li, Runbing Li, Zhanwei Yao, and Sibin Lu. 2019. "A Novel Monitoring Navigation Method for Cold Atom Interference Gyroscope" Sensors 19, no. 2: 222. https://doi.org/10.3390/s19020222

APA Style

Zhang, L., Gao, W., Li, Q., Li, R., Yao, Z., & Lu, S. (2019). A Novel Monitoring Navigation Method for Cold Atom Interference Gyroscope. Sensors, 19(2), 222. https://doi.org/10.3390/s19020222

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