Journal of Sensors

Volume 2018, Article ID 2895187, 6 pages

https://doi.org/10.1155/2018/2895187

## Application of Improved Fast Dynamic Allan Variance for the Characterization of MEMS Gyroscope on UAV

School of Instrument Science and Opto-electronics Engineering, Beihang University, Beijing 100191, China

Correspondence should be addressed to Xueyun Wang; moc.liamg@aaub.nuyeuxgnaw

Received 28 January 2018; Accepted 8 April 2018; Published 8 May 2018

Academic Editor: Paolo Bruschi

Copyright © 2018 Qian Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Microelectromechanical systems (MEMS) are core components in unmanned aerial vehicles (UAV). The precision of MEMS sensors is very important when the GPS signal is invalid. However, the precision and performance of MEMS sensors will be degraded by the changing of environment. Therefore, estimation and identification of the various noise terms existing in MEMS sensors are deemed to be necessary. The Allan variance is a common and standard method to analyze MEMS sensors, but it cannot be used to analyze the dynamic characteristics. The dynamic Allan variance (DAVAR) is a sliding version of the Allan variance. It is a practical tool that could represent the nonstationary behavior of the MEMS signal. As the DAVAR needs to estimate the Allan variance at each time epoch, the computation time grows significantly with the length of the signal series. In this paper, in the case of MEMS gyroscope on UAV, an improved fast DAVAR algorithm based on the choice of relevant time is proposed to shorten the computation time. As an experimental validation, numerical experiments are conducted under the proposed method. The results demonstrate that the improved method could greatly increase the computation speed and does not affect the accuracy of estimation.

#### 1. Introduction

The unmanned aerial vehicles (UAV) have been widely used in civil and military applications, including search and rescue operations, area mapping, weather monitoring, and agricultural operations [1–4]. Whenever the inertial navigation system (INS) of UAV is concerned, cost or weight is always an issue; therefore, the accurate inertial sensors have been constantly excluded. Instead, the microelectromechanical systems (MEMS) have been universally used [5–7], which have the characteristics of lightweight, small mass, less expensive, and lower power requirements [8, 9]. Typically, MEMS sensors have large bias drifts and stochastic errors, which make it difficult to use the MEMS sensors as INS only. Generally, the combination of INS/GPS is used to provide an ideal navigation system with full capability of continuous position, velocity, and attitude output [10–12]. However, the accuracy of the integrated navigation system degrades with time when GPS signals are blocked in environments such as high buildings and indoors. In order to control the simple INS error within a certain range, it is necessary to estimate and identify the various noise terms existing in MEMS sensors.

Allan variance method is a time analysis technique developed by Dr. David Allan to study the characteristic of random noise terms and stability in precision oscillators used in clock application [13]. Allan variance method can be used to determine the characteristics of the underlying random processes which lead to data noises [14, 15], and it is also generally used to identify the errors of inertial sensors (i.e., gyroscopes and accelerometers) [16–19]. The dynamical Allan variance (DAVAR) is a sliding version of Allan variance, which could represent the nonstationary behavior of the signal [20, 21]. For MEMS sensor analyses, DAVAR could track and describe the dynamic characteristics of time series, and it is advantageous to analyze the process of gyroscope errors. The DAVAR is a cluster of Allan variance; therefore, the computational burden is very high because the DAVAR requires the computation of an Allan variance at every time instant [22, 23]. A recursive algorithm for DAVAR is proposed in [24]. In this fast DAVAR algorithm, the relationship of adjacent points of Allan variance has been revealed, then the calculation speed is accelerated. However, the Allan variance also needs to be calculated in many points of relevant times in this algorithm. In this paper, the MEMS gyroscope is regarded as the object of the analysis, and an improved method based on reasonable choice of relevant time is proposed to shorten the computation time further; the results show that the improved fast DAVAR algorithm dramatically reduces the computational time without affecting operation results.

The article is organized as follows: in Section 2, the theoretical basis and implementation process of the Allan variance and DAVAR are presented. In Section 3, the recursive algorithm for the DAVAR is introduced and the existing problems of this algorithm are discussed. Then, the improved fast algorithm of DAVAR is derived step by step. The experimental MEMS gyroscope signals are processed by the improved method in Section 4, followed by conclusion in Section 5.

#### 2. The Principle of Allan Variance and DAVAR

##### 2.1. Allan Variance

The calculation of Allan variance is based on the method of cluster analysis. Assuming that the signal of MEMS gyroscope is acquired at a sampling period , separating the sampling data into clusters and each cluster includes sampling data: where is the angular velocity and its unit is rad/s. The relevant time is defined as . The typical Allan variance can be expressed as where represents the average value of each cluster. By varying the number of samples per cluster, variances are computed at different cluster length by (2). The relationship existing between Allan variance and power spectrum density (PSD) of the intrinsic random process is given by where is the PSD of specified noise , namely, the instantaneous output rate of the gyro, and is the relevant time. Equation (3) indicates that the Allan variance is proportional to the total power output of the random process when passing through a filter with the transfer function of the form . This particular transfer function is the result of the method used to create and operate on the clusters [25].

Allan variance method could be used to analyze the common five basic gyro noise terms, including angle random walk, rate random walk, bias instability, quantization noise, and rate ramp. Each noise term could be accessed through the PSD associated with Allan variance and can be described as follows [26]: where is the quantization noise coefficient, is the angle random walk coefficient, is the bias instability coefficient, is the rate random walk coefficient, and is the drift rate ramp coefficient.

##### 2.2. Allan Variance Estimation Accuracy

In theory, Allan variance is a method which is used to analyze the stochastic process with finite length. Therefore, the estimation error is universal. Due to the finiteness of the divided clusters, the mean squared error of 1 can be obtained by a straightforward calculation as [27]

Combining the definition of (3), it can be seen that the confidence of Allan variance estimation improves as the number of independent clusters increases or the averaging time decreases.

##### 2.3. DAVAR

In the process of discrete calculation, angular increment is always used instead of angular velocity. Therefore, the Allan variance can be rewritten as where is the angular increment at time , which is the integral value of angular velocity in a period of time and can be described as . The unit of angular increment is rad.

The DAVAR is defined as a sliding version of the Allan variance, as Allan variance is a two-dimensional curve of , while the DAVAR is a three-dimensional figure which changes over time and relevant time . DAVAR can be described as where is the discrete-time analysis window.

#### 3. Improved DAVAR Method

As shown in (6), the DAVAR is obtained by computing the Allan variance at each analysis time epoch . With the length of signals increasing, it can result in a large computational burden. Therefore, the algorithm with less time consumption is urgently needed. In [24], the recursive property of the DAVAR was researched. The detailed process is described as follows.

Firstly, (7) can be redefined as where .

Then, the DAVAR at the next time epoch is given by

Comparing (8) and (9), the recursive function can be obtained as where and .

As a conclusion of the above analysis, when the computation time goes from to , the Allan variance at time epoch can be calculated by subtracting from the Allan variance at time epoch and adding to the Allan variance at time epoch . In this way, the computation time can be decreased to a certain extent. However, with the change of , the Allan variances at time epoch 1 which are based on different relevant times also need to be calculated. In this paper, an improved fast DAVAR algorithm based on the selection of relevant time is proposed.

As shown in (2), before calculating the Allan variance, the relevant time is selected based on . In other words, the relevant time is a linear sequence in linear coordinate system. However, for the convenience of intuitively analyzing the different stochastic noise terms, the Allan standard deviation is always plotted on a log-log scale, and the time axis is a logarithmic coordinate. For example, if and the relevant time is selected as linear sequence, the relevant time series distributing in log-log coordinate system is showed in Figure 1.