Next Article in Journal
Enhancement of Sensitivity and Accuracy of Micro/Nano Water Droplets Detection Using Galvanic-Coupled Arrays
Next Article in Special Issue
Direct Wideband Coherent Localization by Distributed Antenna Arrays
Previous Article in Journal
Causality of Biodiversity Loss: Climate, Vegetation, and Urbanization in China and America
Previous Article in Special Issue
An Unconventional Multiple Low-Cost IMU and GPS-Integrated Kinematic Positioning and Navigation Method Based on Singer Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Odometer Velocity and Acceleration Estimation Based on Tracking Differentiator Filter for 3D-Reduced Inertial Sensor System

College of Automation, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Submission received: 10 September 2019 / Revised: 15 October 2019 / Accepted: 15 October 2019 / Published: 17 October 2019
(This article belongs to the Special Issue Multi-Sensor Systems for Positioning and Navigation)

Abstract

:
Velocity information from the odometer is the key information in a reduced inertial sensor system (RISS), and is prone to noise corruption. In order to improve the navigation accuracy and reliability of a 3D RISS, a method based on a tracking differentiator (TD) filter was proposed to track odometer velocity and acceleration. With the TD filter, an input signal and its differential signal are estimated fast and accurately to avoid the noise amplification that is brought by the conventional differential method. The TD filter does not depend on an object model, and has less computational complexity. Moreover, the filter phase lag is decreased by the prediction process with the differential signal of the TD filter. In this study, the numerical simulation experiments indicate that the TD filter can achieve a better performance on random noises and outliers than traditional numerical differentiation. The effectiveness of the TD filter on a 3D RISS is demonstrated using a group of offline data that were obtained from an actual vehicle experiment. We conclude that the TD filter can not only quickly and correctly filter velocity and estimate acceleration from the odometer velocity for a 3D RISS, but can also improve the reliability of the 3D RISS.

1. Introduction

Most current land vehicular navigation is highly dependent on the Global Positioning System (GPS). However, in urban canyons, tunnels, and other GPS-denied environments, GPS service may suffer from possible signal outages, jamming, and multi-path effects. To maintain positioning availability and accuracy in such cases, GPS is augmented with the inertial navigation system (INS). As a standalone approach, INS is inherently immune to external disturbances and is able to provide continuous navigation solutions with short-term accuracy [1]. Therefore, one of the common solutions for vehicular positioning during GPS outages is to augment GPS with INS [2,3].
For low-cost objectives, instead of integrating GPS with a full inertial measurement unit(IMU) containing three accelerometers and three gyroscopes, the reduced inertial sensor system (RISS) has gained more and more attention. Only one azimuth gyroscope and an odometer or wheel encoders is integrated with GPS, referred to as 2D RISS, to provide 2D positioning solutions in planar environments [1,4,5,6,7]. An integrated RISS/GPS module using a particle filter (PF) was proposed in [8] to provide 2D navigation solutions. A 3D RISS [2,9] composed of a 2D RISS and two horizontal accelerometers is another navigation solution suitable for all wheeled moving platforms, and could obtain the pitch and roll angles of a land vehicle. An enhanced version of PF called the Mixture PF was utilized in [10] to perform the tightly coupled integration of a 3D RISS with GPS. Currently, most studies on RISSs focus on the filters of RISS/GPS integration or on inertial sensor error correction techniques because of the low cost and low precision of the utilized inertial sensors [5,11,12].
An odometer is a completely autonomous device, whose measurement error divergence speed is slower than that of the inertial sensors. Therefore, the introduction of odometer measurement information in a pure inertial navigation system could slow down the divergence speed of INS error to a certain extent without decreasing the autonomous characteristics of the system.
In a 3D RISS, the derivative of velocity information from the odometers is used to calculate pitch and roll angles, and the velocity information from the odometer is used to calculate the position together with the inertial sensors. Additionally, during GPS outages, the velocity information from the odometers is also used for the RISS/GPS integrated navigation to correct the performance degradation with a Kalman filter (KF) [13,14,15]. However, in actual roadways, when the vehicle sideslips or jumps off the ground, the output of the odometer will not represent the actual velocity of the vehicle. Meanwhile, the output of the odometer is easily disturbed by random errors and some signal jumps because of the quantization errors and measurement noises of the wheel encoders. In these conditions, the error of the derivative of velocity information obtained in the classical numerical differential method is significantly boosted, which will decrease the navigation accuracy and reliability of the vehicle greatly. If the wrong velocity information is used in integrated RISS/GPS navigation, the error of the KF will diverge rapidly without restraint, especially during GPS outages. Moreover, the precision of the odometer for measuring the velocity of vehicles is low because of the cost of the odometers and the civilian vehicles overall. Therefore, it is necessary to improve the measurement accuracy and reliability of the odometer using error estimation technologies. Recently, many scholars have proposed various algorithms and approaches to improve this problem based on KFs [16,17]. However, KF and its improving algorithms need to build the object models to extract state information, but the actual velocity error model parameters cannot be measured accurately. In [18], a filter was presented that could achieve high-velocity estimation by fusing information from a magnetometer array with visual-inertial navigation systems (VINS), although the complexity and cost of the whole navigation system would be increased.
In engineering applications [19,20], it is common to encounter the requirements for extracting true signals or their derivatives from signals contaminated with noise. There have been various filtering methods proposed by researchers, such as the high-gain observer-based differentiator [21], adaptive filters [22], linear time-derivative trackers [23], robust exact differentiator [19,24], and so on. Actually, it is difficult to guarantee real-time performance because of the large computation complexity of these complicated filtering technology. Tracking Differentiator (TD), originally proposed by Han Wang [25], is mainly used to solve the problem of reasonably extracting continuous and differential signals from measurement signals with discontinuous or random noises. It overcomes the drawbacks of the classical differential algorithm and has a strong ability to suppress noises. Because of its rigid proof in mathematics and good tracking ability for applications in engineering, it has been widely used in various fields [26,27,28]. In [29], a method of extracting the velocity of a moving vehicle from the output signal of an odometer was introduced using TD to improve velocity precision. Compared with numerical differentiation, this method can not only restrain noise, but it also has a strong anti-interference ability, a simple algorithm, convenient parameter adjustment, and good filter effect. On the other hand, the phase lags and the differential output of the TD filter were not taken into account in this study [29].
Inspired by previous studies, this paper attempts to introduce the tracking-differentiator filter to extract proper velocities and their derivatives from odometers with noises, in order to improve the navigation accuracy and reliability of 3D RISSs in GPS-denied environments. In this 3D RISS, the original odometer velocity is filtered using a TD filter, the tracking signal output of the TD filter is used as the velocity, and the differential output of the TD filter is used as the acceleration. Using a group of offline data obtained from an actual vehicle experiment, the effectiveness of the TD filter on the 3D RISS was demonstrated in simulation experiments.

2. 3D RISS Mechanization and Error Analysis

2.1. 3D RISS Mechanization

A 3D RISS is comprised of one azimuth gyroscope providing the azimuth angular rate change, two horizontal accelerometers for calculating the pitch and roll angles in a horizontal plan, and an odometer with a moving velocity in a near-horizontal plan [2,6,10,30]. The 3D RISS mechanization schematic diagram is shown in Figure 1.
During the 3D RISS mechanization process, the pitch angle of the vehicle is calculated by the forward accelerometer output information after error compensation. The pitch angle calculation formula is yielded by Equation (1):
p = sin 1 ( f y a o d g ) ,
where aod is the acceleration of the vehicle in a near-horizontal plan, and g is the Earth’s gravity.
aod is not directly measurable, and is obtained from the derivative of the odometer velocity νod. Generally, aod at each time epoch can be calculated as
a o d = v o d ( k ) v o d ( k 1 ) d t ,
where νod(k) refers to the νod at each time step, and dt is the sample epoch.
After that, the roll angle of the dynamic vehicle is calculated by the transversal accelerometer information fx, the azimuth gyroscope measurement wz, and the odometer velocity information. Therefore, the roll angle calculation formula is yielded by Equation (3):
r = sin 1 ( f x + v o d w z g cos p ) ,
Simultaneously, the azimuth can be derived by Equation (4).
A ˙ = ( w z w i e sin φ v e tan φ R N + h ) ,
After the attitude calculation, the 3D velocity in the local level frame (LLF) can be derived using Equation (5):
v = [ v e v n v u ] = [ v o d sin A cos p v o d cos A cos p v o d sin p ] ,
Finally, the 3D position can be obtained by integration. The step-by-step computation of the 3D position is yield by
h ˙ = v u , φ ˙ = v n R + h , λ ˙ = v e ( R + h ) cos φ ,

2.2. Errors Analysis

2.2.1. Attitude Errors Analysis

Differentiating Equations (1), (3), and (4), respectively, the attitude errors can be yield
δ p = δ f y δ a o d g cos p ( f y a o d ) δ g g 2 cos p ,
δ r = δ f x + δ v o d w z + v o d δ w z g cos p cos r + ( f x + v o d w z ) ( δ g cos p g δ p sin p ) ( g cos p ) 2 cos r ,
δ A ˙ = ( δ w z w i e δ φ cos φ δ v e tan φ + v e δ φ sec 2 φ R N + h + v e tan φ δ h ( R N + h ) 2 ) ,
From Equation (7), the pitch error is determined by the forward accelerometer measurement error δfy, the odometer’s forward accelerate error δaod, and the Earth’s gravity error δg. The roll error is determined by the transversal accelerometer measurement error δfx, the odometer measurement error δvod, the vertical gyroscope measurement error δwz, the pitch error δp, and also the Earth’s gravity error δg. Assuming that the Earth’s radius R is a relative large value, the azimuth error can be simplified as
δ A ˙ = ( δ w z w i e δ ϕ cos ϕ ) ,
Therefore, the azimuth error is mainly determined by the vertical gyro measurement error δwz and the latitude error δφ.

2.2.2. Velocity Errors Analysis:

Differentiating Equation (5), the velocity errors can be yielded:
δ v e = δ v o d sin A cos p + δ A v o d cos A cos p δ p v o d sin A sin p ,
δ v n = δ v o d cos A cos p δ A v o d sin A cos p δ p v o d cos A sin p ,
δ v u = δ v o d sin p + δ p v o d cos p ,

2.2.3. Position Errors Analysis

Differentiating Equation (6), the position errors can be obtained:
δ h ˙ = δ v u ,
δ φ ˙ = δ v n R + h v n δ h ( R + h ) 2 ,
δ λ ˙ = δ v e + δ φ v e tan φ ( R + h ) cos φ v e δ h ( R + h ) 2 cos φ ,
In above analysis, the νod derived from the odometer is at least as important as the data from the gyroscope or accelerometers. The derivative aod of νod is used to calculate the pitch and roll angles. The νod is used to calculate the 3D position together with inertial sensors. Hence, in actual applications, it is of great significant to improve the measurement accuracy and reliability of the odometer velocity.

3. Tracking Differentiator Filter

3.1. TD Filter Principle

In classical control theory, the differentiator is built using a small time-constant inertial unit. The first-order derivate of the input signal U(s) can be obtained by following a linear, time-invariant, continuous-time dynamic system like Equation (17):
Y = s T s + 1 U = 1 T ( 1 1 T s + 1 ) U ,
where Y(s) and U(s) are the output and input, respectively; T is the time constant; and s is the Laplace operator. In fact, when T is small enough, the inertia unit becomes an approximate time-delay unit. That is to say, 1/(Ts + 1) ≈ e-Ts. The inertia unit in Equation (17) could be treated as a time-delay unit with a small time constant. The inverse Laplace transform of Equation (17) is
y ( t ) 1 T ( u ( t ) u ( t T ) ) u ˙ ( t ) ,
When signals u(t) are corrupted by noises, the noises are also amplified 1/T times; so, the differentiator obtained by the classical method is not suitable for most engineering applications. An improved differentiator of Equation (19) was proposed in [25].
Y ( s ) = 1 T 2 T 1 ( 1 T 1 s + 1 1 T 2 s + 1 ) U ( s ) ,
The difference of the two inertia units is used as the differential to depress the noise amplification. In order to obtain the differential by the fastest dynamic part, a resulting control law that drives any initial state point to the origin in the minimum time is introduced to construct the noise-tolerant time optimal control (TOC)-based TD [27].
The double-integral system is defined as
{ x ˙ 1 = x 2 x ˙ 2 = u ,
where |u| ≤ r, and r is a constant constraint of the control input. It was proven in [24] that the resulting feedback control law that drives the state from any initial point to the origin in the shortest time is
u = R sat ( x 1 v + | x 2 | x 2 2 R , δ ) ,
where ν is the desired value for x1.
sat ( A , δ ) = { sign ( A ) , | A | > δ A / δ , | A | δ , δ > 0 ,
where sign(·) is the Sign function, so the TD filter is constructed by using
{ x ˙ 1 = x 2 x ˙ 2 = R sat ( x 1 v + | x 2 | x 2 2 R , δ ) ,
where x1 is the desired trajectory and x2 is its derivative.
Via Euler’s method discretization, the discrete form of TD is given:
{ f h = fhan ( x 1 ( k ) v ( k ) , x 2 ( k ) , r , h 0 ) x 1 ( k + 1 ) = x 1 ( k ) + h x 2 ( k ) x 2 ( k + 1 ) = x 2 ( k ) + h f h ,
where ν is the input signal, x1 is the filter value of ν, x2 is the derivative of x1, r is the tracking velocity of the TD filter, and h is the step size in simulation. The nonlinear switching function fhan (x1-ν, x2, r, h0) is given by Equation (25):
f h = fhan ( x 1 v , x 2 , r , h 0 ) { d = r h 0 d 0 = h 0 d y = x 1 v + h 0 x 2 a 0 = d 2 + 8 r | y | a = { x 2 + ( a 0 d ) sign ( y ) / 2 , | y | > d 0 x 2 + y / h 0 , | y | d 0 fhan = { r sign ( a ) , | a | > d r a / d , | a | d
The state x1 tracks the input signal ν in the maximum velocity r without oscillation due to the function of fhan. The error between x1 and ν goes to zero. The larger the speed factor r is, the faster the signal tracks. However, the larger of the speed factor r is, the stronger the noise amplification is. Moreover, the noise will reduce by adjusting the filter parameter h0 (5~10 times of h). The larger the filter factor h0 is, the better the filtering effect is. However, the larger the filter factor h0 is, the greater the phase loss of the tracking signal. Therefore, in order to obtain a better filtering effect, coordinated adjustment of r and h0 is required.
Generally, there is some phase lag on the output results of filters. Because TD can give the derivative of the input signal, the phase lag of the TD filter could be compensated with the following equations:
{ v 1 ( k ) = v ( k ) + h 1 x 2 ( k ) f h = fhan ( x 1 ( k ) v 1 ( k ) , x 2 ( k ) , r , h 0 ) , x 1 ( k + 1 ) = x 1 ( k ) + h x 2 ( k ) , x 2 ( k + 1 ) = x 2 ( k ) + h f h , ,
where ν1 is the new input signal composed of the original input signal ν, h1 is the forecast time of ν with x2, and x1 is the filter value of ν1. The forecast time h1 is usually 1~1.5 times of h0.

3.2. TD Filter Simulation Examples

3.1.1. Phase Compensation for Signal Filtering

Let y(t) = sin(20πt) be an original input signal. A TD filter is set with the given design parameters: h = 0.005 s, h0 = 5 h, r = 30,000, and h1 = 1.2 h0. The original input signal y(t), the filtering result of TD without phase compensation and the filtering result of TD with phase compensation are shown in Figure 2a. As Figure 2a shows, the phase lag of the TD filter is obviously reduced by the phase compensation. The numerical differential result of y(t), the differential result of TD without phase compensation and the differential result of TD with phase compensation are shown in Figure 2b. As Figure 2b shows, after half a cycle of y(t), the differential result of TD with phase compensation can track the numerical differential result of y(t) very well. So, in the following part of this paper, TD refers in particular to TD with phase compensation.

3.1.2. Noise Reduction

A random noise with uniform distribution in [0,0.1] is added to y(t). The given design parameters of the above TD filter remain unchanged. The original input signal y(t), the noisy signal and the filtering result of TD are shown in Figure 3a. TD can extract the original signal from the noisy signal with a small phase lag. The numerical differential result of original signal and noisy signal, and the differential result of TD are shown in Figure 3b. Compared to the numerical differential result of the original signal, the differential result of TD is obviously more accurate than the numerical differential result of the noisy signal.

3.1.3. Outliers Exclusion

Here, a hypothetical outlier point is magnified 10x from the noisy signal referred to in Section 3.1.2. In addition, the given design parameters of the above TD filter remain unchanged. The original input signal y(t), the noisy signal and the filtering result of TD are shown in Figure 4a. TD cannot only extract the original signal from the noisy signal with a small phase lag, but is also hardly affected by the outlier point. The numerical differential result of the original signal and noisy signal, and the differential result of TD are shown in Figure 3b. The outlier point leads to the remarkable error of the numerical differential result of the noisy signal. However, there is almost no effect caused by the outlier point on the differential result of TD.
As Section 3 demonstrated, TD achieves a better performance on random noises and outliers than traditional numerical differentiation. Using the odometer velocity data of a 3D RISS as the input signal of a TD filter, the filtered value of the odometer velocity data can be obtained from x1 of this TD with noise reduction and outlier exclusion, and the acceleration data can also be obtained from x2 of this TD filter rather than being calculated using the numerical differential method.

4. Simulation Experiments

In this section, a group of offline data prepared for a study of RISSs by the Navigation Instrumentation Research Group in Royal Military College of Canada is used to demonstrate the effect of velocity and acceleration estimation through an actual vehicular experiment based on TD. The 3D RISS mechanization is constructed from a low-grade Xbow IMU with the odometer output at 10 Hz update rates, and high-end Novatel SPAN IMU mechanization results are used as the reference. The odometer velocity (preprocessed by an offline wavelet filter) of this group is used as the original odometer velocity. The original velocity is filtered through a TD filter set with the given design parameters: h = 0.1 s, h0 = 5 h, r = 20,000, and h1 = 1.4 h0. In the 3D RISS with the original odometer velocity filtered by TD, the x1 takes the place of the original velocity νod and the x2 takes the place of aod calculated by Equation (2). The errors in attitude, velocity, position, and 2D trajectory will be plotted for comparison. Finally, error analysis will be discussed.

4.1. Validation of a 3D RISS with Velocity Filtered by TD

The three navigation simulation results of the 3D RISS with the original odometer velocity without a TD filter, the 3D RISS with the original odometer velocity filtered by TD, and the Novatel IMU are compared in Figure 5, Figure 6, Figure 7 and Figure 8. In Figure 5, Figure 6, Figure 7 and Figure 8, “3D-RISS without TD” refers to the simulation results of the 3D RISS with the original odometer velocity without a TD filter, “3D-RISS with TD” refers to the simulation results of the 3D RISS with the original odometer velocity filtered by TD, and “Novatel” refers to the simulation results of the Novatel IMU. The standard deviations of the main navigation errors between the two 3D RISSs and Novatel IMU are listed in Table 1.
Figure 5a depicts a comparison of the tri-axial attitude curves. The two attitude errors between the two types of 3D RISSs and the Novatel IMU are plotted in Figure 5b. The pitch and roll data of the two types of 3D RISSs have a similar variation tendency with the Novatel attitude reference and no divergence trend. However, the azimuth error curves are divergent with time, which is mainly caused by the vertical gyro measurement error δωz. According to the 3D RISS mechanization, the azimuth error divergence cannot be corrected with the odometer velocity. The standard deviations of the pitch and roll errors between the 3D RISS with the original odometer velocity and the Novatel IMU are 0.88° and 0.32°, respectively. The standard deviations of the pitch and roll errors between the 3D RISS with the original odometer velocity filtered by TD and the Novatel IMU are 1.16° and 0.36°, respectively. Because of the filter phase lag, the horizontal attitude error with the TD filter is a little bigger.
Figure 6a depicts a comparison of the tri-axial velocity curves. The two velocity errors between the two types of 3D RISSs and the Novatel IMU are plotted in Figure 6b. As the azimuth error increases, the east velocity VE error and the north velocity VN error of the two RISSs also increases. While the pitch is not divergent, neither is the error of the vertical velocity VU. The results of the VE, VN, and VU of the two 3D RISS do not coincide with each other. The standard deviations of the VE, VN, and VU errors between the 3D RISS with the original odometer velocity and the Novatel IMU are 2.40 m/s, 2.13 m/s, and 0.28m/s, respectively. The standard deviations of the VE, VN, and VU errors between the 3D RISS with the original odometer velocity filtered by TD and the Novatel SPAN IMU are 2.41 m/s, 2.13 m/s and 0.35 m/s, respectively.
Figure 7a depicts a comparison of tri-axial position curves. The two position errors between the two types of 3D RISSs and the Novatel SPAN IMU are plotted in Figure 7b. As the horizontal velocity error increases, the latitude and longitude errors of the two RISSs also increase. The two altitude errors of the two 3D RISSs have the same divergence tendency from the beginning of the 3D RISS mechanization, which is because of the open loop characteristic of the 3D RISS mechanization in the vertical channel. Similar to the velocity results, the positions of the two 3D RISSs coincide with each other as well. The standard deviations of latitude and longitude errors between the 3D RISS with the original odometer velocity and the Novatel SPAN IMU are 0.029107° and 0.052164°, respectively. The standard deviations of latitude and longitude errors between the 3D RISS with the original odometer velocity filtered by TD and the Novatel SPAN IMU are 0.029109° and 0.052155°, respectively. The position results of the 3D RISS with the original odometer velocity filtered by TD are slightly better than those of the 3D RISS with the original odometer velocity.
The three 2D trajectory simulation results of the 3D RISS with the original odometer velocity, the 3D RISS with the original odometer velocity filtered by TD, and the Novatel IMU are compared in Figure 8. The two trajectory errors between the two types of 3D RISSs and the Novatel IMU are almost identical.
As Figure 5, Figure 6, Figure 7 and Figure 8 and Table 1 show, the two navigation errors of the 3D RISS with the original odometer velocity and the 3D RISS with the original odometer velocity filtered by TD occur at the same level. Using the original odometer velocity, TD could quickly and correctly filter the velocity and estimate the acceleration for a 3D RISS.

4.2. Anti-Interference Ability of a 3D RISS with Velocity Filtered by TD

In order to demonstrate the anti-interference effect of TD on a 3D RISS, a white Gaussian noise (with 0.002 variance and a 0 mean value) and a hypothetical outlier point (2 times amplitude and 100 s interval) are added to the original odometer velocity. The white Gaussian noise model is in accordance with our Kalman filter for the RISS/GPS integration, and the hypothetical outlier point is to simulate the odometer velocity error caused by vehicle sideslips or jumps off the ground. The given design parameters of the above TD filter remain unchanged. The results in Section 4.1 of the 3D RISS with the original odometer velocity are used as the reference. The three navigation simulation results of the 3D RISS with a noisy odometer velocity without a TD filter, the 3D RISS with a noisy odometer velocity filtered by TD, and the 3D RISS with the original odometer velocity are compared in Figure 9, Figure 10, Figure 11 and Figure 12. In Figure 9, Figure 10, Figure 11 and Figure 12, “3D RISS with noise” refers to the simulation results of the 3D RISS with a noisy odometer velocity without a TD filter, and “3D RISS with TD” refers to the simulation results of the 3D RISS with a noisy odometer velocity filtered by TD, and “Original 3D RISS” refers to the simulation results of the 3D RISS with the original odometer velocity. The standard deviations of the navigation errors between the two 3D RISSs with a noisy odometer velocity and the reference are listed in Table 2.
Figure 9a depicts a comparison of the tri-axial attitude curves. The two attitude errors between the two 3D RISSs with noisy odometer velocities and the reference are plotted in Figure 9b. Because the aod is in the calculation of the pitch, the pitch errors are most obviously affected in the comparison of attitude curves. The max pitch error between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity reaches 91°. However, the max pitch error between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity filtered by TD is 11°. The standard deviations of the pitch and roll errors between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity are 10.61° and 0.065°, respectively. The standard deviations of the pitch and roll errors between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity filtered by TD are 1.75° and 0.18°, respectively.
Figure 10a depicts a comparison of the tri-axial velocity curves. The two velocity errors between the two 3D RISSs with noisy odometer velocities and the reference are plotted in Figure 10b. The standard deviations of the VE, VN, and VU errors between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity are 24.95 m/s, 24.19 m/s, and 34.79 m/s, respectively. The standard deviations of the VE, VN, and VU errors between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity filtered by TD are 0.26 m/s, 0.25 m/s, and 0.64 m/s, respectively. Due to the excellent filtering performance of TD in velocity estimation, the velocity results of the 3D RISS with a noisy odometer velocity filtered by TD is hardly influenced by the hypothetical noises.
Figure 11a depicts a comparison of the tri-axial position curves. The two position errors between the two 3D RISSs with noisy odometer velocities and the reference are plotted in Figure 11b. The standard deviations of latitude and longitude errors between the 3D RISS with the original odometer velocity and the 3D RISS with a noisy odometer velocity are 0.0021° and 0.0029°, respectively. The standard deviations of latitude and longitude errors between the 3D RISS with a noisy odometer velocity filtered by TD and the 3D RISS with the original odometer velocity are 0.000039° and 0.000077°, respectively. Similarly, the differences of the above two velocity errors are reflected in the two position errors. The three 2D trajectory simulation results of the three types of 3D RISS are compared in Figure 12.
As Figure 9, Figure 10, Figure 11 and Figure 12 and Table 2 show, with the hypothetical significant noises introduced into the original odometer velocity, the divergence speed of navigation errors of the 3D RISS with odometer velocity filtered by a TD filter is much slower than the pure 3D RISS. Thus, the reliability of a 3D RISS is obviously improved by a TD filter used for estimating velocity and acceleration. Additionally, the velocity filtered by TD is favorable to reduce the error divergence risk of the integrated navigation Kalman filter.

5. Conclusions

Compared to a lot of studies on RISSs that have focused on filters of RISS/GPS integration or on the error correction of inertial sensors, there have seldom been studies on odometer velocity for a 3D RISS. However, in 3D RISS mechanization and error analysis, velocity information and its derivatives are at least as important as the information from the other inertial sensors of the RISS. Velocity information from an odometer is prone to noise corruption, which further leads to the noise amplification of acceleration information in a conventional differential method. This paper has presented a solution for odometer velocity and acceleration estimation using a 3D RISS based on a TD filter.
A TD filter does not depend on an object model and has less computation. With a TD filter, an input signal and its differential signal are estimated fast and accurately. Additionally, using the differential signal output by the TD filter, the filter phase lag can be decreased with the prediction method. As Section 3 demonstrated, TD achieves better performance on random noises and outliers than traditional numerical differentiation.
Using a group of offline data obtained from an actual vehicle experiment, the effectiveness of a TD filter on a 3D RISS was demonstrated through simulation experiments. With the odometer velocity data of a 3D RISS as the input signal of a TD filter, the filtered value with noise reduction and the exclusion of outlying odometer velocity data can be obtained from tracking the signal output from this TD, and the acceleration data can be obtained from the derivative output of this TD also, rather than calculated using the numerical differential method. As the results show, a TD filter could not only correctly and quickly filter the velocity and estimate the acceleration from the odometer velocity using a 3D RISS, but could also improve the reliability of a 3D RISS.
The future work will be to transform the TD filter algorithm into a navigational computer program of an actual 3D RISS.

Author Contributions

Conceptualization and methodology, Q.Z. and L.G.; software, Q.Z.; validation, D.X.; writing—original draft preparation, Q.Z.; writing—review and editing, L.G.

Funding

This research was supported by the Fundamental Research Funds for the Central Universities (3072019CF0405, 3072019CFM0402), the National Natural Science Foundation of China (NSFC) (No. 61803118), the Science and Technology Research Program of Chongqing Municipal Education Commission (KJZD-K201804701), and the Post Doc. Foundation of Heilongjiang Province (LBH-Z17053).

Acknowledgments

All the authors are sincerely thankful to the Navigation Instrumentation Research Group in the Royal Military College of Canada for the experiments and the related data support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shen, Z.; Georgy, J.; Korenberg, M.J.; Noureldin, A. Low cost two dimension navigation using an augmented Kalman filter/Fast Orthogonal Search module for the integration of reduced inertial sensor system and Global Positioning System. Transp. Res. C Emerg. Technol. 2011, 19, 1111–1132. [Google Scholar] [CrossRef]
  2. Georgy, J.; Noureldin, A.; Korenberg, M.J.; Bayoumi, M.M. Low-cost three-dimensional navigation solutionfor RISS/GPS integration using mixture particle filter. IEEE Trans. Veh. Technol. 2010, 59, 599–615. [Google Scholar] [CrossRef]
  3. Noureldin, A.; El-Shafie, A.; Bayoumi, M. GPS/INS integration utilizing dynamic neural networks for vehicular navigation. Inf. Fusion 2011, 12, 48–57. [Google Scholar] [CrossRef]
  4. Iqbal, U.; Karamat, T.B.; Okou, A.F.; Noureldin, A. Experimental results on an integrated GPS and multisensory system for land vehicle positioning. Int. J. Navig. Obs. 2009, 2009, 1–18. [Google Scholar]
  5. Georgy, J.; Noureldin, A.; Korenberg, M.J.; Bayoumi, M.M. Modeling the stochastic drift of a MEMS-based gyroscope in gyro/odometer/GPS integrated navigation. IEEE Trans. Intell. Transp. Syst. 2010, 11, 856–872. [Google Scholar] [CrossRef]
  6. Iqbal, U.; Okou, A.F.; Noureldin, A. An integrated reduced inertial sensor system—RISS/GPS for land vehicle. In Proceedings of the Position, Location and Navigation Symposium (2008 IEEE/ION), Monterey, CA, USA, 5–8 May 2008; pp. 1014–1021. [Google Scholar]
  7. Iqbal, U.; Noureldin, A. Integrated Reduced Inertial Sensor System/GPS for Vehicle Navigation: Multi-Sensor Positioning System for Land Applications Involving Single-Axis Gyroscope-Vehicle Odometer and Integrated with GPS, 1st ed.; VDM Verlag Dr. Müller: Saarbrücken, Germany, 2009; pp. 3–10. [Google Scholar]
  8. Georgy, J.; Iqbal, U.; Bayoumi, M.; Noureldin, A. Reduced inertialsensor system (RISS)/GPS integration using particle filtering for landvehicles. In Proceedings of the ION GNSS 2008, Savannah, GA, USA, 16–19 September 2008; pp. 30–37. [Google Scholar]
  9. Han, S.; Wang, J. Land vehicle navigation with the integration of GPS and reduced INS: Performance improvement with velocity aiding. J. Navig. 2010, 63, 153–166. [Google Scholar] [CrossRef]
  10. Georgy, J.; Noureldin, A. Tightly Coupled Low Cost 3D RISS/GPS Integration Using a Mixture Particle Filter for Vehicular Navigation. Sensors 2011, 11, 4244–4276. [Google Scholar] [CrossRef]
  11. Li, X.; Xu, Q.; Li, B.; Song, X. A Highly Reliable and Cost-Efficient Multi-Sensor System for Land Vehicle Positioning. Sensors 2016, 16, 755. [Google Scholar] [CrossRef]
  12. Georgy, J.; Noureldin, A. Nonlinear filtering for tightly coupled RISS/GPS integration. In Proceedings of the ION PLANS 2010, Indian Wells, CA, USA, 4–6 May 2010; pp. 1014–1021. [Google Scholar]
  13. Chang, T.H.; Wang, L.S.; Chang, F.R. A solution to the ill-conditioned GPS positioning problem in an urbanenvironment. IEEE Trans. Intell. Transp. Syst. 2009, 10, 135–145. [Google Scholar] [CrossRef]
  14. Li, Q.; Chen, L.; Li, M.; Shih-Lung, S.; Nuchter, A. A sensor-fusion drivable-region and lane-detection systemfor autonomous vehicle navigation in challenging road scenarios. IEEE Trans. Veh. Technol. 2014, 63, 540–555. [Google Scholar] [CrossRef]
  15. Rose, C.; Britt, J.; Allen, J.; Bevly, D. An integrated vehicle navigation system utilizing lane-detection andlateral position estimation systems in difficult environments for GPS. IEEE Trans. Intell. Transp. Syst. 2014, 15, 2615–2629. [Google Scholar] [CrossRef]
  16. Liu, S.; Zheng, D.; Li, R. Compensation Method for Pipeline Centerline Measurement of in-Line Inspection during Odometer Slips Based on Multi-Sensor Fusion and LSTM Network. Sensors 2019, 19, 3740. [Google Scholar] [CrossRef] [PubMed]
  17. Mirabadi, A.; Khodadadi, A. Slip and Slide Detection and Compensation for Odometery System, UsingAdaptive Fuzzy Kalman Filter. Sens. Lett. 2009, 7, 84–90. [Google Scholar] [CrossRef]
  18. Caruso, D.; Eudes, A.; Sanfourche, M.; Vissière, D.; Le Besnerais, G. A Robust Indoor/Outdoor Navigation Filter Fusing Data from Vision and Magneto-Inertial Measurement Unit. Sensors 2017, 17, 2795. [Google Scholar] [CrossRef] [PubMed]
  19. Levant, A.; Livne, M. Exact differentiation of signals with unbounded higher derivatives. IEEE Trans. Autom. Control 2012, 57, 1076–1080. [Google Scholar] [CrossRef]
  20. Ang, K.; Chong, G.; Li, Y. PID control system analysis, design, and technology. IEEE Trans Control Syst Technol. 2005, 13, 559–576. [Google Scholar]
  21. Ahrens, J.; Khalil, H. High-gain observers in the presence of measurement noise: A switched-gain approach. Automatica 2009, 45, 936–943. [Google Scholar] [CrossRef]
  22. Zhu, Z.; Xiang, G.; Leilei, C. Analysis on the adaptive filter based on LMS algorithm. Optik 2016, 127, 4698–4704. [Google Scholar] [CrossRef]
  23. Ibrir, S. Linear time-derivative trackers. Automatica 2004, 40, 397–405. [Google Scholar] [CrossRef]
  24. Davila, J.; Fridman, L.; Levant, A. Second-order sliding-mode observer for mechanical systems. IEEE Trans Autom. Control 2005, 50, 1785–1789. [Google Scholar] [CrossRef]
  25. Jingqing, H.; Wei, W. Nonlinear tracking-differentiators. J. Syst. Sci. Math. Sci. 1994, 14, 177–183. (In Chinese) [Google Scholar]
  26. Jinping, F.; Wei, W.; Yu, C. An improved tracking-differentiator filter based on Taylor’s formula. Optik 2018, 158, 1026–1033. [Google Scholar]
  27. Hehong, Z.; Yunde, X.; Longhua, S. High-precision tracking differentiator via generalized discrete-time optimal control. ISA Trans. 2019, 5, 1–8. [Google Scholar]
  28. Liqiang, W.; Hao, L.; Jingqing, H. Study of tracking-differentiator on filtering. J. Syst. Simul. 2004, 16, 651–652, 670. (In Chinese) [Google Scholar]
  29. Yan, C.; Yongyuan, Q.; Qi, Z. Research on Application of Tracking Differentiator in Data Processing of Odometer. J. Projectiles Rockets Missiles Guidance 2010, 30, 52–54. [Google Scholar]
  30. Iqbal, U.; Georgy, J.; Abdelfatah, W.F. Enhancing Kalman Filtering–Based Tightly Coupled Navigation Solution Through Remedial Estimates for Pseudorange Measurements Using Parallel Instrum. Sci. Technol. 2012, 40, 530–566. [Google Scholar]
Figure 1. The schematic diagram of 3D RISS mechanization.
Figure 1. The schematic diagram of 3D RISS mechanization.
Sensors 19 04501 g001
Figure 2. Filtering result and differential result of TD.
Figure 2. Filtering result and differential result of TD.
Sensors 19 04501 g002
Figure 3. Filtering result and differential result of TD with noise.
Figure 3. Filtering result and differential result of TD with noise.
Sensors 19 04501 g003
Figure 4. Filtering result and differential result of TD with random noises and outlier point.
Figure 4. Filtering result and differential result of TD with random noises and outlier point.
Sensors 19 04501 g004aSensors 19 04501 g004b
Figure 5. RISS attitude and attitude errors with the original odometer velocity filtered by TD.
Figure 5. RISS attitude and attitude errors with the original odometer velocity filtered by TD.
Sensors 19 04501 g005
Figure 6. RISS velocity and velocity errors with the original odometer velocity filtered by TD.
Figure 6. RISS velocity and velocity errors with the original odometer velocity filtered by TD.
Sensors 19 04501 g006
Figure 7. RISS position and position errors with the original odometer velocity filtered by TD.
Figure 7. RISS position and position errors with the original odometer velocity filtered by TD.
Sensors 19 04501 g007
Figure 8. 2D trajectories comparison between the two RISSs and the reference.
Figure 8. 2D trajectories comparison between the two RISSs and the reference.
Sensors 19 04501 g008
Figure 9. RISS attitude and attitude errors with a noisy odometer velocity filtered by TD.
Figure 9. RISS attitude and attitude errors with a noisy odometer velocity filtered by TD.
Sensors 19 04501 g009
Figure 10. RISS velocity and velocity errors with a noisy odometer velocity filtered by TD.
Figure 10. RISS velocity and velocity errors with a noisy odometer velocity filtered by TD.
Sensors 19 04501 g010
Figure 11. RISS position and position errors with a noisy odometer velocity filtered by TD.
Figure 11. RISS position and position errors with a noisy odometer velocity filtered by TD.
Sensors 19 04501 g011
Figure 12. 2D trajectory comparison between the two RISSs with noises and the reference.
Figure 12. 2D trajectory comparison between the two RISSs with noises and the reference.
Sensors 19 04501 g012
Table 1. Standard deviations of main navigation errors between the two 3D RISSs and the Novatel IMU.
Table 1. Standard deviations of main navigation errors between the two 3D RISSs and the Novatel IMU.
Error3D RISS with Original Odometer Velocity3D RISS with Original Odometer Velocity Filtered by TD
pitch0.88°1.16°
roll0.32°0.36°
VE2.40 m/s2.41 m/s
VN2.13 m/s2.13 m/s
VU0.28 m/s0.35 m/s
latitude0.029107°0.029109°
longitude0.052164°0.052155°
Table 2. Deviations of the main navigation errors between the two 3D RISSs with noisy odometer velocities and the reference.
Table 2. Deviations of the main navigation errors between the two 3D RISSs with noisy odometer velocities and the reference.
Error3D RISS with Noisy Odometer Velocity3D RISS with Noisy Odometer Velocity Filtered by TD
pitch10.61°1.75°
roll0.065°0.18°
VE24.95 m/s0.26 m/s
VN24.19 m/s0.25 m/s
VU34.79 m/s0.64 m/s
latitude0.0021°0.000039°
longitude0.0029°0.000077°

Share and Cite

MDPI and ACS Style

Zhang, Q.; Guan, L.; Xu, D. Odometer Velocity and Acceleration Estimation Based on Tracking Differentiator Filter for 3D-Reduced Inertial Sensor System. Sensors 2019, 19, 4501. https://0-doi-org.brum.beds.ac.uk/10.3390/s19204501

AMA Style

Zhang Q, Guan L, Xu D. Odometer Velocity and Acceleration Estimation Based on Tracking Differentiator Filter for 3D-Reduced Inertial Sensor System. Sensors. 2019; 19(20):4501. https://0-doi-org.brum.beds.ac.uk/10.3390/s19204501

Chicago/Turabian Style

Zhang, Qing, Lianwu Guan, and Dexin Xu. 2019. "Odometer Velocity and Acceleration Estimation Based on Tracking Differentiator Filter for 3D-Reduced Inertial Sensor System" Sensors 19, no. 20: 4501. https://0-doi-org.brum.beds.ac.uk/10.3390/s19204501

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