Next Article in Journal
Influence of Dispersed Oil on the Remote Sensing Reflectance—Field Experiment in the Baltic Sea
Previous Article in Journal
The Simulation of In-Situ Groundwater Detector Response as a Means of Identifying Beta Emitting Radionuclides by Linear Regression Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Improvement of Indoor Positioning Accuracy for UWB Sensors

Department of Electrical Engineering, National Taipei University of Technology, Taipei 10608, Taiwan
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(17), 5731; https://doi.org/10.3390/s21175731
Submission received: 30 May 2021 / Revised: 21 July 2021 / Accepted: 20 August 2021 / Published: 25 August 2021
(This article belongs to the Section Intelligent Sensors)

Abstract

:
Ultra-wideband (UWB) sensors have been widely applied to indoor positioning. The indoor positioning of UWB sensors usually refers to the positioning of the mobile node that interacts with the anchors through radio for calculating the distance between the mobile node and each of the surrounding anchors. The positioning accuracy of the mobile node is affected by the installation positions of surrounding anchors. A mathematical model was proposed in this paper to respectively analyze the mobile node’s 2-dimensional (2D) and 3-dimensional (3D) positioning errors. The factors influencing the mobile node’s positioning errors were explored through the mathematical models. The best installation positions of surrounding anchors were obtained based on the mathematical models. The mobile node’s 2D and 3D positioning errors were reduced based on the anchor positions derived from the mathematical model. Both computer simulations and practical experiments were implemented to justify the results obtained in the mathematical models.

1. Introduction

The Global Positioning System (GPS) is widely used for outdoor positioning and navigation. The GPS receivers receive the GPS signals transmitted from the satellite system for positioning. The received GPS signals are utilized to calculate the user’s 3-dimensional (3D) position. However, the GPS signals transmitted from the satellites are easily blocked by building structures, including roofs, floors, walls, etc. Indoor positioning is not feasible using GPS signals. With the increasing need for indoor positioning for customers in shopping malls, patients in hospitals, containers in terminals, and materials in warehouses [1,2,3,4], the indoor positioning technology industry is growing rapidly. Various indoor positioning systems are already in the market, including ultrasound [5], infrared [6], Bluetooth [7], Zigbee [8], RFID [9], and Wi-Fi [10], and each has its advantages and disadvantages. Ultrasound is susceptible to multipath effects, and its detection range is limited by distance, angle, and obstruction. Infrared is vulnerable to light and obstacles. Bluetooth, Zigbee, and Wi-Fi are based on a received signal strength indicator (RSSI) [8] and are severely affected by barriers, resulting in extensive positioning errors. In addition, RFID is a short-range wireless communication technology; therefore, it is limited by distance. The ultra-wideband (UWB) sensor proposed herein provides reliable, high-precision ranging, covers a relatively wide area, and can bypass certain obstacles [11]. The proposed UWB positioning system comprises multiple anchors and mobile nodes; the measured distances between the anchors and the mobile node are used for positioning.
Localization methods for UWB sensors include trilateration [12], the least-squares (LS) method [13], partial filtering [14], and the extended Kalman filter (EKF) [15]. EKF outperforms other methods regarding positioning accuracy and stability, particularly when the UWB-sensor data have outliers. To further improve positioning accuracy, the UWB mobile node is integrated with an inertia measurement unit (IMU) [16,17,18,19,20,21] to eliminate outliers in the received data. UWB mobile node measures the values such as angle of arrival (AOA) [22], time of arrival (TOA) [23], time of flight (TOF), and time differential of arrival (TDOA) [14], etc., to calculate the position of a mobile node. TOA, TOF, and AOA are based on radio transmission and receiving time differences, which can be used to measure the distance and directional angle. TDOA determines the position of the mobile node based on the time difference between the signals received by different anchors.
Numerous methods have been employed to improve UWB’s positioning accuracy. Multi-sensor fusion is effective for achieving high-precision positioning [24,25,26,27,28]. Combining UWB positioning with an IMU can eliminate drift-free output in UWB-sensor data and correct accumulated IMU errors [29,30]. Both Light-of-sight (LOS) and non-light-of-sight (NLOS) analysis methods can improve the stability and accuracy of UWB positioning systems [31,32,33]. The position accuracy of the mobile node is affected by the installation positions of anchors. The anchors are not all placed in a plane to reduce the mobile node’s 3D position accuracy [14,34,35]. Experimental results in [34] indicated that positioning accuracy differs in the z-axis at different anchor heights, but the reasons were not thoroughly explained. The UWB positioning error was investigated in [36] by using the concept of dilution of precision [37]. A 2-dimensional (2D) positioning method was proposed for automated guided vehicles (AGVs) in [38] to investigate the effects of the distance between neighboring anchors on the positioning accuracy. The mathematical model estimating positioning errors proposed in [38] was specifically for the environment, such as the corridor where the anchors are installed on the wall of the corridor. The heights (positions on the z-axis) of all the anchors on the wall were assumed to be the same. No constraints such as in [38] have been imposed on the installation positions of anchors. In other words, the anchors are not limited to be installed on the wall of a corridor, and the anchors are not necessarily to be installed at the same height on the wall. The mathematical model proposed in [38] was only for calculating the 2D positions of the mobile node. In contrast to the model proposed in [36], the model presented in this paper can be utilized to calculate both 2D and 3D positions. There is no other suitable mathematical model for UWB sensors estimating mobile node’s 3D positioning errors to our best knowledge.
The positioning accuracy of UWB sensors refers to the positioning accuracy of the mobile node. The mobile node’s 2D or 3D position is calculated based on the distance between the mobile node and its surrounding anchors. Therefore, UWB positioning accuracy depends on the installation positions of surrounding anchors. It takes at least three and four anchors for 2D and 3D positioning, respectively. Mathematical models are proposed to analyze both the 2D and 3D positioning accuracies of the UWB mobile node for both the LOS and NLOS conditions. The variance of 2D and 3D positioning errors are mathematically derived. Technically speaking, the mathematical model of 2D positioning errors is a special case of 3D. The best arrangement of anchor installation positions reducing the positioning errors is obtained based on the mathematical model of positioning errors. Computer simulations and practical experiments are then designed to verify anchor installation positions obtained from the mathematical model. The root-mean-square positioning error (RMSPE) of the mobile node is calculated in different arrangements of anchor installation positions. It is shown that both the 2D and 3D RMSPE are reduced if the anchors are installed at the positions suggested in the mathematical model of positioning errors.
The main contributions of this paper are as follows:
  • Mathematical models of both 2D and 3D positioning errors for UWB sensors were derived. To the best of our knowledge, this paper is the first one analyzing suitable UWB anchor installation positions based on a mathematical model of 3D positioning errors.
  • The mathematical models of 2D and 3D positioning errors impose no constraints on anchors’ installation positions. The models are general enough to analyze mobile node’s positioning errors corresponding to any anchor installation positions.
  • Anchor installation positions were suggested based on the mathematical model of 2D and 3D positioning errors for both LOS and NLOS conditions so that the RMSPE can be significantly reduced.
  • Both computer simulations and practical experiments were conducted to verify that the anchor installation positions suggested based on the mathematical model of positioning errors can significantly reduce the RMSPE.
The remainder of the paper is organized as follows. The problem statements and indoor positioning system with UWB sensors are described in Section 2. The mathematical models of 2D and 3D positioning errors are derived in Section 3. The factors affecting the mobile node’s position accuracy are also analyzed. The computer simulations using MATLAB are designed in Section 4 to justify that the anchor installation positions suggested based on the mathematical model do significantly reduce the positioning error. Several practical experiments are further conducted in Section 5 to explain the results obtained through computer simulations practically. Finally, the conclusions are drawn in Section 6.

2. Problem Statement

System Design

The indoor positioning system utilized in this paper was implemented using UWB sensors DecaWave DWM1000 [39], with an optional frequency band range from 3.5 to 6.5 GHz. The data transmission rates included 110 kbps, 850 kbps, and 6.8 Mbps. The anchor was implemented as in Figure 1a using an STM32F0 microcontroller unit (MCU) with the main frequency of 48 MHz and a UWB sensor DWM1000. The mobile node was implemented as in Figure 1b, combining an STM32F4 MCU with the main frequency of 168 MHz, an IMU MPU 9250, and a UWB sensor DWM1000. An embedded system made with a Raspberry Pi 3 was utilized in the mobile node for positioning calculation and the implementation of EKF to improve positioning results. The sampling frequency for positioning was set as 100 Hz.
The positioning of the mobile node was based on measuring the distance between the mobile node and every anchor. The distance measurement relied on the TOF of the radio between the transmitter and the receiver. Two-way radio transmission and receiving were conducted between both sensors. The mobile node was designated as an initiator that initiates the two-way radio communication, while the anchor was set as a responder. The initiator transmits a message through radio to the responder and records the timestamp of transmission. As the responder receives the message, it sends the same message back to the initiator after a preset time delay. The initiator gets the message and records the timestamp of receiving. The TOF is calculated based on the timestamps of transmission and receiving, and the preset time delays at both the responder and the initiator. The calculation of TOF is affected by the clock drift, frequency drift, fading, shadowing, multi-path propagation, etc. Although more delicate ranging methods such as the asymmetric double-sided two-way-ranging (TWR) method [39] or the TWR method integrated with neural network model [40] have been proposed, noise in the measurement of TOF is unavoidable. Let T t o f and T ^ t o f be the ideal TOF and the measured TOF, respectively, and δ T t o f be the noise in TOF measurement. Then,
T ^ t o f = T t o f + δ T t o f .
The measured distance r m between the mobile node and the anchor is calculated as
r m = c ( T t o f + δ T t o f ) = r + ε ,
where c is the speed of radio wave, r is the ideal distance between the mobile node and the anchor, and ε is the noise of measured distance due to the noise of TOF, δ T t o f . The noise ε in (2) is usually modeled as an independent and identical distributed (i.i.d.) random variable. The distribution of the noise ε is usually modeled as the Gaussian distribution [41,42] if the mobile node and the anchor are in the LOS condition. The Gaussian distribution is parameterized with the mean μ G and the variance σ G 2 , i.e., ε N μ G , σ G 2 . The noise modeled in the LOS condition is mainly for the environment that no obstacles are placed to block the radio propagation between the mobile node and the corresponding anchor.
The distribution of the noise ε is modeled as the skew-t distribution [43] for the NLOS condition. The skew-t distribution is parameterized by its location parameter μ S T , spread parameter σ S T 2 , shape parameter δ S T , and degree of freedom ν S T . The probability density function (PDF) of the skew-t distribution for the random variable z is defined as
S T ( z ; μ S T , σ S T 2 , δ S T , ν S T ) = 2 t ( z ; μ S T , δ S T 2 + σ S T 2 , ν S T ) B ( z ˜ ; 0 , 1 , ν S T + 1 ) ,
where t ( · ) denotes the PDF of Student’s t-distribution defined as follows:
t ( z ; μ S T , σ S T 2 , ν S T ) = Γ ( ν S T + 1 2 ) σ S T ν S T π Γ ( ν S T 2 ) ( 1 + ( z μ S T ) 2 ν S T σ S T 2 ) ν S T + 1 2 .
Note that Γ ( · ) in (4) denotes the Gamma function. The function B ( · ) in (3) denotes the cumulative distribution function (CDF) of Student’s t-distribution and the random variable z ˜ in B is defined as:
z ˜ = ( z μ S T ) δ S T σ S T ν S T + 1 ν S T ( δ S T 2 + σ S T 2 ) + ( z μ S T ) 2 .
Therefore, the distance measurement noise ε S T ( z ; μ S T , σ S T 2 , δ S T , ν S T ) for the NLOS condition where skew t-distribution S T ( · ) is defined in (3). The noise modeled in the NLOS condition is mainly for the environment that obstacles are placed between the mobile node and the anchor so that the LOS radio propagation is blocked. The radio multi-path propagation due to moving and/or stationary obstacles usually result in the noise modeled in the NLOS condition [44].

3. Mathematical Model of Positioning Errors

3.1. 3D Positioning

At least four anchors are required for 3D positioning of a mobile node, as shown in Figure 2. Denote q i as the vector containing the coordinate ( x i , y i , z i ) of the ith anchors, i.e., q i = [ x i , y i , z i ] T , i = 1, …,4. The 3D position of the mobile node can be determined by first measuring the distance between the mobile node and every anchor. An EKF is utilized to calculate the position of the mobile node. Denote p and p ^ as the vector containing the coordinate the mobile node’s actual 3D position (x, y, z) and measured 3D position ( x m , y m , z m ) , respectively, i.e., p = [ x , y , z ] T and p ^ = [ x m , y m , z m ] T . The measured distance between the mobile node and the ith anchor is calculated as
r i m = ( x m x i ) 2 + ( y m y i ) 2 + ( z m z i ) 2 ,   i = 1 ,   , 4 .
Similarly, the actual distance between the mobile node and the ith anchor is calculated as
r i = ( x x i ) 2 + ( y y i ) 2 + ( z z i ) 2 ,   i = 1 ,   , 4 .
Referring to (2), the measured distance can be represented as
r i m = r i + ε i ,   i = 1 ,   , 4 .
where ε i denotes the measurement noise of the distance between the mobile node and the ith anchor.
Referring to Figure 2, let k = [ x 2 x 1 , y 2 y 1 , z 2 z 1 ] T , l = [ x 3 x 1 , y 3 y 1 , z 3 z 1 ] T , and m = [ x 4 x 1 , y 4 y 1 , z 4 z 1 ] T be the vectors from the position q 1 to q 2 , q 1 to q 3 , and q 1 , to q 4 , respectively. Note that three points at different places form a plane. For the convenience of analysis, let the vectors k and l be on the same plane.
Let n = k × l , d = l × m , and g = k × m , where
n = ( y 2 y 1 ) ( z 3 z 1 ) ( z 2 z 1 ) ( y 3 y 1 ) ( x 2 x 1 ) ( z 3 z 1 ) + ( z 2 z 1 ) ( x 3 x 1 ) ( x 2 x 1 ) ( y 3 y 1 ) ( y 2 y 1 ) ( x 3 x 1 ) ,
d = ( y 4 y 1 ) ( z 3 z 1 ) ( z 4 z 1 ) ( y 3 y 1 ) ( x 4 x 1 ) ( z 3 z 1 ) + ( z 4 z 1 ) ( x 3 x 1 ) ( x 4 x 1 ) ( y 3 y 1 ) ( y 4 y 1 ) ( x 3 x 1 ) ,
g = ( y 4 y 1 ) ( z 2 z 1 ) ( z 4 z 1 ) ( y 2 y 1 ) ( x 4 x 1 ) ( z 2 z 1 ) + ( z 4 z 1 ) ( x 2 x 1 ) ( x 4 x 1 ) ( y 2 y 1 ) ( y 4 y 1 ) ( x 2 x 1 ) .
With the measured distance r i m defined in (6),
( r 1 m 2 r 2 m 2 ) ( r 1 m 2 r 3 m 2 ) ( r 1 m 2 r 4 m 2 ) = 2 x 2 x 1 x m + 2 y 2 y 1 y m + 2 z 2 z 1 z m + x 1 2 x 2 2 + y 1 2 y 2 2 + z 1 2 z 2 2 2 x 3 x 1 x m + 2 y 3 y 1 y m + 2 z 3 z 1 z m + x 1 2 x 3 2 + y 1 2 y 3 2 + z 1 2 z 3 2 2 x 4 x 1 x m + 2 y 4 y 1 y m + 2 z 4 z 1 z m + x 1 2 x 4 2 + y 1 2 y 4 2 + z 1 2 z 4 2 .
Therefore,
( r 1 m 2 r 2 m 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 m 2 r 3 m 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) ( r 1 m 2 r 4 m 2 ) + ( x 4 2 x 1 2 ) + ( y 4 2 y 1 2 ) + ( z 4 2 z 1 2 ) = 2 x 2 x 1 y 2 y 1 z 2 z 1 x 3 x 1 y 3 y 1 z 3 z 1 x 4 x 1 y 4 y 1 z 4 z 1 x m y m z m .
Referring to (13), let
Δ r m = ( r 1 m 2 r 2 m 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 m 2 r 3 m 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) ( r 1 m 2 r 4 m 2 ) + ( x 4 2 x 1 2 ) + ( y 4 2 y 1 2 ) + ( z 4 2 z 1 2 ) ,
and
H = x 2 x 1 y 2 y 1 z 2 z 1 x 3 x 1 y 3 y 1 z 3 z 1 x 4 x 1 y 4 y 1 z 4 z 1 .
Then, (13) is expressed as
Δ r m = 2 H p ^ ,
where p ^ is the measured 3D position of the mobile node, p ^ = [ x m , y m , z m ] T . If H is nonsingular, i.e., det(H) ≠ 0, the measured position of the mobile node can be calculated as
p ^ = 1 2 H 1 Δ r m .
Since n = k × l , n is perpendicular to the plane formed by three anchor positions q 1 , q 2   and   q 3 . Referring to (14), p ^ is calculated under the condition that det(H) ≠ 0. However, det H = n T m , the condition det(H) ≠ 0 results in the condition that m is not perpendicular to the vector n, i.e., m is not on the plane formed by three anchor positions q 1 , q 2   and   q 3 . In other words, the fourth anchor position q 4 cannot be on the same plane with q 1 , q 2   and   q 3 in order to have a deterministic position p ^ for the mobile node.
Substituting the actual distance r i between the mobile node to the ith anchor into (12)–(14), the actual position p associated with the actual distance can be theoretically determined similar to the measured position in (17) as
p = 1 2 H 1 Δ r ,
where
Δ r = ( r 1 2 r 2 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 2 r 3 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) ( r 1 2 r 4 2 ) + ( x 4 2 x 1 2 ) + ( y 4 2 y 1 2 ) + ( z 4 2 z 1 2 ) .
Denote Δ p as the mobile node’s positioning error,
Δ p p p ^ = 1 2 H 1 Δ r Δ r m = 1 2 H 1 Δ λ .
where Δ λ = Δ r - Δ r m . Referring to (14) and (19), Δ λ is calculated as
Δ λ = 2 r 2 ε 2 2 r 1 ε 1 + ε 2 2 ε 1 2 2 r 3 ε 3 2 r 1 ε 1 + ε 3 2 ε 1 2 2 r 4 ε 4 2 r 1 ε 1 + ε 4 2 ε 1 2 .
Substituting (21) into (20) and going through some derivations yields
Δ p = a 2 M 3 m + a 3 M 2 m + a 4 n n T m ,
where
M 2 = 0 z 2 z 1 ( y 2 y 1 ) ( z 2 z 1 ) 0 x 2 x 1 y 2 y 1 ( x 2 x 1 ) 0 ,
M 3 = 0 ( z 3 z 1 ) y 3 y 1 z 3 z 1 0 ( x 3 x 1 ) ( y 3 y 1 ) x 3 x 1 0 ,
a 2 = 1 2 2 r 2 ε 2 2 r 1 ε 1 + ε 2 2 ε 1 2 ,
a 3 = 1 2 2 r 3 ε 3 2 r 1 ε 1 + ε 3 2 ε 1 2 ,
a 4 = 1 2 2 r 4 ε 4 2 r 1 ε 1 + ε 4 2 ε 1 2 .
Referring to (22),
Δ p T Δ p = a 2 2 m T M 3 T M 3 m ( n T m ) 2 + a 2 a 3 m T M 3 T M 2 m ( n T m ) 2 + a 2 a 4 m T M 3 T n ( n T m ) 2 + a 2 a 3 m T M 2 T M 3 m ( n T m ) 2 + a 3 2 m T M 2 T M 2 m ( n T m ) 2 + a 3 a 4 m T M 2 T n ( n T m ) 2 + a 2 a 4 n T M 3 m ( n T m ) 2 + a 3 a 4 n T M 2 m ( n T m ) 2 + a 4 2 n T n ( n T m ) 2 .
Further derivations in (28) yields
Δ p T Δ p = a 2 2 d 2 ( n T m ) 2 + a 3 2 g 2 ( n T m ) 2 + a 4 2 n 2 ( n T m ) 2 + 2 a 2 a 3 d T g a 2 a 4 d T n + a 3 a 4 g T n ( n T m ) 2 .
Referring to (29), d g d T g , d n d T n and g n g T n . Therefore,
a 2 a 3 d T g a 2 a 4 d T n + a 3 a 4 g T n a 2 a 3 d g + a 2 a 4 d n + a 3 a 4 g n .
Referring to (30), a 2 d a 3 g 2 = a 2 2 d 2 + a 3 3 g 2 2 a 2 a 3 d g 0 implies that ( a 2 2 d 2 + a 3 2 g 2 ) / 2 a 2 a 3 d g . Similarly, ( a 2 2 d 2 + a 4 2 n 2 ) / 2 a 2 a 4 d n and ( a 3 2 g 2 + a 4 2 n 2 ) / 2 a 3 a 4 g n . Therefore, (30) can be further written as
a 2 a 3 d T g a 2 a 4 d T n + a 3 a 4 g T n a 2 2 d 2 + a 3 2 g 2 + a 4 2 n 2 .
Substituting (31) into (29) yields
Δ p T Δ p 3 a 2 2 d 2 ( n T m ) 2 + a 3 2 g 2 ( n T m ) 2 + a 4 2 n 2 ( n T m ) 2 .
Taking the expectation for both sides of (32),
E Δ p T Δ p 3 E a 2 2 d 2 ( n T m ) 2 + E a 3 2 g 2 ( n T m ) 2 + E a 4 2 n 2 ( n T m ) 2 .
Because d = l × m , g = k × m and n = k × l , d 2 = l 2 m 2 sin 2 α , g 2 = k 2 m 2 sin 2 β , and n 2 = k 2 l 2 sin 2 φ , where α , β and φ are the angles between vectors l and m , k and m , k and l , respectively, as shown in Figure 2. The inequality in (33) can be rewritten as
E Δ p T Δ p 3 E a 2 2 l 2 m 2 sin 2 α ( n T m ) 2 + E a 3 2 k 2 m 2 sin 2 β ( n T m ) 2 + E a 4 2 k 2 l 2 sin 2 φ ( n T m ) 2 .
It follows that
E Δ p T Δ p 3 E a 2 2 l 2 m 2 ( n T m ) 2 + E a 3 2 k 2 m 2 ( n T m ) 2 + E a 4 2 k 2 l 2 ( n T m ) 2 .
Because n T m = n m cos ψ and n = k l sin φ , n T m = m k l sin φ cos ψ , where ψ is the angle between m and the normal vector n of the plane containing the three anchors q 1 , q 2 , and q 3 . It follows that
E ( Δ p T Δ p ) 3 sin 2 φ cos 2 ψ E a 2 2 k 2 + E a 3 2 l 2 + E a 4 2 m 2 .
E ( Δ p T Δ p ) in (36) is the variance of the mobile node’s 3D positioning error. The bound for E ( Δ p T Δ p ) in (36) applies for both the LOS and NLOS conditions. Note that E ( a 2 2 ) > 0 , E ( a 3 2 ) > 0 , and E ( a 4 2 ) > 0 in (36). Given that the distance between pairs of anchors, k , l , and m are decided, (36) suggests that arranging the anchor installation positions with φ = 90° and ψ = 0° leads to reasonably small positioning. Moreover, installing the anchors as separated as possible, i.e., k, l, and m being as large as possible, results in reasonably small 3D positioning errors. Note that anchors should be installed at the positions so that k, l, and m are as large as possible, provided that the radio receiving intensity between the mobile node and the anchor is within the rated range.

3.2. 2D Positioning

Essentially, the mathematical model of 2D positioning error is a special case of the one of 3D positioning error. It takes at least three anchors to determine the 2D position of a mobile node. Denote q i as the vector containing the coordinate ( x i , y i , z i ) of the ith anchors, i.e., q i = [ x i , y i , z i ] T , i = 1, …, 3. Referring to Figure 3, let k = [ x 2 x 1 , y 2 y 1 , z 2 z 1 ] T and l = [ x 3 x 1 , y 3 y 1 , z 3 z 1 ] T be the vectors from the position q 1 to q 2 and q 1 to q 3 , respectively. Compared with Figure 2, the fourth anchor q 4 in Figure 2 is not needed for 2D positioning as shown in Figure 3. Denote p′ and p ^ as the vector containing the 2D coordinate the mobile node’s actual position (x, y) and measured position ( x m , y m ) , respectively, i.e., p ^ = [ x , y ] T and p ^ = [ x m , y m ] T . Most of the mathematical analysis for 3D positioning in the previous subsection can be used for the analysis of 2D positioning, except that z = 0 and z m = 0 .
Setting z m = 0 in (12) yields
( r 1 m 2 r 2 m 2 ) ( r 1 m 2 r 3 m 2 ) = 2 x 2 x 1 x m + 2 y 2 y 1 y m + x 1 2 x 2 2 + y 1 2 y 2 2 + z 1 2 z 2 2 2 x 3 x 1 x m + 2 y 3 y 1 y m + x 1 2 x 3 2 + y 1 2 y 3 2 + z 1 2 z 3 2 .
Therefore,
( r 1 m 2 r 2 m 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 m 2 r 3 m 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) = 2 x 2 x 1 y 2 y 1 x 3 x 1 y 3 y 1 x m y m .
Referring to (38), let
Δ r m = ( r 1 m 2 r 2 m 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 m 2 r 3 m 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) ,
and
H = x 2 x 1 y 2 y 1 x 3 x 1 y 3 y 1 .
Then, (38) is expressed as
Δ r m = 2 H p ^ ,
If H′ is nonsingular, i.e., det(H′) ≠ 0, the measured position of the mobile node is calculated as
p ^ = 1 2 H 1 Δ r m .
Referring to (42), p ^ has a deterministic solution under the condition that det(H′) ≠ 0. Referring to (40), let k = [ x 2 x 1 , y 2 y 1 ] T and l = [ x 3 x 1 , y 3 y 1 ] T . If φ is the angle between vectors k and l , det H = k × l = k l sin φ . Note that k′ and l′ are the projection of k and l onto the X-Y plane, respectively, as shown in Figure 4a,b. Therefore, p ^ has a deterministic solution provided that φ 0 , i.e., k′ is not parallel with l′.
Referring to (18), the actual position p′ of the mobile node can be theoretically determined similar to the measured position in (42) as
p = 1 2 H 1 Δ r ,
where
Δ r = ( r 1 2 r 2 2 ) + ( x 2 2 x 1 2 ) + ( y 2 2 y 1 2 ) + ( z 2 2 z 1 2 ) ( r 1 2 r 3 2 ) + ( x 3 2 x 1 2 ) + ( y 3 2 y 1 2 ) + ( z 3 2 z 1 2 ) .
Denote Δ p as the mobile node’s positioning error,
Δ p p p ^ = 1 2 H 1 Δ r Δ r m = 1 2 H 1 Δ λ .
where Δ λ = Δ r - Δ r m and
Δ λ = 2 r 2 ε 2 2 r 1 ε 1 + ε 2 2 ε 1 2 2 r 3 ε 3 2 r 1 ε 1 + ε 3 2 ε 1 2 .
Substituting (46) into (45) and performing some derivations yields
Δ p = a 2 M l a 3 M k k l sin φ ,
where M = 1 0 0 1 , a 2 and a 3 are as shown in (25) and (26), respectively. Then,
Δ p T Δ p = 1 k 2 l 2 sin 2 φ a 2 2 l 2 + a 3 2 k 2 2 a 2 a 3 l T k .
Referring to (48), l k l T k . Therefore,
2 a 2 a 3 l T k 2 a 2 a 3 l k .
Since a 2 l a 3 k 2 = a 2 2 l 2 + a 3 2 k 2 2 a 2 a 3 l k 0 , it implies that a 2 2 l 2 + a 3 2 k 2 2 a 2 a 3 l k . The inequality (49) can be further written as
2 a 2 a 3 l T k a 2 2 l 2 + a 3 2 k 2 .
Substituting (50) into (48) yields
Δ p T Δ p 2 k 2 l 2 sin 2 φ a 2 2 l 2 + a 3 2 k 2 .
Taking the expectation for both sides of (51),
E Δ p T Δ p 2 k 2 l 2 sin 2 φ E a 2 2 l 2 + E a 3 2 k 2 .
It can be further simplified as
E Δ p T Δ p 2 sin 2 φ E ( a 2 2 ) k 2 + E ( a 3 2 ) l 2 .
Similar to the 3D position error, the variance of the mobile node’s 2D positioning error E ( Δ p T Δ p ) in (53) also applies for both the LOS and NLOS conditions. Given that the distance between pairs of anchors k and l are decided, (53) suggests that arranging the anchor positions with φ = 90° leads to reasonably small positioning errors. Moreover, installing the anchors as separated as possible, i.e., k′ and l′ being as large as possible, results in reasonably small 2D positioning errors.

4. Computer Simulations

Computer simulations were implemented to simulate the mobile node’s 2D and 3D positioning. The mobile node was designed to move following a 2D and 3D straight line in a 2D and 3D environment, respectively. A total of 8000 positioning samples were calculated in 2D or 3D simulations, respectively. In other words, the sampling frequency was designed as 100 Hz and 80 s of positioning errors were calculated and recorded for 2D and 3D simulations. The RMSPE was calculated based on the N samples calculated and recorded at the mobile node. Denote Δ p ¯ as the RMSPE and Δ p i as the ith sample of positioning error, i = 1, …, N. Then,
Δ p ¯ = i = 1 N Δ p i T Δ p i N ,
where N = 8000. The mobile node’s RMSPE Δ p ¯ defined in (54) is utilized for the simulation.
Both the LOS and NLOS conditions are to be simulated. The distance measurement noise ε N μ G , σ G 2 for the LOS condition and ε S T ( z ; μ S T , σ S T 2 , δ S T , ν S T ) for the NLOS condition. The parameters of the Gaussian distribution for the LOS condition are set as the mean μ G = 0 and the standard deviation σ G = 0.1 m . The parameters of the skew t-distribution for the NLOS condition are set as μ S T = 0.1 m , σ S T = 0.3 m , δ S T = 3 , and ν = 4 . The variance of 3D position errors defined in (36) can be further simplified if the distance measurement noise ε N 0 , σ G 2 . The terms E ( a 2 2 ) , E ( a 3 2 ) , and E ( a 4 2 ) in (36) can be replaced with E a 2 2 = σ G 2 ( r 1 2 + r 2 2 ) , E a 3 2 = σ G 2 ( r 1 2 + r 3 2 ) , and E a 4 2 = σ G 2 ( r 1 2 + r 4 2 ) for the LOS condition. Similarly, the terms E ( a 2 2 ) and E ( a 3 2 ) in the variance of 2D position errors defined in (53) can also be replaced with E a 2 2 = σ G 2 ( r 1 2 + r 2 2 ) and E a 3 2 = σ G 2 ( r 1 2 + r 3 2 ) for the LOS condition.

4.1. 2D positioning Simulation

In order to simulate the variation of Δ p ¯ with respect to the angle φ between the vectors k’ and l’, the projection of installation positions of three anchors q 1 ,   q 2   and   q 3 onto the X-Y plane are shown in Figure 4a,b, respectively, for the convenience of illustration. Both Figure 4a,b are essentially similar except that the sizes of the simulation environment are different. As shown in Figure 4a,b, 2 anchors q 1 and q 2 were installed at fixed positions (0m, 0m) and (0m, 5m), and (0m, 0m) and (0m, 20m), respectively. The mobile node moves from (0.5m, 0.5m) to (4.5m, 4.5m) in Figure 4a and from (2m, 2m) to (18m, 18m) in Figure 4b for four round trips. The anchor q 3 was simulated to be installed at different positions on the circumference with a radius of 5m and 20m, respectively, in Figure 4a,b, but both centered at (0m, 0m). The angle φ ranged from 0 ° to 160 ° as in Figure 4a,b. The heights of q 1 ,   q 2   and   q 3 on the z-axis were assumed to be 4m, 4.2m, and 4.1m, respectively.
The RMSPE Δ p ¯ defined in (54) with different angles of φ for the LOS and NLOS conditions are listed in Table 1 and Table 2, respectively. It is shown in Table 1 that Δ p ¯ was minimum if φ = 90° in both 5m × 5m and of 20m × 20m environments for the LOS condition. Table 2 shows that Δ p ¯ was minimum if φ = 70° in both 5m × 5m and of 20m × 20m environments for the NLOS condition. However, Δ p ¯ with φ = 90° in both 5m × 5m and of 20m × 20m environments were very close to the minimum values occurring at φ = 70°. The variations of Δ p ¯ with the angle φ for the LOS and NLOS conditions in both Table 1 and Table 2 are illustrated in Figure 5a,b, respectively. The suggestion based on (53) that the anchor installation arrangement with φ = 90° leads to reasonably small positioning error for both the LOS and NLOS conditions is thus justified from both Table 1 and Table 2.
The 2D RMSPE Δ p ¯ is also affected by the anchor distances k   and   l . The variation of Δ p ¯ with different k   and   l for the LOS and NLOS conditions are calculated in Table 3 and Table 4, respectively. Both Table 3 and Table 4 show that the 2D RMSPE Δ p ¯ decreased as k and/or l increased for both the LOS and NLOS conditions. The variations of Δ p ¯ with anchor distances k   and   l for the LOS and NLOS conditions are also illustrated in Figure 6a,b, respectively. The suggestion based on (53) that installing the anchors as separated as possible results in reasonably small positioning error is thus justified.
Since the installation arrangement of anchors suggested in this paper was based on the upper bound of the variance of positioning errors in (53), the accuracy of the upper bound is worth evaluation. Denote the upper bound of variance in (53) for the ith sample as V i u , then
V i u = 2 sin 2 φ E ( ( a 2 , i ) 2 ) k 2 + E ( ( a 3 , i ) 2 ) l 2
where a 2 , i and a 3 , i are the parameters a 2 and a 3 corresponding to the ith sample data. Denote Δ p ¯ u as the upper bound of average positioning error based on (53) for all N samples, then
Δ p ¯ u = i = 1 N V i u N .
The upper bound of average positioning error Δ p ¯ u in (56) was compared with the one proposed in [38]. The theoretical Cramer-Rao lower bound (CRLB) was also calculated and plotted for comparison. Let the upper bound of positioning error for the ith data Δ p i u = V i u . The variations of Δ p i u with samples estimated by our approach and by the approach in [38] are compared in Figure 7a and Figure 8a for the LOS and NLOS condition, respectively. Recall that 8000 samples were recorded in the simulation. Only 20 percent of recorded data are plotted in Figure 7a and Figure 8a for convenience of illustration. Both Figure 7a and Figure 8a show that the upper bound of positioning error Δ p i u was much lower than the ones calculated by the approach in [38], and much closer to the CRLB for the LOS and NLOS conditions, respectively. The variations of upper bound of average positioning error Δ p ¯ u with k for the situation that l = 20m and φ = 90 ° were also compared as in Figure 7b and Figure 8b for the LOS and NLOS conditions, respectively. Figure 7c and Figure 8c show the comparison of Δ p ¯ u with φ if k = l = 20m for the LOS and NLOS conditions, respectively. It is shown in Figure 7b,c and Figure 8b,c that Δ p ¯ u was much lower than the one calculated by the approach in [38], and much closer to the CRLB for both the LOS and NLOS conditions.

4.2. 3D Positioning Simulation

The mobile node’s positioning error defined in (22) was utilized for the 3D simulation. In order to simulate the variations of Δ p ¯ with respect to the angles φ between the vectors k and l and the angle ψ between the vectors m and n, the anchors q 1 ,   q 2 , q 3 , and q 4 were assumed to be installed in a smaller environment 5m × 5m × 5m and a larger environment 20m × 20m × 20m, as shown in Figure 9a,b, respectively. The mobile node moves from (0.5m, 0.5m, 0.5m) to (4.5m, 4.5 m, 4.5m) in Figure 9a and from (2m, 2m, 2m) to (18m, 18m, 18m) in Figure 9b for four round trips. The anchor q 3 and q 4 were both simulated to be installed at different positions on the circumference with a radius of 5m and 20m, respectively, in Figure 9a,b, but both centered at the position of q 1 providing that other anchors were fixed at the positions shown in Figure 9a,b. The positioning errors Δ p ¯ were calculated with angles φ and ψ varying from 0 ° to 160 ° . Note that φ varies by changing the angle between the vectors k and l while ψ varies by changing the angle between the vectors m and n.
The RMSPE Δ p ¯ defined in (54) with different angles of φ for the LOS and NLOS conditions are listed in Table 5 and Table 6, respectively. It is shown in Table 5 that Δ p ¯ was minimum at φ = 80° in the 5m × 5m × 5m environment while Δ p ¯ was minimum at φ = 70° in the 20m × 20m × 20m environment for the LOS condition. However, Δ p ¯ with φ = 90° in both 5m × 5m × 5m and 20m × 20m × 20m environments were very close to the minimum values occurring at φ = 80° and 70° in the 5m × 5m × 5m and 20m × 20m × 20m environment, respectively. Table 6 shows that Δ p ¯ was minimum φ = 90° in both 5m × 5m × 5m and 20m × 20m × 20m environments for the NLOS condition. The variations of Δ p ¯ with the angle φ for the LOS and NLOS conditions in both Table 5 and Table 6 are illustrated in Figure 10a,b, respectively. As for the angle ψ , it is shown in Table 7 that Δ p ¯ was minimum is minimum at ψ = 20° in both 5m × 5m × 5m and 20m × 20m × 20m environments for the LOS condition. However, Δ p ¯ with ψ = 0° in both 5m × 5m × 5m and 20m × 20m × 20m environments were also very close to the minimum values occurring at ψ = 20° in both environments. It is shown in Table 8 that Δ p ¯ was minimum at ψ = 10° in the 5m × 5m × 5m environment while Δ p ¯ was minimum at φ = 20° in the 20m × 20m × 20m environment for the NLOS condition. However, Δ p ¯ with ψ = 0° in both 5m × 5m × 5m and 20m × 20m × 20m environments were very close to the minimum values occurring at ψ = 10° and 20° in the 5m × 5m × 5m and 20m × 20m × 20m environment, respectively. The variations of Δ p ¯ with the angle ψ for the LOS and NLOS conditions in both Table 7 and Table 8 are illustrated in Figure 11a and Figure 11b, respectively. The suggestion based on (36) that the anchor installation arrangement with φ = 90° and ψ = 0° leads to reasonably small positioning errors for both the LOS and NLOS conditions is thus justified from Table 5, Table 6, Table 7 and Table 8.

5. Experiments

After the computer simulations shown in the previous section, an experiment was conducted to further verify the results obtained in 3D positioning for both the LOS and NLOS conditions. The 3D positioning experiment for the LOS condition was conducted in the environment where no obstacle was placed between the mobile node and every anchor. The experiment was implemented in a corridor outside of the laboratory in an environment 6m × 6m × 3.6m. The installation positions of the four anchors are shown in Figure 12a, where anchors q 1 and q 2 were installed on the wall with fixed positions. In contrast, the positions of q 3 and q 4 are adjustable for the positioning evaluations with different angles φ and ψ . The mobile node moves in the experiment environment with 0.52m above the ground. The illustration of the experiment environment and the moving trajectories of adjustable anchors q 3 and q 4 are shown in Figure 12b. The mobile node calculates the 3D positioning with EKF [45] and records the positions with a sampling frequency of 100 Hz. The RMSPEs Δ p ¯ were calculated for every 600 samples.
In the first LOS experiment, the anchors q 1 and q 2 were installed at the positions with coordinates [ 0 , 0 , 2.52 ] T and [ 0 , 5.79 , 2.52 ] T respectively, as shown in Figure 12b. It resulted in k being equal to 5.79m. The installation positions of anchors q 3 were adjusted following the trajectory shown in Figure 12b with φ varying from 10° to 90°. Note that the z coordinate of q 3 was kept constant, 2.52m, while changing the installation positions. Therefore, l   = 5.79m. The installation position of q 4 was also adjusted following the trajectory shown in Figure 12b with ψ varying from 0° to 65°. It is shown in Figure 12b that m remained a constant, 2.49m, while ψ changed from 0° to 65°. The mobile node was parked at the position 2.45 ,   3.74 ,   0.52 T to measure the mobile node’s 3D position for different anchor positions. The RMSPEs Δ p ¯ were calculated with φ varying from 10° to 90° and yet ψ was set to be 0°. Similarly, the RMSPEs Δ p ¯ were calculated with ψ ranging from 0° to 65° and yet φ was set to be 90°. The RMSPEs Δ p ¯ corresponding to different combinations of φ and ψ are shown in Table 9. Table 9 shows that Δ p ¯ was minimum if the installation position of q 3 is adjusted so that φ = 90° under the condition ψ = 0°; or if the installation of q 4 is adjusted so that ψ = 0° under the condition φ = 90°.
In the second LOS experiment, the mobile node was still parked at the position 2.45 ,   3.74 ,   0.52 T , and both q 3 and q 4 were set at the position so that φ = 90° and ψ = 0°. The RMSPEs Δ p ¯ were measured and calculated while changing both k and l from 1.79m to 5.79m, and changing m from 1.19m to 2.19m. Note that Δ p ¯ was calculated by changing only one length of vector at a time with the lengths of two other vectors remaining constants. The experiment results are shown in Table 10. Table 10 shows that the Δ p ¯ decreased if the lengths of 3 vectors k , l , and m were appropriately increased.
The 3D positioning experiment for the NLOS condition was conducted in the same environment as in the LOS experiment, except that several stationary and moving obstacles were placed between the mobile node and the anchors. Figure 13a shows that a stationary obstacle A with the dimension 0.7m × 0.4m × 1.65m was placed between the mobile node and the anchor q 4 while a stationary obstacle B with the dimension 0.6m × 0.25m × 1.65m was placed between the mobile node and the anchor q 2 . In addition to the stationary obstacles, a pedestrian was arranged to walk back and forth following a straight line between the mobile node and the anchor q 2 shown in Figure 13a. The measured distance r i between the mobile node and every ith anchor q i , i = 1, …, 4, for the LOS and NLOS conditions are compared in Figure 13b and Figure 13c, respectively. From the measured distances r 1 , , r 4 compared in Figure 13b,c, it is shown that both stationary and moving obstacles did affect the distance measurement. The same 3D positioning experiment for the NLOS condition was conducted as for the LOS condition. The RMSPEs Δ p ¯ corresponding to different combinations of φ and ψ for the NLOS condition are shown in Table 11. Table 11 shows that Δ p ¯ was minimum if the installation position of q 3 is adjusted so that φ = 90° under the condition ψ = 0°. The RMSPE Δ p ¯ was also minimum if the installation of q 4 is adjusted so that ψ = 0° under the condition φ = 90°. The RMSPEs Δ p ¯ were calculated by changing only one vector length at a time with the lengths of two other vectors remaining constants. The experiment results are shown in Table 12. Table 12 shows that the Δ p ¯ decreased if the lengths of three vectors k , l , and m were appropriately increased.
The suggestion based on (36) that the anchor installation arrangement with φ = 90° and ψ = 0° leads to reasonably small 3D positioning errors for both the LOS and NLOS conditions is thus justified from the experimental results in Table 9 and Table 11. Moreover, the suggestion that installing the anchors as separated as possible leads to reasonably small 3D positioning errors is also justified from the experimental results in Table 10 and Table 12.

6. Conclusions

The upper bounds of variance of both 2D and 3D positioning errors were derived in this paper. The mathematical model for the variance of positioning errors shed light on some tips regarding installation positions of anchors. The positioning errors can be reduced through an appropriate arrangement of anchors’ installation positions according to the mathematical model derived in the paper. Computer simulations and practical experiments were conducted to verify the suggestions of anchor installation positions obtained from the mathematical model.
The positioning accuracy of UWB sensors has attracted a lot of attention lately due to their cost-effective positioning performance. The positioning errors analyzed in the paper were limited to the 2D or 3D positions only. However, the positioning errors can be expanded to include the linear and angular velocities and accelerations if the UWB sensors are combined with the inertial measurement units (IMU). Further works can be performed to derive the mathematical model for the “expanded” positioning errors so that more installation suggestions regarding improving positioning accuracy can be obtained from the model. Further works can also be performed to derive the lower bound of variance of positioning errors based on the proposed mathematical model. The range of positioning errors associated with a combination of anchor installation positions can be estimated with more accuracy if both upper and lower bounds of positioning errors are calculated.

Author Contributions

Conceptualization, L.Y. (Leehter Yao); methodology, L.Y. (Leehter Yao) and Y.-W.W.; software, L.Y. (Lei Yao); validation, L.Y. (Lei Yao); writing—original draft preparation, L.Y. (Lei Yao); writing— review and editing, L.Y. (Leehter Yao); supervision, L.Y. (Leehter Yao); project administration, L.Y. (Leehter Yao); funding acquisition, L.Y. (Leehter Yao). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of Science and Technology, Taiwan, grant number MOST 107-2221-E-027-086-MY3.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, X.; Wang, J.; Liu, C.; Zhang, L.; Li, Z. Integrated wifi/pdr/smartphone using an adaptive system noise extended kalman filter algorithm for indoor localization. ISPRS Int. J. Geo-Inf. 2016, 5, 8. [Google Scholar] [CrossRef] [Green Version]
  2. Zuin, S.; Calzavara, M.; Sgarbossa, F.; Persona, A. Ultra wide band indoor positioning system: Analysis and testing of an IPS technology. IFAC 2018, 51, 1488–1492. [Google Scholar]
  3. Mendoza-Silva, G.M.; Torres-Sospedra, J.; Huerta, J. A more realistic error distance calculation for indoor positioning systems accuracy evaluation. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017; pp. 1–8. [Google Scholar]
  4. Anagnostopoulos, G.G.; Deriaz, M.; Gaspoz, J.; Konstantas, D.; Guessous, I. Navigational needs and requirements of hospital staff: Geneva university hospitals case study. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017; pp. 1–8. [Google Scholar]
  5. Holm, S. Ultrasound positioning based on time-of-flight and signal strength. In Proceedings of the 2012 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sydney, NSW, Australia, 13–15 November 2012; pp. 1–6. [Google Scholar]
  6. Mohebbi, P.; Stroulia, E.; Nikolaidis, I. Sensor-data fusion for multi-person indoor position estimation. Sensors 2017, 17, 2377. [Google Scholar] [CrossRef] [Green Version]
  7. Castillo-Cara, M.; Lovon-Melgarejo, J.; Bravo-Rocca, G.; Orozco-Barbosa, L.; Garcia-Varea, I. An empirical study of the transmission power setting for Bluetooth-based indoor localization mechanisms. Sensors 2017, 17, 1318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Zhou, Y.D.; Wei, M.X.; Zhuang, H. Research on ZigBee indoor technology positioning based on RSSI. Procedia Comput. Sci. 2019, 154, 424–429. [Google Scholar]
  9. Diallo, A.; Lu, Z.; Zhao, X. Wireless indoor localization using passive RFID tags. Procedia Comput. Sci. 2019, 155, 210–217. [Google Scholar] [CrossRef]
  10. Xiao, W.; Ni, W.; Toh, Y.K. Integrated Wi-Fi fingerprinting and inertial sensing for indoor positioning. In Proceedings of the 2011 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Guimaraes, Portugal, 21–23 September 2011; pp. 1–6. [Google Scholar]
  11. Jiménez Ruiz, A.R.; Seco Granja, F. Comparing Ubisense, BeSpoon, and DecaWave UWB position systems: Indoor performance analysis. IEEE Trans. Instrum. Meas. 2017, 66, 2106–2117. [Google Scholar] [CrossRef]
  12. Park, H.; Noh, J.; Cho, S. Three-dimensional positioning system using Bluetooth low-energy beacons. IJDSN 2016, 12. [Google Scholar] [CrossRef] [Green Version]
  13. He, Y.; Bilgic, A. Iterative least squares method for global positioning system. Adv. Radio Sci. 2011, 9, 203–208. [Google Scholar] [CrossRef] [Green Version]
  14. Khalaf-Allah, M. Particle filtering for three-dimensional TDoA-based positioning using four anchor nodes. Sensors 2020, 20, 4516. [Google Scholar] [CrossRef]
  15. Yao, L.; Wu, Y.A.; Yao, L.; Liao, Z.Z. An integrated IMU and UWB sensor based indoor positioning system. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017; pp. 1–8. [Google Scholar]
  16. Zeng, Z.; Liu, S.; Wang, L. UWB/IMU integration approach with NLOS identification and mitigation. In Proceedings of the 2018 52nd Annual Conference on Information Sciences and Systems (CISS), Princeton, NJ, USA, 21–23 March 2018; pp. 1–6. [Google Scholar]
  17. Zeng, Z.; Liu, S.; Wang, L. NLOS detection and mitigation for UWB/IMU fusion system based on EKF and CIR. In Proceedings of the 2018 IEEE 18th International Conference on Communication Technology (ICCT), Chongqing, China, 8–11 October 2018; pp. 376–381. [Google Scholar]
  18. Liu, F.; Li, X.; Wang, J.; Zhang, J. An adaptive UWB/MEMS-IMU complementary Kalman filter for indoor position in NLOS environment. Remote Sens. 2019, 11, 2628. [Google Scholar] [CrossRef] [Green Version]
  19. 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]
  20. 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]
  21. Nakamura, S.; Higashi, Y.; Masuda, A.; Miura, N. A positioning system and position control system of a quad-rotor applying Kalman filter to a UWB module and an IMU. In Proceedings of the 2020 IEEE/SICE International Symposium on System Integration (SII), Honolulu, HI, USA, 12–15 January 2020; pp. 747–752. [Google Scholar]
  22. Jachimczyk, B.; Dziak, D.; Kulesza, W.J. Customization of UWB 3D-RTLS based on the new uncertainty model of the AoA ranging technique. Sensors 2017, 17, 227. [Google Scholar] [CrossRef] [Green Version]
  23. Ferreira, A.G.; Fernandes, D.; Catarino, A.P.; Monteiro, J.L. Performance analysis of ToA-based positioning algorithms for static and dynamic targets with low ranging measurements. Sensors 2017, 17, 1915. [Google Scholar] [CrossRef]
  24. Zhang, K.; Shen, C.; Zhou, Q.; Wang, H.; Gao, Q.; Chen, Y. A combined GPS UWB and MARG position algorithm for indoor and outdoor mixed scenario. Clust. Comput. 2019, 22, 5965–5974. [Google Scholar] [CrossRef]
  25. Ding, G.; Lu, H.; Bai, J.; Qin, X. Development of a high precision UWB/Vision-based AGV and control system. In Proceedings of the 2020 5th International Conference on Control and Robotics Engineering (ICCRE), Osaka, Japan, 24–26 April 2020; pp. 99–103. [Google Scholar]
  26. Song, Y.; Guan, M.; Tay, W.P.; Law, C.L.; Wen, C. UWB/LiDAR fusion for cooperative range-only SLAM. In Proceedings of the 2019 International Conference on Robotics and Automation (ICRA), Montreal, QC, Canada, 20–24 May 2019; pp. 6568–6574. [Google Scholar]
  27. Renaudin, V.; Susi, M.; Lachapelle, G. Step length estimation using handheld inertial sensors. Sensors 2012, 12, 8507–8525. [Google Scholar] [CrossRef]
  28. Song, Y.; Hsu, L.T. Tightly coupled integrated navigation system via factor graph for UAV indoor localization. Aerosp. Sci. Technol. 2021, 108. [Google Scholar] [CrossRef]
  29. Wen, K.; Yu, K.; Li, Y.; Zhang, S.; Zhang, W. A new quaternion Kalman filter based foot-mounted IMU and UWB tightly-coupled method for indoor pedestrian navigation. IEEE Trans. Veh. Tech. 2020, 69, 4340–4352. [Google Scholar] [CrossRef]
  30. Hol, J.D.; Dijkstra, F.; Luinge, H.; Schon, T.B. Tightly coupled UWB/IMU pose estimation. In Proceedings of the 2009 IEEE International Conference on Ultra-Wideband, Vancouver, BC, Canada, 9–11 September 2009; pp. 688–692. [Google Scholar]
  31. Li, X.; Wang, Y.; Khoshelham, K. A robust and adaptive complementary Kalman filter based on mahalanobis distance for ultra wideband/inertial measurement unit fusion positioning. Sensors 2018, 18, 3435. [Google Scholar] [CrossRef] [Green Version]
  32. 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. [Google Scholar] [CrossRef] [Green Version]
  33. Shen, Y.; Win, M.Z. Fundamental limits of wideband localization— Part I: A general framework. IEEE Trans. Inf. Theory 2010, 56, 4956–4980. [Google Scholar] [CrossRef]
  34. Miraglia, G.; Maleki, K.N.; Hook, L.R. Comparison of two sensor data fusion methods in a tightly coupled UWB/IMU 3-D localization system. In Proceedings of the 2017 International Conference on Engineering, Technology and Innovation (ICE/ITMC), Madeira, Portugal, 27–29 June 2017; pp. 611–618. [Google Scholar]
  35. Xu, Y.; Shmaliy, Y.S.; Chen, X.; Li, Y.; Ma, W. Robust inertial navigation system/ultra wide band integrated indoor quadrotor localization employing adaptive interacting multiple model-unbiased finite impulse response/Kalman filter estimator. Aerosp. Sci. Technol. 2020, 98, 105683. [Google Scholar] [CrossRef]
  36. Zhang, R.; Shen, F.; Liang, Y.; Zhao, D. Using UWB aided GNSS/INS integrated navigation to bridge GNSS outages based on optimal anchor distribution strategy. In Proceedings of the 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), Portland, OR, USA, 20–23 April 2020; pp. 1405–1411. [Google Scholar]
  37. Tahsin, M.; Sultana, S.; Reza, T.; Hossam-E-Haider, M. Analysis of DOP and its preciseness in GNSS position estimation. In Proceedings of the 2015 International Conference on Electrical Engineering and Information Communication Technology (ICEEICT), Savar, Bangladesh, 21–23 May 2015; pp. 1–6. [Google Scholar]
  38. Monica, S.; Ferrari, G. UWB-based localization in large indoor scenarios: Optimized placement of anchor nodes. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 987–999. [Google Scholar] [CrossRef]
  39. Decawave. The Implementation of Two-Way Ranging with the DW 1000; Version 2.2; DecaWave Limited: Dublin, Ireland, 2015. [Google Scholar]
  40. Peter, K.; Matjaž, V.; Marko, M. Distance measurements in UWB-radio localization systems corrected with a feedforward neural network model. Sensors 2021, 21, 2294. [Google Scholar]
  41. Larsson, E.G. Cram’er-Rao bound analysis of distribution positioning in sensor networks. IEEE Trans. Signal Process. Lett. 2004, 11, 334–337. [Google Scholar] [CrossRef]
  42. Qi, Y.; Kobayashi, H.; Suda, H. Analysis of wireless geolocation in a non-line-of-sight environment. IEEE Trans. Wireless Commun. 2006, 5, 672–681. [Google Scholar]
  43. Nurminen, H.; Ardeshiri, T.; Piché, R.; Gustafsson, F.A. NLOS-robust TOA positioning filter based on a skew-t measurement noise model. In Proceedings of the 2015 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Banff, AB, Canada, 13–16 October 2015; pp. 1–7. [Google Scholar]
  44. Shen, Y.; Wymeersch, H.; Win, M.Z. Fundamental limits of wideband localization- part II: Cooperative networks. IEEE Trans. Inf. Theory 2010, 56, 4981–5000. [Google Scholar] [CrossRef] [Green Version]
  45. Khan, R.; Khan, S.U.; Khan, S.; Khan, M.U.A. Localization performance evaluation of extended kalman filter in wireless sensors network. Procedia Comput. Sci. 2014, 32, 117–124. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The UWB sensors used in this paper. (a) anchors; (b) mobile node.
Figure 1. The UWB sensors used in this paper. (a) anchors; (b) mobile node.
Sensors 21 05731 g001
Figure 2. Positioning of a mobile node in a 3D space with 4 anchors.
Figure 2. Positioning of a mobile node in a 3D space with 4 anchors.
Sensors 21 05731 g002
Figure 3. Positioning of a mobile node in a 2D space with three anchors.
Figure 3. Positioning of a mobile node in a 2D space with three anchors.
Sensors 21 05731 g003
Figure 4. Projection of installation positions of anchors q 1 ,   q 2   and   q 3 onto the X−Y plane and the moving trajectory of the mobile node. (a) The environment of 5m × 5m. (b) The environment of 20m × 20m.
Figure 4. Projection of installation positions of anchors q 1 ,   q 2   and   q 3 onto the X−Y plane and the moving trajectory of the mobile node. (a) The environment of 5m × 5m. (b) The environment of 20m × 20m.
Sensors 21 05731 g004
Figure 5. Variation of 2D RMSPE Δ p ¯ with angles φ for the LOS and NLOS conditions. (a) The environment of 5m × 5m. (b) The environment of 20m × 20m.
Figure 5. Variation of 2D RMSPE Δ p ¯ with angles φ for the LOS and NLOS conditions. (a) The environment of 5m × 5m. (b) The environment of 20m × 20m.
Sensors 21 05731 g005
Figure 6. Variation of 2D RMSPE Δ p ¯ with different anchor distances k and l . (a) Variation of Δ p ¯ with k provided that l = 20 m ; (b) Variation of Δ p ¯ with l provided that k = 20 m .
Figure 6. Variation of 2D RMSPE Δ p ¯ with different anchor distances k and l . (a) Variation of Δ p ¯ with k provided that l = 20 m ; (b) Variation of Δ p ¯ with l provided that k = 20 m .
Sensors 21 05731 g006
Figure 7. Comparison of upper bound of 2D positioning errors estimated by our approach and the approach in [38] for the LOS condition. (a) Variation of Δ p i u with samples; (b) Variation of Δ p ¯ u with distance k ; (c) Variation of Δ p ¯ u with angle φ .
Figure 7. Comparison of upper bound of 2D positioning errors estimated by our approach and the approach in [38] for the LOS condition. (a) Variation of Δ p i u with samples; (b) Variation of Δ p ¯ u with distance k ; (c) Variation of Δ p ¯ u with angle φ .
Sensors 21 05731 g007
Figure 8. Comparison of upper bound of 2D positioning errors estimated by our approach and the approach in [38] for the NLOS condition. (a) Variation of Δ p i u with samples; (b) Variation of Δ p ¯ u with distance k ; (c) Variation of Δ p ¯ u with angle φ .
Figure 8. Comparison of upper bound of 2D positioning errors estimated by our approach and the approach in [38] for the NLOS condition. (a) Variation of Δ p i u with samples; (b) Variation of Δ p ¯ u with distance k ; (c) Variation of Δ p ¯ u with angle φ .
Sensors 21 05731 g008
Figure 9. Installation positions of anchors q 1 ,   q 2   q 3   and   q 4 and moving trajectory of the mobile node. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Figure 9. Installation positions of anchors q 1 ,   q 2   q 3   and   q 4 and moving trajectory of the mobile node. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Sensors 21 05731 g009
Figure 10. Variation of 3D RMSPE Δ p ¯ with angles φ for the LOS and NLOS conditions. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Figure 10. Variation of 3D RMSPE Δ p ¯ with angles φ for the LOS and NLOS conditions. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Sensors 21 05731 g010
Figure 11. Variation of 3D RMSPE Δ p ¯ with angles ψ for the LOS and NLOS conditions. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Figure 11. Variation of 3D RMSPE Δ p ¯ with angles ψ for the LOS and NLOS conditions. (a) The environment of 5m × 5m × 5m. (b) The environment of 20m × 20m × 20m.
Sensors 21 05731 g011
Figure 12. Experiment of measuring mobile node’s positioning errors. (a) Photograph showing the positions of the four anchors and the mobile node for the LOS condition; (b) Illustration of the experimental environment and the moving trajectory of the adjustable anchors.
Figure 12. Experiment of measuring mobile node’s positioning errors. (a) Photograph showing the positions of the four anchors and the mobile node for the LOS condition; (b) Illustration of the experimental environment and the moving trajectory of the adjustable anchors.
Sensors 21 05731 g012
Figure 13. Experiment of measuring the mobile node’s positioning errors. (a) Photograph showing the positions of the four anchors and obstacles; (b) Measured distances between the mobile node and four anchors for the LOS condition; (c) Measured distances between the mobile node and four anchors for the NLOS condition.
Figure 13. Experiment of measuring the mobile node’s positioning errors. (a) Photograph showing the positions of the four anchors and obstacles; (b) Measured distances between the mobile node and four anchors for the LOS condition; (c) Measured distances between the mobile node and four anchors for the NLOS condition.
Sensors 21 05731 g013
Table 1. The 2D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m and 20m × 20m environments for the LOS condition.
Table 1. The 2D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m and 20m × 20m environments for the LOS condition.
φ   ( d e g ) 102030405060708090100110120130140150160
5m × 5m2.0441.0220.7830.6200.5560.5390.5400.5360.5260.6060.6770.7230.8911.1481.4582.104
20m × 20m0.6540.3250.2130.1760.1550.1500.1500.1490.1490.1840.1950.2270.2700.3330.4480.684
Table 2. The 2D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m and 20m × 20m environments for the NLOS condition.
Table 2. The 2D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m and 20m × 20m environments for the NLOS condition.
φ   ( d e g ) 102030405060708090100110120130140150160
5m × 5m4.1922.0911.5131.3231.1551.1171.1121.1131.1141.2611.4021.4771.7452.2102.8914.219
20m × 20m1.3230.6750.5300.4480.4020.3750.3680.3730.3710.4120.4530.5380.6440.8131.0851.699
Table 3. The 2D RMSPE Δ p ¯ (m) with varying anchor distances k and l for the LOS condition.
Table 3. The 2D RMSPE Δ p ¯ (m) with varying anchor distances k and l for the LOS condition.
Distances (m)48121620
|k’| (m)0.5050.2590.1940.1700.149
|l’| (m)0.5010.2660.1920.1670.149
Table 4. The 2D RMSPE Δ p ¯ (m) with varying anchor distances k and l for the NLOS condition.
Table 4. The 2D RMSPE Δ p ¯ (m) with varying anchor distances k and l for the NLOS condition.
Distances (m)48121620
|k’| (m)1.0390.5780.4630.4140.371
|l’| (m)1.0400.5680.4630.3940.371
Table 5. The 3D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the LOS condition.
Table 5. The 3D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the LOS condition.
φ   ( d e g ) 102030405060708090100110120130140150160
(5m)32.5931.4031.0580.9230.8720.8470.8380.8350.8370.9030.9571.0431.1771.4001.8052.622
(20m)30.7970.4270.3180.2700.2520.2400.2380.2390.2390.2590.2810.3130.3730.4670.6270.961
Table 6. The 3D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the NLOS condition.
Table 6. The 3D RMSPE Δ p ¯ (m) with varying φ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the NLOS condition.
φ   ( d e g ) 102030405060708090100110120130140150160
(5m)34.9362.7252.0761.7941.6761.6391.6211.6141.5911.7191.8282.0392.3682.6493.5355.021
(20m)31.5390.8130.5990.5100.4600.4540.4520.4540.4500.5100.5310.6030.6870.8301.0501.533
Table 7. The 3D RMSPE Δ p ¯ (m) with varying ψ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the LOS condition.
Table 7. The 3D RMSPE Δ p ¯ (m) with varying ψ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the LOS condition.
ψ   ( d e g ) 80706050403020100−10−20−30−40−50−60−70
(5m)32.6791.4371.0820.9250.8480.8300.8220.8340.8280.9100.9941.1131.2741.5592.0022.967
(20m)30.8300.4460.3300.2750.2520.2370.2310.2310.2320.2650.2990.3450.4310.5290.7291.120
Table 8. The 3D RMSPE Δ p ¯ (m) with varying ψ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the NLOS condition.
Table 8. The 3D RMSPE Δ p ¯ (m) with varying ψ in the 5m × 5m × 5m and 20m × 20m × 20m environments for the NLOS condition.
ψ   ( d e g ) 80706050403020100−10−20−30−40−50−60−70
(5m)35.1392.7502.1741.8171.6671.6501.6151.5751.6141.6551.9622.1562.5022.9853.9115.593
(20m)31.4150.7420.5470.4740.4460.4410.4350.4410.4380.5160.5560.6180.7200.9061.1691.780
Table 9. The 3D RMSPE Δ p ¯ with varying angles φ and ψ for the LOS condition.
Table 9. The 3D RMSPE Δ p ¯ with varying angles φ and ψ for the LOS condition.
φ ( deg ) ψ ( deg ) Δ p ¯   ( m )
group 19000.0743
7000.1366
5500.2214
4000.2965
2500.6745
1001.4545
group 29000.0743
90200.1404
90350.2136
90500.3018
90650.5648
Table 10. The 3D RMSPE Δ p ¯ with varying anchor distances k , l , and m for the LOS condition.
Table 10. The 3D RMSPE Δ p ¯ with varying anchor distances k , l , and m for the LOS condition.
|k| (m)|l| (m)|m| (m) Δ p ¯   ( m )
group 15.795.792.190.0743
5.794.792.190.1115
5.793.792.190.1218
5.792.792.190.1807
5.791.792.190.3169
group 25.795.792.190.0743
5.795.791.690.1850
5.795.791.190.2852
group 35.795.792.190.0743
4.795.792.190.1239
3.795.792.190.1544
2.795.792.190.1879
1.795.792.190.2989
Table 11. The 3D RMSPE Δ p ¯ environments for the NLOS condition. φ and ψ for the NLOS condition.
Table 11. The 3D RMSPE Δ p ¯ environments for the NLOS condition. φ and ψ for the NLOS condition.
φ ( deg ) ψ ( deg ) Δ p ¯   ( m )
group 19000.5452
7000.6355
5500.6758
4000.9418
2501.6660
1002.9914
group 29000.5452
90200.5659
90350.7983
90501.0277
90651.1303
Table 12. The 3D RMSPE Δ p ¯ with varying anchor distances k , l , and m for the NLOS condition.
Table 12. The 3D RMSPE Δ p ¯ with varying anchor distances k , l , and m for the NLOS condition.
|k| (m)|l| (m)|m| (m) Δ p ¯   ( m )
group 15.795.792.190.5452
5.794.792.190.7263
5.793.792.190.7575
5.792.792.190.8255
5.791.792.191.0081
group 25.795.792.190.5452
5.795.791.690.8880
5.795.791.191.0658
group 35.795.792.190.5452
4.795.792.190.7302
3.795.792.190.7333
2.795.792.190.7593
1.795.792.190.9898
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yao, L.; Yao, L.; Wu, Y.-W. Analysis and Improvement of Indoor Positioning Accuracy for UWB Sensors. Sensors 2021, 21, 5731. https://doi.org/10.3390/s21175731

AMA Style

Yao L, Yao L, Wu Y-W. Analysis and Improvement of Indoor Positioning Accuracy for UWB Sensors. Sensors. 2021; 21(17):5731. https://doi.org/10.3390/s21175731

Chicago/Turabian Style

Yao, Leehter, Lei Yao, and Yeong-Wei Wu. 2021. "Analysis and Improvement of Indoor Positioning Accuracy for UWB Sensors" Sensors 21, no. 17: 5731. https://doi.org/10.3390/s21175731

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