# A WLS Estimator for Target Localization in a Cooperative Wireless Sensor Network

Slavisa Tomic, Marko Beko, Rui Dinis, Milan Tuba

► To cite this version:

Slavisa Tomic, Marko Beko, Rui Dinis, Milan Tuba. A WLS Estimator for Target Localization in a Cooperative Wireless Sensor Network. 7th Doctoral Conference on Computing, Electrical and Industrial Systems (DoCEIS), Apr 2016, Costa de Caparica, Portugal. pp.273-283, 10.1007/978-3-319-31165-4_27. hal-01438252

HAL Id: hal-01438252

https://inria.hal.science/hal-01438252

Submitted on 17 Jan 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A WLS Estimator for Target Localization in a Cooperative Wireless Sensor Network

Slavisa Tomic<sup>4</sup>, Marko Beko<sup>1,3</sup>, Rui Dinis<sup>2,5</sup>, Milan Tuba<sup>6</sup>

<sup>1</sup> Universidade Lusófona de Humanidades e Tecnologias, Lisbon, Portugal

<sup>2</sup> Instituto de Telecomunicações, Lisbon, Portugal

<sup>3</sup> CTS, UNINOVA - Campus FCT/UNL, Caparica, Portugal

<sup>4</sup> ISR/IST, LARSyS, Lisbon, Portugal

<sup>5</sup> DEE/FCT/UNL, Caparica, Portugal

<sup>6</sup> Megatrend University, Faculty of Computer Science, Belgrade, Serbia
{s.tomic@campus.fc.unl.pt, mbeko@uninova.pt, rdinis@fct.unl.pt, tuba@ieee.org}

**Abstract.** This paper addresses target localization problem in a cooperative 3-D wireless sensor network (WSN). We employ non-traditional methodology which merges distance and angle measurements, respectively withdrawn from the received signal strength (RSS) and angle-of-arrival (AoA) information. Based on RSS measurement model and effortless geometry, a novel non-convex estimator according to the weighted least squares (WLS) criterion is obtained, which closely approximates the maximum likelihood (ML) estimator for small noise. It is shown that the devised estimator is appropriate for distributed implementation. Following the squared range (SR) approach, we propose a suboptimal SR-WLS estimator according to the generalized trust region sub-problem (GTRS) framework, to estimate the locations of all targets in the WSN. According to our simulations, the new estimator has excellent performance in a great variety of considered settings, in which the effectiveness of fusing two radio measurements is confirmed.

**Keywords:** Angle-of-arrival (AoA), generalized trust region sub-problem (GTRS), received signal strength (RSS), wireless sensor network (WSN).

## 1 Introduction

Accurate information about sensor's location is a key component in many practical applications. Wireless localization algorithms usually rely upon range measurements extracted from angle-of-arrival (AoA), received signal strength (RSS), time-of-arrival (ToA) information, or a combination of them [1].

In [2], a hybrid methodology which fuses distance and angle measures has been considered. Both estimators proposed in [2] deal with the non-cooperative target localization problem in a 3-D space: linear least squares (LS) and optimization based. The former one is a fairly simple estimator, and the later one uses Davidson-Fletcher-Powell algorithm [3]. The authors in [4] derived an LS and a maximum likelihood (ML) estimator for a hybrid scheme that fuses RSS difference (RSSD) and AoA measurements. The authors in [4] employed non-linear constrained optimization to estimate the unknown location from multiple RSS and AoA measurements. In [5], a
selective weighted LS (WLS) estimator for mixed RSS/AoA localization problematic
was presented. The target position was established by employing loaded distances from
the two closest anchor measures, together with the serving base station AoA
measurement. Another WLS estimator for a non-cooperative localization problem for
the case where the transmitted power is not known was proposed in [6]. Nonetheless,
similar as the method proposed in [5], the WLS method has been derived without
requiring information about the statistical properties of RSS and AoA measurements,
which might lead to significant performance degradation in practice. Also, the authors
in [6] investigated a small-scale wireless sensor network (WSN) only, with extremely
low noise power.

All of the above approaches examine non-collaborative localization problem only, where the location of a single target, which communicates with anchors exclusively, is established at a time. Contrarily to the above mentioned approaches, here, the target localization problem in an extensive WSN is considered, where the number of anchors is insufficient and the communication range of all sensors is limited (e.g., to prolong sensor's battery life). Hence, only some targets can directly communicate with anchors and sensor cooperation is required in order to acquire adequate amount of information to carry out the localization. Through employing the RSS spreading model and straightforward geometry, we first develop a new local non-convex estimator based on the WLS criterion that closely approximates the local ML one in the case where noise is low. Next, by following the squared range (SR) approach, we propose a suboptimal SR-WLS estimator based on the generalized trust region sub-problem (GTRS) framework, which is possible to solve *exactly* with a bisection method [7]. To the best of the authors' knowledge, distributed localization algorithms for hybrid RSS/AoA systems in cooperative WSNs are yet to be published.

## 2 Relationship to Cyber-Physical Systems

Cyber-physical systems (CPSs) have recently attracted much attention from both the educational and engineering public. These schemes are displaying an enormous potential in interacting with the physical world and creating its management more efficient. Integrating the computation and communication competencies into the modules of the physical surrounding [8] enables this. CPSs represent nowadays the new peer group of networks and embedded structures.

The placement of CPS raises numerous difficulties, where the most significant topic that has activated massive quantity of study is localization [8]. Basically, localization goals at estimating the location of CPS modules. In CPSs where data are closely related with the surroundings and the location at which they are produced, localization is a fundamental task.

Collaboration of a large number of scattered sensors in a WSN can be seen as a CPS.
In such systems, localization is of a crucial importance, since a system may be
configured to react locally to the variations within sensor records; hence, accurate
determination of the location where the deviations arise is the key. Furthermore,
services based on location-awareness represent a key component of countless wireless
structures nowadays. By exploring the synergies between computational and physical
components we can form smart environments which offer improved safety and efficiency in everyday life, e.g. smart parking, assistance for elderly or people with disabilities, monitoring of storage conditions and goods, etc.

### 3 Problem Statement

We examine a large-scale WSN with $N$ anchors and $M$ targets, that may also be viewed as a connected graph, $G(V, E)$, with $|V| = M + N$ vertices and $|E|$ edges (connections), where $|\cdot|$ represents the cardinality of a set. The set of targets and the set of anchors are respectively denoted as $\mathcal{T}$ ($|\mathcal{T}| = M$) and $\mathcal{A}$ ($|\mathcal{A}| = N$), and their locations are denoted by $\mathbf{x}_1, \mathbf{x}_2, \cdots, \mathbf{x}_M$ and $\mathbf{a}_1, \mathbf{a}_2, \cdots, \mathbf{a}_N$, ($\mathbf{x}_i, \mathbf{a}_j \in \mathbb{R}^3, \forall i \in \mathcal{T}$ and $\forall j \in \mathcal{A}$), respectively. Due to limited communication range, $R$, two sensors, $i$ and $j$, can interchange data if and only if they are inside the communication range of each other. The sets of all target/anchor and target/target edges are defined as $\mathcal{E}_{\mathcal{A}} = \{(i, j): \|\mathbf{x}_i - \mathbf{a}_j\| \le R, \forall i \in \mathcal{T}, \forall j \in \mathcal{A}\}$ and $\mathcal{E}_{\mathcal{T}} = \{(i, k): \|\mathbf{x}_i - \mathbf{x}_k\| \le R, \forall i, k \in \mathcal{T}, i \ne k\}$, respectively.

<figure><img src="image_2.png" /><figcaption>Fig. 1. Representation of a target and anchor locations within 3-D.</figcaption></figure>

In favor of ease of expression, we describe a matrix $\mathbf{X} = [\mathbf{x}_1, \cdots, \mathbf{x}_M] (\mathbf{X} \in \mathbb{R}^{3 \times M})$ as the matrix of all unknown target locations. These unknown locations are estimated with the means of a hybrid methodology which merges distance and angle measures.

Throughout this work, we imply that the distance measurements are extracted from the RSS measurements, because ranging based on RSS does not oblige additional hardware [1]. The path loss between two sensors is given by:

$$L_{ij} = L_0 + 10\gamma \log_{10} \frac{\|\mathbf{x}_i - \hat{\mathbf{a}}_j\|}{d_0} + n_{ij}, \forall (i, j) \in \mathcal{E}_{\mathcal{A}} \cup \mathcal{E}_{\mathcal{T}} \quad (1)$$

(see [9],[10]) where $L_0$ represents the path loss value at a reference distance $d_0$ ($\|\mathbf{x}_i - \hat{\mathbf{a}}_j\| \ge d_0$), $\gamma$ is the path loss exponent (PLE), $\hat{\mathbf{a}}_j$ is the j-th neighbor of the i-th target ($\mathbf{a}_j$ if $j \in \mathcal{A}, \mathbf{x}_j$ if $j \in \mathcal{T}$), and $n_{ij}$ is the log-normal shadowing considered as a
Gaussian random variable, i.e., $n_{ij} \sim \mathcal{N}(0, \sigma_{n_{ij}}^2)$. The target/target path loss measures are assumed to be symmetric<fnref n="1" />.

**Fig. 1** gives a representation of a target and anchor locations within 3-D. Looking at **Fig. 1**, $\mathbf{x}_i = [x_{i1}, x_{i2}, x_{i3}]^T$ and $\mathbf{a}_j = [a_{j1}, a_{j2}, a_{j3}]^T$ represent the coordinates of the i-th target (not known) and the coordinates of the j-th anchor (known) respectively, whereas $d_{ij}$, $\phi_{ij}$ and $\alpha_{ij}$ respectively stand for the range, azimuth and elevation angle amongst the i-th target and the j-th anchor. According to the measurement model (1) the ML distance between two sensors is obtained as follows [1]:

$$ \hat{d}_{ij} = d_0 10^{\frac{L_{ij}-L_0}{10\gamma}}, \forall (i, j) \in \mathcal{E}_A \cup \mathcal{E}_T. \quad (2) $$

By applying simple geometry, azimuth and elevation angle measurements can be modeled respectively as [2]:

$$ \phi_{ij} = \tan^{-1} \left( \frac{x_{i2} - a_{j2}}{x_{i1} - a_{j1}} \right) + m_{ij}, \forall (i, j) \in \mathcal{E}_A, \quad (3) $$

and

$$ \alpha_{ij} = \cos^{-1} \left( \frac{x_{i3} - a_{j3}}{\|x_i - a_j\|} \right) + v_{ij}, \forall (i, j) \in \mathcal{E}_A, \quad (4) $$

According to (1), (3) and (4), we can obtain the ML estimator as [11]:

$$ \hat{\mathbf{X}} = \arg \min_{\mathbf{X}} \sum_{i=1}^{3|\mathcal{E}_A|+|\mathcal{E}_T|} \frac{1}{\sigma_i^2} [\theta_i - f_i(\mathbf{X})]^2, \quad (5) $$

where $\theta = [\mathbf{L}^T, \mathbf{\phi}^T, \mathbf{\alpha}^T]^T$ ($\theta \in \mathbb{R}^{3|\mathcal{E}_A|+|\mathcal{E}_T|}$), with $\mathbf{L} = [L_{ij}]^T, \forall (i, j) \in \mathcal{E}_A \cup \mathcal{E}_T, \mathbf{\phi} = [\phi_{ij}]^T, \forall (i, j) \in \mathcal{E}_A, \mathbf{\alpha} = [\alpha_{ij}]^T, \forall (i, j) \in \mathcal{E}_A$, and

$$ f(\mathbf{X}) = \left[ \begin{array}{c} \vdots \\ L_0 + 10\gamma \log_{10} \frac{\|\mathbf{x}_i - \hat{\mathbf{a}}_j\|}{d_0} \\ \vdots \\ \tan^{-1} \left( \frac{x_{i2} - a_{j2}}{x_{i1} - a_{j1}} \right) \\ \vdots \\ \cos^{-1} \left( \frac{x_{i3} - a_{j3}}{\|\mathbf{x}_i - \mathbf{a}_j\|} \right) \\ \vdots \end{array} \right], \sigma = \left[ \begin{array}{c} \vdots \\ \sigma_{n_{ij}} \\ \vdots \\ \sigma_{m_{ij}} \\ \vdots \\ \sigma_{v_{ij}} \\ \vdots \end{array} \right]. $$

The non-convex LS estimator in (5) does not have a closed-form solution. Nonetheless, we will show that it can be transformed into GTRS framework and solved efficiently.

<footnote n="1">¹ This assumption is made without loss of generality; it is readily seen that, if $L_{ij}^T \neq L_{ji}^T, \forall (i, j) \in \mathcal{E}_T$, then it is enough to replace $L_{ij}^T \leftarrow (L_{ij}^T + L_{ji}^T)/2$ and $L_{ji}^T \leftarrow (L_{ij}^T + L_{ji}^T)/2$ when solving the localization problem.</footnote>
## 4 The Proposed SR-WLS Estimator

Assuming that the initial location estimations of the targets, $\hat{\mathbf{X}}^{(0)}$, are given, the ML in (5) can be solved locally per target, resorting merely to the data collected from its neighbors and employing an iterative methodology. Consequently, target $i$ updates its location estimation in each iteration, $t$, by solving the below local ML estimator:

$$ \hat{\mathbf{x}}_{i}^{(t+1)} = \arg \min_{\mathbf{x}_{i}} \sum_{j=1}^{3|\varepsilon_{A_{i}}|+|\varepsilon_{T_{i}}|} \frac{1}{\sigma_{j}^{2}} [\theta_{j} - f_{j}(\mathbf{x}_{i})]^{2}, \forall i \in \mathcal{T}, \quad (6) $$

where $\mathcal{E}_{A_i} = \{j: (i, j) \in \mathcal{E}_A\}$ and $\mathcal{E}_{T_i} = \{k: (i, k) \in \mathcal{E}_T\}$ represent the set of all anchor and all target neighbors of the target $i$ respectively.

Given $\hat{\mathbf{X}}^{(0)}$, for the case where noise power is low, out of (1) we can write:

$$ \lambda_{ij} \|\mathbf{x}_i - \hat{\mathbf{a}}_j\|^2 \approx d_0^2, \forall i \in \mathcal{T}, \forall j \in \mathcal{E}_{A_i} \cup \mathcal{E}_{T_i}, \quad (7) $$

where $\lambda_{ij} = 10^{\frac{L_0-L_{ij}}{5\gamma}}$. Similarly, from (3) and (4) we respectively get:

$$ \mathbf{c}_{ij}^{T}(\mathbf{x}_{i} - \mathbf{a}_{j}) \approx 0, \forall i \in \mathcal{T}, \forall j \in \mathcal{E}_{A_{i}}, \quad (8) $$

and

$$ \mathbf{k}^{T}(\mathbf{x}_{i} - \mathbf{a}_{j}) \approx \|\mathbf{x}_{i} - \mathbf{a}_{j}\| \cos(\alpha_{ij}), \forall i \in \mathcal{T}, \forall j \in \mathcal{E}_{A_{i}}. \quad (9) $$

where $\mathbf{c}_{ij} = [-\sin(\phi_{ij}), \cos(\phi_{ij}), 0]^T$ and $\mathbf{k} = [0,0,1]^T$. To grant more significance to nearby connections, bring into play weights, $\mathbf{w} = [\sqrt{w_{ij}}]$, where $w_{ij} = 1 - \hat{d}_{ij} / \sum_{\forall v \in \mathcal{T}, v \in \mathcal{E}_{A_i} \cup \mathcal{E}_{T_i}} \hat{d}_{ij}$, and substitute $\|\mathbf{x}_i - \hat{\mathbf{a}}_j\|$ in (9) with $\hat{d}_{ij}$ described in (2).

According to the WLS criterion and (7), (8) and (9) every target updates its location through minimization of the below estimator:

$$ \hat{\mathbf{x}}_{i}^{(t+1)} = \arg \min_{\mathbf{x}_{i}} \sum_{j \in \mathcal{E}_{A_{i}} \cup \mathcal{E}_{T_{i}}} w_{ij} (\lambda_{ij} \|\mathbf{x}_{i} - \hat{\mathbf{a}}_{j}\|^2 - d_{0}^{2})^{2} + \sum_{j \in \mathcal{E}_{A_{i}}} w_{ij} (\mathbf{c}_{ij}^{T}(\mathbf{x}_{i} - \mathbf{a}_{j})^{2} + \sum_{j \in \mathcal{E}_{A_{i}}} w_{ij} (\mathbf{k}^{T}(\mathbf{x}_{i} - \mathbf{a}_{j}) - \hat{d}_{ij} \cos(\alpha_{ij}^{A}))^{2}). \quad (10) $$

The above WLS problem is not convex and does not have a closed-form solution. Still, (10) may be expressed as a quadratic programming problem of which the global solution is found effortlessly [7]. Applying the replacement $\mathbf{y}_i = [\mathbf{x}_i^T, \|\mathbf{x}_i\|^2]^T, \forall i \in \mathcal{T}$, equation (10) is modified as:

$$ \mathbf{y}_{i}^{(t+1)} = \arg \min_{\mathbf{y}_{i}} \|\mathbf{W}(\mathbf{A}\mathbf{y}_{i} - \mathbf{b})\|^2 \\ \text{subject to} \\ \mathbf{y}_{i}^{T}\mathbf{D}\mathbf{y}_{i} + 2l^{T}\mathbf{y}_{i} = 0, \quad (11) $$

where $\mathbf{W} = I_3 \otimes \text{diag}(\mathbf{w})$, and $\otimes$ denotes the Kronecker product, $I_q$ represents the $q$ dimensional identity matrix, and
$$A = \begin{bmatrix} \vdots & \vdots \\ -2\lambda_{ij}\hat{a}_j^T & \lambda_{ij} \\ \vdots & \vdots \\ c_{ij}^T & 0 \\ \vdots & \vdots \\ k^T & 0 \\ \vdots & \vdots \end{bmatrix}, b = \begin{bmatrix} d_0^2 - \lambda_{ij} \|\hat{a}_j\|^2 \\ \vdots \\ c_{ij}^T a_j \\ \vdots \\ k^T a_j + \hat{d}_{ij} \cos(\alpha_{ij}) \end{bmatrix}, D = \begin{bmatrix} I_3 & \mathbf{0}_{3 \times 1} \\ \mathbf{0}_{1 \times 3} & 0 \end{bmatrix}, l = \begin{bmatrix} \mathbf{0}_{3 \times 1} \\ \frac{1}{2} \end{bmatrix},$$

i.e., $A \in \mathbb{R}^{3|\varepsilon_{\mathcal{A}i}|+|\varepsilon_{\mathcal{T}i}| \times 4}$, $b \in \mathbb{R}^{3|\varepsilon_{\mathcal{A}i}|+|\varepsilon_{\mathcal{T}i}| \times 1}$, $W \in \mathbb{R}^{3|\varepsilon_{\mathcal{A}i}|+|\varepsilon_{\mathcal{T}i}| \times 3|\varepsilon_{\mathcal{A}i}|+|\varepsilon_{\mathcal{T}i}|}$, and $\mathbf{0}_{p \times q}$ is the $p \times q$-dimensional matrix of all zeros.

Both the objective and the restraint functions in (11) are quadratic. This type of problem is identified as GTRS [7]; it may be solved *exactly* by a bisection method [7].

Assuming that $\mathcal{C}$ is the set of colors of the nodes (in order to coordinate the network [12]), **Algorithm 1** encapsulates the suggested distributed SR-WLS algorithm. Lines 5-7 are carried out concurrently by all targets $i \in \mathcal{C}_c$, which might reduce the realization time of the approach. Information interchange occurs exclusively at Line 7, when targets send their location updates $\hat{x}_i^{(t+1)}$ to their neighbors. Because $\hat{x}_i^{(t+1)} \in \mathbb{R}^3$, one concludes that the new approach obliges a transmission of a maximum $3 \times T_{\max} \times M$ real values. The worst case computational complexity of the proposed method is linear, i.e., $T_{\max} \times M \times O(N_{\max} \times \max_i\{3|\varepsilon_{\mathcal{A}i}| + |\varepsilon_{\mathcal{T}i}|\})$, where $N_{\max}$ denotes the maximum allowed number of iterations of the bisection method. In the further text, we will refer to **Algorithm 1** as "SR-WLS".

<caption>**Algorithm 1** - The proposed distributed SR-WLS algorithm</caption>

Require: $\hat{X}^{(0)}, T_{\max}, \mathcal{C}, \mathbf{a}_j, \forall j \in \mathcal{A}$

1. Set: $t \leftarrow 0$

2. replicate

3. for $c = 1, ..., \mathcal{C}$ do

4. for all $i \in \mathcal{C}_c$ (simultaneously) do

5. Gather $\hat{a}_j, \forall j \in \varepsilon_{\mathcal{A}i} \cup \varepsilon_{\mathcal{T}i}$

6. $\hat{x}_i^{(t+1)} \leftarrow \text{solve (11)}$

7. Broadcast $\hat{x}_i^{(t+1)}$ to $\hat{a}_j, \forall j \in \varepsilon_{\mathcal{A}i} \cup \varepsilon_{\mathcal{T}i}$

8. end for

9. end for

10. $t \leftarrow t + 1$

11. until $t < T_{\max}$

## 5 Performance Evaluation

A set of simulation results is presented here with the intention of assessment of the performance of the new approach in the view of the estimation precision and convergence. All of the considered approaches were solved via MATLAB. In order to demonstrate the benefit of fusing two radio measurements versus traditional
localization systems, we include also the performance results of the proposed method when only RSS measurements are employed, called here “SR-WLS<sub>RSS</sub>”.

In order to produce the radio measures, models (1), (3) and (4) are employed. A random arrangement of sensors within a box with the edge span $B = 20$ m in each Monte Carlo ($M_c$) trial was considered. If not declared otherwise, the PLE was set to $\gamma = 3$, the reference distance to $d_0 = 1$ m, the reference path loss to $L_0 = 40$ dB, and the maximum number of steps in the bisection procedure was $N_{max} = 30$. In practice however, it is almost impossible to perfectly estimate the value of the PLE. Thus, in order to account for a realistic measurement model mismatch and test the robustness of the considered approaches to flawed knowledge of the PLE, the true PLE for each connection was chosen from a uniform distribution on an interval $\gamma_{ij} \in [2.7, 3.3], \forall (i, j) \in E_{cA} \cup E_{cT}$. Finally, the zero estimate of the targets' locations, $\hat{\mathbf{X}}^{(0)}$, was assumed to be at the intersection of the big diagonals of the cube area. The main performance criterion used is the normalized root mean square error (NRMSE), defined as

$$ \text{NRMSE} = \sqrt{\sum_{i=1}^{M_c} \sum_{j=1}^{M} \frac{\|\mathbf{x}_{ij} - \hat{\mathbf{x}}_{ij}\|}{M M_c}}, $$

where $\hat{\mathbf{x}}_{ij}$ represents the estimation of the real location of the j-th target, $\mathbf{x}_{ij}$, in the i-th $M_c$ trial.

<figure><img src="image_3.png" /><figcaption>Fig. 2. NRMSE set against t assessment, when $M = 50$, $R = 6.5$ m, $\sigma_{n_{ij}} = 3$ dB, $\sigma_{m_{ij}} = 6$ deg, $\sigma_{v_{ij}} = 6$ deg, $\gamma_{ij} \in [2.7, 3.3]$, $\gamma = 3$, $B = 20$ m, $L_0 = 40$ dB, $d_0 = 1$ m, $M_c = 500$.</figcaption></figure>

Figs. 2 and 3 illustrate the NRMSE versus t performance when $R = 6.5$ m, $\sigma_{n_{ij}} = 3$ dB, $\sigma_{m_{ij}} = 6$ deg, and $\sigma_{v_{ij}} = 6$ deg, for $M = 50$, $N = 20$ and $N = 30$, and $N = 20$, $M = 50$ and $M = 60$, respectively. From these figures, we can notice that the output from considered approaches improves as t grows, as anticipated. Also, it can be seen from Fig. 2 that the performance of all algorithms improves considerably as more
anchors are included into the network. This behavior is anticipated, because when *N* is increased more reliable information is introduced in the network. Furthermore, **Fig. 3** reveals that both methods need a somewhat elevated number of repetitions to converge for augmented *M*. Nonetheless, the estimation precision of the considered algorithms improves when more targets are added in the network. Moreover, from **Figs. 2** and **3** one can perceive that all major changes in the performance for the considered algorithms occur in the first few iterations ($t \le 15$), and that the performance gain is insignificant afterwards. This result is very important because it shows that our approach necessitates a low number of signal broadcast, which might augment the exploitation productivity of the spectrum, a valuable resource for wireless communications. It also shows that our algorithm is energy efficient; the communication stage is considerably more demanding (in terms of energy) than the data processing one [1],[12].

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** NRMSE versus *t* comparison, when *N* = 20, *R* = 6.5 m, $\sigma_{n_{ij}}$ = 3 dB, $\sigma_{m_{ij}}$ = 6 deg, $\sigma_{v_{ij}}$ = 6 deg, $\gamma_{ij} \in [2.7, 3.3]$, $\gamma = 3$, B = 20 m, $L_0$ = 40 dB, $d_0$ = 1 m, $M_c$ = 500.</figcaption></figure>

**Fig. 4** illustrates the NRMSE versus $\sigma_{n_{ij}}$ (dB) assessment, when $N = 20, M = 50, R = 6.5$ m, and $T_{max} = 30$. In this figure, one can perceive the performance degradation of the proposed algorithm as the quality of the RSS measurement reduces, as projected. Nonetheless, we can see from the figure that the deterioration in the performance is lower than 15 % for the proposed algorithm, which is relatively low for the considered noise range. Finally, **Fig. 4** confirms the effectiveness of measurement fusion, showing a gain of almost 1 m when the hybrid system is employed in comparison with the traditional RSS system.
## 6 Conclusion

Here, the hybrid RSS/AoA target localization problem in a collaborative 3-D WSN was addressed. By fusing information from RSS and AoA, we devised a new problem formulation based on a nonconvex WLS. We showed that the derived novel estimator is suitable for distributed implementation, and we presented a novel non-centralized procedure founded on the GTRS framework. This new algorithm is light in terms of computational complexity, and it can provide accurate localization in a variety of scenarios, having exhibited excellent performance both in view of the estimation precision and convergence. Moreover, the simulation results confirmed the effectiveness of combining RSS and AoA radio measurements in comparison with conventional RSS localization system, displaying a notable enhancement of the estimation precision.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** NRMSE versus $\sigma_{n_{ij}}$ (dB) assessment, for $N = 20$, $M = 50$, $R = 6.5$ m, $\sigma_{m_{ij}} = 1$ deg, $\sigma_{v_{ij}} = 1$ deg, $\gamma_{ij} \in [2.7, 3.3]$, $\gamma = 3$, $T_{max} = 30$, $B = 20$ m, $L_0 = 40$ dB, $d_0 = 1$ m, $M_c = 500$.</figcaption></figure>

**Acknowledgments.** This work was partially supported by Fundação para a Ciência e a Tecnologia under Projects PEst-OE/EEI/UI0066/2014, and PEst-OE/EEI/LA0008/ 2013 (IT pluriannual founding and HETNET), PEst-OE/EEI/UI0066/2011 (UNINOVA pluriannual founding), ADIN PTDC/EEI-TEL/2990/2012, COPWIN PTDC/EEI-TEL/1417/2012 and PTDC/EEITEL/6308/2014-HAMLeT, as well as the grants SFRH/BPD/108232/2015, SFRH/BD/91126/2012 and Ciência 2008 Post-Doctoral Research grant.
References

1. N. Patwari. *Location Estimation in Sensor Networks*. Ph.D. Thesis, University of Michigan, Ann Arbor, MI, USA, 2005.

2. K. Yu, "3-D Localization Error Analysis in Wireless Networks," *IEEE Trans. Wirel. Commun.*, vol. 6, no. 10, pp. 3473--3481, Oct. 2007.

3. R. Fletcher. *Practical Methods of Optimization*. John Wiley & Sons, Chichester, UK, 1987.

4. S. Wang, B. R. Jackson, and R. Inkol, "Hybrid RSS/AOA Emitter Location Estimation Based on Least Squares and Maximum Likelihood Criteria," *IEEE QBSC*, pp. 24--29, Jun. 2012.

5. L. Gazzah, L. Najjar, and H. Besbes, "Selective Hybrid RSS/AOA Weighting Algorithm for NLOS Intra Cell Localization," *IEEE WCNC*, pp. 2546--2551, Apr. 2014.

6. Y. T. Chan, F. Chan, W. Read, B. R. Jackson, and B. H. Lee, "Hybrid Localization of an Emitter by Combining Angle-of-Arrival and Received Signal Strength Measurements," *IEEE CCECE*, pp. 1--5, May. 2014.

7. A. Beck, P. Stoica, and J. Li, "Exact and Approximate Solutions of Source Localization Problems," *IEEE Trans. Sig. Process.*, vol. 56, no. 5, pp. 1770--1778, May 2008.

8. A. Koubaa and M. B. Jamaa, "Taxonomy of Fundamental Concepts of Localization in Cyber-Physical and Sensor Networks," *Technical Report*, CISTER, Research Center in Real-Time & Embedded Computing Systems, Feb. 2013.

9. T. S. Rappaport. *Wireless Communications: Principles and Practice*. Prentice-Hall: Upper Saddle River, NJ, USA, 1996.

10. M. L. Sichitiu and V. Ramadurai, "Localization of Wireless Sensor Networks with a Mobile Beacon," *IEEE MASS*, pp. 174--183, Oct. 2004.

11. S. M. Kay. *Fundamentals of Statistical Signal Processing: Estimation Theory*. Prentice-Hall: Upper Saddle River, NJ, USA, 1993.

12. S. Tomic, M. Beko, and R. Dinis, "Distributed RSS-Based Localization in Wireless Sensor Networks Based on Second-Order Cone Programming," *Sensors*, vol. 14, no. 10, pp. 18410-18432, Oct. 2014.