Temperature drift modeling and compensation of micro-electro-mechanical system gyroscope based on improved support vector machine algorithms

Abstract

This article suggested two methods to compensate for the temperature drift of the micro-electro-mechanical system gyroscopes, which are support vector machine method and C-means support vector machine. The output of X axis which was ranged from −40°C to 60°C based on the micro-electro-mechanical system gyroscope is reduced and analyzed in this article. The results showed the correctness of the two methods. The final results indicate that when the temperature is ranged from −40°C to 60°C, the factor of B is reduced from 0.424 $(° h)$ to 0.02194 $(° h)$ , and when the temperature is ranged from 60°C to −40°C, the factor of B is reduced from 0.1056 $(° h)$ to 0.0329 $(° h)$ , and the temperature drift trend and noise characteristics are improved clearly.

Keywords

Micro-electro-mechanical system gyroscope support vector machine C-means support vector machine temperature drift Allan variance

Introduction

The micro-electro-mechanical system (MEMS) gyroscope is an inertial sensor based on the technology of MEMS, which is used to sense the rotation rate of the carrier. The outstanding characteristics of the MEMS gyroscope are small size, light weight, low power consumption, and batch production. Low-precision MEMS gyroscope products have been very popular in the civil market and are widely used in consumer electronics and automotive fields. At the same time, the MEMS gyroscope is also widely used in the fields of attitude control and vehicle-borne satellite antenna. Most authors pay attention to improve the performance of MEMS gyroscopes. However, two main methods are used to improve their performance.^1,2

On one hand, one way is using hardware methods to reduce the temperature drift of the MEMS gyroscope. Therefore, most hardware methods have been suggested which contain structure and circuit. For instance, Cao et al.³ proposed that a bandwidth expansion method for dual-mass MEMS gyroscopes based on the pole-zero temperature compensation method. Cao et al.⁴ also proposed a sense mode closed-loop method for dual-mass MEMS gyroscopes based on the bi-pole temperature compensation method, and the experiment results show the bias temperature coefficient is decreased to 9.534°/h/°C and 5.991°/h/°C, respectively, based on sense open-loop and closed-loop circuits. Ou et al.⁵ investigated a new stress-released structure to improve the performance of temperature stability of MEMS gyroscope. Li et al.⁶ proposed a new method named frequency mismatch control to improve the temperature stability sensitivity for MEMS gyroscope. The results show that the scale factor temperature coefficient reduces from 1300 to 282.12 ppm/°C. Cao et al.⁷ proposed a new structure, which is a vibration ring gyroscope with a double U-beam and designed using a combination of mathematical analysis and the finite-element method. This new structure can easily improve performance. Although hardware methods have effectively improved the performance of MEMS gyroscopes under different temperature, these methods spend more time and amount of work.

On the other hand, the other way is to establish the mathematical model is to reduce temperature drift that is reported as software compensation. Fu et al.⁸ introduced the principle of MEMS gyroscope and the influence of MEMS gyroscope based on the temperature characteristic. Cao et al.⁹ proposed a method which called radial basis function neural network (RBF NN) to compensate the temperature energy influence. The results indicate that, after RBF NN-GA-KF method compensation, the bias instability of Gyros X, Y, and Z improves from $139 ° / h$ , $154 ° / h$ , and $178 ° / h$ to $2.9 ° / h$ , $3.9 ° / h$ , and $1.6 ° / h$ . Shen et al.¹⁰ proposed a multi-scale parallel model based on sample entropy based bounded ensemble empirical mode decomposition (SE-BEEMD), forward linear prediction (FLP) and extreme learning machine (ELM) for temperature error processing of MEMS gyroscope. The results show that the ELM has the best modeling accuracy and shortest training time, which would be valuable for online temperature drift modeling and compensation. Jianli et al.¹¹ proposed a novel error model based on RBF NN optimized by particle swarm optimization (PSO) to reduce the temperature error. This method is used in ring laser gyroscope (RLG), but it is not mentioned on MEMS gyroscope. Wei et al.¹² describes a novel approach for identifying the drift of RLG based on genetic algorithm and support vector regression. Although the method can easily decrease the drift of RLG based on room temperature, the drift of temperature was not processed in this article. Li et al.¹³ proposed a new algorithm method named support vector machine (SVM) which is used to reduce the temperature error based on the RLG. Li et al.¹⁴ proposed an SVM model to compensate for the temperature drift of the fiber optic gyroscope (FOG). The performance of the FOG was improved effectively. Also the least square support vector machine (LS-SVM) mode has the advantages like training speed, global optimum, high accuracy, and good non-linear properties. Wang and Chen¹⁵ proposed that a method of modeling and compensation for FOGs based on traditional SVM algorithms and the temperature drift of the FOG is greatly reduced. Although the algorithm has been already used to improve the performance of FOG and RLG, it is not mentioned on MEMS gyroscope and the temperature drift of MEMS gyroscope.

In this article, for the sake of improving the performance, a temperature drift model is established, and some methods are suggested to compensate for temperature drift. The algorithm of SVMs which is first used in temperature compensation of MEMS gyroscope is suggested to analyze and reduce the drift. Instead of a single algorithm, a novel improved SVMs algorithm is also suggested to reduce the drift. The Allan variance is suggested to compare the performance of compensation and improve the availability and practicability of the suggest methods.

Model of MEMS gyroscope and temperature

Structure model of MEMS gyroscope

As it is shown in Figure 1, the structure of this MEMS gyroscope can be divided into two modes: drive mode and sense mode, which can be described as a spring-mass-damping second-order system. It can be shown that the drive mode has 1 degree of freedom (DOF) along the X axis, and the sense mode also has 1 DOF along the Y axis. Each axis has two sections, which are stiffness $k$ and damping $d$ . As it is described in Figure 1, Coriolis mass $m$ has 2 DOF in both X axis and Y axis. As usual, MEMS gyroscopes work on the drive mode by drive stiffness $k_{x}$ and drive damping $d_{x}$ , when the Coriolis mass is subjected to an angular velocity of Z axis, the sense stiffness $k_{y}$ and sense damping $d_{y}$ are used by sensing mode to detect angular velocity. The dynamic system can be represented by the simplified equations

m \overset{\cdot\cdot}{x} + d_{x} \overset{\cdot}{x} + k_{x} x = F_{d} \sin (w_{d} t)

(1)

m \overset{\cdot\cdot}{y} + d_{y} \overset{\cdot}{y} + k_{y} y = - 2 m Ω_{z} \overset{\cdot}{x}

(2)

where $m$ is the main mass of this MEMS gyroscope, which used by drive mode and sense mode. $F_{d}$ is the drive force amplitude, $w_{d}$ is the drive force angular frequency, and $Ω_{z}$ is the angular rate around the Z axis.

Figure 1.

The structure of drive and sense mode.

Therefore, the natural frequency and the damping ratio of the drive mode should be described as

w_{n 1} = \sqrt{\frac{k_{x}}{m_{x}}}

(3)

ξ_{1} = \frac{d_{x}}{2 m_{x} w_{n 1}}

(4)

By solving equation (1), the displacement of the MEMS gyroscope in the drive direction can be described as

\begin{matrix} x (t) = & \frac{2 F_{d} ξ_{1} w_{n 1} w_{d} / m_{x}}{{(w_{n 1}^{2} - w_{d}^{2})}^{2} + 4 ξ_{1}^{2} w_{n 1}^{2} w_{d}^{2}} e^{- ξ_{1} w_{n 1} t \cos (w_{n 1} t \sqrt{1 - ξ_{1}^{2}})} \\ + \frac{F_{d} w_{d} (2 ξ_{1}^{2} w_{n 1}^{2} + w_{d}^{2} - w_{n 1}^{2}) / m_{x}}{w_{n 1} \sqrt{1 - ξ_{1}^{2}} [{(w_{n 1}^{2} - w_{d}^{2})}^{2} + 4 ξ_{1}^{2} w_{n 1}^{2} w_{d}^{2}]} e^{- ξ_{1} w_{n 1} t \sin (w_{n 1} t \sqrt{1 - ξ_{1}^{2}})} \\ + \frac{F_{d} / m_{x}}{w_{n 1}^{2} \sqrt{{(1 - \frac{w_{d}^{2}}{w_{n 1}^{2}})}^{2} + {(2 ξ_{1} \frac{w_{d}}{w_{n 1}})}^{2}}} \sin (w_{d} t - φ_{1}) \end{matrix}

(5)

where the phase and amplitude should be determined as

φ_{1} = \arctan \frac{2 ξ_{1} w_{n 1} w_{d}}{w_{n 1}^{2} - w_{d}^{2}}

(6)

A_{1} = \frac{F_{d} / m_{x}}{w_{n 1}^{2} \sqrt{{(1 - \frac{w_{d}^{2}}{w_{n 1}^{2}})}^{2} + {(2 ξ_{1} \frac{w_{d}}{w_{n 1}})}^{2}}}

(7)

When the MEMS gyroscope gets an angular velocity $Ω_{z}$ , the sense mode can be described as

F_{Ω_{Z}} = 2 m_{2} Ω_{z} \times \overset{\cdot}{x} = 2 m_{2} Ω_{z} A_{1} w_{d} \cos (w_{d} t - φ_{1})

(8)

where $F_{Ω_{Z}}$ is the Coriolis force which is generated from the angular velocity $Ω_{z}$ .

So, the equation of motion in sense mode is given as

F_{Ω_{Z}} = 2 m_{2} \overset{\cdot\cdot}{y} + d_{y} \overset{\cdot}{y} + k_{y} y = 2 m_{2} Ω_{Z} A_{1} w_{d} \cos (w_{d} t - φ_{1})

(9)

where the natural frequency and the damping ratio of the sense mode should be described as

w_{n 2} = \sqrt{\frac{k_{2}}{m_{2}}}

(10)

ξ_{2} = \frac{c_{2}}{2 m_{2} w_{n 2}}

(11)

By solving equation (1), the displacement of the MEMS gyroscope in the sense direction can be described as

\begin{matrix} y (t) = & \frac{2 m_{y} Ω_{Z} A_{1} w_{d} [2 ξ_{2} w_{n 2} w_{d} \sin φ_{1} + (w_{n 2}^{2} - w_{d}^{2}) \cos φ_{1}] / m_{y}}{{(w_{n 2}^{2} - w_{d}^{2})}^{2} + 4 ξ_{2}^{2} w_{n 2}^{2} w_{d}^{2}} e^{- ξ_{2} w_{n 2} t} \cos (w_{n 2} t \sqrt{1 - ξ_{2}^{2}}) \\ - \frac{2 m_{y} Ω_{Z} A_{1} w_{d} [ξ_{2} w_{n 2} (w_{n 2}^{2} - 3 w_{d}^{2}) \cos φ_{1} + w_{d} (2 ξ_{2}^{2} w_{n 2}^{2} + w_{n 2}^{2} - w_{d}^{2}) \sin φ_{1}] / m_{y}}{w_{n 2} t \sqrt{1 - ξ_{2}^{2}} [{(w_{n 2}^{2} - w_{d}^{2})}^{2} + 4 ξ_{2}^{2} w_{n 2}^{2} w_{d}^{2}]} e^{- ξ_{2} w_{n 2} t} \sin (w_{n 2} t \sqrt{1 - ξ_{2}^{2}}) \\ + \frac{2 m_{y} Ω_{Z} A_{1} w_{d} / m_{y}}{w_{n 2}^{2} \sqrt{{(1 - \frac{w_{d}^{2}}{w_{n 2}^{2}})}^{2} + {(2 ξ_{2} \frac{w_{d}}{w_{n 2}})}^{2}}} \cos (w_{d} t - φ_{1} - φ_{2}) \end{matrix}

(12)

where the phase and amplitude should be determined as

φ_{2} = \arctan \frac{2 ξ_{2} w_{n 2} w_{d}}{w_{n 2}^{2} - w_{d}^{2}}

(13)

\begin{matrix} A_{2} = & \frac{2 F_{d} w_{d} Ω_{Z}}{m_{x}} \cdot \frac{1}{w_{n 1}^{2} \sqrt{{(1 - \frac{w_{d}^{2}}{w_{n 1}^{2}})}^{2} + {(2 ξ_{1} \frac{w_{d}}{w_{n 1}})}^{2}}} \cdot \\ \frac{1}{w_{n 2}^{2} \sqrt{{(1 - \frac{w_{d}^{2}}{w_{n 2}^{2}})}^{2} + {(2 ξ_{1} \frac{w_{d}}{w_{n 2}})}^{2}}} \end{matrix}

(14)

The temperature model and drift of MEMS gyroscope

Temperature is one of the most significant factors that affects the performance of the MEMS gyroscope. There are two main factors that are to be considered: the size of the structure and the material of the elastic modulus when the temperature is changed. Among the two factors, the size of the structure has less influence on the performance of gyroscopes and the influence of the material of the elastic modulus should be mainly discussed.

The stiffness of the MEMS gyroscope system will be changed when the material of the elastic modulus is changed and the resonant frequency also changes with the stiffness of the system. Therefore, the change of temperature will lead to the drift of resonant frequency. Obviously, elastic modulus varies linearly with the temperature as shown in the following equation

E (T) = E_{0} - E_{0} κ_{ET} (T - T_{0})

(15)

where $E (T)$ is the elastic modulus when the temperature of silicon is $T$ , $E_{0}$ is the elastic modulus when the temperature of silicon is $T_{0}$ , and $κ_{ET}$ is a temperature variation coefficient of the elastic modulus of silicon. The stiffness of the system is proportional to the modulus of elasticity by the simplified equations

K = K_{0} [1 - κ_{ET} (T - T_{0})]

(16)

where $K$ is the stiffness of the system when the temperature is $T$ , $K_{0}$ is the stiffness of the system when the temperature is $T_{0}$ . So the relationship between the resonant frequency and temperature is described by the following equation

\begin{matrix} w_{d} (T) & = \sqrt{\frac{K}{m}} = \sqrt{\frac{K_{0} (1 - κ_{ET} (T - T_{0}))}{m}} \\ = w_{d} (T_{0}) [1 - \frac{1}{2 κ_{ET} (T - T_{0})}] \end{matrix}

(17)

So, we can see from these equations that temperature is one of the most significant factors which affects the performance of the MEMS gyroscope.

The experiment of temperature

The MEMS gyroscope is placed at the inside of a temperature-controlled oven for the sake of outputting the signal and not being affected by temperature variations of environment; the ranges of temperature and the temperature rate are set from −40°C to 60°C and 1°C/min, respectively. First, the initial temperature should be set as 60°C and placed in the temperature for an hour to ensure that the structural temperature is stable at 60°C. Then, the temperature is decreased at a rate of −1°C/min to −40°C and that temperature is also maintained for an hour to ensure that the structural temperature is of −40°C. The temperature-controlled oven should be observed at every 10°C for an hour to collect the data and be ensured that the structural temperature has a good stability and equal to oven temperature. The environment and steps are shown in Figure 2.

Figure 2.

The equipment of temperature experiment.

As for the test, the original output of an axis which described the output data reduce from 60°C to −40°C and rise from −40°C to 60°C is shown in Figure 3.

Figure 3.

The original output data of MEMS gyroscope: (a) $60 ° C$ to $- 40 ° C$ and (b) $- 40 ° C$ to $60 ° C$ .

The model of temperature and algorithms

The drift model of temperature based on algorithms

A novel temperature drift model is established in this article: temperature and temperature variation rate. These two factors could reflect the performance perfectly, therefore, $H$ is defined as the matrix of all two factors

H = {[\begin{matrix} T & \overset{\cdot}{T} \end{matrix}]}^{'}

(18)

where $T$ represents the temperature, and $\overset{\cdot}{T}$ represents the temperature variation rate. For the sake of establishing a model which is used in the new method, a function could be established as a target which can be trained as

f (x) = f {{[\begin{matrix} T & \overset{\cdot}{T} \end{matrix}]}^{'}}

(19)

Algorithm of SVMs

SVMs¹⁶ are based on the theory of statistical learning theory (SLT). SLT has several characteristics such as intuitive geometric explanations and perfect mathematical forms. Since Vapnik put forward the concept of SVM in the 1990s, the basic framework of algorithmic and theoretical concepts has been completely determined, which is based on many kinds of SVMs.

With the majority of data samples based on SVM, we can describe that the linear separable samples are considered as two classes, which are $(x_{1}, y_{1})$ , $(x_{2}, y_{2})$ ,…, $(x_{l}, y_{l})$ $x \in R^{n}$ , $y \in {- 1, + 1}$ ; these linear separable samples are classified by the hyperplane showed on the above equation

w • x + b = 0

(20)

and its standard form can be constrained by the expression of above equation

y_{i} (w • x_{i} + b) ⩾ 1, i = 1, 2, \dots, l

(21)

Figure 4 shows the classification of a sample hyperplane. We can focus on that the distance from the point of the data sample to the hyperplane is described as the equation

d = \frac{| w • x + b |}{‖ w ‖}

(22)

Figure 4.

The classification of sample hyperplane.

This distance is called the classifications distance, and the maximum distance classification is required. In order for the maximum classification interval to be equivalent to the minimum of values of $(1 / 2) | | w | |^{2}$ when the current classification interval is $2 / | | w | |$ , the above equation (21) should be satisfied which is the best classification. Therefore, the original appearance of the SVM can be represented by the following quadratic programming problem equations

\min \frac{1}{2} {‖ w ‖}^{2} + C \sum_{i = 1}^{l} ξ_{i}

(23)

s • t • y_{i} (x_{i} • w + b) ⩾ 1 - ξ_{i}

(24)

where $ξ_{i}$ is a slack variable when $ξ_{i} \geq 0$ , ans $C$ represents a penalty for the error of classification samples. So equation (22) can be represented as the following form of unconstrained Lagrange norm

L (w, b, a) = \frac{1}{2} ‖ w ‖^{2} - \sum_{i = 1}^{l} a_{i} [y_{i} (x_{i} • w + b) - 1]

(25)

where the Lagrange factor is $a_{i}$ , and equation (23) can also be expressed by the above equation

\frac{\partial L}{\partial b} = 0 \to \sum_{i = 1}^{l} a_{i} y_{i} = 0

(26)

\frac{\partial L}{\partial w} = 0 \to w = \sum_{i = 1}^{l} a_{i} x_{i} y_{i} = 0

(27)

We can also express the question of dual optimization by equations (27)–(29)

Q (a) = \sum_{i = 1}^{l} a_{i} - \frac{1}{2} \sum_{i, j = 1}^{l} a_{i} a_{j} y_{i} y_{j} (x_{i}, y_{j})

(28)

where $a_{i}$ is obtained by solving equation (8) because only a small part of $a_{i}$ is not found in the solution, which the corresponding sample is support vector. So the optimal classification function is given as

f (x) = sgn (w • x + b) = sgn (\sum_{i = 1}^{l} a_{i} y_{i} (x_{i} • x)) + b

(29)

where $b$ is threshold of classification which can be expressed by any support vector.

As for non-linear data samples, the non-linear transformation can be transformed into another linear question in high-dimensional space to solve, which needs to be achieved by the kernel function such as $K (x_{i}, x_{j}) = < Φ (x_{i}) • Φ (x_{j}) >$ . At this point, the objective function becomes

Q (a) = \sum_{i = 1}^{l} a_{i} - \frac{1}{2} \sum_{i, j = 1}^{l} a_{i} a_{j} y_{i} y_{j} K (x_{i}, x_{j})

(30)

and the function of classification becomes

f (x) = sgn [\sum_{x \in sv} a_{i} y_{i} K (x_{i} • x) + b]

(31)

In conclusion, these are SVMs. The kernel function is one of the most significant factors in the algorithm of SVM. The present kernel function makes the function solution bypass the eigenspace and get directly in the input space which avoids computational jealousy of non-linear mappings. There are some general kernel functions such as linear kernel function (LKF), polynomial kernel function (PKF), sigmoid kernel function (SKF), radial basis kernel function (RBKF), and spline kernel function (SPKF) which are described as follows

LKF

K (x, x_{i}) = (x, x_{i})

(32)

PKF

K (x, x_{i}) = {[(x, x_{i}) + r]}^{d}

(33)

SKF

\tanh [v • (x • x_{i}) + w]

(34)

RBKF

\exp (- \frac{{‖ x - x_{i} ‖}^{2}}{2 σ^{2}})

(35)

SPKF

1 + 〈 x, x_{i} 〉 + 1 / 2 〈 x, x_{i} 〉 min (x, x_{i}) - 1 / 6 {(x, x_{i})}^{3}

(36)

As is shown on these equations, LKF is a typical global kernel function that can be regarded as classifying the linearly separable data sets. PKF is an extensive applied non-linear mapping function, and the approaching capability of it is determined by the order of the polynomial. PKF is a d-order polynomial classifier which can be described as the above function

f (x) = sgn ({\sum_{SVs} y_{i} a_{i} [x_{i} • x + 1]}^{d} + b)

(37)

SKF actually includes a multi-layer perceptron network with a hidden layer, but the weight and the amount of hidden layer nodes are automatically determined by the algorithm, unlike the traditional perceptron network, which is determined by human experience. RBKF is the most extensive applied localized kernel function and strong in local regression. The smaller the value of $σ$ , the stronger the learning ability of SVM and the larger the confidence interval, but the SVM generalization ability gets worse at the same time. Therefore, the value of $σ$ could be properly selected. SPKF could approximate any non-linear functions theoretically, but the flexibility in practice is not best than RBKF.

Improved algorithm of SVMs

An improved SVM algorithm which is based on fuzzy C-means clustering is suggested to compensate for temperature drift more effectively. The SVM algorithm based on fuzzy C-means clustering includes three main processes which are shown in Figure 5. The processes are to use the fuzzy C-means clustering to preprocess the data of the learning sample. The second process is to describe the application of different classes SVM algorithm which is generated clustering by different classifications. The third process is to calculate the distance from the sample to each clustering center. Finally, the results of the whole system output for each calculation results of SVM algorithm. The steps of C-means support vector machine (C-SVM) method are shown in Figure 5.

Figure 5.

The steps of C-SVM algorithm.

The steps of C-SVM method is as follows:

Step1: Judging whether pre-clustering is necessary according to the size of learning sample data. As the number of clusters is larger than the set value, the number of clusters is calculated and the fuzzy C-means clustering is performed.

Step 2: As for each sample, find out the maximum membership degree of each cluster center and classify it into corresponding classes.

Step 3: After the clustering is completed, the number of samples contained in each class is calculated. And if the number is less than the corresponding value, the relevant class centers and the number of redundant classes will be eliminated. The samples are clustered again according to the new cluster number.

Step 4: Different SVMs are used to learn different classes.

Step 5: A fuzzy classifier is constructed to get the output of SVM. The input for each membership degree of sub-SVM is determined. According to the membership degree, the output of each sub-SVM model is synthesized and the final output is obtained

Y = \sum_{i = 1}^{c} u_{i} Y_{i}

(38)

where $C$ is the number of classes, $Y_{i}$ is the output of the first $i$ SVM, and $u_{i}$ is the test sample points to the membership degree of each cluster center. $u_{i}$ is related to the square of the distance from the sample point to the cluster centers. It can slightly ensure that the SVM of similar classes plays a greater role in the output of classification synthesis

u_{i} = {\begin{matrix} 1 & d (x, c_{i}) \neq 0 \\ \frac{1 / d^{2} (x, c_{i})}{\sum_{j = 1}^{c} 1 / d^{2} (x, c_{j})} & d (x, c_{j}) \neq 0 \end{matrix}

(39)

where $d (x, c_{i})$ expresses the Euclidean distance between the input and the corresponding clustering center, and $c_{i}$ expresses the center of class $i$ .

Computation complexity of SVM and C-SVM

In order to compare the computation complexity of SVM and C-SVM, computation time and space complexity is provided in this article. In order to analyze, we assume that the time spent for all operators is the same and the computation complexity concerns only performance and running hardware. After running the algorithms based on MATLAB, the results are shown in Table 1.

Table 1.

The results of computation time and space complexity based on SVM and C-SVM.

Algorithm	Computation time	Space complexity
SVM	$O ({LOG}_{2}^{N})$	O(N)
C-SVM	$O (2 {NLOG}_{2}^{2 N})$	O(N)

SVM: support vector machine; C-SVM: C-means support vector machine.

Verification and analysis

As for the test, the output data are shown in Figure 6. There are 2200 data points in each temperature experiment, with the majority of data points of the 220 data points are chosen to represent all data points (Figure 6). In order to improve the accuracy of C-SVM and SVM method, we add the output data of Y axis based on another MEMS gyroscope under the temperature experiment. The original data are shown in Figure 7, and Figure 8 includes 220 points of data.

Figure 6.

The output data of MEMS gyroscope: (a) $60 ° C$ to $- 40 ° C$ and (b) $- 40 ° C$ to $60 ° C$ .

Figure 7.

The original output data of another MEMS gyroscope: (a) 60°C to −40°C and (b) −40°C to 60°C.

Figure 8.

The output data of the MEMS gyroscope: (a) 60°C to −40°C and (b) −40°C to 60°C.

From Figures 9 and 10, we can see that the output of X axis varied a lot in the test. For the sake of reducing the temperature drift, SVM, an improved SVM algorithm method and the least square method are used in this article. From Figures 9 and 10, we also can see that the temperature drift model based on SVM algorithm, C-SVM algorithm, and the least square method are established.

Figure 9.

The output of the MEMS gyroscope range from −40°C to 60°C.

Figure 10.

The output of the MEMS gyroscope range from 60°C to −40°C.

From Figures 11 and 12, we can see that the output of Y axis varied a lot in the test. For the sake of reducing the temperature drift, SVM, an improved SVM algorithm method and the least square method are used in this article. From Figures 11 and 12, we also can see that the temperature drift model based on SVM algorithm, C-SVM algorithm, and the least square method are established.

Figure 11.

The output of another MEMS gyroscope range from −40°C to 60°C.

Figure 12.

The output of another MEMS gyroscope range from $60 ° C$ to $- 40 ° C$ .

Based on these temperature drift models, it can be easily recognized that the three curves which are the best algorithm to reduce the temperature drift. Obviously, C-SVM method is suitable to reduce the temperature drift.

Allan variance¹⁷ was proposed to distinguish the various error sources which are used to evaluate the performance of MEMS gyroscopes nowadays. The temperature drift will be divided into mainly five parts and the most important parts which are N (Angle random walk) and B (Bias instability). The method results of Allan variance analysis based on X axis and Y axis are shown in Tables 2 and 3.

Table 2.

The method results of Allan variance analysis based on X axis.

	–40°C to 60°C		60°C to −40°C
	$B (° h)$	$N (° h)$	$B (° h)$	$N (° h)$
Original data	0.424	8.65e–4	0.1056	2.15e–4
SVM data	0.0822	1.65e–4	0.0949	2.06e–4
C-SVM data	0.02194	4.69e–5	0.0329	6.04e–5
Least square	0.3254	6.54e–4	0.1021	1.67e–4

SVM: support vector machine; C-SVM: C-means support vector machine.

Table 3.

The method results of Allan variance analysis based on Y axis.

	–40°C to 60°C		60°C to −40°C
	$B (° h)$	$N (° h)$	$B (° h)$	$N (° h)$
Original data	0.066	1.0473e–4	0.4207	8.89e–4
SVM data	0.0317	6.8850e–5	0.0182	4.53e–5
C-SVM data	0.0145	3.4581e–5	0.0093	2.18e–5
Least square	0.0425	9.2562e–5	0.3703	7.89e–4

SVM: support vector machine; C-SVM: C-means support vector machine.

From Table 2, it can be seen that by using SVM and C-SVM methods, the factors of Allan variance analysis are reduced, which shows that C-SVM method is the better method.^18–24 As for bias instability, when the temperature is ranged from −40°C to 60°C, the factor of B is reduced from 0.424 $(° h)$ to 0.02194 $(° h)$ , and when the temperature is ranged from 60°C to −40°C, the factor of B is reduced from 0.1056 $(° h)$ to 0.0329 $(° h)$ . From Table 3, it also can be seen that by using SVM and C-SVM methods, the factors of Allan variance analysis are reduced, which shows that C-SVM method is the better method. As for bias instability, when the temperature is ranged from −40°C to 60°C, the factor of B is reduced from 0.066 $(° h)$ to 0.0145 $(° h)$ , and when the temperature is ranged from 60°C to −40°C, the factor of B is reduced from 0.4207 $(° h)$ to 0.0093 $(° h)$ .

Conclusion

In this article, two methods are suggested to reduce the temperature drift, which are SVM method and C-SVM. As for bias instability, when the temperature is ranged from −40°C to 60°C, the factor of B is reduced from 0.424 $(° h)$ to 0.02194 $(° h)$ , and when the temperature is ranged from 60°C to −40°C, the factor of B is reduced from 0.1056 $(° h)$ to 0.0329 $(° h)$ . The comparison results between two methods and the original data are shown obviously by Allan variance method.^17–21 Specifically, by using C-SVM, the drift of temperature is reduced more evidently.

Footnotes

Handling Editor: Antonio Lazaro

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by the National Natural Science Foundation of China (grant no. 51705477), and Pre-Research Field Foundation of Equipment Development Department of China (grant no. 61405170104). The research is also supported by program for the Top Young Academic Leaders of Higher Learning Institutions of Shanxi, Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shanxi Province, Shanxi Province Science Foundation for Youths (grant no. 201801D221195), Young Academic Leaders of North University of China (grant no. QX201809), and the Open Fund of State Key Laboratory of Deep Buried Target Damage (grant no. DXMBJJ2017-15).

ORCID iD

Xinwang Wang

References

Shen

Liu

Cao

, et al. Brain-like navigation scheme based on MEMS-INS and place recognition. Appl Sci 2019; 9: 1708.

Xia

Kong

The development of micromachined gyroscope structure and circuitry technology. Sensors 2014; 14: 1394–1473.

Cao

Zhang

Han

, et al. Pole-zero-temperature compensation circuit design and experiment for dual-mass MEMS gyroscope bandwidth expansion. IEEE-ASME T Mech 2019; 24(2): 677–688.

Cao

Liu

Zhang

, et al. Design and experiment of dual-mass MEMS gyroscope sense closed system based on bipole compensation method. IEEE Access 2019; 7: 49111–49124.

Hou

Miao

, et al. A new stress-released structure to improve the temperature stability of the butterfly vibratory gyroscope. Micromachines 2019; 10: 82.

Xiao

, et al. Enhanced temperature stability of sensitivity for MEMS gyroscope based on frequency mismatch control. Microsyst Technol 2017; 23: 3311–3317.

Cao

Liu

Kou

, et al. Design, fabrication and experiment of double U-beam MEMS vibration ring gyroscope. Micromachines 2019; 10: 186.

Chen

, et al. A temperature characteristic research and compensation design for micro-machined gyroscope. Mod Phys Lett B 2017; 31(6): 1750064.

Cao

Zhang

Shen

, et al. Temperature energy influence compensation for MEMS vibration gyroscope based on RBF NN-GA-KF method. Shock Vib 2018; 2018: 2830686.

10.

Shen

Zhang

, et al. Multi-scale parallel temperature error processing for dual-mass MEMS gyroscope. Sensor Actuat A-Phys 2016; 245: 160–168.

11.

Jianli

Feng

Jiancheng

, et al. Temperature error modeling of RLG based on neural network optimized by PSO and regularization. IEEE Sens J 2014; 14(3): 912–919.

12.

Wei

Long

Novel approach for identifying Z-axis drift of RLG based on GA-SVR model. J Syst Eng Electron 2014; 25(1): 115–121.

13.

Zhang

Wei

, et al. Multiple-point temperature gradient algorithm for ring laser gyroscope bias compensation. Sensors 2015; 15: 29910–29922.

14.

Chen

Yuan

, et al. Compensation of FOG temperature drift based on LS-SVM modeling. In: Proceedings of the 35th Chinese control conference, Chengdu, China, 27–29 July 2016. New York: IEEE.

15.

Wang

Chen

Temperature drift modeling and compensation of fiber optical gyroscope based on improved support vector machine and particle swarm optimization algorithms. Appl Opt 2016; 55(23): 6243–6250.

16.

Chang

Lin

CJ.

A library for support vector machines. ACM T Intel Syst Tec 2011; 2(3): 27:1–27:27.

17.

El-Sheimy

Hou

Niu

Analysis and modeling of inertial sensors using Allan variance. IEEE T Instrum Meas 2008; 57(1): 140–149.

18.

Cao

Liu

, et al. An improved interface and noise analysis of a turning fork microgyroscope structure. Mech Syst Signal Pr 2016; 70–71: 1209–1220.

19.

Guo

Tang

, et al. Determine turntable coordinate system considering its non-orthogonality. Rev Sci Instrum. 2019; 90: 033704.

20.

Xiao

Zhou

, et al. Dynamic modeling of the multiring disk resonator gyroscope. Micromachines 2019; 10: 181.

21.

Zhao

Huang

Probabilistic monitoring of sensors in state-space with variational Bayesian inference. IEEE T Ind Electron 2019; 66(3): 2154–2163.

22.

Wang

, et al. Research and application of improved adaptive MOMEDA fault diagnosis method. Measurement 2019; 140: 63–75.

23.

Zhao

Huang

Robust FIR state estimation of dynamic processes corrupted by outliers. IEEE T Ind Inform 2019; 15: 1.

24.

Cao

Shao

, et al. Sensing mode coupling analysis for dual-mass MEMS gyroscope and bandwidth expansion within wide-temperature range. Mech Syst Signal Pr 2018; 98: 448–464.