Next Article in Journal
Supervisory Algorithm for Autonomous Hemodynamic Management Systems
Next Article in Special Issue
Fractional Derivatives Application to Image Fusion Problems
Previous Article in Journal
Research and Development of a COVID-19 Tracking System in Order to Implement Analytical Tools to Reduce the Infection Risk
Previous Article in Special Issue
Triple Estimation of Fractional Variable Order, Parameters, and State Variables Based on the Unscented Fractional Order Kalman Filter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MEMS Accelerometer Noises Analysis Based on Triple Estimation Fractional Order Algorithm

Institute of Control and Industrial Electronics, Warsaw University of Technology, ul. Koszykowa 75, 00-662 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(2), 527; https://doi.org/10.3390/s22020527
Submission received: 12 December 2021 / Revised: 4 January 2022 / Accepted: 7 January 2022 / Published: 11 January 2022
(This article belongs to the Special Issue Fractional Sensor Fusion and Its Applications)

Abstract

:
This paper is devoted to identifying parameters of fractional order noises with application to noises obtained from MEMS accelerometer. The analysis and parameters estimation will be based on the Triple Estimation algorithm, which can simultaneously estimate state, fractional order, and parameter estimates. The capability of the Triple Estimation algorithm to fractional noises estimation will be confirmed by the sets of numerical analyses for fractional constant and variable order systems with Gaussian noise input signal. For experimental data analysis, the MEMS sensor SparkFun MPU9250 Inertial Measurement Unit (IMU) was used with data obtained from the accelerometer in x, y and z-axes. The experimental results clearly show the existence of fractional noise in this MEMS’ noise, which can be essential information in the design of filtering algorithms, for example, in inertial navigation.

1. Introduction

Micromachined Electrical Mechanical Systems (MEMS) are mechanical and electro-mechanical devices made using microfabrication techniques. MEMS technology allows to build miniature (e.g., inside integrated circuits) sensor and actuators which have been employed in many different areas such as medicine [1], biology [2], chemistry [3], aerospace [4], and motor vehicles [5]. A crucial area in which MEMS sensors are used is inertial navigation systems (INS) [6,7] based on double integration of body acceleration processes based on accelerometers and gyroscopes measurements. Due to the double integration action, high accuracy and precision of acceleration measurement are essential because noises (especially biases) are double-integrated and rapidly increase navigation errors. That is why modelling noises, biases, and general dynamics of MEMS sensors are essential. An article [8] uses, for example, an advanced type of recurrent neural network to model some parts of non-modeled MEMS gyroscope dynamics and apply this network into fractional order sliding mode control. In [9], analysis of noises in MEMS sensors is presented. The article [10] presents a mathematical modelling MEMS sensors dynamics, including modelling the noises. Presented in the mentioned work, the model is very complex and underlines the impact of temperature noises. The impact of thermal noise is also presented in [11]. What is important, as it will be discussed later, is that diffusive systems, in particular heat transfer processes, can be very efficiently modelled using fractional order calculus. Moreover, in [12], a fractional order algorithm was effectively used to estimate the bias of the MEMS sensor. That is why in this paper, we will use a fractional order estimation algorithm to identify noises of MEMS’ sensor as a fractional order noise.
The fractional calculus (FC) itself, is an extension of traditional differential and integral calculus. The differential orders in FC can be real or even complex numbers. The fractional derivative appeared for the first time in the correspondence between Leibniz and l’Hôpital in 1695, and thereby, it appeared almost simultaneously with the integer order calculus. The theoretical background for this calculus can be found in already classic works such as [13,14,15,16,17,18], as well as (with some applications) in relatively recently published books [19,20,21].
In contradiction to integer order derivatives, the fractional order derivatives depend not only on local time conditions but also on the whole past of the function [22]. This specific property has an advantage when the analysed dynamics possess a long-term memory nature, and thus, past values of the function are taken into account. The FC was found to be especially efficient in modelling diffusive systems [23,24,25,26]. For example, in the heat transfer process of the solid beam, it is possible to describe dynamics between temperature and heat flux at the desired point as a half order integral. When the heated material is not solid (heterogeneous), the order of the integration can be different by half, as it was presented in [23].
The FC was also recognised as an efficient tool in signal processing to design new types of filters and new tools for signal analysis. Some applications of fractional order calculus to signal processing were presented in [17,27,28,29].
It may happen that using constant order operators applied for some complex dynamic processes may be unsatisfactory, for example, diffusion processes in the porous (inhomogeneous or heterogeneous) environment, where the structure may vary in time [24]. In such a case, the fractional variable order (FVO) operators have to be used. To give a deeper insight into FVO calculus, four switching schemes, equivalent to four definitions of FVO derivatives, are presented in [30,31,32]. The switching strategies introduced, given unambiguously, classify and identify ways of changing the order of derivatives (integrals). Based on those switching schemes, it is possible to categorise fractional order derivatives according to their behaviour and intrinsic properties.
For FVO systems, it was also a generalised Kalman Filter obtaining Variable Order Kalman Filter [33] applied for estimation of fractional variable order state-space systems over a lossy network. A generalisation of the Improved Fractional Kalman Filter (ExFKF) for variable order discrete state-space systems is presented in [34], where the results are used for estimation and smoothing actions.
When the uncorrelated noise (like white noise) passes through a dynamical system, the dynamically correlated noise (coloured noise) is obtained. When the dynamics contain fractional order, the fractional noise is obtained. Article [35] presents an introduction to fractional order noises (the noises obtained by applying uncorrelated white noise to fractional order dynamics).
In [36], fractional signal processing methods were used to electrochemical noise of corrosion processes in stainless steel. It has been shown there that conventional (integer order) methods failed to sufficiently distinguish between electrochemical noise signals obtained from different solutions, and the use of fractional Fourier transforms turned out to be a powerful tool that can better describe the corrosion processes of the electrode.
Article [37] presents estimation schemes for discrete fractional and integer order state-space systems with fractional order coloured noise. Better estimates of the state vector were obtained there thanks to the additional information about noise dynamics used by the proposed estimation algorithm.
Article [38] presents an adaptive filtering approach to filter the noise from MEMS measurements, where an adaptive Kalman filter was derived from the integer order dynamic noise model.
In [39], modelling results of temperature sensor noise were presented. The identification algorithm is based on numerical minimisation of dynamical correlation of identified source noise with the Least Square algorithm. The results presented show that the noise order of the temperature sensor is fractional. The used algorithm assumed a situation in which the only evident fractional order noise is analysed. However, in the real plant, we would rather have a situation when the measured noise is a combination of dynamically correlated and uncorrelated noises.
Interesting results showing the source of noises in MEMS, both of a mechanical and electrical nature, and in particular the impact of thermal-noise on bifurcation MEMS sensors, were presented in [11]. The authors proposed a stochastic integer order model of the sensor that takes these noise sources into account.
In [40], the Triple Estimation algorithm (TEA) for state vector, order, and system’s parameters estimation was proposed and described in detail. In this paper, we will present the TEA’s application to estimate noise parameters for a case when noise is a combination of correlated and uncorrelated dynamically noises. Moreover, we will also present the analysis results of a real plant noise obtained from MEMS sensor.
In [12], Fractional Kalman algorithm in two versions has been used to improve measurement results from MEMS sensors. The results of MEMS noises modelling obtained in our article could explain why Fractional Order Kalman Filter was more efficient in this case.
Thus, the main novelty of this paper is identifying and analysing the accelerometer noises built-in MEMS technology. We have conducted sets of experiments based on real plant data to show the fractional dynamic of the investigated sensor. The paper also describes the identification process of fractional order noises and problems connected with this action. The obtained results present the ability of the TEA to model the noises from real plants. All of these issues make a new contribution to this research topic.
The paper is organised as follows: Section 2 recalls elements of FC and fractional noises. In Section 3, the TEA is presented. Section 4 describes a method for application TEA to fractional noise identification. Finally, Section 5 presents results of using the TEA to model noises obtained from MEMS accelerometer.

2. Fractional Calculus and Fractional Noises

In FC, the three most popular definitions of fractional constant order integral and derivative are used, namely, Grünwald–Letnikov, Riemann–Liouville, and Caputo. These definitions possess different properties and may be applied in various areas of engineering.
In this paper, we use the Grünwald–Letnikov definition, which is usually used in discrete systems, as a base for FVO difference definition. Due to the application nature of this work, we will use a discrete approximation of the Grünwald–Letnikov derivative with a finite (not going to zero) sampling time h. Hence, we have constant order difference definition
Δ k α 0 x k j = 0 k 1 h α ( 1 ) j α j x k j ,
where
α j 1 for j = 0 , α ( α 1 ) ( α j + 1 ) j ! for j > 0 ,
α R is a fractional order and h is a time sampling.
In our paper, we will use the following FVO type of difference:
Δ k α k 0 A x k j = 0 k ( 1 ) j h α k α k j x k j ,
where α k R is FVO.

Fractional Noise

The time-correlated (coloured) noises are the noises that contain a dynamical correlation between the noise’s samples. Such noises can be obtained when some noise (uncorrelated) is passed through dynamical systems. For example, the electromagnetic field noise can induct some current in an electronic circuit, leading to some dynamically correlated noise in voltage because of some dynamic between current and voltage. When the order of the dynamics is an integer, we will have a dynamically correlated integer order noise, which the following relation can describe:
x k + 1 = f x k + ω k ,
where x k is a time-correlated noise, and ω k is an uncorrelated noise for example white Gaussian noise.
When the dynamics of the system are fractional, for example, in temperature transport (for ideal beam temperature is half order integral of heat flux [23]), the uncorrelated heat flux noise can lead to fractional order dynamically correlated noise in temperature. The coloured fractional order noise is given as follows
Δ k + 1 α 0 x k + 1 = f x k + ω k
x k + 1 = h α Δ k + 1 α 0 x k + 1 j = 1 k + 1 ( 1 ) j α j x k j + 1 ,
where x k is a fractional coloured noise, α is an order of the noise, and ω k is an uncorrelated noise.
For the case when the fractional order of the dynamical system changes in time (for example, when the structure of heated medium changes in time [24]), the FVO noise will appear. Depending on the order switching manner, different definitions can describe such dynamics. For example, for A -type definition, we will have the following FVO noise dynamics:
Δ k α k + 1 0 A x k + 1 = f x k + ω k
x k + 1 = h α k + 1 Δ k α k + 1 0 A x k + 1 j = 1 k + 1 ( 1 ) j α k + 1 j x k j + 1 .
Identification of the fractional noise in a real application is a complex process because we do not know the order and system parameters of the noise. We also do not have information about dynamically uncorrelated source noise. In [39], identification algorithm for fractional noise was presented, however, under the assumption that output noise is the only evident fractional order noise. In experimentally obtained noises, we would instead acquire a combination of dynamically correlated and uncorrelated noises
y k = x k + ν k .
That is why in this article, we use a Triple Estimation algorithm to identify parameters of fractional order noises.

3. Triple Estimation Algorithm

The TEA allows estimating state vector, system parameters, and fractional order simultaneously. The main idea of this algorithm is to separate states, parameters, and orders estimation processes. This separation allows a better adjustment of used filters, making it possible to obtain better estimation results. Detailed introduction of TEA was presented in [40].
The TEA will be defined for the following linear Discrete Fractional Variable Order State-Space (DFVOSS) A -type system [41]:
Δ k + 1 α k + 1 0 A x k + 1 = A x k + B u k + ω k ,
x k + 1 = h α k + 1 Δ k + 1 α k + 1 0 A x k + 1 j = 1 k + 1 ( 1 ) j α k + 1 j x k j + 1 ,
y k = C x k + ν k ,
where u k R d is a system input; y k R p is a system output; A R N × N , B R N × d , and C R p × N are the state system, input, and output matrices, respectively; x k R N is a state vector; N is a number of state equations.
In general, the TEA can be treated as a method for simultaneous states, parameters, and order estimation for fractional order systems. Moreover, separation of order and system parameters estimation processes allows better algorithm parameters tuning because we can separately tune parameters for order and system parameters filters.
In the TEA process, the FVO, state variables and parameters estimation is divided into three estimation actions (filters). The first filter, KFx, estimates the state variables vector x ^ k based on order and system parameters estimates from other filters KFo and KFw, respectively. The second one, KFw, estimates the vector of system parameters w ^ k based on state variable and order estimates obtained in the remaining two filters KFx and KFo, respectively. The third filter, KFo, estimates the FVO with the knowledge of state variable and system parameters from filters KFx and KFw, respectively. The scheme of the TEA is given in Figure 1.

3.1. Order Estimation Filter KFo

Because the order estimation problem is highly non-linear (due to ( α k , i j ) ) relations in obtaining state update process), as the KFo filter, the Unscented Fractional Variable Order Kalman Filter is used. The order changing dynamics is assumed to be a constant
α k + 1 = α k + ω k o ,
where ω k o is a noise with variance given by matrix Q k o . The matrix Q k o represents our knowledge of how big fluctuations in time we are assumed. The bigger the value of this matrix, the more the algorithm will spread estimation error to modify the order.
The KFo algorithm equations are given as follows:
α ˜ k = α ^ k 1 ,
P ˜ k o = P ^ k 1 o + Q k 1 o ,
α ˜ k = α ˜ k α ˜ k ± ( L + λ ) P ˜ k o i ,
Δ α ˜ k , i χ ˜ k , i o = A ( w ^ k 1 ) x ^ k 1 + B u k 1 ,
χ ˜ k , i o = h α ˜ k , i Δ α ˜ k , i χ ˜ k , i o j = 1 k ( 1 ) j α ˜ k , i j x ^ k j ,
Y ˜ k , i o = C χ ˜ k , i o ,
y ˜ k o = i = 0 2 L W ( m ) Y ˜ k , i ,
P y k y k o = i = 1 2 L W i ( c ) [ Y ˜ i , k y ˜ k ] [ Y ˜ i , k y ˜ k ] T + R o ,
P α k y k o = i = 1 2 L W i ( c ) [ α ˜ i , k α ˜ k ] [ Y ˜ i , k y ˜ k ] T ,
K k o = P α k y k o ( P y k y k o ) 1 ,
α ^ k = α ˜ k + K k o ( y k y ˜ k o ) ,
P k o = P ^ k o K k o P y k y k o K k o ,
Q k o = ( 1 δ o ) Q k 1 o + δ o ( K k o ) ( y k y ˜ k o ) ( y k y ˜ k o ) T ( K k o ) T ,
where ( ( L + λ ) P k ) i is i-th column of matrix square root (e.g., Cholesky factorisation), L is a dimension of estimated state vector ( 2 L + 1 is a number of sigma points) and coefficients of Unscented transformation W are equal to
W 0 ( m ) = λ / ( L + λ ) ,
W 0 ( c ) = λ / ( L + λ ) + ( 1 A 2 + B ) ,
W i ( m ) = W i ( c ) = 1 / ( 2 ( L + λ ) ) ,
where λ = A 2 ( L + κ ) L , A is a coefficient describing width of point expansion during the transformation (in literature is obtained in the range 1 A 1 e 4 , usually denoted as α , but in this article, because of using order α this notation has been changed); κ is an additional scaling coefficient usually chosen as 3-L; B is a coefficient that corresponds with our knowledge about type of noise, for Gaussian noise is chosen as B = 2 (in literature usually denoted as β ). The δ coefficient is a “forgetting factor” according to Robbins–Monro stochastic approximation scheme for estimating the innovations (see [42] p. 240). The initial values of matrix P 0 o represent our a’priori knowledge about error in choosing initial value of order α 0 (we assume, the initial value is different from the original).

3.2. State Estimation Filter KFx

As the KFx Filter, the Fractional Variable Order Kalman Filter algorithm is used because of the linearity of the state estimation sub-process. This filter is given as follows:
Δ k + 1 α ^ k 0 A x ˜ k + 1 = A ( w ^ k 1 ) x ^ k + B u k ,
x ˜ k + 1 = h α ^ k Δ k + 1 α ^ k 0 A x ˜ k + 1 j = 1 k + 1 ( 1 ) j α ^ k j x ^ k + 1 j ,
P ˜ k = h α ^ k A ( w ^ k 1 ) + α ^ k P k 1 h α ^ k A ( w ^ k 1 ) + α ^ k T
+ Q k 1 + j = 2 k α ^ k j P k j α ^ k j T ,
K k = P ˜ k C T ( C P ˜ k C T + R k ) 1 ,
x ^ k = x ˜ k + K k ( y k C x ˜ k ) ,
P k = ( I K k C ) P ˜ k ,
where initial conditions are
x 0 R N , P 0 = E [ ( x ˜ 0 x 0 ) ( x ˜ 0 x 0 ) T ] ,
and ν k and ω k are assumed to be independent with zero expected value.

3.3. Parameters Estimation Filter KFw

For KFw filter another Unscented Fractional Variable Order Kalman Filter is used. The dynamics of parameter changing is also assumed as constant
w k + 1 = w k + ω k w ,
where ω k w is a noise with variance given by matrix Q k w . The equations of the filter KFw are very similar to filter KFo, and the difference is only in the model replica part:
w ˜ k = w ^ k 1 ,
P ˜ k w = P ^ k 1 w + Q k 1 w ,
W ˜ k = w ˜ k w ˜ k ± ( L + λ ) P ˜ k w i ,
Δ α ^ k 1 χ ˜ k , i w = A ( W ˜ k , i ) x ^ k 1 + B u k 1 ,
χ ˜ k , i w = h α ^ k 1 Δ α ^ k 1 χ ˜ k , i w j = 1 k ( 1 ) j α ^ k 1 j x ^ k j .
Resume, the TEA consists of three sub-filters whose required separate sets of parameters and initial conditions. Parameters of the order estimation filter KFo are denoted with upper index o (e.g., P ˜ k o , Q k 1 o ), parameters of KFw are denoted with upper index w (e.g., P ˜ k w , Q k 1 w ) and parameters of KFx are without upper index. A detailed description of the Triple Estimation Algorithm is presented in [40].

4. Identification and Analysis of Fractional Order Noise Parameters

Before we apply the Triple Estimation algorithm to real plant data (noises estimation of MEMS sensor), we will present the results of some numerical experiments for constant and FVO systems. The one state variable discrete state-space system, used in numerical experiments, is given as follows:
Δ k + 1 α k + 1 0 A x k + 1 = f x k + u k + ω k ,
x k + 1 = h α k + 1 Δ k + 1 α k + 1 0 x k + 1 j = 1 k + 1 ( 1 ) j α k + 1 j x k j + 1 ,
y k = x k + ν k .

4.1. Analysis of Fractional Constant and Variable Order System with Input Signal Known

This section contains a validation of the TEA for analysis and identification of fractional constant and variable order systems. Sets of numerical examples show the capability of the state, parameter and order estimation for known input signal to be a Gaussian noise. The problem formulated in Examples 1–3 is as follows: Estimate the state, order, and parameter of the fractional order system described by (43)–(45) with known input signal being the Gaussian noise. All numerical examples were conducted in Matlab/Simulink environment based on Fractional Variable Order Toolbox [43] with sampling time equals to h = 0.001 s. It is worth noticing that the final results of state, order, and parameter estimation depend on Triple Estimation filters’ parameters (KFx, KFo, and KFw) and should be individually adjusted.
The parameters of TEA applied in Example 1 are as following:
  • Noises parameters
    E [ ω ω T ] = 10 4 ,
    E [ ν ν T ] = 10 3 ,
  • Parameters of KFx filter
    P 0 = 0.01 , Q 0 = 10 4 ,
    x 0 = 0 , R = 10 3 ,
  • Parameters of KFo filter
    P 0 o = 0.01 , Q 0 o = 0.1 ,
    α 0 = 1 , R o = 10 3 , A = 1 , B = 2 , δ o = 0.5 ,
  • Parameters of KFw filter
    P 0 w = 0.01 , Q 0 w = 0.1 ,
    w 0 = 0 , R w = 10 3 , A = 1 , B = 2 , δ w = 0.5 .
An identification of fractional constant order system based on TEA is presented in Example 1.
Example 1.
Let us consider the DFVOSS A -type system given by (43)–(45), where
A = f = 1 , B = 1 , C = 1 , α k = 0.35 , u k N ( 0 , 1 ) .
Numerical results of state, order, and parameter estimation are presented in Figure 2, Figure 3 and Figure 4, respectively. As it can be seen in Figure 2, the state estimation overlaps the original one with high accuracy. Similarly, we can distinguish a high accuracy of order and parameter estimation presented in Figure 3 and Figure 4. The order estimation reaches the original one practically in unnoticeable time. The discrepancy between system parameter and its estimation decays. Due to the initial order value, the estimated order relatively slowly points out the simulated one. The desired simulation value is reached out with decreasing fluctuations at 5 s. In conclusion, the state, order, and parameter estimation using the TEA precisely reflect simulation values. The filter’s parameters can adjust the quality of achieved results individually. It mainly depends on ongoing cases cause the TEA can be treated as a convenient tool for identifying and analysing FO systems and noises. Then the rate of precisely achieved results does not play a significant role. Once identified, FO models can be applied to filter out useless data.
Temperature sensitivity is one of the most significant effects, exerting a huge impact on inertial measurement unit’s noises built-in MEMS technology. Temperature-varying can force the modification in order values during the process, and it is worth to validate the TEA for this case. Therefore, the next examples are devoted to state, order, and parameter estimation of FVO systems. The main difference between these two examples occurs in the order function. In Example 2, the original order function is described by linear function, and in Example 3 reference order function has a parabolic character.
The parameters of TEA applied in Examples 2 and 3 are as follows:
  • Noises parameters
    E [ ω ω T ] = 10 4 ,
    E [ ν ν T ] = 10 5 ,
  • Parameters of KFx filter
    P 0 = 0.01 , Q 0 = 10 4 ,
    x 0 = 0 , R = 10 5 ,
  • Parameters of KFo filter
    P 0 o = 0.1 , Q 0 o = 0.001 ,
    α 0 = 0.2 , R o = 10 5 , A = 1 , B = 2 , δ o = 0.5 ,
  • Parameters of KFw filter
    P 0 w = 0.01 , Q 0 w = 0.001 ,
    w 0 = 0.5 , R w = 10 5 , A = 1 , B = 2 , δ w = 0.5 .
Example 2.
Let us consider the DFVOSS A -type system given by (43)–(45), where
A = f = 1 , B = 1 , C = 1 , u k N ( 0 , 0.1 ) , α k = 0.2 + 0.1 k h f o r   k = 1 , 2 , 3 , . .
The example results of state, order, and parameter estimation are shown in Figure 5, Figure 6 and Figure 7, respectively. In this case, the estimated order precisely overlaps the original one until 4 s, and after that, a minor discrepancy appears. Additionally, there is a small underestimation between the original and algorithm’s parameter. However, both estimated values (order and parameter) lead to high accuracy of state estimation (see Figure 5). The example shows the behaviour of TEA in a wide range of order values while a linear function represents the order function. The estimated order coincides with the simulation one very well. It is a significant advantage of the TEA because the temperature effects influence noise components such as random walk errors or bias instability and directly into the model’s order. It shows that despite unexpected disturbances appearing in the FVO systems, the TEA can be applied for their identification and analysis.
Example 3.
Let us consider the DFVOSS A -type system given by (43)–(45), where
A = f = 1 , B = 1 , C = 1 , u k N ( 0 , 0.1 ) , α k = 0.2 + 0.05 · ( k h ) 2 f o r   k = 1 , 2 , 3 , . .
This time, the results of state, order, and parameter estimation are presented in Figure 8, Figure 9 and Figure 10, respectively. Moreover, the order-varying is forced by parabolic function, and despite of it, high accuracy of state, order, and parameter estimation is achieved. The accelerometers are measurement devices that work well while not being subjected to external forces or disturbances. They are susceptible, and additional actions can imply undesirable errors. It is an essential issue while double integral is needed to specify the position in inertial navigation systems (INS) and any error accumulated inaccuracies very fast. Hence, in this example, we show the situation when the order function changes its value more quickly than in Example 2. It conveys the possibilities of TEA against additional errors.
The results directly coming from Examples 1–3 show that TEA can be successfully used to estimate state, order, and parameter for fractional constant and variable order systems with known input signal. Especially, the capability of FVO systems analysis seems to be very important during noise modelling of MEMS. It is caused by the fact that the temperature effect can change the order of noise. Hence, the time-varying order occurring in the noise model makes it much more flexible and reflects its random walk errors.
The average execution time of TEA equals 495 s for a set of 6000 combined samples of state, order, and parameter estimation, which corresponds to Examples 1 and 2. The time-consuming tests were conducted on PC with Intel Core i7-5500U CPU, 2.4 GHz, RAM 8 GB and Matlab version 2021b 64 bit.

4.2. Identification without Input Signal Knowledge

When an input signal is not measured, the identification process can differ from desired values or order and system parameter. This can be explained by the fact that in practice, the system noise can have some unknown dynamical correlation of some order and parameter. Let us assume the fractional noise system equation in the following form:
Δ 1 α x k + 1 = f 1 x k + ω k ,
where ω k is a system noise also containing the fractional order dynamical correlation described by the following relation:
Δ 2 α ω k + 1 = f 2 ω k + ω k ,
where ω k is assumed to be noise without dynamical correlation.
By combining both equations we obtain
Δ 1 α x k + 1 1 f 2 Δ 2 α ω k + 1 = f 1 x k 1 f 2 ω k .
As we can see, this dynamical correlation can have a direct effect on estimated order and system parameter in the estimation process, which can make obtained estimation results different from those assumed in numerical models, because they take into consideration also the dynamical correlation of the source noise. However, it will not be a problem in estimation of real plant noise because the aim of estimation is to find the most appropriate model with the assumption that the source noise is without dynamical correlation.

5. Identification and Analysis of MEMS Accelerometer’s Noises

This section contains the experimental results of noises modelling for the 3-axes accelerometer being the part of SparkFun MPU9250 Inertial Measurement Unit (IMU) built-in MEMS technology. This is a very popular MEMS, and there exist plenty of robotics projects where it is a base for many Attitude and Heading Reference System (AHRS) and INS implementations. Generally, accelerometers are very sensitive, and they are characterised by high frequencies noises. All of these can have a significant impact on final project results. This is why we decided to apply the TEA for the noises modelling of the accelerometer as part of MPU9250 breakout board.
Overall, the MPU9250 unit is a 9 degree of freedom MEMS with 3 accelerometer’s axes, 3 gyroscope’s axes, and 3 magnetometer’s axes. The MPU9250 breakout board runs on 3.3 VDC and contains I2C and SPI communication protocols.

5.1. Experimental Setup

The experimental approach was divided into two main parts. The first one was devoted to data collection, and the second one was related to their identification based on TEA. Therefore, the first phase can be represented by experimental setup shown in Figure 11. Its main parts are the Arduino Due development board and MPU9250 breakout board. The I2C communication protocol was used for data transmission between them. Additionally, the accelerometer’s range was set to 8 g, and measurement’s data were collected in the Arduino IDE environment for stationary located IMU, with sampling time h = 0.01 s. The lock accelerometer position allows gathering its noises corrupted by constant gravity components. In the second phase of the experiment, the calculated mean value of data for each axis was subtracted from its measurement to have the pure noise corresponding to a particular accelerometer axis. Then, post-processing data for each axis without mean value were adapted as x-axis, y-axis, and z-axis noises in the Matlab environment. Having such prepared noises data and due to the fact that their mathematical models are independent from each other, we decided to apply the TEA separately for each axis. In fact, the state, parameter, and order estimation results were achieved for each axis separately under the same TEA configuration parameters except covariances R, R o , and R w adjusted according to noise measurement corresponding to each axis.
To summarize, in the experimental setup, the problem is formulated as follows: Estimate noise, its order and parameter for x, y and z-axes separately, with no input signal knowledge. Therefore, the real noises data originating from the accelerometer do not take part in the TEA, but they are shown only for validation purposes.
The parameters of TEAs for each acceleromenter’s axis are as the following:
  • Parameters of KFx filter
    P 0 = 0.01 , Q 0 = 0.01 ,
    x 0 = 0 , R = 4 · 10 4   for   x - axis noise 3.9 · 10 4   for   y - axis noise 11 · 10 4   for   z - axis noise
  • Parameters of KFo filter
    P 0 o = 0.01 , Q 0 o = 0.1 ,
    α 0 = 1 , A = 1 , B = 2 , δ o = 0.5 , R o = 4 · 10 4   for   x - axis noise 3.9 · 10 4   for   y - axis noise 11 · 10 4   for   z - axis noise
  • Parameters of KFw filter
    P 0 w = 0.01 , Q 0 w = 0.1 ,
    w 0 = 0 , A = 1 , B = 2 , δ w = 0.5 , R w = 4 · 10 4   for   x - axis noise 3.9 · 10 4   for   y - axis noise 11 · 10 4   for   z - axis noise .

5.2. Experimental Results

In the results of the experimental research, we have obtained the three sets of plots corresponding to x, y, and z-axes of accelerometer’s noises. Therefore, the first set of plots for noise, order, and parameter estimation related to x-axis are presented in Figure 12, Figure 13 and Figure 14, respectively. The estimated order rapidly goes from the initial value to around 0.3 and then, with minor fluctuations, raises to the neighbourhood of 0.4 at time 15 s, while parameter estimation slowly decreases to −4 and reaches this value after 15 s on the plot. We can observe that after 15 s, the order and parameter are relatively stabilised. The plots concerning to y-axis noise together with order and parameter estimation can be found in Figure 15, Figure 16 and Figure 17. This time, the estimated order achieve its stabilised value around 0.3 , and the curve reaches it approximately in 2 s. Then, the fluctuations are relatively slight during the whole estimation process. The estimated parameter goes to −1 rapidly, in approximately 2 s and then with small oscillations tends to a neighbourhood of −1.4. At the end, the experimental results of noise modelling for z-axis are presented in Figure 18. Additionally, its order and parameter estimation are shown in Figure 19 and Figure 20, respectively. For z-axis noise, the estimated order reaches its stabilised value around 0.4 at time 3 s and then tries to keep it. The order curve is relatively slight with minor fluctuations, while the estimated parameter starting with the initial value tends to −2.5, and the plot follows it beginning with 4 s.
Comparing the original noises for x, y, and z-axes to their estimation results presented in Figure 12, Figure 15, and Figure 18 confirm the fractional dynamics of MPU9250 accelerometer’s noises and the high quality of applied method. Moreover, we can see that order values for each axis noise are similar, and they tend to fractional constant values. Hence, we can conclude that the temperature effects did not impact the accelerometer’s random walk errors too much.
The average execution time of TEA equals 122 s for a set of 3000 combined samples of state, order, and parameter estimation. This time corresponds to each accelerometer’s axis, separately.

6. Conclusions

The paper presents the results of the Triple Estimation Algorithm application to estimate fractional order noises. The algorithm can simultaneously estimate state variables, system parameters, and FVO of the noises. Presented numerical experiments show the ability of the algorithm to determine the variable order for the FVO noise. Simultaneous estimation can be an essential feature to test if unknown noise obtained experimentally is stationary or time-varying. As it is known, the low-cost MEMS used in many projects done by hobbyists are temperature sensitive. In this case, the capability of FVO noise analysis can handle its unexpected influence on final results and raise their accuracy. Comparing the real accelerometer noises to their estimates, it could be seen that the TEA was successfully applied to noises modelling of 3-axes accelerometer being the part of MPU9250 breakout board. The experimentally obtained data confirm the fractional character of investigated sensor’s noises, and the estimated order values for x, y and z-axes are similar to each other. Carrying out tests in examples and experiments acknowledges that the TEA can be treated as an effective method for simultaneously estimating states, variable order, and system parameters. Moreover, the main advantage of the proposed algorithm against other existing techniques is that it obtains the unknown system parameters and variable order at once, giving the ability to analyse the time-correlated noise. The obtained results can be extended to examine external factors (like temperature) on the noise order and system parameters. Another possible area of use for the obtained results is the development of estimation algorithms for AHRS or inertial navigation, which will use estimated values of fractional noises. Numerical simulations show that accurate results of state, order, and parameter estimation given by TEA come with the complexity of the algorithm which is relatively high and time consuming.

Author Contributions

Conceptualisation, D.S. and M.M.; data curation, M.M.; formal analysis, D.S. and M.M.; methodology, D.S. and M.M.; project administration, M.M.; software, D.S. and M.M.; validation, M.M.; visualisation, D.S. and M.M.; writing—original draft, D.S., M.M. and W.M.; writing—review & editing, D.S., M.M. and W.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Scientific Council for Automatic Control, Electronics, and Electrical Engineering, Warsaw University of Technology as an internal grant admitted in 2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vasudev, A.; Bhansali, S. Microelectromechanical systems (MEMS) for in vivo applications. In Implantable Sensors Systems for Medical Applications; Woodhead Publishing Series in Biomaterials; Inmann, A., Hodgins, D., Eds.; Woodhead Publishing Ltd.: Swatson, UK, 2013; Volume 52, pp. 331–358. [Google Scholar] [CrossRef]
  2. Panescu, D. MEMS in medicine and biology. IEEE Eng. Med. Biol. Mag. 2006, 25, 19–28. [Google Scholar] [CrossRef]
  3. Deisingh, A. MEMS technology in analytical chemistry. Analyst 2003, 128, 9–11. [Google Scholar] [CrossRef]
  4. Leclerc, J. MEMs for aerospace navigation. IEEE Aerosp. Electron. Syst. Mag. 2007, 22, 31–36. [Google Scholar] [CrossRef]
  5. Classen, J.; Frey, J.; Kuhlmann, B.; Ernst, P. MEMS gyroscopes for automotive applications. In Advanced Microsystems for Automotive Applications 2007; VDI-Buch Series; Valldorf, J., Gessner, W., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 291–306. [Google Scholar]
  6. Goel, A.; Ul Islam, A.; Ansari, A.; Kouba, O.; Bernstein, D.S. An Introduction to Inertial Navigation From the Perspective of State Estimation [Focus on Education]. IEEE Control. Syst. Mag. 2021, 41, 104–128. [Google Scholar] [CrossRef]
  7. Yazdi, N.; Ayazi, F.; Najafi, K. Micromachined inertial sensors. Proc. IEEE 1998, 86, 1640–1659. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, F.; Fei, J.; Xue, Y. Double Recurrent Perturbation Fuzzy Neural Network Fractional-Order Sliding Mode Control of Micro Gyroscope. IEEE Access 2021, 9, 55352–55363. [Google Scholar] [CrossRef]
  9. Mohd-Yasin, F.; Nagel, D.J.; Korman, C.E. Noise in MEMS. Meas. Sci. Technol. 2010, 21. [Google Scholar] [CrossRef]
  10. Krysko, V., Jr.; Awrejcewicz, J.; Yakovleva, T.; Kirichenko, A.; Szymanowska, O.; Krysko, V. Mathematical modeling of MEMS elements subjected to external forces, temperature and noise, taking account of coupling of temperature and deformation fields as well as a nonhomogenous material structure. Commun. Nonlinear Sci. Numer. Simul. 2019, 72, 39–58. [Google Scholar] [CrossRef]
  11. Qiao, Y.; Arabi, M.; Xu, W.; Zhang, H.; Abdel-Rahman, E.M. The impact of thermal-noise on bifurcation MEMS sensors. Mech. Syst. Signal Process. 2021, 161, 107941. [Google Scholar] [CrossRef]
  12. Romanovas, M.; Klingbeil, L.; Traechtler, M.; Manoli, Y. Application of fractional sensor fusion algorithms for inertial MEMS sensing. Math. Model. Anal. 2009, 14, 199–209. [Google Scholar] [CrossRef]
  13. Samko, S.; Kilbas, A.; Maritchev, O. Fractional Integrals and Derivative. Theory and Applications; Gordon & Breach Sci. Publishers: London, UK, 1987. [Google Scholar]
  14. Miller, K.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differenctial Equations; John Wiley & Sons Inc.: New York, NY, USA, 1993. [Google Scholar]
  15. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu, V. Fractional-Order Systems and Controls; Springer: London, UK, 2010. [Google Scholar]
  16. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  17. Magin, R.; Ortigueira, M.D.; Podlubny, I.; Trujillo, J. On the fractional signals and systems. Signal Process. 2011, 91, 350–371. [Google Scholar] [CrossRef]
  18. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, Volume 204 (North-Holland Mathematics Studies); Elsevier Science Inc.: San Diego, CA, USA, 2006. [Google Scholar]
  19. Das, S. Introduction to Fractional Calculus. In Functional Fractional Calculus; Springer: Berlin/Heidelberg, Germany, 2011; pp. 1–50. [Google Scholar] [CrossRef]
  20. Anastassiou, G.A. Generalized Fractional Calculus; Springer: Cham, Switzerland, 2021. [Google Scholar]
  21. Yang, X.J. General Fractional Derivatives: Theory, Methods and Applications; Chapman and Hall/CRC: Boca Raton, FL, USA, 2019. [Google Scholar]
  22. Tarasov, V.E. Generalized Memory: Fractional Calculus Approach. Fractal Fract. 2018, 2, 23. [Google Scholar] [CrossRef] [Green Version]
  23. Sierociuk, D.; Dzielinski, A.; Sarwas, G.; Petras, I.; Podlubny, I.; Skovranek, T. Modelling heat transfer in heterogeneous media using fractional calculus. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2013, 371, 20120146. [Google Scholar] [CrossRef] [Green Version]
  24. Sakrajda, P.; Sierociuk, D. Modeling Heat Transfer Process in Grid-Holes Structure Changed in Time Using Fractional Variable Order Calculus. In Theory and Applications of Non-Integer Order Systems; Babiarz, A., Czornik, A., Klamka, J., Niezabitowski, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 297–306. [Google Scholar]
  25. Reyes-Melo, M.; Martinez-Vega, J.; Guerrero-Salazar, C.; Ortiz-Mendez, U. Application of fractional calculus to modelling of relaxation phenomena of organic dielectric materials. In Proceedings of the 2004 IEEE International Conference on Solid Dielectrics, Toulouse, France, 5–9 July 2004. [Google Scholar]
  26. Dzielinski, A.; Sierociuk, D.; Sarwas, G. Some applications of fractional order calculus. Bull. Pol. Acad. Sci. Tech. Sci. 2010, 58, 583–592. [Google Scholar] [CrossRef] [Green Version]
  27. Ortigueira, M.D.; Valério, D. Fractional Signals and Systems; De Gruyter: Berlin, Germany, 2020. [Google Scholar] [CrossRef]
  28. Sheng, H.; Chen, Y.; Qiu, T. Fractional Processes and Fractional-Order Signal Processing; Springer: London, UK, 2012. [Google Scholar]
  29. Muresan, C.I.; Birs, I.R.; Dulf, E.H.; Copot, D.; Miclea, L. A Review of Recent Advances in Fractional-Order Sensing and Filtering Techniques. Sensors 2021, 21, 5920. [Google Scholar] [CrossRef]
  30. Sierociuk, D.; Malesza, W.; Macias, M. Derivation, interpretation, and analog modelling of fractional variable order derivative definition. Appl. Math. Model. 2015, 39, 3876–3888. [Google Scholar] [CrossRef]
  31. Sierociuk, D.; Malesza, W.; Macias, M. On the Recursive Fractional Variable-Order Derivative: Equivalent Switching Strategy, Duality, and Analog Modeling. Circuits Syst. Signal Process. 2015, 34, 1077–1113. [Google Scholar] [CrossRef] [Green Version]
  32. Macias, M.; Sierociuk, D. An alternative recursive fractional variable-order derivative definition and its analog validation. In Proceedings of the ICFDA’14 International Conference on Fractional Differentiation and Its Applications, Catania, Italy, 23–25 June 2014. [Google Scholar]
  33. Sierociuk, D.; Ziubinski, P. Variable order fractional Kalman filters for estimation over lossy network. Lect. Notes Electr. Eng. 2015, 320, 285–294. [Google Scholar] [CrossRef]
  34. Ziubinski, P.; Sierociuk, D. Improved fractional Kalman filter for variable order systems. In Proceedings of the ICFDA’14 International Conference on Fractional Differentiation and Its Applications, Catania, Italy, 23–25 June 2014. [Google Scholar] [CrossRef]
  35. Wyss, W. Fractional noise. Found. Phys. Lett. 1991, 4, 235–246. [Google Scholar] [CrossRef]
  36. Chen, Y.; Sun, R.; Zhou, A.; Zaveri, N. Fractional order signal processing of electrochemical noises. J. Vib. Control. 2008, 14, 1443–1456. [Google Scholar] [CrossRef]
  37. Sierociuk, D.; Ziubinski, P. Fractional order estimation schemes for fractional and integer order systems with constant and variable fractional order colored noise. Circuits Syst. Signal Process. 2014, 33, 3861–3882. [Google Scholar] [CrossRef] [Green Version]
  38. Bai, Y.; Wang, X.; Jin, X.; Su, T.; Kong, J.; Zhang, B. Adaptive filtering for MEMS gyroscope with dynamic noise model. ISA Trans. 2020, 101, 430–441. [Google Scholar] [CrossRef] [PubMed]
  39. Ziubinski, P.; Sierociuk, D. Fractional order noise identification with application to temperature sensor data. In Proceedings of the 2015 IEEE International Symposium on Circuits and Systems (ISCAS), Lisbon, Portugal, 24–27 May 2015; pp. 2333–2336. [Google Scholar] [CrossRef]
  40. Sierociuk, D.; Macias, M. Triple Estimation of Fractional Variable Order, Parameters, and State Variables Based on the Unscented Fractional Order Kalman Filter. Sensors 2021, 21, 8159. [Google Scholar] [CrossRef] [PubMed]
  41. Sierociuk, D.; Malesza, W. Fractional variable order discrete-time systems, their solutions and properties. Int. J. Syst. Sci. 2017, 48, 3098–3105. [Google Scholar] [CrossRef]
  42. Haykin, S. Kalman Filtering and Neural Networks; John Wiley & Sons, Inc.: New York, NY, USA, 2001. [Google Scholar]
  43. Sierociuk, D. Fractional Variable Order Derivative Simulink Toolkit. MATLAB Central File Exchange. 2012. Available online: http://www.mathworks.com/matlabcentral/fileexchange/38801-fractional-variable-order-derivative-simulink-toolkit (accessed on 10 December 2021).
Figure 1. The Triple Estimation Algorithm scheme.
Figure 1. The Triple Estimation Algorithm scheme.
Sensors 22 00527 g001
Figure 2. Original and estimated state variable from Example 1.
Figure 2. Original and estimated state variable from Example 1.
Sensors 22 00527 g002
Figure 3. Original and estimated order from Example 1.
Figure 3. Original and estimated order from Example 1.
Sensors 22 00527 g003
Figure 4. Original and estimated parameter from Example 1.
Figure 4. Original and estimated parameter from Example 1.
Sensors 22 00527 g004
Figure 5. Original and estimated state variable from Example 2.
Figure 5. Original and estimated state variable from Example 2.
Sensors 22 00527 g005
Figure 6. Original and estimated order from Example 2.
Figure 6. Original and estimated order from Example 2.
Sensors 22 00527 g006
Figure 7. Original and estimated parameter from Example 2.
Figure 7. Original and estimated parameter from Example 2.
Sensors 22 00527 g007
Figure 8. Original and estimated state variable from Example 3.
Figure 8. Original and estimated state variable from Example 3.
Sensors 22 00527 g008
Figure 9. Original and estimated order from Example 3.
Figure 9. Original and estimated order from Example 3.
Sensors 22 00527 g009
Figure 10. Original and estimated parameter from Example 3.
Figure 10. Original and estimated parameter from Example 3.
Sensors 22 00527 g010
Figure 11. The real view of experimental setup with Arduino Due development board and MPU9252 IMU mounted on a shaft of servo motor in lock position.
Figure 11. The real view of experimental setup with Arduino Due development board and MPU9252 IMU mounted on a shaft of servo motor in lock position.
Sensors 22 00527 g011
Figure 12. Original and estimation of x-axis noise.
Figure 12. Original and estimation of x-axis noise.
Sensors 22 00527 g012
Figure 13. Order estimation for x-axis noise.
Figure 13. Order estimation for x-axis noise.
Sensors 22 00527 g013
Figure 14. Parameter estimation for x-axis noise.
Figure 14. Parameter estimation for x-axis noise.
Sensors 22 00527 g014
Figure 15. Original and estimation of y-axis noise.
Figure 15. Original and estimation of y-axis noise.
Sensors 22 00527 g015
Figure 16. Order estimation for y-axis noise.
Figure 16. Order estimation for y-axis noise.
Sensors 22 00527 g016
Figure 17. Parameter estimation for y-axis noise.
Figure 17. Parameter estimation for y-axis noise.
Sensors 22 00527 g017
Figure 18. Original and estimation of z-axis noise.
Figure 18. Original and estimation of z-axis noise.
Sensors 22 00527 g018
Figure 19. Order estimation for z-axis noise.
Figure 19. Order estimation for z-axis noise.
Sensors 22 00527 g019
Figure 20. Parameter estimation for z-axis noise.
Figure 20. Parameter estimation for z-axis noise.
Sensors 22 00527 g020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Macias, M.; Sierociuk, D.; Malesza, W. MEMS Accelerometer Noises Analysis Based on Triple Estimation Fractional Order Algorithm. Sensors 2022, 22, 527. https://doi.org/10.3390/s22020527

AMA Style

Macias M, Sierociuk D, Malesza W. MEMS Accelerometer Noises Analysis Based on Triple Estimation Fractional Order Algorithm. Sensors. 2022; 22(2):527. https://doi.org/10.3390/s22020527

Chicago/Turabian Style

Macias, Michal, Dominik Sierociuk, and Wiktor Malesza. 2022. "MEMS Accelerometer Noises Analysis Based on Triple Estimation Fractional Order Algorithm" Sensors 22, no. 2: 527. https://doi.org/10.3390/s22020527

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