It is difficult to analyze the nonstationary gyro signal in detail for the Allan variance (AV) analysis method. A novel approach in the time-frequency domain for gyro signal characteristics analysis is proposed based on the empirical mode decomposition and Allan variance (EMDAV). The output signal of gyro is decomposed by empirical mode decomposition (EMD) first, and then the decomposed signal is analyzed by AV algorithm. Consequently, the gyro noise characteristics are demonstrated in the time-frequency domain with a three-dimensional (3D) manner. Practical data of fiber optic gyro (FOG) and MEMS gyro are processed by the AV method and the EMDAV algorithm separately. The results indicate that the details of gyro signal characteristics in different frequency bands can be described with the help of EMDAV, and the analysis dimensions are extended compared with the common AV. The proposed EMDAV, as a complementary tool of the AV, which provides a theoretical reference for the gyro signal preprocessing, is a general approach for the analysis and evaluation of gyro performance.
Fundamental Research Funds for the Central UniversitiesNS2015037National Natural Science Foundation of China6153300861374115Jiangsu Key Laboratory of Internet of Things and Control Technologies1. Introduction
As a key device for angular velocity detection, gyro plays an important role in inertial navigation system. However, the errors caused by the influence of gyro working environment will be accumulated with working time, and the navigation accuracy will be deteriorated [1]. Consequently, it is important to analyze the characteristic of gyro signal, which is helpful to identify the performance of the gyros.
The traditional Allan variance (AV) algorithm, which can reflect the characteristics of signal, is widely used to evaluate the performance of gyro [2, 3]. However, it can only achieve the signal characteristics analysis in the time domain. Actually, due to the complexity of practical gyro signal, it is difficult to obtain the signal species and distributions comprehensively in time domain. So it is essential to analyze the gyro signal in the time-frequency domain.
Wavelet transform and empirical mode decomposition (EMD) are commonly adopted to implement the signal time-frequency analysis. Nonetheless, the problems of frequency overlapping, threshold instability, and decomposition and reconstruction algorithm errors existing in the wavelet transform bring some limitations for the gyro signal processing [4]. While being an adaptive local time-frequency analysis method [5], EMD does not possess the above disadvantages. The signal can be decomposed by its characteristic to a finite number of empirical mode functions by EMD, which is commonly applied in signal denoising [6], trend extraction [7], and so on.
In this paper, a gyro signal characteristic analysis approach based on the empirical mode decomposition and Allan variance (EMDAV) is proposed, which employs AV to analyze the decomposed gyro signal of EMD. Compared with the sole time domain characteristic of gyro signal obtained by traditional AV, the proposed method enriches the real gyro performance evaluation approaches, supplies more detailed gyro signal characteristics information in the joint time-frequency domain, and provides the theoretical basis for the denoising threshold determination in signal preprocessing.
2. Allan Variance Analysis
As a common analysis tool in the time domain, AV is widely applied to identify the gyro errors characteristics. Assuming that there is a gyro signal sample sequence of {xt} with the length of N and the sampling frequency of f, its discrete-time AV can be expressed as [8](1)σAVAR2τ=12x-tk+ττ-x-tkτ2,where the operator · denotes ensemble averaging, x-tkτ is the averaged measure of FOG signal at time tk, and τ=N/f is observation interval.
It is assumed that the various error sources are statistically independent; then the AV can be described as the sum of squares of quantization noise, angle random walk, bias instability, rate random walk, and rate ramp [9, 10]:(2)σ2τ=σQN2τ+σARW2τ+σBI2τ+σRRW2τ+σRR2τ=3Q2τ2+N2τ+2B2πln2+K2τ3+R2τ22,where σARW, σRRW, σBI, σQN, and σRR are the variances of the above five kinds of errors separately and Q, N, B, K, and R are the corresponding AV coefficients.
3. EMDAV Analysis
The traditional AV analysis method which presents the gyro noise characteristics in the only time domain cannot reflect its distributions among different frequency bands. However, the noise characteristics vary with the frequency. To reveal more details of the gyro noise characteristic, the EMDAV method is proposed in this paper.
As a method utilized to analyze the nonlinear and nonstationary signal, EMD decomposes the signal on the basis of time scale characteristics. There is no need to set the basic functions in advance for EMD which can decompose any type of signal. Thus EMD has obvious advantages in processing nonlinear and nonstationary signal.
The EMDAV algorithm proposed in the paper is based on the EMD and AV. After decomposing the gyro signal with EMD, the AV can be obtained with the help of intrinsic mode function (IMF). Lastly the AV logarithmic curve over time and frequency can be depicted. The detailed algorithm flow chart is shown in Figure 1.
Flow chart of AV algorithm based on EMD.
The steps of the EMDAV based on gyro signal characteristics analysis method can be summarized as follows:
(1) Find all the local maxima and minima of the original signal xt.
(2) Fit the estimated points of the first step using natural cubic spline to form an upper envelope of emaxt and a lower envelope of emint.
(3) Compute the mean of the two envelopes:(3)mit=emaxt+emint2,where i=1,2,…,n; n is the number of decomposition levels.
(4) Refine the estimate hit of the IMF by subtracting the mean solved in the previous step from the current signal:(4)hit=xit-mit,where(5)xit=xt,i=1,xt-∑j=1i-1hjt,i≥2.(5) Proceed from step (1) again unless a stopping criterion has been fulfilled. The stopping criterion is [11](6)SD=∑t=0Thi-1t-hit2hi-12t,where SD is a screening threshold. The iteration stops when SD is smaller than a constant. The constant value can be set between 0.2 and 0.3 [12].
Eventually, the IMF components of h1t,h2t,…,hnt and the remainder of rnt can be obtained. Then the original signal can be expressed as(7)xt=∑i=1nhit+rnt.(6) Calculate the AV coefficients corresponding to the IMF components and analyze the signal characteristics in different frequency bands.
From the above EMDAV procedures, we can see that the IMF of hit indicates the different frequency components. With the increase of decomposition level, the frequency of decomposed signal decreases gradually. Namely, the first IMF represents the highest frequency component of the original signal. The frequencies of subsequent IMFs reduce almost in the form of the negative power of 2 [13] until the frequency of the remainder is very small. Consequently, each IMF corresponds to the signal component of different frequency.
Compared with the traditional AV analysis method, the characteristics and distributions of gyro signal in different frequency bands can be depicted in detail with EMDAV algorithm. One more piece of dimensional information can be shown in three-dimensional (3D) logarithmic graph of EMDAV compared to that in 2D curve of traditional AV.
4. Experiment Results and Analysis
To verify the algorithm proposed in this paper, AV and EMDAV are compared to analyze the FOG and MEMS gyro signals. Practical fiber optical gyro (FOG) and microelectromechanical system (MEMS) gyro static signals are acquired with the sampling frequency of 100 Hz and the sampling duration of 600 s. The scene of experiment is shown in Figure 2.
FOG and MEMS gyros: (a) FOG on turntable and (b) MEMS gyro.
As we all know, traditional AV coefficients come from original gyro signal, and we can compare them with EMDAV coefficients at the same time. Figure 3 shows the original signals and parts of IMF components decomposed by EMD. It can be seen that the characteristics of FOG and MEMS gyro signals are similar: the first level of decomposed signal belongs to high frequency part, and the frequency decreases with the increasing of decomposition level.
IMFs from EMD of gyro signals: (a) EMD of FOG signal and (b) EMD of MEMS gyro signal.
The original signals and IMF components decomposed by EMD are analyzed with AV, and corresponding coefficients (AV coefficients and EMDAV coefficients) are plotted in a logarithmic graph as shown in Figure 4. From Figure 4 we can find that the 16 curves in the first group from top to bottom are the AV logarithm curves of MEMS gyro original signal (the 0 level, traditional AV) and IMFs (1–15 levels, EMDAV), respectively. Because of the poor accuracy of MEMS gyro, the biggest AV value reaches 105. The 16 profiles in the second group are the respective AV logarithm curves of FOG gyro original signal (the 0 level, traditional AV) and IMFs (EMDAV), where the biggest AV value reaches 1. Taking the AV values at correlation time of 0.01 s for comparison, the maximum is from the original signal for both gyros, by which the first level of decomposed signals followed. The AV values decrease until the decomposition level increases to the last one, whose value corresponds to the minimum.
Logarithmic diagram of AV and EMDAV.
In addition, it can be seen from Table 1 and Figure 4 that the signal distributions of FOG and MEMS gyro are similar; that is to say, with the increase of decomposition level, the proportion of angle random walk is reduced while the proportion of rate random walk is increased. Also it is obvious that angle random walk is dominant in the high frequency signal while rate random walk accounts for a large proportion in the low frequency signal. In consequence, gyro signal can be denoised specifically with the demand of the navigation system.
AV coefficients of gyro signals on some frequencies.
Signal
AV coefficients
N (∘/h1/2)
B (∘/h)
K (∘/h3/2)
R (∘/h2)
Q (μrad)
FOG
AV (original signal)
1.40e-3
5.71e-1
1.22e1
4.43e1
3.55e-2
EMDAV (signals decomposed by EMD)
Level 1
1.60e-3
4.64e-1
9.84e0
3.56e1
3.76e-2
Level 4
4.15e-4
1.75e-1
3.95e0
1.44e1
7.40e-3
Level 7
1.33e-4
6.66e-2
1.20e0
4.25e0
1.90e-3
Level 10
2.19e-5
1.02e-2
3.20e-1
1.24e0
3.17e-4
Level 15
4.32e-7
2.10e-4
2.36e-2
7.72e-2
6.15e-6
MEMS gyro
AV (original signal)
5.91e-1
3.10e1
6.20e2
2.24e3
1.39e0
EMDAV (signals decomposed by EMD)
Level 1
1.37e-1
1.41e1
2.10e2
7.15e2
8.78e0
Level 4
1.01e-1
9.85e1
2.16e3
7.85e3
2.37e0
Level 7
9.93e-2
4.92e1
8.50e2
2.99e3
1.40e0
Level 10
1.59e-2
7.36e0
2.54e2
9.81e2
2.30e-1
Level 15
1.30e-3
6.12e-1
6.54e1
1.80e2
1.81e-2
Figure 5 shows the relationship between AV and EMDAV coefficients, where the zero level means the AV coefficient of original gyro signals. As mentioned before, the frequency of each IMF is decreased in the form of negative power of 2 gradually, which indicates the fact that the higher the EMD decomposition levels the smaller the signal frequency. It also shows that the EMDAV coefficients trend to decrease with the decrease of signal frequency.
Relationship between AV coefficients (level 0) and EMDAV coefficients: (a) FOG and (b) MEMS gyro.
From the comparison of Figures 5(a) and 5(b), the coefficients of MEMS gyro are much larger than that of FOG; namely, the accuracy of MEMS gyro is worse than FOG. As for FOG (in Figure 5(a)), the EMDAV coefficients in the 5th level are larger, which means that the angle random walk noise, bias instability noise, rate random walk noise, and rate ramp noise are more dominant. The EMDAV coefficients of the 8th level and higher level are very small, and the EMD denoising method will be effective when the decomposition levels are larger than 7. As for MEMS gyro (in Figure 5(b)), the EMDAV coefficients are dominant in the 2nd level and are largest for 4 parameters, and reference decomposition level of 9 will be applicable for EMD denoising method. Therefore the algorithm proposed in this paper can provide guidance for gyro denoising.
The relationships among AV values of IMFs, correlation time, and EMD decomposition levels are shown in the AV and EMDAV 3D logarithmic graph of Figure 6, where the top curves are AV values of MEME gyro while the bottom are of FOG. Traditional AV logarithmic diagram is the 2D curve on the 0 level, and EMDAV logarithmic diagram is the 3D view on 1–15 levels. It can be found that though the AV values of MEMS gyro are larger than that of FOG, the signal magnitudes of the both gyros tend to reduce with the decrease of frequency. The 3D figures of both FOG and MEMS gyro signals show the special characteristics of the FOG and MEMS gyro performance. The AV values decrease with the increase of EMD levels and correlation time. The ridges, which also can be seen from Figure 4, are the collections of the points with a zero slope in logarithmic diagram of every IMF signal. Obviously the ridge of FOG signal is consistent with that of MEMS gyro signal.
AV and EMDAV 3D logarithmic graph: (a) front view and (b) back view.
The contrasts between AV and EMDAV are listed in Table 2. As can be seen from the table, EMDAV can analyze the both stationary and nonstationary signals with the analysis results showing the 3D relationship (served as the 3D fingerprint of a gyro) among AV, time, and frequency. In addition, various kinds of noises can be identified to provide more details about the signal characteristics of gyro. EMDAV can obtain more referable parameters than AV and provide a reference for signal denoising.
Performance contrasts between AV and EMDAV.
AV
EMDAV
Signal
Stationary signal
Stationary and nonstationary signal
Result presentation (fingerprint)
2D
3D
Domain of signal analysis
Time domain
Time-frequency domain
Ability of signal identification at different frequencies
No
Yes
Numbers of obtained parameters can be referenced
Less
More
Reference for signal denoising
No
Yes
5. Conclusion
An AV analysis algorithm based on EMD for gyro error is proposed in this paper. The AV value is calculated by the signal decomposed based on the characteristic time scale, and the signal characteristics of gyro error in different frequency bands are revealed. Theoretical analysis and experiments of practical FOG and MEMS gyro data indicate the following:
Compared with the characteristic of gyro error in the sole time domain given by the traditional AV, the EMDAV algorithm can describe the distributions on different kinds of gyro signals in time-frequency domain.
The EMDAV algorithm extends the dimensions of gyro signal characteristics analysis and hence obtains more parameters than traditional AV.
The EMDAV algorithm is able to determine the decomposition level for gyro signal denoising and provide a theoretical guidance to the gyro signal denoising and other preprocessing algorithms.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research has been mainly funded by the Fundamental Research Funds for the Central Universities (NS2015037) and partly funded by the National Natural Science Foundation of China (Grants nos. 61533008 and 61374115 to Jianye Liu) and Jiangsu Key Laboratory of Internet of Things and Control Technologies.
LaiJ.LvP.LiuJ.JiangB.Noncommutativity error analysis of strapdown inertial navigation system under the vibration in UAVs20129, article 13410.5772/527582-s2.0-84868016909PeesapatiR.SabatS. L.KarthikK. P.NayakJ.GiribabuN.Efficient hybrid Kalman filter for denoising fiber optic gyroscope signal2013124204549455610.1016/j.ijleo.2013.02.0132-s2.0-84883241502ChikovaniV. V.Performance parameters comparison of ring laser, coriolis vibratory and fiber-optic gyros based on Allan variance analysisProceedings of the IEEE 2nd International Conference on Actual Problems of Unmanned Air Vehicles Developments (APUAVD ’13)October 2013Kiev, UkraineIEEE15315610.1109/apuavd.2013.67053122-s2.0-84893663609QuC.-S.LuT.-Z.TanY.A modified empirical mode decomposition method with applications to signal de-noising2010361677310.3724/sp.j.1004.2010.000672-s2.0-77649195680GanY.SuiL.XiaoG.DuanY.EMD de-noising theory considering static and dynamic conditions and its applications in INS201324550751510.1007/978-3-642-37407-4-482-s2.0-84880735546ZhangY.WangS.XiaD.EMD-based denoising methods in the MEMS gyroscope de-driftProceedings of the 5th IEEE International Conference on Nano/Micro Engineered and Molecular Systems (NEMS ’10)January 2010Xiamen, China59159410.1109/nems.2010.55924712-s2.0-78649256457TianH.YuanR.ZhangC.JinJ.SongN.LinS.Nonlinear empirical mode predictive drift extraction on fiber optical gyroscope2011122541542110.1016/j.ijleo.2010.02.0242-s2.0-78751566761JainS.NandyS.ChakrabortyG.KumarC. S.RayR.ShomeS. N.Error modeling of various sensors for robotics application using allan variance techniqueProceedings of the IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC ’11)September 2011Xi’an, ChinaIEEE1410.1109/icspcc.2011.60617702-s2.0-81055149783ZhangN.LiX.Research on theoretical improvement of dynamic Allan variance and its application20113111707510.3788/aos201131.11060032-s2.0-82755188143LiJ. C.GaoF. Q.WangG. L.GaoW.ZhuW.WangM.Analysis of dynamic Allan variance for fiber optic gyro under vibration and variable temperature conditions20134091710.3788/cjl201340.09080042-s2.0-84885457181DangS.TianW.QianF.EMD- and LWT-based stochastic noise eliminating method for fiber optic gyro201144102190219310.1016/j.measurement.2011.07.0162-s2.0-84860403558HuangN. E.ShenZ.LongS. R.WuM. C.ShihH. H.ZhengQ.YenN.-C.TungC. C.LiuH. H.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis1998454197190399510.1098/rspa.1998.0193MR1631591GanY.SuiL. F.De-noising method for gyro signal based on EMD20114067457502-s2.0-84857412562