Next Article in Journal
Advances in Integration, Wearable Applications, and Artificial Intelligence of Biomedical Microfluidics Systems
Next Article in Special Issue
A Self-Oscillating Driving Circuit for Low-Q MEMS Vibratory Gyroscopes
Previous Article in Journal
Path Planning and Navigation of Miniature Serpentine Robot for Bronchoscopy Application
Previous Article in Special Issue
High-G MEMS Accelerometer Calibration Denoising Method Based on EMD and Time-Frequency Peak Filtering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature Drift Compensation for Four-Mass Vibration MEMS Gyroscope Based on EMD and Hybrid Filtering Fusion Method

1
Shanxi Software Engineering Technology Research Center, Taiyuan 030051, China
2
School of Software, North University of China, Taiyuan 030051, China
3
School of Instrument and Electronics, North University of China, Taiyuan 030051, China
4
Beijing Institute of Aerospace Control Devices, Beijing 100039, China
5
School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
6
Key Laboratory of Instrumentation Science & Dynamic Measurement, Ministry of Education, North University of China, Taiyuan 030051, China
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(5), 971; https://doi.org/10.3390/mi14050971
Submission received: 19 April 2023 / Revised: 26 April 2023 / Accepted: 27 April 2023 / Published: 28 April 2023
(This article belongs to the Special Issue MEMS Inertial Device)

Abstract

:
This paper presents an improved empirical modal decomposition (EMD) method to eliminate the influence of the external environment, accurately compensate for the temperature drift of MEMS gyroscopes, and improve their accuracy. This new fusion algorithm combines empirical mode decomposition (EMD), a radial basis function neural network (RBF NN), a genetic algorithm (GA), and a Kalman filter (KF). First, the working principle of a newly designed four-mass vibration MEMS gyroscope (FMVMG) structure is given. The specific dimensions of the FMVMG are also given through calculation. Second, finite element analysis is carried out. The simulation results show that the FMVMG has two working modes: a driving mode and a sensing mode. The resonant frequency of the driving mode is 30,740 Hz, and the resonant frequency of the sensing mode is 30,886 Hz. The frequency separation between the two modes is 146 Hz. Moreover, a temperature experiment is performed to record the output value of the FMVMG, and the proposed fusion algorithm is used to analyse and optimise the output value of the FMVMG. The processing results show that the EMD-based RBF NN+GA+KF fusion algorithm can compensate for the temperature drift of the FMVMG effectively. The final result indicates that the random walk is reduced from 99.608°/h/Hz1/2 to 0.967814°/h/Hz1/2, and the bias stability is decreased from 34.66°/h to 3.589°/h. This result shows that the algorithm has strong adaptability to temperature changes, and its performance is significantly better than that of an RBF NN and EMD in compensating for the FMVMG temperature drift and eliminating the effect of temperature changes.

1. Introduction

MEMS (micro-electro-mechanical systems) gyroscopes are widely used in both civil and military fields. They are mainly used in fields such as automobile navigation, inertial navigation, attitude determination, robotics, micro-signal detection, aerospace and spatial positioning, and application electronics [1,2,3,4,5,6]. However, temperature variation can seriously reduce the precision of the MEMS gyroscope when used in a changing environment. Therefore, it is necessary to study how to eliminate the influence of temperature drift on MEMS gyroscopes. In general, there are two main methods for the temperature compensation of MEMS gyroscopes: hardware compensation methods and software compensation methods.
The hardware method is mainly used to optimise the structure of the MEMS gyroscope. For example, Yu proposed to use a polarising resonator to improve the thermal stability of a resonator fibre optic gyro [7]. Gao proposed that by calculating the optimal layer number of the fibre coil, the temperature drift of the interferometric fibre optic gyroscope can be effectively reduced [8]. Wang proposed that the use of a symmetrical push-pull structure can eliminate the effects of temperature drift [9]. Yan proposed to use a single-polarisation fibre, and the resonant fibre optic gyro was found to exhibit good stability over a wide temperature range [10]. Ling proposed a novel winding idea of dicyclic fibre coils for suppressing the thermal-induced bias drift error of interferometric fibre optic gyroscopes [11]. Qian proposed a hybrid resonator consisting of a polymer-based long-range surface plasmon polariton (LRSPP) waveguide coupler and a silica fibre to reduce the effects of temperature noise on the resonant fibre optic gyroscope [12]. Zhang proposed an improved QAD fibre coil, which helped to reduce the thermal-induced drift error of a fibre optic gyroscope and improved its precision [13]. Zhang proposed a special resonator structure that can generate oriented thermal stress to reduce the temperature drift of the vibrating axis of the resonator [14]. To some extent, hardware compensation methods have the disadvantage of increasing the size and power consumption and cost of MEMS gyroscopes, so these methods are not suitable for some applications that require low cost and low power.
The second solution is a software compensation method that takes as its input the relationship between the temperature of the MEMS gyroscope and the output signal and then establishes a temperature drift model to compensate for the effects of temperature drift. For example, Fontanella proposed a new method for accurate thermal compensation of the inertial device gyro bias using an augmented state Kalman filter [15]. Cao proposed a fusion algorithm based on an RBF NN-GA-KF, which effectively reduced the effect of the temperature drift on the MEMS vibration gyroscope [16]. Baxpameeb proposed an improved method for real-time thermal drift compensation of fibre optic gyros, through which the fibre gyro temperature drift was effectively compensated [17]. Narasimhappa proposed an ARMA model based on an adaptive unscented fading Kalman filter to reduce the drift of a fibre-optic gyroscope [18]. Zhang proposed a novel wavelet threshold denoising method, which effectively suppresses the environmental drift of fibre optic gyroscopes [19]. Song proposed a new modelling and compensation method combining an artificial fish swarm algorithm (AFSA) and a back-propagation (BP) neural network to improve the output precision of fibre optic gyroscopes [20]. Han proposed a simplified model to compensate for the output of ring laser gyros, which effectively improved the output accuracy of the gyroscopes [21]. Ahtohoba used factors such as the temperature change rate of the environment; an error model of the fibre optic gyroscope was established, and the output precision of the gyro was ultimately improved [22]. Zha established an IUKF neural network model for temperature drift compensation of the fibre optic gyro [23]. Prikhodko proposed to compensate high-Q MEMS gyroscopes by using temperature self-sensing technology, and the compensated gyroscope has high precision [24]. Chen proposed a fibre optic gyro adaptive positive linear prediction denoising algorithm based on the ambient temperature rate of change and a temperature drift modelling compensation concept, which was used to correct errors caused by drastic changes in the ambient temperature [25]. Wang proposed to establish an HRG temperature compensation model based on the natural frequency of the resonator and then used this model to compensate for the temperature drift of the gyroscope; this method can be used over a wide temperature range [26]. Cai established a parallel processing model based on multi-objective particle swarm optimisation based on a variational modal decomposition-time-frequency peak filter (MOVMD-TFPF) and Beetle antennae search algorithm-Elman neural network (BAS-Elman NN) to eliminate the noise and temperature drift in a Micro-Electro-Mechanical Systems (MEMS) gyroscope’s output signal [27]. Zhou proposed the improved VMD and TFPF hybrid denoising algorithm, which combines variational modal decomposition (VMD) and time-frequency peak filtering (TFPF), which has a smaller signal distortion and stronger denoising ability [28]. In recent years, the academic circle has shown great interest in the research of a four-mass gyroscope. Zhou proposed a novel Center Support Quadruple Mass Gyroscope, which combines the advantages of tuning fork gyros and micro hemisphere resonant gyros and is expected to achieve performance breakthroughs of flat structures [29]. Trusov proposed a new tuning fork architecture addressing the limitations of the conventional designs [30]. The software method has the advantages of a simple structure, low cost, and convenient implementation, and it is convenient for improving the precision of the MEMS gyroscope.
In this paper, an improved RBF NN filtering algorithm based on EMD is proposed, and this method has been successfully applied to eliminate the temperature drift of the MEMS gyro. Different from the other MEMS gyroscope temperature compensation algorithms, this paper proposes a new fusion algorithm (EMD-based RBF NN+GA+KF fusion algorithm) to eliminate the influence of temperature drift on MEMS gyroscopes and make the MEMS gyroscope more accurate.
The structure of the MEMS gyroscope is introduced in Section 2; the fusion algorithm is shown in Section 3; Section 4 describes the temperature experiment; and Section 5 gives the results of the data processing and a comparison of the various algorithms. The conclusion of this paper is given in Section 6.

2. Structure of the Four-Mass Vibration MEMS Gyroscope

The original signal is generated by a four-mass vibration MEMS gyroscope (FMVMG). The schematic diagram of the FMVMG is shown in Figure 1a. The FMVMG is composed of support anchor points, quadruple mass blocks, a drive comb, a drive sense comb, a sense comb, a force rebalances comb, and multiple folding beams. This gyroscope has two operating modes: a drive mode and a sense mode. Both modes can be considered as a “mass-spring-damped” second-order vibration system. The basic mechanical equivalent model is shown in Figure 1b.
Firstly, after the preliminary design and optimisation, the geometric shape of the FMVMG was determined, and the specific parameters are shown in Table 1. Then, the modal component in ANSYS software was used to analyse the modal of the FMVMG, and the modal analysis mainly included the natural frequency, modal shape, and vibration stability of the MEMS gyroscope structure. Finally, the driving and sensing mode shapes of the FMVMG are shown in Figure 2.
As shown in Figure 2, the vibration mode of the driving mode is the motion of the Coriolis mass block on the X-axis, while the vibration mode of the sensing mode is the motion of the Coriolis mass block on the Y-axis. In both the driving and sensing modes, the motion of the adjacent mass blocks is reversed. Among them, the resonant frequency of the driving mode is 30,740 Hz, and the resonant frequency of the sensing mode is 30,886 Hz due to the highest mechanical sensitivity when the natural frequencies of the driving and sensing modes of the MEMS gyroscope are equal, but the bandwidth is small. The frequency difference of this structure is 146Hz, which not only ensures the mechanical sensitivity of the gyroscope but also gives it a certain bandwidth.
The gyroscope control and detection system are shown in Figure 3. In the drive loop, the drive frame displacement, x(t), is detected by the drive sensing combs and picked up by the differential amplifier ①. Then, the signal phase is delayed by 90° (through ②) to satisfy the phase requirement of the AC drive signal, VdacSin(ωdt). After that, VdacSin(ωdt) is processed by a full-wave rectifier ③ and a low pass filter ④. Afterwards, Vdac is compared (in ⑤) with the reference voltage, Vref ⑥. Next, the drive PI controller ⑦ generates the control signal, which is modulated by VdacSin(ωdt), and then the signal is superposed (through ⑩) by VDC ⑨ to stimulation drive mode.
The sensing system employs an open loop, which utilises the same interface as the drive circuit. Firstly, the left and right masses’ sensing signals are detected separately with the differential detection amplifier ⑪. In addition, the output signals are processed by a second differential amplifier ⑫ to generate a signal, Vstotal. Then, Vstotal is demodulated by signal VdacSin(ωdt) (in ⑬). After that, the demodulated signal, Vdem, passes through the low pass filter ⑭ so the sensing mode’s movement signal, VOopen, can be obtained.
The monitoring system consists of three PCBs. PCBI contains the connection circuit and is connected to the structure chip, and PCB II is the drive circuit. PCB I and PCB III are nested on the top and back sides, respectively. PCB III contains the detection loop, and the output signal is from PCB III. The prototype of the FMVMG is shown in Figure 4.

3. Model and Algorithm

3.1. SE-EMD Algorithm

Empirical mode decomposition is a new adaptive signal time-frequency processing method that is especially suitable for the analysis and processing of nonlinear non-stationary signals. In essence, this method is also a process of smoothing non-stationary signals by generating a series of data sequences, and their feature scales are different. Each of these data sequences is called an intrinsic mode function (IMF). The specific steps of the EMD are as follows:
(1) Take out all the local extremum points of the signal and use the cubic spline curve to connect all the local maxima to form the upper envelope. Similarly, the minimum values are also connected to form the lower envelope. In addition, find the average value of the upper and lower envelopes, denoted as m1, and h1 can be obtained as follows:
x ( t ) m 1 = h 1
where x(t) is the original data, and if h1 meets the IMF component conditions, then h1 can be seen as the first IMF component of x(t).
(2) If h1 is not the IMF, set h1 as the original data, continue with step 1, obtain the average value m11, and evaluate h11 = h1m11; if h11 is still not an IMF component, repeat the loop k times and obtain h1(k−1)m1k = h1k until h1k meets the conditions of the IMF. Note that c1 = h1k. Then, c1 is the first of the signals x(t) that meets the IMF condition.
(3) Separate c1 from x(t) and obtain Equation (2):
r 1 = x ( t ) c 1
(4) Substitute the value solved by Equation (2) into Step 1 and Step 2, obtain the second IMF-compliant component, and repeat this step n times to obtain the nth IMF-compliant components of the signal, namely:
r 1 c 2 = r 2 r n 1 c n = r n
(5) If the resulting rn is a monotonic function and no additional IMF-compliant components can be derived from it, the steps are not repeated. The final result of EMD decomposition is shown in Equation (4):
x ( t ) = i = 1 n c i + r n
where c1, c2, c3, …, cn are the individual IMF components, and rn represents the average trend of the signal.
After EMD decomposition, the output signal of the gyroscope is decomposed into many IMF components. To simplify the analysis, the IMFs are divided into three component types: a noise component, a mixed component, and a drift component. The sample entropy (SE) is introduced here to realise this work. The sample entropy is a measure of the complexity of a time series. The larger the value of SE is, the more complex and irregular the data. SE can be defined as follows:
S E ( m , r ) = lim N { l n [ A m r B m r ] }
where m is the embedding dimension, referred to as the length of the data to be compared, and r is the tolerance value.

3.2. The Algorithm of RBF NN+GA+KF

The radial basis function (RBF) is a non-negative nonlinear function characterised by a local distribution of the centre point and a symmetric attenuation of the radial centre point. The smaller the distance between the input point and the centre of the radial basis function, the stronger the output signal.
The basis function of the RBF NN is shown in Equation (6):
α j ( x ) = ψ j ( x c j / σ j )
where x is an n-dimensional input vector, cj is the centre of the jth hidden node, σj is the breadth of the hidden layer neuron, ||·|| is the Euclidean distance, and ψj has a maximum value at cj, which is a radially symmetric function. When ||x−cj|| increases, ψj decreases quickly.
To establish an RBF neural network model, an RBF is introduced into the hidden layer of the neural network. The RBF NN structure is shown in Figure 5, which is similar to the other forward neural networks. The RBF neural network has three layers: the input layer, the hidden layer, and the output layer. The input layer consists of signal source nodes. The number of elements in the hidden layer is determined according to the specific problem, and the hidden layer centre transform function is a radial basis function. The output layer is the response to the input mode.
Here, the Gaussian basis function is chosen as the transfer function of the RBF neurons, which is shown in Equation (7):
a j ( x ) = ψ j ( x c j / σ j ) = e x c j 2 2 σ j 2
The output of the RBF NN is shown in Equation (8):
y i = j = 1 m w i j α j ( x ) , i = 1 , 2 , , r
where Wij is the connection weight between the ith hidden layer node and the jth output layer node. The learning algorithm is shown in Equation (9):
w i j ( l + 1 ) = w i j ( l ) + β [ y i d y i ( l ) ] α j ( x ) / α T ( x ) α ( x )
where β is the learning rate. When β falls within the range of 0 to 2, the algorithm offers good convergence performance.
During the training process of the RBF NN, the local amplification of the RBF neural network occurs. If x is far from cj, the output of αj(x) is also nearly 0 when passing through the linear nerve of the second layer. When x and cj become very close to 0, the output of the layer is almost equal to 1, and the output value is almost equal to the weight value on the second layer when passing through the second layer
The RBF neural network has the characteristics of a simple structure, fast operation speed, and local function approximation. In general, the more hidden layer neuron nodes of the RBF neural network, the stronger the computing power and mapping ability are and the better the function is. The approximation ability can approach a complex function curve with arbitrary precision.
However, the performance of the RBF NN model depends mainly on the length of the training sample set and the expansion constant of the radial basis function (breadth σj). The length of the training sample set affects the performance of the network and the real-time performance of the training directly; if the expansion constant, σj, is too small, the network will over-fit, and, conversely, if it is too large, the network will not have the ability to fit. Therefore, genetic algorithms are introduced to solve this problem.
The genetic algorithm is an adaptive global probability search algorithm that simulates the inheritance and survival of creatures in their natural environment. The algorithm mainly includes five basic steps: coding, initial population generation, fitness setting function, genetic operator design, and operational parameter determination. The process of a GA is shown in Figure 6. In combining the GA with the RBF neural network, use the GA to solve the parameter selection problem of the RBF NN model. The specific algorithm is as follows: set the maximum algebra, GMAX, of the genetics and combine each body in the population as RBF NN training parameters to find the corresponding output error. After training the algebra maximum, stop the calculation.
The Kalman filter algorithm is an optimal estimation of the autoregressive data processing algorithm. It is widely used in the elimination of gyro noise because of its high estimation accuracy, fast convergence speed, strong adaptability, and simple calculation. Since the original data contain the influence of noise, the optimal choice of data processing can be regarded as a filtering process. The Kalman filtering process can be approximated as the following five equations:
X ^ k / k 1 = Φ k , k 1 X ^ k 1
X ^ k = X ^ k / k 1 + K k ( Z k H k X ^ k / k 1 )
K k = P k / k 1 H k T ( H k P k / k 1 H k T + R ) 1
P k / k 1 = Φ k , k 1 P k 1 Φ k , k 1 T + Γ k 1 Q Γ k 1 T
P k = ( I K k H k ) P k / k 1
where Hk is a measurement matrix, Pk/k−1 is a matrix of the prior estimate error covariance, Pk is a matrix of the estimation error variance, Kk is a filter gain matrix, R is the covariance of the measurement noise, and Q is the covariance of the process noise. According to the above equations, it can be found that the state estimation of X ^ k ( k = 1 , 2 , ) at time point k can be calculated by the measurement vector of Zk at time point k if the initial value of X ^ 0 and P0 are determined.
In this section, the RBF NN, GA, and KF are combined, and the RBF NN+GA+KF fusion algorithm is proposed. The algorithm has the advantages of online real-time use, training capability, and fast learning. The fusion algorithm of RBF NN+GA+KF is shown in Figure 7.

3.3. RBF NN+GA+KF Denoising Based on SE-EMD

After being decomposed by EMD, the original data are decomposed into multiple IMF components. Each IMF component is classified by using the sample entropy (SE). Generally, these IMF components can be divided into three groups: the noise component, the mixture component, and the drift component. The noise has the characteristics of high frequency and low drift. The high-frequency noise component of the IMFS can be removed directly. However, the mixed component contains both the useful signal and the equivalent noise signal, so the RBF NN+GA+KF is applied to the filtering of this component, which is beneficial for protecting the useful signal of the gyroscope and removing the noise signal effectively. The processing of the drift component can be divided into two steps. First, the output of the temperature sensor is decomposed by EMD to obtain the trend term of the temperature, and then the trend term of the temperature is used to compensate for the drift component of the gyroscope. Finally, the signal reconstruction of the MEMS gyroscope is performed by using the mixed component after denoising and the drift component after compensation by the trend term of the temperature. The specific flow chart is shown in Figure 8.

4. Temperature Experiment Proposal

Through temperature experiments, the temperature characteristics of the MEMS gyroscope were tested. In the temperature experiment, the temperature-controlled oven can accurately provide a temperature range of −40 °C to +150 °C, as shown in Figure 9. The gyroscope output is collected by dedicated data acquisition software. Power is provided by a Gwinstek GPS-4303C DC power supply. The real-time temperature in the gyroscope metal casing is obtained by a temperature sensor, whose value is synchronised with the gyroscope output. The indoor temperature is 25 °C.
Firstly, the gyroscope is fixed to the static plane to avoid the influence of motion. Secondly, the output lines are connected to the Gwinstek DC power supply and the laptop. Then, the oven temperature is set at −30 °C to +60 °C. Finally, the power is switched on, and the original output signal of the gyroscope is recorded. The data collection process is continuous, and the temperature experiment results are shown in Figure 10. It can be seen that the gyroscope output changes more obviously when the temperature changes. This phenomenon indicates that the temperature affects the precision of the gyroscope output seriously. Therefore, it is necessary to establish a temperature drift model of the MEMS gyroscope.

5. Verification and Analysis

Figure 11 shows the results of the EMD processing of the gyroscope temperature drift error, which is decomposed into 7 IMF components and a residual component. These IMF components are characterised by non-stationary and complex components. For the convenience of processing, the sample entropy (SE) method is used to classify the IMF components of each layer, which can reduce the cumulative error when processing a single IMF reconstruction. Figure 12 shows the feature extraction based on the SE. It can be seen from Figure 13 that these IMF components can be roughly divided into three groups: IMF1, IMF2 as a group; IMF3, IMF4 as a group; and IMF5 to IMF7 as a group. These three groups are the three signal components mentioned above: the noise component (S1), the mixed component (S2), and the drift component (S3). The three signals are processed separately below.
According to the temperature error handling method proposed in Figure 8, first, the noise component (S1) is removed due to the high-frequency and low-drift characteristics of the noise, and the mixed component (S2) is denoised by the RBF NN+GA+KF algorithm. The denoised results are shown in Figure 14a. Second, after the temperature sensor output signal is decomposed by EMD, as shown in Figure 13, the temperature trend is obtained. Then, the temperature component is used to compensate for the drift component of the gyroscope output. Finally, the S2 and S3 component signals after processing are reconstructed, and the reconstruction results are shown in Figure 14b.
It can be seen from Figure 14 that the EMD improved algorithm proposed in this paper can effectively eliminate the error caused by temperature drift.
To compare the temperature compensation results quantitatively, the Allen analysis of variance was introduced. The Allen analysis is an analysis method that is based on the time domain. The main feature of the Allan variance method is that it can easily characterise and identify various error sources and their contribution to the overall noise statistics, and it also has the advantages of easy calculation and easy separation. Figure 15 shows the different compensation results for the Allan standard deviation curve, and Table 2 gives the quantitative results. It can be concluded from Table 2 that the proposed EMD-based RBF NN+GA+KF fusion method has the best performance, and the bias stability is increased from 34.66°/h to 3.589°/h. The angular random walk is reduced from 99.608°/h/Hz1/2 to 0.967814°/h/Hz1/2.

6. Conclusions

This paper studies the design and analysis of a new four-mass vibration MEMS gyroscope (FMVMG). The temperature compensation algorithm of the MEMS gyroscope is also studied in detail. First, the structure is simulated by finite element analysis software. The simulation results show that the FMVMG has two working modes: a driving mode and a sensing mode. The resonant frequency of the driving mode is 9.6096 kHz, the resonant frequency of the sensing mode is 9.6154 kHz, and the frequency interval between those two modes is 5.8 Hz. In the temperature compensation, a new fusion algorithm based on empirical mode decomposition (EMD) and combining a radial basis function neural network (RBF NN), genetic algorithm (GA), and Kalman filter (KF), is proposed. First, the gyroscope output is decomposed into three components by EMD and the sample entropy (SE), which are the noise component, the mixed component, and the drift component. Second, to process these three components, the specific steps are as follows: the high-frequency noise component is directly removed, and the RBF NN+GA+KF algorithms are used to denoise the mixed components. Third, the trend component of the temperature sensor output is used to compensate for the drift component of the MEMS gyroscope. Finally, the signal is reconstructed to obtain the final compensated signal. The results show that this method not only removes the noise signal of the gyroscope effectively but also protects the effective signal of the gyroscope. Finally, the Allan variance coefficient shows the comparison of the method with the original data. The data show that if the gyroscope output is compensated by the EMD and RBF NN + GA with the KF fusion method, the bias stability is increased from 34.66°/h to 3.589°/h. The angular random walk is reduced from 99.608°/h/Hz1/2 to 0.967814°/h/Hz1/2.

Author Contributions

Z.L. and Y.C. conceived and designed the experiments; K.C., H.C., Y.G. and G.W. performed the experiments; J.Y. and Z.L. analysed the data; Z.L. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key Research and Development Program of China (No. 2022YFB3205000), the National Natural Science Foundation of China, the NSAF (Grant No. U2230206), the Technology Field Fund of Basic Strengthening Plan of China (2020JCJQJJ409, 2021-JCJQ-JJ-0315), the National Defense Basic Scientific Research Program (WDZC20190303), and the Pre-Research Field Foundation of Equipment Development Department of China (No. 61405170104 and No. 80917010501). The research is also supported by the Fundamental Research Program of Shanxi Province (20210302123020 and 20210302123062), the Shanxi Province Key Laboratory of Quantum Sensing and Precision Measurement (201905D121001), the Key Research and Development (R&D) Projects of Shanxi Province (202003D111004), the Beijing Key Laboratory of High Dynamic Navigation Technology Open Founding (HDN2021102), and the Fund for Shanxi “1331 Project” Key Subjects Construction.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pribanić, T.; Petković, T.; Đonlić, M. 3D registration based on the direction sensor measurements. Pattern Recognit. 2019, 88, 532–546. [Google Scholar] [CrossRef]
  2. Davidson, P.; Virekunnas, H.; Sharma, D.; Piché, R.; Cronin, N. Continuous Analysis of Running Mechanics by Means of an Integrated INS/GPSDevice. Sensors 2019, 19, 1480. [Google Scholar] [CrossRef]
  3. Guo, Y.; Wang, P.; Ma, G.; Wang, L. Envelope oriented singularity robust steering law of control moment gyros for spacecraft attitude maneuver. Trans. Inst. Meas. Control 2019, 41, 954. [Google Scholar] [CrossRef]
  4. Li, Y.; Hu, H.; Zhou, G. Using Data Augmentation in Continuous Authentication on Smartphones. IEEE Internet Things J. 2018, 6, 628–640. [Google Scholar] [CrossRef]
  5. Jauhiainen, M.; Puustinen, J.; Mehrang, S.; Ruokolainen, J.; Holm, A.; Vehkaoja, A.; Nieminen, H. Identification of Motor Symptoms Related to Parkinson Disease Using Motion-Tracking Sensors at Home (KAVELI): Protocol for an Observational Case-Control Study. JMIR Res. Protoc. 2019, 8, e12808. [Google Scholar] [CrossRef] [PubMed]
  6. Voicu, R.A.; Dobre, C.; Bajenaru, L.; Ciobanu, R.I. Human Physical Activity Recognition Using Smartphone Sensors. Sensors 2019, 19, 458. [Google Scholar] [CrossRef] [PubMed]
  7. Yu, X.; Ma, H.; Jin, Z. Improving thermal stability of a resonator fiber optic gyro employing a polarizing resonator. Opt. Express 2013, 21, 358. [Google Scholar] [CrossRef] [PubMed]
  8. Gao, Z.; Zhang, Y.; Gao, W. Method to determine the optimal layer number for the quadrupolar fiber coil. Opt. Eng. 2014, 53, 084106. [Google Scholar] [CrossRef]
  9. Wang, J.; Yu, Y.; Chen, Y.; Luo, H.; Meng, Z. Research of a double fiber Bragg gratings vibration sensor with temperature and cross axis insensitive. Opt.—Int. J. Light Electron Opt. 2015, 126, 749–753. [Google Scholar] [CrossRef]
  10. Yan, Y.; Ma, H.; Jin, Z. Reducing polarization-fluctuation induced drift in resonant fiber optic gyro by using single-polarization fiber. Opt. Express 2015, 23, 2002–2009. [Google Scholar] [CrossRef]
  11. Ling, W.; Li, X.; Xu, Z.; Wei, Y. A Dicyclic Method for Suppressing the Thermal-induced Bias Drift of I-FOGs. IEEE Photonics Technol. Lett. 2015, 28, 272–275. [Google Scholar] [CrossRef]
  12. Qian, G.; Fu, X.C.; Zhang, L.J.; Tang, J.; Liu, Y.R.; Zhang, X.Y.; Zhang, T. Hybrid fiber resonator employing LRSPP waveguide coupler for gyroscope. Sci. Rep. 2017, 7, 41146. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Z.; Yu, F. Analysis for the thermal performance of a modified quadrupolar fiber coil. Opt. Eng. 2018, 57, 017109. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Wu, Y.; Xu, X.; Xi, X.; Wu, X. Research on the method to improve the vibration stability of vibratory cylinder gyroscopes under temperature variation. Int. J. Precis. Eng. Manuf. 2017, 18, 1813–1819. [Google Scholar] [CrossRef]
  15. Fontanella, R.; Accardo, D.; Moriello, R.S.L.; Angrisani, L.; Simone, D.D. An Innovative Strategy for Accurate Thermal Compensation of Gyro Bias in Inertial Units by Exploiting a Novel Augmented Kalman Filter. Sensors 2018, 18, 1457. [Google Scholar] [CrossRef] [PubMed]
  16. Cao, H.; Zhang, Y.; Shen, C.; Liu, Y.; Wang, X. Temperature Energy Influence Compensation for MEMS Vibration Gyroscope Based on RBF NN-GA-KF Method. Shock. Vib. 2018, 2018, 2830686. [Google Scholar] [CrossRef]
  17. Vahrameev, E.I.; Galyagin, K.S.; Oshivalov, M.A.; Savin, M.A. Method of numerical prediction and correction of thermal drift of the fiber-opticgyro. Izv. Vuzov. Priborostr. 2017, 60, 32–38. [Google Scholar] [CrossRef]
  18. Narasimhappa, M.; Nayak, J.; Terra, M.H.; Sabat, S.L. ARMA model based adaptive unscented fading Kalman filter for reducing drift of fiber optic gyroscope. Sens. Actuators A Phys. 2016, 251, 42–51. [Google Scholar] [CrossRef]
  19. Zhang, Q.; Wang, L.; Gao, P.; Liu, Z. An Innovative Wavelet Threshold Denoising Method for Environmental Drift of Fiber Optic Gyro. Math. Probl. Eng. 2016, 2016, 9017481. [Google Scholar] [CrossRef]
  20. Song, R.; Chen, X.; Shen, C.; Zhang, H. Modeling FOG Drift Using Back-Propagation Neural Network Optimized by Artificial Fish Swarm Algorithm. J. Sens. 2014, 2014, 276043. [Google Scholar] [CrossRef]
  21. Han, X.; Hu, S.M.; Luo, H. A Simplified Model of the Compensation Method for the Thermal Bias of a Ring Laser Gyro. Lasers Eng. 2014, 27, 119–126. [Google Scholar]
  22. Antonova, M.V.; Matveev, V.A. Model of error of A fiber-opyic gyro exposed to thermal and magnetic fields. Her. Bauman Mosc. State Tech. Univ. Ser. Instrum. Eng. 2014, 3, 73–80. [Google Scholar]
  23. Zha, F.; Xu, J.; Li, J.; He, H. IUKF neural network modeling for FOG temperature drift. J. Syst. Eng. Electron. 2013, 24, 838–844. [Google Scholar] [CrossRef]
  24. Prikhodko, I.; Trusov, A.; Shkel, A.M. Compensation of drifts in high-Q MEMS gyroscopes using temperature self-sensing. Sens. Actuators A: Phys. 2013, 201, 517–524. [Google Scholar] [CrossRef]
  25. Chen, X.; Shen, C. Study on error calibration of fiber optic gyroscope under intense ambient temperature variation. Appl. Opt. 2012, 51, 3755. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Wu, W.; Fang, Z.; Luo, B.; Li, Y.; Jiang, Q. Temperature Drift Compensation for Hemispherical Resonator Gyro Based on Natural Frequency. Sensors 2012, 12, 6434–6446. [Google Scholar] [CrossRef] [PubMed]
  27. Cai, Q.; Zhao, F.; Kang, Q.; Luo, Z.; Hu, D.; Liu, J.; Cao, H. A Novel Parallel Processing Model for Noise Reduction and Temperature Compensation of MEMS Gyroscope. Micromachines 2021, 12, 1285. [Google Scholar] [CrossRef] [PubMed]
  28. Zhou, Y.; Cao, H.; Guo, T. A Hybrid Algorithm for Noise Suppression of MEMS Accelerometer Based on the Improved VMD and TFPF. Micromachines 2022, 13, 891. [Google Scholar] [CrossRef]
  29. Zhou, B.; Zhang, T.; Yin, P.; Chen, Z.; Song, M.; Zhang, R. Innovation of flat gyro: Center Support Quadruple Mass Gyroscope. In Proceedings of the 2016 IEEE International Symposium on Inertial Sensors and Systems, Laguna Beach, CA, USA, 22–25 February 2016. [Google Scholar]
  30. Trusov, A.A.; Schofield, A.R.; Shkel, A.M. Gyroscope architecture with structurally forced anti-phase drive-mode and linearly coupled anti-phase sense-mode. In Proceedings of the TRANSDUCERS 2009-2009 International Solid-State Sensors, Actuators and Microsystems Conference, Denver, CO, USA, 21–25 June 2009. [Google Scholar]
Figure 1. (a) Structure of FMVMG, (b) The lumped model of FMVMG.
Figure 1. (a) Structure of FMVMG, (b) The lumped model of FMVMG.
Micromachines 14 00971 g001
Figure 2. (a) Driving mode, (b) Sensing mode.
Figure 2. (a) Driving mode, (b) Sensing mode.
Micromachines 14 00971 g002
Figure 3. FMVMG monitoring system.
Figure 3. FMVMG monitoring system.
Micromachines 14 00971 g003
Figure 4. FMVMG prototype.
Figure 4. FMVMG prototype.
Micromachines 14 00971 g004
Figure 5. Model of a radial basis function neural network.
Figure 5. Model of a radial basis function neural network.
Micromachines 14 00971 g005
Figure 6. The process of the RBF NN based on a GA (RBF NN+GA).
Figure 6. The process of the RBF NN based on a GA (RBF NN+GA).
Micromachines 14 00971 g006
Figure 7. The process of RBF NN+GA+KF.
Figure 7. The process of RBF NN+GA+KF.
Micromachines 14 00971 g007
Figure 8. The process of RBF NN+GA+KF denoising based on SE-EMD.
Figure 8. The process of RBF NN+GA+KF denoising based on SE-EMD.
Micromachines 14 00971 g008
Figure 9. (a) Temperature experiment equipment. (b) Inside the oven.
Figure 9. (a) Temperature experiment equipment. (b) Inside the oven.
Micromachines 14 00971 g009
Figure 10. Temperature experiment results in a temperature-controlled oven.
Figure 10. Temperature experiment results in a temperature-controlled oven.
Micromachines 14 00971 g010
Figure 11. IMFs decomposed by EMD (gyroscope output).
Figure 11. IMFs decomposed by EMD (gyroscope output).
Micromachines 14 00971 g011
Figure 12. Feature component extraction based on the SE.
Figure 12. Feature component extraction based on the SE.
Micromachines 14 00971 g012
Figure 13. IMFs decomposed by EMD (temperature).
Figure 13. IMFs decomposed by EMD (temperature).
Micromachines 14 00971 g013
Figure 14. (a) Mixed component after denoising, (b) The final reconstructed signal.
Figure 14. (a) Mixed component after denoising, (b) The final reconstructed signal.
Micromachines 14 00971 g014
Figure 15. Allan variance analysis of the gyroscope’s output after error processing.
Figure 15. Allan variance analysis of the gyroscope’s output after error processing.
Micromachines 14 00971 g015
Table 1. MEMS gyroscope structure mechanical values.
Table 1. MEMS gyroscope structure mechanical values.
ParameterValue
Elastic modulus (E)169 GPa
Poisson’s ratio (µ)0.27
Density (ρ)2328.3 kg/m³
Prototype length3015 μm
Prototype width2331 μm
Prototype height60 μm
Table 2. Allan variance analysis.
Table 2. Allan variance analysis.
DenoisingTemperature Compensation
Original DataRBF NNRBF NN+GARBF NN+GA+KF
B(°/h)N(°/h/Hz1/2)B(°/h)N(°/h/Hz1/2)B(°/h)N(°/h/Hz1/2)B(°/h)N(°/h/Hz1/2)
34.6699.60828.4299.603723.6599.60113.5890.967814
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Z.; Cui, Y.; Gu, Y.; Wang, G.; Yang, J.; Chen, K.; Cao, H. Temperature Drift Compensation for Four-Mass Vibration MEMS Gyroscope Based on EMD and Hybrid Filtering Fusion Method. Micromachines 2023, 14, 971. https://doi.org/10.3390/mi14050971

AMA Style

Li Z, Cui Y, Gu Y, Wang G, Yang J, Chen K, Cao H. Temperature Drift Compensation for Four-Mass Vibration MEMS Gyroscope Based on EMD and Hybrid Filtering Fusion Method. Micromachines. 2023; 14(5):971. https://doi.org/10.3390/mi14050971

Chicago/Turabian Style

Li, Zhong, Yuchen Cui, Yikuan Gu, Guodong Wang, Jian Yang, Kai Chen, and Huiliang Cao. 2023. "Temperature Drift Compensation for Four-Mass Vibration MEMS Gyroscope Based on EMD and Hybrid Filtering Fusion Method" Micromachines 14, no. 5: 971. https://doi.org/10.3390/mi14050971

APA Style

Li, Z., Cui, Y., Gu, Y., Wang, G., Yang, J., Chen, K., & Cao, H. (2023). Temperature Drift Compensation for Four-Mass Vibration MEMS Gyroscope Based on EMD and Hybrid Filtering Fusion Method. Micromachines, 14(5), 971. https://doi.org/10.3390/mi14050971

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