Unbalanced Hybrid AOA/RSSI Localization for Simplified Wireless Sensor Networks

Abstract

Source positioning using hybrid angle-of-arrival (AOA) estimation and received signal strength indicator (RSSI) is attractive because no synchronization is required among unknown nodes and anchors. Conventionally, hybrid AOA/RSSI localization combines the same number of these measurements to estimate the agents’ locations. However, since AOA estimation requires anchors to be equipped with large antenna arrays and complicated signal processing, this conventional combination makes the wireless sensor network (WSN) complicated. This paper proposes an unbalanced integration of the two measurements, called 1AOA/nRSSI, to simplify the WSN. Instead of using many anchors with large antenna arrays, the proposed method only requires one master anchor to provide one AOA estimation, while other anchors are simple single-antenna transceivers. By simply transforming the 1AOA/1RSSI information into two corresponding virtual anchors, the problem of integrating one AOA and N RSSI measurements is solved using the least square and subspace methods. The solutions are then evaluated to characterize the impact of angular and distance measurement errors. Simulation results show that the proposed network achieves the same level of precision as in a fully hybrid nAOA/nRSSI network with a slightly higher number of simple anchors.

1. Introduction

Global positioning system (GPS) has made our daily lives much more convenient due to the tracking function available in mobile phones, cars, watches, etc. With GPS receivers equipped in such devices, their movements in outdoor areas can be accurately located and tracked, providing vast benefits in navigation and management. However, since at least four visible GPS satellites are required to locate a receiver’s position, the performance of GPS-based systems becomes inaccurate in indoor environments, dense urban areas, or deep forests [1]. To improve accuracy, the combination between GPS and opportunistic terrestrial signals, such as Wi-Fi, has been explored in the literature, e.g., [2]. In emergency scenarios, such as fire, post-tsunami, post-earthquake, and military missions however, these signals may be unavailable, sabotaged, or even denied. Hence, a wireless sensor network (WSN) must be deployed to provide the localization services for first responders.

Typically, a WSN includes a number of anchors, whose locations are known, and other unknown nodes, called agents. In a range-based WSN, anchors receive signals from unknown nodes, perform ranging measurements, and send information to a center node which estimates agents’ locations [3]. There are many different techniques to obtain ranging information based on time-of-arrival (TOA), time-different-of-arrival (TDOA), time-of-flight (TOF), angle-of-arrival (AOA), and received signal strength indicator (RSSI). The more accurate the ranging information is, the more precise the locations of agents can be estimated. However, the accuracy of ranging measurements depends on the complexity of the anchors. For example, distance estimation based on TOA, TDOA, and TOF can be very accurate, but it requires strict synchronization among anchors and agents [4]. Therefore, a fundamental problem of localization in WSNs is to precisely locate the agents, while the complexity and power consumption of the anchors are minimized.

Among various approaches in the literature, hybrid AOA/RSSI methods have been proven to be very effective [5]. Since RSSI is always available in receivers, using RSSI for localization is one of the most popular methods [6]. Using path loss models, if we know the transmitted power of an agent, its distances to anchors can be determined from RSSI. These distances can be used in a lateration method to draw multiple circles, then their intersections are used to localize the agent. However, as RSSI is affected by multipath, noise, and shadowing, precision of standalone RSSI positioning is low. Meanwhile, if there exist line-of-sight paths from agents to anchors, AOA estimation can provide the directions toward the unknown nodes. Combination of RSSI and AOA, therefore, can significantly reduce the estimation error [7].

An example of combining RSSI and AOA localization can be found in Guideloc systems proposed in [8] for disaster rescue applications. In this system, unmanned aerial vehicles (UAVs) are used to search the positions of the targets. The RSSI and AOA received by the UAVs are combined to control the trajectories of the UAVs to fly over the targets. Thereby, the targets’ locations can be derived from the GPS coordinates of the UAVs. Since both RSSI and AOA measurements are utilized, the search time can be reduced significantly.

In other applications such as localization and tracking of first responders, RSSI and AOA information is collected by the anchors and processed at a data fusion center. In this case, target positioning using hybrid AOA/RSSI data can be cast as a nonconvex optimization problem [9]. Solutions of this problem can be categorized into optimization-based and least square (LS) estimators [10]. Examples of the former approaches can be found in [10,11] using a maximum likelihood (ML) estimator. While the ML algorithm achieves great accuracy, it suffers from severe computational complexity and the global solution is not guaranteed [6]. To reduce complexity while obtaining a comparable accuracy to ML, semidefinite programming and second-order cone programming relaxation techniques have been proposed in [12,13,14,15]. However, the computational complexity of these algorithms grows significantly with the network size [5]. In comparison, the LS estimators are less accurate, but they are more computationally efficient. Therefore, great efforts have been devoted to improve their accuracy by applying a weighting matrix to the LS solution, so called weighted LS (WLS). In [9,16,17,18], a weighting covariance matrix is derived based on noise statistics. However, since noise statistics are not available, they need to be estimated first. As a result, they are not suitable for real-time applications due to the delay occurring in the noise estimation process. Another way to construct the weighting matrix is to utilize the distances among the anchors and the agents as proposed in [19]. It was proved that both RSSI and AOA short-range measurements are more reliable than the long ones. Hence, emphasis should be given to the shorter links. As shown in [5], this WLS estimator achieves remarkable precision, while its complexity is comparable to LS method.

The above-mentioned hybrid AOA/RSSI localization methods assume that AOA and RSSI measurements are available at all anchors. However, AOA estimation is a challenging task which normally requires an antenna array and complicated signal processing. For example, in [20], two perpendicularly orientated directional antennas are employed and the RSSI from these antennas are fitted into parabolic functions to estimate AOA. Otherwise, a large antenna array with multiple signal classification algorithms must be equipped to obtain AOA information [6]. Obviously, such a WSN requiring AOA estimation at every anchor leads to significant hardware complexity and power consumption. Especially in applications such as emergency responses and disaster rescues, which require flexible and fast deployment of WSNs, anchors can be carried by UAVs [8,21]. Therefore, requiring every anchor to be equipped with an antenna array is not practical.

In this paper, an unbalanced hybrid AOA/RSSI localization method, called 1AOA/nRSSI, is proposed to simplify the complexity of the WSN. In particular, LS algorithms are employed to estimate the unknown nodes based on 1 AOA and n RSSI measurements(n≥3). As a result, the WSN in this case only needs a master anchor which is equipped with an antenna array to perform AOA estimation, while other anchors are simple single-antenna transceivers to acquire RSSIs. Simulations are then conducted to show that, with a slightly higher number of simple anchors, the root mean square error of the proposed method is equivalent or superior to that obtained by a fully equipped WSN performing the 3AOA/3RSSI WLS algorithm. Meanwhile, under a significant shadowing environment, the proposed 1AOA/nRSSI method approaches the precision of a 3AOA/3RSSI network if more simple anchors are deployed. The contributions of this paper are summarized as follows. Firstly, unbalanced hybrid AOA/RSSI localization is proposed to reduce the burden of AOA estimation in all anchors of the WSN. Secondly, by simply transforming one AOA information into two virtual anchors with corresponding ranging measurements, we present two algorithms, namely, LS and subspace estimators. The impact of measurement errors on the performance of these two algorithms is evaluated. Finally, simulation results show that the proposed 1AOA/nRSSI approach is promising in low shadowing environments, for example, from a UAV to the ground agents.

The remainder of this paper is organized as follows. Section 2 introduces system models and some related works. Section 3 presents the proposed 1AOA/nRSSI estimator with two different algorithms. Section 4 provides the evaluation results of the two algorithms under different simulation scenarios. Finally, Section 5 concludes the paper. Note that a list of abbreviations used in this paper is given at the end of the paper.

2. Related Works

Consider a WSN including N anchors whose locations are known at_ai=[_aix,_aiy],i=1,…,N, and M agents whose actual locations are_xm=[_xmx,_xmy],m=1,…,M. In this paper, we only present the results for a two-dimensional map for the ease of illustration. At the anchor_ai, the RSSI_Piand AOA_αimeasured from the m-th agent are expressed as

_Pi=_P0−10γ_log10(∥_xm−_ai∥_r0)+_niαi=^tan−1(_xmy−_aiyxmx−_aix)+_mifori=1,…,N,

where_P0is the received power at the reference distance_r0,γis the path-loss exponent (PLE),∥.∥denotes the Euclidean norm,_niis the log-normal shadowing, and_miis the angle estimation noise. Both_niand_miare modeled as zero-mean Gaussian random variables with standard deviations_σni (dB) and_σmi (degree), respectively—i.e.,_ni∼N(0,_σni2),_mi∼N(0,_σmi2). The problem here is to estimate the location of the agent_xmbased on_Piand_αi. Some typical localization methods are presented as follows.

2.1. RSSI Localization

RSSI is used to determine the distance_rifrom the agent to the i-th anchor. Assuming PLE is known at the anchor, the distance_rican be determined by

_ri=_r0^{10_P0−_Pi+_ni10γ}.

The localization problem becomes solving a set of nonlinear equations of two variables_xmx,_xmyexpressed as

^{(_xmx−_aix)2}+^{(_xmy−_aiy)2}=_ri2,fori=1,…,N.

The agent’s location can be found with as few as two anchors using the bilateration method [6], which plots two circles centered at the two anchors and uses one of their intersections as the agent’s position. However, due to the noise effect, the two circles may not intersect and the location cannot be found. Therefore, normally, at least three anchors are required, so Equation (3) becomes an overdetermined set of nonlinear equations. To solve this problem, it can be firstly transformed into linear equations and then solved by a linear LS solution. As presented in [22], choosing the r-th anchor as a reference anchor, and then subtracting the remaining equations in Equation (3) to the r-th one, we have a set of(N−1)linear equations as

(_aix−_arx)_xmx+(_aiy−_ary)_xmy=12_aix2+_aiy2−_arx2−_ary2−_ri2+_rr2,fori=1,…,N,i≠r.

The agent’s location, denoted for this RSSI localization method as_xmR , can be found by solving Equation (4) based on the LS solution as

_x^mR=^{(HTH)−1 HT}b,

where

H=_a1x−_arxa1y−_ary⋮⋮_aNx−_arxaNy−_ary,b=12_K12−_Kr2−_r12+_rr2⋮_KN2−_Kr2−_rN2+_rr2,

and_Ki2=_aix2+_aiy2,i=1,…,N . The linear LS solution is simple, but its accuracy is relatively low compared to nonlinear LS method [6]. However, since a nonlinear LS is an optimization problem, its global solution may not be guaranteed while it requires a very high level of complexity.

To achieve a better accuracy at a reasonable cost of complexity, subspace localization using TOA [23] and RSSI measurements [24] is presented as follows. Assume that the WSN is fully cooperative, i.e., all the anchors and agents are in the communication ranges, or the anchors know the positions of each other. Define anN×2matrixXas

X=_a1x−_xmxa1y−_xmy⋮⋮_aNx−_xmxaNy−_xmy.

FromX, construct a multidimensional similarity matrix, denoted byD, as

D=X^XT.

Here,_Di,kelement is

_Di,k=(_aix−_xmx)(_akx−_xmx)+(_aiy−_xmy)(_aky−_xmy)=12(_di2+_dk2−_dik2),

where_diand_dkare the real distances between the i-th and k-th anchors to the m-th agent, respectively, while_dik2=^{∥_ai−_ak∥2},i,k=1,…,N. Applying the eigenvalue decomposition ofD, we have

D=UΛ^UT,

whereΛ=diag{_λ1,_λ2,…,_λN}is the diagonal matrix of eigenvalues ofDwith_λ1≥_λ2≥…≥_λN, andUis an orthonomal matrix whose columns are the corresponding eigenvectors_un,n=1,…,N. Since the rank ofDis 2,_λ3=_λ4=…=_λN=0 , therefore, Equation (9) can be written as

D=_{Us Λs1/2} ^{(_{Us Λs1/2})T}=_{Us Λs1/2}Ω^{(_{Us Λs1/2}Ω)T},

where_Us=[_u1u2],_Λs=diag{_λ1,_λ2}, andΩis an unknown rotation matrix satisfyingΩ^ΩT=_IN, where_IN denotes the identity matrix of order N. From Equations (7) and (10), we have

X=_{Us Λs1/2}Ω.

The unknown rotation matrixΩcan be determined based on the LS rule as

Ω=^{[(_{Us Λs1/2})T(_{Us Λs1/2})]−1 (_{Us Λs1/2})T}X.

Substituting Equation (12) into Equation (11), we obtain

X=_{Us UsT}X.

This means that the agents’ position_xmcan be extracted from the signal subspace eigenvectors. SinceDis not available, its approximation version, denoted byD^, can be constructed by

_D^i,k=12(_ri2+_rk2−_dik2),

where_riand_rk are calculated based on Equation (2). Hence, Equations (12) and (13) become approximations. As a result, we have

X≈_{U^s U^sT}X,

where_U^sis the signal subspace ofD^. DenotingA=_a1xa1y⋮⋮_aNxaNy,Xcan be written asX=A−_{1N xm}, where_1Nis theN×1 vector of all ones. Substituting this expression to Equation (15) and rearranging it, we have

(_IN−_{U^s U^sT})_{1N xm}≈(_IN−_{U^s U^sT})A.

Note that_IN−_{U^s U^sT}=_{U^n U^nT}, where_U^nis the noise subspace ofD^ . The LS solution of Equation (16) gives the agent’s estimated location, denoted by_xmS, as

_xmS=_{1NT U^n U^nT}A_{1NT U^n U^nT 1N}.

Figure 1 compares the root mean square errors (RMSEs) of the subspace method and the LS counterpart under different levels of shadowing and number of anchors. We can see that under the same levels of RSSI errors, the subspace method outperforms the LS method, especially when a larger number of anchors is deployed.

2.2. Hybrid AOA/RSSI Localization

Due to measurement errors, standalone RSSI localization cannot precisely locate agents. Therefore, combining the range measurement with an AOA estimation can significantly improve the accuracy. As shown in Figure 2, using both RSSI and AOA information, possible location of the agent (red star) is within the shaded area, which is smaller than that defined by standalone RSSI (3 circles) and AOA (3 lines). With both AOA and RSSI available at each anchor, we have a new set of linear equations as

_x^mx=_aix+_ricos(_αi)_x^my=_aiy+_risin(_αi),i=1,…,N.

Using the LS method, the solution for this problem is found [25]:

_x^mnAnR=^{(STS)−1 ST}u,

whereS=_1N0N0N1N;_1Nand_0Nare column vectors of N ones and N zeros, respectively; andu=^{[_a1x+_r1cos(_α1),⋯,_aNx+_rNcos(_αN),_a1y+_r1sin(_α1),⋯,_aNy+_rNsin(_αN)]T}. The accuracy of this method can be improved by applying a weight_wi=1−_{ri∑i=1N ri} on different links [19]. This is because a closer link normally has a smaller error than the longer one. The weighted least square (WLS) solution is expressed as

_x^mnAnR=^{(ST WTS)−1 ST WT}u,

whereW=_I2⊗diag{w},_I2denotes an identity matrix of order 2, ⊗ is Kronecker product, andw=_w1,⋯,_wN. Due to its simplicity and high accuracy, this WLS solution, referred to as “nAnR-WLS”, will be used as a benchmark in this paper.

3. The Proposed Method

Geometrically, conventional hybrid AOA/RSSI localization aims to reduce the area where the agent possibly resides, as shown in Figure 2. However, not all of the AOA measurements contribute to the improvement of the accuracy, especially, when AOA estimation suffers from high measurement errors. Since a high-resolution AOA estimation requires complicated signal processing and a large antenna array, conventional hybrid AOA/RSSI leads to complex anchors and so does the WSN.

To overcome this problem, we propose an unbalanced hybrid method, called 1AOA/nRSSI, as follows. Instead of estimating AOA from all anchors, we only use one AOA estimation to be combined with N RSSI measurements. As a result, in the WSN with the 1AOA/nRSSI method, only one master anchor, which is equipped with a large antenna array and powerful signal processing, is employed to obtain the AOA information. Other anchors are simple single-antenna transceivers to measure RSSI. In case a significant error is encountered in RSSI measurements, more simple anchors can be deployed to obtain the same reduction of the potential area as in the case of a fully hybrid AOA/RSSI method. This reduction is shown in Figure 3a.

There are many different ways of integrating one AOA and N RSSI measurements. Unlike conventional hybrid nAOA/nRSSI that enjoys N linear equations as in Equation (18), in the case of 1AOA/nRSSI, we only have one linear equation which is 1AOA/1RSSI, while N RSSI measurements are nonlinear ones. Therefore, N ranging measurements can be transformed into angular information, as proposed in [26]. In this transformation, given two anchors_a1,_a2with the ranges to the agent of_r1,_r2, draw a line that passes the intersections of the two circles centered at_a1and_a2with the radii_r1and_r2, respectively. This line is thus perpendicular to the line connecting these two anchors. Thereby, the vertices of the polygonal area in which the agent resides can be found from the ranging measurements_r1,_r2, and_d12. Then, the LS solution is applied to(N+1)linear equations. However, when the shadowing effect is significant, ranging-to-angular transformation leads to cumulative errors, yielding a lower precision.

We propose a simple method to integrate 1AOA with N RSSI measurements as follows. Without loss of generality, we assume that the master anchor is at the origin. From the AOA_α1and the range_r1estimated at the master anchor, we can determine two new virtual anchors, denoted as_av1,_av2, corresponding to the ranging information_rv1,_rv2 , respectively. As illustrated in Figure 3b, the virtual anchors’ locations and_rv1,_rv2can be calculated as

_av1=[_r1cos(_α1)0]_av2=[0_r1sin(_α1)]_rv1=_r1sin(_α1)_rv2=_r1cos(_α1).

We now have a set of(N+2) anchors and the corresponding ranging information. From this set, LS and subspace solutions presented in Section 2.1 can be used to locate the agents. We define these methods as “1AnR-LS” and “1AnR-Subspace” in the rest of the paper.

One problem of the LS method is to select the reference anchor (cf. Equation (4)). As presented in [27], the anchor that is closest to the agent should be chosen. However, in this set of(N+2)anchors, the master anchor at the origin is used three times, i.e., once as the real anchor_a1with the range_r1, and twice as virtual anchors_aviwith the ranges_rvi,i=1,2 . This means that the performance of the 1AOA/nRSSI localization relies significantly on the master anchor, especially when the number of anchors is small. As a result, choosing the master anchor as a reference node is a good decision in case of a small number of available anchors. This interpretation will be confirmed later in our simulations. The solution that employs the master anchor as the reference is denoted as 1AnR-LS. Its performance will be compared to that of 1AnR-LS-minRef [27]—which is the LS solution using the closest anchor as the reference—in the next section.

4. Performance Evaluation and Discussions 4.1. Simulation Setup

Simulations are conducted using MATLAB to evaluate and compare the precision of the proposed approaches and that of the conventional nAOA/nRSSI. A WSN with N anchors is randomly deployed in an area of 100 m × 100 m to localize 2000 agents whose locations are also random in the area. Among N anchors, the first one is fixed at the origin (_a1=[00]) which provides AOA estimation for the 1AOA/nRSSI method, while other anchors can make both AOA and RSSI measurements for the conventional AOA/RSSI positioning. To evaluate the accuracy of different localization methods, we use the expectation of RMSE, defined as

RMSE=E1M∑m=1M^(Δ_rm)2,

whereE{.}stands for expectation over Monte Carlo simulations, andΔ_rm=^{(_xmx−_x^mx)2}+^{(_xmy−_x^my)2} is the distance from the estimated location to the real location of the m-th agent. In the following simulations, 1000 Monte Carlo runs will be performed and 2000 agents will be estimated in each run. Since anchors are assumed to be deployed by UAVs, the air-to-ground channel model in the urban area evaluated in [28] is adopted in this paper. As shown in [28], with the transmitted signal frequencies from 1 GHz to 2 GHz, the PLE range is from 3.76 to 4.77. Further, we assume that all anchors have fixed positions so that the RSSI values are obtained after averaging many measurements to reduce the errors caused by multipath and shadowing. Therefore, without loss of generality, the averaged PLEγ=4is adopted at all anchors for simulations. The log-normal shadowing effect on RSSI at each anchor is generated independently in each iteration and dynamically for different agents. However, the log-normal shadowing variance_σni is assumed to be the same at all anchors. The angle estimation noise variances_σmi are set to be the same at all anchors in case nAOA/nRSSI is used.

4.2. Impact of AOA and RSSI Measurement Errors

We firstly evaluate the impact of angle estimation errors on the performance of the above localization methods. In this simulation, the AOA estimation error_σmi is set in the range of0to10degrees, while the log-normal shadowing is_σni =0dB.

Figure 4 shows the RMSE obtained by different methods withN=10anchors. We can see that—if perfect ranging measurement is obtained, i.e.,_σni =0dB—AOA estimation errors cause nAnR-WLS to be worse than both 1AnR algorithms. The reason is that when AOA estimations suffer from significant errors, such as under a non-line-of-sight (NLOS) environment, the more AOA information is involved and the worse the level of precision that can be obtained.

We then evaluate the impact of RSSI measurement errors on the performance of those methods. In the second simulation,_σni is set in the range from 0.1 to 1 dB, while_σmi is 0 degrees. Over the same number of runs and agents, the RMSE is presented in Figure 5. We can see that 1AnR-LS is much more sensitive to the range estimation error than the subspace one.

In the last two simulations, the RMSE is calculated under both AOA and RSSI measurement errors. The RMSEs are also compared to the Cramer–Rao lower bounds (CRLB) for both conventional nAOA/nRSSI and proposed 1AOA/nRSSI estimators, denoted by CRLB-nAnR and CRLB-1AnR, respectively. Detailed derivations of these bounds are presented in Appendix A.

Considering the AOA error_σmi =5degrees and the number of available anchors from 3 onward, the RMSEs for_σni =0.3dB and_σni =1 dB are depicted in Figure 6 and Figure 7, respectively. Under a small shadowing effect (_σni =0.3dB), the two 1AOA/nRSSI methods can achieve the same accuracy as the 3AOA/3RSSI approach with a slightly higher number of simple anchors. In particular, both 1AnR methods surpass the accuracy of 3A3R-WLS with 5 anchors. When more than 10 simple anchors are employed, the subspace method outperforms the fully nAnR-WLS. It is worth noting that the fully 3A3R-WLS and nAnR-WLS always require an antenna array and a more complex signal processing capability (thus power supply) at each anchor, while the 1AnR methods only require them at the master anchor.

The dashed blue line is 1AnR-LS-minRef, in which the reference anchor is chosen to have the shortest estimated range to the agent. In Figure 6, 1AnR-LS with the master anchor (at the origin) as the reference is superior to 1AnR-LS-minRef. However, in case of a more severe shadowing environment in Figure 7, 1AnR-LS-minRef outperforms 1AnR-LS. The reason is that, in the former case, the master anchor has a significant impact on the performance of the LS method as it is used three times compared to a small total number of anchors, as explained in Section 3. Meanwhile, in the latter case with more severe distance measurement errors, the reference anchor found from a higher number of simple available anchors is more reliable than the master one.

Under the more severe shadowing effect (_σni =1 dB) as shown in Figure 7, the performance of the 1AOA/nRSSI LS estimator is far worse than the subspace counterpart. 1AnR-Subspace needs up to 12 simple anchors to achieve the same accuracy as 3A3R-WLS; while in case of 1AnR-LS, the required number of simple anchors will be much higher (N>30) to reach the accuracy of 3A3R-WLS.

These simulation results confirm that the proposed 1AnR-Subspace and 1AnR-LS perform well in low-shadowing-effect environments with a slightly higher number of simple anchors (N=5) to estimate the agents’ location. In practice, one quick deployment method for these simple anchors is to fly some UAVs in certain trajectories through the expected anchor positions. The deployment of UAVs also allows the links from the anchors to the ground agents to have low shadowing. Meanwhile, in more severe shadowing environments, 1AnR-Subspace is more efficient as it requires a significantly lower number of anchors compared to 1AnR-LS, while its computational complexity is slightly higher than the LS counterpart.

It is worth noting that in both cases of the shadowing effects, CRLB-nAnR and CRLB-1AnR are almost the same. This means that the 1AOA/nRSSI estimator can ultimately achieve the same precision as the conventional nAOA/nRSSI—for example, when the number of anchors in both methods is very large—while the proposed method significantly reduces both hardware and computational complexities.

These simulation results imply that there is still room to improve the performance of the 1AOA/nRSSI. For example, more complicated algorithms may be used to achieve the same accuracy as the nAOA/nRSSI method, while the additional number of anchors remains small. In addition, as observed from Figure 6 and Figure 7, the accuracy can be improved by selecting an optimal position of the master anchor. As discussed in Section 2.2, the anchor that has the closest distance to the agent should be emphasized because it has smaller estimation errors. However, since all anchors are deployed at their fixed positions, while targets can appear anywhere in the area, it is impossible to select the master anchor with the shortest range to every target. To demonstrate the advantage of the optimal master anchor placement, we simulate a simple case of the master anchor placed at the center of the area. Since targets are uniformly distributed in the area, the average distance from the targets in the area to the master anchor is the shortest. Figure 8 depicts the RMSE for the same setup in Figure 6 but the master anchor is located at the center of the area, rather than at the origin. Clearly, under the same level of measurement errors, both 1AnR methods outperform the nAnR-WLS one whenN≥6compared toN≥10 , as shown in Figure 6.

Furthermore, the accuracy of the 1AOA/nRSSI estimator can be improved if the master anchor provides a more accurate AOA estimation. This means that instead of equipping all anchors with AOA measurements as in the conventional nAOA/nRSSI method, we can invest in the master anchor with a significantly large antenna array and sufficient signal processing capability, so that the 1AOA/nRSSI estimator has superior angle information than that in the nAOA/nRSSI. Hence, unbalanced hybrid AOA/RSSI localization allows a more flexible deployment of WSN. 5. Conclusions

This paper proposes an unbalanced hybrid AOA/RSSI localization method to simplify the WSN. By requiring only one AOA information to combine with RSSI measurements, the large antenna array and complicated signal processing required at all anchors nodes are not necessary except for one master anchor. Two different 1AOA/nRSSI localization methods are presented and evaluated under different scenarios. The simulation results show that the 1AOA/nRSSI LS estimator is suitable when more accurate RSSI measurements are available, while the subspace method is more suitable in a shadowing environment with an acceptable increase in the number of anchors and computational complexity. Possible future directions include improving the precision of 1AOA/nRSSI localization algorithms using different approaches listed in Section 4 and adopting an analog least mean square loop presented in [29,30,31,32] for a more accurate ranging estimation. Additionally, we might consider orthogonal frequency division multiplexing (OFDM) for correlated multipath fading channels [33] in the localization and positioning contexts.

Figure 1. Comparison of the subspace and least square (LS) solutions.

Figure 2. Hybrid nAOA/nRSSI localization. AOA-angle-of-arrival; RSSI-received signal strength indicator.

Figure 3. Unbalanced hybrid AOA/RSSI localization: (a) 1AOA/4RSSI and (b) Angular-to-ranging transformation.

Figure 4. Root mean square error (RMSE) versus AOA estimation errors;σni=0dB.

Figure 5. RMSE versus RSSI estimation errors;σmi=0degrees.

Figure 6. RMSE forσni=0.3dB,σmi=5degrees.

Figure 7. RMSE forσni=1dB,σmi=5degrees.

Figure 8. RMSE forσni=0.3dB,σmi=5degrees, with the master anchor at the centre of the area.

Author Contributions

Conceptualization, A.T.L. and L.C.T.; methodology, A.T.L. and L.C.T.; software, A.T.L.; validation, L.C.T. and X.H.; formal analysis, L.C.T. and A.T.L.; investigation, L.C.T. and A.T.L.; writing original draft preparation, A.T.L.; writing, review and editing, A.T.L., L.C.T., X.H., and S.L.P.; supervision, L.C.T. and X.H.; project administration, L.C.T., X.H., C.R., E.D.; funding acquisition, L.C.T., X.H., C.R., E.D., S.L.P., A.B., and D.F. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the New South Wales (NSW) Defence Innovation Network (DIN) and the NSW State Government, Australia, for the financial support of this project through the DIN Pilot Project Grant (Project ID: 888-006-985, Funding Years: 2019-2020).

Conflicts of Interest

The authors declare no conflict of interest. The approval for publications has been granted by the sponsor of this project.

Abbreviations

The following abbreviations are used in this manuscript:

AOA	Angle-of-Arrival
CRLB	Cramer-Rao Lower Bound
FIM	Fisher Information Matrix
GPS	Global Positioning System
LS	Least Square
NLOS	Non-Line-of-sight
ML	Maximum Likelihood
OFDM	Orthogonal-Frequency Division Multiplexing
PDF	Probability Density Function
PLE	Path-loss Exponent
RMSE	Root Mean Square Error
RSSI	Received Signal Strength Indicator
TOA	Time-of-Arrival
TDOA	Time-Different-of-Arrival
TOF	Time-of-Flight
UAV	Unmanned Aerial Vehicle
WLS	Weighted Least Square
WSN	Wireless Sensor Network

Appendix A. CRLB Derivation

Cramer-Rao Lower Bound (CRLB) is an important benchmark to compare different unbiased estimators. CRLB is calculated from a Fisher information matrix (FIM) which contains the expected value of the observed ranging information. Define a generalized measurement vectorθas

θ=f(x)+n

wheref(.)is a nonlinear function of the agents' positionsxandnis an additive zero-mean noise vector. Since hybrid AOA/RSSI information is in use,θincludes both RSSI and AOA measurements, i.e.,θ=^{[_θrss,1,...,_θrss,N,_α1,...,_αL]T}whereL=Nfor the case of the conventional nAOA/nRSSI andL=1for the case of 1AOA/nRSSI._θrss,i,i=1,...,N , is the RSSI represented in the Equation (A1) form, i.e.,

_θrss,i=_Pi-_P0=-10γ_log10(_di)+_ni.

Substituting_di=^(_xx-_aix)2+^(_xy-_aiy)2 into Equation (A2) and utilizing_αi calculated from the second equation in Equation (1),f(x)is written as

f(x)=-10γ_log10^(_xx-_a1x)2+^(_xy-_a1y)2⋮-10γ_log10^(_xx-_aNx)2+^{(_xy-_aNy)2tan-1}(_xy-_a1yxx-_a1x)⋮^tan-1(_xy-_aLyxx-_aLx).

The FIM, denoted byI(x), is the expected values of the second-order derivatives of the logarithm of the measurement probability density function (PDF) with respect tox, i.e.,

I(x)=E^∂2lnp(θ)∂x∂^xT

wherep(θ)is the PDF of the measurementsθ . As proved in [34], whennis zero-mean Gaussian distributed,I(x)can be also calculated by

I(x)=^{∂f(x)∂xT C-1}∂f(x)∂x

whereC=diag{_σn12,...,_σnN2,_σm12,...,_σmL2}is the covariance matrix of the noise vectorn.

Taking partial derivations off(x)over the two dimensions x and y, denoted as∂f(x)∂_xxand∂f(x)∂_xyrespectively, we have

∂f(x)∂_xx=^{-η_xx-_a1xd1,...,-η_xx-_aNxdN,-_xy-_a1yd12,...,-_xy-_aNydN2T}∂f(x)∂_xy=^{-η_xy-_a1yd1,...,-η_xy-_aNydN,_xx-_a1xd12,...,_xx-_aNxdN2T}

whereη=10γln10. Hence,I(x)is found:

I(x)=∑i=1N^{η2 (_xx-_aix)2}_{σni2 di2}+∑j=1L^(_xy-_ajy)2_{σmj2 dj4}∑i=1N^η2(_xx-_aix)(_xy-_aiy)_{σni2 di2}-∑j=1L(_xx-_ajx)(_xy-_ajy)_{σmj2 dj4}∑i=1N^η2(_xx-_aix)(_xy-_aiy)_{σni2 di2}-∑j=1L(_xx-_ajx)(_xy-_ajy)_{σmj2 dj4}∑i=1N^{η2 (_xy-_aiy)2}_{σni2 di2}+∑j=1L^(_xx-_ajx)2_{σmj2 dj4}.

Finally, the CRLB is obtained by

CRLB=E1M∑m=1M_{[^I-1(x)]1,1(m)}+_{[^I-1(x)]2,2(m)}2

whereE{.}denotes expectation over all Monte Carlo iterations.

References

1. Gonçalves Ferreira, A.F.G.; Fernandes, D.M.A.; Catarino, A.P.; Monteiro, J.L. Localization and positioning systems for emergency responders: A survey. IEEE Commun. Surv. Tutor. 2017, 19, 2836-2870.

2. Femminella, M.; Reali, G. A zero-configuration tracking system for first responders networks. IEEE Syst. J. 2017, 11, 2917-2928.

3. Han, G.; Xu, H.; Duong, T.Q.; Jiang, J.; Hara, T. Localization algorithms of Wireless Sensor Networks: A survey. Telecommun. Syst. 2013, 52, 2419-2436.

4. Win, M.Z.; Shen, Y.; Dai, W. A theoretical foundation of network localization and navigation. Proc. IEEE 2018, 106, 1136-1165.

5. Tomic, S.; Beko, M.; Dinis, R.; Bernardo, L. On target localization using combined RSS and AoA measurements. Sensors 2018, 18, 1266.

6. Zekavat, R.; Buehrer, R.M. Handbook of Position Location: Theory, Practice and Advances, 1st ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2011.

7. Nguyen, N.M.; Tran, L.C.; Safaei, F.; Phung, S.L.; Vial, P.; Huynh, N.; Cox, A.; Harada, T.; Barthelemy, J. Performance evaluation of non-GPS based localization techniques under shadowing effects. Sensors 2019, 19, 2633.

8. Wang, A.; Ji, X.; Wu, D.; Bai, X.; Ding, N.; Pang, J.; Chen, S.; Chen, X.; Fang, D. GuideLoc: UAV-assisted multitarget localization system for disaster rescue. Mob. Inf. Syst. 2017, 2017, 1267608.

9. Kang, S.; Kim, T.; Chung, W. Hybrid RSS/AOA localization using approximated weighted least square in Wireless Sensor Networks. Sensors 2020, 20, 1159.

10. Yu, K. 3-D localization error analysis in wireless networks. IEEE Trans. Wirel. Commun. 2007.

11. 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 26th Biennial Symposium on Communications (QBSC), Kingston, ON, Canada, 28-29 May 2012; pp. 24-29.

12. Yu, K.; Sharp, I.; Guo, Y.J. Ground-Based Wireless Positioning; IEEE Press Series on Digital & Mobile Communication; Wiley: Chichester, West Sussex, UK, 2009.

13. Biswas, P.; Lian, T.C.; Wang, T.C.; Ye, Y. Semidefinite programming based algorithms for sensor network localization. ACM Trans. Sens. Netw. 2006, 2, 188-220.

14. Qi, H.; Mo, L.; Wu, X. SDP relaxation methods for RSS/AOA-based localization in sensor networks. IEEE Access 2020, 8, 55113-55124.

15. Qi, Q.; Li, Y.; Wu, Y.; Wang, Y.; Yue, Y.; Wang, X. RSS-AOA-based localization via mixed semi-definite and second-order cone relaxation in 3-D Wireless Sensor Networks. IEEE Access 2019, 7, 117768-117779.

16. Chan, Y.T.; Chan, F.; Read, W.; Jackson, B.R.; Lee, B.H. Hybrid localization of an emitter by combining angle-of-arrival and receive 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.

17. Khan, M.; Salman, N.; Kemp, A.; Mihaylova, L. Localisation of sensor nodes with hybrid measurements in Wireless Sensor Networks. Sensors 2016, 16, 1143.

18. L. N. Nguyen, T.; D. Vy, T.; Shin, Y. An efficient hybrid RSS-AoA localization for 3D Wireless Sensor Networks. Sensors 2019, 19, 2121.

19. 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.

20. Jiang, J.R.; Lin, C.M.; Lin, F.Y.; Huang, S.T. ALRD: AoA Localization with RSSI differences of directional antennas for Wireless Sensor Networks. Int. J. Distrib. Sens. Netw. 2013, 9, 529489.

21. Popescu, D.; Stoican, F.; Stamatescu, G.; Chenaru, O.; Ichim, L. A survey of collaborative UAV-WSN systems for efficient monitoring. Sensors 2019, 19, 4690.

22. Wang, Y. Linear least squares localization in sensor networks. EURASIP J. Wirel. Commun. Netw. 2015, 2015, 51.

23. So, H.C.; Chan, F.K.W. A generalized subspace approach for mobile positioning with time-of-arrival measurements. IEEE Trans. Signal Process. 2007, 55, 5103-5107.

24. Lin, L.; So, H.C.; Chan, F.K.W. Multidimensional scaling approach for node localization using received signal strength measurements. Digit. Signal Process. 2014, 34, 39-47.

25. Sayed, A.H.; Tarighat, A.; Khajehnouri, N. Network-based wireless location finding. Signal Process. Mag. IEEE 2005, 22, 24-40.

26. Rupp, M.; Schwarz, S. An LS localisation method for massive MIMO transmission systems. In Proceedings of the 2019 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brighton, UK, 12-17 May 2019.

27. Guvenc, I.; Gezici, S.; Watanabe, F.; Inamura, H. Enhancements to linear least squares localization through reference selection and ML estimation. In Proceedings of the 2008 IEEE Wireless Communications and Networking Conference, Las Vegas, NV, USA, 31 March-3 April 2008; pp. 284-289.

28. Feng, Q.; McGeehan, J.; Tameh, E.K.; Nix, A.R. Path loss models for air-to-ground radio channels in urban environments. In Proceedings of the IEEE 63rd Vehicular Technology Conference, Melbourne, VIC, Australia, 7-10 May 2006; Volume 6, pp. 2901-2905.

29. Huang, X.; Guo, Y.J. Radio frequency self-interference cancellation with analog least mean-square loop. IEEE Trans. Microw. Theory Tech. 2017, 65, 3336-3350.

30. Le, A.T.; Tran, L.C.; Huang, X. Cyclostationary analysis of analog least mean square loop for self-interference cancellation in in-band full-duplex systems. IEEE Commun. Lett. 2017, 21, 2738-2741.

31. Le, A.T.; Tran, L.C.; Huang, X.; Guo, Y.J.; Vardaxoglou, J.C. Frequency domain characterization and performance bounds of ALMS loop for RF self-interference cancellation. IEEE Trans. Commun. 2019, 67, 682-692.

32. Le, A.T.; Tran, L.C.; Huang, X.; Guo, Y.J. Analog least mean square loop for self-interference cancellation: A practical perspective. Sensors 2020, 20, 270.

33. Le, N.P.; Tran, L.C.; Safaei, F. Energy-efficiency analysis of per-subcarrier antenna selection with peak-power reduction in MIMO-OFDM wireless systems. Int. J. Antennas Propag. 2014, 2014, 313195.

34. Chan, Y.T.; Ho, K.C. A simple and efficient estimator for hyperbolic location. IEEE Trans. Signal Process. 1994, 42, 1905-1915.