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

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., ) with smaller error than positioning approaches with an individual sensor’s data. It is given by , 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 and 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 , where k denotes the time index, , and are the velocity, the offset of acceleration and angular velocity. Consider the following system model: It characterizes the target’s kinematics and sensor measurements, where is the state vector with transition function , is measurement vector obtained by the measurement function . is the input or control vector, and are the process and measurement noise, respectively.

In INS/UWB integrated positioning system, is determined by the kinematics of INS, and 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 and the rotation matrix jointly.

Figure 1

The transformation between the navigation frame and the body frame.

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]: where is the skew symmetric operator, and 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 . It consists of the nominal state and the error state , which satisfies where . ⊕ 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 . Correspondingly, the transition of the nominal state and error state has the following forms [32]: where is the gravity acceleration vector in {N}, and are the acceleration and angular velocity measurements, respectively. and are the measurement noise of the accelerometer and gyroscope, respectively. and 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 , where 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 , then the measurement distance from the UWB tag to the n-th anchor is written as

Similarly, if the true position of tag is , the true distance can be defined as

The positioning problem of the tag can be formulated as finding the solution which minimizes the the error between the measurement distance and true distance, i.e., Note that N may be greater than the dimension of the unknown position vector. From Equation (6), we have

Let the second row, third row, …, and N-th row subtract the first row, respectively. A linear equation is obtained as where

Accordingly, Equation (10) can be solved by using the least squares method:

Meanwhile, the velocity of UWB tag is also obtained as

Suppose that the position and velocity vector of INS at time k are and respectively, then the measurement equation of the integrated positioning system is where is the measurement noise vector, is its covariance matrix, and 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, , , 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.

Figure 2

The framework of robust inertial navigation systems (INS)/ultra wideband (UWB) integrated positioning system.

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 where denotes the nonlinear state transition function. It can be derived by the Euler integrals of Equation (4). As the control inputs, 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

Note that since the mean of the error state is initialized to zero, the linear Equation (17) always returns zero. is the error state transition matrix, it can be obtained by the discretization of Equation (5): where is the noise driven matrix and 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 is defined as and its theoretical covariance is

The theoretical covariance equals the truth when the filter is strictly converged. Correspondingly, the convergence criterion is derived as and the theoretical estimate of MNC matrix is

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

In this equation, is the modified innovation contribution weight (MICW) where . reduces to the traditional ICW when . 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 is a regulatory factor of MICW which can be obtained by the fuzzy inference system.

The estimated MNC at time k (i.e., .

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 was increased, and (b) the regulatory factor (i.e., ) of MICW was obtained through fuzzy inference for improving the estimation ability to .

We define the adaptive modified innovation contribution weight (AMICW) as the product of regulatory factor and MICW (i.e., ), 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 .

Figure 3

The changes of innovation contribution weight (ICW) and adaptive modified innovation contribution weight (AMICW).

To obtain the regulatory factor , we construct a fuzzy reference system. As the input of a fuzzy reference system, the difference coefficient between the estimated innovation covariance and the theoretical covariance should be calculated first, i.e., where is calculated through Equation (23), and is the trace of a matrix. Then, the difference coefficient serves as the input of the fuzzy reference system, and the corresponding output can be obtained: where is regarded as a nonlinear function. Because of the unknown a priori information of measurement noise, an initial estimate of is obtained by the simplified Sage–Husa adaptive filter at the beginning (i.e., ). When , MICW tends to be stable. At this stage, the regulatory factor was inferred by the fuzzy inference system and plays a dominant role in adaptive tuning.

The first adaption for MICW relies on the forgetting factor b, and it is adjusted adaptively in the second time through the regulatory factor .

The inference rules of the above fuzzy logic system are written as

According to Equation (29), it can be seen that intuitively reflects the error between the estimated and theoretical value of innovation covariance, which characterizes the change degree of MNC. The larger the 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 . Conversely, the smaller the , the smaller the contribution of innovation to the estimated MNC. In this case, the present should be estimated by more information of the previous (i.e., slightly decreasing ). Thus should be increased to enhance the innovation contribution for the estimate of . The input and output membership function of fuzzy inference system are shown in Figure 4.

Figure 4

The membership function of input and output. (a) The input of fuzzy inference system. (b) The output of fuzzy inference system.

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 and the theoretical expectation of is

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 . Suppose that the i-th diagonal element of and are and , respectively. If , then where is the predetermined sensitivity threshold of outliers recognition. Thus the outlier is corrected as where . From Equation (34) we can see that the larger (i.e., the deviation between measurement and prediction) is, the smaller .

To approximate the true in Equation (32), the estimated innovation covariance and the expectation need to be calculated. The value 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.

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 where 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 of the error state can be calculated, i.e.,

The estimated error state vector and its corresponding covariance are updated as

Thus the error state is injected into the nominal state, and the true state is estimated by Equation (3). Particularly, Afterwards, the estimator (i.e., ) 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.

Figure 5

The flowchart of Sage–Husa fuzzy adaptive error state Kalman filter.

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 (), (), (), () (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 . The interference noise is random with uniform distributed (i.e., ). 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 .

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.

In simulations, the initial position of this robot is located at () 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 (), 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.

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.

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.

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.

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.

Table 1

The root mean square error (RMSE) with different positioning approaches.

Positioning Approach	RMSE (m)
	East	North	Absolute
UWB	0.5661	0.5707	0.8038
ESKF	0.2802	0.2719	0.3904
SHAF	0.2357	0.2159	0.3196
SHFAF	0.1045	0.0991	0.1440

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.

Figure 9

The hardware communication of INS/UWB integrated system.

Table 2

The technical specifications of sensors.

Sensors		Specification
IMU	Gyroscope	Range: ±2000^∘/secResolution: 0.07^∘/sec/LSBGyro rate noise: 0.03 dps/HzData bits: 16 bits
	Accelerometer	Range: ±16 gResolution: 3.9 mg/LSBImpact resistance: 10,000 gData bits: 12 bits
	Mgnetometer	Range: ±8 gaussResolution: 5 milli-gaussHysteresis: ±25 ppmData bits: 16 bits
UWB		Communication rate: 110 kbit/s, 850 kbit/s, 6.8 Mbit/sFrequency: 3.5 GHz∼6.5 GHzTransmit power: −35 dbm/MHz∼62 dbm/MHzCommunication distance: 30 m (Indoor), 50 m (Unobstructed)Ranging error: ±10 cm (Typical), ±30 cm (General hindered)

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 (), (), () and (), 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.

Figure 10

The INS/UWB integrated positioning experimental platform.

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 () 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 rad/s and 0.27 rad/s, respectively. Finally, the robot stopped at point (). 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 () and passed through (), (), and (). Finally, this robot returned to the starting point.

Figure 11

The reference trajectory of integrated positioning experiments.

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.

Figure 12

The positioning results of Trajectory A.

Figure 13

The positioning results of Trajectory B.

Figure 14

The absolute positioning errors. (a) Trajectory A. (b) Trajectory B.

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.

Figure 15

The distribution of positioning errors. (a) Trajectory A. (b) Trajectory B.

Figure 16

The estimated attitude angles with proposed approach. (a) Trajectory A. (b) Trajectory B.

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.