Research on Wireless Sensor Network Coverage Path Optimization Based on Biogeography-Based Optimization Algorithm

. Path selection is one of the key technologies of wireless sensor network (WSN). A reasonable choice of coverage path can improve the service quality of WSN and extend the life cycle of WSN. Biogeography-based optimization (BBO) is widely used in the ﬁeld of cluster intelligent optimization because its search method has a better incentive mechanism for population evolution. In this paper, the move-in and move-out operation and mutation operation of the BBO algorithm enable WSN to ﬁnd an eﬃcient routing path. In this paper, simulation experiments are carried out in two scenarios of regular deployment and random deployment of WSN nodes. The experimental results show that the quality of the WSN coverage path solution optimized by the BBO algorithm in the two scenarios is better than that of the particle swarm algorithm and genetic algorithm.


Introduction
With the gradual deepening research of quality of service (QoS) issues in WSN, we are paying more attention to the coverage path optimization of wireless sensor network [1,2].As the energy of sensor nodes in WSN is limited and cannot be maintained in time, establishing an efficient path connection method can not only effectively save network energy on the one hand but also improve the quality of service of WSN on the other hand.e main goal of WSN coverage path optimization is to establish an efficient and reliable path from the sensor sending node to the receiving node and maximize the network life of the WSN as much as possible [3,4].
In this paper, through the research and improvement of the biogeography-based optimization algorithm, it will make full use of the ability of the biogeography-based optimization algorithm to sense interactive data in multidimensional and high-dimensional problems to achieve the optimal construction of the wireless sensor network coverage path.e main contributions of this article are as follows: (1) Analyze the main problems of wireless sensor network coverage path optimization, explain the impact of coverage path optimization on wireless sensor network deployment performance, and give a specific optimization model.(2) Explain in detail the main principles, mathematical model, migration, and mutation operation steps of the biogeography-based optimization algorithm, and apply it to the WSN coverage path optimization problem.rough the comparison with the coverage path optimization experiments of the genetic algorithm and particle swarm optimization (PSO) algorithm, the effectiveness of the biogeography-based optimization algorithm for wireless sensor network coverage path optimization is further verified, and the local convergence problem in the iterative optimization process is solved and also the robustness of the communication link after WSN networking is enhanced.
e rest of this paper is organized as follows.Section 2 explains the structure and principle of the WSN coverage path optimization model.Section 3 explains the main mechanism and principle of the BBO algorithm.Section 4 gives the design of the WSN coverage path optimization algorithm based on the BBO algorithm.Section 5 explains the simulation experiments conducted under two WSN node deployment conditions to verify the effectiveness of the BBO algorithm in the WSN coverage path optimization problem.

WSN Coverage Path Optimization Model
In WSN structure, sensor nodes are divided into source node, sink node, and target node.
e sensor node that collects monitoring data is the source node, whose data are summarized to the sink node through a certain multihop path and then sent to the target node (or user) through the Internet [5][6][7][8][9].In order to improve the efficiency of data transmission and balance the energy consumption of the entire network, all sensor nodes of WSN are grouped [10][11][12][13].Each group of sensor nodes is called a cluster.A cluster head is selected from the cluster as the sink node of the cluster to process the data of the cluster and communicate with other clusters and users.e cluster head node of each cluster is the target node on the communication path of the cluster, which is not fixed, but dynamically adjusted according to some energy consumption mechanism [14][15][16][17].As the energy of each node will gradually decrease as the time goes on, the coverage path in each cluster of WSN needs to be adjusted timely according to the energy situation of the whole network.Because of the complicated topological structure of WSN, the plane structure of WSN will be simplified in order to facilitate research without losing the effectiveness of the algorithm [18,19].e WSN cluster is simplified to a rectangular area with a known size in a twodimensional space.e number and location of sensor nodes in the area are known [20][21][22].e sensor nodes are homogeneous nodes (the sensing radius, communication radius, and initial energy are the same), and these nodes are treated as a point within the region and are numbered in an orderly manner (1, 2, . . ., n). erefore, an undirected graph can be used to describe the WSN system structure.e undirected weight graph of any cluster G � 〈V, E〉, where is the node set in the cluster, V Hea d is the cluster head node of G, and V i (i � 1, 2, . . ., n) is other sensor nodes in G.
e communication link set between nodes is denoted by E � (e 1 , e 2 , . . ., e m ), which can be written as follows: where λ 0 is the effective transmission distance between each sensor node.e energy consumption set of the link between nodes j and k is C(j, k) � C(k, j). e judgment value of whether the communication link between node j and node k is connected is as follows: In summary, the mathematical model of WSN coverage path optimization is as follows: Equation ( 4) can ensure the shortest path between the source node and the target node [23,24].
e node selection mechanism of WSN coverage path optimization problem is different from the simple routing node selection of traditional wireless network [25][26][27].e limit energy of sensor nodes in WSN, the more difficult positioning of the nodes, and the uncertainty of the monitoring area environment make it necessary for the WSN to take into account a lot of factors when selecting the path, thus using more complex routing mechanisms [28,29].
Aiming at this more complex routing mechanism, this paper uses the BBO algorithm with an adaptive mechanism to optimize the WSN coverage path.

BBO Algorithm
e BBO algorithm is also an intelligent algorithm based on biological populations, which is similar to PSO.Biological populations live in different habitats, and each habitat has its own fitness index (HSI) [30,31].Factors affecting HSI include rainfall, vegetation diversity, and climate conditions, which form a set of fitness vectors (SIV) [32,33].e population distribution and individual migration of habitats are greatly affected by HSI.Habitats with higher HSI hold more populations, and the ecological load tends to be saturated, so the emigration rate is high and the immigration rate is low [34,35].Habitats with lower HSI hold fewer populations, so the emigration rate is low and the emigration rate is high.Of course, as the population moves out and the biodiversity of the habitat declines, HSI will also decrease.In the same way, as populations move in and the biodiversity of the habitat increases, HSI will also increase.In addition, it is also necessary to consider that organisms will have a certain probability of mutation, which will also affect HSI.
In the BBO algorithm, the habitat suitability vector SIV corresponds to the solution of the optimization problem [36,37].Habitat migration and migration operations correspond to the information sharing of the algorithm.e program with higher HSI operates according to its own migration rate and shares the information with the program with lower HSI.e program with lower HSI starts from higher.e HSI scheme has received many new problemsolving factors; the existence of mutation factors can help the 2 Complexity diversity of the population to a certain extent, which increases the search ability of the algorithm in the optimization problem space.

Migration Model.
e migration model of a certain habitat population in biogeography is shown in Figure 1. e abscissa represents the number of populations in the habitat, and the ordinate represents the rate of habitat migration.e emigration rate and immigration rate of the population are denoted by μ(s) and λ(s) respectively.
Assume that μ(s) and λ(s) are the immigration rate and the emigration rate when the number of habitat populations is s.During the time period from t to t + Δt, the probability that the habitat accommodates s populations at the same time is as follows.
Suppose μ(s) and λ(s) are the emigration rate and the emigration rate when the number of habitat population is s, respectively.In the time period from t to t + Δt, the probability that the habitat can accommodate s populations at the same time is as follows: In order for the habitat to have s populations at any time t + Δt, the following three conditions must be met: (1) At time t, there are s populations in the habitat, and no populations will be migrated in or out during the period from t to t + Δt (2) At time t, there are s − 1 populations in the habitat, and only one population migrates in the time period from t to t + Δt (3) At time t, there are s + 1 populations in the habitat, and only one population migrates out from t to t + Δt For the convenience of calculation, it is assumed that Δt is so small that the migration in or out of the population within the time period can be ignored.en, let Δt ⟶ 0 and find the limit of equation ( 5) as follows: In which, n � S max , P � [P 0 , . . ., P n ] T , and we can get the matrix form of _ P s as follows: In equation (7), the matrix A is as follows: According to the migration model function shown in Figure 1, if a given population size is k, the emigration rate and immigration rate corresponding to this population size can be calculated as follows: Considering the special case when E � I, In [3], the author gives the calculation of special cases.When the eigenvalue of A ′ is 0, the corresponding eigenvector is as follows:

Complexity
in which ceil ( * ) is the function that takes the upper limit integer of the formal parameter.At the same time, it can be also inferred that when the maximum population number in the habitat is n, the vector of probability composition corresponding to different population numbers is as follows: 3.2.Migration Operation.Migration operation is an important part of the BBO algorithm.rough the mutual migration of populations between various habitats, the algorithm makes a wide-area search in the solution space [38].
is feature is similar to the genetic strategy in the genetic algorithm.ere is a mapping relationship between HSI and population size in the BBO algorithm.Before the migration operation is performed, all habitats are sorted according to their HSI values.s 1 represents the habitat with lower HSI, and s 2 represents the habitat with higher HSI; then, λ(s 1 ) > λ(s 2 ) and μ(s 1 ) < μ(s 2 ).In the process of migration operation, the probability of habitat change is proportional to its migration rate and the probability of change of the habitat of the new population is proportional to its migration rate.
It is assumed that the BBO algorithm is consisted of n habitats, each habitat has a D-dimensional fitness variable, and the vector x i � (x i1 , x i2 , . . ., x i D ) represents the potential solutions in the D-dimensional search space, i � (1, 2, . . ., n).First, the algorithm determines whether habitat i is selected for the introduction of new populations according to the global migration rate P mod ∈ [0, 1].Under a certain probability, if the habitat is selected, the suitability variable x ij of habitat i, (j � 1, 2, . . ., D) will be determined whether to be changed or not according to its immigration rate λ(s i ).If the suitability variable of habitat i is selected, the selection is made according to the emigration rate μ(s k ) of other habitats k (k ≠ i).If habitat k is selected as the immigration source of habitat i, then the suitability variable x kl , (l � 1, 2, . . ., D) of habitat k will replace the suitability variable x ij of habitat i.

Mutation Operation.
e mutation operation ensures that the potential solutions of the BBO algorithm have a certain diversity and will not be degenerated in the iterative process.e main problem of mutation operation lies in how to give the corresponding mutation probability according to the existing population number in the habitat.e probability of the population quantity in the habitat has been given by equation ( 6), which represents the pre-existing possibility of the solution to the problem.If the population probability of a habitat is low, the probability of an optimal solution is low too, but if the habitat is mutated, then the optimal or better solution may appear.In contrast, if the population probability of a habitat is high, the probability of a new solution mutated is very low.Based on this, the equation for mutation operation can be obtained as follows: in which m max is the maximum mutation rate of the userdefined population and P max is the migration probability when the population size of the habitat reaches the maximum.

Design of WSN Coverage Path Optimization Algorithm
Since the BBO algorithm is an iterative algorithm, a fitness function is needed to judge the degree of superiority of each iteration of the algorithm.For the fitness function of the sensor nodes in the process of WSN path optimization, not only the energy of the node must be considered but also the distance to the neighbor node and the energy consumption during communication.In order to further simplify the algorithm, the distance between nodes can be used to measure the energy consumption during communication.e specific fitness function is as follows: where t is the sensor node number, n is the total number of sensor nodes, c t is the current energy value of sensor node t, c i is the remaining energy of sensor node i for sleep, d i is the Euclidean distance between sensor node i and current node t, ω 1 is the energy weighting coefficient of sensor node t, ω 2 is the energy weighting coefficient of neighboring nodes, and ω 1 +ω 2 � 1 and 0 < ω 1 , ω 2 < 1.
4.1.Setting Parameters.N sensor nodes are deployed in the monitoring area, and the coordinates, radius, initial energy, and other information of the sensor nodes are known.Assuming that the total number of habitats is n, the fitness variable of each habitat has D dimensions, the maximum can habitat accommodate the population is S max , the maximum mutation rate of the population is m max , the migration probability when the habitat population reaches the maximum is P max , the maximum number of iterations of the algorithm is MaxNum, and the number of elite individuals retained in each generation is z.

Algorithm Flow.
e implementation process of the WSN coverage path optimization algorithm is as follows: Step 1: e fitness of energy consumption of each habitat is calculated according to equations ( 3) and ( 14), and it is taken as the adaptation index HSI of each habitat.

Complexity
Step 2: Habitats are ranked according to HSI value from small to large.
Step 3: e habitat data of group Z with the largest value of HSI were saved as the elite group.
Step 4: e corresponding immigration rate λ i and emigration rate μ i are calculated according to the population number of each habitat.According to the pseudocode of the above migration operation, each habitat is migrated.
Step 5: e mutation rate of each habitat is calculated according to equation (13), and then mutation operation is carried out on each habitat according to their respective mutation rate.
Step 6: e HSI of each habitat is calculated, and the habitats are sorted from small to large.
Step 7: e retained data of z elite groups are assigned to the habitat group z with the lowest HSI at this time.
Step 8: e HSI of each habitat was calculated again and sorted from small to large, and the optimal value was recorded.
Step 9: If the termination condition of the algorithm is met (the maximum number of iterations is MaxNum, or the global optimal value is found), the algorithm is terminated; otherwise, it returns to Step (4) to continue execution.

Simulation Results and Analysis
e simulation experiment is carried out in Matlab2019 (b), and the BBO algorithm is compared with the PSO algorithm and genetic algorithm.
e experiment is divided into group A and group B. In the experiment of group A, 20 sensor nodes are regularly deployed in a two-dimensional plane and the data of sensor nodes are known, which is shown in Figure 2. Group B randomly deployed 20 sensor nodes in a two-dimensional plane, and the data of sensor nodes were also known, which is shown in Figure 3.In the two sets of experiments, the number of iterations of each algorithm is 200.In order to make the final results without loss of generality, 100 experiments were carried out for the three algorithms in the two sets of experiments, and the final results were averaged.
e WSN coverage paths were optimized by the BBO algorithm for the sets of experiments which are shown in Figures 4 and 5, respectively.e comparison results of the three algorithms in the two sets of experiments are shown in Tables 1 and 2 and Figures 6 and 7.
It can be seen from Figure 6 and Table 1 that in the experiment of group A, the distribution of communication links between the sensor nodes after the rule deployment is relatively even, the optimization result of the PSO algorithm has good adaptability and success rate, while the BBO algorithm further increased the fitness to 0.9902, and no nodes that exceeded the threshold appeared.Figure 5 shows the WSN coverage path optimization results of the BBO algorithm in the A group experiment.
rough the above experimental results of Figure 7 and Table 2, it can be observed that the relative positions deployed by the sensor nodes are distributed irregularly, and the construction of the coverage path is more complicated.
e success rate of the genetic algorithm and the PSO algorithm in the coverage path optimization and the fitness is not ideal enough, and there are many nodes that exceed the threshold.Although there are 2 nodes exceeding the threshold in the optimization results of the BBO algorithm, its success rate and adaptability are better than those of the genetic algorithm and the PSO algorithm.Figure 5 also shows the WSN coverage path of the BBO algorithm in the B group experiment.As a result of the optimization, the    6 Complexity consumption after the network is stabilized.Future research will be inclined to further improve the convergence speed of the algorithm on the premise of ensuring the optimization of the coverage path.

Figure 2 :
Figure 2: Initial deployment of experimental WSN in group A.

Figure 3 :
Figure 3: Initial deployment of experimental WSN in group B.

Figure 4 :Figure 6 :Figure 7 :
Figure 4: Optimization of WSN coverage path based on the BBO algorithm in group A.

Figure 5 :
Figure 5: WSN coverage path optimization based on the BBO algorithm in group B.