1. Introduction
Recent advances in Micro Electro Mechanical Systems (MEMS) based inertial sensors are quite significant in that they promise to be smaller and cheaper systems. MEMS is the integration of mechanical elements, sensors, actuators, and electronics on a common silicon substrate through the utilization of microfabrication technology (
Figure 1). MEMS are expected to revolutionize many product categories by bringing together silicon based microelectronics with micromachining technology, and enabling complete systems-on-a-chip to be realized [
1].
Among the applications of MEMS technology, MEMS- based inertial sensors such as MEMS based gyroscopes and MEMS based accelerometers have been adopted as aiding sensors to improve navigation information continuity. While the MEMS based gyroscope is a relatively new technology and still an ongoing research activity for commercial uses, the MEMS based accelerometer has been widely utilized in a variety of developments and its success is significant.
In the navigation field, it has been reported that many efforts have been made to make the smaller and less expensive navigation devices available for more users. MEMS based inertial sensors have recently drawn great attention as aiding GPS outages with low inherent cost, small size, low power consumption, and solid reliability. However, their performance is still considered poor in accuracy for certain applications. In this paper, the MEMS based inertial sensors' performance will be investigated and a land vehicle navigation system prototype will be developed by integrating MEMS based inertial sensors, a low cost GPS receiver and a digital compass.
In spite of the smaller size and cost effectiveness of MEMS based inertial sensors, the error behaviour of MEMS based inertial sensors must be appropriately treated in order to turn the raw sensor measurements into reliable and useful data for vehicle position determination. When we confine the scope of application of MEMS based inertial sensors to aiding GPS solutions to a relatively short period of time, some deterministic error sources (zero-offset bias and 1st order scale factor) and stochastic variation (random noise) can be considered as the main concerns to be discussed among the different types of error sources for MEMS based inertial sensors. Besides, the understanding of their stochastic variations is of significant importance for the development of optimal estimation algorithms.
In the subsequent sections, the conventional error model of inertial sensors will be simplified considering MEMS based sensor design and a short time period usage assumption. The deterministic error sources will be estimated by using multi-position test which is well described in the reference [
2], and the stochastic variation will be modeled by 1
st order Gauss-Markov (GM), which has been widely used in navigation field, and a higher order AutoRegressive (AR) model introduced in [
3]. The deterministic error sources (zero-offset bias and 1
st order scale factor) of MEMS based inertial sensors estimated by using multi-position testing in the laboratory will be referenced to initial measurement in kinematic environments. For the stochastic variation, not only the conventional 1
st order GM model but also a higher order AR model will be used in optimal estimation algorithm (i.e. Kalman filter) to quantify the effect of precise stochastic modeling method for MEMS based inertial sensor applications in kinematic environments. When the performance of MEMS-based inertial sensors are admissible for a certain application such as land vehicle navigation, a continuous integrated navigation system will be available by integrating GPS with cheaper and smaller inertial sensors in urban canyons with GPS signal blockages.
2. Accelerometer/Gyroscope Error Model and Stochastic Modeling
There are two major aspects that should be considered in the error analysis of any MEMS-based sensor: (1) error analysis to identify deterministic error and non-deterministic (stochastic) error sources; and (2) the development of stochastic modeling methods used to characterize the random part of the sensor output.
2.1. Error Sources and Error Models
Current commercial accelerometers/gyroscopes are mainly classified as either mechanical or solid-state. As mentioned before, all accelerometers/gyroscopes are suffering from a variety of error sources which are slightly different depending upon different types of the manufacturing principles. The error equation of conventional mechanical inertial sensors from the reference [
2] will be first introduced and the error equation will then be simplified according to the tolerance of a specific application such as land vehicle navigation system and MEMS technology.
Conventionally, the measurement in the X-axis provided by accelerometer (
ãx ) can be expressed in terms of the applied acceleration acting along its sensitive axis (
ax ) and the accelerations acting along the pendulum and hinge axes,
ay and
az respectively, by the equation [
2]:
where
Sx is the scale factor error, usually expressed in polynomial form to include non-linear effects,
My ,
Mz are the cross-axis coupling factors,
Bf is the measurement zero-offset bias,
Bv is the vibro-pendulous error coefficient, and
nx is the random noise.
For an accelerometer based on MEMS technology and non-pendulous design, it is reasonable to expect that the cross-axis coupling factors and vibro-pendulous error would be insignificant because most MEMS accelerometers are assembled as three single-axis accelerometers so that they have low cross-axis coupling factors [
4]. Then, the conventional error model can be simplified as below,
As indicated by
Equation (2), the zero-offset bias and the 1
st order scale factor are the main concerns for the deterministic error sources and the last term is the stochastic variation of the sensor output. The Y-axis and Z-axis measurements can be expressed in the same way.
Similarly, current commercial gyroscopes utilize different development principles, resulting in various types of gyroscopes with distinct characteristics for each one. Accordingly, assuming the acceleration sensitive errors are negligible, the measured angular rate can be modeled for many applications as [
2]
where
Sx is the scale factor which may be expressed as a polynomial in
ωz to represent scale factor non-linearity,
My ,
Mz are the cross-axis coupling factors,
Bf is the measurement zero-offset bias,
nz is the random noise. Using the same assumption in the accelerometer case,
Equation (3) can be simplified as
in which only the zero-offset bias and 1
st order scale factor are included with significant contribution to the deterministic error sources.
Equations (2) and
(4) will be used to estimate the deterministic error sources (zero-offset bias and 1
st order scale factor) by using multi-position testing.
2.2. Stochastic Modeling
Considering only linear stationary stochastic processes, one way to specify a random process is to describe in detail the conceptual chance experiment giving rise to the process [
5]. As it can be seen, many signals are quite different, even with same mean and variance values, so it is clear that more information than just mean and variance is needed in order to describe the random process more precisely. It has been seen that the autocorrelation function, denoted as
RX (
τ) in the sequel, is an important descriptor of a random process that is relatively easy to obtain because it depends on only the second-order probability density for the process [
5]. Thus, if
RX (
τ) is known a-prior or if we can estimate
RX (
τ) from observational data, then we can use this information to help “identify” which of the special models (if any) would fit the process under study [
6].
If a stationary Gaussian process
Xt has an exponential Autocorrelation, it is called a Gauss-Markov (GM) process. As shown in
Figure 2, the autocorrelation and spectral functions of 1
st order GM process have the form of
The mean-square value and time constant for the process are given by the σ2 and 1/β parameters, respectively. The Gauss-Markov process is a very useful process in applied work because (i) it can fit a large number of physical processes with reasonable accuracy, and (ii) it has a relatively simple mathematical description. In positioning and navigation fields, 1st order Gauss-Markov process has been frequently used to describe the stochastic behaviors due to its simple estimation.
However, quite often even with simple representation of time-correlated signal behaviour, the estimated autocorrelation and its FFT transform are quite different from 1
st order Gauss-Markov process shown in
Figure 2. Considering the very noisy measurements and poor performance of MEMS based inertial sensors, the more precise and appropriate stochastic modeling is desirable. In this paper, a higher order AutoRegressive model will be applied and compared with 1
st order Gauss-Markov process.
Some of special discrete parameter stochastic models which can provide us with a structure for fitting models to practical data have been studied in probability and mathematical statistics theory such as the purely random noise {
εt}, the AutoRegressive model (AR), the MovingAverage model (MA), the mixed AutoRegressive/MovingAverage model (ARMA), the harmonic model, just to mention a few. Among the discrete parameter stochastic models, AR model has been widely adopted to describe random noise output of physical systems in many fields because the relatively simple parameter estimation and the value of {
εt } is “drawn into” the process {
Xt } and thus influences all future values,
Xt,
Xt+1,
Xt+2, ...... resulting in its autocorrelation function “dies out gradually” [
6]. The attempt to apply AR model for inertial sensors was first introduced in [
3] and its mathematical description can be found in [
3][
6][
7][
8].
With the standard form of AR model,
where
b0 is the variance
of
εt which represents a random process. For a variety of methods to estimate the parameters
ak the reader is referred to the literature ([
6], [
7]) and this paper will apply one of them. Also, if the process is both causal and stable, then all the poles of
H(
z) must lie inside the unit circle of the z-plane because the Region of Convergence (ROC) is of the form | z | > r
max, and since the unit circle is included in the ROC, one must have r
max < 1, where r
max equals the largest magnitude of any of the poles of
H(
z) [
9].
The most considered AR model based parameter estimation methods are the Yule-Walker, Burg and Unconstrained Least-Squares methods [
7]. For a large dataset, the results of these three different methods provide a fairly close estimation results of the parameters. However, there are still some different characteristics for each method to be noticed. The authors chose the Burg method, which is quite popularly used and estimates the reflection coefficients by minimizing both forward and backward prediction errors in the least square sense with the constraint that the AR parameters satisfy the Levinson-Durbin recursion. Unless
a priori information about the order of an AR Model was given, the order of the AR model is an unknown, so it needs to be estimated. Several method to determine the order of the AR Model have been reported in literatures such as Final Prediction Error (FPE), Akaike's Information Criterion (AIC), Minimum Description Length (MDL), Investigation of Residual Variance (IRV) and etc. One of the common ways to determine the order of the AR model is to investigate the residual variance for different orders. Assuming the true model is of finite order, as the estimated order is getting close to the true model, the residual variance wouldn't reduce significantly [
7]. It should be kept in mind that a higher order AR model would increase Kalman Filter error states. As a result, it would increase the computational loads and might result in unstable solutions [
3].
4. Multi-Position Testing
RGA300CA system consists of a MEMS angular rate gyro and a triaxial silicon MEMS accelerometer. The triaxial accelerometer is a bulk-micromachined capacitive accelerometer whose input range is ±2g.
First, RGA300CA has been tested in the rotation panel connected to SmartMotor from Animatics Corporation. The testing rate table was carefully leveled relative to the local gravity vector. Once an accelerometer was attached to the testing rate table properly, the accelerometer output was collected with constant speed of rotation. The actual measurements of X/Y/Z axes of the accelerometer were compared with the reference acceleration determined by the testing rate table orientation. Using the error model described in the previous section and the “best fit” line regression method, the deterministic error sources have been estimated and used for calibration. Both of accelerometer and rotation panel were turned on and off every time with about 10 minutes apart and were warmed up for about 5 minutes before each datalogging [
11]. The room temperature (about 21°C) was also monitored. The local gravity value (9.8080.m/s
2) in the Multi-Sensor Lab at the University of Calgary has been used for the reference gravity value. One set of the accelerometer measurements coinciding z-axis with rotation axis of rate table is shown in
Figure 6.
For the sensor, 20 Hz logging rate and approximate bandwidth 10 Hz were used with data logging systems as GYRO-VIEW from Crossbow for RGA300CA. Also, the model SM2330SQ version 4.11 motor was used with SMI 1.310 windows S/W from Animatics. The RGA300CA and the rotation panel were connected to separate computers with RS-232 port cables and the output of the sensor measurements was saved in text file format.
Due to the instability of the rotating motor and the initial alignments, the angular velocity should be calculated for every run. The zero-offset bias (
Bf ) and 1
st order scale factor (
Sx ) can be obtained by using Least Squares method with rotational measurements and reference gravity value. The bias and the scale factor stability results are given in
Figures 7 and
8 with their mean & standard deviation shown in
Table 2.
As discussed in the previous section, the deterministic error sources of RGA300CA Yaw rate gyroscope are zero-offset bias and 1
st order scale factor. Analogous to the accelerometer case, the simplified form of error
Equation (4) without any modification will be used to analyze the actual gyroscope's Yaw rate measurements. This time, the rotational table has been precisely leveled out horizontally to provide the reference angular rate which is supposed to be correspondent to Yaw rate of gyroscope assuming that Earth rotation rate effect is nominal.
During a typical test schedule, the rotation rate of the rate table is stepped through a series of angular rates starting from zero deg/s recording data at each stage. The rotation speed is kept constant for a period at each step and the sensor outputs are allowed to stabilize, before recording the output signals. The applied angular rate is varied in incremental steps between the maximum and minimum desired rotation rates. At each step, the signals from gyroscope are recorded when the sensor is in equilibrium [
2].
In this experiment, the applied rotation rate has been increased from 0 deg/s to 80 deg/s and then, decreased until negative 80 deg/s. After that, it resumed to increase from -80 deg/s to 0 deg/s. For each rotation rate steps (10 deg/s), the dwell time consists of stabilization time (about 10 seconds) and sample time (about 10 seconds). 33 subsets of data have been recorded based on the same scheme and combined together to compose a series of measurements. One of the testing results has been illustrated in
Figure 9.
Those recorded data has been averaged out to provide a list of measurements resulting in measurement matrix in Least Squares estimation scheme. Accordingly, two parameters (zero-offset bias and 1
st order scale factor) could be estimated by simple Least Squares process with 33 measurements. The same test has been performed ten times with approximately 30 minutes interval. The results of the ten tests with their mean & standard deviation for Yaw-rate are shown in
Table 3. The mean values of
Table 2 and
Table 3 will be referred in the initial static leveling in kinematic testing.
5. Stochastic Modeling of RGA300CA
Based on the discussion in previous sections, the stochastic variation (random noise) of the experimental output of accelerometer and gyroscope inside RGA300CA will be analyzed and modeled appropriately. Since the usage of 1st Gauss-Markov model is well known in many literatures, only the stochastic modeling by AutoRegressive model will be described here. AR model parameters will be estimated by using Burg method and corresponding order of AR model will be approximated by investigating the residual variance in accordance to different orders considering the increase of state vector in Kalman Filter error state.
Shown in
Figures 10 and
11 are the normalized autocorrelation functions of a triaxial accelerometer and yaw rate gyroscope outputs in well-leveled static mode with about 3 hours warming time. It is clear that the temperature variation of the sensor unit affects the sensor measurements significantly. It is well indicated in many literatures that the temperature is the main concern of sensor output stability. That is why measurements after three hours were being used for the analysis. Therefore, the relatively stable parts of the original accelerometer/gyroscope measurements were only used and their trends were removed.
As expected, the RGA300CA measurements have shown very short correlation times in
Figures 10 and
11. For the gyroscope, its stochastic variation has behaved as quite purely random process (white noise). On the other hand,
Figure 10 has indicated that there exists time correlated behaviour in the stochastic variation of the accelerometer measurements and it should be modeled precisely.
There are two main steps involved in AR modeling, namely, parameter estimation and order determination. Once the three parameter estimation methods in the previous section were performed to estimate the parameters using the sample dataset (about 8 hours with a sampling rate of 20Hz), they have provided very close results from one to the other. Therefore, in spite of some distinct characteristics, any methods could be used in this experiment.
In the testing, the Burg method has been applied. To assess the proper determination of the order for AR model, the estimated residual variance
in accordance to different orders was chosen to be analyzed. In order to avoid abrupt increase in the error states of the Kalman filter due to the increase of the order of the AR model, an appropriate order ought to be determined when the variance plot starts to be leveled out in
Figure 12.
Based on the results shown in
Figures 12 and
13, a 4
th order AR model for the stochastic variation of the accelerometer has been chosen and its parameters were estimated by the Burg method. A total of 10 sample datasets have been used for the estimation of the AR model parameters and the mean values will be used for bias modeling in optimal estimation algorithm for kinematic testing. Therefore, the stochastic variation of accelerometer output will be modeled by 1
st order GM and 4
th AR models and the stochastic variation of gyroscope output will be modeled by purely random process (white noise) based on
Figure 11 and
13. Now the error state vector includes 14 states, namely 2 for position, 2 for velocity, 1 for misalignment, 8 for 4
th order AR model of X/Y-axis of accelerometer biases and 1 for white noise for yaw rate bias. The corresponding dynamic model driven by the white noise
w and the measurement model equations have the expressions as follows:
Again, the parameter estimation results in the multi-position testing have been used to construct the state transition matrix and the covariance matrix (Q ) associated with wk in the dynamic model while the covariance matrix (R) associated with vk in the measurement model has been constructed by using RMS values of GPS trajectory accuracy compared with PPP solutions.
6. Kinematic Testing and Results
To fulfill the testing assumptions described previously the kinematic testing has been conducted in one of the parking lots in the University of Calgary which is a relatively flat and open area. The vehicle was driven at speeds of 10 to 30 km/hr with six major turns. Around the corners, the speed was reduced and then, was accelerated along the straight path comparatively. The same driving testing was repeated 10 times with the same routine of data collection in the same area.
As illustrated and also described previously, the data collection consists of two parts, namely, static mode and kinematic mode after warming up and compass calibration periods. The static mode dataset for about five minutes was referenced with zero-offset bias estimation which was explained in multi-position testing. With an initial position from PPP processing of GPS measurements and initial heading corrected by CGRF, a Kalman filter error estimation has been conducted which is composed of a dynamic model using measurements of X/Y axes of accelerometer, Yaw rate gyroscope of RGA300CA and a measurement model using measurements of GPS-9543 module. The trajectories of the integrated system using 4th order AR model and 1st order GM model were first generated with the 1-sec update interval and were then compared with the PPP solution trajectory and GPS-9543 solution trajectory. After that, the system trajectories were generated with 5-sec, 10-sec, 20-sec, 30-sec, 60-sec update intervals.
In the 1-sec update case, the system trajectories by both 4
th AR model and 1
st order GM model processes have indicated slightly better performance than the GPS-9543 solution. In the 5-sec, 10-sec, 20-sec update cases, the system solutions are showing that the position errors have increased significantly between updates. The biggest position errors have often occurred in the corner sections. This can be because the Yaw rate error is significant in position determination at the corners.
Figures 14 and
15 illustrate the horizontal position trajectory for 1-sec update and 5-sec update intervals, respectively.
Shown in
Figures 16,
17,
18 and
19 are northing and easting position errors from different updating intervals comparing to PPP solutions. It is noticed that the solutions of 4
th order AR model and 1
st order GM model are varying as the updating intervals have increased.
It is clear that the position error in each channel tends to increase without GPS position updates and settle down with updates. Also, the kinematic mode position error is much bigger than the static mode position error in all update interval cases. In the 10-sec updates, the maximum position error has reached to about 60 meters and even worse since 100 meter position error has been shown in the 20 sec updates. The numerical result of position errors of the kinematic testing is summarized in
Table 4.
Based on the position error plots and numerical values in
Table 4, the estimation with 4
th order AR model has produced better results than the estimation with 1
st order GM model in kinematic testing. However, the improvement is relatively smaller than the one in static testing [
8]. It could be explained that the time-correlated behavior in the stochastic variation of MEMS inertial sensors are quite subtle and more unmodeled deterministic error sources are involved in the kinematic environments. This however is dependent on the testing equipments actually used in the testing.
The kinematic testing described here has made very important assumptions mentioned earlier. Therefore, the results with those assumptions and data processing method presented above should be understood very carefully.
6. Conclusions
The main purpose of this paper was to characterize the error behaviour of MEMS based inertial sensors and to analyze the performance of an integrated system prototype of MEMS based inertial sensors and a low-cost GPS receiver with digital compass. Along with deterministic error sources (zero-offset bias and 1st order scale factor), the stochastic variation modeling of MEMS based inertial sensors was focused in this paper due to the assumption of short period of aiding time to GPS. The 1st order Gauss Markov and the 4th order AR models have been used to model the stochastic variation of the MEMS accelerometers of RGA300CA. The static and kinematic tests have shown a better performance using the 4th order AR model. However the test results have also demonstrated that the quality performance of MEMS based inertial sensors (RGA300CA) is not yet quite acceptable to aid GPS system for land vehicle application for longer period of GPS signal outrage. As the performance of MEMS sensors continues to improve, it is recommended to test the latest MEMS sensors available in the market to develop GPS/MEMS based inertial sensor integrated navigation system for continuous navigation solutions with the procedures described here.