Abstract

Localization is crucial to wireless sensor networks. Among the recently proposed localization algorithms, the mobile anchor-assisted localization (MAL) algorithm seems promising. A MAL algorithm using a single mobile anchor has low energy consumption but a high localization error. Conversely, a MAL algorithm with three or more mobile anchors has minor localization errors but high energy consumption. By balancing energy consumption and localization accuracy, our study developed a localization algorithm assisted by two mobile anchors. A mobile anchor traverses the network along a double anchor SCAN (DASCAN) path, which divides the deployment region into grids and requires the two mobile anchors to traverse different horizontal lines in a zigzag pattern. Sensor nodes estimate their locations using a multiple-disturbance strategy grey wolf optimization (MDS-GWO) algorithm, which improves optimization by introducing a nonlinearly decreasing weight, a random perturbation of grey wolves and a mirror grey wolf. Using MATLAB, DASCAN was compared with GTURN, GSCAN, PP-MMAN, H-Curves, M-Curves, and SCAN paths by their energy consumption and localization rates. The localization error of MDS-GWO was compared with trilateration, PSO, WOA, and GWO. The impacts of radio irregularity, radio radius, and fading effect on MDS-GWO with different paths were also analyzed. The simulation results showed that the energy consumption of DASCAN was, on average, 30.1% less than GSCAN, GTURN, and PP-MMAN, but they had almost the same localization accuracy. The energy consumption of DASCAN was an average of 18.67% more than M-Curves, H-Curves, and SCAN, but the localization error of DASCAN was average of 32.3% less than SCAN, H-Curves, and M-Curves. The localization error of MDS-GWO was average of 25.5% less than trilateration, PSO, WOA, and GWO. Moreover, the performance of the proposed algorithm was less affected by different setups than the compared methods.

1. Introduction

The wireless sensor network (WSN) consists of a large amount of sensor nodes deployed in a given region of interest (ROI). It is widely applied to environmental surveys, habitat monitoring, medical diagnosis, and disaster rescue. Node localization is a key technology of WSNs. However, equipping each sensor node with a global navigation satellite system (GNSS) device is impractical due to cost and energy consumption constraints. Only a few sensor nodes, called anchors, know their positions. Other sensor nodes, unknown nodes, are localized with the help of beacons broadcast by anchors.

More anchors result in higher localization accuracy, but the anchors are more expensive than unknown nodes, so the mobile anchor-assisted localization (MAL) algorithms have attracted considerable interest. MAL algorithms require the mobile anchor to traverse the ROI along a given path and broadcast beacons periodically. The locations where the mobile anchor broadcasts beacons are called virtual anchors. If the MAL algorithm only uses one mobile anchor, obtaining high localization accuracy with a short path is difficult. If the MAL algorithm uses three or more mobile anchors, localization becomes more accurate, but at the cost of high energy consumption. This paper proposes a MAL algorithm that uses two mobile anchors, balancing energy consumption and localization accuracy. Our algorithm uses a double anchor scan (DASCAN) as the moving path and a multiple disturbance strategy for the grey wolf optimizer (MDS-GWO) to estimate the locations of unknown nodes. The contributions of our study are: (1)It presents the DASCAN path for two mobile anchors, reducing the number of beacons and path length. This allows the two anchors to scan adjacent rows of ROI simultaneously so that three neighboring virtual anchors form an equilateral triangle.(2)It presents an improved GWO algorithm with multiple disturbance strategies to estimate the locations of unknown nodes, which can improve localization accuracy.(3)It reports extensive experiments that compared the performance of the proposed algorithm with similar current algorithms. DASCAN was compared with GTURN, GSCAN, PP-MMAN, H-Curves, M-Curves, and SCAN paths in terms of energy consumption and localization rate. MDS-GWO was compared with trilateration, particle swarm optimization (PSO), whale optimization algorithm (WOA), and grey wolf optimizer (GWO) regarding localization errors. The impacts of radio irregularity, radio radius, and fading effect on MDS-GWO with different paths were also analyzed.

The rest of this paper is organized as follows. Section 2 introduces related work, and Section 3 presents the problem statements and performance metrics. Section 4 describes the proposed algorithm, including the DASCAN path and MDS-GWO algorithm. Section 5 summarizes the simulations and analysis. Section 6 concludes the paper. Table 1 lists the acronyms used in this paper and their definitions.

The MAL algorithm needs to define a moving trajectory for a mobile anchor, and the location algorithm to estimate the location of unknown nodes. The path-planning algorithm may be dynamic or static [1]. The dynamic form determines the path according to the distribution of unknown nodes and environments, which may not cover all unknown nodes and requires additional hardware support. The static form determines the path in advance and requires the mobile anchor to move along a given path during localization, which is cheaper than the dynamic approach. Song et al. [2] proposed a path for self-adaptive anchors based on the Gaussian–Markov model, and they applied an alternating minimization algorithm to estimate the positions of unknown nodes. However, the random path may not cover the whole ROI. SCAN, Double-SCAN, and Hilbert proposed in [3] are simple, but their parameters must be carefully designed to avoid collinearity. H-curve [4] is an H-shaped path that generates collinear beacons. SLMAT (a mobile anchor node based on trilateration and scan) [5] ensures that each unknown node is covered by an equilateral triangle formed by beacons. ∑-SCAN [6] combines SCAN and zigzag paths to achieve high localization accuracy and cover the ROI with a short path. M-curve [7] is an M-shaped path and applies dolphin optimization algorithm to localize unknown nodes, but it may not localize unknown nodes near the borders of ROI. The static search-and-decide (SSD) and dynamic search-and-decide (DSD) [8] paths are, respectively, static and dynamic, which consist of a search phase and a decision phase. In the first phase, SSD visits a subset of virtual anchors to determine the grids occupied by unknown nodes. In the decision phase, the mobile anchor revisits the grids containing unknown nodes for localization. DSD differs from SSD in the second phase: It generates anchors based on perpendicular bisectors. The localization accuracy of these MAL algorithms with a single mobile anchor is low, but they consume comparatively little energy. GSCAN and GTURN [9] use three mobile anchors. The mobile anchors in GSCAN repeatedly broadcast beacons at the same location; GTURN avoids this problem using a longer path. However, GSCAN and GTURE need to use two boundary strategies to guarantee that they can cover the whole ROI. PP-MMAN [10] also uses three mobile anchors, but it requires them to move horizontally simultaneously to overcome the drawbacks of GSCAN and GTURN. These MAL algorithms using multiple mobile anchors improve localization accuracy at the cost of high energy consumption. The enhanced RSSI-based tree-climbing mechanism (ERTC) [11] requires the mobile anchor to be equipped with both omnidirectional and directional antennae, where the omnidirectional antenna broadcasts the message and directional antennas to receive messages of sensor nodes. It identifies the trajectory of the mobile anchor with the virtual force of unknown nodes in the network, and it uses a circumcenter algorithm to localize unknown nodes. It is more complicated than the other MAL algorithms and needs additional hardware support.

The location of unknown nodes can be estimated as an optimization problem. Nowadays, swarm intelligence optimization has been widely applied in localization due to its low complexity and easy implementation. H-Best PSO [12] is an improved PSO algorithm to estimate the locations of unknown nodes where the mobile anchor traverses along a Hilbert curve. Song et al. [13] proposed a localization algorithm based on glowworm swarm optimization of a hybrid chaotic strategy to control the moving distance of each firefly by chaos mutation and chaotic inertial weight when the firefly falls into a local optimum. The algorithm proposed in [14] reduces the hop distance error by leading to the average hop distance error correction value. It applies differential evolution to optimize the localization result of the unknown node. Considering the impact of obstacles, the algorithm proposed in [15] divides the anchors into multiple groups. If the anchors of a group cannot accurately localize the unknown node, the other anchors in the nearby groups are used to assist localization. In [16], WOA is applied to optimize the relationship between RSSI (received signal strength indicator) and signal transmission distance to improve localization accuracy.

Because the traversal of mobile anchors forms some specific curves, the location of unknown nodes can be estimated by some geometry principles. For example, PI (perpendicular intersection) [17] is based on the perpendicular intersection principle, which utilizes the geometric relationship of a perpendicular intersection to compute the locations of unknown nodes. The algorithm proposed in [18] is based on the principle that the perpendicular bisector of a chord of a circle passes through the center of the circle, so the unknown node is localized by two chords constructed by a mobile anchor. These methods are easy to implement, but they heavily depend on the mobile anchor’s path and the density of virtual anchors.

The MAL algorithm with one mobile anchor has difficulty obtaining high localization accuracy with low energy consumption, while the MAL algorithm using three or more mobile anchors has high localization accuracy but high energy consumption. To balance energy consumption and localization accuracy, this paper proposes a MAL algorithm with two mobile anchors. DASCAN produces a path that is as short as possible and provides sufficient virtual anchors to each unknown node. Compared with the other evolutionary algorithms, GWO has high optimization efficiency and is simple to implement because of fewer control parameters. MDS-GWO makes the best of these advantages and improves its ability to jump out of local optima.

3. Problem Statements and Performance Metrics

3.1. Problem Statement

A WSN consists of a set of unknown nodes randomly deployed in an ROI with height and length . The two mobile anchors are denoted by and . Let be the actual location of , and be its estimated location, denoted by . The proposed algorithm applies RSSI to measure the distance between the mobile anchor and the unknown node. The RSSI-based ranging technique is based on where is the received power at distance , denotes the received power at reference distance , is the path-loss exponent, and is a log-normal random variable with variance to account for fading.

Ideally, the communication range of a sensor node is a circle with radius . However, signal propagation is easily affected by the environment. This paper uses the degree of irregularity (DOI) to represent communication irregularity. DOI is defined as: where represents the DOI in the -th direction and satisfies , and is a random number in (0,1).

Let the neighboring virtual anchors of unknown node be where is the number of neighboring virtual anchors. The objective of localization is formulated as where is the location of , and is the measured distance between and . is the unknown quantity to be found, and its solution is .

3.2. Performance Metrics

The localization error of is defined as

An unknown node in two-dimensional space can be localized only if it receives at least three noncollinear beacons. Let be the number of localized unknown nodes, and be the localization rate, that is

The average localization error is

In addition, the root mean square error (RMSE) of localization error is

The Cramer–Rao lower bound (CRLB) provides the lower bound on the variance of any unbiased estimator. This paper uses the CRLB introduced in [19] as the standard benchmark where is a fixed constant. Let be the union set of unknown nodes and virtual anchors, where is the set of unknown nodes and is the set of virtual anchors. The real position of is , so for and for Let be the CRLB of error variance of estimating ; it is proved in [19] that where

In the abovementioned equations, is the set of sensor nodes and virtual anchors that make a pairwise observation with . is defined as

Based on the abovementioned equations, the RMSE of localization error variance is

The energy consumption of the mobile anchors depends on the path length and number of virtual anchors: where and are, respectively, the energy consumption of broadcasting a virtual anchor and moving per meter; M and L are, respectively, the numbers of virtual anchors and path length.

Table 2 presents a list of symbols used in this paper.

4. Proposed Methodology

4.1. DASCAN Path

As Figure 1 shows, ROI is divided into horizontal rows, numbered from bottom to top, and the gap between neighboring rows is . In addition, ROI is divided into vertical columns horizontally, numbered from left to right, and the gap between two neighboring columns is .

The virtual anchors on even and odd rows are, respectively, indicated by and . The bottom left point of ROI is taken as the coordinate origin, so the locations of virtual anchors are: (1), where is an even number, and.(2), where is an odd number, and .

Because the virtual anchors deployed as equilateral triangles can achieve better localization accuracy [9, 10], DASCAN requires and to be equilateral triangles, so . The points in Figure 1 are virtual anchors. Algorithm 1 shows the process of obtaining DASCAN, where the trajectory of the two mobile anchors is carefully designed to be as short as possible.

Input: .
Output: : The lists of virtual anchors of and .
1.  ,
2.  ,,,
3.  While ()
4.   If
5.    
6.    
7.    
8.   Else if
9.    
10.    
11.    
12.   End if
13.   
14.  End while
15.  If
16.   If
17.    
18.   Else if
19.    
20.    
21.   End if
22.  Else if
23.   If
24.    
25.   Else if
26.    
27.    
28.   End if
29.  End if

In Algorithm 1, and are, respectively, the lists of virtual anchors to be traversed by and . Lines 1 and 2 initialize the two lists and parameters. Lines 3–14 generate the two lists when . The variable decides the traversal direction. and move from left to right if ; they move from right to left if . Lines 15–29 decide and when . If , traverses row and traverses row . If , traverses the top row. Moreover, the traversal direction also depends on .

According to different and after the while (line 14), DASCAN has four types of paths, as Figure 2 shows.

4.2. Location Estimation Using Improved GWO Algorithm

GWO has a higher optimization efficiency and is easier to implement than other evolutionary algorithms. Its main drawback is that, once it falls into a local optimum, it is difficult for the algorithm to jump out. Thus, we propose the MDS-GWO algorithm to estimate the positions of unknown nodes.

4.2.1. Improved GWO Algorithm

GWO is a swarm intelligence optimization algorithm that imitates a wolf group’s hunting behavior. In GWO, the best, the second-best, and the third-best wolves are, respectively, denoted by , , and . The rest of the wolves are denoted by . When hunting, directs to surround, hunt, and attack their prey. In addition, and assist the work of then capture the prey.

The mathematical model of encircling prey is as follows: where is the current iteration number, is the position vector of prey, and is the position vector of a grey wolf. The coefficient vectors and are where and are random vectors in (0,1), and each element of is where is the maximum iteration number.

The mathematical model of hunting prey is as follows: where , , and , respectively, represent the position vectors of , , and in the current iteration; , , and are the distances between the current grey wolf and the three best wolves.

Inspired by [20, 21], this paper introduces multiple disturbance strategies into GWO to conquer the drawback of GWO, and this improved GWO algorithm is referred to as MDS-GWO. The improvements are: (1)Nonlinear decrease of. The convergence speed of the GWO algorithm depends on the speed at which decreases from 2 to 0. Using Equation (16), decreases at the same speed, which cannot balance the global and local search, so this paper updates by

Based on Equation (20), decreases quickly in the early stage to make GWO converge fast, and decreases slowly later to improve the ability of global search while preventing GWO from falling into a local optimum. (2)Random perturbation. Furthermore, random perturbation is introduced to GWO to enhance the ability of global search. For each wolf of , a new individual is introduced aswhere and are, respectively, the upper and lower bounds of the search space, and is a random number in [0,1], so the grey wolf is updated as where and are the fitness values of , respectively. (3)Mirror grey wolf. When a grey wolf surrounds the prey, a mirror wolf is produced at a position symmetrical to the prey, which is:where is the position of the mirror wolf, and and are random vectors in [0,1]. The grey wolf is updated as

Finally, the grey wolf is updated according to the value of , defined by where is a predefined constant, and is a random number in [0,1]. In each iteration, if , the grey wolf has a 50% chance to update its position using Equation (22); If , the grey wolf has a 50% chance to update its position using Equation (24). In other cases, the grey wolf updates its position using Equation (19).

4.2.2. Location Estimation Using MDS-GWO

The virtual anchors of DASCAN form equilateral triangles, so the neighboring virtual anchors of an unknown node are likely to form equilateral triangles. Therefore, an unknown node can choose three neighboring virtual anchors forming an equilateral triangle to get an initial estimation. The location estimation of unknown node consists of three steps: (1)Use an RSSI-weighted centroid algorithm to estimate the approximate location (2)Determine the search space of MDS-GWO for according to (3)Use MDS-GWO to estimate the coordinates of

Each unknown node maintains a virtual anchor list, , which contains the locations and RSSIs of neighboring virtual anchors. The virtual anchors with the first- and second-largest RSSIs are, respectively, denoted by , and . If a virtual anchor in can form an equilateral triangle with and , it is chosen as . Otherwise, the virtual anchor with the third-largest RSSI is chosen as . is calculated as where . The search space is .

Algorithm 2 presents the details of MDS-GWO-based location estimation. Note that this algorithm is run by each unknown node simultaneously.

Input: , ,
Output:
1.  the virtual anchor with the largest RSSI in
2.  the virtual anchor with the second-largest RSSI in
3.  If ( satisfying is an equilateral triangle)
4.   
5.  Else
6.   the virtual anchor with the third-largest RSSI in
7.  End if
8.  Get by eq. (26)
9.  
10.  Initialize a grey wolf population in
11.  Initialize , , and by Eq. (14), Eq. (15), Eq. (20), and Eq. (25)
12.  Calculate the fitness of each wolf by Eq. (3)
13.  the best wolf
14.  the second-best wolf
15.  the third-best wolf
16.  
17.  While
18.   For each wolf
19.    Get by Eq. (19)
20.    a random number in [0,1]
21.    If
22.     Get by Eq. (21)
23.     Get the final by Eq. (22)
24.    Else if
25.     Get by Eq. (23)
26.     Get the final by Eq. (24)
27.    End if
28.   End for
29.   Update , , and by Eq. (14), Eq. (15), Eq. (20), and Eq. (25)
30.   Calculate the fitness of each wolf by Eq. (3)
31.   Update , and
32.   
33.  End while
34.  Return

In this algorithm, lines 1–7 got three virtual anchors , , and , and lines 8 and 9 decided the search space based on these three virtual anchors. Lines 10–16 initialize the wolf population and necessary parameters of MDS-GWO, and lines 17–33 are the main process of location estimation by MDS-GWO. When it has not reached , lines 18–28 update each grey wolf according to different conditions, and lines 29–32 update the parameters necessary to prepare for the next iteration. Finally, line 34 returns the position of wolf as the estimated location of unknown node .

5. Results and Discussion

This section reports the performance of the proposed algorithm compared with the other algorithms. It also presents the impact of parameters on localization performance in terms of energy consumption, localization rate, and localization error. The simulation parameters are shown in Table 3. The other parameters are the same as in the corresponding references.

5.1. Comparisons of Different Paths

This section compares the different paths of mobile anchors of GTURN, GSCAN, PP-MMAN, H-Curves, M-Curves, SCAN, and DASCAN. (1)Energy Consumption of Mobile Anchors. The energy consumption of a mobile anchor depends on the path length and the number of virtual anchors. Figure 3 presents the lengths of all paths, which become shorter as increases. On average, GTURN is the longest path, followed by GSCAN, PP-MMAN, DASCAN, H-Curves, M-Curves, and SCAN. GTURN, GSCAN, and PP-MMAN use three mobile anchors, which are longer than the other paths. DASCAN is shorter than the other multianchor paths and longer than one-anchor paths. On average, DASCAN is 53%, 36.5%, and 14% shorter than GTURN, GSCAN, and PP-MMAN, respectively; it is 14%, 15.5%, and 34% longer than M-Curves, H-Curves, and SCAN, respectively.

DASCAN uses two mobile anchors with a carefully designed path, while GSCAN, GTURN, and PP-MMAN use three mobile anchors and boundary compensations, so DASCAN is shorter than GSCAN, GTURN, and PP-MMAN. By contrast, M-Curves, H-Curves, and SCAN use a single mobile anchor to traverse the ROI without repeated scans, so they are shorter than DASCAN.

Figure 4 shows the number of virtual anchors. All paths generate fewer virtual anchors as increases. Generally, GSCAN generates the most virtual anchors, followed by GTURN, PP-MMAN, DASCAN, M-Curves, H-Curves, and SCAN. The three mobile anchors of GSCAN repeatedly broadcast beacons at the same locations, and GSCAN provides the second-longest path, so it has the most virtual anchors. However, SCAN traverses the ROI along straight lines, and its path is the shortest, so it has the fewest virtual anchors. The proposed DASCAN has fewer virtual anchors than the other multianchor paths and more virtual anchors than one-anchor paths. The number of virtual anchors of DASCAN is, respectively, 11% and 53% fewer than GTURN and GSCAN, and it is, respectively, 12.3% and 32% more than H-Curves and SCAN. The number of virtual anchors of DASCAN is slightly less than PP-MMAN and slightly more than M-Curves.

The boundary compensation methods of GSCAN, GTURN, and PP-MMAN require these paths to generate more virtual anchors, so GSCAN has redundant virtual anchors. By contrast, DASCAN has no duplicate virtual anchors or boundary compensation, so it has fewer virtual anchors than GSCAN, GTURN, and PP-MMAN. SCAN traverses ROI with straight lines and generates virtual anchors periodically, while H-Curve and M-Curve generate virtual anchors at each turn of the path. DASCAN makes slightly more turns than M-Curve, so it has more virtual anchors than SCAN, H-Curve, and M-Curve.

Figure 5 shows the energy consumptions. GSCAN consumes the most energy, followed by GTURN, PP-MMAN, DASCAN, M-Curves, H-Curves, and SCAN. Since GSCAN and GTURN are the longest paths and have the most virtual anchors, they consume the most energy. DASCAN consumes less energy than the other multianchor paths and consumes more energy than one-anchor paths. On average, the energy consumption of DASCAN is, respectively, 42.2%, 38.8%, and 9.3% less than GSCAN, GTURN, and PP-MMAN, and 9%, 14%, and 33% more than M-Curves, H-Curves, and SCAN.

Although DASCAN consumes less energy than GTURN, GSCAN, and PP-MMAN, their localization rates and error levels are almost identical. The localization rate and accuracy of DASCAN are higher than H-Curves, M-Curves, and SCAN, although it consumes more energy. (2)Localization Rates. In this section,  m. As Figure 6 shows, the localization rates of all paths decrease as DOI increases because the larger DOI causes the unknown nodes to receive fewer virtual anchors. The localization rates of all paths increase as increases because a larger causes the unknown node to receive more virtual anchors.

GSCAN, GTURN, PP-MMAN, and DASCAN use multiple anchors, and each anchor has a radio range, so they are less affected by DOI than other paths using a single mobile anchor. The average localization rate differences of DASCAN, GSCAN, GTURN, and PP-MMAN are less than 2%, and their average localization rates are higher than H-Curves, M-Curves, and SCAN. When , the average difference in GSCAN, GTURN, PP-MMAN, and DASCAN is only 0.56%, but the average localization rate of DASCAN is, respectively, 13%, 38.3%, and 43.4% higher than M-Curves, H-Curves, and SCAN.

5.2. Comparisons of Location Estimation Algorithms

This section compares the location estimation algorithms in terms of the localization error. These algorithms, with different paths, included trilateration, PSO, WOA, GWO, and MDS-GWO. In this section, and are used. The average localization error over all is shown in Table 4. Figure 7 shows the results for different values.

The localization errors of trilateration and PSO were the two largest. Those of GWO and WOA were the same. With each path, MDS-GWO achieved the highest localization accuracy. In all cases, trilateration with SCAN had the lowest localization accuracy. The localization errors of MDS-GWO with GTURN, GSCAN, PP-MMAN, and DASCAN were almost identical and smaller than MDS-GWO with the other paths. Considering all paths, the average localization errors of trilateration, PSO, WOA, GWO, and MDS-GWO were, respectively, 0.259, 0.240, 0.213, 0.202, and 0.183 m. The localization error of MDS-GWO was, respectively, 42%, 31.7%, 16.9%, and 11.4% less than trilateration, PSO, WOA, and GWO. Considering MDS-GWO with different paths, the localization accuracy of DASCAN was almost the same as GTURN, GSCAN, PP-MMAN, and the localization error of DASCAN was, respectively, 70.6%, 16.3%, and 10% less than SCAN, H-Curves, and M-Curves. The virtual anchors of DASCAN, GSCAN, GTURN, and PP-MMAN formed equilateral triangles. The localization rates of these paths were almost the same as mentioned above, so their localization errors were almost the same. SCAN, H-Curves, and M-Curves use only one mobile anchor, which is heavily affected by DOI, so their localization errors were higher than the other paths.

Figure 8 shows the RMSE of localization algorithms when m,  l, , and the mobile anchors use DASCAN as the moving path. MDS-GWO had the minimal RMSE, followed by GWO, WOA, PSO, and trilateration. On average, the RMSE of MDS-GWO was, respectively, 8.2%, 19%, 35.1%, and 38% less than GWO, WOA, PSO, and trilateration. Figure 9 shows the CDF of localization errors. MDS-GWO localized all sensor nodes with an error of less than 0.3 m, and more than 90% of the unknown nodes were localized with an error of less than 0.1 m. Comparatively, GWO and WOA, respectively, localize about 85% and 80% of unknown nodes with localization errors less than 0.1 m. Furthermore, the RMSE of MDS-GWO was only 0.03 larger than . MDS-GWO made the most of GWO’s advantages and applied multiple strategies to jump out of the local optimum, so we concluded that MDS-GWO was superior to the other localization methods.

5.3. Impacts of Different Parameters on Localization Errors

This section analyzes the impact of DOI, , and on localization errors of MDS-GWO over different paths. (1)Impact of DOI. In this part, m and  m. Figure 10 shows that SCAN had the most significant localization error under all DOI, followed by H-Curves and M-Curves. With the increase in DOI, the localization errors of all paths increased. For every 0.05 increase in DOI, the localization error of SCAN increased by 0.02 m on average, which is the largest among all paths. DASCAN was less affected by DOI than H-Curves, M-Curves, and SCAN, and almost the same as GTURN, GSCAN, and PP-MMAN. Under all DOI, the localization error of DASCAN was the same as GTURN, GSCAN, and PP-MMAN and less than H-Curves, M-Curves, and SCAN. The localization error of DASCAN was, respectively, 0.03, 0.02, and 0.07 m less than H-Curves, M-Curves, and SCAN.(2)Impact of. In this part,  m, m, and . As Figure 11 shows localization errors of all paths increased with . The localization error of SCAN was the largest in all cases because unknowns may not receive sufficient beacons to be localized, and many virtual anchors of SCAN were collinear. The localization errors of GTURN, GSCAN, PP-MMAN, and DASCAN were almost the same. For every 5 m increase in , the localization errors of M-Curves and SCAN increased by 0.013 m and 0.15 m, respectively, while those of GTURN, GSCAN, PP-MMAN, DASCAN, and H-Curves grew by less than 0.01 m. Of all paths, DASCAN was the least affected by . Under all values, the localization error of DASCAN was slightly smaller than GTURN, GSCAN, and PP-MMAN and was significantly less than H-Curves, M-Curves, and SCAN. On average, the localization error of DASCAN was, respectively, 0.05, 0.03, and 0.08 m less than H-Curves, M-Curves, and SCAN.(3)Impact of. In this part,  m,  m, , and . Figure 12 shows that the localization error of SCAN was the largest, followed by H-Curves and M-Curves. For every two increases of , the localization errors of H-Curves, M-Curves, and SCAN, respectively, increased 0.23, 0.218, and 0.27 m on average, and the localization errors of GTURN, GSCAN, PP-MMAN, and DASCAN increased by 0.2 m. DASCAN was less affected by than H-Curves, M-Curves, and SCAN, and it was almost the same as GTURN, GSCAN, and PP-MMAN. On average, the localization error of DASCAN was, respectively, 0.1, 0.05, and 0.2 m less than H-Curves, M-Curves, and SCAN.(4)Analysis of RMSE and CDF. In this part,  m,  m, and . Figure 13 shows that the RMSE of seven paths with DASCAN had the least RMSE, followed by GSCAN, GTURN, PP-MMAN, M-Curves, H-Curves, and SCAN. The RMSE of DASCAN was, respectively, 20.3%, 21.7%, and 29% less than H-Curves, M-Curves, and SCAN, and it was only 0.015 higher than . As Figure 14 shows, the CDF of DASCAN, GTURN, and GSCAN increased faster than the other paths, which means they could localize all unknown nodes with higher accuracy than the other paths. After a localization error of 0.11 m, the CDF of DASCAN was always larger than the other paths. For example, DASCAN localized more than 89% of unknown nodes with a localization error of less than 0.15 m, which was the largest among all paths.

6. Conclusions

This paper proposes a localization algorithm using two mobile anchors to balance the localization accuracy and energy consumption. The study developed a specified path, DASCAN, which uses two mobile anchors to traverse different rows after dividing the ROI into regular grids. The neighboring virtual anchors generated by the two mobile anchors form equilateral triangles to provide noncollinear beacons for each unknown node to be localized. The paper also proposes an improved GWO algorithm (MDS-GWO) to estimate the positions of unknown nodes. MDS-GWO introduces multiple disturbance strategies into GWO to improve its ability to jump out of local optima.

Data Availability

Data sets are not used in this article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by the National Key R&D Program of China (Grant No. 2018YFC1406203), National Natural Science Foundation of China (Grant No. 62072287), Natural Science Foundation of Shandong Province, China (Grant No. ZR2021LZH004), and Science and Technology Support Plan of Youth Innovation Team of Shandong Higher School, China (Grant No. 2019KNJ024).