Next Article in Journal
Global-and-Local Context Network for Semantic Segmentation of Street View Images
Previous Article in Journal
C-MHAD: Continuous Multimodal Human Action Dataset of Simultaneous Video and Inertial Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros

1
College of Automation, Harbin Engineering University, Harbin 150001, China
2
Southwest China Research Institute of Electronic Equipment, Chengdu 610036, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(10), 2906; https://doi.org/10.3390/s20102906
Submission received: 18 April 2020 / Revised: 16 May 2020 / Accepted: 18 May 2020 / Published: 21 May 2020
(This article belongs to the Section Physical Sensors)

Abstract

:
The conventional temperature drift error (TDE) compensation model cannot decouple temperature dependence of Si-based materials because temperature correlated quantities (TCQ) have not been obtained comprehensively, and Micro-Electro-Mechanical System gyros’ (MEMS-gyros’) environmental adaptability is reduced in diverse, complicated conditions. The study presents modification of TDE compensation model of MEMS-gyros based on microstructure thermal effect analysis (MTEA). First, Si-based materials’ temperature dependence was studied in microstructure with thermal expansion effect and TCQ that determines the structural deformation were extracted to modify the conventional model, including temperature variation and its square. Second, a precise TDE test method was formed by analyzing heat conduction process between MEMS-gyros and thermal chamber, and temperature experiments were designed and conducted. Third, the modified model’s parameters were identified based on radical basis function artificial neural network (RBF ANN) and its performance was evaluated. Last, the conventional and modified models were compared in performance. The experimental results show MEMS-gyros’ bias stability was up to 10% of the conventional model, the temperature dependence of Si-based materials was decoupled better by the modified one and the environmental adaptability of MEMS-gyros was improved to expand their application in diverse complicated conditions.

1. Introduction

With the progress of science and technology, humans’ willingness to explore space and develop resources is increasing [1,2,3,4,5,6,7]. As is known, there are many rare resources in deep space, and there may even be an adaptive environment supporting life on the earth, and humans are attracted to go into action. Due to the extremely harsh environment in space, humans are unable to reach far. Instead, unmanned intelligent devices have been designed and used, such as all-weather Unmanned Aerial Vehicle (UAV) monitoring systems, lunar and Mars rovers, micro satellites, etc. [8,9,10,11,12]. Considering navigating or guiding unmanned intelligent devices to the target location safely, it is necessary to stabilize their attitudes. Gyros are indispensable to such navigation systems. Currently, unmanned intelligent devices have characteristics such as small size, low power consumption, incredible environmental adaptability—all of which require gyros to meet or even exceed stringent requirements. Base on this, resonator gyros have been invented and used. By measuring Coriolis force, carriers’ angular velocities are obtained, and their attitudes are measured. Table 1 shows the comparison between the resonator gyros.
Considering the accuracy, cost, size, stability and reliability comprehensively, Micro-Electro-Mechanical System gyros (MEMS-gyros) are easy to apply as well as precise to navigate and guide, and this is the best choice to unmanned intelligent devices. They are used to measure the attitude of unmanned intelligent devices’, which plays an important role in their stability and safety. For example, when unmanned rovers cruise on a planet with a complex terrain, measuring their stability in real time can predict if overturning and sinking is imminent, and alert them to take solutions to avoid the coming risk, such engaging an emergency brake or changing course. In addition, it can warn unmanned aerial vehicles of unsafe flight situations to avoid crash. However, MEMS-gyros are fabricated with temperature dependent Si-based materials, and their physical properties change as ambient temperature varies. Ambient temperature in space is about −180 °C–130 °C, and inevitably TDE come into play to reduce MEMS-gyros’ stability. Taking the accuracy of MEMS-gyro into account (±0.00875°/s) for example, when ambient temperature varies by 10 °C, its TDE is about 0.7°/s [13]. This results in attitude errors, velocity errors and heading errors accumulating over time, and some wrong references are given to unmanned intelligent devices. Maybe, they take some improper solutions to aggravate the dangerous situations. Hence, TDE restrict MEMS-gyros use in diverse complicated conditions; eliminating TDE plays a very important role in unmanned intelligent devices’ stability and safety. Decoupling temperature dependence of Si-based materials becomes the key to eliminate TDE and improve the environmental adaptability of MEMS-gyros [14,15].
To solve these problems, Masako Tanaka et al. analyzed the key factors affecting MEMS-gyros’ stability in depth and determined that structural consistency influenced driving frequency of driving circuit and combs stability of sensing circuit [16]. Liu et al. studied the effect of temperature variation on MEMS-gyros and especially the process that structural consistency changed the driving frequency [17]. They showed MEMS-gyros’ structural consistency lied in temperature dependence of Si-based materials in essence whose permanence is induced by constant ambient temperature. Hence, temperature control is optional to maintain structural consistency [18]. However, it needs high-power equipment to stabilize ambient temperature at the target, which conflicts with MEMS-gyros’ merits mentioned. Instead, digital compensation correction is an alternative, like Least Square Method (LSM), Kalman filter, artificial neural network. Jiaying Du et al. compared LSM, Kalman filter and digital filter in TDE compensation and studied their merits of their accuracy in real time [19]. TDE are estimated by them with ambient temperature. Kalman filter has highest accuracy, but worst real time due to its high dimensional matrixes and complex functions. To improve its real time further, Rita Fontanella et al. built a compensation model with augmented state Kalman filter. Its accuracy is maintained and its real time is improved by optimizing internal structure [20]. Markedly, Kalman filter depends on hardware resources to support its high accuracy and inevitably increases costs. Instead, LSM has the merits of estimating TDE accurately and quickly [21,22,23,24]. Igor P. Prikhodko et al. established LSM linear model with quality factor and ambient temperature [25]. Using quality factor, ambient temperature is estimated approximately and scale factor—as well as bias—is estimated to compensate TDE in real time. Scale factor and bias are up to 700 ppm and 2°/h. However, LSM is weakly able to restrain random error, so it needs a restraining method. Fuchao Liu and Hao Zheng et al. proposed an integrated model based on LSM and Kalman filter [26,27]. Kalman filter restrains random error in sample data and LSM estimates TDE with ambient temperature. Considering that Kalman filter only restrains random error, but does not calculate TDE, its real time is harmed, but in an allowable range. Improving the accuracy of LSM relies on increasing the fitting order and the higher the fitting order is the greater accuracy appears. However, it is unilateral that greater accuracy only relies on increasing the fitting order, which leads to over-compensation even real time reduction. Based on a number of test data, Bourgeteau et al. concluded complex nonlinearity appeared between ambient temperature and TDE, and describing it accurately was a prerequisite to TDE compensation accuracy. It also showed LSM could not increase TDE compensation accuracy further [28]. BP-ANN was introduced to describe the nonlinearity accurately [29,30,31]. It takes ambient temperature as the input and TDE as the output. By training BP-ANN, TDE compensation model is established after the compensation requirements are met. However, there may be local minimums which worsens its generalization ability and compensation accuracy in global scope. More important, TCQ exciting TDE should be figured out clearly and completely, which offers precise references to describe TDE. Jan K. Bekkeng et al. researched TCQ in depth and revealed that ambient temperature variation was a key factor to the structural deformation of MEMS-gyros [32]. Using ambient temperature variation, MEMS-gyros’ bias stability was improved and increased by more than one order of magnitude.
Improving the performance of TDE compensation is determined by three factors: TCQ exciting TDE, accurately describing the complex nonlinearity between TCQ and TDE, testing TDE precisely to identify the compensation model’s parameters. Hence, structural deformation of Si-based materials in MEMS-gyros is analyzed precisely with thermal expansion theory to extract TCQ comprehensively. With heat conduction analysis, TDE precise test method forms and its key parameters are deduced, including temperature jump interval and its period, and temperature experiments are conducted. A modified RBF ANN-based TDE compensation model is established and its parameters are identified precisely. It can estimate TDE more accurately and decouple Si-based materials’ temperature dependence effectively, and that increases MEMS-gyros’ stability and improves the environmental adaptability. It is significant to expand MEMS-gyros’ application in diverse complicated conditions and guarantees the safety and stability of unmanned intelligent devices and similar smart systems.
This article is organized as follows: in Section 2, TCQ of MEMS-gyros are extracted to form a modified TDE compensation model. Section 3 demonstrates the entire scheme for the modified model based on input-modified RBF ANN and shows the collected data utilized to train the input-modified RBF ANN and the experimental results of the MEMS-gyros-compensated temperature experiments. Section 4 evaluates TDE compensation performances of the proposed scheme compared to the previous conventional schemes. Section 5 presents the conclusions and benefits of the novel method.

2. Modification of TDE Compensation Model for MEMS-Gyros

2.1. Conventional TDE Compensation Model

MEMS-gyros are a kind of miniaturized devices manufactured with Si-based materials. They are mainly composed of the mass, the driving circuit, the sensing circuit and the substrate. With a series of procedure, including design, process, manufacture, measure and control, all of components are combined and integrated as a micromachining unit. Figure 1 shows the hardware design diagram of MEMS-gyros and their system schematic diagram [17].
Where k x and k y are stiffness coefficients of springs in driving direction x-axis and sensing direction y-axis, c x and c y are damping coefficients in x-axis and y-axis. With the help of the driving circuit, the mass m vibrates along x-axis under sinusoidal voltage with frequency ω d . When MEMS-gyros rotate along z-axis at angular velocity Ω , Coriolis force causes an displacement along y-axis. Technically, combs of sensing circuit can be abstracted as plate capacitors composed of moving plates and fixed plates. Any carrier’s angular velocity is obtained by measuring the capacitance variation in y-axis. Figure 2 shows the state diagram of combs before and after clockwise or counter clockwise rotation.
Based on plate capacitors’ definition, C 3 between 1# fixed plate and 2# fixed plate can be obtained:
C 3 = | C 1 C 2 | = | ε 4 π k S 0 ( d 0 + Δ d ) ε 4 π k S 0 ( d 0 Δ d ) | = 2 Δ C
Hence, angular velocity of the carrier can be obtained by measuring the capacitance variation. However, due to temperature dependence of Si-based materials, the stiffness of the driving circuit and the sensing circuit changes with ambient temperature. In particular, the stiffness is critical to the resonant frequencies in driving direction x-axis. According to elastic modulus formula, the elastic modulus variation of Si-based materials with ambient temperature can be described as follows:
E ( T ) = E ( T 0 ) [ 1 k ( T T 0 ) ]
where E ( T ) and E ( T 0 ) are the elastic modulus of Si-based materials at temperature T and T 0 separately, k is the elastic coefficient of Si-based materials and k = 7 0   ppm , T 0 = 300   K . Based on that, the spring stiffness coefficient is expressed as follows:
K ( T ) = K ( T 0 ) [ 1 k ( T T 0 ) ]
where K ( T ) and K ( T 0 ) are the stiffness coefficients at temperature T and T 0 . Hence, resonant frequency in x-axis ω x ( T ) can be induced as follows:
ω x ( T ) = K ( T 0 ) [ 1 k ( T T 0 ) ] / m x
where m x is the effective mass in driving mode. The resonant amplitude u max ( T ) and its phase α x ( T ) can be deduced as follows [17]:
u max ( T ) = F e m x { [ ω x ( T 0 ) ] 2 [ 1 k ( T T 0 ) / 2 ] 2 ω d 2 } 2 + { ω x ( T 0 ) [ 1 k ( T T 0 ) / 2 ] ω d / Q x } 2 α x ( T ) = arctan ω x ( T 0 ) [ 1 k ( T T 0 ) / 2 ] ω d Q x ω x ( T 0 ) 2 [ 1 k ( T T 0 ) / 2 ] 2 ω d 2
where F e is the amplitude of the external driving force in x-axis, Q x = k x / c x ω x . From Equation (5), the resonant amplitude and its phase are determined by ambient temperature T and its reference T 0 , which results in the resonant frequency instability as ambient temperature variation Δ T and Δ T = T 1 T 0 . Unavoidably, it introduces TDE in the output of MEMS-gyros. Hence, the conventional TDE compensation model considers ambient temperature variation as TCQ, and it is shown as follows:
Δ E M E M S = f ( Δ T )

2.2. Modified TDE Compensation Model

As it is known, although the stability of MEMS-gyros is directly determined by the resonant frequency in the driving circuit, it is not the only key factor to TDE. Furthermore, it cannot be ignored that the stiffness of the sensing circuit changes with ambient temperature as well, and the measuring error of the capacitance in sensing circuit should be reconsidered as another key factor. Taking the combs in sensing circuit for example, due to temperature dependence of Si-based materials, structural deformation of the combs appears in three-dimensional space as ambient temperature varies and the structural consistency changes. Assuming that ambient temperature is T and the angular velocity of carriers is ω , MEMS-gyros under diverse conditions were simulated as follows.
•  T = T 0 and ω = 0
When ambient temperature is T 0 constantly, the internal structure of MEMS-gyros is stable. When the carriers does not rotate and ω = 0 , moving plates and fixed plates are in a balanced state in ideal case. Figure 3 shows the combs of moving plates and fixed plates when ω = 0 .
Where a 0 is the thickness of the combs of fixed plates, b 0 is the length of the overlap between moving plates and fixed plates, c 0 is the width of the overlap between moving plates and fixed plates, d 0 is the distance between the combs of moving plates and fixed plates, e 0 is the thickness of the combs of moving plates. As shown in Figure 3, assuming that the combs are fabricated in an ideal case, then the capacitance measured in sensing circuit is as follows:
C 3 = | C 1 C 2 | = | ε 4 π k b 0 c 0 d 0 ε 4 π k b 0 c 0 d 0 | = 0
•  T = T 0 and ω = ω 0 ( ω 0 0 )
When ambient temperature stays still T 0 and the carrier rotates at the angular velocity of ω = ω 0 , Coriolis force acts on the mass m and a displacement appears in y-axis. Moving plates as well as fixed plates are in an unbalanced state, which changes the capacitance in the sensing circuit. Figure 4 shows the change of the combs of moving plates and fixed plates when ω = ω 0 .
From Figure 4, the combs of moving plates displace from the balanced state and according to Equation (1) the capacitance measured in the sensing circuit is expressed as follows:
C 3 = | C 1 C 2 | = | ε 4 π k b 0 c 0 d 0 Δ d ε 4 π k b 0 c 0 d 0 + Δ d | = 2 × ε b 0 c 0 4 π k Δ d ( d 0 Δ d ) ( d 0 + Δ d )
•  T = T 1 ( T 1 T 0 ) and ω = ω 0 ( ω 0 0 )
When T = T 1 , the combs deform in three dimensions as ambient temperature varies because of temperature dependence of Si-based materials, including expanding or contracting. Figure 5 shows the deformations before and after temperature variation when ω = 0 .
Where a 1 is the thickness of the combs of fixed plates after deformation, b 1 is the length of the overlap between moving plates and fixed plates after deformation, c 1 is the width of the overlap between moving plates and fixed plates after deformation, d 1 is the distance between the combs of moving plates and fixed plates after deformation, e 1 is the thickness of the combs of moving plates after deformation, Δ a is the thickness variation of the combs of fixed plates after deformation, Δ e is the thickness variation of the combs of moving plates after deformation. According to thermal expansion theory, three-dimensional sizes of the combs shown in Figure 5 after deforming are shown as follows:
{ a 1 = a 0 [ α T ( T 1 T 0 ) + 1 ] = a 0 ( α T Δ T + 1 ) c 1 = c 0 [ α T ( T 1 T 0 ) + 1 ] = c 0 ( α T Δ T + 1 ) e 1 = e 0 [ α T ( T 1 T 0 ) + 1 ] = e 0 ( α T Δ T + 1 )
where α T is the thermal expansion coefficient of Si-based materials, Δ T is ambient temperature variation and Δ T = T 1 T 0 . Under the excitation of TCQ, the combs of moving plates and fixed plates deform, respectively at the same time. In Figure 5b, assuming that the combs of moving plates deform in transverse left direction and its transverse length expands or contracts by Δ e , the combs of fixed plates deform in transverse right direction and its transverse length also expands or contracts by Δ a because fixed plates and moving plates have the same degrees of freedom for deformation. Hence, the total transverse length of the overlap of combs b 1 can be described as follows:
b 1 = 2 b 0 [ α T ( T 1 T 0 ) + 1 ] = 2 b 0 ( α T Δ T + 1 )
In addition, from Figure 5a, the transverse extension of the combs of fixed plates Δ a and the transverse extension of the combs of moving plates Δ e are, respectively shown as follows:
Δ a = a 1 a 0 = a 0 α T Δ T Δ e = e 1 e 0 = e 0 α T Δ T
Because the combs deform in transverse left direction and in transverse right direction, the distance between the combs of moving plates and fixed plates after deformation d 1 is expressed:
d 1 = d 0 Δ a 2 Δ e 2 = d 0 α T Δ T ( a 0 + e 0 ) 2
Given that the fabricated MEMS-gyros have firm structure, the mass has the same displacement under the same Coriolis force when the carriers rotate at the same angular velocity ω = ω 0 . That means the relative distance between the mass’ center when the carriers rotate at ω = 0 and ω = ω 0 still stay constant. Figure 6 shows the deformation of the combs when ω = ω 0 and T = T 1 .
According to Equation (1), the capacitance measured in the sensing circuit is shown as follows:
C 3 = | C 1 C 2 | = | ε 4 π k b 1 c 1 d 1 Δ d ε 4 π k b 1 c 1 d 0 + Δ d | = | 2 × ε b 1 c 1 4 π k Δ d ( d 1 Δ d ) ( d 1 + Δ d ) |
Then, from Equations (9), (10) and (12), we obtain:
C 3 = | 2 ε Δ d 2 b 0 c 0 4 π k × ( α T Δ T + 1 ) 2 d 1 2 Δ d 2 | = | 2 ε Δ d 2 b 0 c 0 4 π k × 4 ( α T Δ T + 1 ) 2 4 d 0 2 + α T 2 Δ T 2 ( a 0 + e 0 ) 2 4 d 0 α T Δ T ( a 0 + e 0 ) ( 2 Δ d ) 2 |
According to Equations (7) and (14), then the capacitance error Δ C E can be described:
Δ C E = C 3 C 3 = 2 ε b 0 c 0 Δ d 4 π k × 2 ( α T Δ T + 1 ) 2 d 1 2 Δ d 2 2 ε b 0 c 0 4 π k Δ d ( d 0 Δ d ) ( d 0 + Δ d ) = 2 ε Δ d b 0 c 0 4 π k ( d 0 2 Δ d 2 ) × α T 2 [ 8 d 0 2 ( a 0 + e 0 ) 2 8 Δ d 2 ] Δ T 2 + 4 α T [ 4 d 0 2 + ( a 0 + e 0 ) d 0 4 Δ d 2 ] Δ T + 4 d 0 2 4 Δ d 2 [ α T 2 ( a 0 + e 0 ) 2 Δ T 2 4 α T ( a 0 + e 0 ) d 0 Δ T + 4 d 0 2 4 Δ d 2 ]
From Equation (15), the capacitance error in sensing circuit is relevant to ambient temperature variation Δ T and its square Δ T 2 . Moreover, from Equation (6), ambient temperature variation is the key factor to the resonant frequency and important to TDE. Hence, Δ T and Δ T 2 are the critical references to compensate TDE accurately. Hence, the modified TDE compensation model is established:
Δ E M E M S = f ( Δ T , Δ T 2 )

3. Design of Modified TDE Compensation Model

Implementation of the modified TDE compensation model for MEMS-gyros is a priori process, and a mathematical model with perfect structure and clear parameters is gradually established by analyzing the sufficient experimental data which are tested and obtained beforehand. Hence, precise test for TDE is necessary and essential to the modified TDE compensation model, and the parameters identification are an important guarantee to its implementation.

3.1. Design of Precise Test for Temperature Drift Error

According to Equation (16), the modified TDE compensation model uses Δ T and Δ T 2 as the model inputs and TDE as the model output. Its accuracy depends on how precisely TDE can be accurately described with TCQ and that illustrates testing TDE and TCQ precisely is a prerequisite to the accuracy of the modified TDE compensation model. As it is known, TDE Δ E M E M S consists of bias error, trend error and random error and they are expressed as follows:
Δ E M E M S = E b i a s + E t r e n d + E r a n d o m
where E b i a s is bias error which is a type of fixed deviation between the measured value of the angular velocity of the carriers and its theoretical value, E t r e n d is trend error which is a type of linear or nonlinear deviation as ambient temperature varies, E r a n d o m is random error which is a series of small random fluctuations of some related factors and mutually compensated in the long term. E b i a s and E t r e n d account for most part of TDE, and it is possible to estimate and compensate accurately them due to their definite law and expression. When MEMS-gyros are manufactured, their environmental adaptability are unalterable. From the datasheets, their TDE can be grossly described as follows:
Δ E = α Δ T + β Δ T
where Δ E are roughly estimated values of TDE, α is MEMS-gyros’ character named “Zero-rate level change vs. temperature”, β is MEMS-gyros’ character named “Sensitivity change vs. temperature”, Δ T = T T 0 , and T 0 is the referenced ambient temperature. Theoretically, Δ E stands for part of TDE and is smaller than the true value of TDE in amplitude, so we can obtain as follows:
Δ E Δ E M E M S
Because MEMS-gyros’ sensitivity determines the measured minimum of the angular velocity of the carriers Δ E S , it is very possible that TDE appears greater than MEMS-gyros’ sensitivity when ambient temperature jumps rapidly, even completely submerging the actual angular velocity of the carriers. In that case, the amplitude of E r a n d o m also increases as TDE and greater than the actual angular velocity reference, then the deviation of the angular velocity reference is introduced. Based on that, it will make the angular velocity of the carrier be measured inaccurately and reduce the credibility of the angular velocity reference from the rate table. In order to test TDE accurately, TDE should be less than Δ E S , which is Δ E M E M S Δ E S . Hence, according to Equation (18), we obtain:
Δ T Δ E S | α | + | β |
Hence, a precise test for TDE should meet the equation given in Equation (20), which defines ambient temperature variation in precise test for TDE. In that case, MEMS-gyros are installed on a precise rate table in thermal chamber. At present, the thermal chambers adopt the structural design of front-door opening and temperature control unit (TCU) arranged on left and right sides, and a precise rate table locates in the center of the thermal chamber. Moreover, it adopts closed insulation design to prevent heat leakage to form high or low temperature condition. Figure 7 shows the schematic diagram of MEMS-gyros installed on the rate table inside the thermal chamber.
Considering the universality and the reliability of the test results, the following key operations should be paid more attention:
• Heat conduction measures
In order to control temperature gradient effect perfectly to improve the real time in the test, heat conduction measures should be taken to ensure ambient temperature in the thermal chamber is completely the same as that of MEMS-gyros, which reduces heat conduction delay effect.
• Precise temperature measurement system
To obtain TCQ more accurately, precise temperature measurement system is utilized, and temperature sensors are installed closely at the surface of MEMS-gyros. The measurement accuracy of precise temperature measurement system should be more than 2 times more precise than ambient temperature variation and its measurement frequency should be higher than the output frequency of MEMS-gyros to reserve the margin for the accuracy of the test results.
• Reasonable temperature control sequence
According to heat conduction theory, it takes some time to transfer heat from place A to place B . Assuming that temperature at place A varies in a sequence T A = [ T , T , T ] , after transferring in time sequence t = [ t 1 , t 2 ] , temperature at place B varies in a sequence T B = [ T , T , T ] . That points out it takes t 1 and t 2 to transfer the heat from place A to place B completely. During the period t 1 , temperature at place B is between T and T . At that time, if temperature at place A varies to T , temperature at place B will finally vary to T . However, temperature at place B never stays stable at T for a while. If that appears in precise test for TDE, TDE of MEMS-gyros is not tested accurately. Hence, a reasonable temperature control sequence that keeps transferring heat completely and stalely is an important guarantee for testing TDE accurately.
As shown in Figure 7, the internal space of the thermal chamber is artificially divided into two independent spaces with identical physical characteristics, 1#VR and 2#VR. TCUs control ambient temperature in the thermal chamber through its inner wall. Given that two independent spaces are cubes with size of L × L 1 × L 2 mm, heat from TCUs uniformly transfers to the joint of the independent spaces along length L which is perpendicular to the inner wall. The farther the location is away from the inner wall, the longer heat from TCUs transfer. The joint of independent spaces is the last area where ambient temperature stays stable. From Thermal Conductivity Formula, we can obtain:
k = Q t s L A Δ T
where Q is the conducted heat, t s is the time for heat conduction, L is the length of heat conduction, A is the section area of heat conduction, Δ T is ambient temperature variation. From the calculation formula of specific heat capacity, the heat heating the independent spaces uniformly can be expressed:
Q = C m Δ T
where C is specific heat capacity of air inside thermal chamber, m is the total mass of air in a closed state. Substituting Equation (22) into Equation (21), a new equation can be obtained as follows:
k = C m Δ T t s L A Δ T = C m t s L A
The time for heat conduction uniformly in the independent spaces can be deduced:
t s = C m k L A = C ρ v k L A = C ρ L A k L A = C ρ L 2 k
where ρ is air density in the thermal chamber. From Equation (24), the time for heat conduction from the inner wall of the thermal chamber to the center of the rate table can be calculated precisely. To guarantee that both of independent spaces are heated uniformly, the period t p when the current temperature control target goes to the next one can be expressed:
t s t p
Hence, Equation (20) is taken as a reference to temperature jump interval and Equation (25) is taken as a reference to temperature jump period. Based on the discussion, L3GD20 manufactured by ST company is tested for TDE. From the datasheet, Δ E S = 8.75   mdps / digit , α = ± 0.03   dps / ° C , its temperature range is 40   ° C ~ 85   ° C . After dimensional transformation, β = 2 % × F S / [ ( 85 ° C ) ( 40   ° C ) ] = 0.04   dps / ° C . Using Equation (20), temperature jump interval Δ T is found as 0.125   ° C . Given that accurately testing TDE and simplifying test steps, it is essential to set Δ T = 0.1   ° C . Using thermal chamber SET-Z-021 to test L3GD20, C = 1.005   kJ / ( kg × K ) , k = 0.0267   W / m ° C , L = 0.6   m , ρ = 1.293   kg / m 3 . From Equation (24), t s = 17.5209   s . The temperature period jump of 17.5209 s is taken for TCUs to vary temperature jump interval 0.1 °C, which ensures that the heat transfers from TCU to the center area uniformly and stably. To simplify the test steps, set t P = 20   s . L3GD20 is chosen randomly and tested five times, and temperature is monitored using precise temperature measurement system with the accuracy of ±0.03 °C as well as the frequency of 10 Hz [18]. Based on that, a temperature experiment was designed whose flow chart is shown in Figure 8 and its steps are described as follows:
1
MEMS-gyro was closely attached to a metal shell as a module with thermal silicone grease and installed on the precise rate table. The temperature sensor of precise temperature measurement system was attached closely to the metal shell and measures the temperature of MEMS-gyro T t 1 . Get PC ready to record the data from MEMS-gyro in real time;
2
Start the rate table at the target ω s , while get PC ready to record the data from MEMS-gyro D t 1 ;
3
Ambient temperature in thermal chamber goes down to −40 °C. After the data from MEMS-gyro and precise temperature measurement system are stable, start recording T t 1 and D t 1 ;
4
Ambient temperature in the thermal chamber goes up to 85 °C at a heating rate of 18 °C/h and stop the experiment when the data from MEMS-gyro and precise temperature measurement system keep stable for 1 h. Meanwhile, all the data during this period are recorded;
5
Repeat step (2) to step (4) 5 times and randomly select one group as the test data.

3.2. Parameter Identification for Modified TDE Compensation Model

Under the premise of establishing modified TDE compensation model and testing TDE and TCQ, the accuracy of compensating TDE depends on parameter identification based on TDE and TCQ. To test MEMS-gyros conveniently, set ω s = 0 and the output reference was 0°/s. Figure 9 shows one group of the experimental data of L3GD20 and its ambient temperature. Ambient temperature goes up from −40 °C to 85 °C, and its initial value keeps stable for a while which was set as the reference temperature. Temperature variation and its square as well as its TDE are also shown in Figure 9.
As shown in Figure 9, when ambient temperature T varies, ambient temperature variation Δ T and its square Δ T 2 vary in similar trend. With the help of Δ T and Δ T 2 , TDE have approximate trend. Based on that, it concludes that there is a complex nonlinearity among Δ T , Δ T 2 and TDE. In addition, on the basis of Equation (16), it is necessary and essential to apply a nonlinear model with multiple inputs and multiple outputs which has high accuracy and remarkable real time to fit the complex nonlinearity. RBF ANN uses neurons as the basic computing units and neural layers as the basic computing framework. Neurons are distributed in different neural layers, including an input layer, a hidden layer and an output layer. The inputs are calculated and transmitted by neurons in three layers and the kernel functions in hidden layers. Finally, the outputs approximately approach to the targets. Figure 10 shows the structure of RBF ANN.
Where X i ( i = 1 N ) is the i th input of RBF ANN, and Y i ( i = 1 M ) is the i th output of RBF ANN, and I i ( i = 1 N ) is the i th neuron in input layer, and H i ( i = 1 K ) is the i th neuron in hidden layer, and O i ( i = 1 M ) is the i th neuron in output layer. Usually, Gaussian function is chosen as kernel function and the inputs are divided into several groups by kernel function, which is described by:
ϕ j ( x ) = e x c j 2 / 2 σ j 2                         j = ( 1 , 2 , , K )
The most key point for RBF ANN is fixing the center c j and the width σ j of kernel function. As shown in Figure 10, input layer has N neurons and hidden layer has K neurons, as well output layer has M neurons. From Equation (26), ϕ j ( x ) is the output of the j th neuron in the hidden layer, c j is the center vector of kernel function of the j th neuron in the hidden layer, x is a N-dimensional input vector, σ j is the width of Gauss function of the j th neuron in the hidden layer, x c j is the distance between input vector and the center vector of Gauss function. Hence, the output of RBF ANN is shown:
y j = j = 1 K W i j ϕ j ( x )                         j = ( 1 , 2 , , M )
where y j is the output of the j th neuron in the output layer, W i j is the weight between the j th neuron in the output layer and the j th neuron in the hidden layer. The sample set is divided to several groups by RBF ANN with Equation (26), and the outputs of neurons in the hidden layer is obtained. Using Equation (27), RBF ANN calculates the actual output with the weights and the outputs of neurons in the hidden layer. Comparing the actual output with the targets, the difference decides if W i j , c j and σ j will be adjusted, and their adjusted magnitudes are shown as Δ W i j , Δ c j and Δ σ j respectively:
Δ W i j = η w ( d i y i ) ϕ j ( x ) Δ c j = η c i = 1 M [ ( d i y i ) W i j ] ( x c j ) σ j 2 ϕ j ( x ) Δ σ j = η σ i = 1 M [ ( d i y i ) W i j ] x c j 2 σ j 3 ϕ j ( x )
When RBF ANN is being trained, W i j , c j and σ j will be adjusted separately with Equation (28). After being adjusted for several times, the actual outputs of RBF ANN meet the design requirement [18]. In addition, there are two advantages about RBF ANN as follows:
  • RBF ANN can avoid local minimums. Owing that RBF ANN works on the basis of Gaussian functions, the current results are optimal in global scope, even in complex conditions like some flat areas where error gradient approximate to zero.
  • According to Kolmogorov theorem, a three-layer forward network can approach any continuous function with any desired accuracy [18]. RBF ANN has the typical structure of input layer, hidden layer and output layer and is able to realize the nonlinearity in any accuracy. In addition—considering real time and the universality—the structure of the modified TDE compensation models should be as simple as possible. Hence, RBF ANN is good at improving the real time and the universality of TDE compensation models for MEMS-gyros.
Therefore, RBF ANN is the best choice to describe the nonlinearity among Δ T , Δ T 2 and TDE accurately. Hence, Equation (16) can also be deduced as follows:
Δ E M E M S = A N N R B F ( Δ T , Δ T 2 )
The parameter of the modified TDE compensation model should be identified as follows:
1
Two temperature experiments are carried out. The experimental data from one temperature experiments group are in training sample set and the other data are in verification sample set.
2
The sample set of TDE is obtained by subtracting the reference outputs of MEMS-gyros from their actual outputs in the training sample set. The sample set of Δ T is obtained by subtracting the reference temperature of MEMS-gyros from their temperature in the training sample set, and Δ T 2 is obtained by multiplying itself.
3
RBF ANN is trained with Δ T and Δ T 2 as the inputs and TDE as the output. Training will not stop until the differences between the outputs of RBF ANN and the corresponding TDE meet the design requirements.
4
The compensated results are obtained from subtraction between the outputs of RBF ANN and the corresponding outputs of MEMS-gyros.
Based on all the steps above, Equation (29) will be trained with the experimental data shown in Figure 9, and the parameters will be identified accurately. Then, the modified TDE compensation model will be checked again with verification sample set. Figure 11 shows the primary outputs of MEMS-gyros shown in Figure 9 and its compensated outputs.
According to Figure 11, the modified model can estimate TDE accurately and MEMS-gyros run stably while ambient temperature varies from −40 °C to 85 °C. The experimental results means ambient temperature almost has no impact on MEMS-gyros. Usually, MEMS-gyros’ performance is evaluated with the indexes, bias stability, angle random walk, angle rate random walk, quantization noise and rate ramp. Bias stability not only shows the dispersion degree between MEMS-gyros’ output and its reference, but also can illustrate the dispersion degree between their outputs before and after TDE compensation. Hence, bias stability is applied to evaluate the accuracy of the modified TDE compensation model and the stability of MEMS-gyro after compensation, which is shown as follows:
B S = M S E ( x x )
where x is the evaluated sample, x is the reference of the evaluated sample, M S E is mean square error algorithm, and B S is bias stability between x and x . Bias stability is a remarkably intuitive index indicating the dispersion degree between the evaluated sample and its reference, which also shows the fluctuation of MEMS-gyro output after compensation. The smaller the bias stability, the smaller the dispersion degree between the evaluated sample and its reference is and the more accurately the modified model can estimate TDE. To illustrate the credibility and repetitiveness of the modified model, bias stabilities in five experiments before and after compensation are shown in Table 2.
From Table 2, the modified TDE compensation model effectively improve their stabilities and bias stabilities after compensation is about 4 orders of magnitude higher than before compensation. Hence, MEMS-gyro is unaffected by ambient temperature and the modified compensation model is able to decouple temperature dependence of Si-based materials remarkably.

4. Comparison of Test Results Before and After Modification

To verify the performance and universality of the modified TDE compensation model further, MEMS-gyro I3G4250D manufactured by ST company is selected instead as the test object. Based on temperature experiment and parameter identification, the conventional and the modified models of MEMS-gyros in x-axis, y-axis and z-axis were established, and their compensation performances were verified and compared. In order to guarantee the universality of test results, the referenced angular velocity in x-axis, y-axis, z-axis were randomly set as ω r e f x = 10   dps , ω r e f y = 5   dps , ω r e f z = 20   dps . Figure 12 shows the comparison of test results.
According to Equation (30), the evaluated formulas for the primary data, the data compensated by the conventional model and the data compensated by the modified model are shown as follows:
B S 1 = M S E ( k k r e f i ) B S 2 = M S E [ k ( Δ T ) k r e f i ] B S 3 = M S E [ k ( Δ T , Δ T 2 ) k r e f i ]
where k is the primary results of MEMS-gyros in x-axis, y-axis and z-axis; k ( Δ T ) is the test results compensated by the conventional model; k ( Δ T , Δ T 2 ) is the test results compensated by the modified model of MEMS-gyros; k r e f i ( i = x , y , z ) is the referenced angular velocity in x-axis, y-axis and z-axis. Bias stabilities in five experiments are shown as follows. To demonstrate the improvement between the primary data and the compensated results, performance improvement index is shown as below:
P i = B S i + 1 B S i         ( i = 1 , 2 )
According to Figure 12, the modified models estimate and compensate TDE more accurately, which keeps MEMS-gyros running stably for a long time while ambient temperature varies. From Table 3, Table 4, Table 5, Table 6 and Table 7, bias stabilities of the modified models were significantly smaller than the conventional model, which is increased to about 10% of bias stability of the conventional model. Hence, the modified TDE compensation model decouple temperature dependence of Si-based materials more remarkably, which achieves the purpose of improving the environmental adaptability.

5. Conclusions

In this study, a modification of an MTEA-based temperature drift error compensation model of MEMS-gyros was presented. Using microstructure thermal effect analysis, the novel TCQ (temperature variation and its square) were extracted. Then, two key parameters for TDE precise test method, temperature jump interval and its period, were deduced with heat conduction analysis. The modified TDE compensation models were built based on input-modified RBF ANN and their performances were verified and compared. The experimental results show MEMS-gyros run stably and accurately while ambient temperature varies and bias stability was increased by more than one order of magnitude. Temperature dependence of Si-based materials was decoupled completely and the environmental adaptability of MEMS-gyros was improved even in diverse complicated conditions.

Author Contributions

Conceptualization, B.Q., F.W. and F.L.; formal analysis, B.Q.; funding acquisition, J.C.; investigation, B.Q.; methodology, B.Q.; project administration, F.W.; supervision, F.W. and F.L.; validation, F.W.; visualization, J.C.; writing—original draft, B.Q.; writing—review & editing, F.W. and F.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant No.61633008, No.61773132 and No.61803115, in part by the 7th Generation Ultra Deep Water Drilling Unit Innovation Project sponsored by the Chinese Ministry of Industry and Information Technology, in part by the Heilongjiang Province Science Fund for Distinguished Young Scholars under Grant JC2018019 and in part by the Fundamental Research Funds for Central Universities under Grant HEUCFP201768 and 3072020CFT0403.

Acknowledgments

The authors thank the National Natural Science Foundation of China, the 7th Generation Ultra Deep Water Drilling Unit Innovation Project, the Heilongjiang Province Science Fund for Distinguished Young Scholars and the Fundamental Research Funds for Central Universities for financial support and also thank all the editors and anonymous reviewers for improving this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vetrella, A.R.; Fasano, G.; Accardo, D.; Moccia, A. Differential GNSS and vision-based tracking to improve navigation performance in cooperative multi-UAV systems. Sensors 2016, 16, 2164. [Google Scholar] [CrossRef] [PubMed]
  2. Renga, A.; Graziano, M.D.; D’Errico, M.; Moccia, A.; Menichino, F.; Vetrella, S.; Accardo, D.; Corraro, F.; Cuciniello, G.; Nebula, F. Galileo-based space-airborne bistatic SAR for UAS navigation. Aerosp. Sci. Technol. 2013, 27, 193–200. [Google Scholar] [CrossRef]
  3. Majdik, A.L.; Verda, D.; Albers-Schoenberg, Y.; Scaramuzza, D. Air-ground matching: Appearance-based GPS-denied urban localization of micro aerial vehicles. J. Robot. 2015, 32, 1015–1039. [Google Scholar] [CrossRef] [Green Version]
  4. Sanchez, J.; Denis, F.; Checchin, P.; Dupont, F.; Trassoudaine, L. Global Registration of 3D LiDAR Point Clouds Based on Scene Features: Application to Structured Environments. Remote. Sens. 2017, 9, 1014. [Google Scholar] [CrossRef] [Green Version]
  5. Hernández-Aceituno, J.; Arnay, R.; Toledo, J.; Acosta, L. Using Kinect on an Autonomous Vehicle for Outdoors Obstacle Detection. IEEE Sens. J. 2016, 16, 3603–3610. [Google Scholar] [CrossRef]
  6. Iuzzolino, M.; Accardo, D.; Rufino, G.; Oliva, E.; Tozzi, A.; Schipani, P. A cubesat payload for exoplanet detection. Sensors 2017, 17, 493. [Google Scholar] [CrossRef] [Green Version]
  7. Munguia, R.; Grau, A. A Practical Method for Implementing an Attitude and Heading Reference System. Int. J. Adv. Robot. Syst. 2014, 11, 62. [Google Scholar] [CrossRef] [Green Version]
  8. Chang, H.W.; Georgy, J.; El-Sheimy, N. Improved cycling navigation using inertial sensors measurements from portable devices with arbitrary orientation. IEEE Trans. Instrum. Meas. 2015, 64, 2012–2019. [Google Scholar] [CrossRef]
  9. Munguia, R.; Urzua, S.; Bolea, Y.; Grau, A. Vision-Based slam system for unmanned aerial vehicles. Sensors 2016, 16, 372. [Google Scholar] [CrossRef] [Green Version]
  10. Sheng, H.L.; Zhang, T.H. MEMS-based low-cost strap-down AHRS research. Eng. Meas. 2015, 59, 63–72. [Google Scholar] [CrossRef]
  11. Dan, L.H.; Ying, S.X.; Sheng, L.X. Application of Strongly Tracking Kalman Filter in MEMS Gyroscope Bias Compensation. In Proceedings of the 2017 6th International Conference on Advanced Materials and Computer Science, Zhengzhou, China, 29–30 April 2017; pp. 23–29. [Google Scholar]
  12. DeBruin, J. Control systems for Mobile Satcom antennas-Establishing and maintaining high-bandwidth satellite links during vehicle motion. IEEE Control Syst. Mag. 2008, 28, 86–101. [Google Scholar]
  13. Qi, B.; Cheng, J.H.; Zhao, L. A Novel Temperature Drift Error Model for MEMS Capacitive Accelerometer. In Proceedings of the 2017 IEEE 2nd Advanced Information Technology, Electronic and Automation Control Conference, Chongqing, China, 24–26 March 2017; pp. 182–186. [Google Scholar]
  14. Yang, C.; Li, H.S. Digital control system for the MEMS tuning fork gyroscope based on synchronous integral demodulator. IEEE Sens. J. 2015, 15, 5755–5764. [Google Scholar] [CrossRef]
  15. Zhang, L.; Masek, V.; Sanatdoost, N. Structural optimization of Z-axis tuning-fork MEMS gyroscopes for enhancing reliability and resolution. Microsyst. Technol. 2014, 21, 1187–1201. [Google Scholar] [CrossRef]
  16. Tanaka, M. An Overview of Quartz MEMS Devices. In Proceedings of the 2010 IEEE International Frequency Control Symposium, Newport Beach, CA, USA, 1–4 June 2010; pp. 162–167. [Google Scholar]
  17. Liu, D.C.; Chi, X.Z.; Cui, J.; Lin, L.T. Research on temperature dependent characteristics and compensation methods for digital gyroscope. In Proceedings of the 3rd International Conference on Sensing Technology, Taiwan, China, 30 November–3 December 2008; pp. 273–277. [Google Scholar]
  18. Cheng, J.H.; Qi, B.; Chen, D.D.; Landry, R.J. Modification of an RBF ANN-Based Temperature Compensation Model of Interferometric Fiber Optical Gyroscopes. Sensors 2015, 15, 11189–11207. [Google Scholar] [CrossRef] [Green Version]
  19. Du, J.Y.; Gerdtman, C.; Linden, M. Signal processing algorithms for temperature drift in a MEMS-gyro-based head mouse. In Proceedings of the 21ST International Conference on Systems, Signals and Image Processing, Dubrovnik, Croatia, 12–14 May 2014; pp. 123–126. [Google Scholar]
  20. Fontanella, R.; Accardo, D.; Lo Moriello, R.S.; Angrisani, L.; De Simone, 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] [Green Version]
  21. Gunhan, Y.; Unsal, D. Polynomial degree determination for temperature dependent error compensation of inertial sensors. In Proceedings of the 2014 IEEE/ION Position, Location and Navigation Symposium-PLANS 2014, Monterey, CA, USA, 5–8 May 2014; pp. 1209–1212. [Google Scholar]
  22. Zhang, B.Q.; Chu, H.R.; Sun, T.T.; Guo, L.H. Thermal calibration of a tri-axial MEMS gyroscope based on parameter-interpolation method. Sens. A Phys. 2017, 261, 103–116. [Google Scholar] [CrossRef]
  23. Chen, X.Y.; Xu, C.Y. Neural network modeling for FOG zero point drift based on forward linear prediction algorithm. J. Chin. Inert. Technol. 2007, 15, 334–337. [Google Scholar]
  24. Doostdar, P.; Keighobadi, J. Design and implementation of SMO for a nonlinear MIMO AHRS. Mech. Syst. Process. 2012, 32, 94–115. [Google Scholar] [CrossRef]
  25. Prikhodko, I.P.; Trusov, A.A.; Shkel, A.M. Compensation of drifts in high-Q MEMS gyroscopes using temperature self-sensing. Sens. A Phys. 2013, 201, 517–524. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, F.C.; Su, Z.; Li, Q.; Li, C.; Zhao, H. Error Characteristics and Compensation Methods of MIMU with Non-centroid Configurations. In Proceedings of the 37th Chinese Control Conference, Wuhan, China, 25–27 July 2018; pp. 4877–4882. [Google Scholar]
  27. Zheng, H.; Yuan, D.; Yang, J.B.; Chen, F. Hardware Design of Virtual Gyroscope Based on Array Technology. In Proceedings of the 2017 2nd International Conference on Frontiers of Sensors Technologies, Shenzhen, China, 14–16 April 2017; pp. 24–27. [Google Scholar]
  28. Bourgeteau, B.; Levy, R.; Lavenus, P.; Le Traon, O. Multiphysical Finite Element Modeling of a Quartz Micro-Resonator Thermal Sensitivity. In Proceedings of the 2015 Symposium on Design Test Integration and Packaging of MEMS and MOEMS, Montpellier, France, 27–30 April 2015. [Google Scholar]
  29. Fontanella, R.; Accardo, D.; Lo Moriello, R.S. MEMS gyros temperature calibration through artificial neural networks. Sens. A Phys. 2018, 279, 553–565. [Google Scholar] [CrossRef]
  30. Xia, D.Z.; Chen, S.L.; Wang, S.R.; Li, H.S. Microgyroscope temperature effects and compensation-control methods. Sensors 2009, 9, 8349–8376. [Google Scholar] [CrossRef] [PubMed]
  31. Fontanella, R.; Accardo, D.; Lo Moriello, R.S. An Extensive Analysis for the use of Back Propagation Neural Networks to Perform the Calibration of MEMS Gyro Bias Thermal Drift. In Proceedings of the 2016 IEEE/ION Position, Location and Navigation Symposium (PLANS), Savannah, GA, USA, 11–14 April 2016; pp. 672–680. [Google Scholar]
  32. Bekkeng, J.K. Calibration of a Novel MEMS Inertial Reference Unit. IEEE Trans. Instrum. Meas. 2009, 58, 1967–1974. [Google Scholar] [CrossRef]
Figure 1. Functional schematic diagram of MEMS-gyros. (a) hardware design; (b) system schematic.
Figure 1. Functional schematic diagram of MEMS-gyros. (a) hardware design; (b) system schematic.
Sensors 20 02906 g001
Figure 2. State diagram of the combs before and after rotation. (a) clockwise; (b) counter-clockwise.
Figure 2. State diagram of the combs before and after rotation. (a) clockwise; (b) counter-clockwise.
Sensors 20 02906 g002
Figure 3. Combs of moving plates and fixed plates when T = T 0 and ω = 0 . (a) top view; (b) side view.
Figure 3. Combs of moving plates and fixed plates when T = T 0 and ω = 0 . (a) top view; (b) side view.
Sensors 20 02906 g003
Figure 4. Combs of moving plates and fixed plates when T = T 0 and ω = ω 0 . (a) top view; (b) side view.
Figure 4. Combs of moving plates and fixed plates when T = T 0 and ω = ω 0 . (a) top view; (b) side view.
Sensors 20 02906 g004
Figure 5. Deformation of the combs when T = T 1 and ω = 0 . (a) top view; (b) side view.
Figure 5. Deformation of the combs when T = T 1 and ω = 0 . (a) top view; (b) side view.
Sensors 20 02906 g005
Figure 6. Deformation of the combs when T = T 1 and ω = ω 0 . (a) top view; (b) side view.
Figure 6. Deformation of the combs when T = T 1 and ω = ω 0 . (a) top view; (b) side view.
Sensors 20 02906 g006
Figure 7. Schematic diagram of MEMS-gyros installed on the rate table inside the thermal chamber.
Figure 7. Schematic diagram of MEMS-gyros installed on the rate table inside the thermal chamber.
Sensors 20 02906 g007
Figure 8. Flow chart of temperature experiment.
Figure 8. Flow chart of temperature experiment.
Sensors 20 02906 g008
Figure 9. Experimental results of L3GD20. (a) ambient temperature and temperature variation as well as its square; (b) outputs of L3GD20 and its temperature drift error (TDE).
Figure 9. Experimental results of L3GD20. (a) ambient temperature and temperature variation as well as its square; (b) outputs of L3GD20 and its temperature drift error (TDE).
Sensors 20 02906 g009
Figure 10. Structure of a radical basis function artificial neural network (RBF ANN).
Figure 10. Structure of a radical basis function artificial neural network (RBF ANN).
Sensors 20 02906 g010
Figure 11. Primary outputs of MEMS-gyros and its compensated outputs.
Figure 11. Primary outputs of MEMS-gyros and its compensated outputs.
Sensors 20 02906 g011
Figure 12. Comparison of the test results in five experiments. (a) 1st experiment; (b) 2nd experiment; (c) 3rd experiment; (d) 4th experiment; (e) 5th experiment.
Figure 12. Comparison of the test results in five experiments. (a) 1st experiment; (b) 2nd experiment; (c) 3rd experiment; (d) 4th experiment; (e) 5th experiment.
Sensors 20 02906 g012aSensors 20 02906 g012b
Table 1. Comparison between resonator gyros.
Table 1. Comparison between resonator gyros.
.PrincipleMeritsDemerits
Tuning fork
gyro
Coriolis effect and
momentum conservation
Low cost
Simple implementation
Low precision
Piezoelectric
gyro
Coriolis effect and
Piezoelectric effect
Long service life
High reliability
Worse environmental adaptability
MEMS-gyroCoriolis effect and
capacitor
Low cost and small size
Higher precision
Bad environmental adaptability
HR-gyroCoriolis effect and
precession effect of standing wave in radial vibration
Highest precision
High bandwidth
Strong overload
High cost
Complex implementation
Table 2. Bias stabilities before and after compensation.
Table 2. Bias stabilities before and after compensation.
BS1
(The Primary Data)
BS2
(After Compensation)
Improvement
(BS2/BS1)
1st experiment30.04027.882 × 10−32.6239 × 10−4
2nd experiment29.38375.641 × 10−31.9199 × 10−4
3rd experiment28.38165.615 × 10−31.9782 × 10−4
4th experiment29.25485.630 × 10−31.9244 × 10−4
5th experiment29.42015.758 × 10−31.9570 × 10−4
Table 3. Bias stabilities of the experimental data in 1st experiment.
Table 3. Bias stabilities of the experimental data in 1st experiment.
BS1BS2BS3P1P2P3 = P2/P1
x-axis15.37110.00553.885 × 1043.594 × 1042.528 × 1057.03%
y-axis29.26210.00776.201 × 1042.617 × 10−42.119 × 10−58.09%
z-axis1.53360.00754.571 × 1044.889 × 1032.980 × 1046.09%
Table 4. Bias stabilities of the experimental data in 2nd experiment.
Table 4. Bias stabilities of the experimental data in 2nd experiment.
BS1BS2BS3P1P2P3 = P2/P1
x-axis15.06190.01671.716 × 1031.111 × 1031.139 × 10410.26%
y-axis28.69820.01881.879 × 1036.558 × 10−46.549 × 10−59.99%
z-axis1.59750.01771.807 × 1031.109 × 1021.131 × 10310.19%
Table 5. Bias stabilities of the experimental data in 3rd experiment.
Table 5. Bias stabilities of the experimental data in 3rd experiment.
BS1BS2BS3P1P2P3 = P2/P1
x-axis15.70150.01741.654 × 1031.109 × 1031.054 × 1049.49%
y-axis29.10820.01801.737 × 1036.183 × 10−45.969 × 10−59.65%
z-axis1.54720.01841.904 × 1031.190 × 1021.231 × 10310.34%
Table 6. Bias stabilities of the experimental data in 4th experiment.
Table 6. Bias stabilities of the experimental data in 4th experiment.
BS1BS2BS3P1P2P3 = P2/P1
x-axis15.65230.01791.870 × 1031.146 × 1031.195 × 10410.42%
y-axis29.01140.01871.813 × 1036.449 × 10−46.249 × 10−59.69%
z-axis1.40890.01801.740 × 1031.275 × 1021.235 × 1039.69%
Table 7. Bias stabilities of the experimental data in 5th experiment.
Table 7. Bias stabilities of the experimental data in 5th experiment.
BS1BS2BS3P1P2P3 = P2/P1
x-axis15.37460.01721.749 × 1031.117 × 1031.138 × 10410.18%
y-axis29.43750.01861.757 × 1036.319 × 10−45.970 × 10−59.45%
z-axis1.58110.01852.038 × 1031.171 × 1021.289 × 10311.01%

Share and Cite

MDPI and ACS Style

Qi, B.; Wen, F.; Liu, F.; Cheng, J. Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros. Sensors 2020, 20, 2906. https://doi.org/10.3390/s20102906

AMA Style

Qi B, Wen F, Liu F, Cheng J. Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros. Sensors. 2020; 20(10):2906. https://doi.org/10.3390/s20102906

Chicago/Turabian Style

Qi, Bing, Fuzhong Wen, Fanming Liu, and Jianhua Cheng. 2020. "Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros" Sensors 20, no. 10: 2906. https://doi.org/10.3390/s20102906

APA Style

Qi, B., Wen, F., Liu, F., & Cheng, J. (2020). Modification of MTEA-Based Temperature Drift Error Compensation Model for MEMS-Gyros. Sensors, 20(10), 2906. https://doi.org/10.3390/s20102906

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