UWB Base Station Cluster Localization for Unmanned Ground Vehicle Guidance

In this paper, we seek to provide unmanned ground vehicles with positioning service using ultrawideband (UWB) technology, a high-accuracy positioning approach. UWB is chosen for two distinct reasons. First, it does not rely on global navigation satellite systems like GPS, making it able to be applied indoors or in an environment where GPS signal is unstable. Second, it is immune to interference from other signals and accurate enough to guide unmanned ground vehicles moving precisely in a complex environment within a narrow road. In this paper, three UWB base stations are aggregated as a group in a 2D space for localization. A large number of tests are performed with a UWB base station cluster in order to validate its positioning performance. Based on the experiment results, we further develop a dynamic particle swarm optimization-based algorithm and a genetic algorithm to deploy multiple clusters of UWB base stations to cover an area of interest. (e performance of the proposed algorithms has been tested through a series of simulations. Finally, experiments using unmanned ground vehicles are carried out to validate the localization performance. (e results confirm that the robots can follow complex paths accurately with the proposed UWB-based positioning system.


Introduction
Accurate positioning is essential for robot control while challenging to achieve, especially in an unstructured environment. Many applications, such as multirobot formation control and robot path planning, strongly rely on positioning accuracy either indoors or outdoors. Depending on the environment, positioning technologies can be roughly classified into outdoor localization, such as GPS, and indoor positionings, among which we can find Bluetooth, WiFi, Zigbee, and UWB technologies. As UWB-based localization can achieve higher accuracy positioning compared with other techniques, it is chosen in this paper to guide unmanned ground vehicles (UGVs).
UWB has been a particularly attractive localization technology, especially for indoor positioning. Normally, UWB base stations, also known as beacons or anchor nodes, are scattered in the area of interest. en, a UWB tag or agent needs to communicate with several base stations to obtain its position [1,2]. One advantage of this scheme is an area of interest can be covered with fewer base stations, reducing the cost. Its drawback is that it is challenging to achieve time synchronization among anchor nodes. us, complex communication protocols need to be designed to eliminate clock offset among anchors, and multiple communication packages are required between tags and anchors for distance measurements or for localization. As a result, this strategy works well for fewer tags but is unable to provide real-time positioning for a tag swarm. Besides, it is also hard to obtain accurate positions of anchors during deployment or installation when it has a long distance between two anchors, which can introduce positioning errors for the entire system.
In response to these problems, we apply UWB beacon clusters for localization. To be specific, a group of base stations can be regarded as an independent positioning device, which is able to provide highly accurate position information for tags within its coverage range. is is very useful for swarm robotic control with the leader-follower structure [3][4][5] because one cluster can be equipped on the leader robot, providing highly accurate positioning information for its followers. Since the beacons in one cluster are near each other, the installation accuracy can be guaranteed. Furthermore, it is feasible to share the same clock source among beacons in the same group, achieving clock synchronization. In this way, it is promising to realize accurate positioning with less communication, so that more tags can be supported simultaneously.
In order to demonstrate our idea, a relatively simple situation is first considered in this paper, where three base stations are combined as a group to provide 2D petitioning for tags. Plenty of experiments are carried out to test the performance of the UWB cluster to determine the proper distances among anchor nodes. ree anchors are installed as equilateral triangles with different side lengths. In order to guarantee the system performance and improve efficiency, we develop a test platform that can adjust the side length continuously and an app used to record the positions of tags automatically. Based on the experiment results, an optimal distance between a pair of base stations and the range a beacon group can cover are obtained. Secondly, two strategies are explored to deploy UWB clusters to cover an area of interest. e main challenges for these strategies are the irregularity of UWB cluster coverage, the randomness of the shape of the area of interest, and uncertainty in the number of UWB clusters. us, we develop two algorithms to determine the best solution to deploy beacon clusters, namely, a dynamic particle swarm optimization-(DPSO-) based method and a genetic algorithm-(GA-) based method. Finally, the experiments of UGVs are performed based on UWB positioning to validate the performance of the localization system. e results show that UGVs can follow multiple complex paths. e maximum error is within 30 cm, and the average error is around 6.7 cm. Furthermore, one beacon group can simultaneously support multiple UGVs. e main contributions of this paper are as follows: (i) An idea is proposed to combine several UWB base stations as one positioning cluster, which can potentially support more tags simultaneously, laying the ground for future application in swarm robotics (ii) Experiments are carried out to test the positioning performance of a UWB group, concluding the tendency of positioning accuracy with the coverage range and the change of distances between base stations, which can be used as reference guidelines for the future application of UWB clusters (iii) Two algorithms are proposed to deploy UWB clusters to cover an area of interest, addressing the issues of the irregularity of UWB cluster coverage, the randomness of the shape of the area of interest, and uncertainty in the number of UWB clusters e remaining of the paper is organized as follows. In Section 2, current works related to positioning for robots are presented. In Section 3, the developed test platform is introduced and the experiment results are provided. In Section 4, two proposed algorithms for deploying a set of UWB groups into a random area of interest are explained. In Section 5, we introduce the design of a UGV prototype and experiment results of UGVs following different trajectories under the positioning of a UWB cluster. Finally, in Section 6, concluding remarks are provided.

Related Work
Positioning plays an essential role in robot control. Various positioning methods have been proposed so far, with different features and suitable applications. e most widely used positioning technique for robots in an outdoor environment is GPS [6,7]. It is low in cost and applicable to various seniors. Two main drawbacks constrain its further applications, i.e., its dependency on satellites and its low accuracy, with positioning error reaching meters. ese two drawbacks encourage the study and application of other positioning techniques on robots.
Inertial positioning has a long history, and lots of great achievements have been made so far in this field. is approach estimates positions according to acceleration and is widely used in many applications such as underwater positioning since it does not rely on any external information. However, the integral of acceleration introduces cumulative errors and finally affects the estimation accuracy. So far, various hybrid positioning methods have been proposed by combining inertial positioning with other methods, such as GPS [17,18], acoustic ranging [19], and Doppler speedometer [20]. e point cloud-based positioning includes a family of methods using kinds of sensors to detect environments around and represent sensing data as a point cloud. From the point cloud, the environment around a robot can be reconstructed. By comparing the environment with a digital map, the robot can recognize its position. Point clouds can be generated using various sensors suitable for different environments. Laser radar [11] can reach an accuracy of millimeters but is expensive. Although millimeter-wave radar costs less, it correspondingly provides few data points in space [12]. An infrared (IR) array can be equipped around a robot to provide a set of range information [21]. e provided data points are so rare that they can only be used for positioning in a structured environment. Ultrasonic sensors are used in a similar way as an IR array. Due to the low cost, they are equipped with various robots, for example, the AmigoBot [22]. From the perspective that positioning is achieved via comparing a reconstructed digital environment with a preinput map, vision-based localization can also be regarded as a special point cloud-based positioning, treating pixels as data points in the cloud. It is a rapidly developed field and various schemes have been proposed, such as positioning by recognizing landmarks [23][24][25] or detecting motion by comparing successive images [26]. Novel algorithms are mainly based on deep learning, enabling robots to construct a 3D map with only one camera [13]. However, vision-based positioning heavily relies on a preconstructed map, making it hard to fit for dynamic environments. Furthermore, it is prone to be deceived by similar scenes.
Indoor positioning includes a family of positioning technologies suitable for GPS-denied areas, which include WiFi, Bluetooth, FRID, magnetic field, and UWB. Usually, these methods rely on beacons or base stations. For example, the WiFi-based positioning requires to predeploy a set of APs in the environment [27], and the magnetic field-based positioning method requires to preembed magnetic nails on the floor [28,29]. At present, WiFi-based and Bluetoothbased methods mainly use the attenuation of signals to estimate distance, probably resulting in high positioning errors reaching 10 meters. e positioning accuracy may be enough for guiding customers in malls, but not enough for robot guidance. FRID-based positioning is realized by recognizing a tag when the tag is near to a reader. However, the distance between them needs to be close enough, constraining its application. Magnetic field-based positioning is precise enough to guide robots, making it suitable for automatic parking and automatic storage [30,31]. However, the deployment of magnetic nails needs to be completed during the construction of the building. UWB-based positioning is achieved by measuring the distances between a base station and a tag using the time of flight, two-way ranging, or time of arrival, etc. e centimeter-level accuracy can be achieved by using UWB technology [32][33][34].
In this paper, we seek for a positioning method suitable for robots moving in a complex area.
is requires high positioning accuracy and working with fewer constraints. ese requirements leave us limited choices. GPS, WiFi, Bluetooth, and RFID are not suitable due to large positioning errors. Furthermore, GPS is prone to fail indoors. e point cloud-based localization schemes such as IR, laser, and ultrasonic either are expensive or require pregiven maps for localization. Image-based positioning methods suffer from the same problem. As a result, UWB positioning is chosen as the solution considering the balance among performance, environment adaptation, and cost. Besides high positioning accuracy, UWB technology also has the advantages of high transmission rate, low power consumption, and strong anti-interference, making it suitable for guiding robots. In this paper, the positioning is realized based on ranging between base stations and a tag using twoway ranging. For more details about current research of UWB-based positioning, refer to [33] and other related literature studies.

Performance of UWB Station Cluster
As mentioned in Section 1, three UWB base stations are combined as one group in 2D space, providing positioning information for tags within the group's coverage region, i.e., the overlap coverage region of the three base stations. To avoid bias caused by asymmetry, three base stations are initially placed as an equilateral triangle. Besides, we try to make side lengths of a triangle small so that the UWB cluster can be made as one small device. erefore, the effect of side lengths on the localization accuracy needs to be determined by checking the positioning accuracy of tags in different locations.

Test Settings.
e test platform is shown in Figure 1, which consists of two tracks with an angle of 60°. e sliders can move continuously along guide rails. TT DC motor is installed on each slider, which is reserved for automatic testing in the future. One UWB beacon is installed at the intersection of the extension of two rails, which is defined as the origin (0, 0) in the Cartesian coordinate system of the UWB cluster. e other two UWB beacons are mounted on top of the two sliders.
In order to facilitate the testing process including data recording and data analysis, a Windows app is developed as shown in Figure 2. e app can record UWB tag distances and positions, show tag trajectories, and control the movement of a UGV. It can also be used to update positions of UWB base stations.

Test
Results. In this part, UWB-based ranging and positioning tests are carried out to determine the proper side edges of one UWB cluster. Due to the limitations of the physical environments, the maximum distance between a UWB tag and a base station is around 8 m in a room and around 38 m in a lobby. In order to achieve the best performance, the tests are conducted under line-of-sight propagation conditions. e absolute positions of sampling points are depicted in Figures 3 and 4 , respectively. For each point, more than 50 samples are collected. e mean distance errors between a UWB tag and base stations are shown in Figures 5 and 6 , respectively. ese experiments are carried out in a room and in a lobby in order to determine the ranging performance of UWB. In a room, the mean distance errors are less than 20 centimeters in most cases. Within 40 meters, the maximum mean distance error is less than 40 centimeters. e distance errors grow with the increase of the distance between a station and a tag, indicating by Figures 5 and 6. is implies that the positioning error will be large when a tag is far from the three base stations.
Meanwhile, mean positioning errors are measured in the same environments. e corresponding results are depicted in Figures 7 and 8 , respectively. At present, two side lengths, i.e., 50 cm and 60 cm, are verified. From Figure 7, it can be seen that the performance with two side lengths is similar.
en, more experiments are conducted in a lobby with long distance, whose results are shown in Figure 8. e mean position error increases approximately in proportion to the distance from the origin in a long-term perspective. e mean position errors are no more than 60 centimeters in most cases within a 10-meter distance. erefore, 50 cm side length is chosen for a UWB cluster and a 10-meter distance is determined as the coverage radius of a UWB base station.

Region Coverage with Multiple UWB Clusters
Based on the experiment results in Section 3, both the proper side lengths of the triangle composed of three UWB base stations and the range covered by each UWB cluster can be determined.
In general, for a complex area of interest, one group of UWB base stations is unable to cover the whole area. en, we need to find the best solutions to deploy base station groups. ere are three main challenges: (1) e region covered by all the three base stations are the overlap area of three discs, which is irregular, making it hard to determine if a subarea is covered by the group (2) e area of interest is of random shape with inaccessible region (3) Because of (2), the number of UWB groups necessary to cover the area of interest cannot be determined in advance To solve this problem, we first approximate the overlap region of the three base stations as a disc, as explained in Section 4.1. en in Sections 4.2 and 4.3, two algorithms are proposed to deploy UWB base station groups to cover a randomly generated area. Finally, in Section 4.4, the performance of the two algorithms is compared.

Approximate the Covered Region as a Circle.
In this paper, UWB positioning is achieved through ranging between base stations and tags. In 2D space, if one robot wants to determine its position, one option is to obtain distance measurements to at least three base stations that are not in a line. In our case, three base stations are located near to each other, forming a group. To determine positions with the group, a robot needs to be within the coverage range of each of them, i.e., it should be within the overlap region of the coverage range of all the three base stations, as shown in Figure 9. We call this overlap region S and obviously it is irregular. To simplify the problem, S needs to be approximated into a regular shape, and there are two choices. We can either use an inscribed triangle S t or a circumcircle S c to approximate S, and the one whose area is closest to that of S is the solution. us, it is required to calculate the area of S, S t , and S c . To simplify the expression, in the following, S, S t , and S c are used to represent the area of the corresponding region.  In our case, three base stations are located at (x i , y i ), i � 1, 2, 3, composing an equilateral triangle. As all the three base stations are identical in hardware configuration, the coverage range is supposed to be the same, which is a disc with a radius of R. With the information, S can be determined. e three vertices of S are also the vertices of S t , which can be obtained with the following procedure.
Define the vertices as x i,j k , where i, j, k ∈ 1, 2, 3 { } and i, j, k are different, representing the vertex who is a cross point of disks centered by (x i , y i ) and (x j , y j ), and this cross point is within the range of the third disk centered by (x k , y k ). en, we have (1) With the expressions above, we can determine the three vertices. Furthermore, we can obtain any edge length of S t by computing the distance between the corresponding two vertices. For example, As the center of discs is placed as an equilateral triangle, thus the three vertices are equivalent, making S t also an equivalent triangle, whose edges have us, For S c , its center can be obtained with e radius of S c can be measured as the distance between (x c , y c ) and a random vertex (x So the area of S c can be obtained as follows: e overlap region S is a bit tricky, as shown in Figure 10. It consists of the inscribed triangle S t and three arch form areas, represented as S k arch , where k indicates the center of the circle corresponding to the arch and k � 1, 2, 3. en, e area of S k arch can be obtained with where S k fan is the corresponding fan centered by (x k , y k ) and S triangle is the triangle with three vertices of (x k , y k ), (x k,j i , y k,j i ), and(x k,i j , y k,i j ), as shown in Figure 10. To calculate the area of S fan , it requires its radius R and its span angle θ. R is obtained by experiments, and θ is lacking.   The side length is 50 cm The side length is 60 cm To calculate θ, we first get two vectors starting from the center of the disk (x k , y k ) and pointing at the other two points as Since e area of the fan can be obtained as where S triangle is an isosceles triangle with two edges being R and the base line is the chord of the arch, whose length is en the area of the triangle can be calculated with Heron's formula. Half of the perimeter is given by en the area is

Mathematical Problems in Engineering
With S k fan and S k triangle , As the three arch form areas are equivalent, the area of S then is where k could be 1, 2, or 3. Compared with S, S t is a bit smaller while S c is a bit larger. Define the rate of approximate error as e change of approximation error is shown in Figure 11. If the distances among base stations are small, approximate the overlap region as a disc is a better choice. In our case, according to Section 3, three base stations are placed as an equilateral triangle, whose side length is 50 cm. e error rate is 1.9023%. is means the disk can approximate the actual overlap region with a minor error. en methods to place multiple such disks to cover an area of a random shape can be considered in the following.

Deployment with Dynamic Particle Swarm Optimization.
We intend to place multiple groups of UWB base stations to provide positioning services for robots. In Section 4.1, it is shown that the irregular overlap region of the three base stations can be approximated as a disk reasonably. en the problem becomes where to place these disks and how many such disks are needed. erefore, the deployment scheme is a trade-off between coverage rate and cost.
Inspired by [35], a PSO-based algorithm is proposed in this paper that can dynamically adjust the number of UWB clusters. e algorithm is called dynamic PSO, or short as DPSO.
Generally, DPSO consists of two loops. e outer loop adjusts the number of station groups, while the inner loop adjusts the positions of station groups following the idea of PSO. Assume that in the i out iteration of the outer loop, N out i station groups are used. en, for its inner loops, the particles are encoded as the following sequence: where (x j , y j ) ∈ S and (x j , y j ), j � 1, 2, . . . , N out i , indicates the position of the j th station group. e fitness function is set to maximize the coverage rate where S P is the valid area covered by the station groups and S is the area of the valid area in the map. "Valid" refers to only the region accessible for the robots counts. To measure S P and S in an irregular area, we mesh the map into a set of grids and estimate areas by counting the number of grids. A swarm of particles is generated randomly in the init process. en the particles search in the solution space guided by some historical best solutions randomly, as introduced in [36]. Besides, the particle P should not move out of the valid area S, so that a filter is designed. According to the idea above, the pseudocode of DPSO is presented in Algorithm 1. Key functions of Algorithm 1 are detailed in Algorithms 2 and 3, respectively. A general flowchart of Algorithm 1 is shown in Figure 12.
us, Algorithm 1 has a complexity of In Algorithm 1, the inputs are the valid area S and the radius R of the approximate coverage disk. e output is a set S best , consisting of best particles with different numbers of base station clusters. e parameters for the algorithm are set in lines 1-2. In line 4, the range of numbers of clusters to be explored is calculated. e floor is the number of groups needed to stuff the area with disks, which is S/(πR 2 ), and the ceil is the number of groups needed to stuff the area with the inscribed square of the disk, which is S/( � 2 √ R) 2 . e loop in lines 5-20 is the outer loop, which increases the number of groups by one each time, expending the length of the particle. In each iteration, the particle swarm P, the velocity set V, the historical best solution for each particle p best , and the best solution for all particles g best are reset. en, within the inner loop in lines 10-18, the velocity and position of each particle are updated randomly, affected by g best and p best . en g best and p best are updated.
In line 6, we use function GeneratePos to obtain a random particle, as in Algorithm 2. e challenge here is that S is of random shape, and there is a high chance that the particle obtained with a general random number generation algorithm will go out of range. To solve this problem, in GeneratePos, we first scatter S into a set of grids and then stack the grids up into a sequence. e grids in the sequence are then labeled successively. We generate i out random numbers within the range of [1, N s ], where N s is the number of the grids. ese random numbers are the indexes of grids chosen as elements in the particle. e complexity of Algorithm 2 is O(N p ).
In line 7, the function GenerateVel() initializes the velocity of a particle with a zero sequence having the same size of the particle.
In line 14, we use function Constrain to ensure particles will not move out of S, as in Algorithm 3. According to the encoding rule, the particle represents the coordinates of the station groups. If a station group moves out of S, we just ignore the update of its position and corresponding velocity this time. Meanwhile, the updates of velocities and positions for other station groups are not affected. Algorithm 3 has a complexity of O(N g ).
In lines 9, 15, and 17, the fitness function f(·) needs to be calculated to choose the best particle. As S is irregular, a meshing scheme is applied to estimate the coverage rate: 8 Mathematical Problems in Engineering where N covered is the number of grids covered by at least one of the station groups and N S is the number of grids in S. For each grid, we calculate the distance from its center to all the center of station groups as dis � d i |i � 1, . . . , i out . If min(dis) < R, we treat the grid being covered. In this way, N covered can be obtained. Running this algorithm, we notice that there is a chance that the best particle with more station groups performs worse than that with fewer station groups, as shown in Figure 13. is is because in line 6 of Algorithm 1, each time the number of station groups is changed, the particle swarm is regenerated randomly. e red dots in Figure 13(a) indicate that, with 48 stations, the coverage rate is higher than that of 49 and 50 station groups. Clearly, the solution in Figure 13(b) performs better, with fewer station groups.
ere is a chance that the init swarm is rather bad, and we do not iterate enough times in the inner loop to obtain a good solution. It is not easy to check and decide when to increase the iteration number for the program automatically. With DPSO, the tendency is that more station groups will increase coverage rate, but not stable.
To solve this problem without increasing the iteration number, the DPSO algorithm is improved into the Elite DPSO algorithm, which ensures that the more UWB clusters are used, the better solution can be obtained. Elite DPSO is almost the same as DPSO, but line 6 in Algorithm 1 is changed as P ≔ P random ∪ P elite , P random ≔ p k |p k ←GeneratePos, k � 2, . . . , N swarm , is means that each time the group number is increased, N swarm − 1 particles are randomly regenerated, as is done in DPSO. A special particle called elite particle is also generated with function GenerateElite(g best ) based on the best particle g best obtained in the last iteration of the outer loop. It should perform better than g best . Two kinds of GenerateElite(g best ) are explored.
One way to generate elite particle is to add a cluster with a random position in S to g best . In the worst case, the newly added cluster overlaps with another cluster in g best , resulting in f(p elite ) � f(g best ). Otherwise, p elite covers more region, ensuring that f(p elite ) ≥ f(g best ). is scheme is named as Random Elite DPSO.
Another way to generate elite particle is to add a cluster far from other clusters in g best , encouraging station groups to move to a relatively empty region. Assume that the newly added cluster is located at (x, y) and other clusters in g best are located at (x i , y i ), i � 1, . . . , i out . en, the optimal position can be obtained with

Error of S c
Error of S t Figure 11: Approximation errors of S c and S t with the increase of distances among base stations. When the distances among stations are less than 1 m, the approximation error of S c is less than 5%, while that of S t reaches 60%.
Apparently, d(x, y) is a convex function, so its maximum value is obtained with (x, y) on the border of the region. All grids at the border of S are exhausted to obtain (x, y). is method is named Sparse Elite DPSO because it intends to encourage the locations of station groups to be sparse. e performance of both schemes is compared, showing in Figure 14. From the figures, the Sparse Elite DPSO method performs slightly better than the Random Elite DPSO method. However, considering the computation cost on average, the Random Elite DPSO is a more reasonable choice.

Deployment with Genetic Algorithm.
e DPSO methods proposed in Section 4.2 explore a set of group numbers and provide the best solutions correspondingly. However, they are unable to determine the best group number automatically. A GA-based algorithm is then developed, which can provide the best number of clusters needed to cover a region and the deployment scheme, indicating positions of clusters.
We scatter the area of interest with random shape into N S grids and use the center of the grids to represent them.
us, these grids are labeled as g 1 , g 2 , . . . , g N S . en in GA, the chromosome is encoded as an element vector with the length of N S . Each element is either 0 or 1. If the i th element is 1, then a cluster is deployed at the center of grid g i . Otherwise, a cluster will not be deployed. is implies that where S g1 indicates the sum of the area only covered by one cluster, N covered1 is the number of grids whose center is covered only by one cluster, and N S stands for the total number of grids. e ratio of N covered1 and N S is used to estimate the fitness function value. Defining the fitness function in this way encourages clusters to cover more grids, which satisfies the target that the whole area of interest to be covered. Meanwhile, this definition can prevent the algorithm from adding too many clusters to achieve a high coverage rate because the increase of overlap rate results in a decrease of fitness value. In this way, the algorithm can adjust the number of clusters automatically. en a swarm of chromosomes is generated. e fitness function is used to evaluate and choose the chromosomes. A cross operator is used to create a new chromosome, and a mutation operator is applied to add randomness to the chromosomes, as the idea introduced in [37]. e algorithm is shown in Algorithm 4. Details of functions in Algorithm 4 are further introduced in Algorithms 5-9 . A general flowchart of Algorithm 4 is shown in Figure 15.
In Algorithm 4, necessary parameters are set in lines 1-3. en the area S is scattered into grids, helping to calculate the fitness function value and enable us to decode the best chromosome. In line 5, the swarm is generated using a function GenerateChromosome(), and the chromosome with the biggest fitness value is chosen as the current best chromosome, which is presented in detail in Algorithm 5. en the loop enters lines 7-14. In each iteration, the fitness values of chromosomes in the swarm are first calculated.
en go through selection operation as introduced in Algorithm 6, cross operation as introduced in Algorithm 7, and Input: valid area S, radius of approximate disk R Output: set S best consists of best particles with different lengths (1) Set parameters: swarm size N swarm , iteration number N in (2) Set parameters for PSO: ω, c 1 , c 2 (3) S best ←ϕ (4) Group number range: (16) end for (17) g best ←arg max f(p), p ∈ p 1 best , . . . , p N swarm best (18) end for (19) S best ←S best ∪ g best (20) end for ALGORITHM 1: Dynamic PSO.

Mathematical Problems in Engineering 13
Function GenerateChromosome() randomly generates a sequence whose length is the same as the number of grids in S, with each element corresponding to a grid. e sequence only consists of 0 and 1, with 1 indicating a cluster to be deployed at the center of a grid. To generate such a sequence, i.e., a chromosome for the GA, we first Input: N S , p init Output: random chromosome c (1) function GENERATECHROMOSOME (N S , p init ) (2) Generate N S random numbers obey normal distribution: c � [c 1 , . . . , c N S ] (3) for i � 1, 2, . . . , N S do (4) if c i > p init then end if (9) end for (10) end function ALGORITHM 5: Function GenerateChromosome().

Input: S, c best
Output: positions for station groups P (1) function DECODE (S, c best ) end if (7) end for (8)  generate N s random values obeying normal distribution, as line 2 in Algorithm 5. If the value is larger than parameter p init , it is set as 1, and 0 otherwise, as lines 4-8.

Algorithm 5 has a complexity of O(N s ).
It is obvious that the choice of p init affects the number of clusters deployed initially and affects the optimization process indirectly. us, setting a proper p init is important. Our suggestion of setting p init is to follow where P is a probability distribution and S c is the area covered by one cluster. In this paper, the normal distribution is chosen for P. en, A proper p init can be obtained by solving this equation. Function Selection() presented in Algorithm 6 takes current swarm C and the fitness values of these chromosomes F as inputs. With a roulette scheme, half of the chromosomes are selected into the next generation, which means N selection � (1/2)N swarm . e roulette scheme is achieved by creating a floor sequence R floor and ceiling sequence R ceil based on the ratio of a chromosome's fitness value to the sum of the fitness values, shown in lines 2-9. R floor and R ceil form a sequence of slots. A random number is generated and falls into a slot. e chromosome corresponding to the slot is selected, as in lines 12-13. Algorithm 6 has a complexity of O(N swarm ).
e selection operation presented in Algorithm 6 is a general operation in GA and it works well in various scenarios. However, it performs poorly in our case, as shown in Figure 16. e reason is that the fitness values of a swarm are similar, resulting in a roulette scheme unable to select proper chromosomes for the next generation effectively. is problem is solved by adopting an elite chromosome scheme. e idea is to choose some best chromosomes and put them directly into the next generation. en, lines 11-14 in Algorithm 6 are replaced with the following: (i) Sort elements in C as descending order of elements in F (ii) C new ←c 1 , c 2 , . . . , c N selection ∈ C N selection is the number of chromosomes directly passed to the next generation. To add more randomness, we set N selection � (1/4)N swarm . From Figure 16, it can be seen that this elite chromosome scheme can ensure that the fitness value never decreases. Besides, the algorithm can automatically adjust the number of clusters. Comparing with the roulette scheme, the elite chromosome scheme can achieve a higher fitness value and a coverage rate with fewer clusters.
is point is clearly shown in Figure 17, where we exhibit the deployment of clusters according to solutions obtained with GA using two different selection operators. In both cases, 35 clusters are deployed. e solution with the elite scheme resulting in clusters spreading more evenly and covering more area than that of the roulette scheme.
Function Cross() in Algorithm 7 generates the rest chromosomes. It randomly picks two chromosomes as parents, exchanges their elements after a random point, and generates two children (lines 2-6).
en, in lines 7-8, a random child is chosen and added to the swarm. Algorithm 7 has a complexity of O(N swarm ).
Function Mutation() in Algorithm 8 randomly generates N S numbers. If the number is larger than the mutation parameter μ, the corresponding element in the chromosome is inverse from 0 to 1 or from 1 to 0. Attention should be paid that μ should be set as a relatively small number, as a big μ introduces 1s rapidly during iterations as shown in Figure 18, overwhelming the algorithm's ability to adjust cluster numbers. Algorithm 8 has a complexity of O(N s × N swarm ).
Function Decode() in Algorithm 9 decodes the best chromosome into a set of coordinates for cluster deployment. Algorithm 9 has a complexity of O(N s ).

Comparison of the Two Schemes.
Experiments have been carried out to test the algorithms provided in Sections 4.2 and 4.3. e results in a rectangle region and a randomly generated irregular region are presented here, as shown in Figure 19. For the PSO, the swarm size is set as 20 and the iteration number is 50. e station group number increased from the number of clusters necessary to just fulfill the area ignoring the shape, which is 21. e solution of GA contains 29 clusters. To make it comparable, in DPSO, we also set the maximum number of clusters as 29. It turns out that the coverage can reach 87.98%. For the Elite GA, the swarm size is set as 50 and iteration number 200. Parameter μ is set as 0.00001 and p init is automatically set with the method introduced in Section 4.3. e solution provided by the Elite GA contains 29 clusters, with a coverage rate reaching 89.87%, slightly better than that of the Elite DPSO. Taking the randomness of both algorithms into consideration, we cannot claim the Elite GA is better than the Elite DPSO. In fact, the algorithms have different features, and some discussion is provided in the following.
DPSO can provide a set of particles corresponding to a different number of groups. us, users can choose the best one according to the cost and the coverage rate. GA automatically decides how many clusters are needed, and users can directly use the scheme generated. However, users are unable to assign the number of clusters.
In DPSO, base station groups move continually in the area of interest, while, in GA, base station groups jump among grids. is means DPSO has more flexibility than GA. If we are seeking for a more precise solution with GA, the area of interest has to be meshed into more grids, increasing the length of the chromosome and resulting in more computation. For example, if a region of 100 × 100 m 2 is scattered into 1 × 1 m 2 grids, the length of the chromosome is (100 × 100/1 × 1) � 10 4 . In case the region is rather complex, we have to divide it into grids of 0.1 × 0.1 m 2 . en, the length of the chromosome becomes (100 × 100/0.1 × 0.1) � 10 6 , increasing the computation dramatically.
In conclusion, both methods can provide a scheme to deploy base station clusters, and the choice depends on the users' preference.

Experiment Validation of UGV Navigation Based on UWB Cluster Positioning
In Section 3, the performance of a UWB cluster is tested and validated, and then in Section 4, the deployment methods of UWB clusters are proposed and validated by simulations. In this section, the outcomes of the previous sections are applied to UGVs. Based on the experiments in Section 3, 50 cm side length is chosen for a UWB cluster and 10-meter distance is determined as the coverage radius of a UWB base station. UGVs are then set out to follow different paths, showing that under the positioning of UWB groups, errors of UGV trajectories are constrained around 25 cm. With UWB clusters, UGVs can move smoothly to their destinations.

Design of UGVs.
In order to validate the navigation performance based on the UWB cluster positioning system, UGV prototypes are fabricated. e control kernel hardware is composed of an MCU, a UWB tag, an IMU, a WiFi module, and actuators. e structure is shown in Figure 20. e MCU acts as the "brain" of a UGV, receiving input information and making control decisions. A 32-bit MCU,    STM32H750VB, is adopted, which includes an Arm Cortex-M7 with 128 Kbytes of flash memory, 1 MB RAM, 480 MHz CPU, L1 cache, large set of peripherals. e UWB tag uses DWM1000 of Decawave. For IMU, MPU-9250 is applied, which has a 9-axis motion processing unit including a gyroscope, an accelerometer, and a digital compass. e WiFi connection is established using a low-cost WiFi microchip, ESP8266, with a full TCP/IP stack. e operation of DC motors is driven by A4950s, full-bridge DMOS PWM motor drivers. e UGV prototype is shown in Figure 21.

UGV Automated Driving under UWB Positioning.
In this part, experiments are carried out with UGVs using a UWB cluster positioning system. First, UGV automated driving is realized in an 8 m×9 m room, as shown in Figure 22. A UGV is programmed to follow designated trajectories shown in Figure 23. It can complete fixed-point back and forth trajectory, square trajectory, and circular trajectory. For each test, target positions or waypoints are generated. e UGV checks whether the current position is within a circle around these positions to guarantee that it follows the designated paths.
e live demonstration video is provided in supplemental materials. e test results of a circular trajectory are depicted in Figure 24 and the corresponding positioning errors are provided in Figure 25. In all the experiments, the position errors are within 30 cm as shown in Figure 26, which are consistent with test results in Section 3.2.

Conclusion
In this paper, the concept of UWB clusters composing of multiple UWB base stations as one device is proposed to provide positioning service for robots. Tests in Section 3 show that this positioning system works well. Within 8 meters, the mean distance errors are less than 20 centimeters in most cases. Within 10 meters, the mean position errors are no more than 60 centimeters in most cases. By approximating the overlap area covered by all the three base stations in one group as a circle, two algorithms are proposed to obtain the best deployment scheme in an area of interest with random shape. Both algorithms can help determine the number of groups necessary. Furthermore, the UWB positioning system is verified by conducting multiple tests with UGVs. e results confirm that UGVs can follow complex trajectories precisely. e UWB cluster's ability to support multiple robots highlights its potential application in swarm robotics. Furthermore, combining multiple base stations as one positioning device, it can be equipped on a larger vehicle such as an unmanned surface vehicle or a heavy drone, providing high-accuracy positioning service for a swarm of smaller robots. Note that, in this paper, we mainly focus on the feasibility of the proposed idea, and the performance improvement of the positioning system is reserved as future work.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.