Estimation errors for Model I and Model II in different simulation runs of the one dimensional localization problem.

## 1. Introduction

In the past years, the development in micro electromechanical systems (MEMS), radio frequency (RF), integrated circuit (IC), etc., greatly enhanced the advancement of wireless sensor networks (WSNs). As an ubiquitous sensing technology, WSNs find more and more applications, such as structural monitoring [34], precision agriculture [3], gas-leak localization [14], volcano monitoring [33], robot navigation [4, 15], health monitoring [20], to name a few. For most existing applications of WSNs, the location information is crucial. For example, in the structural monitoring application, we can conclude that the structure is out of condition if fault is detected by one or more sensors in the network of sensors mounted everywhere on the structure. However, we are unable to accurately report the faulty position without localization capability of the WSN. In contrast to other type of networks, e.g., Internet, a prominent difference is that WSNs are location-based networks. Therefore, the design of localization hardware and localization algorithms is an important procedure in the development of a WSN system.

There are mainly two classes of localization approaches for WSNs: one is pre-localization and the other one is self-localization. The pre-localization method measures the position of sensors in the deployment stage. After the deployment and position measurement, the position is stored in the memory of the sensor. For this method, any movement of the sensors will result in errors in the location information. Differently, the self-localization method computes the locations of each sensor based on real-time measurements and therefore is robust to the variance of the environment. With GPS devices embedded, sensors are enabled with self-localization capability. However, the relatively high cost of GPS devices often makes it not practical to apply GPS to all sensors in a network. Instead, the strategy with a portion of sensors equipped with GPS as beacons and using triangulation or trilateration to iteratively determine the positions of blind sensors based on the distance or angle measurements between neighboring sensors provides a less expensive way for self-localization [13, 16, 26]. Although many GPS devices are saved, as a tradeoff, the sensors are required to have the ability to measure the distance or the relative angle to its neighbor, which may result in costs for extra hardware. Without introducing extra hardware, received signal strength (RSS) based distance measurement method [17, 30], relying on the estimated distance according to the signal strength received from the neighboring sensor, provides a promising direction for self-localization. Another promising self-localization method is range-free localization, which even does not require the information on the signal strength received from the neighbor but the connectivity information, i.e., a sensor only need to know who is its neighbor. This technology implies that localization can be a by-product of communication since connectivity information can be obtained in communication. For example, if Sensor

Dynamic models gained great success in realtime signal processing [28], robotics [12, 22], online optimization [29], etc.. In this chapter, we overview our previous work on dynamic model based range-free localization [10, 11, 25]. Particularly, we will examine two dynamic models for the real time localization of WSNs. The models are described by nonlinear ordinary differential equations (ODEs). The state value of the ODEs converges to the expected position estimation of sensors. Both of the two models find feasible solutions to the formulated optimization problem. Particularly, the second model, by exploiting heuristic information, has a tendency to converge to better solutions in the sense of localization error. The real time processing ability of the models allows possible movement of the sensor nodes, which often happens in mobile sensor networks [23]. Besides the real time localization capability, another prominent feature of the proposed models is that both of them are completely distributed, i.e., each sensor in the network only need to exchange information with its neighbor and thus no message passing is needed in the network. This advantage makes the proposed algorithms scalable to large scale networks involving thousands of sensors or more.

The remainder of this paper is organized as follows. In Section 2, some preliminaries on range-free localization of WSNs are presented. In Section 3, we formulate the localization problem from an optimization perspective. Two dynamic models are presented in Section 4 to solve the formulated optimization problems. In Section 5, simulations are performed to demonstrate the effectiveness of our method. Section 6 concludes this paper.

## 2. Preliminaries

In this chapter, we assume that all sensors are equipped with communication modules and the locations of beacon sensors are known. Fig. 1 sketches the connectivity topology of a WSN consisting of beacon sensors and blind sensors. In the network, The beacon sensors are those with known locations. The locations can be obtained either by GPS or by pre-deployment. The blind sensors are those without pre-known positions. A sensor can communicate with other sensors within the signal coverage area. The communication links and sensors therefore form a network with sensors as nodes and communication links as edges.

The signal strength at a given distance from the emitter varies due to propagation conditions, material coverage, antenna configurations and battery conditions [31] and the calculated distance according to the received signal strength often has a large error [8, 18]. Nevertheless, the nominal maximum range, which is measured under ideal conditions in open environments without obstacles along the signal propagation route, without material coverage, with a proper configuration of the antenna and with a full power of the battery, etc., gives an upper bound on the distance between the emitter and receiver pair. As depicted in Fig. 2. Subfigure (a) depicts the ideal open environment, where the communication radius, denoted by

## 3. Problem formulation

In this section, we present the mathematical formulation of the problem.

### 3.1. Nonlinear inequality problem formulation

As discussed in Section 2, the position of beacon sensors are known and the distance between two neighbor sensor (in the sense of communication) is less than

where

**Remark 1** *There is no explicit objective function but inequality and equality constraints in problem (1). The solution to this problem is generally not unique. We are more concerned with finding a feasible solution in real time instead of finding all the feasible solutions. Based on this consideration, we explore finding a feasible solution to problem (1) in real time via a dynamic model.*

### 3.2. Optimization problem formulation

To find a feasible solution of problem (1) numerically, we first transform the problem into an optimization problem and employ dynamic evolutions to solve it.

The solution of problem (1) is identical to the one of the following normal optimization with an explicit objective function,

where

## 4. Solving the problem via nonlinear dynamic evolution

In this section, we present two ODE models, both of which are able to solve the range-free localization problem (2). As the solution to the problem is generally not unique. Property employment of heuristic information may improve the solution performance. Based on the feasible solution obtained by the first ODE model, the second ODE model proposed in this chapter indeed realizes the improvement in performance.

### 4.1. Model I

The partial sub-gradient relative to

where

In the ODE model, each blind sensor is associated with a dynamic module. The modules interact with their neighbor modules and all the modules together perform the localization task and solve the problem (1). The dynamic evolution of

Notably, the ODE model (3) is a distributed one. Communication only happens between neighboring sensors. No routing or cross-hop communication is required for the implementation of the ODE model. The distributed nature of the model thoroughly reduces the communication burden and makes the method scalable to a network with a large number of sensors involved.

About the ODE model I (3), we have the following theorem,

**Theorem 1 [25]** *The ODE model I (3) withε1>0, wijfor all possible iandj, asymptotically converges to a feasible solution xi* (for all iin the blind sensor set) of problem (1).*

The proof of this theorem is based on Lyapunov stability theory. Interested readers are refereed to our previous work [10, 25] for a detailed proof. This theorem reveals that the ultimate output of the ODE model I is a feasible solution to problem (1).

### 4.2. Model II

Section 4.1 provides an ODE model to find a feasible solution of the problem. The model presented in this part is also a dynamic ODE model. Different from Model I, which is initialized randomly and does not use any heuristic information, Model II is initialized with the ultimate output of Model I with

where

We use the following gradient based dynamics to solve (5):

where

The ODE model (6) is a distributed one since the update of

About the initialization of the ODE model, we have the following remark,

**Remark 2** *The ODE model (6) is initialized with the ultimate output of the ODE model I (3) with Rreplaced by R-δin (4) withδ≪R. The goal is to ensure the ultimate output of Model I strictly locates inside the open set formed by (1), which is necessary for the barrier term in (5) to restrict the solution always stays inside the feasible region.*

According to Theorem 1, the ODE model I with

with which we conclude that

**Theorem 2 [11]** The ODE model II (6) with

## 5. Simulations

In this section, simulations are used to verify the two ODE models in both the one dimensional space and the two dimensional space.

### 5.1. Range-free localization in one dimensional spaces

In this part, we investigate the range-free localization of sensors in a network deployed in a one dimensional topology.

#### 5.1.1. Background

There are a bunch of applications which deploys sensors along an one dimensional line. For example, WSNs for highway monitoring [1, 9, 21] are often deployed along the highway direction and thus form a one dimensional deployment topology, as sketched in Fig. 3. Other applications, such as WSNs for bridge health monitoring [32] and WSNs along a tunnel [5] for traffic safety, can also be put into the category of one dimensional localization problem.

#### 5.1.2. Simulation setup and simulation results

We consider a wireless sensor network with one dimensional deployment. There are

For the dynamic models, the state values of Model I are randomly initialized. We choose

Estimation Error of Model I | Estimation Error of Model II | |||

E1 | E2 | E1 | E2 | |

1 | 41.2948 | 84.0870 | 9.9419 | 23.7890 |

2 | 44.5053 | 84.4876 | 10.2123 | 27.2300 |

3 | 43.7765 | 80.0258 | 13.2497 | 26.1105 |

4 | 44.3971 | 83.3457 | 11.1692 | 26.9130 |

5 | 40.9815 | 83.3549 | 12.4231 | 28.3951 |

6 | 40.7992 | 79.4815 | 10.8819 | 25.6199 |

7 | 39.5168 | 86.9493 | 11.4397 | 26.0199 |

8 | 37.4373 | 65.4550 | 10.4166 | 24.0404 |

9 | 42.9898 | 83.3865 | 11.4330 | 26.6570 |

10 | 44.0201 | 78.4496 | 13.2752 | 28.7782 |

### 5.2. Range-free localization in two dimensional Spaces

In this part, we investigate the range-free localization of sensors in a network deployed in a two dimensional topology.

#### 5.2.1. Background

Most existing literatures deal with the general localization problem in two dimensional space. The one dimensional sensor localization problem investigated in the last section falls into this category by fixing the value of sensor positions along one dimension. Localization in applications, such as wildlife monitoring [7], WSN aided robot navigation [6, 19] and animal tracking [27], can be abstracted as two dimensional localization problems.

#### 5.2.2. Simulation setup and simulation results

In the simulation, we consider a

As to the dynamic models, we choose the scaling factors

Fig. 9 and Fig. 10 show the time histories of the position estimations by Model I along x-direction and that along y-direction respectively. From the figures we can observe that after a short period of transient, the estimation results converge to constant values. The time histories of the position estimations along x-direction and y-direction estimations by Model II are plotted in Fig. 11 and Fig. 12, respectively. Compared to the transient of Model I, the change of state values in Fig. 11 and Fig. 12 are much milder. The adjustment of values refine the initial estimation with the tendency to even distributions under the communication connectivity constraints. With time elapses, the estimation results by Model II converge. It is worth noting that the ultimate values by Model II shown in Fig 11 and Fig. 12 are not exactly uniformly distributed. This is due to the compromise of the heuristic information driving to even distribution in space and the inequality constraints imposed by communication connectivity. The final position estimations by Model I and Model II are shown in Fig. 13 and Fig. 14, respectively. It can be observed that the estimated results shown in both figures are within the area covered by the circle centered at the true position with a radius

Estimation Error of Model I | Estimation Error of Model II | |||

E1 | E2 | E1 | E2 | |

1 | 19.2071 | 34.9071 | 10.9093 | 24.0850 |

2 | 20.6657 | 42.7582 | 10.2647 | 18.3666 |

3 | 21.2873 | 37.9441 | 10.2106 | 18.2110 |

4 | 20.0403 | 41.3617 | 9.9603 | 16.8761 |

5 | 20.6692 | 39.0069 | 9.1838 | 17.2634 |

6 | 18.1865 | 40.7181 | 9.5004 | 21.8774 |

7 | 21.0053 | 37.7784 | 10.2448 | 18.6334 |

8 | 23.1289 | 44.2787 | 12.3591 | 25.5854 |

9 | 23.5922 | 40.3665 | 10.6922 | 21.0049 |

10 | 19.5115 | 41.8718 | 11.6586 | 24.0909 |

### 5.3. Discussions

In the above two subsections, we considered the range-free localization problem in one dimensional case and two dimensional case respectively. In some applications of WSNs, higher dimensional cases [24] (see Fig. 15 for the sketch of a typical three dimensional one after normalization along three axial directions), such as building monitoring [35], underwater acoustic sensor networks[2], may be encountered. For example, the localization problem of sensors for building monitoring is actually defined in a three dimensional space since sensors are deployed in two dimensions on each floor and the whole network constructed by sensors on different floors forms a three dimensional one. Also, sensors in the underwater acoustic sensor networks are often deployed at different depth, with different longitude and latitude and thus form a three dimensional sensor network.

The presented two models in this chapter admit the higher dimensional localization problems as we did not specify the number of dimensions in the problem formulation and the model works for all possible dimensions.

As demonstrated in simulations, Model II outperforms Model I in the sense of estimation error for the cases with the same simulation setup. However, it is notable that Model II requires an extra dynamic model for the initialization, which is at the cost of implementation complexity and longer computation time. Fortunately, the dynamic models can be implemented with either digital or analog devices and thus the computation can be completed in a very short time. For example, the simulation examples for the two dimensional localization problem showed that it takes a time of

## 6. Conclusions and future work

In this chapter, we overviewed our recent work on range-free localization of sensors in WSNs via dynamic models. The range-free localization problem is formulated as two different optimization problems, each of which corresponds to a dynamic model, namely Model I and Model II, for the solution. Simulations in both one dimensional case and two dimensional case are performed and the two models are compared in both sceneries. The simulation results demonstrate effectiveness of the dynamic models.

Compared with conventional range-free localization algorithms, a prominent advantage of ODE based solutions is that the models are implementable by parallel hardware. As a promising direction for WSN localization, the following aspects are of fundamentally importance in both theory and practise and are open to all researchers,

Is there better heuristic information applicable to improve the performance?

Can Model II be modeled from a probabilistic perspective and therefore to evaluate its performance improvement relative to Model I in quantity?

Recent advancement in ODE shows that a properly designed ODE model can receive finite time convergence to its steady state value. Can the results applied to the range-free WSN localization to gain a theoretically finite time convergence?

Hardware implementation with VLSI, FPGA, etc. is a bright direction in practise.

Can the presented models be used in range-base WSN localization, if not applicable directly, is the results helpful to supply a good "warm start", which accelerates the convergence?

The presented models are essentially ODE models for solving nonlinear inequalities defined on a network. Are the models extendable to other network applications, such as robotic networks, smart grids, the internet of things, etc?

### Acknowledgement

The authors would like to acknowledge the constant motivation by the fact presented in the movie Gattaca that never saving anything for the swim back enables Vincent to swim longer and faster than his genetically superior brother. The authors would like to acknowledge the support by the National Science Foundation of China under Grant NSFC: 61105090.