An Efficient Hybrid RSS-AoA Localization for 3D Wireless Sensor Networks

Abstract

Wireless sensor networks (WSNs) enable many applications such as intelligent control, prediction, tracking, and other communication network services, which are integrated into many technologies of the Internet-of-Things. The conventional localization frameworks may not function well in practical environments since they were designed either for two-dimensional space only, or have high computational costs, or are sensitive to measurement errors. In order to build an accurate and efficient localization scheme, we consider in this paper a hybrid received signal strength and angle-of-arrival localization in three-dimensional WSNs, where sensors are randomly deployed with the transmit power and the path loss exponent unknown. Moreover, in order to avoid the difficulty of solving the conventional maximum-likelihood estimator due to its non-convex and highly complex natures, we derive a weighted least squares estimate to estimate jointly the location of the unknown node and the two aforementioned channel components through some suitable approximations. Simulation results confirm the effectiveness of the proposed method.

1. Introduction

Localization plays a beneficial role in enabling various applications in wireless sensor networks (WSNs), especially for location-aware applications. By computing the locations of wireless nodes in the networks, it is possible to reveal the knowledge for the services that the users have used. Regarding the WSN characteristics (e.g., short-range communication, limited battery and capacity, etc.), one big challenge is to solve the localization problem with the constraints in energy consumption and communication bandwidth [1,2]. For example, developing a smart parking system for vehicles produces a 2D/3D plan to find the nearest parking slot. In order to delivery a better parking experience, other Internet-of-Things (IoT) functionalities such as mobile location service should be integrated to achieve a real-time availability status. Detailed surveys on location-based applications can be found in [3,4,5,6].

Recently, localization based on hybrid measurements has received much attention due to its accuracy and wide range of applications in WSNs such as surveillance, navigation, and tracking. A range-based hybrid scheme can be any combination of time-of-arrival (ToA), time-difference-of-arrival (TDoA), angle-of-arrival (AoA), and received signal strength (RSS) [7,8,9]. Among them, the hybrid RSS-AoA localization can significantly improve the estimation accuracy and have lower costs [10] because the AoA measurement can be adopted to achieve self-localization, while RSS measurement is obtained easily with low accuracy in harsh environments. For instance, it has been shown in [6] that a minimum of two AoA measurements from two different anchors is enough to obtain a location estimation of a target. Moreover, the hybrid RSS-AoA localization can avoid the strict requirement of precise timing synchronization in ToA and TDoA measurements. Thus, the localization problem with hybrid RSS-AoA measurements needs to be carefully investigated.

Despite the recent interest in the hybrid RSS-AoA approach, there remain two key challenges for practical implementation. First, it is assumed that the receiver is able to observe and extract hybrid measurement from the same received signal. However, this assumption may not hold in practice as practical sensors are not yet able to measure both measurements directly. Second, the conventional receiver architecture designed for RSS and AoA may operate with different power sensitivities at the receiver. Thus, for a system that involves hybrid measurements, the receiver architecture should be optimized according to a specific application purpose. Regarding the RSS, all sensor nodes are supposed to have radios. The distance information thus can be obtained from the RSS measurements. For the AoA, an antenna array [11,12,13,14,15] or an omnidirectional antenna [16,17] is implemented at the anchor nodes to measure AoA measurements from different normal sensors. In such settings, all normal sensors are equipped with a digital compass [18,19] to obtain the orientation information and to reduce energy consumption. When both the distance and the AoA measurements are available at each anchor node, the target location can be uniquely determined by using a certain number of anchors [11,12,13]. In order to quickly access range information, a quantized RSS method is implemented to estimate the distance. In summary, the main contributions of this paper are summarized as follows.

(i)

We introduce a 3D hybrid localization that combines the RSS and the AoA measurements, where the channel parameters, i.e., the transmit power and the path loss exponent (PLE) are unknown.

(ii)

We propose a fast and accurate localization scheme where an unbiased solution is developed based on the hybrid measurements. In particular, we transform the original maximum likelihood estimator (MLE), which is non-convex and difficult to solve, into an alternative form after some suitable approximations, then apply a certain numerical optimization.

(iii)

We derive the Cramer–Rao lower bound (CRLB) and geometric dilution of precision (GDoP) analysis for the proposed weighted least squares (WLS).

(iv)

In the simulation results, we compare our work with other existing localization approaches under some likely practical conditions in terms of root mean squared error and bias.

The remainder of the paper is organized as follows. In Section 2, we introduce a 3D localization problem including the system description and related works. In Section 3, we present step-by-step procedures on our solution for improving the localization accuracy and reducing the run time. Simulation results are demonstrated in Section 4, followed by conclusions and future works in Section 5.

2. System Description and Related Works 2.1. Notations

The following notations are used throughout the paper.E(·)andVar(·)denote the expectation and the variance operators, respectively. We denoteA=[_aij]as a matrix whose entries are_aij,_Inas then×nidentity matrix and0as all-zero vectors or matrices.^(·)Tand^(·)−1correspondingly denote the transpose and the inverse operations.diag(·)denotes the diagonal matrix.N(μ,^σ2)denotes the Gaussian distribution with a meanμand a variance^σ2.||·_||1and||·||denote the_ℓ1and the_ℓ2-norm, respectively.Unif[a,b]denotes the uniform distribution on the interval[a,b]._[y][i:j]denotes theithto thejthelements of a vectory.

2.2. Problem Formulation

We assume there are N anchor nodes and one unknown node in a 3D WSN, where the locations of the anchor nodes are denoted by_a1,⋯,_aN∈^R3(_ai=^{[_xi,_yi,_zi]T},i=1,⋯,N)and the location of the unknown node is denoted byx=^[x,y,z]T , as illustrated in Figure 1. We also assume the anchors are equipped with hybrid RSS-AoA antennas. Thus, the combined measurements are given by:

_Pi=_P0−10β_{log10di d0}+_mi,_ϕi=^tan−1y−_yix−_xi+_ni,_αi=^cos−1z−_zidi+_vi,

where_Piis the received signal power (in dB) of the signal transmitted by the target node at theithanchor,_ϕiand_αiare the azimuth and the elevation angles, respectively,_di=||x−_ai||is the distance between theithanchor node and the target node,_P0is the power measured at a reference distance_d0,βis the PLE, and_mi∼N(0,_σmi2),_ni∼N(0,_σni2),_vi∼N(0,_σvi2)represents the log-normal shadowing and the measurement errors for the azimuth and the elevation angles, respectively. Without loss of generality, we assume_d0=1m. Note that_P0is relevant to the transmit power_PT. Our goal is to determinex based on the measurements (1), where_P0andβare treated as unknown parameters.

Given the observation vectorr=^{[PT,ϕT,αT]T}∈^R3N, whereP=[_Pi],ϕ=[_ϕi]andα=[_αi](i=1,⋯,N), the joint probability density function (PDF) is given as:

p(r|θ)=∏i=13N12π_σi2exp−12_σi2^{[_ri−_fi(θ)]2},

whereθ=^{[x,y,z,_P0,β]T},_fi(θ)=_P0−10β_{log10di d0},_fi+N(θ)=^tan−1y−_yix−_xi,_fi+2N(θ)=^cos−1z−_zidi,_σi2=_σmi2,_σi+N2=_σni2,_σi+2N2=_σvi2(i=1,⋯,N). The MLE to maximize the prior PDF is obtained by solving the following non-convex optimization problem.

θ^=argminθ∑i=13N1_σi2^{[_ri−_fi(θ)]2}.

This estimator is evidently non-convex and difficult to solve, especially when both_P0andβare unknown. In the next sections, we will derive a suboptimal technique to solve this problem.

2.3. Strategy for Parameter Estimation

In this section, we consider a general case where both transmit power_P0and PLEβ are unknown, then transform (3) to an alternative form that is easier to solve. For instance, when we need to estimate the target location of a mobile nodein an asynchronous system, it is often desirable to develop a distributed algorithm that is independent of any timing parameters. Furthermore, the system may need to be frequently recalibrated for optimal measuring. Thus, the joint estimate of the channel parameters (e.g., PLE,_P0), as well as the target location can be computed only from the received signals. We observe that ifx is known, (3) becomes linear with respect to_P0andβ. Thus, if we obtain an estimate for the target locationx, we can also obtain the estimates for_P0andβ. For a given maximum number of iterations M, the estimation procedures are described as follows.

1. Step 1: Set up initial guessing intervals for both_P0andβ, and pick_P^0∈[_P0min,_P0max]andβ^∈[_βmin,_βmax].

2. Step 2: Solve (3) using(_P^0,β^)to obtain an estimate forx^.

3. Step 3: Use^{x^(k−1)}to obtain the linear solutions for(_P^0(k),^{β^(k)})by the following.

(P^0(k),β^(k))=argmin(P0,β)∑i=1N1σi2Pi−P0+βti(k−1)2,

where_ti(k−1)=10_log10||^{x^(k−1)}−_ai||. If either_P^0(k)or^{β^(k)}does not lie in the predefined interval, i.e.,_P^0(k)∉[_P0min,_P0max]or^{β^(k)}∉[_βmin,_βmax], then stop the iteration. Otherwise, go to the next step.

4. Step 4: Use the update(_P^0(k),^{β^(k)}) to solve (3) again and get a new update forx,^{x^(k)}. Then, check the stopping criteria: stop ifk>Mor|f(^θ(k))−f(^θ(k−1))|/f(^θ(k−1))≤ϵ, otherwise,k=k+1, and go to Step 2.

2.4. Related Works

A typical wireless localization process follows three main step: range identification, location estimation, and location refinement. Due to hardware characteristics, ranging is not always available for wireless devices. Rather, it has been extracted from different characteristics of radio signals. Many localization algorithms utilized hybrid RSS-AoA measurements [14,15,19,20,21,22,23] because they are more accurate than other methods. The tradeoff is that using hybrid system provides more information on sensor location while increasing network complexity and implementation costs. Thus, according to a specific application purpose, we can select an appropriate localization scheme with a given level of complexity and accuracy.

There are some related literature works on hybrid RSS-AoA measurements [14,15,19,22,23,24]. According to the capabilities of diverse hardware, we may classify the conventional hybrid RSS-AoA-based localization into either of two main categories: nonlinear minimization localization or convex approximation localization. In the first category, supposing that the measurements from the reference node are only corrupted by zero-mean Gaussian noise and the reference nodes are identical or much closer to each other than the target, the basic idea is solving (3) with the given information. In particular, this nonlinear minimization problem can be solved by the Newton–Gauss iteration. For instance, ref. [19] proposed the least squares (LS) and optimization estimators using the distance and the AoA information. In [20], hybrid RSS-AoA LS and maximum likelihood (ML) estimators were derived for the emitter geolocation. The method of combining ML-LS implicitly assumes that a rough estimate range between anchor-target can be obtained. The cost function given by (3) weakly depends on this value; thus, the roughness will not significantly affect the solution. It has been reported in [22] that the accuracy of such a scheme suffers from the environmental dynamics. On the other hand, departing from the first category, most of the works in the second category simplify the location region update due to the target range and the network communication cost. The objective function (3) can be solved by itself or approximated by an alternative form. The approximation for (3) often makes use of the convex constraint. One typical example is [22] where the geometric constraints (distance and angle measurements) between target and anchor node are represented in linear inequalities. All of them are combined to form a single semidefinite programming (SDP) problem. For each iteration, the bounding region for each node is updated accordingly. In [14], the authors presented a hybrid weighted LS (WLS) approach to intra-cell target localization in non-line-of-sight environments, while another closed-form WLS was developed in [15]. The common characteristics of those approaches are their concise problem formulation with clear model representation and refined mathematical solutions, while likely precluding themselves in practice due to their complexities. For example, the complexity of solving the SDP in [22] is at leastO(_Ncv3), where_Ncv is the number of convex constraints needed to describe the network connectivity and measurement information. Other approaches such as SDP and second-order cone programming were proposed in [22,23]. However, the aforementioned works were designed either on 2D space [14], or had a high computational complexity [22,23], or only worked well with low noise [14,19,24]. Moreover, the RSS signals mainly depend on the transmit power of the target node_P0(e.g., its battery) and the PLEβ . Only a few studies have considered the transmit power in (1) as an unknown parameter [15,20,22,23], while in practice, these two parameters may vary significantly in different times and places. To the best of our knowledge, there is no existing work that derives a solution for hybrid RSS-AoA in 3D space when both transmit power and PLE are unknown variables. Therefore, in the paper, we derive a suboptimal solution for this case.

Moreover, given the combined RSS-AoA measurements, our localization scheme is proposed based on the following assumptions. The first one is that the AoAs can be measured and calibrated from the omnidirectional antennas. The ambiguities in the hybrid receiver architecture and its sensitivity are neglected in this paper. We also do not address the noise uncertainty issue over fading channel factors through the localization problem formulation.

3. The Proposed Weighted Least Squares Approach 3.1. An Inspiring Example

In this section, we give an example for the system utilizing hybrid measurements as illustrated in Figure 2. In particular, we consider an intelligent vehicle parking system that consists of the following three parts: (i) a vehicle carries a mobile phone, which is able to connect with the Internet; (ii) a sensor system consists of several sensors placed at the parking slots to detect vehicle presence/absence and many anchor nodes and relay nodes located over the parking area. (iii) a base station receives the localization request message from the mobile phone and provides a possible parking slot. The relay node collects the sensor readings and transmits to the base station for further processing. The anchor node is placed at a known location and helps to cover the network connectivity. As shown in Figure 2a, a typical parking site is installed in several layers to optimize the parking space by saving on the depth and the height. Thus, it may have a lower maintenance and operational cost and a lower construction cost compared to the conventional ones. Each layer represents a parking level, which usually consists of three main components: anchor nodes, relay nodes, and sensors, as depicted in Figure 2b. Given this typical setting, a localization scheme is proposed based on two assumptions. The first one is that the layout of the parking system is organized in an optimal model, i.e., it satisfies the coverage, connectivity, and geolocation constraints, known as the optimal sensor placement model. Thus, the system is able to function effectively in terms of efficient data collection and financial cost. Second, the available of RSS holds over the whole network.

In order to apply the cooperative hybrid localization scheme, we performed two procedures: region identification and location refinement. The basic idea is that we first determined which area may have available slots by quantizing the RSS readings. Second, by utilizing the AoA measurements, we estimated the slot location. After that, a refinement method was applied for better accuracy.

3.2. Issues and Motivations

The aforementioned works [14,15,19,20,21,22,23] faced several issues. First, due to the nonlinear nature of the localization problem, the location estimation via the classical method such as the MLE requires high complexity and computation cost. On the other hand, low complexity approaches like the LS or WLS methods are based on the low noise assumption, which is hard to achieve in real-world applications. Thus, constructing a localization scheme that achieves very accurate positioning in a short time for both low and high noise conditions, as well as under different channel settings is needed. Second, it can be seen that the measurement noises have a large impact on the solution accuracy. Consequently, these effects must take part in the localization procedure. Thus, unlike the conventional approaches, we investigated the noise covariance of hybrid measurements and produced an unbiased estimate under a realistic assumption of unknown transmit power and PLE. Third, the current setting for the anchor nodes is seriously restricted to the ideal scenarios, that is the anchors are positioned at the vertices of a cube region in which the target is inside the convex hull formed by them; thus, no GDoP effect can be examined. For a more practical scenario, we considered a random network in the simulation setting, where the target may fall outside the convex hull defined by the anchors.

3.3. Quantized RSS-Based Ranging

We considered a realistic modification of (1), i.e., instead of using the analog version of the RSS measurements, we took the quantized values of them. In this case, we no longer have access to the actual RSS reading, but only partial information, related to the ability of the device to communicate with others. In particular, by using an S-level quantization operatorQ(·), we obtain:

_Qi=Q(_Pi)=0,_L0≤_Pi≤_L11,_L1≤_Pi≤_L2⋮⋮S−1,_LS−1≤_Pi≤_LS,

where_L0,⋯,_LSare quantized levels,_L0=−∞dBm and_LS=∞dBm. Let_d^ibe the distance at which the mean received power is_L¯i, i.e., the estimated distance between the targetxand the anchor_ai; we can compute_d^ifrom_L¯ias:

_d^i=_d0 ^{10_P0−_L¯i10β}.

Figure 3 shows an example of distance estimation in indoor (PLE = 1.5 and additional reflection coefficients) and outdoor (PLE = 3.5) environments. Other parameters are given in Section 4. In each run, the estimated range was obtain from (6) and the quantization process. The coordinate of a single target was randomly generated for the given area. It has been shown in [3] that the relation between the average received power and the distance is proportional. The estimated range was obtained from the raw measurements over 100 realizations. We observed that in the outdoor environment the quantized scheme performed more stably and had smaller errors compared to the indoor space. This is because the wireless signals in the indoor environment suffered from multipath and reflection effects from walls and obstacles. In addition, the results provided a clue that could help in choosing an appropriate quantization level. A good accuracy of about 1 m average error within a range of 10 m can be achieved, which is acceptable for ranging with inexpensive devices.

3.4. Location Estimation

In the presence of both range and angle measurements, we can obtain the target location as:

x=_xi+_d^isin_α^icos_{ϕ^i δi},y=_yi+_d^isin_α^isin_{ϕ^i δi},z=_zi+_d^icos_{α^i δi},

where_δiis the bias-reducing constant for the RSS-AoA measurements given by:

_δi=exp_σni22+_σvi22−_σni22^(βγ)2.

Hence, (7) can be rewritten in a matrix form asAx+q=b^, whereA=diag{_eN,_eN,_eN}∈^R3N×3,b^=^{[_b^x,_b^y,_b^z]T}∈^R3N×1with three elements as:

_b^x=^{[_x1+_d^1sin_α^1cos_{ϕ^1 δ1},⋯,_xN+_d^Nsin_α^Ncos_{ϕ^N δN}]T},_b^y=^{[_y1+_d^1sin_α^1sin_{ϕ^1 δ1},⋯,_yN+_d^Nsin_α^Nsin_{ϕ^N δN}]T},_b^z=^{[_z1+_d^1cos_{α^1 δ1},⋯,_zN+_d^Ncos_{α^N δN}]T},

andqis the noise vector, which has a zero mean vector and covariance_Σx=E[(b^−b)^{(b^−b)T}]. Thus, the WLS solution is obtained by minimizing the following cost function.

ϵ(x)=^{(b^−Ax)T} _Σx−1(b−Ax),

which yields the intermediate solution:

x^=^{A† b†},

where^A†=^{[AT _Σx^−1A]−1 AT}and^b†=_Σx^−1b^.

3.5. Linear Approximation and Channel Factor Refinement

From Section 3.4, the target node should obtain its initial location estimate. In some cases, it might not be well localized since it does not have enough anchor nodes. In order to have knowledge on the location estimate and internal relay parameters of its neighbors, the location estimate of the target is refined by applying a simple LS algorithm. Thus, the square root term in (6) can be linearized by using the Taylor series approximation. Following Section 2.3, first we pick initial_P^0∈[_P0min,_P0max]andβ^∈[_βmin,_βmax] . The right-handed side of (6) can be approximated using the first-order Taylor expansion as:

_d^i=^(x−_xi)2+^(y−_yi)2+^(z−_zi)2≈_d˜i+∂_d^i∂xδx+∂_d^i∂yδy+∂_d^i∂zδz≈_d˜i+x˜−_xid˜i(x−x˜)+y˜−_yid˜i(y−y˜)+z˜−_zid˜i(z−z˜).

Here,^{(_xi,_yi,_zi)T}is the known coordinate of either an anchor or a known neighbor node and_d˜i=^(x˜−_xi)2+^(y˜−_yi)2+^(x˜−_zi)2. Thus, the location estimate is updated byx^←x^+δx,y^←y^+δyandz^←z^+δz.

Next, the unknown PLE can be expressed asβ^=_β01+β^−_β0β0=_β0[1+δ], where_β0≠0is the tuning parameter, which can be set in the simulation, andδ=_β0−β^_β0. For small|δ|and small_P^0−_Pi10_β0δln10, by using the facts that1/(1+δ)≈1−δand^u−δ≈1−δlnu, let_λi=^{10_P^0−_Pi10_β0}, and we perform:

_d^i≈^{10_P^0−_Pi10_β0(1+δ)}η1+ln1010β^_mi≈_λiη(1−δln_λi)1+ln1010β^_mi,

whereη≜^{10(_P0−_P^0)/(10β^)}. In order to take advantage of neighboring anchor nodes, each weight factor is defined as_wi=1−_Pi/_{∑i=1N Pi}.

In the next step, we try to extract the channel factors from the raw sensor data based on (4). By denoting_d¯i=||x¯−_ai||(i.e., the new estimated distance between the target node and theithanchor node), we first rewrite the path-loss model as:

_Pi=_P0−10β_{log10 d¯i}+_mi.

Here, an LS method is applied to compute the estimates of the transmit power_P0and the PLEβ. The solution is obtained by:

[_P¯0,β¯]=^{(_{A1T A1})−1} _{A1T b1},

where_A1=1−10_{log10 d¯1}⋮⋮1−10_{log10 d¯N},_b1=^{[_P1,⋯,_PN]T}. Then, we exploit[_P¯0,β¯] to calculate the distance (12) and use its estimated value to solve the WLS problem given by (10). Note that in Step 1 given in Section 2.3, the initialization state sets up the guessing interval of channel factors based on the mean RSS data between a reference node and anchors in the environment.

Remark 1

(Complexity analysis). Given a maximum number of iterations M, the complexity of solving MLE, LS, WLS, SOCP, mixed SDP/SOCP, and the proposed method is summarized in Table 1.

3.6. Cramer–Rao Lower Bound for the Proposed Hybrid RSS-AoA Localization

In order to provide a limit of estimator performance, we computed the CRLB to reflect the effect of model parameters on the estimation accuracy. Letθ^be an unbiased estimator of a parameterθ; the CRLB for theithelement ofθ^is defined by:

Var(_θ^i)=_{[^J−1(θ)]i,i}≜CRLB(_θ^i).

Letp(r|θ)denote the conditional PDF of the observation datar=^{[PT,ϕT,αT]T}parameterized byθ. The entries of the corresponding Fisher information matrixJ(θ)are derived as:

_[J(x)]i,j=−E∂lnp(r|x)∂_xi∂lnp(r|x)∂_xj.

When applying the quantization operation (5), the conditional PDF of_ri=_Pi∈[_Ls,_Ls−1](i=1,⋯,N) is given by:

p(_ri|θ)=_{∫LsLs−1 pi}(_Li)dP=Φ(_Dis−1−_Dis),

whereΦ(t)=_∫−∞t12π^e−v2/2dv,Φ(_Di0)=1,Φ(_DiS)=0, and_Dis=β_σmi2ln(_di/_d^i)ln(_di/_d0),i=1,⋯,N.

Let_d2,i=^(x−_xi)2+^(y−_yi)2, then we have by some manipulations,

∂_fi(x)∂x=^{10_τi(x−_xi)_di2ln10,−(y−_yi)_d2,i2,(x−_xi)(z−_zi)_{di2 d2,i}T},

∂_fi(x)∂y=^{10_τi(y−_yi)_di2ln10,(x−_xi)_d2,i2,(y−_yi)(z−_zi)_{di2 d2,i}T},

∂_fi(x)∂z=^{10_τi(z−_zi)_di2ln10,0,−_{d2,i di2}T},

∂_fi(x)∂_P0=^1,0,0T,∂_fi(x)∂β=^{−10_{log10 di},0,0T},

where:

_τi=12π∑i=1S^{e−12(_Dis−1)2}β_σmi −_Dis−^{e−12(_Dis)2}β_σmi −_DisΦ(_Dis−1)−Φ(_Dis).

By using the Taylor series, we construct the relationship betweenxandx^as:

x≈x^+∑i=13N∂x∂_fi[_ri−_fi(x)],

where∂x∂u=^{(ATWA)−1 AT}W∂b∂u. Assuming that the measurements are independent and differentiable, the approximated variances of the estimator errors are given by:

Var(_[x^]j)≈∑i=1N^_∂x∂Pij2 _σmi2+^{_∂x∂ϕij2} _σni2+^{_∂x∂αij2} _σvi2,j=1,2,3.

4. Simulation Results

Computer simulations were conducted to show the performance of the proposed approach compared to the MLE (solving by the Gauss–Newton method [25]), the SOCP (using the CVX package [26]), and the WLS method [21]. Simulation parameters have been adopted, and we kept the same set to achieve fairness of comparison among the aforementioned works. Furthermore, the true values of the unknown parameters were used as the starting points. In particular, the following parameters were used:_d0=1m,_P0∼Unif[−15,−5]dBm,β∼Unif[2,5], the standard deviationof the RSS was_σmi =_σ1, and that of the AoA was_σni =_σvi =_σ2 . We used N anchors to locate one unknown sensor node. Those anchor nodes were picked randomly in a cube-shaped region, which had an edge length of 100 m. The unknown node was deployed randomly in this region. In terms of sensor deployment settings, the distance between the normal nodes and the anchor nodes should meet the connectivity condition of the networks. The maximum communication among nodes was set at 25 m. We used (1) to generate hybrid RSS-AoA measurements.

In order to evaluate the performance, the root mean squared error (RMSE) and the bias were considered as follows.

RMSE(x^)=1_Mc∑m=1_Mc||^x(m)−^{x^(m) ||2},

Bias(x^)=1_Mc∑m=1_Mc||^x(m)−^{x^(m)} _||1,

where_Mc=5000is the number of Monte Carlo runs. In fact, small RMSE and bias values can produce the estimated solutions that are close to the traceable actual values. Given a localization algorithm, the RMSE and the bias show how well the estimated target location matches the actual one. More specifically, besides the computation cost, a good localization algorithm must ensure a low RMSE and a stable bias under the same settings scenario. In this section, we will evaluate and compare those aforementioned algorithms under different noise conditions to make accurate validation.

In order to first find an appropriate threshold for the RSS measurements, we investigated the effect of the RSS quantization level on the RMSE performance. Without loss of generality, we set_σ1=_σ22andN=10 . One of the basic choices in selecting the quantization level is the number of anchors. As Figure 4 illustrated, the higher quantized level S resulted in a lower RMSE value. However, the RMSE performance was not improved by increasing the quantization level S. This is because when S becomes large, there is not much difference among the quantized intervals, and boosting error can increase with S, especially when the measurements are noisy. We also observed that the value of quantized level had a significant impact on the performance of the estimated range (6). A small value S reduced the localization expense, but degraded the localization accuracy. When S approached infinity, the estimator in (10) was actually reduced to the conventional weight least squares method with the bias-reducing constant (8). Thus, with a refinement scheme to exploit the above tradeoff, the proposed localization scheme worked well under the condition settings. With an acceptable RMSE, we choseS=8for the next simulation results.

In Figure 5, we evaluate the RMSE with different numbers of anchor nodes N where_σ12=_σ22=0(i.e., noiseless measurements). When the number of anchor nodes N increased, the RMSEs of all estimators decreased significantly and converged to the CRLB. With a fixed number of anchor nodes, the proposed estimator achieved an acceptable RMSE and smaller bias compared to the others.

Figure 6 illustrates the RMSE and the bias versus_σ1when_σ2=1deg andN=13. We observed that the proposed scheme achieved a comparable RMSE and smaller bias when_σ1was small compared to other methods. When_σ1 became large, even our solution had a lower RMSE, but higher bias. This phenomenon may occur from the approximations of the linearization process explained in Section 3.5.

Figure 7 shows the RMSE and the bias versus_σ2when_σ1=3dBm. From this figure, we observe that the proposed method still had a good performance when the noise became large. The bias performance was comparable to the MLE and SOCP and much lower than the WLS. However, it was also acceptable since the proposed method had a lower complexity than the MLE and SOCP estimator.

Remark 2.

In some approaches, the accuracy may be defined as the expected distance between the actual positionxand the estimated onex^, i.e.,E(||x^−x||) . Thus, the evaluation metric can be indicated as the percentage of the results satisfying a predefined accuracy requirement. In this section, we only considered two metrics in (25) and (26) to validate the localization effectiveness. Based on the observation from the simulation results, there were some major factors that caused large errors during the localization process. Assuming that information on the anchor positions was obtained accurately and reliably, the error source may be caused by the distance estimation error from (6), the multipath fading, and the AoA measurement noises. The same phenomenon occurred in other localization schemes.

Remark 3.

It has been seen that that localization system runs based on an event-driven mechanism, i.e., a sensor only collects the RSS-AoA measurements when it receives a localization request. The quantized RSS helps the sensors to know quickly their ranges to the neighbor anchors/relay nodes and to indicate the AoA of the incoming signals. Until the localization packet (sensor ID, RSS-AoA measurements) reaches the neighbor anchor nodes, they collect and deliver back to the BS their locations. Finally, the localization algorithm can be executed by the BS to obtain the sensor location on the received packets. The unknown node only sends one message for localization instead of continuously reporting the measurements to the neighbor anchors; thus, the number of messages required for localization per unknown node can be reduced.

In order to measure the uncertainty of the location estimatex^relative to its meanE(x^) , we considered the GDoP that measured the ratio of the RMSE of a location fix to the statistical RMSE of the ranging measurements as [3]:

GDoP(x^)=1_σrE^{x^−E(x^)T}x^−E(x^).

We now investigate the effects of the GDoP in the network localization. For a given boundξ, we calculated the probability of GDoP smaller thanξ, which means the probability of acceptable estimationP{x^:GDoP(x^)≤ξ} . The exponential increase in the detection probability as we took more measurements from the anchor nodes is illustrated in Figure 8, which plots the performance predicted by (27) for a range ofξwith the number of anchor nodes N. The results are consistent and should approach one exponentially fast as N andξare increased.

5. Conclusions

In this paper, we introduced the localization problem on RSS-AoA measurements in 3D WSNs, then derived a suboptimal solution that can obtain a good accuracy under some mild conditions. Simulation results showed that the proposed method outperformed the MLE, SOCP, and WLS solutions in terms of the RMSE and the complexity and had an acceptable bias performance.

Some effectiveness of the proposed scheme is addressed as follows. First, as we presented in Section 2.4 and Section 3.2, we considered two main problems that had not been solved yet: (i) how to handle the nonconvexity of network deployment with low computational complexity; (ii) how to extract quickly the location information to keep lower maintenance. It is not easy to solve both problems in 2D space, but much more difficult in 3D space. Moreover, since the conventional approaches were often conducted in a distributed or an iterative manner, error during estimation process leads to inaccurate location estimates. Through mathematical analysis and simulation, the proposed approach has demonstrated that it is an effective technique to tackle those issues. Second, with the event-based localization protocol, a sensor node can save its energy and does not require to wake up all the time. Thus, the proposed scheme is suitable for location-aware applications.

In the future work, we are interested in device-free localization algorithms where the “to-be-located target” cannot join the online localization process. It needs to be detected and localized by using the changes of the environment (e.g., RSS changes) collected at scattered devices. Another direction of future research with good potential is investigating different radio propagation under multipath, shadowing affects caused by reflection and scattering from obstacles, and interference with other devices that operate at the same frequency, etc. The location of unknown nodes cannot be uniquely determined due to ranging errors from location ambiguities in a practical scenario. Moreover, the newly-estimated locations may become the references to help with locating other sensor nodes in the localization process for the multiple target localization problem.

[Image omitted. See PDF.]

[Image omitted. See PDF.]; (b) Infrastructure at a single layer."]

[Image omitted. See PDF.]

[Image omitted. See PDF.]

[Image omitted. See PDF.]

[Image omitted. See PDF.]

[Image omitted. See PDF.]

[Image omitted. See PDF.]

Method	Complexity
MLE [20]	O(M∗N3)
ine LS [19]	O(KmaxN), whereKmax is the maximum number of steps in the bisection procedure
ine WLS [14]	O(MN+KN) , where K is the maximum number of steps in the bisection procedure [25]
ine SOCP [21]	O(MSOCP∗(N+4)3.5), whereMSOCPis related to a search interval[p,q]and the solution precisionϵ:MSOCPis the minimum integer that satisfiesM>log2((q−p)/ϵ) .
ine Mixed SDP/SOCP [23]	At leastO(3(81(N+0.5)2))
ine Proposed method	O(MN+KmaxN), whereKmax is the maximum number of iterations for solving (14).

Author Contributions

The work was realized with the collaboration of all the authors. T.L.N.N. and T.D.V. contributed to the main results and code implementation. Y.S. provided the funding, supervised the research, and reviewed the draft of the paper. All authors discussed the results, approved the final version, and wrote the paper.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (2016R1A2B2014497).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

3D three-dimensional

AoA angle-of-arrival

CRLB Cramer–Rao lower bound

GDoP geometric dilution of precision

IoT Internet-of-things

LS least squares

ML maximum likelihood

MLE maximum likelihood estimator

PDF probability density function

PLE path loss exponent

SDP semi-definite programming

SOCP second-order convex programming

RMSE root mean square error

ToA time-of-arrival

TDoA time-different-of-arrival

WLS weighted least squares

WSN wireless sensor network

References

1. Masazade, E.; Niu, R.; Varshney, P.K.; Keskinoz, M. Energy aware iterative source localization for wireless sensor networks. IEEE Trans. Signal Process. 2010, 58, 4824–4835.

2. Niu, R.; Vempaty, A.; Varshney, P.K. Received-signal-strength-based localization in wireless sensor networks. Proc. IEEE 2018, 106, 166–1182.

3. Dargie, W.; Poellabauer, C. Fundamentals of Wireless Sensor Networks: Theory and Practice; John Wiley & Sons: Hoboken, NJ, USA, 2011; Chapter 2.

4. International Electrotechnical Commission. Internet of Things: Wireless Sensor Networks. Available online: https://www.iec.ch/whitepaper/pdf/iecWP-internetofthings-LR-en.pdf (accessed on 6 May 2019).

5. Atzori, L.; Iera, A.; Morabito, G. The internet of things: A survey. Comput. Netw. 2010, 54, 2787–2805.

6. Nguyen, T.L.N.; Shin, Y. A new approach for positioning based on AOA measurements. In Proceedings of the 2013 International Conference on Computing, Management and Telecommunications (ComManTel), Ho Chi Minh City, Vietnam, 21–24 January 2013; pp. 208–211.

7. Gu¨venc, I.; Chong, C. A survey on TOA based wireless localization and NLOS mitigation techniques. IEEE Commun. Surv. Tutor. 2009, 11, 107–124.

8. Xu, J.; Ma, M.; Law, C.L. AOA cooperative position location. In Proceedings of the IEEE GLOBECOM 2008—2008 IEEE Global Telecommunications Conference, New Orleans, LO, USA, 30 November–4 December 2008; pp. 1–5.

9. Wang, G.; Chen, H.; Li, Y.; Jin, M. On received-signal-strength based localization with unknown transmit power and path loss exponent. IEEE Wirel. Commun. Lett. 2012, 1, 536–539.

10. Dai, W.; Shen, Y.; Win, M.Z. Distributed power allocation for cooperative wireless network localization. IEEE J. Sel. Areas Commun. 2015, 33, 28–40.

11. Dogancay, K. 3D pseudolinear target motion analysis from angle measurements. IEEE Trans. Signal Process. 2015, 63, 1570–1580.

12. Badriasl, L.; Dogancay, K. Three-dimensional target motion analysis using azimuth/elevation angles. IEEE Aerosp. Electron. Syst. Mag. 2014, 50, 3178–3194.

13. Shao, H.J.; Zhang, X.P.; Wang, Z. Efficient closed-form algorithms for AOA based self-localization of sensor nodes using auxiliary variables. IEEE Trans. Signal Process. 2014, 62, 2580–2594.

14. Gazzah, L.; Najjar, L.; Besbes, H. Selective hybrid RSS/AOA weighting algorithm for NLOS intra cell localization. In Proceedings of the 2014 IEEE Wireless Communications and Networking Conference (WCNC), Istanbul, Turkey, 6–9 April 2014; pp. 2546–2551.

15. Chan, Y.T.; Chan, F.; Read, W.; Jackson, B.R.; Lee, B.H. Hybrid localization of an emitter by combining angle-of-arrival and received signal strength measurements. In Proceedings of the 2014 IEEE 27th Canadian Conference on Electrical and Computer Engineering (CCECE), Toronto, ON, Canada, 4–7 May 2014; pp. 1–5.

16. P, R.; Sichitiu, M.L. Angle of arrival localization for wireless sensor networks. In Proceedings of the 2006 3rd Annual IEEE Communications Society on Sensor and Ad Hoc Communications and Networks, Reston, VA, USA, 28 September 2006; pp. 374–382.

17. Nasipuri, A.; Li, K. A directionality based location discovery scheme for wireless sensor networks. In Proceedings of the WSNA ’02 Proceedings of the 1st ACM International Workshop on Wireless Sensor Networks and Applications, Atlanta, GA, USA, 28 September 2002; pp. 105–111.

18. Kim, E.; Yoo, S.; Kim, H.; Kim, H. Enhancing base station connectivity with directional information of movement. ICT Express 2016, 2, 122–125. [Green Version]

19. Yu, K. 3-D localization error analysis in wireless networks. IEEE Trans. Wirel. Commun. 2007, 6, 3473–3481.

20. Wang, S.; Jackson, B.R.; Inkol, R. Hybrid RSS/AOA emitter location estimation based on least squares and maximum likelihood criteria. In Proceedings of the 2012 26th Biennial Symposium on Communications (QBSC), Kingston, ON, Canada, 28–29 May 2012; pp. 24–29.

21. Tomic, S.; Beko, M.; Dinis, R.; Montezuma, P. A closed-form solution for RSS/AoA target localization by spherical coordinates conversion. IEEE Wirel. Commun. Lett. 2016, 5, 680–683.

22. Tomic, S.; Marikj, M.; Beko, M.; Dinis, R.; Orfao, N. Hybrid RSS-AoA technique for 3-D node localization in wireless sensor networks. In Proceedings of the 2015 International Wireless Communications and Mobile Computing Conference (IWCMC), Dubrovnik, Croatia, 24–28 August 2015; pp. 1277–1282.

23. Tomic, S.; Beko, M.; Dinis, R. 3-D target localization in wireless sensor network using RSS and AoA measurement. IEEE Trans. Veh. Technol. 2017, 66, 3197–3210.

24. Multi-Level Solutions for Car Parking. Available online: https://www.nbmcw.com/product/car-parking/29700-multilevel-solutions-for-car-parking.html (accessed on 6 May 2019).

25. Cheney, W.; Kincaid, D. Numerical Mathematics and Computing, 6th ed.; Brooks/Cole: Boston, MA, USA, 2008; Chapter 6.

26. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming. Available online: http://cvxr.com/cvx (accessed on 6 May 2019).