CN111447627B - WSN node positioning method based on differential evolution genetic algorithm - Google Patents

WSN node positioning method based on differential evolution genetic algorithm Download PDF

Info

Publication number
CN111447627B
CN111447627B CN202010180046.1A CN202010180046A CN111447627B CN 111447627 B CN111447627 B CN 111447627B CN 202010180046 A CN202010180046 A CN 202010180046A CN 111447627 B CN111447627 B CN 111447627B
Authority
CN
China
Prior art keywords
vector
positioning
steps
formula
population
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
Application number
CN202010180046.1A
Other languages
Chinese (zh)
Other versions
CN111447627A (en
Inventor
吴斌
金洁丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Post And Telecommunication College
Original Assignee
Zhejiang Post And Telecommunication College
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Post And Telecommunication College filed Critical Zhejiang Post And Telecommunication College
Priority to CN202010180046.1A priority Critical patent/CN111447627B/en
Publication of CN111447627A publication Critical patent/CN111447627A/en
Application granted granted Critical
Publication of CN111447627B publication Critical patent/CN111447627B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/22Traffic simulation tools or models
    • H04W16/225Traffic simulation tools or models for indoor or short range network
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention discloses a WSN node positioning method based on a differential evolution genetic algorithm. The method comprises the following steps: step 1, establishing a node positioning optimization model and designing a positioning optimization objective function; step 2, initializing a population; step 3, generating a variation vector by utilizing differential variation operation after the population is subjected to initialization assignment; step 4, after the mutation operation is carried out, a target vector and the mutation vector need to be subjected to binomial intersection to generate a final test vector; step 5, one-to-one tournament selection is executed between the test vector and the corresponding target vector; step 6, updating the individuals by using the success probability; and 7, repeating the steps 3 to 6, and outputting the optimal position value when the maximum iteration times are reached. On the basis of the basic GA algorithm, a differential evolution mechanism is provided, so that the iteration times and the positioning error of the algorithm are considered; the method can give consideration to both network real-time performance and high-precision positioning performance, and has effectiveness and practicability.

Description

WSN node positioning method based on differential evolution genetic algorithm
Technical Field
The invention belongs to the field of node positioning in a wireless sensor network, and particularly relates to a WSN node positioning method using a differential evolution genetic algorithm.
Background
At present, algorithms for solving WSN network positioning are endless, and the existing algorithms can be divided into 2 types according to whether the algorithms measure the distance or not, and the positioning algorithms are based on ranging and do not need ranging. The node location algorithm based on ranging mainly includes received signal strength indication (rssi), time difference of arrival (TDOA), and angle of arrival (AOA). The positioning algorithm without distance measurement mainly includes a centroid algorithm, an (approximate point-in-three-alignment, APIT) algorithm, a (DV-Houpo optimization) DV-Hop algorithm and an Amorphous algorithm. The former is to measure the actual position of the distance estimation node by using some tools, and the distance is not required to be measured when the position of the positioning node is estimated without distance measurement, so that the positioning error based on the distance measurement is far smaller than the positioning algorithm of the node without the distance measurement. Since the former has higher positioning accuracy, the node positioning method based on ranging is mainly considered herein. In the method based on the distance measurement, the RSSI algorithm is widely used because it does not need an additional large amount of hardware and has high positioning accuracy. However, due to the existence of errors in the ranging, a plurality of experts and scholars model the positioning problem of the wireless sensor into an optimization problem for solving the minimum ranging error to reduce the errors, and an intelligent optimization algorithm is adopted to solve the optimization problem. However, single intelligent optimization such as genetic algorithm, particle swarm, differential evolution algorithm and the like often has the defects of poor global optimization capability, low iteration speed, high positioning error and the like, so that the research of an intelligent algorithm with low positioning error and high iteration speed is very important.
Disclosure of Invention
Aiming at the defects that a single intelligent algorithm is easy to fall into a local minimum value, low in optimizing speed, low in precision and the like, a difference optimization mechanism is provided on the basis of a basic genetic algorithm, so that the iteration rate and the positioning performance of the algorithm are considered.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
step 1, establishing a node positioning optimization model and designing a positioning optimization objective function;
step 2, initializing a population;
step 3, generating a variation vector by utilizing differential variation operation after the population is subjected to initialized assignment;
step 4, after performing mutation operation, the target vector x needs to be transformed i,g And the variation vector v i,g Generating final test vector u by carrying out binomial intersection i,g =[u i1,g , u i2,g ,…,u iD,g ];
Step 5, one-to-one tournament selection is executed between the test vectors and the corresponding target vectors;
step 6, updating the individuals by using the success probability;
and 7, repeating the steps 3 to 6, and outputting the optimal position value when the maximum iteration times are reached.
The invention has the beneficial effects that:
1. a WSN node positioning method based on a differential evolution genetic algorithm is provided. The method can give consideration to both network real-time performance and high-precision positioning performance, and has effectiveness and practicability.
2. On the basis of a basic GA algorithm, a differential evolution mechanism is provided, so that the iteration times and the positioning errors of the algorithm are considered.
Drawings
FIG. 1 is a schematic view of the present invention.
Detailed Description
The following describes the implementation steps of the present invention in further detail with reference to the attached drawings.
As shown in fig. 1, a method for positioning a WSN node based on a differential evolution genetic algorithm specifically includes the following steps:
step 1, establishing a node positioning optimization model and designing a positioning optimization objective function;
step 2, initializing a population;
step 3, generating a variation vector by utilizing differential variation operation after the population is subjected to initialization assignment;
step 4, after performing mutation operation, the target vector x needs to be transformed i,g And the variation vector v i,g Generating final test vector u by carrying out binomial intersection i,g =[u i1,g , u i2,g ,…,u iD,g ];
Step 5, one-to-one tournament selection is executed between the test vector and the corresponding target vector;
step 6, updating the individuals by using the success probability;
and 7, repeating the steps 3 to 6, and outputting the optimal position value when the maximum iteration times are reached.
Establishing a node positioning optimization model and designing a positioning optimization objective function in the step 1, specifically, the following steps are performed:
setting anchor P 1 (x 1 ,y 1 ),P 2 (x 2 ,y 2 ),…,P n (x n ,y n ) And the actual true lengths of the unknown nodes P (x, y) are r 1 ,r 2 ,…,r n The difference of the distance measurement errors is epsilon 12 ,…,ε n Then, | r is satisfied i -d i |<ε i Where i =1,2, ..., n, then the unknown node P (x, y) constraint:
Figure BDA0002412182790000031
solving for (x, y), such that:
Figure BDA0002412182790000032
the node location problem is converted into an optimization problem for solving the minimum value of the ranging error.
The initialization population in step 2 is specifically as follows:
default values are given to the chromosome individuals within the limit range, function values Fitness (i) of all chromosomes in the total group are calculated according to a formula (2), and the individual p with the highest Fitness value in the current algebraic individuals is determined g And the position p with the highest personal cutoff algebraic fitness value h
After the population in the step 3 is subjected to initialization assignment, generating a variation vector v by utilizing differential variation operation i,g The method comprises the following steps:
v i,g =x r1,g +F*(x r2,g -x r3,g ) (3)
as shown in FIG. 1, in the formula (3), r 1 ,r 2 ,r 3 ∈[1,2,....NP]The scaling factor F is a constant in the interval of (0, 1) and is used for controlling the size of the differential vector; x is the number of r1,g 、x r2,g And x r3,g Are a vector randomly selected from the population.
After performing mutation operation as described in step 4, the target vector x is required to be transformed i,g And the variation vector v i,g Generating final test vector u by carrying out binomial intersection i,g =[u i1,g , u i2,g ,…,u iD,g ]The method comprises the following steps:
Figure BDA0002412182790000033
in the formula, j rand Is an integer randomly selected from the set {1,2, \ 8230;, D }, to ensure the variation vector v i,g At least one dimension of information is preserved; the crossover probability CR is a constant within the interval (0, 1).
Step 5, one-to-one tournament selection is performed between the test vector and the corresponding target vector, which specifically comprises the following steps:
for test vector u i,g And a target vector x i,g Comparing the target function values of the test vectors, and if the test vectors have more optimal target functions, replacing the target vectors with the test vectors; otherwise, the target vector remains unchanged; taking the objective function minimization as an example, the following formula is used for representation:
Figure BDA0002412182790000041
updating the individuals by using the success probability in the step 6 specifically comprises the following steps:
7-1, calculating chromosomes of each individual by using the formula (6);
x′ i =x i +N(0,σ) (6)
where N (0, σ) represents the mean value of 0, σ is a variable, and x i Denotes the individual before update, x' i Representing the updated individual;
7-2 updating the adaptive mutation probability P s
P s The updating method comprises the following steps:
P s =(1-β)P s +β (7)
wherein beta is-1/N;
7-3 update sigma, and three update modes are provided;
if P s Less than 1/3, update σ using equation (8) if P s If the ratio is more than 1/3, updating by using the formula (9), otherwise, keeping unchanged;
σ=C*σ 0 (8)
σ=C/σ 0 (9)。
and 7, repeating the steps 3 to 6, and outputting the optimal position value when the maximum iteration number is reached, wherein the method specifically comprises the following steps:
and (3) judging whether the maximum iteration number is reached, if the maximum iteration number is not reached, executing the step 3-step 6, guiding chromosome evolution by using the formula (2) as a fitness function, and obtaining a code with the maximum function value as the best estimation of the position when the maximum iteration number is reached.

Claims (2)

1. A WSN node positioning method based on a differential evolution genetic algorithm is characterized by comprising the following steps:
step 1, establishing a node positioning optimization model and designing a positioning optimization objective function;
step 2, initializing a population;
step 3, generating a variation vector by utilizing differential variation operation after the population is subjected to initialization assignment;
step 4, after performing mutation operation, the target vector x needs to be transformed i,g And the variation vector v i,g Generating final test vector u by carrying out binomial intersection i,g =[u i1,g ,u i2,g ,…,u iD,g ];
Step 5, one-to-one tournament selection is executed between the test vector and the corresponding target vector;
step 6, updating the individuals by using the success probability;
step 7, repeating the steps 3 to 6, and outputting an optimal position value when the maximum iteration times are reached;
establishing a node positioning optimization model and designing a positioning optimization objective function in the step 1, specifically, the following steps are performed:
anchor node P 1 (x 1 ,y 1 ),P 2 (x 2 ,y 2 ),…,P n (x n ,y n ) And the actual true lengths of the unknown nodes P (x, y) are r 1 ,r 2 ,…,r n The difference of the distance measurement errors is epsilon 12 ,…,ε n Then | r is satisfied i -d i |<ε i Where i =1,2, \8230, n, then unknown node P (x, y) constraints:
Figure FDA0004071745090000011
(x, y) is solved such that:
Figure FDA0004071745090000012
converting the node positioning problem into an optimization problem for solving the minimum value of the ranging error, and taking a formula (2) as a fitness function;
the initialization population in step 2 is specifically as follows:
assigning default values to the chromosome individuals within the limit range, calculating function values Fitness (i) of all chromosomes in the total group according to a formula (2), and determining the individual p with the highest Fitness value in the current algebraic individuals g And the position p with the highest individual cutoff algebraic fitness value h
After the population stated in step 3 is subjected to initialization assignment, generating a variation vector v by utilizing differential variation operation i,g The method comprises the following steps:
v i,g =x r1,g +F*(x r2,g -x r3,g ) (3)
in the formula (3), r 1 ,r 2 ,r 3∈[1,2,…NP] The scaling factor F is a constant in the interval of (0, 1) and is used for controlling the size of the differential vector; x is the number of r1,g 、x r2,g And x r3,g All are a vector randomly selected from the population; after performing mutation operation as described in step 4, the target vector x is required to be transformed i,g And the variation vector v i,g Generating final test vector u by carrying out binomial intersection i,g =[u i1,g ,u i2,g ,…,u iD,g ]The method comprises the following steps:
Figure DA00040717450938732231
in the formula, j rand Is randomly selected from the set {1,2, ..., D }An integer is selected to ensure the variation vector v i,g At least one dimension of information is preserved; the crossover probability CR is a constant within the interval (0, 1);
step 5, one-to-one tournament selection is performed between the test vectors and the corresponding target vectors, which is specifically as follows:
for test vector u i,g And a target vector x i,g Respectively solving the fitness values by using a formula (2) for comparison, and replacing the target vector with the test vector if the test vector has a better value; otherwise, the target vector remains unchanged; expressed by the following formula:
Figure FDA0004071745090000022
updating the individuals by using the success probability in the step 6 specifically comprises the following steps:
7-1, calculating chromosomes of each individual by using the formula (6);
x′ i =x i +N′(0,σ) (6)
where N' (0, σ) represents mean 0, σ as variable, x i Denotes the individual before update, x' i Representing the updated individual;
7-2 updating the adaptive mutation probability P s
P s The updating method comprises the following steps:
P s =(1-β)P s +β (7)
wherein beta is-1/N;
7-3 update sigma, and three update modes are provided;
if P s Less than 1/3, update σ using equation (8) if P s If the ratio is more than 1/3, updating by using a formula (9), otherwise, keeping unchanged;
σ=C*σ 0 (8)
σ=C/σ 0 (9)。
2. the method for positioning a WSN node based on a differential evolution genetic algorithm as claimed in claim 1, wherein the step 7 is repeated from step 3 to step 6, and when the maximum number of iterations is reached, the optimal position value is output, specifically as follows:
and (3) judging whether the maximum iteration number is reached, if the maximum iteration number is not reached, executing the step 3-step 6, guiding chromosome evolution by using the formula (2) as a fitness function, and obtaining a code with the maximum function value as the best estimation of the position when the maximum iteration number is reached.
CN202010180046.1A 2020-03-16 2020-03-16 WSN node positioning method based on differential evolution genetic algorithm Active CN111447627B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010180046.1A CN111447627B (en) 2020-03-16 2020-03-16 WSN node positioning method based on differential evolution genetic algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010180046.1A CN111447627B (en) 2020-03-16 2020-03-16 WSN node positioning method based on differential evolution genetic algorithm

Publications (2)

Publication Number Publication Date
CN111447627A CN111447627A (en) 2020-07-24
CN111447627B true CN111447627B (en) 2023-04-18

Family

ID=71648700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010180046.1A Active CN111447627B (en) 2020-03-16 2020-03-16 WSN node positioning method based on differential evolution genetic algorithm

Country Status (1)

Country Link
CN (1) CN111447627B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115372903B (en) * 2022-07-28 2024-03-08 国网江苏省电力有限公司电力科学研究院 Near-ultrasonic-based long and narrow pipeline personnel positioning method and device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201717125D0 (en) * 2016-11-28 2017-11-29 National Univ Of Defense Technology Differential evolution method oriented to agile satellite multi-target task planning
CN109379780A (en) * 2018-10-23 2019-02-22 华南理工大学 Wireless sensor network locating method based on adaptive differential evolution algorithm
CN109461471A (en) * 2018-08-30 2019-03-12 浙江工业大学 A kind of adaptive Advances in protein structure prediction based on championship mechanism
CN109640261A (en) * 2018-11-22 2019-04-16 天津理工大学 A kind of location algorithm of the heterogeneous wireless sensor network based on support vector regression
CN109688537A (en) * 2018-12-21 2019-04-26 延安大学 Forest fire monitoring Internet of things node localization method based on differential evolution algorithm
CN110839245A (en) * 2019-10-31 2020-02-25 杭州电子科技大学 Wireless sensor network node deployment method applied to indoor positioning

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201717125D0 (en) * 2016-11-28 2017-11-29 National Univ Of Defense Technology Differential evolution method oriented to agile satellite multi-target task planning
CN109461471A (en) * 2018-08-30 2019-03-12 浙江工业大学 A kind of adaptive Advances in protein structure prediction based on championship mechanism
CN109379780A (en) * 2018-10-23 2019-02-22 华南理工大学 Wireless sensor network locating method based on adaptive differential evolution algorithm
CN109640261A (en) * 2018-11-22 2019-04-16 天津理工大学 A kind of location algorithm of the heterogeneous wireless sensor network based on support vector regression
CN109688537A (en) * 2018-12-21 2019-04-26 延安大学 Forest fire monitoring Internet of things node localization method based on differential evolution algorithm
CN110839245A (en) * 2019-10-31 2020-02-25 杭州电子科技大学 Wireless sensor network node deployment method applied to indoor positioning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑建国等.混合粒子群和差分进化的定位算法.计算机测量与控制.2019,第第27卷卷(第第27卷期),192-195. *

Also Published As

Publication number Publication date
CN111447627A (en) 2020-07-24

Similar Documents

Publication Publication Date Title
CN109597864B (en) Method and system for real-time positioning and map construction of ellipsoid boundary Kalman filtering
CN109581281B (en) Moving target positioning method based on arrival time difference and arrival frequency difference
CN109379780B (en) Wireless sensor network positioning method based on self-adaptive differential evolution algorithm
CN105203106A (en) WMPS network layout optimizing method based on simulated annealing particle swarm optimization
CN110225454B (en) Confidence transfer distributed type volume Kalman filtering cooperative positioning method
CN109039720B (en) Parameter estimation method for binary sensor network
CN111447627B (en) WSN node positioning method based on differential evolution genetic algorithm
CN113326912A (en) Information sharing Harris eagle optimization-based ultra-wideband positioning method
CN101917762A (en) Node positioning method of particle swarm optimization sensor with penalty function
CN115907304A (en) Evaluation index weight optimization method based on inter-zone analytic method
Wang et al. Research on hybrid PSODE with three populations based on multiple differential evolutionary models
CN115775038A (en) Short-term load prediction method based on IGWO optimization LSSVM
CN109253727B (en) Positioning method based on improved iteration volume particle filter algorithm
CN113993152A (en) Communication base station flow prediction method based on improved Husky algorithm optimization support vector machine regression machine
CN113432608A (en) Generalized high-order CKF algorithm based on maximum correlation entropy and suitable for INS/CNS integrated navigation system
CN112529328A (en) Product performance prediction method and system
CN115550837B (en) DV-Hop positioning method based on chaotic mapping and gray wolf algorithm optimization
CN116865255A (en) Short-term wind power prediction method based on improved entropy weight method and SECEEMD
CN111259340A (en) Saturated load prediction method based on logistic regression
CN113361146B (en) Improved particle swarm optimization-based manganese-copper shunt structure parameter optimization method
Chen et al. Improved DV-Hop Node location Optimization Algorithm Based on Adaptive Particle Swarm
CN114637956A (en) Novel double-Kalman filtering method
CN109753695B (en) Method for expanding characteristics of engine parts
CN106961656B (en) Wireless sensor network data prediction method
CN112037289A (en) Off-line parameter calibration method based on genetic algorithm

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
GR01 Patent grant