Abstract

Target localization plays an important role in the application of radar, sonar, and wireless sensor networks. In order to improve the localization performance using only two stations, a hybrid localization method based on angle of arrival (AOA) and time difference of arrival (TDOA) measurements is proposed in this paper. Firstly, the optimization model for localization based on AOA and TDOA are built, respectively, in the sensor network. Secondly,the majorization-minimization (MM) method is employed to create surrogate functions for solving the multiple objective optimization problem. Next, the hybrid localization problem is solved by the projected gradient decent (PGD) method. Finally, the Cramer–Rao lower bound (CRLB) for the joint AOA and TDOA method is derived for the comparison. Simulations proved that the proposed method has improved localization performance using AOA and TDOA measurements from only two base stations.

1. Introduction

Source localization in the sensor network is a fundamental problem which has received an upsurge of attention in recent years. A number of separated sensors in the network will be employed to measure the emitted or reflected signal from the target [1]. Many algorithms have been presented in literature to determine the source position, based on the time of arrival (TOA), time difference of arrival (TDOA), angle of arrival (AOA), and received signal strength (RSS) method.

Localization problem based on the TOA method in the quasi-synchronous network is addressed in [2], in which a two-step linear algorithm is presented for estimating the passive object position. TOA-based localization problem in an adverse environment is solved in [3], which employs two convex relaxation methods for mitigating non-line-of-sight (NLOS) effect. Using a single transmitter and multiple receivers, the optimal geometries for TOA localization configuration in the two-dimensional environment are analyzed in [4]. With the help of floor plan information, Liu et al. [5] propose an effective geometrical localization method based on TOA in the non-line-of-sight (NLOS) environment without considering multiple-bound scattering signal paths.

A coarse position estimation algorithm-based TDOA is presented in [6], which is robust with regard to the initial position without redundant receivers. As a further step, a low-complexity nonlinear expectation maximization localization algorithm is presented in [7]. Target location is estimated by the TDOA data in [8], where empirical mode decomposition is used to decompose seismic signals, and multiple features of intrinsic mode functions are extracted. The work of [9] transforms a TDOA model into a TOA model, by presenting a semidefinite programming method. The TDOA method is also employed in passive radar source localization [10], taking advantage of the characteristics of the same sensor in different positions.

For obtaining the best localization performance, the angular separation requirements for AOA sensors are established in [11]. A closed-form solution for 3D localization using AOAs is presented in [12] to handle the effect of sensor position errors in a wireless sensor network. In the AOA localization system, an evaluation function by frame theory is derived and proved in [13], in order to solve the sensor augmentation problem.

Localization algorithm based on the RSS method in NLOS environment is presented in [14], by correcting the NLOS measurements. A novel fingerprinting paradigm is introduced for RSS localization in [15], which obtains a significant gain for multiple sources and multipath environments. In [16], RSS localization methods based on convex optimization are presented to solve the noncooperative and cooperative problems in sensor networks.

However, localization based on the TOA method needs an accurate time synchronization. TDOA-based approach requires at least four located sensors. RSS algorithms are based on known exact received signal without multipath effect, and the AOA method can perform well only when the target is not far away from the sensors [17]. There are many advantages with the hybrid localization method [17, 18], which has already attracted many attentions. Using the hybrid-bearing AOA and TDOA sensors [19], the constrained optimization method is employed to estimate the maximum likelihood values. The RSS-assisted TOA-based localization method is presented in [20], which improves the efficiency and localization accuracy of the indoor geolocation system. In adverse NLOS environment, a new mitigation method is proposed using the combined RSS and TOA measurements in [21]. Combining TDOA and RSS measurements, a two-step position estimator is presented in [22] for the visible light positioning system. The hybrid AOA- and RSS-based localization method is presented in [23] for 3D wireless sensor network without central processor. Using the same hybrid method, a simple closed-form solution method is proposed by using the spherical coordinate conversion in [24]. Using the hybrid measurements of TDOA and AOA methods, the mean square error matrix is derived under small error condition in [25], and a closed-form solution based on two base stations (CFS-2BSs) is proposed. In order to reduce the sensors for target localization, we propose a hybrid TDOA/AOA localization approach only with two sensors, taking advantage of the complementary property between the TDOA- and AOA-based localization methods.

This paper investigates the localization problem based on AOA and TDOA measurements from two stations. The main contributions of this paper are summarized as follows:(1)A localization model using only two stations is developed, which is converted to be a multiple optimization problem.(2)The MM method is used to solve this optimization problem, which maximizes the similarity between unknown and known information for AOA and TDOA measurements.(3)A single objective optimization model is constructed after the derivation based on the MM. The PGD method is employed to solve this problem effectively.(4)Simulations prove the superior localization performance using only two stations in the sensor network.

This paper is organized as follows. Section 2 describes the basic TDOA localization model. Section 3 proposes the hybrid localization method based on AOA and TDOA data from two base stations. CRLB is derived for evaluating the localization performance in Section 4. Section 5 shows the simulation results based on the proposed method. Section 6 gives the conclusion of this paper.

Notations. Bold uppercase (e.g., H) and lowercase (e.g., b) letters represent the matrices and vectors, respectively. The notations and stand for transpose and Hermitian of their argument, respectively. denotes the norm of a vector. The gradient of at is denoted by .

2. Basic Localization Methods Based on TDOA

Let denotes the position of th base station and represents the unknown position of the master station. The distance between the th base station and the master station can be written aswhere is the number of base stations. Time difference from the th and 1th base station can be defined aswhere c is the propagation speed. Distance difference can be derived as

Taylor iteration method and Chan’s method are proposed, respectively, as the classical localization methods based on the time difference.

An iterative method is employed to solve the hyperbolic equations for estimating the MS position in the Taylor iteration method, where a local linear least square (LS) solution is selected based as the deviation of each iteration. The position of master station can be written aswhere is initial coordinate values for the iteration and and can be derived aswhere is the covariance matrix for the time difference estimation and and are written as, respectively, as follows:

Besides, Chan’s method is proposed to obtain a noniterative solution. Let be the solution, which can be derived aswhere

3. Hybrid Localization Method Based on AOA and TDOA

3.1. Multiple Objectives’ Design

The unknown location of the target from the noisy measurements of angle of arrival (AOA) can be estimated, as there is a nonlinear relationship between the AOA and target location. Only two base stations with two angles are considered, , which is assumed to be known in Figure 1.

In Figure 1, and are the azimuth angles of the base stations. can be written as follows:

Based on the AOA measurements of every BS, the first fitness function is designed as follows:where the vector is designed as follows:

Using the TDOA data from the two base stations, the sum of squared range difference errors is used to design the second fitness function, which is defined aswhere

3.2. Hybrid Optimization Model

In order to solve the hybrid optimization problems based on (10) and (12), the minorization-maximization method is employed to create the surrogate functions here. Let , where can be majorized by , which is derived as follows:where is the solution at the th iteration, and

Let ; then,

The gradient of can be written as

Let can be majorized by , which can be written as follows:where

Based on (14) and (18), the final objective function for the localization using two stations is designed aswhere is the scalarization coefficient, and

Then, we use the PGD method to solve the optimization problem. can be calculated iteratively bywhere is the Euclidean projection. The gradient of can be written as

In conclusion, the hybrid localization procedure based on the TDOA and AOA data is described in Algorithm 1.

(i)Input: Information of AOA and TDOA, Number of iterations
(ii)Result: MS position
(1) Set ; start from one random position
 (2) repeat
  (3) Construct surrogate functions by using formulas from (22) to (25).
  (4) Design objective function (20)
  (5) Set and to some initial vector, e.g.,
  (6) repeat
   (7) Compute using (27)
   (8)
  until;
  (9) Obtain
  (10)
  until;

4. Derivation of Cramer–Rao Lower Bound

Only Gaussian range sampling error is assumed in the localization here, so Fisher’s information matrix [26] is then given bywhere is defined as the inverse matrix of :where and are, respectively, the element of and  = 1, 2 and  = 1, 2. The localization CRLB based on the AOA method can be computed as follows:

Based on the AOA measurements, (31) can be derived aswherewhere is the Gaussian error variance of and denotes the angle between the MS and th BS with respect to a reference direction.

Using the same method, CRLB based on the TDOA method can be calculated aswherewhere , where is supposed to contain a Gaussian error with average and variance of 0 and , respectively.

The answer to the multiobjective optimization problem is a set of solutions that defines the best trade-off between competing objectives. Here, the log-likelihood function for the multiobjective optimization problem can be written aswhere and .

The CRLB for the localization based AOA and TDOA using two base stations can be derived as

5. Simulations

To verify the effectiveness of the proposed approach, computer simulations are conducted for comparisons with other methods. As it is shown in Figure 2, the coordinates of the BSs are, respectively, set to , , , and ; the MS location is chosen randomly according to a uniform distribution defined over the region formed by the BSs.

In our proposed localization algorithm, and are selected as the two stations. Then, 15 TDOA and AOA measurements from the selected BSs are obtained. The root of mean squared error (RMSE) of MS estimation from the th BS is given bywhere denotes the true position of MS.

5.1. Comparison with the Classical TDOA Localization Methods Using More BSs

As we know, it is necessary for the classical TDOA-based methods to employ more than 3 stations realize the localization. To start with the comparison of the localization methods, the results of the classical source localization method are presented in Figure 3. Based on Taylor iteration method and Chan’s method, it can be observed that the estimated value is much closer to the sensor position whenever more BSs are employed. This is because Taylor iteration and Chan’s methods apply the iterative and noniterative solutions to achieve their optimum performance, respectively. Thus, they are considerably affected by the accuracy of TDOA measurements. Under the same TDOA measurement condition, it can be observed that the proposed algorithm has much better localization accuracy than the classical methods using three BSs, while it uses only two BSs. In some cases, e.g., Chan’s method, using three BSs, a superiority of more than one order of magnitude is observable. Since increasing the number of BSs is expected to improve accuracy, reformulating the proposed method to a higher number of BSs is expected to result in even higher accuracy. However, the reformulation is beyond the topic of this article and can be addressed in future investigations.

5.2. Comparison with Single Objective Optimization-Based Methods

In order to prove the effectiveness of the proposed method using only two stations, source localization results using 3 BSs are presented here based on optimization methods, including weighted least squares (WLS-O-3BSs) [7], and optimization using time difference of arrival (TDOA-O-3BSs) [21]. They are then compared with the proposed methods. Figure 4 depicts localization error for a various number of testing points. It can be observed that the results obtained in this work outperform others, being able to achieve a lower localization error. Specifically, an improvement of the order of one degree of magnitude is observable in terms of peak localization error, i.e., the highest deviation in estimation.

Using two stations, the TDOA (TDOA-O-2BSs), AOA (AOA-O-2BSs) based on equations (10) and (12), and CFS-2BSs in [25] are realized, respectively. The results are demonstrated in Figure 5. Due to the lack of measurement information, there are large errors in the TDOA and AOA measurements, and the single objective function for optimization cannot perform very well in such an environment. CFS-2BSs, in [25], proves its good localization performance through the simulations, but its performance depends on the covariance matrix of the measurement noise. Interestingly, the proposed method shows very low localization error in a broad range for testing point numbers. When all methods use 2 base stations, the achieved accuracy improvement is about two orders of magnitude compared to TDOA and less than one order of magnitude compared to the AOA method in terms of peak localization error.

5.3. Comparison with Other Multiobjective Optimization-Based Methods

There are many optimization methods for the multiple objectives, such as the conventional multiobjective particle swarm optimization (MOPSO) [27] and nondominated sorting genetic algorithm (NSGAII) [28, 29], which can be used to solve the localization problems using two stations. Here, the proposed algorithm is compared with the MOPSO and NSGAII algorithms to demonstrate its superiority over a range of state-of-the-art investigations. In this regard, the population is set to 100 for all three optimizers, and the maximum number of iterations is set to 200. Four simulations are conducted using different optimization methods for the same multiple objectives in equations (10) and 12), where their results are demonstrated in Figure 6. Besides, the percentages of crossover and mutation populations and the mutation rate of NSGAII are set to 0.7, 0.4, and 0.02, respectively. The root CRLB in equations (34) and (38) is also computed and depicted in the same figure, labeled as AOA-CRLB and TOA-AOA-CRLB, respectively. Under the same conditions, the localization RMSE of the proposed method is smaller than other methods, being significantly close to CRLB.

Based on the theory of projected gradient descent, our proposed converges at rate , where k is the iteration number. The runtime comparison results using the same computer are depicted in Figures 7 and 8. The experiments confirm the effectiveness of the proposed MM-based method in solving optimization and localization problems under the uncertainties of MS positions. This is because the nondominated sorting GA (NSGAII) must sort out all the individuals according to the Pareto-domination relationship and selects the individuals with better ranks to form the next-generation population, but with a lower convergence speed [30]. For conventional MOPSO, the neighborhood of one particle is formed by some nearest particles according to the objective space of one objective, but the best particle in the neighborhood cannot be determined very well by the fitness of the other objective. The CFS-2BSs method proposed by constructing relationships between the hybrid measurements and the unknown source position is not the iterative implementation method, which takes less computations.

6. Conclusion

In this paper, a significant increase in the positioning accuracy is achieved by using the MM and PGD methods, and the TDOA and AOA information together, while decreasing the number of required base stations. Accordingly, the localization problem is formulated for a two station scenario and is then derived by the MM method and then solved by the PGD method. Then, theoretical CRLB analysis is realized to establish a basis for the accuracy comparison. The extensive performance comparison with several state-of-the-art localization methods and some other optimization methods reveal the superiority of the proposed method. Specifically, they show that the proposed algorithm can achieve a superior localization accuracy using only two base stations compared to other recent methods required three or more base stations.

Data Availability

We note that there are no data-sharing issues since all of the numerical information is produced by solving the equations of the proposed algorithm, which are realized by MATLAB software in the paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Fund (61671221) in China.