Next Article in Journal
Graphene-Based Ammonia Sensors Functionalised with Sub-Monolayer V2O5: A Comparative Study of Chemical Vapour Deposited and Epitaxial Graphene
Next Article in Special Issue
Solving Monocular Visual Odometry Scale Factor with Adaptive Step Length Estimates for Pedestrians Using Handheld Devices
Previous Article in Journal
A Human Activity Recognition Algorithm Based on Stacking Denoising Autoencoder and LightGBM
Previous Article in Special Issue
Characterization of Multipath Effects in Indoor Positioning Systems by AoA and PoA Based on Optical Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Approach to Robust INS/UWB Integrated Positioning for Autonomous Indoor Mobile Robots

1
School of Information Engineering, Henan University of Science and Technology, Luoyang 471023, China
2
School of Communication and Information Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(4), 950; https://doi.org/10.3390/s19040950
Submission received: 15 January 2019 / Revised: 6 February 2019 / Accepted: 19 February 2019 / Published: 23 February 2019
(This article belongs to the Special Issue Sensor Fusion and Novel Technologies in Positioning and Navigation)

Abstract

:
The key to successful positioning of autonomous mobile robots in complicated indoor environments lies in the strong anti-interference of the positioning system and accurate measurements from sensors. Inertial navigation systems (INS) are widely used for indoor mobile robots because they are not susceptible to external interferences and work properly, but the positioning errors may be accumulated over time. Thus ultra wideband (UWB) is usually adopted to compensate the accumulated errors due to its high ranging precision. Unfortunately, UWB is easily affected by the multipath effects and non-line-of-sight (NLOS) factor in complex indoor environments, which may degrade the positioning performance. To solve above problems, this paper proposes an effective system framework of INS/UWB integrated positioning for autonomous indoor mobile robots, in which our modeling approach is simple to implement and a Sage–Husa fuzzy adaptive filter (SHFAF) is proposed. Due to the favorable property (i.e., self-adaptive adjustment) of SHFAF, the difficult problem of time-varying noise in complex indoor environments is considered and solved explicitly. Moreover, outliers can be detected and corrected by the proposed sliding window estimation with fading coefficients. This facilitates the positioning performance improvement for indoor mobile robots. The benefits of what we propose are illustrated by not only simulations but more importantly experimental results.

1. Introduction

The positioning technology plays an important role in navigation for mobile robots and can be classified into the absolute positioning and relative positioning [1,2,3]. The former generally depends on reference nodes with known absolute locations, which is adopted by global navigation satellite system (GNSS) [4]. However, GNSS signals are often unreliable or even unavailable because they are difficult to be detected by mobile robots all the time in indoor environments [5]. Moreover, the stability and accuracy of GNSS signals are not guaranteed. For these reasons, ultra wideband (UWB) is widely applied for indoor positioning in recent years due to the high resolution, strong anti-jamming performance, low transmit power, and so on [6,7], but it may be affected by multipath effects and the non-line-of-sight (NLOS) in complicated environments, leading to large positioning errors [8]. Unlike the absolute positioning, the relative positioning technology tracks target state incrementally by estimating changes of relative pose. So it is also called dead reckoning (DR) [9]. As a typical application of the relative positioning, inertial navigation system (INS) is usually applied for trajectory tracking of moving objects owing its favorable properties (e.g., high frequency, fast response and high short-term positioning accuracy) [10]. Unfortunately, the positioning errors of INS are accumulated gradually with time and distance, which makes the estimated trajectory deviate from the true one [11]. Thus, INS is not suitable for accurate long-distance indoor navigation and positioning.
With the development of indoor positioning technologies, the INS/UWB integrated positioning system has been extensively studied and widely applied [12,13]. The inertial measurement unit (IMU) is used to obtain the continuous and regular dynamic information. Meanwhile, the UWB location module outputs the discrete positions with noncumulative error [14]. These two heterogeneous data complement and benefit each other in a fusion system, which overcomes the limitation of single sensor in information utilization. However, in practical applications of INS/UWB integrated positioning, UWB data are usually affected by the time-varying noise in complicated indoor environments [15]. Moreover, outliers may appear in measurements. As a result, not only is the optimal estimation lost, but also the divergence of filter occurs [16], which causes large positioning errors especially when measured data is greatly influenced by external environment.
In practical integrated positioning systems, most kinematics/measurement models are complex and nonlinear. Although nonlinear filters such as the extended Kalman filter (EKF) [17,18], unscented Kalman filter (UKF) [19] and cubature Kalman filter (CKF) [20,21] are widely applied to fuse heterogeneous data, their filtering performance may be reduced when there exists the time-varying measurement noise (caused by multipath effects and the NLOS factor) in complicated indoor environments. This is because these filters usually assume the measurement noise is time invariance, which is clearly not in accordance with practical statistics. In view of this, several adaptive filters have been proposed to improve the integrated positioning performance by estimating the statistical characteristics of noise, e.g., innovation adaptive estimation (IAE) and multiple model adaptive estimation (MMAE) [22]. The IAE method may estimate the noise covariance adaptively by analyzing the filter’s innovation sequence and the attenuation factor was introduced [23]. In [24], the fuzzy inference system was used to output the adaptation value of the measurement noise covariance (MNC), which contributes to achieving adaptive estimation for time-varying noise. [25] proposed an enhanced innovation adaptive estimation Kalman filter (EIAE-KF). It can be switched between the prediction model and the update model based on an innovation sequence for getting the optimal estimation. To deal with the influences of unexpected dynamic errors, [26] presented the strong tracking filter (STF) in integrated positioning system. In [27], a simplified Sage–Husa adaptive filter (SHAF) is used to weaken the effect of time-varying noise on filter performance. Although these approaches may improve the integrated positioning performance to some extent, they do not consider the case that outliers occur in measurements.
Thus, [28] established a fuzzy logic system to recognize outliers, and [29] used the chi-square test to judge whether the measured data were outliers or not. Ref. [30] constructed the probability density function (PDF) of different status based on the Kullback–Leibler (KL) divergence theory, and the outlier detection experiments with UWB radar were carried out as well. In [31], based on the orthogonal properties of innovation sequence, an approach was proposed to detect outliers. Thus measurements containing outliers are corrected by using weight coefficients.
Note that above approaches need to assume that the noise covariance changes little or the number of outliers is small. However, these assumptions do not necessarily coincide with the reality, because phenomenons such as wireless signal interference, pedestrians walking and obstacles block may actually exist in complicated indoor environments. These factors make the statistical characteristics of measurement noise change dramatically, but more importantly, more outliers appear in measured data and degrade the performance.
To solve these problems simultaneously, this paper proposes a Sage–Husa fuzzy adaptive filter (SHFAF) for INS/UWB integrated positioning. The error state is modeled for improving the system robustness instead of using the nominal state directly [32]. The IAE method is also considered, and the regulatory factor of the modified innovation contribution weight (ICW) is obtained through fuzzy inference system, so that the MNC can be estimated accurately in real time. Moreover, the sliding window estimation with fading coefficients is introduced to approximate the true innovation covariance. More importantly, it is applied to outliers detection and correction for providing precise positioning results based on the innovation orthogonality.
Compared with existing approaches, this work has the following innovative aspects:
  • An effective framework of INS/UWB positioning system is proposed for autonomous indoor mobile robots in which the proposed modeling approach is simple to implement. This facilitates the design of a Sage–Husa fuzzy adaptive filter and accordingly improves the positioning accuracy and robustness.
  • In our approach, a modified innovation contribution weight is obtained by the fuzzy inference system, which serves our Sage–Husa fuzzy adaptive filter. Thanks to its favorable property (i.e., self-adaptive adjustment), the measurement noise covariance can be estimated accurately.
  • To reduce the adverse effect of outliers in a complicated indoor environment, a sliding window estimation with fading coefficients is presented for outlier detection and correction based on innovation orthogonality criterion.
  • Instead of merely using numerical simulations, we adopt the Gazebo simulation engine to evaluate the positioning performance. The measurement data are generated in real time and complicated factors (e.g., kinematics and interferences) in real scenarios are considered explicitly. More importantly, in this paper, practical experiments are also performed to verify the effectiveness of what we propose.
This paper is organized as follows. Section 2 formulates the INS/UWB integrated positioning problem. Besides, the kinematic and measurement model of INS/UWB are presented, respectively. In Section 3, an effective system framework of INS/UWB integrated positioning is established, where the Sage–Husa fuzzy adaptive filter is easily obtained to improve the positioning performance, the difficult but important problems of time-varying noise and outliers are considered and solved. In Section 4, simulation and experimental results demonstrate the effectiveness of the proposed INS/UWB integrated positioning approach. The last section concludes the paper.

2. System Model

2.1. Problem Formulation

For the problem of INS/UWB integrated positioning, the key is how to integrate two heterogeneous sensors’ data effectively and reasonably according to the prior information such as kinematic model, measurement model and noise characteristics of moving target, so as to achieve the estimated pose state (i.e., x p o s e ) with smaller error than positioning approaches with an individual sensor’s data. It is given by x p o s e = [ p T , q T ] T , which consists of the position and attitude (i.e., the orientation of a moving target). This paper focuses on the target pose state estimation in three-dimensional (3D) cartesian coordinates, in which p and q denote position and attitude in 3D space, respectively. Unlike in the 2D space, the pose in 3D space needs more variables to be described. In addition, the attitude vector and the dynamics parameters are nonlinearly coupled (to be demonstrated later), which is difficult to be handled.
To improve the pose estimation performance, other indirect state variables related to the pose are also considered. They are included into the whole state vector in x k = [ p k T , v k T , q k T , a b , k T , ω b , k T ] T , where k denotes the time index, v , a b and ω b are the velocity, the offset of acceleration and angular velocity. Consider the following system model:
x k = f ( x k 1 , u k 1 , w k 1 ) z k = h ( x k , v k )
It characterizes the target’s kinematics and sensor measurements, where x k is the state vector with transition function f ( · ) , z k is measurement vector obtained by the measurement function h ( · ) . u is the input or control vector, w and v are the process and measurement noise, respectively.
In INS/UWB integrated positioning system, f ( x k 1 , u k 1 , w k 1 ) is determined by the kinematics of INS, and h ( x k , v k ) is obtained by UWB measurement equations. The positioning systems of INS and UWB can be integrated in a unified framework and described by Equation (1). Besides, based on the error state, the kinematic model of INS and UWB measurement model will be established and introduced below.

2.2. Kinematic Model

For INS/UWB integrated positioning, the dynamics of a moving target should be analyzed to obtain its state transition equation, which facilitates the kinematic model. Based on the kinematic model, the relative position, velocity and attitude can be incrementally estimated by using the inertial data (i.e., acceleration and angular velocity). In INS, there are two reference frames: body reference frame {B} and navigation frame {N}. The transformation between these reference coordinates can be shown in Figure 1, in which the transformation from {B} to {N} is described by the translation vector d B N and the rotation matrix C B N jointly.
Compared with Euler angles, the quaternion may not be affected by the gimbal lock and has one-to-one correspondence with the rotation [33]. Thus it is adopted for attitude representation in state vector. Besides, the rotation matrix is also used in this paper, which can be transformed by the quaternion with the following equation [34]:
C = ( 2 q 0 2 1 ) I 3 2 q 0 [ q v ] × + 2 q v q v T ,
where [ · ] × is the skew symmetric operator, q 0 and q v = [ q 1 , q 2 , q 3 ] T are the real and imaginary part of the rotation quaternion, respectively.
The joint state vector of INS/UWB positioning system in the navigation coordinates {N} is x = [ p T , v T , q T , a b T , ω b T ] T . It consists of the nominal state x ^ and the error state δ x , which satisfies
x = x ^ δ x ,
where δ x = [ δ p T , δ v T , δ θ T , δ a b T , δ ω b T ] T . ⊕ denotes a generic composition (quaternion products for the quaternion states and sums for other states). Note that δ θ is the angle error vector respecting to the three axes of body coordinates {B}. Since δ θ is very small, the rotation error of quaternion can be approximately expressed as δ q = [ 1 , δ θ T / 2 ] T . Correspondingly, the transition of the nominal state and error state has the following forms [32]:
x ^ ˙ = p ^ ˙ v ^ ˙ q ^ ˙ a ^ ˙ b ω ^ ˙ b = v ^ C ^ ( a m a ^ b ) + g 1 2 q ^ ( ω m ω ^ b ) 0 0
δ x ˙ = δ p ˙ δ v ˙ δ θ ˙ δ a ˙ b δ ω ˙ b = δ v C ^ [ a m a ^ b ] × δ θ C ^ δ a b C ^ a n [ ω m ω ^ b ] × δ θ ω w ω n a w ω w ,
where g is the gravity acceleration vector in {N}, a m and ω m are the acceleration and angular velocity measurements, respectively. a n and ω n are the measurement noise of the accelerometer and gyroscope, respectively. a w and ω w are the random walk noise of biases.

2.3. Measurement Model

Generally, the processed UWB measurements provide an absolute reference for the integrated positioning system and are used to correct the accumulative error of INS, which are modeled to the position of UWB tag. To improve the positioning accuracy, several UWB anchors may be used sometimes. However, overdetermined equations also may appear in the process of dealing with range measurements obtained by using time of arrival (TOA). So in this paper, the least squares (LS) method was used for more accurate position estimation of UWB tag, which is demonstrated as follows.
Assuming that the position vector of UWB anchors in navigation coordinates is p n = [ x n , y n , z n ] T , where n { 1 , 2 , 3 , , N } denotes the anchor’s serial number, and N is the total number of anchors. If the position vector of the UWB tag at time k is p m , k = [ x m , k , y m , k , z m , k ] T , then the measurement distance from the UWB tag to the n-th anchor is written as
d n , k = p n p m , k 2 = ( x n x m , k ) 2 + ( y n y m , k ) 2 + ( z n z m , k ) 2 .
Similarly, if the true position of tag is p t , k , the true distance can be defined as
d ¯ n , k = p n p t , k 2 = ( x n x t , k ) 2 + ( y n y t , k ) 2 + ( z n z t , k ) 2 .
The positioning problem of the tag can be formulated as finding the solution p ^ m , k which minimizes the the error between the measurement distance and true distance, i.e.,
p ^ m , k = arg min p m , k n = 1 N ( d n , k d ¯ n , k ) 2
Note that N may be greater than the dimension of the unknown position vector. From Equation (6), we have
d 1 , k 2 = ( x 1 x m , k ) 2 + ( y 1 y m , k ) 2 + ( z 1 z m , k ) 2 d 2 , k 2 = ( x 2 x m , k ) 2 + ( y 2 y m , k ) 2 + ( z 2 z m , k ) 2 d N , k 2 = ( x N x m , k ) 2 + ( y N y m , k ) 2 + ( z N z m , k ) 2 .
Let the second row, third row, …, and N-th row subtract the first row, respectively. A linear equation is obtained as
2 G p m , k = b k ,
where
G = x 1 x 2 y 1 y 2 z 1 z 2 x 1 x N y 1 y N z 1 z N
b k = d 2 , k 2 d 1 , k 2 + x 1 2 x 2 2 + y 1 2 y 2 2 + z 1 2 z 2 2 d N , k 2 d 1 , k 2 + x 1 2 x N 2 + y 1 2 y N 2 + z 1 2 z N 2 .
Accordingly, Equation (10) can be solved by using the least squares method:
p m , k = 1 2 ( G T G ) 1 G T b k .
Meanwhile, the velocity of UWB tag is also obtained as
v m , k = Δ p k Δ t k = p m , k p m , k 1 t k t k 1 .
Suppose that the position and velocity vector of INS at time k are p I N S , k and v I N S , k respectively, then the measurement equation of the integrated positioning system is
z k = δ p k δ v k = p m , k p I N S , k v m , k v I N S , k = H δ x k + v k ,
where v N { 0 , R } is the measurement noise vector, R is its covariance matrix, and H = I 3 0 0 0 0 0 I 3 0 0 0 is the measurement matrix of the error state.

3. INS/UWB Integration Strategy with Sage–Husa Fuzzy Adaptive Filter

To achieve precise positions of moving targets using two different types of data (measured from IMU and UWB), this paper proposes a framework of INS/UWB integrated positioning system. As shown in the framework of Figure 2, ω , α , φ 0 are measured by IMU, which serve as the inputs of INS. Thus the position, velocity and attitude of a robot can be obtained by using recursive integration, in which the measurements of triaxial magnetometer are used for its orientation initialization. Absolute measurements (i.e., position and velocity) in navigation coordinates are obtained from UWB positioning system by using the TOA method based on the distances between UWB anchors and tag. In this INS/UWB integrated positioning system, UWB measurements are used to correct the INS output. However, since indoor environments are often complicated (dynamic, occlusion factors, etc.), the UWB measurements are usually affected by the multipath effects and NLOS factor. This may bring the outliers and therefore cause the performance degradation or even filter divergence.
In view of the above, this paper proposes a Sage–Husa fuzzy adaptive filter (SHFAF) based on the error state. It has a favorable property that the nonlinear nominal state estimation problem can be transformed into a optimal linear one [35]. More importantly, the estimation accuracy and robustness improved with low computational complexity.
As shown in Figure 2, SHFAF included the following modules: outlier detection and correction, MNC adaptive estimator, and error state Kalman filter. The module of outliers detection and correction aims to recognize outliers and reduce the deviation degree of them by pre-processing the raw measurements. In MNC adaptive estimator, the time-varying MNC can be estimated accurately based on the innovation adaptive estimation. Moreover, the corrected measurements and their estimated noise covariance are as the inputs of error state Kalman filter. Finally, the compensation values calculated by SHFAF were injected into the outputs of INS, and the pose estimation of robot can be achieved. The implementation of proposed SHFAF is explained in detail below.

3.1. Prediction

From Equation (3), the whole state vector is composed of the nominal state and error state, which can be predicted by the transition equation recursively. The nominal state prediction is written as
x ^ k | k 1 = f ( x ^ k 1 , u k 1 ) ,
where f ( · ) denotes the nonlinear state transition function. It can be derived by the Euler integrals of Equation (4). As the control inputs, u k = [ a m , k T , ω m , k T ] T are measurements of the triaxial accelerometer and gyroscope. The process noise in error state transition is accumulated over time, and the uncertainty of error state can be reflected by a covariance matrix. Correspondingly, the prediction of the error state and its covariance are
δ x k | k 1 = F δ x , k 1 ( x ^ k 1 , u k 1 ) · δ x k 1
P k | k 1 = F δ x , k 1 P k 1 F δ x , k 1 T + Γ n Q n Γ n T .
Note that since the mean of the error state δ x is initialized to zero, the linear Equation (17) always returns zero. F δ x , k 1 is the error state transition matrix, it can be obtained by the discretization of Equation (5):
F δ x , k 1 = I 3 I 3 Δ t 0 0 0 0 I 3 C ^ k 1 [ a m , k 1 a ^ b , k 1 ] × Δ t C ^ k 1 Δ t 0 0 0 C ^ k 1 { ( ω m , k 1 ω ^ b , k 1 ) Δ t } 0 I 3 Δ t 0 0 0 I 3 0 0 0 0 0 I 3
Γ n = 0 0 0 0 I 3 0 0 0 0 I 3 0 0 0 0 I 3 0 0 0 0 I 3
Q n = ( σ a n Δ t ) 2 I 3 0 0 0 0 ( σ ω n Δ t ) 2 I 3 0 0 0 0 ( σ a ω Δ t ) 2 I 3 0 0 0 0 ( σ ω ω Δ t ) 2 I 3 ,
where Γ n is the noise driven matrix and Q n is the process noise covariance matrix.

3.2. Covariance Adaptive Estimator of Measurement Noise

The adaptive Kalman filter is usually applied in complex and unknown noise circumstances, because it estimates the noise covariance at each recursive step. In the traditional Kalman filter framework, however, the noise covariance is a constant matrix. Thus the innovation ε k is defined as
ε k = z k H δ x k | k 1 ,
and its theoretical covariance is
S k = H P k | k 1 H T + R k .
The theoretical covariance equals the truth when the filter is strictly converged. Correspondingly, the convergence criterion is derived as
ε k ε k T = H P k | k 1 H T + R k ,
and the theoretical estimate of MNC matrix R ^ k is
R ^ k = ε k ε k T H P k | k 1 H T .
When the filter runs, the MNC matrix is estimated recursively, and smooth operation is necessary to ensure its continuity. To synthesize the above analysis, the Sage–Husa fuzzy adaptive estimator can be represented by
R k = ( 1 s k α d k ) R k 1 + s k α d k ( ε k ε k T H P k | k 1 H T ) .
In this equation, d k is the modified innovation contribution weight (MICW)
d k = ( λ b ) / ( λ b k + 1 ) ,
where λ 1 . d k reduces to the traditional ICW when λ = 1 . The value λ is used to increase the asymptotic stability value of ICW which provides a wide range for the output of the fuzzy inference system, b is a forgetting factor (usually between 0.95 and 0.99), and s k is a regulatory factor of MICW which can be obtained by the fuzzy inference system.
Remark 1.
The estimated MNC at time k (i.e., R k ) is weighted by two parts; one is the covariance of previous sampling time (i.e., R k 1 ), the other part is the difference between the estimated and theoretical covariance of innovation at time k. Since the traditional ICW is very small after stabilization, the information provided by innovation is little as well. When time-varying noise is observed, the changes of MNC may not be estimated in time. Therefore, the traditional ICW is modified as Equation (27), and it will be analyzed in detail later.
In this paper, we have made two improvements to ICW; (a) λ in Equation (27) was introduced to increase the steady-state value of ICW, and the contribution of innovation to R k was increased, and (b) the regulatory factor (i.e., s k ) of MICW was obtained through fuzzy inference for improving the estimation ability to R k .
We define the adaptive modified innovation contribution weight (AMICW) as the product of regulatory factor and MICW (i.e., s k d k ), then the changes of ICW and AMICW are compared in Figure 3. It can be seen that AMICW always maintained a larger scale and can be adjusted adaptively according to the changes of external noise characteristics. In comparison, the value of ICW was very small and basically invariant after several iterations, which made little contribution to the estimate of R k .
To obtain the regulatory factor s k , we construct a fuzzy reference system. As the input of a fuzzy reference system, the difference coefficient r k between the estimated innovation covariance S ^ k and the theoretical covariance S k should be calculated first, i.e.,
S ^ k = E ( ε k ε k T )
r k = Tr ( S ^ k ) Tr ( S k ) 1 ,
where S k is calculated through Equation (23), and T r ( · ) is the trace of a matrix. Then, the difference coefficient r k serves as the input of the fuzzy reference system, and the corresponding output s k can be obtained:
s k = 1 , k k s f u z z y ( r k ) , k > k s ,
where f u z z y ( · ) is regarded as a nonlinear function. Because of the unknown a priori information of measurement noise, an initial estimate of R is obtained by the simplified Sage–Husa adaptive filter at the beginning (i.e., k k s ). When k > k t , MICW tends to be stable. At this stage, the regulatory factor s k was inferred by the fuzzy inference system and plays a dominant role in adaptive tuning.
Remark 2.
The first adaption for MICW relies on the forgetting factor b, and it is adjusted adaptively in the second time through the regulatory factor s k . As a scaling of s k , α is chosen depending on the practical situations. A larger α approaches a true R by using fewer iterations, but this way leads to unstable estimates (e.g., the estimation oscillates around the ground truth); a smaller α may smooth the estimator, while more steps are necessary to reach the ground truth.
The inference rules of the above fuzzy logic system are written as
I f r k Equal , t h e n s k Equal I f r k More , t h e n s k More I f r k Less , t h e n s k Less
According to Equation (29), it can be seen that r k intuitively reflects the error between the estimated and theoretical value of innovation covariance, which characterizes the change degree of MNC. The larger the r k is, the greater the error between the estimated and theoretical value is. This means the external noise statistical characteristics have changed a lot, and correspondingly more extra information is required to correct R k . Conversely, the smaller the r k , the smaller the contribution of innovation to the estimated MNC. In this case, the present R k should be estimated by more information of the previous R k (i.e., slightly decreasing s k ). Thus s k should be increased to enhance the innovation contribution for the estimate of R k . The input and output membership function of fuzzy inference system are shown in Figure 4.
Through two stages of adaptive estimation, the MNC can be accurately estimated in the whole process, and the performance of the filter can be improved.

3.3. Outlier Detection and Correction

UWB may affected by unknown and uncertain disturbances such as multipath effects and NLOS factor. In this case, the outliers exist in measurements and are far away from the ground truth. This phenomenon whose adverse impact on the filter performance is difficult to be reduced only by IAE method. To avoid large estimation error and even the filter divergence, this paper proposes approach to the outliers detection and correction based on the innovation orthogonal criterion. According to the orthogonality of innovation, we have
E ( z k z k T ) = E ( ε k ε k T ) + E ( z k | k 1 z k | k 1 T ) ,
and the theoretical expectation of z k z k T is
D k = H P k | k 1 H T + R k + H δ x k | k 1 δ x k | k 1 T H T .
It can be judged whether the measurement at time k is an outlier or not by calculating the ratio between the estimated and theoretical expectation of z k z k T . Suppose that the i-th diagonal element of E ( z k z k T ) and D k are G i , k and D i , k , respectively. If M i , k = G i , k / D i , k , then
f r i , k = 1 1 / M i , k   M i , k ξ i M i , k > ξ i ,
where ξ i is the predetermined sensitivity threshold of outliers recognition. Thus the outlier is corrected as
z k = f r k × z k ,
where f r k = d i a g [ f r 1 , k , f r 2 , k , , f r r a n k ( R ) , k ] . From Equation (34) we can see that the larger M i , k (i.e., the deviation between measurement and prediction) is, the smaller f r i , k .
To approximate the true E ( z k z k T ) in Equation (32), the estimated innovation covariance S ^ k and the expectation E ( z k | k 1 z k | k 1 T ) need to be calculated. The value S ^ k is usually generated by sliding window estimation, in which the innovation sequence consisting of the recent N innovations is used to approximate the true covariance at the current time. However, the statistical characteristics of innovations obtained by the traditional sliding window estimation may only reflect their average status. Once outliers appear, the covariance may not vary greatly (especially when the size of sliding window is large), causing the outliers not to be captured correctly and timely.
Remark 3.
Reducing the width of sliding window can enhance the tracking sensitivity to the innovation covariance, but the estimate accuracy of innovation covariance may be degraded if the window is getting too narrow.
To improve the ability to capture outliers, the coefficient series with the ‘forgotten’ properties are weighted to the innovation sequence. This makes new data play a major role (i.e., given by a larger weight), and the innovations far away from filtering moment are conferred smaller weights. As a result, the sensitivity to outliers can be improved while maintaining the width of the sliding window unchanged. The sliding window estimation is modified as
S ^ k = E ( ε k ε k T ) = j = k l + 1 k σ j ε j ε j T ,
where σ j = a k j ( 1 a ) / ( 1 a l ) is the fading coefficient of the j-th innovation covariance in a sliding window with the width of l, and a characterizes the fading rate (usually taken to be between 0.95 to 0.99).

3.4. Measurement Update

After the above steps, the Kalman gain K of the error state δ x k can be calculated, i.e.,
K k = P k | k 1 H T ( H P k | k 1 H T + R k ) 1 .
The estimated error state vector and its corresponding covariance are updated as
δ x k = K k ε k
P k = ( I K k H ) P k | k 1 .
Thus the error state is injected into the nominal state, and the true state is estimated by Equation (3). Particularly,
p ^ k v ^ k q ^ k a ^ b , k ω ^ b , k = p ^ k 1 + δ p k v ^ k 1 + δ v k q ^ k 1 q { δ θ k } a ^ b , k 1 + δ a b , k ω ^ b , k 1 + δ ω b , k .
Afterwards, the estimator (i.e., x ^ k ) is considered as the nominal state for next recursion. Note that the estimate of the error state should be reset to zero after nominal states have been updated. The flow chart of our algorithm is shown in Figure 5.

4. Experiment Analysis and Evaluation

In this section, the proposed INS/UWB integrated positioning approach is evaluated and analyzed to illustrate its effectiveness; a) two simulation scenarios (only the time-varying noise and the time-varying noise with outliers in measurements) are provided to test the convergence and robustness of our approach and b) the experimental platform of INS/UWB integrated positioning system is established, where the practical positioning experiments in indoor environments with two different motion trajectories of a mobile robot are performed to verify the effectiveness of the proposed approach in real applications.
Note that several different positioning approaches are compared in our simulations and practical experiments, including TOA based UWB positioning (see Section 2.3), ESKF [32], SHAF [36] and SHFAF.

4.1. Simulation Experiments and Analysis

In this part, simulations are performed by the Gazebo physical simulation engine. The east–north–up coordinate system is adopted for our navigation, and four UWB anchors are used and located in ( 5 , 1 ), ( 1 , 1 ), ( 1 , 5 ), ( 5 , 5 ) (see Figure 6a). UWB tag is fixed on the geometric center of a robot body. The UWB measurement noises are generally assumed to consist of the ranging noise and the time-varying noise [37]. The later is usually caused by NLOS scenarios, multipath effects and other uncertain factors. The ranging noise is zero-mean Gaussian white noise and its variance is 0.1 2 . The interference noise is random with uniform distributed (i.e., u ( 0.2 , 0.2 ) ). The update frequency of UWB data is 5 Hz, and the TOA-based positioning is used to get the tag position. The IMU and UWB tag are placed in the same location for collecting the acceleration and angular velocity of the robot in simulation environments. The update frequency of IMU data is 100 Hz, and Gaussian noise parameter is 0.06 2 .
In simulations, the initial position of this robot is located at ( 0 , 0 ) in the navigation coordinates, it moves along the diagonal line of 4 m × 4 m square with constant speed 0.4 m/s. Note that there is a relative transformation between the initial body coordinates and the navigation coordinates, which can be seen in Figure 6b. As the robot moving on, IMU and UWB measurements are recorded and processed simultaneously. Meanwhile, several different positioning approaches are compared for performance evaluation.
To verify the ability of self-adaptation and outliers suppression of the proposed INS/UWB integrated positioning approach, two different scenarios are designed as follows.
Scenario 1: The measurement noises are composed of the Gaussian white noise and random walk noise (without outliers). The aim of this experiment is to test the convergence and adaptive capacity of the presented approach. The filter is initialized at ( 0.57 , 0 ), i.e., it starts from a larger error point due to lack of information. In addition, the MNC is initialized with a large matrix, which means that there exists great uncertainty in the initial stage.
Simulation results for this scenario are show in Figure 7, where the absolute error is chosen as a measurable indicator for position accuracy. Figure 7a gives the estimated trajectories without considering outliers. Clearly, SHFAF outperforms ESKF because it has a stronger adaptive ability to the time-varying noise and its performance is not affected by the initial condition of MNC. This leads to poor positioning accuracy, and even worse than the results of UWB-only positioning at the initial stage. From Figure 7b we can see that the absolute positioning error of SHFAF is less than 0.2 m after 2 s, while ESKF needs about 10 s to achieve such an effect. It can be concluded that SHFAF has faster convergence speed and better self-adaptive ability than ESKF.
Scenario 2: Unlike in Scenario 1, outliers exist in measurements with Gaussian white noise and random walk noise. This scenario is designed to evaluate outlier suppression and adaptive ability of the proposed approach comprehensively.
Figure 8 shows that the position is not estimated accurately only by using UWB when outliers exist. Meanwhile, the positioning error can be reduced effectively by using INS/UWB integrated approach. That is, the position estimation of SHFAF is little affected because the outliers processing and MNC adaptive estimator are implemented together. SHFAF has good accuracy and its positioning errors are mainly below 0.2 m. Conversely, ESKF causes large errors when meeting outliers because it lacks of the ability of adaptive estimation to MNC. Compared with ESKF, SHAF has the self-adaptive ability to the time-varying noise so that smaller positioning error can be achieved. Though SHAF outperforms ESKF, it is still difficult to apply it in complicated environments because of low positioning accuracy.
To evaluate the performance of the proposed SHFAF from a quantitative view, the root mean square error (RMSE) of robot positions in the east, north and absolute distance are calculated, which are listed in Table 1. Clearly, SHFAF has the lowest positioning error among all four approaches in this table. Simulation results demonstrate the benefit of the proposed positioning approach—SHFAF has strong robustness and good self-adaptive ability.

4.2. Experimental Analysis and Performance Evaluation

To further evaluate the performance of the proposed positioning approach in real scenarios, we build a INS/UWB integrated positioning platform (see Figure 9), and the technical specifications of sensors are listed in Table 2. Note that the inertial data are collected by 9-DOF Razor IMU of Sparkfun. The IMU is integrated with a three-axis gyroscope, accelerometer and magnetometer, which serves to measure the nine-DOF inertial information. Besides, this paper used the TREK1000 indoor positioning suite of Decawave based on the DW1000 chip. It consists of four UWB anchors and one UWB tag. The 3D position of a tag can be solved by this system which is considered as an independent positioning system. Communications between UWB modules are built on channel 2 in the frequency range from 3.74 GHz to 4.24 GHz, while the router works on the 2.4 GHz frequency band. So in this case, these two wireless networks hardly interfere with each other.
A large experimental site is obviously beneficial, so we make the robot cover space as large as possible under the limited experimental conditions. The configuration of the experimental site refers to [11,38], which is shown in Figure 10. The navigation coordinates is established in a rectangular area (i.e., 6.4 m × 4.8 m) which is surrounded by four UWB anchors. They are located at ( 0 , 0.5 ), ( 6.4 , 0.5 ), ( 0 , 4.3 ) and ( 6.4 , 4.3 ), respectively. All anchors and tag work in the same frequency band and communicate with each other. The sampling frequency of IMU data is 50 Hz, and the update frequency of UWB positioning data is 5 Hz.
In a practical experiment, we designed two different motion trajectories to test our positioning approach (see Figure 11). In trajectory A, a robot started from point ( 0 , 2 ) and headed to the north. This trajectory was an S-shaped curve consisting of two semicircles with 1.5 m radius. The linear velocity was 0.405 m/s, and the angular velocity of the two semicircular trajectories was 0.27 rad/s and 0.27 rad/s, respectively. Finally, the robot stopped at point ( 6 , 2 ). Unlike trajectory A, trajectory B was a 4 m × 4 m square in which a robot moved around the edge of the square counterclockwise, i.e., it started at ( 1 , 0 ) and passed through ( 5 , 0 ), ( 5 , 4 ), and ( 1 , 4 ). Finally, this robot returned to the starting point.
As shown in Figure 12 and Figure 13, the positioning results are not consistent with the reference trajectory in the case of existing outliers. This is because UWB measurements are affected by the multipath effects and NLOS factor. However, SHFAF has the best positioning performance among several contrast approaches, which can be shown in Figure 14. Clearly, the estimate of SHFAF is closer to the truth than that of ESKF and SHAF. ESKF is greatly affected by outliers and time-varying noise, leading to a large positioning error. Although SHAF has better positioning performance than ESKF, there are still large errors of position estimation (especially in positions with abnormal measurements). Note that the experimental results of integrated positioning approaches (i.e., ESKF, SHAF and SHFAF) are better than the individual approaches including INS and UWB. In addition, the performances of these integrated positioning approaches are consistent with the simulation results in Section 4.1.
The absolute positioning error range of trajectory A and B are shown in Figure 15. The position estimation errors less than 0.2 m of SHFAF account for 88.6% of the total sampling points in Figure 15a, while SHAF and ESKF account for 69.3% and 54.3%, respectively. From Figure 15b, we can see that the positioning results of SHFAF, SHAF and ESKF with the error less than 0.2 m stand at 88.2%, 76.8% and 61.1%, respectively. It can be concluded that the error distribution of SHFAF is more concentrated than ESKF and SHAF in a small area of less than 0.2 m. In addition, the proposed integrated positioning approach can also estimate the attitude of the robot because of using the 6-DOF kinematic model. The estimated attitude angles in this experiment are shown in Figure 16.
In summary, results and performance analyses demonstrate that SHFAF has good self-adaptive ability and robustness for applying it in the INS/UWB integrated positioning system, which satisfies the requirement on the positioning accuracy for a mobile robot in complicated indoor environments.

5. Conclusions

This paper focuses on the INS/UWB integrated positioning for autonomous mobile robots in complex indoor environments. To deal with the time-varying noise and outliers in UWB measurements caused by the multipath effects and NLOS factor, a robust INS/UWB integrated positioning approach is proposed based on the Sage–Husa fuzzy adaptive filter. The difficult but important problem of time-varying noise is considered explicitly in SHFAF. So in this paper, the regulatory factor is obtained by the fuzzy inference system to adjust innovation contribution weight of SHFAF adaptively, which facilitates the accurate estimation of measurement noise covariance. Specifically for outliers detection and correction, we propose an effective sliding window estimation with fading coefficients based on innovation orthogonality. As a result, the capture ability to outliers is enhanced and the accuracy of covariance estimation is guaranteed. The effectiveness of what we proposed is demonstrated through simulation and experimental results, which achieves excellent positioning performance with strong robustness for mobile robots in complicated indoor environments. Because of the limited experimental conditions, we can only evaluate the proposed INS/UWB integrated positioning approach in a regular sized laboratory. Actually, a lager experimental site should also be taken into consideration. In this context, more work remains to be done.

Author Contributions

J.L. conceived of the idea and developed the proposed approaches. L.S. gave advice on the research and helped edit the paper. J.P. and Z.H. improved the quality of the manuscript and the completed revision.

Funding

This research received no external funding.

Acknowledgments

This work is supported in part by a grant for the National Natural Science Foundation of China NSFC (Nos. U1504619, 61671139, 61473115), the National Defense Basic Scientific Research Program of China (No. JCKY2018419C001), the Scientific Technology Program of Henan Province, China (No. 182102110397), and the Science and Technology Programme Foundation for the Innovative Team of Henan Province University Grant (No. 18IRTSTHN011).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mostafa, M.; Zahran, S.; Moussa, A.; El-Sheimy, N.; Sesay, A. Radar and visual odometry integrated system aided navigation for UAVS in GNSS denied environment. Sensors 2018, 18, 2776. [Google Scholar] [CrossRef] [PubMed]
  2. Nascimento, P.P.L.L.; Kimura, B.Y.L.; Guidoni, D.L.; Villas, L.A. An Integrated Dead Reckoning with Cooperative Positioning Solution to Assist GPS NLOS Using Vehicular Communications. Sensors 2018, 18, 2895. [Google Scholar] [CrossRef] [PubMed]
  3. Song, Y.; Nuske, S.; Scherer, S. A multi-sensor fusion MAV state estimation from long-range stereo, IMU, GPS and barometric sensors. Sensors 2017, 17, 11. [Google Scholar] [CrossRef] [PubMed]
  4. Ziebold, R.; Medina, D.; Romanovas, M.; Lass, C.; Gewies, S. Performance characterization of GNSS/IMU/DVL integration under real maritime jamming conditions. Sensors 2018, 18, 2954. [Google Scholar] [CrossRef] [PubMed]
  5. Li, L.; Yang, M.; Wang, C.; Wang, B. Hybrid Filtering Framework Based Robust Localization for Industrial Vehicles. IEEE Trans. Ind. Inform. 2018, 14, 941–950. [Google Scholar] [CrossRef]
  6. Huang, Z.; Zhu, J.; Yang, L.; Xue, B.; Wu, J.; Zhao, Z. Accurate 3-D Position and Orientation Method for Indoor Mobile Robot Navigation Based on Photoelectric Scanning. IEEE Trans. Instrum. Meas. 2015, 64, 2518–2529. [Google Scholar] [CrossRef]
  7. Zhao, S.; Huang, B.; Liu, F. Localization of Indoor Mobile Robot Using Minimum Variance Unbiased FIR Filter. IEEE Trans. Autom. Sci. Eng. 2018, 15, 410–419. [Google Scholar] [CrossRef]
  8. Li, Z.; Tian, Z.; Zhou, M.; Zhang, Z.; Jin, Y. Awareness of Line-of-Sight Propagation for Indoor Localization Using Hopkins Statistic. IEEE Sens. J. 2018, 18, 3864–3874. [Google Scholar] [CrossRef]
  9. Sabet, M.T.; Daniali, H.R.; Fathi, A.R.; Alizadeh, E. Experimental analysis of a low-cost dead reckoning navigation system for a land vehicle using a robust AHRS. Robot. Auton. Syst. 2017, 95, 37–51. [Google Scholar] [CrossRef]
  10. Huang, Y.; Zhang, Y.; Wang, X. Kalman-Filtering-Based In-Motion Coarse Alignment for Odometer-Aided SINS. IEEE Trans. Instrum. Meas. 2017, 66, 3364–3377. [Google Scholar] [CrossRef]
  11. Wang, C.; Li, K.; Liang, G.; Chen, H.; Huang, S.; Wu, X. A heterogeneous sensing system-based method for unmanned aerial vehicle indoor positioning. Sensors 2017, 17, 1842. [Google Scholar] [CrossRef] [PubMed]
  12. Hellmers, H.; Kasmi, Z.; Norrdine, A.; Eichhorn, A. Accurate 3D positioning for a mobile platform in non-line-of-sight scenarios based on IMU/magnetometer sensor fusion. Sensors 2018, 18, 126. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, Y.; Ahn, C.K.; Shmaliy, Y.S.; Chen, X.; Li, Y. Adaptive robust INS/UWB-integrated human tracking using UFIR filter bank. Measurement 2018, 123, 1–7. [Google Scholar] [CrossRef]
  14. Fan, Q.; Sun, B.; Sun, Y.; Zhuang, X. Performance Enhancement of MEMS-Based INS/UWB Integration for Indoor Navigation Applications. IEEE Sens. J. 2017, 17, 3116–3130. [Google Scholar] [CrossRef]
  15. Aditya, S.; Molisch, A.F.; Behairy, H.M. A Survey on the Impact of Multipath on Wideband Time-of-Arrival Based Localization. Proc. IEEE 2018, 106, 1183–1203. [Google Scholar] [CrossRef]
  16. Davari, N.; Gholami, A.; Shabani, M. Multirate Adaptive Kalman Filter for Marine Integrated Navigation System. J. Navig. 2017, 70, 628–647. [Google Scholar] [CrossRef]
  17. Atia, M.M.; Liu, S.; Nematallah, H.; Karamat, T.B.; Noureldin, A. Integrated indoor navigation system for ground vehicles with automatic 3-D alignment and position initialization. IEEE Trans. Veh. Technol. 2015, 64, 1279–1292. [Google Scholar] [CrossRef]
  18. Deng, Z.A.; Wang, G.; Qin, D.; Na, Z.; Cui, Y.; Chen, J. Continuous indoor positioning fusing WiFi, smartphone sensors and landmarks. Sensors 2016, 16, 1427. [Google Scholar] [CrossRef] [PubMed]
  19. Jayasiri, A.; Nandan, A.; Imtiaz, S.; Spencer, D.; Islam, S.; Ahmed, S. Dynamic Positioning of Vessels Using a UKF-Based Observer and an NMPC-Based Controller. IEEE Trans. Autom. Sci. Eng. 2017, 14, 1778–1785. [Google Scholar] [CrossRef]
  20. Liu, J.; Cai, B.; Wang, J. Cooperative Localization of Connected Vehicles: Integrating GNSS With DSRC Using a Robust Cubature Kalman Filter. IEEE Trans. Intell. Transp. 2016, 18, 1–15. [Google Scholar] [CrossRef]
  21. Cui, B.; Chen, X.; Tang, X. Improved Cubature Kalman Filter for GNSS/INS Based on Transformation of Posterior Sigma-Points Error. IEEE Trans. Signal Process. 2017, 65, 2975–2987. [Google Scholar] [CrossRef]
  22. Kottath, R.; Poddar, S.; Das, A.; Kumar, V. Window based Multiple Model Adaptive Estimation for Navigational Framework. Aerosp. Sci. Technol. 2016, 50, 88–95. [Google Scholar] [CrossRef]
  23. Liu, Y.; Fan, X.; Lv, C.; Wu, J.; Li, L.; Ding, D. An innovative information fusion method with adaptive Kalman filter for integrated INS/GPS navigation of autonomous vehicles. Mech. Syst. Signal Process. 2018, 100, 605–616. [Google Scholar] [CrossRef]
  24. Al-Sharman, M.K.; Emran, B.J.; Jaradat, M.A.; Najjaran, H.; Al-Husari, R.; Zweiri, Y. Precision landing using an adaptive fuzzy multi-sensor data fusion architecture. Appl. Soft Comput. 2018, 69, 149–164. [Google Scholar] [CrossRef]
  25. Fang, X.; Nan, L.; Jiang, Z.; Chen, L. Noise-aware fingerprint localization algorithm for wireless sensor network based on adaptive fingerprint Kalman filter. Comput. Netw. 2017, 124, 97–107. [Google Scholar] [CrossRef]
  26. Zhong, M.; Guo, J.; Zhou, D. Adaptive In-Flight Alignment of INS/GPS Systems for Aerial Mapping. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 1184–1196. [Google Scholar] [CrossRef]
  27. Han, H.; Xu, T.; Wang, J. Tightly coupled integration of GPS ambiguity fixed precise point positioning and MEMS-INS through a troposphere-constrained adaptive kalman filter. Sensors 2016, 16, 1057. [Google Scholar] [CrossRef] [PubMed]
  28. Nyqvist, H.E.; Skoglund, M.A.; Hendeby, G.; Gustafsson, F. Pose estimation using monocular vision and inertial sensors aided with ultra wide band. In Proceedings of the 2015 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Banff, AB, Canada, 13–16 October 2015. [Google Scholar]
  29. Li, Z.; Chang, G.; Gao, J.; Wang, J.; Hernandez, A. GPS/UWB/MEMS-IMU tightly coupled navigation with improved robust Kalman filter. Adv. Space Res. 2016, 58, 2424–2434. [Google Scholar] [CrossRef]
  30. Wang, W.; Zhang, B.; Wang, D.; Jiang, Y.; Qin, S.; Xue, L. Anomaly detection based on probability density function with Kullback–Leibler divergence. Signal Process. 2016, 126, 12–17. [Google Scholar] [CrossRef]
  31. Wang, L.; Li, S. Enhanced Multi-sensor Data Fusion Methodology based on Multiple Model Estimation for Integrated Navigation System. Int. J. Control Autom. 2018, 16, 295–305. [Google Scholar] [CrossRef]
  32. Zhen, W.; Zeng, S.; Soberer, S. Robust localization and localizability estimation with a rotating laser scanner. In Proceedings of the 2017 IEEE International Conference on Robotics and Automation (ICRA), Singapore, 29 May–3 June 2017; pp. 6240–6245. [Google Scholar]
  33. Herda, L.; Urtasun, R.; Fua, P.; Hanson, A. Automatic Determination of Shoulder Joint Limits Using Quaternion Field Boundaries. Int. J. Robot. Res. 2003, 22, 419–436. [Google Scholar] [CrossRef]
  34. Fang, W.; Zheng, L.; Deng, H. A motion tracking method by combining the IMU and camera in mobile devices. In Proceedings of the 2016 International Conference on Sensing Technology (ICST), Nanjing, China, 11–13 November 2016. [Google Scholar]
  35. Widy, A.; Woo, K.T. Robust attitude estimation method for underwater vehicles with external and internal magnetic noise rejection using Adaptive Indirect Kalman Filter. In Proceedings of the 2017 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Vancouver, BC, Canada, 24–28 September 2017; pp. 2595–2600. [Google Scholar]
  36. Sun, J.; Xu, X.; Liu, Y.; Zhang, T.; Li, Y. FOG random drift signal denoising based on the improved AR model and modified Sage–Husa adaptive Kalman filter. Sensors 2016, 16, 1073. [Google Scholar] [CrossRef] [PubMed]
  37. Ke, W.; Wu, L. Mobile location with NLOS identification and mitigation based on modified Kalman filtering. Sensors 2011, 11, 1641–1656. [Google Scholar] [CrossRef] [PubMed]
  38. Fan, Q.; Sun, B.; Sun, Y.; Wu, Y.; Zhuang, X. Data Fusion for Indoor Mobile Robot Positioning Based on Tightly Coupled INS/UWB. J. Navig. 2017, 70, 1079–1097. [Google Scholar] [CrossRef]
Figure 1. The transformation between the navigation frame and the body frame.
Figure 1. The transformation between the navigation frame and the body frame.
Sensors 19 00950 g001
Figure 2. The framework of robust inertial navigation systems (INS)/ultra wideband (UWB) integrated positioning system.
Figure 2. The framework of robust inertial navigation systems (INS)/ultra wideband (UWB) integrated positioning system.
Sensors 19 00950 g002
Figure 3. The changes of innovation contribution weight (ICW) and adaptive modified innovation contribution weight (AMICW).
Figure 3. The changes of innovation contribution weight (ICW) and adaptive modified innovation contribution weight (AMICW).
Sensors 19 00950 g003
Figure 4. The membership function of input and output. (a) The input of fuzzy inference system. (b) The output of fuzzy inference system.
Figure 4. The membership function of input and output. (a) The input of fuzzy inference system. (b) The output of fuzzy inference system.
Sensors 19 00950 g004
Figure 5. The flowchart of Sage–Husa fuzzy adaptive error state Kalman filter.
Figure 5. The flowchart of Sage–Husa fuzzy adaptive error state Kalman filter.
Sensors 19 00950 g005
Figure 6. The simulation experiment in Gazebo. (a) The initial pose, desired trajectory of mobile robot and distribution of anchors. (b) The navigation and body coordinates in simulations.
Figure 6. The simulation experiment in Gazebo. (a) The initial pose, desired trajectory of mobile robot and distribution of anchors. (b) The navigation and body coordinates in simulations.
Sensors 19 00950 g006
Figure 7. The simulation results in Scenario 1. (a) The estimated position trajectory in which outliers are not considered, but time-varying noise and initial error are set up. (b) The absolute distance error of positioning results.
Figure 7. The simulation results in Scenario 1. (a) The estimated position trajectory in which outliers are not considered, but time-varying noise and initial error are set up. (b) The absolute distance error of positioning results.
Sensors 19 00950 g007
Figure 8. The simulation results in Scenario 2. (a) The positioning results that outliers with large errors exist in measurement data. (b) The absolute distance error of different integrated positioning approaches.
Figure 8. The simulation results in Scenario 2. (a) The positioning results that outliers with large errors exist in measurement data. (b) The absolute distance error of different integrated positioning approaches.
Sensors 19 00950 g008
Figure 9. The hardware communication of INS/UWB integrated system.
Figure 9. The hardware communication of INS/UWB integrated system.
Sensors 19 00950 g009
Figure 10. The INS/UWB integrated positioning experimental platform.
Figure 10. The INS/UWB integrated positioning experimental platform.
Sensors 19 00950 g010
Figure 11. The reference trajectory of integrated positioning experiments.
Figure 11. The reference trajectory of integrated positioning experiments.
Sensors 19 00950 g011
Figure 12. The positioning results of Trajectory A.
Figure 12. The positioning results of Trajectory A.
Sensors 19 00950 g012
Figure 13. The positioning results of Trajectory B.
Figure 13. The positioning results of Trajectory B.
Sensors 19 00950 g013
Figure 14. The absolute positioning errors. (a) Trajectory A. (b) Trajectory B.
Figure 14. The absolute positioning errors. (a) Trajectory A. (b) Trajectory B.
Sensors 19 00950 g014
Figure 15. The distribution of positioning errors. (a) Trajectory A. (b) Trajectory B.
Figure 15. The distribution of positioning errors. (a) Trajectory A. (b) Trajectory B.
Sensors 19 00950 g015
Figure 16. The estimated attitude angles with proposed approach. (a) Trajectory A. (b) Trajectory B.
Figure 16. The estimated attitude angles with proposed approach. (a) Trajectory A. (b) Trajectory B.
Sensors 19 00950 g016
Table 1. The root mean square error (RMSE) with different positioning approaches.
Table 1. The root mean square error (RMSE) with different positioning approaches.
Positioning ApproachRMSE (m)
EastNorthAbsolute
UWB0.56610.57070.8038
ESKF0.28020.27190.3904
SHAF0.23570.21590.3196
SHFAF0.10450.09910.1440
Table 2. The technical specifications of sensors.
Table 2. The technical specifications of sensors.
SensorsSpecification
IMUGyroscopeRange: ±2000 /sec
Resolution: 0.07 /sec/LSB
Gyro rate noise: 0.03 dps/ H z
Data bits: 16 bits
AccelerometerRange: ±16 g
Resolution: 3.9 mg/LSB
Impact resistance: 10,000 g
Data bits: 12 bits
MgnetometerRange: ±8 gauss
Resolution: 5 milli-gauss
Hysteresis: ±25 ppm
Data bits: 16 bits
UWBCommunication rate: 110 kbit/s, 850 kbit/s, 6.8 Mbit/s
Frequency: 3.5 GHz∼6.5 GHz
Transmit power: 35 dbm/MHz∼62 dbm/MHz
Communication distance: 30 m (Indoor), 50 m (Unobstructed)
Ranging error: ±10 cm (Typical), ±30 cm (General hindered)

Share and Cite

MDPI and ACS Style

Liu, J.; Pu, J.; Sun, L.; He, Z. An Approach to Robust INS/UWB Integrated Positioning for Autonomous Indoor Mobile Robots. Sensors 2019, 19, 950. https://doi.org/10.3390/s19040950

AMA Style

Liu J, Pu J, Sun L, He Z. An Approach to Robust INS/UWB Integrated Positioning for Autonomous Indoor Mobile Robots. Sensors. 2019; 19(4):950. https://doi.org/10.3390/s19040950

Chicago/Turabian Style

Liu, Jianfeng, Jiexin Pu, Lifan Sun, and Zishu He. 2019. "An Approach to Robust INS/UWB Integrated Positioning for Autonomous Indoor Mobile Robots" Sensors 19, no. 4: 950. https://doi.org/10.3390/s19040950

APA Style

Liu, J., Pu, J., Sun, L., & He, Z. (2019). An Approach to Robust INS/UWB Integrated Positioning for Autonomous Indoor Mobile Robots. Sensors, 19(4), 950. https://doi.org/10.3390/s19040950

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