A Virtual Force Algorithm-Lévy-Embedded Grey Wolf Optimization Algorithm for Wireless Sensor Network Coverage Optimization

The random placement of a large-scale sensor network in an outdoor environment often causes low coverage. In order to effectively improve the coverage of a wireless sensor network in the monitoring area, a coverage optimization algorithm for wireless sensor networks with a Virtual Force-Lévy-embedded Grey Wolf Optimization (VFLGWO) algorithm is proposed. The simulation results show that the VFLGWO algorithm has a better optimization effect on the coverage rate, uniformity, and average moving distance of sensor nodes than a wireless sensor network coverage optimization algorithm using Lévy-embedded Grey Wolf Optimizer, Cuckoo Search algorithm, and Chaotic Particle Swarm Optimization. The VFLGWO algorithm has good adaptability with respect to changes of the number of sensor nodes and the size of the monitoring area.


Introduction
Wireless sensor networks (WSNs) were originally designed for military applications, however currently WSNs are also widely used in civilian applications, including vehicle tracking, forest monitoring, seismic observation, building monitoring, and water resource monitoring [1,2]. The coverage rate is an important measure of WSN performance. How to use the minimum number of sensor nodes to monitor the target area has been a question of particular interest regarding WSN technology. Wireless sensor nodes are typically placed randomly in the monitored area, which results in a non-uniform distribution of sensor nodes and a low coverage rate of the monitored area. Therefore, it is of great significance to improve the coverage rate of WSNs in the monitored area by adjusting the position of sensor nodes.
In 1999, the concept of swarm intelligence was formally proposed. Algorithms and solutions for a distributed problem inspired by the collective behavior of social insects and animal groups are part of swarm intelligence research. This type of algorithm is called a swarm intelligence optimization algorithm. Among intelligent algorithms, the Taguchi-teaching learning-based optimization algorithm and Jaya algorithm are also widely used [3,4]. In 2009, Yang and Deb of Cambridge University proposed a new biological heuristic algorithm, namely the Cuckoo Search algorithm [5]. This algorithm is widely used because of its simple structure, few control parameters, and strong search capability. Thanh-Phong Dao [6] proposed a hybrid Taguchi-Cuckoo Search (HTCS) to be applied to compatible focus positioning platforms (CFPP). The proposed HTCS approach can effectively optimize multiple objectives for the CFPP, and would be a useful technique for related optimization problems.

1.
This paper applies the LGWO algorithm to wireless sensor node deployment.

2.
For wireless sensor node deployment, this paper proposes the VFLGWO algorithm, which achieves a higher coverage rate.

3.
This paper designs sensor node matching algorithms to make the average moving distance of sensor nodes shorter.
The remainder of this paper is structured as follows. Section 2 describes the related work. Section 3 describes the proposed VFLGWO based sensor deployment scheme. The simulation results and comparison of the CPSO algorithm, CS algorithm, and LGWO algorithm are presented in Section 4, and conclusions are drawn in Section 5.

Research Situation
In Guanglin Zhang's work [8], Guanglin Zhang and Shan You studied the coverage of WSNs. They employed a Voronoi diagram to address the coverage problem. The authors transformed the network area coverage problem into cell coverage problems with the Voronoi diagram, which only needs to optimize local coverage for each cell in a decentralized way.
The research on using swarm intelligence optimization algorithms to solve WSN coverage optimization problems is gradually increasing.
Ning-ning Qin and Leopoldo Eduardo Cárdenas-Barrón [9,10] applied differential evolution and its improved algorithm to complete wireless sensor node deployment. Differential evolution has the characteristics of a simple structure, easy implementation, fast convergence, and strong robustness. Jia, Jie [11] adopted a genetic algorithm to solve the WSN coverage optimization problem. Genetic algorithms have a better ability to conduct a global search, but their implementation is relatively complex and their convergence speed is slow. The PSO algorithm [12][13][14] has been applied to the deployment of wireless sensor networks by scholars. The PSO is easy to implement and can effectively balance the directionality, diversity, and balance in the search. In other research [15], the PSO algorithm was employed to solve the optimization problem for a compliant rotary joint for an upper limb assistive device, and obtained a CRJ with a smaller error. However, the PSO algorithm has some shortcomings, including relying too much on search parameters, easily falling into local extremum, and premature convergence. In Celal Ozturk's work [16], an artificial Sensors 2019, 19, 2735 3 of 20 bee colony algorithm was used to maximize area coverage in WSNs. The artificial bee colony algorithm has the advantage of easy realization, and the shortcoming of easily falling into the local solution. Some scholars have adopted the Chaotic Artificial Fish Swarm algorithm to solve this problem [17]. The advantage of the Artificial Fish Swarm algorithm is that it has better robustness. The algorithm adopts an adaptive step size, which effectively increases the coverage of the algorithm.

An Overview of the GWO Algorithm
Inspired by wolf hunting behavior, Mirjalili et al. proposed the Grey Wolf Optimization (GWO) [18], which is easy to implement because of its few control parameters. The GWO algorithm is inspired by the social hierarchy and hunting strategies of grey wolves in the wild. Each grey wolf represents a potential solution to the population. The position of α wolf, which is the leading wolf, is the best position. The position of the β wolf in the second class of the wolf rank, and the δ wolf responsible for reconnaissance, vigilance, warding and guarding, are the second-best solution and the third best solution, respectively. The other candidate solutions are the lower-level ω wolf positions. In the conventional GWO, the motion of wolves is described as follows.

Population Initialization
Since the performance of the GWO is less affected by the initial value of the population, it is initialized in the algorithm by randomly generating: where X is the grey wolf population, i ∈ [1, 2, 3 . . . , N] and j ∈ [1, 2, 3 . . . , sizepop], N is the number of grey wolf populations, sizepop is the population dimension; lb and ub are the lower bound and upper bound of the search interval respectively; and U is a random uniform distribution function.

Encircling Prey
As mentioned above, grey wolves encircle prey during the hunt. In order to mathematically model encircling behavior, the following equations are proposed: where → D is the distance between the prey and grey wolf; t indicates the current iteration,

Population Location Update
By calculating the objective function value of the GWO algorithm, the best solution, the second best solution, and the third best solution are set to α wolf, β wolf, and δ wolf. The position of other grey wolves is determined by the position of α wolf, β wolf and δ wolf, as shown in Equations (6)- (8). After generating a new group, boundary control is performed on the elements in the population to complete an iteration. The above process is repeated until the algorithm termination condition is satisfied, and finally the optimal solution is outputted.

Lévy-Embedded Grey Wolf Optimization Algorithm
The GWO algorithm has the advantages of not relying on parameter settings and convenient implementation, however the GWO algorithm is still prone to prematurely falling into local extremum when solving complex optimization problems; that is, the phenomenon of premature convergence occurs. In order to alleviate the above problems, Heidari proposed a Lévy-embedded Grey Wolf Optimization (LGWO) algorithm based on the GWO algorithm combined with Lévy Flight [7]. The Lévy Flight can assist the GWO in searching based on deeper searching patterns. Using this concept, it can be ensured that the GWO can handle global searching more efficiently. The LGWO algorithm mainly improves the GWO algorithm in three aspects.

1.
The role of δ wolves in the social hierarchy will be played by other wolves.

2.
The Lévy Flight concept is embedded into the modified GWO.

3.
The greedy selection (GS) strategy is employed for the LGWO.
The specific process of the LGWO algorithm is summarized as follows, and the pseudo code of the LGWO algorithm is described in detail, as shown in Table 1. 1.
Population initialization. The process is the same as Equation (1)  Population search. The process is the same as Equations (2) and (3) of the GWO algorithm in Section 2.2.2.

3.
By calculating the objective function value of the LGWO algorithm, the best solution α wolf and second-best solution β wolf are obtained. The positions of other grey wolves are determined by the position of α wolf and β wolf, as shown in Equation (9): where where u and v follow a normal distribution: In the LGWO Algorithm 1 proposed by Heidari, the parameter β is a random number of [0,2].
Algorithm 1 Pseudo code for the LGWO algorithm 1.
Determine the initial swarm size N and number of iterations T 2.
Randomly generate the initial population of wolves 3. Initialize a, p, Compute the fitness of each wolf

5.
Set to → X α be the best wolf 6. Set to → X β be the second best wolf 7.
For each wolf 9.
Update the position of current wolves by Equation (13)  After generating a new group, Heidari [7] selected whether to retain the updated grey wolf through the greedy selection strategy depicted by Equation (11), and completes an iteration. The above process is repeated until the algorithm termination condition is satisfied, and finally the optimal solution is outputted.
where r new and p are random values inside (0, 1).
Additionally, referring to the greedy selection (GS) strategy from the DE algorithm, the concept of "survival of the fittest" is utilized here with the probability p. According to this strategy, new superior positions in each generation can continue to be more enriched for the next generations and the worse ones are disregarded, so the LGWO algorithm has a stronger search capability.

Virtual Force Algorithm
The Virtual Force Algorithm (VFA) was originally an algorithm that allowed mobile robots to evade obstacles in an unknown environment [19][20][21]. The VFA has been applied to WSN coverage optimization, and can quickly disperse sensor nodes. The VFA treats the sensor nodes as charged particles. The forces between the sensor nodes are presented by attractive and repulsive patterns. When two nodes are close enough, smaller than the threshold D th , the force is in a repulsive pattern, which intends to separate them. When the distance between two nodes is larger than the threshold D th , the force becomes an attractive pattern, which draws them closer. Assuming that the set of sensor nodes in the region is S = {s 1 , s 2 , s 3 , . . . , s l }, then the distance d ij between the i, j sensor nodes is s i (x i , y i ) and s j (x j , y j ): Then, the force of the sensor node s j on the sensor node s i is: where ω α and ω γ are the virtual force coefficients, ω γ ω α ; D th is the threshold for generating gravitation and repulsive force, and the value of D th is √ 3R s [20]; and α ij is the angle of the line between sensor node i and sensor node j with respect to the y-axis. If d ij < D th , it is considered that sensor node i is too close to sensor node j, and the force is in a repulsive pattern; if d ij > D th , the force becomes an attractive pattern. The resultant force of the sensor node i is → F i , the resultant force in the horizontal direction is → F ix , the resultant force in the vertical direction is → F iy , and the position of the node is updated by Equations (14) and (15).
where x(k) and y(k) are the horizontal and vertical coordinates of sensor node s i at time k, respectively; and MaxStep is the step size for the wireless sensor node to move.

Improved Virtual Force Algorithm
Generally, the communication radius R of the sensor node is defined as twice the perceived radius R s . Actually, when d ij is larger than R, they cannot exchange information with each other, they cannot calculate or obtain the distance, and the VFA cannot be used in this case. The purpose of applying the improved virtual force algorithm in this paper is to apply it to wireless sensor network node deployment, and specifically to apply it to wireless sensor network nodes, which not only have sensor sensing capabilities, but also have communication capabilities. Therefore, the distance between sensor nodes must be less than the communication radius, and the sensor node distribution should not be too dense.
In the improved virtual force algorithm applied in this paper, when the d ij is greater than R, there is no virtual force between the sensor nodes. Equation (13) of the sensor node s j to the sensor node s i is modified: The physical meaning of the variables and parameters in Equation (16) is the same as Equation (13), and the updated equation of the sensor node position is the same as Equations (14) and (15).

Problem Description
All sensor nodes are mobile.

2.
All sensor nodes are omnidirectional sensors, and the perceptual model is a probabilistic detection model. 3.
All sensor node location information is known by GPS. 4.
All sensor nodes are able to move to the scheduled position, where the scheduled position is within their mobility range.

Coverage Rate
The sensor node in this paper adopts the probabilistic detection model. Compared with the disk model, this model can objectively reflect the real network deployment environment, as shown in Figure 1. In Figure 1, the inside of the solid circle is the determined area, and the area between the solid circle and the outermost dotted circle is an uncertain area. The outside of the outermost dotted circle is an unmeasured area. Assuming that the monitoring area is a two-dimensional planar area, each sensor node can be deployed anywhere within the area. If the coordinates of the sensor node s i in the sensor network are (x i , y i ) and the coordinates of the monitoring target point p (pixel point) in the region are (x p , y p ), then the monitoring probability C p (s i , p) of the sensor node s i to the target point p is set as follows: where d(s i ,p) is the Euclidean distance between the sensor node s i and the target point p; R s is the distance radius that the sensor node can perceive; r e (0 < r e < R s ) is the measurement reliability parameter of the sensor node; α 1 = r e − R s + d(s i ,p), α 2 = r e + R s − d(s i ,p), and α 1 > 0, α 2 > 0; λ 1 , λ 2 , β 1 , β 2 are the measurement parameters related to the characteristics of the sensor node. The parameters of the work by Author [22] are used to set λ 1 = 1, λ 2 = 0, β 1 = 1, β 2 = 1.5.

Problem Description
3.1.1. Network Model 1. All sensor nodes are mobile. 2. All sensor nodes are omnidirectional sensors, and the perceptual model is a probabilistic detection model. 3. All sensor node location information is known by GPS. 4. All sensor nodes are able to move to the scheduled position, where the scheduled position is within their mobility range.

Coverage Rate
The sensor node in this paper adopts the probabilistic detection model. Compared with the disk model, this model can objectively reflect the real network deployment environment, as shown in Figure 1. In Figure 1, the inside of the solid circle is the determined area, and the area between the solid circle and the outermost dotted circle is an uncertain area. The outside of the outermost dotted circle is an unmeasured area. Assuming that the monitoring area is a two-dimensional planar area, each sensor node can be deployed anywhere within the area. If the coordinates of the sensor node si in the sensor network are (xi, yi) and the coordinates of the monitoring target point p (pixel point) in the region are (xp, yp), then the monitoring probability Cp (si, p) of the sensor node si to the target point p is set as follows: where d(si,p) is the Euclidean distance between the sensor node si and the target point p; Rs is the distance radius that the sensor node can perceive; re(0 < re < Rs) is the measurement reliability parameter of the sensor node; α1 = re − Rs + d(si,p), α2 = re + Rs − d(si,p), and α1 > 0, α2 > 0; λ1, λ2, β1, β2 are the measurement parameters related to the characteristics of the sensor node. The parameters of the work by Author [22] are used to set λ1 = 1, λ2 = 0, β1 = 1, β2 = 1.5. Let the probability threshold of target point p be detected as Cth, and the joint monitoring probability Cp(sall, p) of the sensor node to target point p in the whole monitoring area is: Let the probability threshold of target point p be detected as C th , and the joint monitoring probability C p (s all , p) of the sensor node to target point p in the whole monitoring area is: where s all represents the set of sensor nodes that measure target point p, and l is the number of sensor nodes in the monitoring area.
Assuming that the monitoring area is a rectangle of m × n (m 2 ), the area to be measured is divided into m × n identical grids of equal size and area of 1 m 2 . Then, the grids are simplified to pixel points with a discrete precision of 1. In this paper, the WSN node coverage rate is defined as the ratio of the number of grids covered by Equation (19) to the total number of grids in the monitoring area, namely: Then, the coverage optimization problem is briefly described as follows: Step 1. Calculate the probability of monitoring one-pixel point by one sensor node using Equation (17).
Step 2. Calculate the joint monitoring probability of all sensor nodes to one-pixel point using Equation (18).
Step 3. The coverage rate of the monitored area is calculated using Equation (19). Equation (19) is used as the objective function of the coverage optimization algorithm.

Uniformity of Sensor Node Distribution
Uniformity is an important indicator for measuring the uniform distribution of sensor nodes in wireless sensor networks. It is used to describe the uniformity of the sensor node distribution. It is usually expressed by the mean of the sum of the standard deviation of all sensor nodes. The wireless sensor network node uniformity indicator E is calculated as shown in Equation (20): where l is the total number of sensor nodes, k is the number of neighboring nodes of sensor node i, D i,j is the Euclidean distance between adjacent nodes i, j, and M i is the average of the Euclidean distance between the sensor node i and all its neighboring sensor nodes. It can be found from Equation (20) that the closer the values of D i,j and M i are, the more uniform the distribution is. Therefore, the smaller the E value, the better the uniformity of the sensor node distribution.

Average Moving Distance of Sensor Nodes
The average moving distance of sensor nodes refers to the average distance that each sensor node in the wireless sensor network moves from the initial position to the best position with the largest final coverage rate. It is assumed that the initial position information of each sensor node is known. Through the MATLAB simulation, the VFLGWO algorithm is used to solve the optimal position of each sensor node with the largest coverage rate. Each sensor node is directly moved from the initial position to the optimal position in a straight line. Obviously, for wireless sensor nodes, the shorter the average moving distance is, the smaller the energy consumption is.

The Modified LGWO Algorithm
The LGWO algorithm applies GS strategy to allow the algorithm to retain a better solution for each iteration. For the LGWO algorithm, f ( → X) represents a convergence curve and the smaller the value, the better. However, for wireless sensor node deployment, f ( → X) represents f itness( → X) (the coverage rate). The larger the f itness( → X) value that is obtained, the larger the corresponding coverage rate. Therefore, the modified LGWO proposed in this paper is the same as the LGWO, except that Equation (11) is modified to Equation (21). When the new value f itness( → X) obtained by a new iteration is not greater than the previous generation value and the satisfaction probability condition is reached, the optimal solution is not updated; otherwise, the obtained superior solution is substituted for the previous generation solution.

Combination of the Improved Virtual Force Algorithm and the Modified LGWO Algorithm
The performance of the LGWO is significantly better than that of the GWO. However, when the modified LGWO algorithm is applied to WSN node deployment, the problem of non-uniform distribution of sensor nodes still exists, while the Virtual Force algorithm can make the sensor nodes more uniform. Therefore, combined with the modified LGWO algorithm, the improved Virtual Force algorithm is introduced. In this paper, the Virtual Force-Lévy-embedded Grey Wolf Optimization (VFLGWO) algorithm is proposed.
In the modified LGWO algorithm, the position of a wolf represents a possible solution to the optimization problem. Therefore, the deployment of the sensor nodes in the detection area refers to a wolf (a solution) in the algorithm.
The procedure of deploying wireless sensor nodes by the VFLGWO algorithm can be summarized as follows: N sensor nodes are randomly generated to be the initial position of the real sensor nodes. Initialize each wolf (each set of solutions) by the modified LGWO algorithm as the virtual sensor node position.
Enter the iteration, update the position of each wolf using Equations (9) and (21). Then, apply the improved Virtual Force algorithm to update the sensor nodes' positions for each wolf, calculate the virtual force using Equation (20), and adjust the position of the nodes for each wolf using Equations (14) and (15). Apply the survival of the fittest selection rule using Equation (22), and retain the better solution obtained by the improved virtual force algorithm: Through the superior solution of the t th generation → X (t), whether to update α wolf or β wolf is determined; the VFLGWO algorithm completes an iteration. If the number of iterations does not reach the specified number of iterations T, perform the above iterative process again; otherwise, the process ends. The specific implementation of the VFLGWO algorithm will be described in Section 3.3.

Procedure for the VFLGWO Algorithm
Application of the VFLGWO algorithm in WSN node deployment is shown in Figure 2. The steps of the VFLGWO algorithm are described as follows: 1.
Determine the size of the monitoring area and the number of sensor nodes N to be deployed, and randomly deploy the sensor nodes in the monitoring area at the real sensor node positions, the set of which is S 0 = {s 1 , s 2 , s 3 , . . . , s N }.

3.
Generate the initialized grey wolf population → X(t) randomly and calculate the coverage rate fitness t of each wolf (each set of solutions) using Equations (17)- (19).

4.
Set the solution with the largest coverage rate to α wolf, set the solution with the second largest coverage rate to β wolf, and record the coverage rate fitness (α) and fitness (β) of α wolf and β wolf.

5.
Enter the iteration, set t = 1, calculate the position of each grey wolf after the tth iteration using Equation (11), calculate the coverage fitness t of each wolf, judge the better grey wolf using Equation (21), and then update a, p, Apply the improved virtual force algorithm to calculate the virtual force of each sensor node of each wolf using Equation (16), adjust the position of each sensor node of each wolf using Equations (14) and (15), calculate the coverage fitness t ' of each wolf, and obtain better grey wolves → X (t) using Equation (22). 7.
Compare the fitnesst' value of the tth iteration with the fitness(α) and fitness(β) obtained by the previous iteration to determine whether to update the α wolf or the β wolf. 8.
Determine whether the number of iterations reaches the specified number of iterations T, "No" then t = t + 1, return 5); "Yes" then the iteration ends. 9.
Output the optimal solution α wolf as the optimal target position of the real sensor nodes, recorded 5. Enter the iteration, set t = 1, calculate the position of each grey wolf after the tth iteration using Equation (11), calculate the coverage fitnesst of each wolf, judge the better grey wolf using Equation (21), and then update a, p, ， .
6. Apply the improved virtual force algorithm to calculate the virtual force of each sensor node of each wolf using Equation (16), adjust the position of each sensor node of each wolf using Equations (14) and (15), calculate the coverage fitnesst' of each wolf, and obtain better grey wolves ′( ) using Equation (22). 7. Compare the fitnesst' value of the tth iteration with the fitness(α) and fitness(β) obtained by the previous iteration to determine whether to update the α wolf or the β wolf. 8. Determine whether the number of iterations reaches the specified number of iterations T, "No" then t = t + 1, return 5); "Yes" then the iteration ends. 9. Output the optimal solution α wolf as the optimal target position of the real sensor nodes, recorded as 0  There are currently two options for WSN node deployment: (1) give the target sensor node position without considering how the sensor nodes move; (2) match the same sequence numbers of the sensor nodes before and after optimization and design the corresponding relationship of node movements.
The pseudo code of the sensor nodes matching algorithm proposed in this paper is shown in Algorithm 2. The initial position set of N sensor nodes in the region is S 0 = {s 1 , s 2 , s 3 , . . . , s N }, and the target position for obtaining the optimal solution is S 0 = {S 1 , S 2 , S 3 , . . . , S N }. The matching of S 0 and S 0 is achieved by the sensor nodes matching algorithm. Algorithm 2 Pseudo code for sensor nodes matching algorithm 1.
Input initial position set S 0 , optimal solution position set S 0 , number of nodes N 2.
For i from 1 to N do 3.
For j from 1 to N do 4.
Calculate the Euclidean distance between the S 0 node position s i and the S 0 node position s j and record 5. The VFLGWO algorithm proposed in this paper applies the LGWO algorithm with a stronger search ability to the deployment of wireless sensor nodes, and combines the improved Virtual Force algorithm to achieve a wireless sensor network with a higher coverage rate and more uniform distribution. In order to achieve the goal of making the average moving distance of the sensor nodes shorter, the randomly generated real node positions are taken as the initial positions of the sensor nodes, and the sensor node position of the optimal solution obtained by the VFLGWO algorithm is taken as the final position of the real node. The sensor nodes matching algorithm proposed in this paper is used to complete the WSN node deployment. The VFLGWO algorithm proposed in this paper executes the improved Virtual Force algorithm again after the LGWO search phase. The time complexity of the improved Virtual Force algorithm is O (N 2 ). Therefore, according to the VFLGWO algorithm process, it needs to be updated to N times per iteration. The time complexity is O (m × N 2 ). In addition, the T set in this paper is much larger than N. Therefore, the time complexity of the sensor nodes matching algorithm proposed in this paper is negligible. The time complexity of the total coverage strategy of the VFLGWO algorithm is O (T ×m × N 2 ).

Algorithm Time Complexity Analysis
Although the VFLGWO algorithm increases the amount of computation, the coverage rate and uniformity obtained by the VFLGWO algorithm are significantly higher than other algorithms. The position adjustment of the sensor nodes after random delivery is not very time-sensitive. The position-adjusted coverage rate is significant, therefore it makes sense to increase the amount of computation to achieve higher wireless sensor network deployment performance.

VFLGWO Algorithm Simulation
Firstly, the relevant parameters of the VFLGWO algorithm are discussed through simulation experiments. Different simulation experiments are designed to test the performance of the VFLGWO algorithm. The performance indexes of the simulation test in this paper include: coverage rate, uniformity, average moving distance of sensor nodes (m), and running time (s). When testing the average moving distance of sensor nodes, the VFLGWO algorithm applies the sensor nodes matching algorithm proposed in this paper. Other comparison algorithms still use the conventional method in which the S 0 and S 0 are matched by the sequence number of sensor nodes.
In this paper, the experiments are performed using MATLAB 2014. The simulation environment is an area of 50 m × 50 m. The number of randomly distributed mobile sensor nodes is N = 50, and the sensing radius is R s = 5 m for all sensor nodes; the communication radius is R = 2R s = 10 m, the detection probability threshold is C th = 0.8, the measurement reliability parameter is r e = 2.5 m, the number of basic iterations is T = 3000 generations, and the population number is sizepop = 30 (the algorithm uses the cell array to store the sensor node position coordinates, and each wolf contains 50 sensor node position coordinates).

Determination of Gravity and Repulsion Parameters
According to the gravity and repulsion parameters of the VFLGWO algorithm, six groups of experiments are used to obtain relatively optimal ω α and ω γ . Each group of experiments is performed 20 times independently, and the average values of each indicator are obtained. The control variable method is applied. Firstly, the MaxStep is unified at 1.2 m, and the experimental value of ω α is unchanged at 1. The values of ω γ are 600, 800, 1000, 1500 and 2000, respectively. Since this article mainly discusses how to deploy sensor nodes to maximize the coverage rate of the monitoring area, coverage rate is taken as the first indicator. The four optimization indicators calculated by the VFLGWO algorithm are listed in Table 1. As can be seen from Table 1, when the value of ω γ is 1000, the coverage rate is the largest. Therefore, in the simulation experiment of the proposed algorithm, ω α is set to 1 and ω γ is set to 1000.  The experimental results are shown in Table 2. When the MaxStep is 1.2 m, the coverage rate is the largest. Therefore, set the MaxStep to 1.2 m.

VFLGWO Algorithm Effectiveness Simulation Experiment
In order to verify the performance of the VFLGWO algorithm, the VFLGWO is compared with the wireless sensor network coverage optimization algorithm using Chaotic Particle Swarm (CPSO), Cuckoo Search (CS), and LGWO.

Comparison Algorithm Parameter Settings
The parameters of the CPSO, CS, and LGWO are set as shown in Table 3 [7]. The other parameters, such as the size of the population sizepop and the number of basic iterations T, are the same as those of the VFLGWO algorithm.  Figure 3 shows the initial positions of the sensor nodes. Figure 4a-d are the development of the best solutions (each large circle represents a schematic representation of the sensing area of the sensor nodes at its center, with a perceived radius R s = 5 m) after one optimization by the CPSO, CS, LGWO, and VFLGWO algorithms. The coverage rate and uniformity of the CPSO, CS, LGWO, and VFLGWO algorithms for this experiment are listed in Table 4. It can be seen from Figure 4 that the sensor nodes have the most uniform distribution and the largest coverage rate after the VFLGWO algorithm is optimized. The experimental data in Table 4 also confirmed this conclusion.

Multiple Performance Index of Four Algorithms
In this paper, 20 independent experiments of 3000 generations are performed on the four algorithms. The average values of the coverage, uniformity, average moving distance of the sensor nodes, and running time are listed in Table 5. In order to facilitate visual comparison, the data of Table 5 is plotted as a histogram, as shown in Figure 5. It can be seen from Table 5 that the coverage rate of the VFLGWO is improved by 4.05%, 4.87% and 9.94% compared with the other three

Multiple Performance Index of Four Algorithms
In this paper, 20 independent experiments of 3000 generations are performed on the four algorithms. The average values of the coverage, uniformity, average moving distance of the sensor nodes, and running time are listed in Table 5. In order to facilitate visual comparison, the data of Table 5 is plotted as a histogram, as shown in Figure 5. It can be seen from Table 5 that the coverage rate of the VFLGWO is improved by 4.05%, 4.87% and 9.94% compared with the other three

Multiple Performance Index of Four Algorithms
In this paper, 20 independent experiments of 3000 generations are performed on the four algorithms. The average values of the coverage, uniformity, average moving distance of the sensor nodes, and running time are listed in Table 5. In order to facilitate visual comparison, the data of Table 5 is plotted as a histogram, as shown in Figure 5. It can be seen from Table 5 that the coverage rate of the VFLGWO is improved by 4.05%, 4.87% and 9.94% compared with the other three algorithms, respectively; the uniformity of the VFLGWO algorithm is optimal. The average moving distance of the sensor nodes of the VFLGWO algorithm is significantly shorter than the average moving distance of the sensor nodes of the other three algorithms. In terms of running time of the VFLGWO algorithm, it is about twice that of the other algorithms. algorithms, respectively; the uniformity of the VFLGWO algorithm is optimal. The average moving distance of the sensor nodes of the VFLGWO algorithm is significantly shorter than the average moving distance of the sensor nodes of the other three algorithms. In terms of running time of the VFLGWO algorithm, it is about twice that of the other algorithms.  5. Comparison of the four algorithms. Figure 6 shows the network coverage rate curves of the four algorithms, CPSO, CS, LGWO and VFLGWO; running 3000 generations. It can be seen from Figure 6 that the coverage rate of the VFLGWO algorithm is significantly larger than that of the other three algorithms.   Figure 6 shows the network coverage rate curves of the four algorithms, CPSO, CS, LGWO and VFLGWO; running 3000 generations. It can be seen from Figure 6 that the coverage rate of the VFLGWO algorithm is significantly larger than that of the other three algorithms. algorithms, respectively; the uniformity of the VFLGWO algorithm is optimal. The average moving distance of the sensor nodes of the VFLGWO algorithm is significantly shorter than the average moving distance of the sensor nodes of the other three algorithms. In terms of running time of the VFLGWO algorithm, it is about twice that of the other algorithms.  5. Comparison of the four algorithms. Figure 6 shows the network coverage rate curves of the four algorithms, CPSO, CS, LGWO and VFLGWO; running 3000 generations. It can be seen from Figure 6 that the coverage rate of the VFLGWO algorithm is significantly larger than that of the other three algorithms.

Different Numbers of Sensors in the Same Area
The performance index of the coverage rate, uniformity, and average moving distance of the number of sensor nodes placed in the 50 m × 50 m monitoring area from 40 to 60 (each group adds five sensor nodes) are listed in Tables 6-8. In order to facilitate visual comparison, the data of Tables 6-8 are plotted as graphs, as shown in Figure 7. It can be seen from Figure 7 that when the size of the monitoring area is the same, the coverage rate of the four algorithms increases as the number of sensor nodes increases. This is because the number of sensor nodes per unit area in the monitoring area increases, and the area coverage rate naturally increases. At the same time, the coverage rate and uniformity of the VFLGWO algorithm are better than the other three algorithms, and the average moving distance of the sensor nodes of the VFLGWO algorithm is significantly shorter than the other three algorithms. The performance index of the coverage rate, uniformity, and average moving distance of the number of sensor nodes placed in the 50 m × 50 m monitoring area from 40 to 60 (each group adds five sensor nodes) are listed in Tables 6-8. In order to facilitate visual comparison, the data of Table  6-8 are plotted as graphs, as shown in Figure 7. It can be seen from Figure 7 that when the size of the monitoring area is the same, the coverage rate of the four algorithms increases as the number of sensor nodes increases. This is because the number of sensor nodes per unit area in the monitoring area increases, and the area coverage rate naturally increases. At the same time, the coverage rate and uniformity of the VFLGWO algorithm are better than the other three algorithms, and the average moving distance of the sensor nodes of the VFLGWO algorithm is significantly shorter than the other three algorithms.   Tables 9-11. In order to facilitate visual comparison, the data of Tables 9-11 are plotted as histograms, as shown in Figure 8. It can be seen from Figure 8 that the VFLGWO algorithm is superior to the other three algorithms in coverage rate, uniformity, and average moving distance of the sensor nodes in different size monitoring areas.   Tables 9-11. In order to facilitate visual comparison, the data of Tables 9-11 are plotted as histograms, as shown in Figure 8. It can be seen from Figure 8 that the VFLGWO algorithm is superior to the other three algorithms in coverage rate, uniformity, and average moving distance of the sensor nodes in different size monitoring areas.  In summary, coverage rate, uniformity, and average moving distance of sensor nodes of the VFLGWO algorithm achieve a better performance compared to the other three algorithms. The sensor nodes adjusted by the VFLGWO algorithm are more uniformly distributed and have shorter sensor node movement distances. This paper also performs the indicators of the VFLGWO algorithm when the working environment changes. The experimental results of the VFLGWO algorithm are still better than the experimental results of the CPSO, CS and LGWO algorithms.

Conclusions
In this paper, the modified LGWO algorithm is applied to WSN node deployment. In order to obtain a better performance index, a coverage optimization algorithm for wireless sensor networks with the Virtual Force-Lévy-embedded Grey Wolf Optimization (VFLGWO) algorithm is proposed to solve the coverage optimization problem. The VFLGWO algorithm can effectively make sensor node distribution more uniform, the coverage rate higher in the same environment, the average moving distance of the sensor nodes shorter, and it can ensure the expected design goal is obtained. The simulation results show that the wireless sensor network coverage rate, uniformity, and average moving distance of the VFLGWO algorithm are significantly better than those of the CPSO, CS and LGWO, and the VFLGWO algorithm has good environmental adaptability.
After obtaining the optimal position of the wireless sensor nodes through the VFLGWO algorithm, the sensor nodes matching algorithm proposed in this paper can greatly shorten the average moving distance of the sensor nodes. The running time of the VFLGWO algorithm is relatively long, however because there is no need for frequent scheduling during the deployment of wireless sensor nodes, it will not affect the deployment of WSN nodes-and the coverage rate, In summary, coverage rate, uniformity, and average moving distance of sensor nodes of the VFLGWO algorithm achieve a better performance compared to the other three algorithms. The sensor nodes adjusted by the VFLGWO algorithm are more uniformly distributed and have shorter sensor node movement distances. This paper also performs the indicators of the VFLGWO algorithm when the working environment changes. The experimental results of the VFLGWO algorithm are still better than the experimental results of the CPSO, CS and LGWO algorithms.

Conclusions
In this paper, the modified LGWO algorithm is applied to WSN node deployment. In order to obtain a better performance index, a coverage optimization algorithm for wireless sensor networks with the Virtual Force-Lévy-embedded Grey Wolf Optimization (VFLGWO) algorithm is proposed to solve the coverage optimization problem. The VFLGWO algorithm can effectively make sensor node distribution more uniform, the coverage rate higher in the same environment, the average moving distance of the sensor nodes shorter, and it can ensure the expected design goal is obtained. The simulation results show that the wireless sensor network coverage rate, uniformity, and average moving distance of the VFLGWO algorithm are significantly better than those of the CPSO, CS and LGWO, and the VFLGWO algorithm has good environmental adaptability.
After obtaining the optimal position of the wireless sensor nodes through the VFLGWO algorithm, the sensor nodes matching algorithm proposed in this paper can greatly shorten the average moving distance of the sensor nodes. The running time of the VFLGWO algorithm is relatively long, however because there is no need for frequent scheduling during the deployment of wireless sensor nodes, it will not affect the deployment of WSN nodes-and the coverage rate, uniformity, and average