A Novel Ultra ‐ Wideband Double Difference Indoor Positioning Method with Additional Baseline Constraint

: Ultra ‐ wideband (UWB) technology is suitable for indoor positioning owing to its high reso ‐ lution and penetration. However, the current UWB positioning methods not only fail to fully analyze errors, but do not have the ability to eliminate gross and large random errors. In this article, the errors of UWB indoor positioning are analyzed comprehensively, and the basic function model is given. An indoor positioning method based on a double difference UWB with ranging observations is proposed and realized. In the proposed method, two UWB rover stations and a common base station are intro ‐ duced, and the known baseline length between two rovers is used as the constraint condition for qual ‐ ity control. The observations and coordinate estimations are constrained by the prior and posteriori, respectively, and the weight of ranging observations with large residuals is reduced. Two groups of static experiments are designed. After adopting the proposed method, the plane error of one rover is 3.4 cm and 2.1 cm, and plane error of another rover is 3.3 cm and 2.0 cm, respectively. The positioning precision is improved by more than 80% compared with the traditional method. In the dynamic ex ‐ periment, the coordinates of the starting and ending point obtained by the proposed method are basically consistent with the truth value, and the positioning results are close to the reference tra ‐ jectory. The experimental results show that the proposed method can eliminate systematic and large random errors and improve the positioning precision effectively.


Introduction
Nowadays, global satellite navigation systems (GNSS) are widely used for outdoor positioning and navigation and verified in many high-precision positioning fields. Common GNSS high-precision positioning technologies include real-time kinematic (RTK) and precise point positioning (PPP). They are the representatives of outdoor dynamic high-precision positioning solutions and, nowadays, are fully realized. However, GNSS signal is unreliable in GNSS-denied conditions due to occlusion, such as in indoor environments [1]. Therefore, a variety of indoor positioning technologies are gradually researched. For example, indoor positioning technologies with an accuracy of 3-5 m, such as Wi-Fi, Bluetooth, ZigBee and geomagnetism, can be used in emergency security, intelligent warehousing, precision marketing and mobile health. While audio, UWB, pseudolite and infrared positioning technologies with high-precision positioning accuracy are commonly used in aircraft manufacturing, large facilities monitoring, mine or tunnel engineering geodesy and robot intelligent positioning navigation. Among them, UWB technology is considered to be suitable for highprecision indoor positioning because of its advantages of fast transmission, strong penetration and anti-multipath fading ability [2].
The research into navigation and positioning generally focuses on the anchors arrangement, positioning algorithm and data fusion, error analysis and quality control. Similarly, most of the research on UWB indoor positioning focuses on the anchors arrangement, positioning algorithm, error modeling including non-line of sight (NLOS) and noise.
With regard to the anchors arrangement, the geometric dilution of precision (GDOP), Cramer-Rao lower bound (CRLB) and Root mean square error (RMSE) are usually used as objective functions. Especially the RMSE and GDOP are widely utilized as metrics to evaluate the optimal anchors arrangement [3]. Lin et al. [4] used a multi-lateration method to deal with the anchors arrangement problem, and the minimal position dilution of precision (PDOP) was the goal of optimization. Monica et al. [5] proposed an approach of the optimized placement of anchor nodes using the criterion of minimizing the average RMSE. Chen et al. [6] used the genetic algorithm (GA) to optimize locations of the anchor nodes, and the optimization objective was to minimize the average RMSE. Based on the Chan algorithm, Li et al. [7] analyzed the base station arrangement method in combination with the system construction site. Their optimization objective function was similar to that in the previous literature. Pan et al. [8] considered CRLB and coverage degree criterion as the objective function of the differential evolution (DE) algorithm, simultaneously. Wang et al. [9] proposed a search method for the optimal arrangement of the reference nodes based on a genetic algorithm. Wang et al. [10] proved that the dilution of precision DOP could not be used as the decisive metric for the configuration and optimization of the anchors. In fact, it was difficult to find a ubiquitous arrangement optimization principle due to the multipath and NLOS problems in the indoor environment. The anchors arrangement can only be adjusted according to the actual situation.
In terms of a UWB positioning algorithm, the Chan algorithm and nonlinear least squares estimation are often used for positioning the calculation if no extra nodes are used [2]. Some errors can also be corrected in advance or estimated according to the environmentrelated parameters [11]. These solutions are widely developed and discussed. In order to expand the availability of positioning, other positioning technologies are integrated with UWB technology. Geng et al. [12] fused UWB and ZigBee technologies using maximum likelihood estimation. Graichen et al. [13] applied occupancy grid maps (OGMs) to assist UWB positioning, and the proposed method increased the positioning accuracy. Krapež et al. [14] adopted additional calibration modules (CMs) in a 3D UWB positioning system. The results showed that the calibration before positioning could effectively improve the positioning accuracy. Wang et al. [15] proposed several positioning algorithms for different indoor scenes. The improved UWB location algorithm could meet a variety of positioning requirements. The most common method was to integrate UWB with inertial navigation system (INS). Li et al. [16] introduced map matching to optimize PPP/INS/UWB tight coupling positioning, which could effectively improve the positioning accuracy. Nguyen et al. [17] fused the inertial measurement unit (IMU) and UWB with camera. UWB ranging observations were focused on and the elimination of time offset was studied. The experimental results were clearly better than the traditional methods. In addition, Li et al. [18] proposed applying a particle filter to IMU/UWB positioning. Li et al. [19] used an extended Kalman filter (EKF) to improve the positioning accuracy of UWB and adopted an error-state Kalman filter (ESKF) to fuse UWB and IMU. Guo et al. [20] realized an indoor positioning method based on pedestrian dead reckoning (PDR) and UWB. The average error could reach up to a decimeter. Wang et al. [21] proposed a zero-velocity update (ZUPT)/UWB fusion algorithm based on graph optimization, with an accuracy of 0.4 m. The same idea was also presented by Zhang et al. [22]. Lutz et al. [23] proposed to combine UWB and a visual-inertial navigation system (VINS) to provide accurate pose estimations. Additionally, the drift in VINS could be alleviated. Actually, INS/UWB integration could improve positioning accuracy and availability, and maintain the continuity of positioning results when the UWB signal was missing. However, INS could not be calibrated when the UWB positioning accuracy was not guaranteed. Therefore, the error analysis and quality control of the UWB positioning results was very crucial.
Many studies focused on positioning error analysis and quality control. Bharadwaj et al. [24] studied the effect of human position on line-of-sight (LOS) and NLOS observations of UWB. Cao et al. [25] calibrated the UWB in the LOS scene and used three estimation methods to improve the positioning accuracy in the NLOS scene. Yang et al. [26] introduced the concepts of "virtual inertial point" and "environmental factor" to compensate the NLOS error in UWB indoor positioning. Zhang et al. [27] proposed a phase-difference-of-arrival (PDOA)-assisted UWB positioning method for the NLOS environment. In addition, neural networks were used to classify and recognize LOS/NLOS signals [28][29][30][31]. For example, Cui et al. [32] proposed a LOS/NLOS identification method based on a Morlet wave transform and convolutional neural networks. As for the processing of UWB error (including NLOS), a robust and adaptive Kalman filter was often utilized [33][34][35]. In addition, Xia et al. [36] applied a particle swarm optimization (PSO) algorithm for UWB positioning. The traditional PSO algorithm was used to obtain the initial value, and all the ranging information was applied to refine the positioning result, which effectively reduced the ranging error. However, the signal error was not classified, and for the analyses in these methods, the quality was not been controlled either. In fact, the analysis of UWB positioning errors was of great significance to understand its essence.
In addition to the theoretical research, UWB is widely used for the real-time high-precision positioning in indoor industrial engineering. For example, the status of people and equipment in coal operation are available at any time in the UWB mine personnel positioning system. Early warning and rescue measures can be provided. Furthermore, UWB technology is commercially applied to high-precision indoor and outdoor cable channel inspection. The Local Sense UWB precise positioning system developed by Tsinghua University has a positioning accuracy of 10 cm [2]. It solved the problem of the inability to precisely position components, equipment, vehicles and personnel in the industrial field, such as electric power, tunnels, and in the coal and chemical industry. The UWB system developed by Zebra Company of the United States achieved a sub-meter positioning precision and hwas applied in warehousing and logistics [37]. Ubisense Company in the UK launched the UWBbased transmitting and receiving module. The system measures the arrival time difference and the angle of arrival of the signal, and obtains the three-dimensional coordinates of the target according to the positioning algorithm. The error of the three-dimensional coordinates can be controlled within 15 cm with the Ubisense system [37]. DW1000, launched by the Decawave (now Qorvo) company, is a widely used UWB product [38] which is applied to many UWB-integrated development companies and scientific research institutions in colleges and universities. Wen et al. [39] also put forward a new strategy combining the downlink time difference of arrival (TDOA) based on real-time UWB in tunnels and highway scenarios. While the accuracy only reached 1 decimeter to 1 m. Hasan et al. [40] used UWB to compensate the position of automated guided vehicles (AGV) for the distribution of goods, scheduling of machine processes and coordination. Lawrence et al. [41] proposed using three UWB tags to determine the positions and attitudes of the mobile laser scanner for building information modeling (BIM) 3D data collection.
Although UWB positioning systems in industrial engineering are widely applied, there are still many details that can be further optimized, especially in precision, error analysis and quality control. Current research tends to focus on either error modelling and correction, or error elimination using a robust method. However, it is essential to classify errors to find main factors affecting the positioning. In this article, the main error that affects the accuracy is obtained through the exploration of UWB errors. Combined with the actual ranging observations, it is verified that the ranging errors are mainly composed of systematic and random errors. Regarding GNSS double difference technology, time-independent error components in UWB ranging observations are processed by a differential method, and a positioning algorithm based on double difference (DD) UWB ranging observations is proposed and deduced. Then, we introduce two rover stations to form a known baseline, which is used as an additional condition to constrain the original observations and coordinate the estimation values from the perspective of prior and posteriori. Therefore, the weight of ranging observations with large residual errors is reduced. Finally, through the experiments, the proposed method is verified to eliminate the systematic errors and large random errors effectively. The proposed method has the ability to control the quality of positioning the results and guarantees a precision at the centimeter level. For the convenience of description, the proposed method is named BC-DUWB (baseline constrained double difference UWB positioning).

Two-Way Time of Flight Ranging
The UWB devices used in this article are PulsON410 and 440 (hereinafter referred to as P410 and P440) manufactured by Time Domain Company. Range observations between the anchor node and the tag are obtained using the Two-Way Time of Flight (TW-TOF) ranging approach. Firstly, a positioning tag send a request to an anchor node. Then, the anchor node returns a response to the positioning tag after an observation is obtained. Finally, the tag takes an observation, incorporates the tag observation and a range observation is calculated after the tag receives the response. When combined, the speed of light, time to send the request, time to return the response, time to receive the response, and the centimeter-level distance observation can be calculated. TW-TOF ranging method does not require time synchronization between anchor nodes. The schematic diagram of TW-TOF is shown in Figure 1. The distance P between the anchor node and positioning tag can be expressed by the following formula: where c is the speed of light; Ta1 is the time that the anchor node sends the request to the positioning tag; Ta2 is the time that the anchor node receives the response from the positioning tag; Tb1 is the time that the positioning tag receives the request from the anchor node; and Tb2 is the time that the positioning tag sends the request to the anchor node.

Error Analysis of UWB Ranging
Although the clock drift error is not considered in the TW-TOF ranging method, many errors in UWB ranging derive from a lack of discussion. In this article, UWB ranging errors are classified into the following three categories in Table 1 and analyzed comprehensively. The phase-center deviation is caused by the non-coincidence of the geometric center and the phase center, which vary with the intensity and direction of the signal. Actually, there is an electrical delay between the antenna surface and internal sampler of the electronical device. Additionally, each device used in this paper is pre-programmed with an internal estimate to represent the electrical delay. The electrical delay is time-independent and only related to electronic devices and antenna types. It is also defined as the antenna delay in some articles [23]. The base station coordinate error is small, since the base station coordinate is generally obtained by other high-precision measurement technologies. In this article, a total station is used for the base station coordinate measurement. Total station is an optical instrument commonly used in construction, surveying and civil engineering. The measurement accuracy can reach mm level. For ease of understanding, the coordinate error of the base station is analogous to the orbit error in GNSS. Multipath and NLOS propagation errors are two difficult problems to be solved in positioning, especially in the complex indoor environment. UWB has a strong penetration ability in most cases, and the penetration error is related to the thickness of obstacles, such as wood and walls. In this article, the penetration error is regarded as NLOS error. Furthermore, irrelevant radio equipment should be placed away from anchor nodes and positioning stations during the measurement to avoid the electromagnetic interference.

UWB Positioning Principle
According to the summary of UWB ranging errors, the ranging observations of UWB devices can be expressed by the following function model, which is proposed by referring to the GNSS error model [42]: where k is the current epoch, r is the UWB positioning tag, and s is the UWB anchor node; chor node coordinate deviation may not be the same at different epochs. The truth value of geometric distance can be expressed as: where ( , , ) When ranging observations of the tag to n (n > 3) anchor nodes are measured synchronously, assuming 's rk  is a ranging observation at k epoch without considering the random error, the equations are as follows: where s = 1, 2, …, n. where: The linearized observation equation can be obtained by expanding Taylor series and taking the first term: where Ek is the truncation error and all random errors. Formula (4) can be written as follows: Combined with Formulas (2) and (9), it can be written in the form of matrix: The corrections of the unknown parameters can be obtained as: Given initial coordinate 0 , the coordinates of the positioning tags can be obtained. Initial coordinates are updated by using the coordinate corrections, and satisfactory results can be obtained after 2~3 Gauss-Newton iterations.

Double Difference (DD) UWB Positioning
By analyzing error terms in Formula (2), DD can be considered to eliminate the electronic delay. The antenna phase center deviation and anchor node coordinate deviation can also be reduced. In this article, a stationary tag with known coordinates is called the base station, and a tag with unknown coordinates to be located is called the rover station or rover. In addition, one of the anchor nodes should be selected as the reference node before DD according to the elevation angle or signal-to-noise ratio (SNR). If the ranging observations of Formula (2) are used for single difference (SD), which means the ranging observations between one positioning tag and two anchor nodes are used to calculate the difference simultaneously, Formula (2) can be written as follows:  (12) where r is the rover station tag; k is the current epoch; and Δ is the single difference operator; s = 1, 2.
Through SD of anchor nodes, the electronic delay at the rover station can be eliminated, and the antenna phase center deviation of the rover can be weakened. Similarly, the single difference observations of base station can be written as follows: where b represents base station. The basic expression of UWB double difference observations is as follows: where  represents the double difference operator.
Since the antenna phase center deviation and anchor node coordinates errors are greatly reduced and can be classified as noise, the Formula (14)   where k E  represents the truncation error and all random errors, and Coordinates corrections can be obtained by the least square method, so as to obtain the three-dimensional coordinates of the rover station.

Positioning Quality Control Based on Baseline Constraint
The proposed UWB double difference positioning method is similar to the GNSS RTK principle, but different in the process of error decomposition and difference. Theoretically, time-independent errors in UWB ranging observations can be eliminated or weakened using DD, but the noise is amplified at the same time. Moreover, DD cannot identify gross errors and large random errors. Although there is no systematic deviation in the positioning results, there will still be large outliers. In order to solve this problem and control the quality of positioning results, we introduce the baseline constraint on the basis of UWB double difference positioning, which is defined as BC-DUWB. In BC-DUWB, two or more rover stations are used for DD positioning using a common base station. The baseline length of the two rover stations is known because of the fixed spatial relationship. These known baseline constraints are attached to the positioning process from the two perspectives of priori and posteriori.
The first one is the prior baseline constraint. If two rovers are ranging with a common anchor node (assuming node number is s), the anchor node s and two rovers will form a triangle with known sides. The absolute value of the difference between two rovers and anchor node s should be less than the known baseline length. Although UWB devices have electrical delay and other errors, the systematic error would be close if same types of devices are adopted. Considering the nominal accuracy of UWB ranging, the quality control model of prior baseline constraint can be written as follows: where W s represents the weight of the ranging observations between the anchor node s and the rover station at the current time; 1 2 , s s r r P P represent the ranging observations between two rovers and anchor node s, respectively; lr represents the known baseline length; σ0 represents the nominal ranging accuracy of P440; s r e  represents the difference in systematic errors including electronic delay.
The main purpose of the prior baseline constraint is to identify the larger gross errors in the ranging observations, and set the weight of the observations with gross errors to 0. Figure 2 shows the schematic diagram of prior baseline constraint: The second one is the posterior baseline residual constraint. After obtaining the coordinates of the two rovers, the baseline length lc can be calculated. The difference between the calculated value lc and the truth value lr is defined as the baseline residual. The baseline residual can be used to further optimize the positioning results. If the baseline residual is greater than 3σ, the positioning result is considered unreliable. The anchor node with large ranging errors is found by eliminating the anchor nodes one by one. σ can be obtained through a large number of experiments, and can also be set according to the actual positioning requirements. For example, if the user needs the positioning accuracy to be within 5 cm, the σ can be set to 5 cm. However, this value should not be less than 2.3 cm of p440 nominal positioning accuracy.

Flow Chart of BC-DUWB Algorithm
The whole process of the method proposed in Section 2.2 is represented by a flow chart as follows Figure  Combined with the description of the previous section, the flow of the proposed method is shown in Figure 3. Firstly, the reference anchor node is defined according to the elevation angle after obtaining the observations of the base station and two rover stations. The baseline constraint is then performed using the method proposed in Section 2.2.3. When the reference anchor has no fault, the DD positioning can be carried out combined with the priori baseline constraint. Finally, the posteriori baseline constraint can be used to further improve the precision of positioning results after obtaining the coordinates.

Experimental Design
Seven P410 and P440 UWB ranging devices are selected as anchor nodes. The types of devices numbered 101, 102, 103 is P440 and those of devices numbered 107, 108, 109, 110 are P410. The devices are arranged in a laboratory of about 14 m × 14 m × 4 m (length, width and height). The identifier of the base station is 104 (P440) and the identifier of the rover station is 105 and 106 (P410). Many parameters of P410 and P440 are similar, and some of the main parameters are listed in Table 2. The location of the anchor nodes is shown in Figures 4 and 5. Positioning tags should not be placed on the same plane as the anchor nodes. It is generally considered that the elevation angle from the positioning tag to the anchor node is greater than 30 degrees [43]. In addition, the minimum GDOP principle is also used to place anchor nodes [3], but whether this method is suitable for UWB remains to be discussed. Thus, the principle of UWB anchors distribution in this article mainly refers to the PulsON440 user guide [43].  In addition to the ranging experiment in Section 3.2, two groups of static and one group of dynamic positioning experiments are carried out, named Test 1, Test 2 and Test 3. In static experiments, the rover and base station are all fixed and truth values are measured by total station. In the dynamic experiment, two rover stations are placed on the trolley. The trajectory of the trolley is approximately a rectangle. Truth values of the starting point, the corners of the rectangle and the end point are measured by the total station. However, since the corner of the rectangle may not be exactly reached when turning, the trajectory around the corner becomes a curve.

Error Analysis of UWB Ranging Observations
The purpose of the ranging experiment is to verify that UWB ranging errors are composed of time-independent systematic errors, gross errors and random errors. In this experiment, two UWB anchor nodes are randomly selected for ranging, and the total station is used to obtain the truth distance. Figure 6 shows the ranging observations between the anchor node #104 and #107/#108, respectively. Table 3 lists the ranging observations, standard deviations and truth values.  According to Figure 6 and Table 3, it can be considered that ranging error is composed of systematic, gross and random errors. The gross error is generally more than three times of the mean square error, the random error is small, and the system deviation is about 10 cm. In addition, the ranging observations of P440 and P410 are greater than the truth values, and systematic differences between different nodes are relatively consistent. This provides a favorable condition for the identification of large gross errors with prior baseline constraints.

Static Experiments
In the static experiments, Test 1 and Test 2, two rover stations are placed somewhere in the laboratory, and the base station is fixed in a place with a good observation condition where the ranging signals of seven anchor nodes can be received at the same time. The true locations of the base stations and rover stations in Test 1 and Test 2 are shown in Figure 7.  In order to simulate the positioning effect under different observation conditions, we walk around the rover stations to add interference and occlusion in Test 1. Instead, the interference and occlusion are not allowed in Test 2. It can be considered that the observation condition of Test 1 is worse than that of Test 2.

Test 1 with Interference and Occlusion
Referring to the description in Section 2, when the baseline constraint is carried out, the difference in ranging observations between an anchor node and two rover stations needs to be calculated. Taking two anchor nodes (#109 and #110) as examples, the difference and the known baseline length are shown in Figure 8. According to Figure 8a, the ranging observation difference between the two rover stations and the anchor node #109 is less than the known baseline length, which indicates that there is no gross errors in the observations of #109. However, there are a few outliers in the figure which cannot be identified by the prior baseline constraint but can be eliminated by the posterior baseline constraint. As for Figure 8b, there are many outliers in the ranging observations difference between the two rover stations and the anchor node #110. The observations difference in many epochs is much larger than the known baseline length, which indicates that there are a lot of gross errors in #110 observations. Therefore, the weight of the observations of #110 is set to 0. For the observations of the other anchor nodes, the same way is used for the prior quality control. Finally, to ensure the positioning precision, epochs with a large random error are further eliminated through the posterior baseline residual constraint. The positioning results of the common method, DD UWB positioning method without additional constraints (DUWB), and DD UWB positioning method with additional baseline constraints (BC-DUWB), are plotted in Figure 9. According to Figure 9a,b, the positioning results of the common method contain gross errors and also have obvious systematic deviations from the truth value. The systematic error of the positioning results is eliminated in DUWB, but the influence of gross error is further amplified. Since UWB signals cannot penetrate the human body, the observations in some directions will be lost due to human interference. This makes the positioning results present a great deviation. Different from the previous two methods, the positioning results are very close to the truth value in BC-DUWB, since systematic and large gross errors are eliminated, and only some small random errors exist. Tables 4 and 5 verify this conclusion through STD and RMSE. When the data are of poor quality, the positioning precision is low because DUWB cannot deal with gross errors. In Test 1, the RMSEs of the two rover stations calculated by BC-DUWB are 3.4 cm and 3.3 cm, respectively. The positioning precision of the proposed method is improved by about 81% compared with the common method. The STDs of the two rovers are 2.5 cm and 2.7 cm, respectively, which indicates that the positioning results are stable and the accuracy is high.

Test 2 without Interference and Occlusion
Different from Test 1, the observation environment of Test 2 is better as it experiences no interference. Taking two anchor nodes (#107 and #108) as examples, the known baseline length of the two rover stations, as well as the ranging observation differences between #107 and the two rover stations are shown in Figure 10a. The similar results of #108 are shown in Figure 10b. According to Figure 10a,b, there are a few outliers in the ranging observations difference between the two rovers and the anchor nodes #107/#108. The observation differences in some epochs are much larger than the known baseline length, which indicates that there are a few gross errors in the observations of #107 and #108. In these epochs, the observations of #107 and #108 should be reduced to 0. Similarly, for the observations of other anchor nodes, the same way is used for the prior quality control. In order to ensure the positioning precision, large random errors are further eliminated through the posterior baseline residual. For comparison, the positioning results of the common method, DUWB positioning method and BC-DUWB positioning method are plotted in Figure 11. According to Figure 11a, the positioning results of rover #105 using the common method only contains a small amount of gross error due to the good quality of observation data in Test 2. However, there is still an obvious systematic deviation from the truth value. Similarly, the systematic error is eliminated in DUWB, but the gross error is still unsolved. The positioning results are very close to the truth value in BC-DUWB, since the systematic and large gross errors are all eliminated, and only some small random errors exist. As for Figure 11b, due to the original observations that rover #106 does not contain gross errors, the positioning results obtained by the common method only contain the systematic deviation. In this case, the positioning results in DUWB and BC-DUWB are basically the same, since the systematic error can be eliminated in both methods. Tables 6 and 7 verify this conclusion from STD and RMSE. When the quality of observations is good, the STD of the positioning results of all methods is small and the accuracy improvement is not clear. The reason for this is that there are a few gross errors in the good quality data and the effect of the additional baseline constraint is not clear. In this test, the RMSEs of two rovers calculated by BC-DUWB are 2.1 cm and 2.0 cm, respectively; the accuracy is higher than that of the common method by 84% and 81%, respectively. The standard deviations of the positioning results of the two rovers are 1.5 cm and 1.7 cm, respectively, slightly worse than for common methods due to DD amplifying the noise.  Test 1 and Test 2 show that the baseline constraint has little significance when the quality of observations is good. However, when the data quality is poor, the positioning results deviation will become larger by using DD. Therefore, the baseline constraint is necessary to ensure the positioning precision.

Dynamic Experiments
The dynamic experiment, Test 3, is carried out to verify the dynamic positioning effect of the proposed method. In Test 3, two rovers are placed on the trolley and the trajectory of the trolley is roughly a rectangle. Since it is difficult to obtain the truth value in the process of dynamic positioning, the coordinates of each corner of the rectangle are measured by the total station before the experiment. Furthermore, the coordinates of starting and ending point are measured as well. When the trolley is pushed to move, the actual trajectory at the corner appears as a curve since it may not arrive exactly at the corner during turning. The location of the two rover stations and the trolley in Test 3 is shown in Figure 12.  According to Figure 13, there is a significant deviation between the results and the reference trajectory when using the common method. Compared with #106, more gross errors appear in the observations of rover station #105, which are most likely due to the greater impact of occlusion on the movement. As a result, reliable positioning results cannot be obtained in all epochs. After applying DD, although the systematic deviation is eliminated in DUWB, there are still many gross errors, especially in rover station #105. The results of BC-DUWB proposed in this paper are the closest to the reference trajectory. Furthermore, the positioning results of the starting point and ending point are completely consistent with the truth values. In addition, the results of BC-DUWB do not contain large random errors. In order to further compare the differences of the positioning results, the positioning results of #105 and #106 obtained by these three methods are plotted in Figure  14a and Figure 14b, respectively. The GDOP and the number of available anchors during the dynamic test are shown in Figure 14c,d.   Figure 14c,d shows that when the number of available anchors decreases, GDOP increases rapidly. In addition, small GDOP during the experiment indicates that the current distribution of UWB anchors does not affect the positioning precision.

Conclusions and Discussion
When UWB technology is used for indoor ranging and positioning, the three-dimensional coordinates of the rover stations are generally taken as the parameters to be estimated. The Chan algorithm and nonlinear least squares estimation are often used for positioning calculations. In this article, the UWB error is analyzed and an UWB indoor double difference positioning algorithm with a baseline constraint is proposed. The main contents and contributions are as follows: (1) The error is classified and verified, and the function model is established.
Current research focuses more on error modeling and elimination, but lacks comprehensive error decomposition and analysis. Based on PulsON UWB devices, this article focuses on exploring UWB ranging errors, and proposes dividing the errors into three parts: anchor nodes, tags and propagation-path related errors. The corresponding function model is also given in this article. Furthermore, UWB ranging errors are classified into systematic errors represented by electrical delay, antenna phase center deviation and anchor node coordinate error, and random errors and gross errors are represented by noise and multipath.
(2) A double difference method is proposed to eliminate some systematic errors.
As for systematic deviation, a double difference positioning method based on UWB TW-TOF ranging observations is proposed and the formulas are derived to prove the feasibility of eliminating this systematic deviation theoretically.
(3) A strategy based on known baseline length constraint is proposed to identify and eliminate large random errors.
For large random and gross errors, two rover stations are introduced and the DD UWB positioning method with additional baseline constraints is proposed for positioning quality control.
The experimental results show that the proposed method can eliminate the systematic error effectively and reach the high accuracy of the positioning results. After adopting (1) The positioning precision of two rover stations in static experiments is more than 80% higher than that of the common method. When the data quality is poor, the positioning precision can reach about 3.4 cm. The positioning precision can reach 2.0 cm when the data quality is good, which further illustrates the robustness of the proposed method. (2) In the dynamic experiment, the results of BC-DUWB are not only closer to the reference trajectory than adopting the common method, but the coordinates of the starting point and ending point are almost consistent with the truth value. Furthermore, the gross error and large random error are also eliminated, so that the maximum deviation is less than 5 cm.
However, there are still some deficiencies in this article, such as the fact that the principle of anchor nodes arrangement has not been studied. Therefore, various arrangement methods will be compared in the near future. In addition, the trolley velocity in the dynamic experiment is slow, and lacks tests for fast-moving scenes. The dynamic analysis model needs to be improved to adapt to the fast-moving application scenarios. UWB has a strong penetration ability in some cases, such as in wood and stones, and the penetration error can be corrected by modeling. However, it has a weak ability to penetrate water and metal. Therefore, of the expected errors in this paper, other UWB errors including penetration error and NLOS will be the next research focus.These factors will affect the validity of the results. Since the main purpose of this paper is to verify the feasibility of the method, there is no modeling of penetration errors or NLOS in the algorithm. For the applicability of the algorithm, the algorithm should be optimized to adapt to the NLOS environment in the future. Furthermore, the anchor nodes will be placed in different rooms and conduct networking experiments in practical environments, such as large scene factories and underground projects.
Author Contributions: All authors contributed to the manuscript and discussed the results; Yinzhi Zhao developed the original idea and designed the software; Lixian Cai conceived the experiments and finished the data collection; Jingui Zou performed data processing and contributed to the final revision of the paper; Jiming Guo provided critical comments and performed data acquisition; Gege Huang assisted in data analysis and contributed to the data expression. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data that support the findings of this study are available from the first and corresponding author, upon reasonable request.