CN116828484B - Sensor network coverage optimization method based on improved species life-to-life optimization algorithm - Google Patents
Sensor network coverage optimization method based on improved species life-to-life optimization algorithm Download PDFInfo
- Publication number
- CN116828484B CN116828484B CN202311059852.3A CN202311059852A CN116828484B CN 116828484 B CN116828484 B CN 116828484B CN 202311059852 A CN202311059852 A CN 202311059852A CN 116828484 B CN116828484 B CN 116828484B
- Authority
- CN
- China
- Prior art keywords
- species
- particle
- dimension
- max
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 20
- 241000894007 species Species 0.000 claims abstract description 148
- 239000002245 particle Substances 0.000 claims description 100
- 230000006870 function Effects 0.000 claims description 39
- 238000001212 derivatisation Methods 0.000 claims description 24
- 238000004880 explosion Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 230000008602 contraction Effects 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 claims description 6
- 230000008033 biological extinction Effects 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 4
- 241000854291 Dianthus carthusianorum Species 0.000 claims description 3
- 206010054949 Metaplasia Diseases 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000015689 metaplastic ossification Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
- H04W16/18—Network planning tools
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
- H04W16/22—Traffic simulation tools or models
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W84/00—Network topologies
- H04W84/18—Self-organising networks, e.g. ad-hoc networks or sensor networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
A method for optimizing sensor network coverage using an improved specie life-to-death algorithm is provided, the method comprising the steps of: establishing a sensor probability coverage model; clustering; optimizing the position deployment of the rest dynamic sensor nodes by adopting an SED-PSO algorithm; and outputting an optimal position deployment result of the sensor node. The method is used for optimizing the coverage rate of the sensor network, can effectively realize the coverage optimization of the sensor network based on the probability coverage model, and has the advantages of high convergence rate, high stability and high solution quality of the algorithm.
Description
Technical Field
The invention relates to a sensor network coverage optimization technology, in particular to a method for optimizing the coverage rate of a sensor network by utilizing an improved species life-to-death algorithm.
Background
Coverage rate is a measure of detection of a monitoring area of a sensor network and reflects the quality of perceived service that the sensor network can provide. At present, the optimization method for the sensor network coverage is mainly based on a binary coverage model (a disc model or a 0-1 model), and two types of sensor network coverage methods are formed: firstly, a sensor network coverage method based on geometry is provided, and the method solves the optimal position deployment of nodes through the geometrical relationship among the nodes; secondly, a sensor network coverage method based on a group intelligent optimization algorithm is adopted, so that complicated geometric deduction is avoided by adopting the group intelligent optimization algorithm, but the problems of easy local optimum, premature convergence, large iteration times and the like exist.
Because the coverage model of the sensor is not a theoretical disc coverage model, but a probability coverage model, but the coverage optimization of the sensor network based on the probability coverage model is often complex and difficult to apply a geometric method, the position deployment of the sensor nodes is simpler, more convenient and easier by adopting a group intelligent optimization algorithm for the sensor network coverage optimization aiming at the probability coverage model conforming to the actual detection of the sensor.
The species life-to-death algorithm (system engineering and electronic technology, 2018, 40 (4): 941-947) is a new evolutionary computing (evolutionary computing, EC) algorithm which is proposed by the inventor according to the large explosion and large extinction phenomena of the species which have been found in the biological evolution history and by referring to the thought of the species catastrophe evolution theory, realizes optimization by executing the large explosion and large extinction on the species, achieves the purpose of balancing the global optimizing capability and the local optimizing capability of the algorithm by introducing strategies such as main branch transfer, species derivative capability contraction and the like, and has higher quality optimal solution and higher stability.
Disclosure of Invention
In order to solve the problem that the geometrical optimization method is difficult to adopt in the sensor network coverage optimization based on the probability coverage model, and simultaneously in order to further improve the optimization performance of the species life-to-death (species explode AND DERACINATE, SED) algorithm, the invention provides a sensor network coverage optimization method based on an improved species life-to-death optimization algorithm, which specifically comprises the following steps:
the first step: establishing a sensor probability coverage model shown in a formula (1);
(1)
Wherein, Representing the perceived probability of the jth sensor for an object at a distance R, the sensor perceived radius R s, the uncertainty perceived range delta x R s, delta being a coefficient less than 1,/>Representing the distance between the target T and the sensor S j, x t、yt representing the coordinates at the target point,/>, respectively、/>Respectively representing coordinates of the sensor S j; lambda 1、λ2、β1、β2 is a parameter related to the sensor characteristics, alpha 1、α2 is a first and second input parameter, respectively, the values of which are related to r s and d (S j, T), the relation being
(2)
And a second step of: clustering; clustering the static sensors by adopting a K-means algorithm, and dividing all sensor nodes into 4 clusters; a dynamic sensor is respectively arranged at the central position of each cluster and the central position of the optimization area to serve as a cluster head and a network central point;
And a third step of: optimizing the position deployment of the rest dynamic sensor nodes by adopting an SED-PSO algorithm; the SED-PSO algorithm comprises the following steps:
(1) Setting parameters of a species life-to-death algorithm, comprising: the number S of surviving species, the derivatization capability A of the species, the generation number G of the species, the explosion multiple M of the species, the contraction coefficient alpha and the iteration number FE max of the algorithm;
(2) Setting PSO algorithm parameters, including: population size P, first learning factor c 1 and second learning factor c 2, inertial weight factor omega, maximum positive movement speed v max and maximum negative movement speed v -max of particles, and iterative times iter of algorithm;
(3) Establishing an adaptability function; the fitness function of the sensor network coverage based on the improved species life-to-life optimization algorithm is the reciprocal of coverage;
(4) Initializing a species life-to-death algorithm in a randomizing mode according to the range of the optimized region, and calculating the fitness value of each species; each dimension of each species initializes the formula:
(4)
(5)
Wherein X s,d represents the value of the d dimension of the S-th species, a s represents the derivatization capacity of the S-th species, s=1, 2, …, S represents the S-th surviving species; d=1, 2, …, D representing the D-th dimension of the optimization problem, D representing the dimension of the optimization problem, rand representing the random number between (0, 1); [ L min, Lmax ] represents a search interval, and L min, Lmax is the lower limit and the upper limit of the search interval respectively; a epsilon [ A min, Amax],Amax ] represents the maximum species derivatization capacity, and A min is more than or equal to 0 represents the minimum species derivatization capacity;
(5) Performing an outbreak on each surviving species; the generation of new species at each major burst is shown in the following formula;
(6)
(7)
Wherein, Representing the new species derived from the s-th species in the kth optimization at the j-th outbreak,/>The s-th surviving species expressed in the kth optimization, rand is AND/>Random vector between (0, 1) of the same dimension,/>Representing the derivatization ability of the s-th species in the kth optimization,/>Representing the derivatization capacity of a new species derivatized by a jth outbreak of an s-th species in a kth optimization, j=1, 2, …, M representing a jth outbreak;
(6) Judging and processing out-of-range of new species generated by each burst; each dimension of the species generated by the outbreak should have a value within the interval [ L min, Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range species with a random number in a [ L min, Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
(7) Performing deep iteration on new species generated by each burst by adopting a PSO algorithm;
① Based on new species generated by explosion, generating N particles, wherein the generation formula of each dimension of each particle is as follows:
(8)
(9)
Wherein, Value representing the d-th dimension of the i-th particle,/>A speed value representing the d-th dimension of the i-th particle, i=1, 2, …, N representing the i-th particle; /(I)Representation/>D-th dimensional value of (2);
② Carrying out boundary crossing judgment and treatment on each particle; each dimension of each particle has a value within the interval [ L min, Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range particles with a random number in a [ L min, Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
③ According to the established fitness function, calculating fitness values of N particles, selecting a global optimal solution g, and enabling the particle self optimal solution to be the current value of the particle;
④ Updating each dimension of each particle velocity according to the following formula;
(10)
wherein ω is inertial weight; d-th dimension of speed at t-th iteration of the ith particle,/> Is the d dimension of the self history optimal position of the ith particle at the t time of iteration,/>Is the d dimension of the position of the ith particle at the t iteration,The d dimension of the optimal position of the group at the t-th iteration; and c 1 and c 2 are acceleration factors, which are constants greater than 0, and are used for controlling the influence of the particle history optimal position and the overall optimal particle position of the population on the current particle speed. r 1 and r 2 are random numbers between (0, 1);
Judging whether the speed is in the interval of [ - (L max-Lmin),( Lmax - Lmin) ] or not, if not, replacing the speed by using a random function to generate random numbers in the interval of [ - (L max - Lmin), (Lmax - Lmin) ] and if so, not processing;
updating the position of the particles according to the following formula;
(11)
Judging whether each dimension of each updated particle is in the [ L min, Lmax ] interval, if not, replacing the updated dimension by a random number in the [ L min, Lmax ] interval generated by a random function; if the interval is within the interval, no processing is needed;
⑤ After the particle speed and the position are updated each time, calculating an fitness function value of the updated particle, if the fitness function value is smaller than the fitness function value of the historical optimal solution of the particle, replacing the historical optimal solution of the particle with the position of the updated particle, otherwise, keeping the historical optimal solution of the particle unchanged;
After all particles complete the position updating of the iteration, selecting the particle with the smallest fitness function value of all updated particles, comparing the fitness value of the particle with the fitness value of the global optimal solution, and if the fitness value of the particle is smaller, replacing the global optimal solution with the position of the particle; otherwise, the global optimal solution remains unchanged;
⑥ Let t=t+1 and repeat ③-⑤ steps until iter iterations are completed;
⑦ Replacing the value of the new species generated by the current burst operation with the iterated global optimal solution g iter, namely: ;
(8) Calculating and comparing fitness function values of the species after and before the explosion, and marking the new species as the main species if the fitness value of the new species is smaller, otherwise, the old species as the main species; only the main species can perform the large burst operation; if a major transfer occurs, the generation number G of the species is increased by 1, and the derivatization capacity of the major species is updated to the derivatization capacity of the new species as shown in the following formula;
(12)
(13)
(9) Repeating the steps (5) to (8) until M bursts are completed;
(10) Sequencing the fitness function values of all new species generated by M bursts and existing species before the bursts, selecting S species with the optimal fitness value as metaplasia species for the next large burst and large extinction, and contracting the derivative capacity of surviving species according to the following formula, wherein all other species are eliminated;
(14)
Wherein, Is the coefficient of contraction;
(11) Repeating the steps (5) to (10) until a termination condition is reached (the termination condition is generally that the maximum iteration number is reached or the required optimizing precision is reached);
fourth step: and outputting an optimal position deployment result of the sensor node.
In one embodiment of the invention, in the third step (3), the coverage is calculated as follows: dividing a monitoring area into [100×100] grids according to granularity g=40m, and calculating the proportion of the grids with the sensor network sensing probability of more than p th =0.6 in the center of the grids to all the grids; the sensing probability of the sensor network is as follows
(3)。
In one embodiment of the invention, the parameters are shown in Table 1;
TABLE 1 scene and algorithm parameters
Scene and algorithm | Parameters (parameters) |
Monitoring area: 4km×4km | Lmin=0,Lmax=4km; |
Number of nodes: 38 sensor nodes, wherein 20 sensor nodes are static sensors, and the rest are dynamic sensors; | D=(38-20-4-1)×2=13×2=26; |
Probability coverage model parameters of (1) | rs=1.2km,δ=0.6,λ1=1,λ2=0.6,β1=3,β2=2; |
SED-PSO algorithm | M=20,S=30,α=0.05,FEmax=2000,N=2,ω=0.724,c1=c2=2,iter=5; |
According to the invention, an SED algorithm and a Particle Swarm Optimization (PSO) algorithm are combined, an improved SED algorithm fused with the PSO algorithm is provided, the SED algorithm is recorded as an SED-PSO algorithm, and a sensor network coverage optimization method based on an improved species life-to-life optimization algorithm is provided based on a sensor probability coverage model.
Drawings
FIG. 1 shows a schematic flow diagram of a sensor network coverage optimization method based on an improved species life-to-life optimization algorithm;
FIG. 2 shows a sensor node location profile after optimization using the SED-PSO algorithm;
FIG. 3 shows sensor network coverage and coverage after optimization using the SED-PSO algorithm;
FIG. 4 shows the average coverage over the number of iterations for 30 Monte Carlo (Monte-Carlo) trials;
FIG. 5 shows the results of optimization of 30 Monte-Carlo experiments.
Detailed Description
The present invention is described in detail below with reference to the accompanying drawings.
The detailed steps of the sensor network coverage optimization method based on the improved species life-to-death optimization algorithm are as follows.
The first step: a probabilistic coverage model of the sensor is established as shown in formula (1) (the coverage model can also be established according to the actual sensing condition of the sensor).
(1)
Wherein,Representing the perceived probability of the jth sensor for an object at a distance R, the sensor perceived radius R s, the uncertainty perceived range delta x R s, delta being a coefficient less than 1, beta 1、β2 being a parameter related to the sensor characteristics, alpha 1、α2 being the first and second input parameters, respectively, the values/>Representing the distance between the target T and the sensor S j, x t、yt representing the coordinates at the target point,/>, respectively、/>The coordinates of the sensor S j are indicated respectively. Lambda 1、λ2, related to r s and d (S j, T), the relation of which is as follows
(2)
And a second step of: clustering. And clustering the static sensors by adopting a K-means (K-means) algorithm, and dividing all sensor nodes into 4 clusters. And arranging a dynamic sensor at the central position of each cluster and the central position of the optimized area as a cluster head and a network central point.
And a third step of: and optimizing the position deployment of the rest dynamic sensor nodes by adopting an SED-PSO algorithm. The SED-PSO algorithm comprises the following steps:
(1) Setting parameters of a species life-stop algorithm, comprising: number of surviving species S, derivatization ability of species a, number of generation of species G, number of outbreaks of species M, coefficient of contraction α, number of iterations of algorithm FE max.
(2) Setting parameters of a PSO algorithm, including: the population size P, the first and second learning factors c 1 and c 2, the inertial weight factor ω, the maximum positive velocity v max and the maximum negative velocity v -max of the particles, the number of iterations iter of the algorithm.
(3) And establishing an adaptability function. The fitness function of the sensor network coverage based on the improved species life-stop optimization algorithm is the reciprocal of the coverage rate. The coverage is calculated as follows: dividing a monitoring area into [100×100] grids according to granularity g=40m, and calculating the proportion of the grids with the sensor network sensing probability of more than p th =0.6 at the center of the grids to all the grids. The sensing probability of the sensor network is as follows
(3)
(4) Initializing a species life-to-death algorithm in a randomization mode according to the range of the optimized region, and calculating the fitness value of each species; each dimension of each species initializes the formula:
(4)
(5)
wherein X s,d represents the value of the d dimension of the S-th species, a s represents the derivatization capacity of the S-th species, s=1, 2, …, S represents the S-th surviving species; d=1, 2, …, D representing the D-th dimension of the optimization problem, D representing the dimension of the optimization problem, rand representing the random number between (0, 1); [ L min, Lmax ] represents a search interval, and L min, Lmax is the lower limit and the upper limit of the search interval, respectively. A.epsilon.A min, Amax],Amax represents the maximum species derivatization capacity and A min.gtoreq.0 represents the minimum species derivatization capacity.
(5) Performing an outbreak on each surviving species; the generation of new species at each major burst is shown in the following formula.
(6)
(7)
Wherein,Representing the new species derived from the s-th species in the kth optimization at the j-th outbreak,/>The s-th surviving species expressed in the kth optimization, rand is AND/>Random vector between (0, 1) of the same dimension,/>Representing the derivatization ability of the s-th species in the kth optimization,/>Represents the derivatization ability of the new species derivatized by the jth outbreak of the jth species in the kth optimization, j=1, 2, …, M, representing the jth outbreak.
(6) Judging and processing out-of-range of new species generated by each burst; each dimension of the species generated by the outbreak should have a value within the interval [ L min, Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range species with a random number in a [ L min, Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
(7) Performing deep iteration on new species generated by each burst by adopting a PSO algorithm;
① Based on new species generated by explosion, generating N particles, wherein the generation formula of each dimension of each particle is as follows:
(8)
(9)
Wherein, Value representing the d-th dimension of the i-th particle,/>A speed value representing the d-th dimension of the i-th particle, i=1, 2, …, N representing the i-th particle; /(I)Representation/>D-th dimensional value of (c).
② Carrying out boundary crossing judgment and treatment on each particle; each dimension of each particle has a value within the interval [ L min, Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range particles with a random number in a [ L min, Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
③ According to the established fitness function, calculating fitness values of N particles, selecting a global optimal solution g, and enabling the particle self optimal solution to be the current value of the particle.
④ Updating each dimension of each particle velocity according to the following formula;
(10)
wherein ω is inertial weight; d-th dimension of speed at t-th iteration of the ith particle,/> Is the d dimension of the self history optimal position of the ith particle at the t time of iteration,/>Is the d dimension of the position of the ith particle at the t iteration,The d dimension of the optimal position of the group at the t-th iteration; and c 1 and c 2 are acceleration factors, which are constants greater than 0, and are used for controlling the influence of the particle history optimal position and the overall optimal particle position of the population on the current particle speed. r 1 and r 2 are random numbers between (0, 1).
Judging whether the speed is in the interval of [ - (L max-Lmin),( Lmax - Lmin) ] or not, if not, replacing the speed by using a random function to generate random numbers in the interval of [ - (L max - Lmin), (Lmax - Lmin) ] and if so, not processing;
updating the position of the particles according to the following formula;
(11)
Judging whether each dimension of each updated particle is in the [ L min, Lmax ] interval, if not, replacing the updated dimension by a random number in the [ L min, Lmax ] interval generated by a random function; if the interval is within the interval, no processing is needed;
⑤ After the particle speed and the position are updated each time, calculating an fitness function value of the updated particle, if the fitness function value is smaller than the fitness function value of the historical optimal solution of the particle, replacing the historical optimal solution of the particle with the position of the updated particle, otherwise, keeping the historical optimal solution of the particle unchanged;
After all particles complete the position updating of the iteration, selecting the particle with the smallest fitness function value of all updated particles, comparing the fitness value of the particle with the fitness value of the global optimal solution, and if the fitness value of the particle is smaller, replacing the global optimal solution with the position of the particle; otherwise, the global optimal solution remains unchanged;
⑥ Let t=t+1 and repeat ③-⑤ steps until iter iterations are completed;
⑦ Replacing the value of the new species generated by the current burst operation with the iterated global optimal solution g iter, namely: ;
(8) Calculating and comparing fitness function values of the species after and before the explosion, and marking the new species as the main species if the fitness value of the new species is smaller, otherwise, the old species as the main species; only the main species can perform the large burst operation; if a major transfer occurs, the generation number G of the species is increased by 1, and the derivatization capacity of the major species is updated to the derivatization capacity of the new species as shown in the following formula;
(12)
(13)
(9) Repeating the steps (5) to (8) until M bursts are completed;
(10) Sequencing the fitness function values of all new species generated by M bursts and existing species before the bursts, selecting S species with the optimal fitness value as metaplasia species for the next large burst and large extinction, and contracting the derivative capacity of surviving species according to the following formula, wherein all other species are eliminated;
(14)
Wherein, Is the coefficient of contraction.
(11) Repeating steps (5) to 10 until a termination condition is reached (the termination condition is typically the maximum number of iterations or the required optimization accuracy is reached).
Fourth step: and outputting an optimal position deployment result of the sensor node.
In the present invention, the parameters in the specific embodiments are shown in table 1.
TABLE 1 scene and algorithm parameters
Scene and algorithm | Parameters (parameters) |
Monitoring area: 4km×4km | Lmin=0,Lmax=4km。 |
Number of nodes: 38 sensor nodes, wherein 20 sensor nodes are static sensors, and the rest are dynamic sensors. | D=(38-20-4-1)×2=13×2=26。 |
Probability coverage model parameters of (1) | rs=1.2km,δ=0.6,λ1=1,λ2=0.6,β1=3,β2=2。 |
SED-PSO algorithm | M=20,S=30,α=0.05,FEmax=2000,N=2,ω=0.724,c1=c2=2,iter=5。 |
Claims (3)
1. A sensor network coverage optimization method based on an improved species life-to-death optimization algorithm is characterized by comprising the following steps:
the first step: establishing a sensor probability coverage model shown in a formula (1);
Wherein, Representing the perceived probability of the jth sensor for an object at a distance R, the sensor perceived radius R s, the uncertainty perceived range delta x R s, delta being a coefficient less than 1,/>Representing the distance between the target T and the sensor S j, x t、yt representing the coordinates at the target point, respectively, and x s,j、ys,j representing the coordinates of the sensor S j, respectively; lambda 1、λ2、β1、β2 is a parameter related to the sensor characteristics, alpha 1、α2 is a first and second input parameter, respectively, the values of which are related to r s and d (S j, T), the relation being
And a second step of: clustering; clustering the static sensors by adopting a K-means algorithm, and dividing all sensor nodes into 4 clusters; a dynamic sensor is respectively arranged at the central position of each cluster and the central position of the optimization area to serve as a cluster head and a network central point;
And a third step of: optimizing the position deployment of the rest dynamic sensor nodes by adopting an SED-PSO algorithm; the SED-PSO algorithm comprises the following steps:
(1) Setting parameters of a species life-to-death algorithm, comprising: the number S of surviving species, the derivatization capability A of the species, the generation number G of the species, the explosion multiple M of the species, the contraction coefficient alpha and the iteration number FE max of the algorithm;
(2) Setting PSO algorithm parameters, including: population size P, first learning factor c 1 and second learning factor c 2, inertial weight factor omega, maximum positive movement speed v max and maximum negative movement speed v -max of particles, and iterative times iter of algorithm;
(3) Establishing an adaptability function; the fitness function of the sensor network coverage based on the improved species life-to-life optimization algorithm is the reciprocal of coverage;
(4) Initializing a species life-to-death algorithm in a randomizing mode according to the range of the optimized region, and calculating the fitness value of each species; each dimension of each species initializes the formula:
Xs,d=Lmax+(Lmin-Lmax)×rand (4)
As=Amax-Amin (5)
Wherein X s,d represents the value of the d dimension of the S-th species, a s represents the derivatization capacity of the S-th species, s=1, 2, …, S represents the S-th surviving species; d=1, 2, …, D representing the D-th dimension of the optimization problem, D representing the dimension of the optimization problem, rand representing the random number between (0, 1); [ L min,Lmax ] represents a search interval, and L min,Lmax is the lower limit and the upper limit of the search interval respectively; a epsilon [ A min,Amax],Amax ] represents the maximum species derivatization capacity, and A min is more than or equal to 0 represents the minimum species derivatization capacity;
(5) Performing an outbreak on each surviving species; the generation of new species at each major burst is shown in the following formula;
Wherein, Representing the new species derived from the s-th species in the kth optimization at the j-th outbreak,/>The s-th surviving species expressed in the kth optimization, rand is AND/>Random vector between (0, 1) of the same dimension,/>Representing the derivatization ability of the s-th species in the kth optimization,/>Representing the derivatization capacity of a new species derivatized by a jth outbreak of an s-th species in a kth optimization, j=1, 2, …, M, representing a jth outbreak;
(6) Judging and processing out-of-range of new species generated by each burst; each dimension of the species generated by the outbreak should have a value within the interval [ L min,Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range species with a random number in a [ L min,Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
(7) The PSO algorithm is adopted to carry out deep iteration on new species generated by each burst:
① Based on new species generated by explosion, generating N particles, wherein the generation formula of each dimension of each particle is as follows:
vi,d=(Lmax-Lmin)×(rand-0.5) (9)
where x i,d represents the value of the d-th dimension of the i-th particle, v i,d represents the speed value of the d-th dimension of the i-th particle, i=1, 2, …, N represents the i-th particle; representation/> D-th dimensional value of (2);
② Carrying out boundary crossing judgment and treatment on each particle; each dimension of each particle has a value within the interval [ L min,Lmax ], and if not, the dimension is beyond the boundary; replacing the dimension value of the out-of-range particles with a random number in a [ L min,Lmax ] interval generated by a random function; if the boundary is not crossed, no treatment is needed;
③ According to the established fitness function, calculating fitness values of N particles, selecting a global optimal solution g, and enabling the particle self optimal solution to be the current value of the particle;
④ Updating each dimension of each particle velocity according to the following formula;
wherein ω is inertial weight; d-th dimension of speed at t-th iteration of the ith particle,/> Is the d dimension of the self history optimal position of the ith particle at the t time of iteration,/>Is the d dimension of the position of the ith particle at the t-th iteration,/>The d dimension of the optimal position of the group at the t-th iteration; c 1 and c 2 are acceleration factors, which are constants greater than 0, and are used for controlling the influence of the particle history optimal position and the population overall optimal particle position on the current particle speed, and r 1 and r 2 are random numbers between (0, 1);
Judging whether the speed is in the interval of [ - (L max-Lmin),(Lmax-Lmin) ] or not, if not, replacing the speed by using a random function to generate random numbers in the interval of [ - (L max-Lmin),(Lmax-Lmin) ] and if so, not processing;
updating the position of the particles according to the following formula;
Judging whether each dimension of each updated particle is in the [ L min,Lmax ] interval, if not, replacing the updated dimension by a random number in the [ L min,Lmax ] interval generated by a random function; if the interval is within the interval, no processing is needed;
⑤ After the particle speed and the position are updated each time, calculating an fitness function value of the updated particle, if the fitness function value is smaller than the fitness function value of the historical optimal solution of the particle, replacing the historical optimal solution of the particle with the position of the updated particle, otherwise, keeping the historical optimal solution of the particle unchanged;
After all particles complete the position updating of the iteration, selecting the particle with the smallest fitness function value of all updated particles, comparing the fitness value of the particle with the fitness value of the global optimal solution, and if the fitness value of the particle is smaller, replacing the global optimal solution with the position of the particle; otherwise, the global optimal solution remains unchanged;
⑥ Let t=t+1 and repeat ③-⑤ steps until iter iterations are completed;
⑦ Replacing the value of the new species generated by the current burst operation with the iterated global optimal solution g iter, namely:
(8) Calculating and comparing fitness function values of the species after and before the explosion, and marking the new species as the main species if the fitness value of the new species is smaller, otherwise, the old species as the main species; only the main species can perform the large burst operation; if a major transfer occurs, the generation number G of the species is increased by 1, and the derivatization capacity of the major species is updated to the derivatization capacity of the new species as shown in the following formula;
(9) Repeating the steps (5) to (8) until M bursts are completed;
(10) Sequencing the fitness function values of all new species generated by M bursts and existing species before the bursts, selecting S species with the optimal fitness value as metaplasia species for the next large burst and large extinction, and contracting the derivative capacity of surviving species according to the following formula, wherein all other species are eliminated;
wherein α is the coefficient of contraction;
(11) Repeating the steps 5 to 10 until the termination condition is reached;
fourth step: and outputting an optimal position deployment result of the sensor node.
2. The sensor network coverage optimization method based on the improved species life-to-death optimization algorithm according to claim 1, wherein in the third step (3), the coverage rate is calculated as follows, the monitoring area is divided into [100×100] grids according to the granularity g=40m, and the proportion of the grids with the sensor network sensing probability larger than p th =0.6 at the center of the grids to all grids is calculated; the sensing probability of the sensor network is as follows
3. The method of claim 1, wherein in step (11), the termination condition is that the maximum number of iterations is reached or that the required optimizing accuracy is reached.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311059852.3A CN116828484B (en) | 2023-08-22 | 2023-08-22 | Sensor network coverage optimization method based on improved species life-to-life optimization algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311059852.3A CN116828484B (en) | 2023-08-22 | 2023-08-22 | Sensor network coverage optimization method based on improved species life-to-life optimization algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116828484A CN116828484A (en) | 2023-09-29 |
CN116828484B true CN116828484B (en) | 2024-06-07 |
Family
ID=88126027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311059852.3A Active CN116828484B (en) | 2023-08-22 | 2023-08-22 | Sensor network coverage optimization method based on improved species life-to-life optimization algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116828484B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103945395A (en) * | 2014-02-27 | 2014-07-23 | 北京理工大学 | Particle-swarm-based rapid optimization deployment method for wireless network sensors |
WO2018072351A1 (en) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | Method for optimizing support vector machine on basis of particle swarm optimization algorithm |
CN112163665A (en) * | 2020-10-09 | 2021-01-01 | 中国人民解放军空军工程大学 | PSO algorithm-based species birth and death optimization method |
CN114827929A (en) * | 2022-03-16 | 2022-07-29 | 西安电子科技大学 | Multi-sensor network collaborative deployment method and device based on improved particle swarm optimization |
CN115278808A (en) * | 2022-06-22 | 2022-11-01 | 浙江工商大学 | SDWSN-oriented distributed efficient entropy energy-saving clustering routing method |
CN115396905A (en) * | 2021-05-25 | 2022-11-25 | 鲁东大学 | Wireless sensor network coverage optimization method based on improved genetic algorithm |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105704732B (en) * | 2014-11-27 | 2019-01-25 | 中国科学院沈阳自动化研究所 | Relay node robust covering method towards double-layer structure wireless sensor network |
-
2023
- 2023-08-22 CN CN202311059852.3A patent/CN116828484B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103945395A (en) * | 2014-02-27 | 2014-07-23 | 北京理工大学 | Particle-swarm-based rapid optimization deployment method for wireless network sensors |
WO2018072351A1 (en) * | 2016-10-20 | 2018-04-26 | 北京工业大学 | Method for optimizing support vector machine on basis of particle swarm optimization algorithm |
CN112163665A (en) * | 2020-10-09 | 2021-01-01 | 中国人民解放军空军工程大学 | PSO algorithm-based species birth and death optimization method |
CN115396905A (en) * | 2021-05-25 | 2022-11-25 | 鲁东大学 | Wireless sensor network coverage optimization method based on improved genetic algorithm |
CN114827929A (en) * | 2022-03-16 | 2022-07-29 | 西安电子科技大学 | Multi-sensor network collaborative deployment method and device based on improved particle swarm optimization |
CN115278808A (en) * | 2022-06-22 | 2022-11-01 | 浙江工商大学 | SDWSN-oriented distributed efficient entropy energy-saving clustering routing method |
Non-Patent Citations (10)
Title |
---|
A PSO Based Routing with Novel Fitness Function for Improving Lifetime of WSNs;Damodar Reddy Edla;《Wireless Personal Communications》;20180929;全文 * |
A wireless sensor network coverage optimization algorithm based on particle swarm optimization and Voronoi diagram;Nor Azlina Bt. Ab Aziz;《2009 International Conference on Networking, Sensing and Control》;20090329;全文 * |
Coverage optimization of sensor network based on improved particle swarm optimization;Yang Yongjian;《Systems Engineering and Electronics》;20170713;全文 * |
Data Aggregation in WSN Based on SOFM Neural Network;Yang Yongjian;《Chinese Journal of Sensors and Actuators》;20140731;全文 * |
En Wang ; Yongjian Yang ; Jie Wu.Energy Efficient Phone-to-Phone Communication Based on WiFi Hotspot in PSN.《2015 24th International Conference on Computer Communication and Networks (ICCCN)》.2015,全文. * |
Particle Swarm Optimization in Wireless-Sensor Networks: A Brief Survey;Raghavendra V. Kulkarni;《IEEE Transactions on Systems, Man, and Cybernetics, Part C》;20200723;全文 * |
Species explode and deracinate algorithm;Y. Yang;《Systems Engineering and Electronics》;20180401;全文 * |
基于IEEE 802.15.4的星形拓扑无线传感器网络性能分析;CAI Jie, GU Ming;《清华大学学报(自然科学版) 》;20151231;全文 * |
物种生灭算法;杨永建;系统工程与电子技术;20181231;全文 * |
物种生灭算法的改进策略;杨永建;《计算机工程与应用 》;20191231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116828484A (en) | 2023-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107808212B (en) | Solar energy collection power prediction method based on grey neural network | |
CN112698634B (en) | Event trigger-based traffic intelligent system fixed time dichotomy consistency method | |
CN106162663B (en) | A kind of sensing node covering method based on improvement ant colony algorithm | |
CN112887994B (en) | Wireless sensor network optimization method based on improved binary particle swarm and application | |
CN110030843B (en) | Heat accumulating type aluminum melting furnace parameter optimization setting method based on improved whale optimization algorithm | |
Zhang et al. | Energy‐efficient adaptive dynamic sensor scheduling for target monitoring in wireless sensor networks | |
CN113905389A (en) | Wireless sensor network coverage method based on particle swarm optimization imperial butterfly algorithm | |
CN115866621B (en) | Wireless sensor network coverage method based on whale algorithm | |
CN111553469A (en) | Wireless sensor network data fusion method, device and storage medium | |
CN114580786A (en) | Arrangement optimization method for monitoring points of urban water supply pipe network | |
Azharuddin et al. | A GA-based approach for fault tolerant relay node placement in wireless sensor networks | |
CN112580855A (en) | Cable group steady-state temperature rise prediction method based on self-adaptive variation PSO-BP neural network | |
CN116828484B (en) | Sensor network coverage optimization method based on improved species life-to-life optimization algorithm | |
CN108805253B (en) | PM2.5 concentration prediction method | |
CN113115342B (en) | WSNs deployment method and system of virtual force-guided sparrow search algorithm | |
CN113095477B (en) | Wind power prediction method based on DE-BP neural network | |
CN107092779B (en) | A kind of sensor actuator method for optimizing position based on minimum across displacement entropy | |
CN107396375B (en) | Quantum brain storm-based optimal target coverage method for heterogeneous sensor network | |
CN114237282B (en) | Unmanned aerial vehicle flight path intelligent planning method for intelligent industrial park monitoring | |
Ou et al. | Modeling wireless sensor network coverage based on random radius | |
CN115759979A (en) | Process intelligent processing method and system based on RPA and process mining | |
CN112434857B (en) | Unmanned aerial vehicle airport network site selection arrangement method based on GA-PSO optimization | |
CN114449610A (en) | Wireless sensing node clustering and cluster head election method based on network division | |
CN111010704B (en) | Underwater wireless sensor network data prediction optimization method based on exponential smoothing | |
CN114048668B (en) | Comprehensive energy-oriented wireless sensor wheel value method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |