Temperature Energy Influence Compensation for MEMS Vibration Gyroscope Based on RBF NN-GA-KF Method

-is paper proposed three methods to compensate the temperature energy influence drift of the MEMS vibration gyroscope, including radial basis function neural network (RBF NN), RBF NN based on genetic algorithm (GA), and RBF NN based on GA with Kalman filter (KF). -ree-axis MEMS vibration gyroscope (Gyro X, Gyro Y, and Gyro Z) output data are compensated and analyzed in this paper. -e experimental results proved the correctness of these three methods, and MEMS vibration gyroscope temperature energy influence drift is compensated effectively. -e 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, respectively. And the angle random walk of Gyros X, Y, and Z was improved from 3.03°/h, 4.55°/h, and 5.89°/h to 1.58°/h, 2.58°/h, and 0.71°/h, respectively, and the drift trend and noise characteristic are optimized obviously.

e MEMS vibration gyroscope is one of the most important applications of MEMS vibration devices, and temperature energy influence drift restrains its high-precision application [11].Plenty of literatures focus on improving the temperature characteristics of the MEMS vibration gyroscope.However, there are mainly two methods to improve the temperature compensation of the MEMS vibration gyroscope.
One is employing hardware to improve the MEMS vibration gyroscope temperature performance, and a large number of methods of hardware compensation have been reported which include architecture improvement and circuit control.For example, Cao and Li investigated how the Gyro structure mechanical model is affected by temperature and proposed methods to improve the silicon structure temperature robustness [12].Liu et al. proposed robust structural design for the MEMS vibration gyroscope to compensate the temperature energy influence drift [13].Cao et al. established a new silicon structure equivalent electric model in order to improve the accuracy of performance under higher temperature [14].Fu et al. proposed a simpler and more applicable compensation method which is provided with a new circuit to reduce the influence of damping coefficient and resonating frequency [15].Yang et al. proposed an architecture that utilizes the on-chip temperature sensor of the micro-Gyro to achieve the onchip temperature compensation and used the integrated serpentine microheater to realize the on-chip temperature control of the micro-Gyro [16].Cao et al. proposed a temperature compensation bandwidth expanding method [17].Although hardware compensation is used in many MEMS vibration gyroscopes, more and more money and time is wasted.
e other is to establish the mathematical model to acquire temperature compensation which are called software compensations.Software compensations spend less money and time to improve the measurement precision, and its algorithm can be adjusted easily to improve the accuracy of the MEMS vibration gyroscope.Wang et al. proposed the relationship between temperature and frequency on the Hemispherical Resonator Gyro, and a method by changing the frequency is provided to compensate temperature energy influence drift [18].Aggarwal et al. proposed an efficient thermal variation model for low-cost MEMS inertial sensors and units in order to compensate the temperature [19].Song et al. proposed a new hybrid algorithm back propagation (BP) NN optimized by the artificial fish swarm algorithm (AFSA), and it is used to describe the temperature energy influence drift characteristic of FOG [20].Chen and Shen proposed a novel method based on the forward linear prediction (FLP) algorithm to decrease noise and a novel temperature compensation by using the ambient temperature change rate which can reduce the influence of intense temperature variation effectively [21].Feng et al. proposed an adaptive algorithm of strong tracking KF to compensate the bias drift of the MEMS vibration gyroscope in different temperatures and successfully improved the accuracy of the MEMS vibration gyroscope [22].Chong et al. proposed a modeling method based on Elman NN and GA with multiple temperature variable inputs in [23], and the temperature model was established based on temperature, temperature variation rate, and the coupling term.In [24], Wei et al. proposed GA and successfully used in the modeling of ring laser Gyro for temperature energy influence drift, and GA was employed to select the optimal parameters of support vector regression.Ding et al. proposed a modified RBF NN method of temperature compensation for laser Gyro and improved the accuracy of laser Gyro under different temperatures [25], and this method can quickly and accurately identify the effect of temperature on laser gyro zero bias.Cheng et al. proposed a modification of an RBF ANN based on temperature compensation models for IFOGs and improved the accuracy of temperature compensation [26], and three temperature-relevant terms were extracted including temperature of fiber loops, temperature variation of fiber loops, and temperature product term of fiber loops.Xia et al. proposed a temperature prediction and control based on BP NN and fuzzy-PID control method and successfully used to reduce the zero rate output [27].Software compensation based on MEMS vibration gyroscopes is more convenient than hardware compensation.
In this paper, in order to improve the accuracy of the MEMS vibration gyroscope, a temperature energy influence drift model is established, and three methods are proposed to compensate temperature energy influence drift.Instead of a single algorithm, RBF NN, RBF NN based on GA and fusion algorithms of RBF NN based on GA with KF which is first used in temperature compensation of MEMS vibration gyroscope are proposed to analyze and compensate the temperature energy influence drift.

Temperature Energy Influence Drift Model.
A new temperature energy influence drift model of MEMS is established in this paper: temperature, temperature variation rate, and temperature product term are considered to combine with the temperature energy influence drift model.
ese three factors can easily reflect the performance of the MEMS vibration gyroscope; therefore, D is defined as the matrix of three factors: where T is the temperature, _ T is the temperature variation rate, T × ΔT is the product term.After that, these three factors are all the inputs of the drift model as shown in Figure 1.In order to establish a complex model as inputs to NN, a function is constructed as a target function to be trained, as shown in the following equation: (2)

e Algorithm of RBF NN.
RBF is a kind of nonnegative and nonlinear functions which has features of local distribution and radial symmetry attenuation of the center point.Local distribution is described that the meaningful nonzero response will be achieved by the RBF of hidden unit when the input is placed in a small designated area.
Radial symmetry attenuation of the center point is described as the center of RBF with the equal input of the radial distance, and then, the equal output is produced by the RBF of hidden layer points, and the closer the distance between input and the center of RBF is, the bigger response of the hidden node will be.By using the hidden layer of neural work into RBF, the RBF NN is established.RBF NN, a typical multilayer feedforward neural network, is shown as Figure 2 which has a multilayer structure of sensationassociation-reaction that includes the input layer, the hidden layer, and the output layer.Compared with other neural networks, there are plenty of advantages which belong to the RBF NN.For example, in the RBF NN, the training time can be shorter, the functional approximation can be optimal, and the RBF NN can be approximated to arbitrary continuous function with arbitrary precision [26].e basis function of RBF NN is shown in the following equation: where c j is the center point of the basis function in j, σ j is the breadth of the center point around the basis function which is a freedom factor, ‖x − c j ‖ is the distance between x and c j which is the norm of x − c j , and ψ j is a function of radial symmetry which has a unique maximal value at c j , with the increasing of ‖x − c j ‖ and ψ j is quickly decreased to zero.For the determined output of x ∈ R n , only a little processing element is activated, where the center is closed to x. e transfer function of RBF neurons has many different forms; in this paper, Gaussian basis function is shown in the following equation: Shock and Vibration As is shown in Figure 3, RBF NN is formed with two parts: one part achieve the nonlinear mapping method from x to a j (x) � ψ j (‖x − c j ‖/σ j ) and the other part achieve the linear mapping from α j (x) to y, which is shown as ( e learning algorithm of connection weight is shown in the following equation: where β is the learning rate, and β is decided from 0 to 2 in order to erform the convergence of iterative learning algorithms. In the training process of NN, α j (x) is nearly equal to 0 when x is far from c j , and the output of α j (x) is also nearly equal to 0 when pass by the linear neural of the second layer.e corresponding weight value will be influenced when α j (x) is more than a value such as 0.05, the output of the layer is nearly equal to 1 when the distance between x and c j is small, the output value will nearly be equal to the weight value on the second layer when pass by the second layer.After these trainings, the advantage of the quick study speed of the local network is provided with RBF NN.

e Algorithm of RBF NN Based on GA.
GA is an intelligent method based on group optimization, which is followed by the rule of survival of the fittest to look for optimal solution and has advantages of global searching ability, being independent of initial value, and fast convergence rate.
e GA mainly includes five basic steps: parameter coding, initial population, fitness function choice, genetic operators, and trainable connection setting.e process of GA is shown in Figure 4, and the steps can be described as follows: (1) e initial set of RBF NN (initial population) is produced by a random number generator and consists of the GA parameters.

Input layer
Hidden layer Output layer

RBF NN Based on GA with KF.
e process of KF is shown in Figure 5, it indicates that the process of KF needs two steps: time upgrading and measurement upgrading, and the process of calculating can be divided into two steps which are plus loop and filter loop.In other words, the process is recursive.
erefore, the process of KF can be described as the following five equations: where H k is the measurement matrix, K k is the filter gain matrix, P k/k−1 is the matrix of prior estimate error covariance, P k is the matrix of estimation error variance, R is the covariance of measurement noise covariance, and Q is the covariance of process noise.From what have been discussed above, as long as the initial value of  X 0 and P 0 can be determined, based on the measurement vector of Z k at k moment, the state estimation of  X k (k � 1, 2, . ..) at k moment can be calculated.
As for two methods RBF NN and RBF NN + GA, a new method which is a fusion of RBF NN and GA with KF (RBF NN + GA + KF) is carried out in this section.e fusion algorithm of RBF NN + GA + KF not only provided the advantage of online real time but also provided the advantage of studying ability.
e flow chart of the fusion algorithm of RBF NN + GA + KF is shown in Figure 3. e steps of the fusion algorithm of RBF NN + GA with KF can be described as follows: (1) e initial parameters are produced by temperature experiment and set as input data to KF (2) KF model is established by the time series analysis method, and the output data of KF is ready for correcting the output data based on the method of RBF NN + GA.

Experiment Proposal
e structure and shell of LPMS-USBAL2 MEMS vibration gyroscope that 10 Hz of sampling rate is a system-inpackage featuring a 3D Gyro, which the bias stability is 0.05 °/s for Gyros.e main test device is a temperaturecontrolled oven, which can provide the ranges of temperature from −50 °C to +150 °C.
e equipment of the experiment setup is shown in Figure 6, which was used to collect above data.In order to test the temperature characteristic of the MEMS vibration gyroscope, a temperature experiment is set up.
e MEMS is kept on the temperature-controlled oven in order to output the signal of MEMS vibration gyroscopes and not influenced by outside temperature variation; then, the ranges of temperature and the temperature rate are set from −40 °C to 60 °C and 1 °C/min, respectively.Firstly, the initial temperature should be set as 60 °C and keep the temperature for an hour in order to make sure that the inside temperature of the MEMS vibration gyroscope is stabled at 60 °C.Secondly, the temperature is reduced at the speed of −1 °C/min into −40 °C and the temperature is kept for an hour to make sure that the inside temperature of the MEMS vibration gyroscope is −40 °C.e temperature-  Shock and Vibration controlled oven should be stayed at each 10 °C for an hour in order to collect the output of the MEMS vibration gyroscope and make sure that the inside temperature of the MEMS vibration gyroscope is stable and equal to oven temperature.

Verification and Analysis
Based on the temperature experiment of the MEMS vibration gyroscope, the output data of Gyro X, Gyro Y, and Gyro Z are shown in Figure 7.
From Figure 7, it can be observed that three MEMS vibration gyroscopes output signals varied about 3 °/s during −40 °C to 60 °C temperature range.In order to overcome the drift of temperature, three methods are used to compensate the temperature energy influence drift.From Figures 8-10, the temperature energy influence drift models based on RBF NN, RBF NN + GA and RBF NN + GA + KF are established.
Based on the established model, the compensation drift of three methods can be calculated.From Figures 8-10, it can be seen that the modal curves of three methods follow the original signal well, and the noise performance of RBF NN + GA + KF is the best.And based on this phenomenon, RBF NN + GA + KF method is used for compensation.
Figure 11 shows the compensation results of RBF NN + GA + KF method, and comparing with output data in Figure 7, RBF NN + GA + KF method successfully decreased the temperature influence drift.In order to have a quantitative evaluation of compensation result of the MEMS vibration gyroscope's temperature energy influence drift, the Allan variance method was applied to analyze the drift compensation result.e Allan variance curves are shown in Figure 12, and angle random walk (N) and bias instability (B) values are listed in Table 1.
From Table 1 and Figure 12, it can be seen that based on three methods, factors in Allan variance of B and N are declined, and the RBF NN + GA + KF method is better than the others.As for bias instability, Gyro X is from 139 °/h to 2.9 °/h, Gyro Y is from 154 °/h to 3.9 °/h, and Gyro Z is from 178 °/h to 1.6 °/h.As for angle random walk, Gyro X is from 3.03 °/h 1/2 to 1.58 °/h 1/2 , Gyro Y is from 4.55 °/h 1/2 to 2.58 °/h 1/2 , Gyro Z is from 5.89 °/h 1/2 to 0.71 °/h 1/2 .

Conclusion
In this paper, the MEMS vibration gyroscope temperature energy influence drift model is investigated, and three methods are proposed to compensate the temperature energy influence drift.RBF NN, RBF NN + GA, and RBF NN + GA + KF are proposed to compensate the temperature energy influence drift by using the temperature

( 2 )Figure 1 :
Figure 1: e framework of the temperature energy influence drift model.

Figure 2 :Figure 3 :
Figure 2: Model of radial basis function neural network.

( 3 )
Temperature energy influence drift model is established.(4) By the temperature energy influence drift model in step (3), RBF NN + GA steps are worked which is shown in Figure 3. (5) e original data after processing are output by RBF NN + GA and the correction of KF.

Figure 4 :
Figure 4: e process of RBF NN based on GA (RBF NN + GA).

Figure 6 :
Figure 6: e equipment of temperature test.

Figure 8 :Figure 7 :
Figure 8: e three temperature energy influence drift models of Gyro X.

6
Shock and Vibration energy influence drift model which is established in this paper.Firstly, RBF NN is established to improve the accuracy of temperature energy influence drift, and RBF NN + GA is used to optimize parameters of RBF NN to improve the precision.Secondly, a novel algorithm RBF NN + GA + KF is proposed to improve the accuracy of these two methods and the original output.en, a novel temperature energy influence drift model is established by temperature, temperature variation rate, and temperature product term.Finally, the comparison results among three methods and the original data are shown by Allan variance coefficients obviously.Specifically, by using RBF NN + GA

Figure 11 :Figure 10 :
Figure 11: e compensation results of Gyro X, Gyro Y, and Gyro Z based on three methods.