UWB and MEMS IMU Integrated Positioning Algorithm for a Work-Tool Tracking System

Abstract

In this paper, we address a system that can accurately locate and monitor work tools in a complex assembly process, such as automotive production. Our positioning monitoring system is positioned by a combined sensor of the UWB module and the MEMS IMU (inertial measuring unit) sensor based on the extended Kalman filter. The MEMS IMU sensor provides the positioning calibration information. The proposed method incorporates IMU and UWB positioning to compensate for errors that can only occur in UWB positioning through the extended Kalman filter (EKT). This EKT is improved by the error dynamic equation derived from the sparse state-space matrix. Also, the proposed method computes the transmission time and distance between the tag and anchor of the UWB module by the TWR (two-way range) system. The tag of a mobile node, which is attached to a moving tool, measures the position of the work tool and transmits the position coordinate data to the anchor. Here, the proposed method uses the trilateration localization method by the confidence distance compensation to prevent the distance error by obstacles and changes in the indoor environment. Experimental results verified that the proposed method confirms whether a specific tool is accurately used according to the prescribed regulations and has more positioning accuracy than the conventional methods.

1. Introduction

The automobile industry is the largest manufacturing industry in the world with revenues of over 1 trillion US dollars and more than 10 million employees handling roughly 20,000 components in complex assembly processes. The automobile assembly process is labor intensive, in which approximately 3000 units of interior material, dashboard, seat, windshield, engine, transmission, and axle are assembled in a painted automobile body, thus having the lowest automation rate compared to other processes [1].

The current automobile assembly process is not monitored in real time, which causes difficulty in identifying a missing operation in each process. Therefore, if the location of the tool currently being worked on can be determined, it can improve the productivity and quality in the manufacturing process by identifying missing or redundant tasks to prevent manufacturing detects [2,3].

Herein, we propose a position tracking system of a tool using a UWB (ultra-wide band) and a MEMS (microelectro-mechanical system) IMU (inertial measurement unit) sensor for realtime monitoring of the position of tools. Recently, extensive research has been conducted on real-time position recognition systems using the UWB technology in which Wi-Fi uses a bandwidth of approximately 20 MHz, whereas UWB uses a wide bandwidth of several gigahertz; hence, UWB requires ten times less power than Wi-Fi but has a short transmission distance of 10 m [4,5]. Furthermore, an IMU is a device for measuring the speed, direction, gravitational force, and acceleration of a moving object by using an accelerometer, an angular velocity meter, a geomagnetic meter, and an altimeter, respectively [6].

We have previously studied the Kalman filter-based positioning system using MEMS IMU and UWB sensors [7]. However, any work tools are positioned out of 200 mm because of the linear system of the Kalman filter. This system should be improved by an advanced Kalman filter, such as the extended Kalman filter, complementary Kalman filter and the accurate TWR system, to prevent the distance by obstacles in the indoor environment. In this study, a UWB communication module was used to measure the position of a work tool, and the MEMS IMU sensor was used to provide positioning calibration data for reducing the measurement error of the UWB. The output data of the IMU sensor or UWB contain noises in addition to positioning data because of electrical or mechanical reasons.

Our study has the following features about UWB and MEMS IMU integrated positioning algorithm for accurate work-tool tracking. Firstly, the accurate position, velocity, and attitude are measured through the MEMS IMU sensor derived from the inertial navigation formula in 3 coordinate frames of indoor positioning frame, local navigation frame, EGI frame. Secondly, the distance between an anchor node and a tag node by the TWR method is computed by the confidence-distance-compensated trilateration localization method to prevent the distance error by obstacles or changes in the indoor environment. Thirdly, the improved EKR (extended Kalman filter) by the error dynamic equation is derived from 8 × 8 sparse state-space matrix used for the integrated positioning algorithms of MEMS IMU and UWB sensors.

Any work tools that have deviated from the set operation path can be identified in real time using the position tracking system of a tool proposed in this study, thus enabling any omitted operations to be identified; the positioning accuracy was found to be 200 mm or less through multiple experiments, which indicates that the location of the tools can be accurately identified when the operation points are spaced more than 200 mm apart.

In Section 2, the proposed position tracking system based on a composite sensor is introduced. In Section 3, the validity of the proposed position tracking system was verified through the basic unit testing results. In Section 4, the experiment results for evaluating the performance of the proposed position tracking system are presented.

2. Proposed Position Tracking System Based on Composite Sensors

The UWB wireless technology is applied in telecommunication using a significantly wide frequency band over several gigahertz in baseband without using a carrier wave. This system does not entail mutual interference with existing mobile communication systems by using a very narrow pulse of several nanoseconds based on a very low spectrum power, such as the noise of an existing wireless system. The proposed position tracking system based on a composite sensor consists of a positioning tag attached on a moving object for measuring the object position and an anchor installed on the predetermined position for UWB communication with the tag, as shown in Figure 1a.

2.1. Positioning Tag and Anchor Based on UWB and MEMS IMU Sensor

Herein, a position tracking system is proposed in which the UWB communication module and the MEMS IMU sensor are integrated, as shown Figure 1b; the positioning tag consists of the UWB module [8], MEMS IMU sensor, MCU, and battery module. The MEMS IMU sensor is removed while either Wi-Fi or Ethernet module is included for the anchor, because it is installed on a fixed position and does not need to measure the position data [9,10]. In particular, the UWB modules [8] of the tag and anchor transmit signals between each other to measure the transmission time for estimating the distance between the tag and anchor. The MEMS IMU sensor provides the positioning calibration information on the basis of the accelerometer, angular velocity meter, geomagnetic meter, and altimeter. The tag, which is a moving node, is attached to a moving tool to measure the three-dimensional position of the tool and to deliver the coordinate information to the anchor, while notifying the operator as to whether the tool is performing the correct operation through alarming via an LED lamp. The anchor is installed in a fixed position and determines the distance between the tags and collects position data and alarm information from multiple tags. Figure 1c shows the UWB sensor board of the PCB assembly that we implemented and the design of the tag and anchor.

2.2. MEMS IMU Positioning Algorithm

The MEMS IMU sensor is an inertial sensor with Gyroscope, Accelerometer, and Magnetometer for measurements on three axes (Roll $\varphi$, Pitch $\theta$, and Yaw $\psi$). We apply the MEMS IMU positioning algorithm by Yao [11]. The errors of speed and acceleration data measured by an IMU sensor are complemented by a UWB communication system, which makes the position measurement accurate. A tag sensor of UWB, accelerometers and gyros are mounted on the platform body of B frame and 3 anchors of UWB are located on any positions.

Given that $p_p\left(t_0\right)$ and $v_p\left(t_0\right)$ are the initial platform position and velocity to E frame, the position $p_p\left(t\right)$ and the velocity $v_p\left(t\right)$ for the platform motion to an indoor positioning coordinate frame (E frame) can be obtained by the computable accelerations $a_m^N$ for a local navigation coordinate frame (N frame).

$$p_p\left(t\right)=p_p\left(t_0\right)+{{\displaystyle \int }}_{t_0}^t{v}_p\left(\tau \right)d\tau \mathrm{and}{v}_p\left(t\right)=v_p\left(t_0\right)+{{\displaystyle \int }}_{t_0}^t{a}_m^N\left(\tau \right)d\tau$$

(1)

Here, $a_m^N$ is computed by the following inertial navigation formula;

$$a_m^N={\widehat{C}}_B^N{a}_m^B-{\omega }_{I\to N}^N{v}_m^N+g^N$$

(2)

where $a_m^B$ is the measured body accelerations and $g^N$ is the gravity vector. $v_m^N$ is the computed translational velocity of the platform and ${\omega }_{I\to N}^N$ is the angular velocity vector representing the transport rate vector from the EGI (Earth-Center-Inertial) frame (I frame) to the navigation coordinate frame (N frame).

$${\omega }_{I\to N}^N=\left[\begin{array}{ccc}0& -{\omega }_z& {\omega }_y\\ {\omega }_z& 0& -{\omega }_x\\ -{\omega }_y& {\omega }_x& 0\end{array}\right]$$

(3)

$\left({\omega }_x,{\omega }_y,{\omega }_z\right)$ is the measured gyro output. ${\widehat{C}}_B^N$ is a 3 × 3 direction cosine matrix computed from $C_I^N$ and ${\widehat{C}}_B^I$,

$${\widehat{C}}_B^N=\left[C_I^N\right]\left[{\widehat{C}}_B^I\right]$$

(4)

which $C_I^N$ represents the attitude of I frame for N frame and $C_B^I$ represents the transported B frame for I frame.

$$C_I^N=\left[C_z\left(\pi \right)C_y\left(-lat+\frac{\pi }{2}\right)C_z\left(lon\right)\right]\left[C_z\left({\theta }_e\right)\right].$$

(5)

$${\widehat{C}}_B^I=C_B^I{\mathit{\omega }}_{I\to N}^N=C_B^I\left[\begin{array}{ccc}0& -{\omega }_z& {\omega }_y\\ {\omega }_z& 0& -{\omega }_x\\ -{\omega }_y& {\omega }_x& 0\end{array}\right].$$

(6)

where

$$C_y\left(\varphi \right)=\left[\begin{array}{ccc}\cos \left(\varphi \right)& 0& \sin \left(\varphi \right)\\ 0& 1& 0\\ -\sin \left(\varphi \right)& 0& \cos \left(\varphi \right)\end{array}\right]\mathrm{and} C_z\left(\psi \right)=\left[\begin{array}{ccc}\cos \left(\psi \right)& -\sin \left(\psi \right)& 0\\ \sin \left(\psi \right)& \cos \left(\psi \right)& 0\\ 0& 0& 1\end{array}\right].$$

(7)

Here, ${\theta }_e={\theta }_e\left(t_0\right)+{\omega }_e\left(t-t_0\right)$ by the earth’s rotation rate ${\omega }_e$. ($lat, lon$) are the latitude and longitude of motion platform. By using these equations, we compute Euler angles of roll $\varphi$, pitch $\theta$, and yaw $\psi$ and the position ${\widehat{p}}^N$ and velocity ${\widehat{v}}^N$ on N frame, as shown in Figure 2.

2.3. UWB Positioning Algorithm

The two-way ranging (TWR) method [12,13,14] was developed to calculate the distance between two nodes by measuring the turnaround time of a pulse transmitted between the two nodes in which an Anchor node (Node 1) sends the pulse to a Tag node (Node 2) that transmits the pulse back to Node 1, and then by subtracting the processing time at Node 2, as shown in Figure 2.

In general, the moving distance of light is calculated by multiplying the speed of light $c=3\times {10}^8 \mathrm{m}/\mathrm{s}$ and the measured time; thus, the distance between an anchor and a tag is expressed as $d$ and

$$d=c\times TOF$$

(8)

where $TOF$ (Time of Flight) is

$$TOF=\frac{t_{\mathrm{total}}-t_{\mathrm{r}}}{2}=\frac{\left(t_2-t_1\right)-t_{\mathrm{r}}}{2}.$$

(9)

as shown in Figure 3. Here, $TOF$ represents the time required for a signal to be transmitted between an anchor. It is calculated by $TOF$ and is calculated by time stamps $t_1$ of Tx and $t_2$ of Rx and the response time $t_{\mathrm{r}}$. In addition, $t_{\mathrm{r}}$ refers to the time required for a tag (Node 2) that has received the pulse signal to transmit the signal back. The TWR method has the advantage of simplifying the system through accurate positioning without synchronizing the nodes.

Such wireless communication-based distance measurement is difficult to accurately recognize positions when the distance errors occur due to obstacles or changes in the indoor environment. Typical trilateration localization methods to reduce the distance measurement error can be used by topology [15], probability [16], confidence [17], and geographical similarity [18]. We used the confidence-based distance compensation trilateration localization method [17] to calculate the distance between 3 anchors and a target tag.

Firstly, we selected the fixed positions of 3 anchors: $\left\{\left(x_i,y_i\right)\right|i\in \left[1, 3\right]\}$. The distances $\left\{d_i\right|i\in \left[1, 3\right]\}$ between 3 anchors and a tag can be computed by TOFs in Equation (1) and they satisfy the condition; ${\left(x_i-x_t\right)}^2+{\left(y_i-y_t\right)}^2=d_i^2$ for all $i\in \left[1, 3\right]$. Here, the intersection points $\left(x_{ij}^{+},y_{ij}^{+}\right),\left(x_{ij}^{-},y_{ij}^{-}\right)$ for the $i$th anchor and $j$th anchor can be derived from the simultaneous equation.

$$\{\begin{array}{c}{\left(x_i-x_t\right)}^2+{\left(y_i-y_t\right)}^2=d_i^2\\ {\left(x_j-x_t\right)}^2+{\left(y_j-y_t\right)}^2=d_j^2\end{array}$$

(10)

and

$$\{\begin{array}{c}{x}_{ij}^{+}=\frac{ab+{\epsilon }_i}{1+a^2}, x_{ij}^{-}=\frac{ab-{\epsilon }_i}{1+a^2}\\ y_{ij}^{+}=\frac{b+a{\epsilon }_i}{1+a^2}, y_{ij}^{-}=\frac{b-a{\epsilon }_i}{1+a^2}\end{array}$$

(11)

where ${\epsilon }_i=\sqrt{a^2{b}^2-\left(1+a^2\right)\left(b^2-d_i^2\right)}$. $a$ and $b$ are defined by

$$a=\frac{x_j-x_i}{y_j-y_i}$$

(12)

$$b=\frac{{\left(x_j-x_i\right)}^2+{\left(y_j-y_i\right)}^2-\left(d_j^2+d_i^2\right)}{2\left(y_j-y_i\right)}.$$

(13)

The distance $d$ of an anchor and a tag is updated to the new distance $d^{*}$ by the Cram$\acute{\mathrm{e}}$r–Rao lower bound (CRLB) $J$ of TOF estimate and the confidence interval $△d$.

$$d^{*}=d+△d$$

(14)

where

$$△d=1.96J/\sqrt{N},\left(N=3\right)$$

(15)

Given the positions of 3 anchors $\left\{\left(x_i,y_i\right),\left(x_j,y_j\right),\left(x_k,y_k\right)\right\}$, radii $\left\{r_i,r_j,r_k\right\}$ of circles of them, and intersection points $(x_{ij}^{+},y_{ij}^{+}),(x_{ij}^{-},y_{ij}^{-})$ of $i,j$ anchors in Equation (4), the optimal intersection point $(x_{ij}^{*},y_{ij}^{*})$ can be selected by the following inequality equation.

$$\{\begin{array}{c}(x_{ij}^{+},y_{ij}^{+})=(x_{ij}^{*},y_{ij}^{*}), \mathrm{if}{(d_k-\sqrt{{(x_{ij}^{+}-x_k)}^2+{(y_{ij}^{+}-y_k)}^2})}^2<{(d_k-\sqrt{{(x_{ij}^{-}-x_k)}^2+{(y_{ij}^{-}-y_k)}^2})}^2\\ (x_{ij}^{-},y_{ij}^{-})=(x_{ij}^{*},y_{ij}^{*}),\mathrm{otherwise}\end{array}$$

(16)

The optimal intersection points $\{(x_{ij}^{*},y_{ij}^{*}),(x_{ik}^{*},y_{ik}^{*}),(x_{jk}^{*},y_{jk}^{*})\}$ for 3 anchors are obtained by the above method, then the position of target tag node is calculated by the centroid-based trilateration.

$$\{\begin{array}{c}{\widehat{x}}_t=\frac{1}{3}\left(x_{ij}^{*}+x_{ik}^{*}+x_{jk}^{*}\right)\\ {\widehat{y}}_t=\frac{1}{3}\left(y_{ij}^{*}+y_{ik}^{*}+y_{jk}^{*}\right)\end{array}$$

(17)

By changing the position of 3 anchors, the position of the target tag node can be computed by the above method.

2.4. EKT-Based MEMS IMU and UWB Integrated Positioning Algorithm

Integrated positioning algorithms of MEMS IMU and UWB sensors have been studied by using various Kalman filters, such as the extended Kalman filter [11,19] and complementary Kalman filter [20,21]. We applied Yao’s integrated positioning algorithm with the extended Kalman filter (EKT) [11]. The Kalman filter of a recursive filter estimates the state of a linear dynamic system containing noises that can recursively process the measurement data containing measurement noises, thus being widely used across various fields, including computer vision, robotics, and radar. As shown in Figure 3, the extended Kalman filter consists of five steps, largely comprising a step for predicting the current state and a step for estimating a more accurate value, including the measurement value [6].

We consider a 2D positioning frame for a platform that moves on an x-y plane with a body-mounted IMU and an UWB tag sensor. This means that two orthogonal accelerometers measure platform body acceleration on the x and y axes on the fixed E frame, and measure body angular velocity on the I frame with one gyro along the z axis. The state variables consist of position errors, velocity errors, and attitude angle errors and additional state variables for IMU sensors consists biases ($b$) of accelerations and gyros and misalignment errors. Among them, 8 state error variables ${\mathbf{ϵ}}_z$ can be defined [11], as follows;

$${\mathbf{ϵ}}_z={\left[e_1,e_2,e_3,e_4,e_5,e_6,e_7,e_8\right]}^T$$

(18)

where

$$\begin{array}{cccc}{e}_1=x_p-{\widehat{x}}_p,& e_2=v_{xp}-{\widehat{v}}_{xp},& e_3=b_x-{\widehat{b}}_x,& e_4=y_p-{\widehat{y}}_p\end{array}$$

(19)

$$\begin{array}{cccc}{e}_5=v_{yp}-{\widehat{v}}_p,& e_6=b_y-{\widehat{b}}_y,& e_7=\psi -{\widehat{\psi }}_x,& e_8=b_x-{\widehat{b}}_x\end{array}.$$

(20)

By applying the small perturbations in the IMU positioning system, the error dynamic equation can be derived by the state-space form.

$${\widehat{x}}_z=F_s\left(t\right){\mathbf{ϵ}}_z+\mathbf{\eta }$$

(21)

Here, the 8 × 8 sparse state-space matrix $F_s\left(t\right)$ on any time can be defined by

$$F_s\left(t\right)=\left[\begin{array}{cccccccc}0& 1& 0& 0& 0& 0& 0& 0\\ 0& 0& -\cos \left(\psi \right)& 0& -{\omega }_z& \sin \left(\psi \right)& f_1(\psi ,{\mathbf{a}}^B,\mathbf{b})& {\widehat{v}}_{yp}\\ 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 1& 0& 0& 0\\ 0& -{\omega }_z& -\sin \left(\psi \right)& 0& 0& -\cos \left(\psi \right)& f_2(\psi ,a_x^B,a_y^B)& -{\widehat{v}}_{xp}\\ 0& 0& 0& 0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0& 0& 0& -1\\ 0& 0& 0& 0& 0& 0& 0& 0\end{array}\right],$$

(22)

where $f_1\left(\psi ,a^B,b\right)$ and $f_2\left(\psi ,a^B,b\right)$ are computed by the error of body accelerations $a^B$ and biases $b$ to the yaw $\psi$ on z axis as the following.

$$\left[\begin{array}{c}{f}_1\left(\psi ,a^B,b\right)\\ f_2\left(\psi ,a^B,b\right)\end{array}\right]=\left[\begin{array}{cc}-\sin \left(\psi \right)& \cos \left(\psi \right)\\ \cos \left(\psi \right)& -\sin \left(\psi \right)\end{array}\right]\left[\begin{array}{c}{a}_x^B-{\widehat{b}}_x\\ a_y^B-{\widehat{b}}_y\end{array}\right]$$

(23)

Next, the state-space parameter $\mathsf{\eta }$ is a vector for gyro measurement defined by

$$\mathsf{\eta }={\left[0,{\eta }_1,b_{ax},0,{\eta }_2,b_{ay},g_{ang},g_{rate}\right]}^T$$

(24)

where $\left(g_{ang},g_{rate}\right)$ are gyro angle random walk (ARW) and gyro rate random walk (RRW) and $\left(b_{ax},b_{ay}\right)$ are accelerometer biases stabilities. $\left({\eta }_1,{\eta }_2\right)$ are computed by the accelerometer noises $(n_{ax}$,$n_{ay})$ to the yaw $\psi$ on the z axis, as follows.

$$\left[\begin{array}{c}{\eta }_1\\ {\eta }_2\end{array}\right]=-\left[\begin{array}{cc}\cos \left(\psi \right)& -\sin \left(\psi \right)\\ \sin \left(\psi \right)& \cos \left(\psi \right)\end{array}\right]\left[\begin{array}{c}{n}_{ax}\\ n_{ay}\end{array}\right]+\left[\begin{array}{c}{\widehat{v}}_{yp}{g}_{ang}\\ -{\widehat{v}}_{xp}{g}_{ang}\end{array}\right]$$

(25)

$({\widehat{v}}_{xp},{\widehat{v}}_{yp})$ are $({{\displaystyle \int }}_0^t{a}_x^N\left(\tau \right)d\tau ,{{\displaystyle \int }}_0^t{a}_y^N\left(\tau \right)d\tau )$ and $({\widehat{x}}_{xp},{\widehat{x}}_{yp})=({{\displaystyle \int }}_0^t{\widehat{v}}_{xp}\left(\tau \right)d\tau ,{{\displaystyle \int }}_0^t{\widehat{v}}_{yp}\left(\tau \right)d\tau )$. $\psi ={{\displaystyle \int }}_0^t{\widehat{\omega }}_z\left(\tau \right)d\tau$. $\left(a_x^N,a_y^N\right)$ are derived from Equation (2).

$$\left[\begin{array}{c}{a}_x^N\\ a_y^N\end{array}\right]=\left[\begin{array}{cc}\cos \left(\psi \right)& -\sin \left(\psi \right)\\ \sin \left(\psi \right)& \cos \left(\psi \right)\end{array}\right]\left[\begin{array}{c}{a}_x^B-{\widehat{b}}_x\\ a_y^B-{\widehat{b}}_y\end{array}\right]-\left[\begin{array}{c}{\widehat{v}}_{xp}\\ -{\widehat{\omega }}_z{\widehat{v}}_{yp}\end{array}\right]$$

(26)

${\widehat{\omega }}_z$ is the error of the angle velocity to the bias on the z axis; ${\widehat{\omega }}_z={\omega }_z-{\widehat{b}}_z$.

The range measurement $z_z$ and estimated rand measurement ${\widehat{z}}_z$ for 3 anchors of UWB can be defined by the Euclidean distance of anchors and a target tag platform including an anchor height $h$ and the measurement noises $n={\left[n_1,n_2,n_3\right]}^T$ as follows:

$$z_z\left(t\right)=z+n\mathrm{for}z={\left[z_1,z_2,z_3\right]}^T$$

(27)

$${\widehat{z}}_z\left(t\right)=\widehat{z}+n \mathrm{for} \widehat{z}={\left[{\widehat{z}}_1,{\widehat{z}}_1,{\widehat{z}}_1\right]}^T$$

(28)

where

$$\left[\begin{array}{c}{z}_i\\ {\widehat{z}}_i\end{array}\right]=\left[\begin{array}{c}\sqrt{{\left(x_i-x_p\left(t\right)\right)}^2+{\left(y_i-y_p\left(t\right)\right)}^2+h^2}\\ \sqrt{{\left(x_i-{\widehat{x}}_p\left(t\right)\right)}^2+{\left(y_i-{\widehat{y}}_p\left(t\right)\right)}^2+h^2}\end{array}\right]+n_i$$

(29)

The extended Kalman filter residual μ, the error between two range measurements $z_z\left(t\right),{\widehat{z}}_z\left(t\right)$, can be approximated by the gradient matrix $H_z=\frac{\partial g\left(x\right)}{\partial x}$ and the state variables $x_z$ as follows:

$$\mathit{\mu }=z_z\left(t\right)-{\widehat{z}}_z\left(t\right)\approx \left[H_z\right]x_z+n$$

(30)

The position and velocity of platform are estimated from the platform and navigation state variable $x_z$ from IMU sensors and the range measurements $z_z$ of 3 anchors of UWB sensors, as shown in Figure 3.

Let us take a closer look at the positioning process by EKT (Extended Kalman filter) that linearizes the state equation, as shown in Figure 4. EKT assumes the differentiability of a state transition function instead of the model linearity. There is a closed-form expression that the estimated state ${\widehat{x}}_{k+1}$ can be computed by a nonlinear function $f$ of the previous state ${\widehat{x}}_k$, control $u_k$, noise $w_k$, and time $t$.

$${\widehat{x}}_{k+1}=f\left({\widehat{x}}_k,u_k, t\right)+w_k$$

(31)

The Jacobians of the predicted state to the previous state and the noise take simple forms when the noise goes linearly into the state update equation. In EKT, the measurement ${\widehat{z}}_k$ can be computed by a nonlinear function $h$ of the state $x_k$ and the measurement noise $v_k$.

$${\widehat{z}}_{k+1}=h\left({\widehat{x}}_k\right)+v_k$$

(32)

Here, $f,h$ are differentiable. In a sensor module, the position data ${\mathrm{z}}_k$ measured at $t$ time contain noise components, which interfere with signal analysis. Therefore, to predict the estimated value ${\widehat{x}}_k$ for which noise components have been removed, the near-optimal Kalman gain $K_k$ is computed to predict ${\widehat{x}}_k$ in the EKT.

$$K_k=P_{k|k-1}{H}_k^T{S}_k^{-1}$$

(33)

$S_k$ is the residual covariance; $S_k=H_k{P}_{k|k-1}{H}_k^T+R_k$. The updated state estimate ${\widehat{x}}_{k|k}$ and covariance estimate can be expressed as

$${\widehat{x}}_{k|k}={\widehat{x}}_{k|k-1}+K_k{\tilde{y}}_k,$$

(34)

$$P_{k|k}=\left(\mathrm{I}-K_k{H}_k\right)P_{k|k-1}.$$

(35)

where ${\tilde{y}}_k$ is the measurement residual; ${\tilde{y}}_k=z_k-h\left({\widehat{x}}_{k|k-1}\right)$.

Here, ${\widehat{x}}_{k|k-1}$ and $P_{k|k-1}$ represent the predicted state estimate and the predicted covariance estimate, respectively, as follows:

$${\widehat{x}}_{k|k-1}=f\left({\widehat{x}}_{k-1|k-1},u_k\right)$$

(36)

$$P_{k|k-1}=F_k{P}_{k-1|k-1}{F}_k^T+Q_k$$

(37)

$F_k$ and $H_k$ are the state transition and state variable matrices defined to be Jacobians. $Q_k$ and $R_k$ are the covariance matrix of the system and the measurement noise, respectively.

The final estimated value ${\widehat{x}}_{k|k}$ is predicted by multiplying an appropriate weight on the basis of $K_k$ with the predicted value ${\widehat{x}}_{k|k-1}$ and the measured value $z_k$, followed by summing the two values, while $K_k$ is inversely proportional to the measured noise component, $R$. Specifically, as the measured noise component becomes greater, $K_k$ decreases, and the rate at which $z_k$ is reflected in ${\widehat{x}}_k$ is decreased. It can be seen that the predicted error covariance, $P_{k|k-1}$, is proportional to the system noise component, $Q_k$, and that $K_k$ is proportional to $P_{k|k-1}$. Therefore, $K_t$ is proportional to $Q_k$. In other words, an increase in $Q_k$ signifies that the reliability of ${\widehat{x}}_t^{\prime }$ is decreased, in which the components of the measured value are considerably reflected when predicting ${\widehat{x}}_{k|k-1}$ [6].

3. Experimental Results by Unit Testing

In this chapter, for evaluating the basic performance of the proposed positioning system, a basic telecommunication test, a LOS (line of sight) test, a non-LOS (NLOS) test, and a stopping test in an anechoic chamber were conducted indoors and outdoors.

3.1. Basic Telemmounication Test

A basic telecommunication test was conducted in various places including indoor, hallway, and outdoor areas, as shown in Figure 5. The experiment examined whether multiple anchors installed in a fixed position of each place effectively monitor the coordinate values of immovable or moving tags. The experimental result showed that the packet with accurate position coordinate values is received normally.

3.2. LOS Test and Non-LOS Test

An LOS test was conducted after verifying the communication function between the anchor and node. Our LOS test installed three anchors in an indoor space without obstacles and we repeatedly measured the coordinates of 16 fixed positions for a tag along the x-, y-axes, as shown in Figure 6. The results of the repeated LOS and NLOS tests are presented in

Table 1. Indoor LOS and NLOS measurement data.

Measurement Error [mm]	LOS		NLOS
	x Axis	y Axis	x Axis	y Axis
Average	68.01	−29.56	−43.25	−45
Minimum	−194	−194	289	107
Maximum	195	188	−382	−237

. All of the errors in indoor spaces were within ±200 mm.

After the LOS test, an NLOS test was conducted to inspect the positioning performance in an indoor space with obstacles. In this experiment, the position coordinates of four tags in the indoor space were repeatedly measured. From Table 1, the measurement errors in the NLOS environment were slightly higher than those in the LOS environment, as predicted. Average measurement errors of LOS and NLOS tests are out 0 mm. However, a lot of tests for various indoors allows the average errors to approach 0 mm.

3.3. Anechoic Chamber Test

Lastly, the measurement test was conducted in an anechoic chamber. As shown in Figure 7, four anchors were installed at a 5 m interval in the length, width, and height directions, and the position of tags were repeatedly measured from four locations in the chamber.

The results of the stopping test in an anechoic chamber are presented in

Table 2. Measurement error for an immovable target.

Reference Values	Reference 1			Reference 2			Reference 3			Reference 4
	x	y	z	x	y	z	x	y	z	x	Y	z
Ref. coordinates (mm)	0	0	1800	0	5000	1800	5000	0	1800	5000	5000	1800
Max. deviation measurement value (mm)	−43	−58	1893	75	4943	1943	4947	−62	1965	4948	4949	1624
Measurement absolute error (mm)	43	58	93	75	57	143	53	62	165	52	51	176

. The tags were installed to be fixed on a tripod. Measurement errors occurred on average 55.75 mm and 57 mm in the x- and y-axes directions, but on average 144.25 mm in the z-axis direction, slightly more than the x- and y-axes direction.

4. Experimental Result of the Position Tracking System and Discussion

Following the measurement testing at a fixed position, the experimental results of a moving object are explained in this chapter. This experiment was conducted in an indoor anechoic chamber and the actual working environment.

4.1. Position Tracking Experiment in Anechoic Chamber

For this experiment, four anchors were installed at a 5 m interval in the length, width, and height directions, as shown in Figure 8. The position was repeatedly measured by controlling the movement of tags in a cubical space, during which the tags were controlled by being installed on a triaxial Cartesian robot. The results of measuring the position of tags moving, as shown in Figure 9, along the triaxial robot, are presented in

Table 3. Measurement data of orthogonal reciprocating motion using testbed.

Axial Direction		Ref. Coordinates (mm)	Mean Value (mm)	Mean Error (mm)
x axis	Min.	2800	2760	40
	Max.	3190	3191	41
y axis	Min.	2710	2730	30
	Max.	3090	3053	37
z axis	Min.	1500	1321	179
	Max.	1890	1985	95

. In this experiment, the position tracking was accurate in the x-axis and y-axis directions, but the positioning error in the z-axis direction was large.

Therefore, a moving average (MA) filter, Kalman filter [7], and our extended Kalman filters (EKT) were applied to the positioning data for improving the accuracy of positioning coordinates. Specifically, the MA, Kalman filter, and extended Kalman filter were applied to the measurement data after placing the tag on the reference coordinate (2700, 2700) on the fixed z axis. The measurements were taken for approximately 380 times to obtain reliable data; the results are shown in Figure 10.

The figure shows the results of applying the MA, Kalman filter [7], and our EKT to the measurement data of moving nodes in an indoor space where multiple anchors are installed, and it can be observed that the position measurement error is significantly reduced when both types of filters are applied. The statistical results of the 380 measurements, as shown in Figure 10, are presented in

Table 4. Measurement data from sensor units and filtering results.

Min/Max	Measurement Value (mm)		Moving Average (mm)		Kalman Filter [7] (mm)		Proposed Extended Kalman Filter (mm)
	x	y	X	y	x	y	x	Y
Min.	2687.5	2673.3	2691.2	2679.2	2693.4	2683.2	2696.6	2686.8
Max.	2722.1	2723.2	2709.6	2710.9	2708.4	2709.4	2707.4	2705.2
|Max.−Min.|	34.6	49.9	18.4	31.7	15.0	26.2	10.8	18.4

; From these results, we found that our ETK-based method had less errors of 7.6 mm and 13.3 mm on the x and y axes compared to the MA method and also had less errors of 4.2 mm and 7.8 mm on the x and y axes compared to the Kalman filter method [7]. According to the numerical data, the application of the extended Kalman filter produced a smaller positioning error, which also can be confirmed through the graphical data.

4.2. Position Tracking Experiment in Anechoic Chamber

Lastly, an experiment was conducted to examine the positioning accuracy for a tool in the actual automobile assembly process at an automobile repair shop, which resembles the real automobile assembly site. The fixed anchors were installed on four posts of the lift, and the tags were attached to the tool to measure its moving path, as shown in Figure 11. For displaying the current position of the tool operating in an automobile model, an application was developed for operation monitoring to display the position of the tool, as shown in Figure 12.

4.3. DAA(Detection and Avoid) Experiment

The UWB central frequency band used practically in the development product is 4492.8 MHz, which is included in the technical criteria of the existing frequency band and the changed frequency band. However, changes in radio frequency and UWB-related laws act as unexpected variables in development. The mean power density, including antenna absolute gain, is −70 dBm/MHz or less. IAT (Interference Avoidance Technology) reduces the amount to less than −70 dBm/MHz within two seconds and avoids it within 2 s if it detects signals from other radio stations greater than −61 dBm during operation.

Spectroscopic analysis shows that the BW (Band Width) was wider than the set 499.2 MHz, as well as signals in addition to BW. If signals were emitted from BW other than 3.735–4.8 GHz BW, the UWB radio frequency band was used by adjusting the register value of DW1000. After BW adjustment, we conducted a wire test and anaerobic radiation test through SMA connector wire. The final test conducted through a spectrum analyzer confirmed that it was within the frequency range in the regulation.

5. Conclusions

In this study, the UWB and MEMS IMU integrated positioning method was proposed for a system for accurately identifying and monitoring the position of work tools which are used in the assembly process of automobile production. Our method measures the accurate position, velocity, and attitude through a MEMS IMU sensor that is derived from the inertial navigation formula in 3 coordinate frames of indoor positioning frame (E-frame), local navigation frame (N frame), EGI frame (I-frame). Then, we computed the distance between an anchor node and a tag node by the TWR method. Here, we used the confidence distance compensation trilateration localization method to prevent the distance error by obstacles or changes in the indoor environment. Our method used the improved EKR by the error dynamic equation derived from 8 × 8 sparse state-space matrix, for the integrated positioning algorithms of MEMS IMU and UWB sensors. We compared our method to the previous UWB and MEMS IMU positioning method-based moving averaging filter and Kalman filter [7]. The results verified that our method has less errors (4.2 mm~13.3 mm) than the previous methods. Furthermore, we verified that the bandwidth is wider than 499.2 MHz and it is within the regularized frequency range, through DAA experiments. Also, we developed to monitor the accurate position of a work tool operating in an automobile model. Recently, the wireless positioning in the IoT has been researched much for indoor location identification [22,23]. Thus, we intend to study the applicability of the wireless positioning in the IoT for the future works.