International Journal of Distributed Sensor Networks

Volume 2013 (2013), Article ID 794805, 11 pages

http://dx.doi.org/10.1155/2013/794805

Research Article

On the Joint Time Synchronization and Source Localization Using TOA Measurements

Ming Sun¹ and Le Yang²

¹Department of Electrical and Computer Engineering, University of Missouri, Columbia, MO 65211, USA

²School of Internet of Things (IoT) Engineering, Jiangnan University, Wuxi, Jiangsu 214122, China

Received 1 October 2012; Revised 1 January 2013; Accepted 7 January 2013

Academic Editor: Long Cheng

Copyright © 2013 Ming Sun and Le Yang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper considers the problem of estimating the clock bias and the position of an unknown source using time of arrival (TOA) measurements obtained at a sensor array to achieve time synchronization and source localization. The study starts with deriving the localization mean square error (MSE) for the case where we pretend that the source clock bias is absent and apply TOA positioning to find the source position. An upper bound on the clock bias, over which we shall obtain a higher localization MSE than that from jointly identifying the clock bias with the source position, is established. Motivated by the MSE analysis, this paper proceeds to develop a new efficient solution for joint synchronization and source localization. The new method is in closed-form, computationally attractive, and more importantly; it is shown analytically to attain the CRLB accuracy under small Gaussian TOA measurement noise. Computer simulations are conducted to corroborate the theoretical development and illustrate the good performance of the proposed algorithm.

1. Introduction

Recently, source localization using wireless sensor networks (WSNs) has attracted great research efforts [1–4]. To determine the source position, positioning parameters such as time of arrival (TOA), time difference of arrival (TDOA), received signal strength (RSS), and angle of arrival (AOA) are commonly used. In this paper, we shall consider determining the source position using TOAs of the source signal received at a wireless sensor network, which is essential for emerging applications including logistics, search and rescue, medical service, smart homes, environmental monitoring, and surveillance [5]. For this problem, a large number of algorithms such as those in [6–14] have been developed in literatures. Nevertheless, almost all of them assumed that the source and the sensors are synchronized ahead of the localization process. In other words, the clock biases of the source and the sensors are known a priori so that the obtained TOA measurements can be transformed into range measurement for source localization via TOA positioning.

In practice, the source clock bias may not be readily available and the TOA measurements are therefore subject to unknown source clock bias. In this case, we may simply ignore the source clock bias and still locate the source through TOA positioning. It can be expected by intuition that the obtained source localization accuracy would be worse than the case where the source clock bias is known a priori. However, the amount of performance degradation would be small if the source clock bias is negligible compared to the source-sensor distances. On the other hand, another way to handle the presence of the unknown source clock bias is to estimate it jointly with the source position. This would also yield a source localization accuracy poorer than that when the source clock bias is known a priori, due to the increase in the number of unknowns to be identified. Clearly, a theoretical analysis is needed here to compare the performance of the above two methods in terms of their localization accuracy. The obtained results would reveal the sensitivity of TOA positioning to the presence of the unknown source clock bias in the TOA measurement. More importantly, they would establish the conditions under which joint time synchronization and source localization should be preferred in order to achieve better source localization accuracy. To the best of our knowledge, this has not been addressed in existing literatures.

When the source clock bias cannot be ignored and/or it is also of interest besides the source position, an algorithm that can estimate both the source position and clock bias from the TOA measurements is definitely needed. For example, this could occur when a sensor is newly added to the WSN and it needs to be located and synchronized. Conventionally, synchronization and source localization are treated separately. In particular, synchronization is usually achieved first via applying one of the available protocol-domain techniques, such as [15–17], and the localization task can then be accomplished by executing various algorithms, such as those developed in [7–9, 12–14]. A more recent trend is to perform joint synchronization and source localization owing to their close relationship [18–22]. For this problem, given the statistical model of the TOA measurement noise, the maximum likelihood estimator (MLE) can be developed. However, since the cost function is nonconvex, MLE often resorts to the iterative search with a good initial guess for finding a globally optimal solution. This approach possesses the local convergence and even divergence problems. Another disadvantage is the high computational complexity due to possibly large number of iterations. Therefore, a closed-form solution is highly desirable. Bancroft's closed-form solution [23] is the first linearization-based algorithm for 1D timing (i.e., clock bias estimation) and 3D positioning. Zhu and Ding extended Bancroft's method in [21] to the case where the number of TOA measurements is more than the dimension of the unknowns (i.e., the source position plus the clock bias). However, this approach is in general not efficient. In other words, it can attain the Cramér-Rao Lower Bound (CRLB), the best accuracy for an unbiased estimator, in certain scenarios only.

The contribution of this paper is twofold. Firstly, we conduct a mean square error (MSE) analysis to derive the source localization MSE when the source clock bias in the TOA measurements is neglected and TOA positioning is performed to locate the source. For this purpose, a first-order analysis on a pseudo-MLE source location estimator is utilized, which indicates that the obtained results are valid for small TOA noise and source clock bias. The obtained source localization MSE would be partially dependent on the value of the source clock bias and as a result, it reflects the sensitivity of the TOA positioning accuracy to the presence of source clock bias. Moreover, an upper bound on the absolute value of the source clock bias is derived, over which joint synchronization and source localization would provide better localization accuracy in terms of smaller localization MSE than TOA positioning after ignoring the source clock bias. Secondly, this paper proposes a new efficient solution in closed-form for joint synchronization and localization. The newly developed method consists of two processing steps, where the first step locates the source and the second step estimates the clock bias. Compared with the iterative MLE method, the proposed solution is computationally more attractive, since it does not possess local convergence or divergence problem. Moreover, in contrast to the closed-form Bancroft's method [23] and its generalized version [21], the new approach is shown analytically to be able to attain the CRLB accuracy for both source location and source clock under small TOA noise. We support our theoretical developments using extensive computer simulations.

It is worthwhile to point out that the theoretical approach adopted in this paper to perform the MSE analysis and develop the proposed algorithm can be applied in a straightforward manner to the case when the TOAs are deduced from twoway/round-trip time of flight (TOF) measurements. TOF is usually considered as an inexpensive but effective alternative to obtain TOAs without using the cumbersome and costly source-sensor synchronization [24–26]. In this case, however, the turn-around time at the source needs to be either estimated [26] or small enough to be negligible [27, 28] so that TOA positioning can be utilized. With the above observation in mind, it is not difficult to show that the turn-around time in the TOF-deduced TOA measurements can be considered mathematically as an equivalent source clock bias. As a result, the MSE analysis and the proposed algorithm presented in this paper can be applied without much modifications to investigate the effect of ignoring the turn-around time and perform joint estimation of the turn-around time and the source position.

The structure of this paper is as follows. Section 2 formulates the joint time synchronization and source localization problem under Gaussian TOA noise model. Besides, the corresponding CRLB is established and analyzed in detail to motivate the MSE analysis and the new estimation algorithm development. Section 3 derives and verifies by computer simulations the localization MSE when the source clock bias in the received TOAs is ignored and the source position is found via TOA positioning. Section 4 presents the proposed closed-form solution for joint synchronization and source localization. Theoretical performance analysis that establishes the efficiency of the new method under small TOA measurement error is also conducted. Section 5 gives the simulation results to illustrate the good performance of the newly developed algorithm, and Section 6 concludes the paper and discusses the future work.

2. Problem Formulation and CRLB

In this section, we shall first formulate the joint time synchronization and source localization problem in consideration. The CRLB for the unknowns is then derived and analysed to motivate the MSE analysis and the algorithm development presented in the following two sections.

2.1. Signal Model

We shall consider jointly synchronizing and locating a single source in a 2D plane. The extensions of the theoretical developments in this paper to the more general case of synchronizing and localizing a source in the 3D space are straightforward.

The source is located at an unknown position . It has an unknown clock bias of seconds with respect to the reference time (i.e., the true time). Mathematically, the local time at the source would be if the reference time is . The value of can be negative or positive. Besides, is unknown but deterministic.

A sensor array composed of sensors at accurately known positions , , is used to identify and by measuring the TOAs of the source signal. Without loss of generality, we assume that the source starts emitting signals at its local time 0, which corresponds to the reference time . The source signal reaches sensor at the reference time , where is the true distance between the source and sensor , denotes the 2-norm and is the signal propagation speed. Let be the known clock bias of sensor with respect to the reference time. The TOA measurement obtained at sensor would be [21], after taking into account the measurement error ,

Multiplying both sides of the above equation with the signal propagation speed yields the TOA equation

where , , , and . Note that the scaled clock biases in (2), namely and , now have the units of meters. Noting that estimating the original source clock bias is equivalent to identifying its scaled version , we shall focus on determining in the rest of this paper to simplify the presentation.

Following the noise model in [21], we assume that the TOA measurement errors in (2) are independently and identically distributed (i.i.d.) Gaussian random variables with zero mean and variance . This model has been commonly adopted in literatures such as [10–14] for developing TOA-based localization algorithms and/or studying their performance via computer simulations. Nevertheless, it is a simplification of the real scenario where besides zero-mean Gaussian random errors due to the additive Gaussian noise at every sensor [10], the TOA measurements are also subject to errors owing to the multipath and nonline-of-sight (NLOS) propagation of the source signal as shown in recent experimental studies [29–31]. The fundamental limits of localization accuracy under these realistic factors was established in [32]. Extending our work to more realistic TOA noise models is beyond the scope of this paper, but it is under investigation.

By introducing the vectors,

we can express the measured TOA in (2) using the following signal model in matrix form:

The problem of interest is to estimate the source position vector and the source clock bias given and .

2.2. CRLB Analysis

The CRLB for the composite unknown vector under the joint time synchronization and localization scenario presented in the previous subsection is [21]

where is a unit vector from to . The CRLB of the source clock bias , denoted by CRLB, is given by the upper left component of CRLB while the lower right block of CRLB is the CRLB of the source position , denoted by CRLB. Applying the partitioned matrix inversion formula [33], we have

We proceed to evaluate (6b) more carefully to gain more insights. Collecting , , yields the partial derivative matrix given as

CRLB can then be expressed as

where , is an identity matrix and is an column vector whose elements are all equal to one as defined in (4). The matrix can be rewritten as the following matrix product:

where is an identity matrix and is an column vector whose elements are all equal to one. Applying the matrix inversion lemma [33] to the matrix in the middle of the right hand sides of (9) yields

Clearly, is an matrix with all its diagonal elements equal to 2 and others equal to 1. Substituting (10) back into (9) and putting the result into (8), we have

where is an matrix whose th row, , is equal to

Comparing (11) with from [34], we notice that within the considered time synchronization and source localization framework, the CRLB of the source position is equal to that of locating the source from time difference of arrival (TDOA) measurements , where . In other words, the TDOA measurements are indeed generated by subtracting the TOA obtained at sensor from the remaining TOAs. From the definition of the TOA measurements given in (2), all the TOAs are linearly related to the unknown source clock bias and hence, the subtraction operation eliminates the presence of in the obtained TDOAs . This enables TDOA positioning of the source and also forms the basis of the joint synchronization and localization algorithm proposed in Section 4 of this paper.

It is worthwhile to point out that we choose to produce the TDOAs through subtracting the TOA measured at sensor 1 from the rest TOAs for two reasons. First, this facilitates the comparison of the CRLB result derived in (11) with from [34], where the TDOA measurements were also obtained with respect to sensor 1. Second, subtracting the TOA at any sensor other than sensor 1 from the remaining TOAs would yield a different set of TDOA measurements but with the same quality as in terms of the source localization CRLB. This can be verified by applying the fact from Section 2 that the TOA measurement noises are i.i.d. and they have the same variance. Nevertheless, the above observation would become invalid if are no longer i.i.d. Investigating the approach that takes into account the statistical information on to produce the optimal TDOA measurement set is an important topic subject to future researches.

When the clock bias of the source, , is available a priori, but the source position remains unknown, the joint synchronization and localization task reduces to the classic problem of TOA positioning using the measurements in (2). The best possible localization accuracy for this case becomes [21]

and given in (6b) are the CRLBs of the source location estimation under two different scenarios, namely, when the source block bias is known and when is unknown. They are both matrices and have the same units, which makes subtracting in (13) from CRLB in (6b) feasible. In fact, it can shown that CRLB- is a positive semidefinite matrix, which indicates that compared with TOA positioning, the problem of joint time synchronization and source localization considered in this paper has a worse localization accuracy in general. This is expected, because we need to estimate one more unknown, specifically the source clock bias , from the same set of TOA measurements.

On the other hand, carefully examining (2) reveals that if the clock bias is negligible compared to the source-sensor distance , , or equivalently , and the source position is of primary interest, we may simply ignore the presence of and apply TOA positioning technique to the TOA measurements in (2) for identifying . In this way, the obtained location estimate will be biased, but it could have a smaller localization MSE than that of the location estimate from an efficient estimator for joint time synchronization and source localization. We shall theoretically illustrate the above observations in the following section.

3. Localization MSE Analysis

In this section, we shall derive the localization MSE when identifying the source position from the TOA measurements in (2) via TOA positioning and pretending that the source clock bias is zero. The obtained results will be contrasted with the best achievable localization MSE under the framework of joint estimation of and , which is equal to the trace of the CRLB given in (6b). The theoretical developments are based on an estimator that utilizes the Taylor-series linearization and estimates the source location from via the classic weighted least squares (WLS) technique. The reason behind the use of a WLS source location estimate is that it enables obtaining its localization MSE in closed form in terms of the TOA measurement noise power and the source clock bias to gain more insights. The localization MSE result in this section is valid for the scenario where the TOA measurement noise and the source clock bias are both small, due to the use of first-order approximation in the considered location estimator. But it applies to any TOA positioning technique that achieves the CRLB accuracy (13) when the source clock bias is known a priori. We shall corroborate the analytical results using numerical examples in Section 3.1 where the maximum-likelihood (ML) estimator for TOA positioning that attains (13) asymptotically when the source clock bias is known is simulated.

The MSE analysis starts with approximating the TOA measurement in (2) via applying the Taylor-series expansion up to the linear term to the source-sensor distance using an initial source position guess . Mathematically, we have

It can be deduced from (14) that the identification of the source position is equivalent to find . For this purpose, we rearrange (14) and stack the result for to arrive at

where and are defined in (4). Pretending that and utilizing being zero-mean Gaussian random vector, the maximum likelihood (ML) estimate for can be derived using the WLS technique [35]. Adding the result to the initial solution guess for yields the final source position estimate

Substituting (15) into the above equation and subtracting the true source position from both sides give the estimation error, when the clock bias is ignored

Assuming that the initial source position guess is sufficiently close to the true value , we can ignore the error in and approximate it with defined in (7) whose row vectors are . Therefore, (17) becomes

It can been observed that is a biased estimate of and the estimation bias is equal to , which is proportional to the value of the source clock bias . Moreover, the localization MSE can be derived by multiplying both sides of (18) with the transpose of and taking expectation as well as the matrix trace. We have

where is the trace of a matrix.

On the other hand, the localization MSE of an efficient estimator for jointly identifying and can be derived by applying the matrix inversion lemma [33] to (8) and evaluating its trace. Subtracting the result from (19) yields

This indicates that ignoring the clock bias may still lead to a location estimate with a smaller localization MSE than that of jointly identifying it together with the source position , if the following inequality holds:

The term on the right hand side of (21) is dependent on the TOA measurement noise power and the localization geometry. In particular, it increases with . Furthermore, it would also become relatively large when the source lies away from the sensor array. This is because in this case, the matrix inverse in the denominator, which is in fact the scaled version of the best TOA localization accuracy CRLB when the clock bias is known (see (13)), would increase. On the other hand, the term on the right hand side of (21) would decrease if is close to a zero vector. This could occur when the source is inside the sensor array.

3.1. Numerical Examples

We shall verify the theoretical analysis results via computer simulations. The considered scenario is depicted in Figure 1, where sensors are uniformly deployed along a circle centered at the origin and having a radius of 20m. The source is located either at m which is outside the sensor array or at m which is inside the sensor array.

Figure 1: Joint synchronization and source localization scenario. The sensors are denoted by circle symbols while the upper triangle symbols represent the two true source positions considered.

The MLE is realized to identify the source position and its localization MSE is defined as , where is the source location estimate in the th ensemble run and is the total number of ensemble runs. The noisy TOA measurements fed to the MLE in each ensemble run are generated by adding to the true values zero-mean Gaussian noise with variance set to be −20 in log scale. The source position is found via maximizing the pseudo-ML cost function with respect to the source position , where and are both defined in (4), collects the TOA measurements (having clock bias) and is functionally dependent on . In this way, the MLE performs the TOA-localization task by pretending that the clock bias in is absent. The maximization is conducted via applying the iterative Taylor-series method [36] with the initial solution guess produced by adding to the true source position zero-mean Gaussian random vector with a covariance matrix of 2CRLB, where CRLB is defined in (6b). This ensures that the MLE would normally need several iterations to converge to a solution.

Two sets of simulation results are generated, one for the source outside the sensor array and the other for the source inside the sensor array. They are plotted as function of the source clock bias in Figures 2 and 3. Also included in the figures are the theoretical localization MSE when ignoring the clock bias, which is given in (19). We plot as well the traces of CRLB in (6b) and CRLB in (13). They denote the lowest possible localization MSEs when the clock bias is jointly estimated with the source position and when is known a priori.

Figure 2: Localization accuracy for the source outside the sensor array when ignoring the presence of clock bias. from (13), from (6b), MSE from (19), star symbol: source localization MSE from the simulated pseudo-ML estimator.

Figure 3: Localization accuracy for the source inside the sensor array when ignoring the presence of clock bias. from (13), from (6b), MSE from (19), star symbol: source localization MSE from the simulated pseudo-ML estimator.

We can see from Figure 2 that when the source clock bias is zero, the localization MSE of the simulated ML estimator reaches CRLB as expected. As increases, the localization performance of the considered MLE degrades as predicted in (19). The simulation MSE matches the theoretical results very well, which justifies the validity of the localization MSE analysis. Furthermore, the localization MSE of the MLE remains to be smaller than that of jointly identifying the source clock bias and the source position, until the clock bias reaches a certain level (in this simulation, around  m). This is consistent with the analysis in (21). In Figure 2, the boundary condition (i.e., the crossing point of curves and ) that makes (21) achieve equality is also verified through simulation. Curves and are indeed the localization MSEs when the source is located jointly with the source clock bias and when the source position is found via ignoring the source clock bias and applying TOA positioning. The existence of a crossing point indicates that the two techniques can yield the same performance in terms of localization MSE for a particular value of the source clock bias.

The observations from Figure 3 where the source is inside the sensor array are very similar to those obtained from Figure 2 where the source is outside the sensor array. An important difference is that in this case, the crossing point of curves and appears when the source clock bias is equal to m. This observation is consistent with the analysis under (21). The crossing point being much closer to the origin compared with Figure 2 indicates that when the source lies inside the sensor array, ignoring the source clock bias and locating the source via TOA positioning is more likely to produce worse localization performance than joint clock bias and source position estimation, considering the presence of clock asynchronization (see, e.g., [32]). It is worthwhile to point out that the above observation is obtained without taking into account other realistic factors such as the multipath effect in signal propagation. We are currently extending our theoretical developments to more realistic signal models.

It can also be seen from Figure 3 that when the clock bias is as low as 0.05 m, which corresponds to around 0.17 nanoseconds, if the signal propagates at the speed of light, ignoring leads to an amount of more than 10 dB degradation in source localization MSE, compared to the joint time synchronization and localization. With those observations in mind and noting that existing closed-form joint synchronization and localization techniques, such as the one proposed in [21], are not able to provide the CRLB accuracy, we are motivated to develop in the next section a novel algorithm in closed-form that can identify and the source position efficiently.

4. Algorithm and Performance Analysis

In this section, we shall develop an efficient algorithm for jointly estimating the clock bias and the source position from the TOA measurements in (4). The algorithm is in closed-form and has low computational complexity. It consists of two processing steps, where Step 1 locates the source and Step 2 estimates . We shall show analytically that under small TOA measurement noise, the proposed algorithm reaches the CRLB accuracy for both the source position and the clock bias.

Step 1. The development of Step 1 processing is motivated by the CRLB analysis presented in Section 2.2. In particular, we have shown there that the CRLB of the source position, when jointly identified with the clock bias, is identical to that under TDOA positioning, where the TDOAs are produced by subtracting the TOA measurement obtained at sensor 1 from the remaining TOA measurements. Following this observation, we locate the source position as follows. First, we obtain the TDOA measurements via

From (2), the measurement errors in are . Collecting forms the column vector and its error vector is correspondingly

where is the TOA measurement noise vector defined in (4). is a zero-mean random vector because have zero mean and its covariance matrix is , where is defined in (10). Next, the source position is estimated from using the closed-form twostage algorithm developed in [34]. We shall briefly summarize the computations. Interested readers can find the algorithm details in [34].

Stage 1 of Step 1 processing estimates together with an extra variable that is the distance between the source and sensor 1. The functional relationship between them is ignored in Stage 1. Define the unknown vector as and its weighted least squares (WLS) estimate is found via

where the weighting matrix is equal to , and , , are the ranges between the source and sensor (see (2)). The matrix and the column vector both have rows and their th rows are and . The estimation error is equal to which is approximately zero-mean (under small noise condition and being approximately noise-free) and has a covariance matrix .

Stage 2 of Step 1 processing refines the source position estimate from Stage 1, which is , by exploring the relation between the unknowns and . The unknown vector of this stage is , where is the Schur product (element by element multiplication). Its estimate is again found via applying the WLS technique, which is given by

where the weighting matrix is and . and are equal to and , where is a identity matrix and is a column vector of ones. The estimation error in is with a covariance matrix of .

The source position estimate is obtained through mapping using

where is a sign function equal to 1 for and 0 otherwise. The source localization error would be , where . The covariance matrix of can be shown to be cov.

Step 2. Step 2 estimates the clock bias . In particular, under the joint synchronization and localization problem formulated in Section 2.1 (see (4)), the maximum likelihood (ML) estimates of and the source position can be found via solving the following minimization problem:

where the fact that the TOA measurement noise vector is a zero-mean Gaussian random vector has been applied. Differentiating the cost function in (27) with respect to and setting the result to zero yield the ML estimator for the clock bias . Mathematically, we have

is defined in (4) and it is unknown because the true source position is not available. We replace it with that has the same functional form as except that has been replaced by Step 1 output . Hence, the proposed algorithm estimates the clock bias via

We now summarize the processing required to estimate the source position and the clock bias . The proposed algorithm accomplishes the joint synchronization and source localization task by evaluating sequentially (24), (25), (26), and (29). Note that in computing and using (24) and (25), the true source position is needed to produce the weighting matrices and . To bypass these difficulties, when evaluating (24), we first set to be an identity matrix of appropriate dimension and compute to find . Then, (24) is calculated again with in being replaced by that is just obtained. In evaluating (25), the source position estimate from Stage 1 of Step 1 processing is utilized to generate .

4.1. Performance Analysis

We shall establish the efficiency of the algorithm proposed in the previous subsection under the condition that the TOA measurement noise are sufficiently small. Mathematically, we need to show that the algorithm output, namely the source position estimate and the clock bias estimate , is unbiased and their covariance matrices are approximately equal to the corresponding CRLBs, that is,

is the covariance matrix of given under (26) and CRLB is the CRLB of the source position defined in (6b). The validity of (30) can be verified by following the performance analysis procedure adopted in [34]. The unbiasedness of can be shown as follows. The estimation error in can be written as, from the definitions of , and given below (26), (25), and (24),

where is the error in the TDOA measurement vector produced from the TOA measurements , which is defined in (23) and shown to have zero mean. Under the small TOA measurement noise condition, we can ignore the noise in and as a result, is linearly proportional to and is zero-mean, which establishes that the source location estimate from the proposed algorithm is unbiased.

We proceed to prove (31), where is the covariance matrix of the clock bias estimate in (29) and CRLB is the CRLB defined in (6a). As pointed out above (29), the clock bias is identified within the proposed algorithm using an ML estimator. We have shown above that the source position estimate is an efficient estimate and as such, it can be expected from the property of the ML estimator [35] that the clock bias estimate would have an accuracy approximately equal to its CRLB. To derive , the estimation error in , denoted by , needs to be found. For this purpose, we expand in (29) at the true source position using the Taylor-series expansion up to the linear term, substitute (4), and subtract from both sides of (29) to arrive at

where is defined in (7). It can be shown by putting (32) into (33) that the clock bias estimate from the proposed algorithm is also zero-mean since its estimation error is linearly proportional to the zero-mean TOA measurement noise vector , under small TOA noise condition.

can be obtained via squaring both sides of (33) and taking expectation. We have, after some simplifications,

where (30) has been substituted and . We shall show that , where is an vector of zeros. In particular, after putting (32) and (23),

Putting the above result into (34) and comparing (34) with (6a) after applying the matrix inversion lemma [33] and , which is

would yield (31). This completes the establishment of the efficiency of the proposed algorithm.

5. Simulations

We shall demonstrate the estimation performance of the algorithm developed in the previous section for identifying the position and the clock bias of the unknown source in a twostep manner via computer simulations. The simulation scenario is the same as in Section 3.1, which is depicted in Figure 1.

The newly proposed algorithm in Section 4 is applied to estimate the source position as well as the clock bias fixed at m. The accuracy for source localization and synchronization is quantified using and , where is the total number of ensemble runs and and are the source position and clock bias estimates in ensemble . In each ensemble run, the erroneous TOA measurements are produced by adding to the true values zero-mean Gaussian noise with variance . For the purpose of comparison, we also realize the algorithm developed in [21] for joint source localization and synchronization.

For each of the two sources shown in Figure 1, we generate two figures for plotting the synchronization and source localization MSEs as function of the TOA measurement noise power , respectively. The corresponding CRLBs, specifically CRLB and CRLB from (6a) and (6b), are included in the figures as performance benchmark.

Figures 4 and 5 plot the clock bias and source position estimation MSEs for the source outside the sensor array. It can be seen from the figures that the proposed algorithm is able to attain the CRLB accuracy for the source position and clock bias before the TOA measurement noise power reaches in log scale. This is consistent with the theoretical performance analysis presented in Section 4.1 that the proposed algorithm is approximately efficient for accurate TOA measurements. On the other hand, the previously developed method from [21] cannot reach the CRLB accuracy, and its synchronization and localization MSEs are higher than the CRLB by an amount of more than 17 dB when lies in the range from −40 to −25 in log scale. As the TOA measurement noise power further increases over −10 in log scale, the performance improvement due to the use of the newly proposed method is even more significant.

Figure 4: Synchronization accuracy for the source outside the sensor array. Solid line: from (6a), star symbol: synchronization MSE from the proposed method, right triangle symbol: synchronization MSE from the algorithm in [21].

Figure 5: Localization accuracy for the source outside the sensor array. Solid line: from (6b), star symbol: source localization MSE from the proposed method, right triangle symbol: source localization MSE from the algorithm in [21].

Figures 6 and 7 depict the estimation MSEs for the clock bias and the position of the source inside the sensor array. Comparing with Figures 4 and 5 immediately reveals that the estimation performance is much better in this case, mainly because the localization geometry is improved. Moreover, in contrast to the case where the source is outside the sensor array, the proposed method and the algorithm from [21] yield similar estimation accuracy that matches the CRLBs for both the clock bias and the source position. This again verifies the performance analysis results in Section 4.1 on the approximate efficiency of the proposed solution. Interestingly, the method from [21] suffers from the threshold effect later than the new algorithm in this simulation. However, its better performance does not persist as the source location moves outsides the sensor array, as shown in Figures 4 and 5. In contrast, the proposed algorithm can attain the CRLB performance for both cases where the source lies within and outside the sensor array, when the TOA measurement noise power is not sufficiently large.

Figure 6: Synchronization accuracy for the source inside the sensor array. Solid line: from (6a), star symbol: synchronization MSE from the proposed method, right triangle symbol: synchronization MSE from the algorithm in [21].

Figure 7: Localization accuracy for the source inside the sensor array. Solid line: from (6b), star symbol: source localization MSE from the proposed method, right triangle symbol: source localization MSE from the algorithm in [21].

6. Conclusions and Future Work

In this paper, the TOA-based joint synchronization and source localization problem was considered. Effects of neglecting the presence of the source clock bias in TOA measurements on source location estimation accuracy was first investigated. For this purpose, an MSE analysis was performed for the case where the source is localized via TOA positioning when assuming the source clock bias does not exist, but in fact it is nonzero. Comparing the obtained source localization MSE with that from joint estimating the source position and clock bias, we derived a condition under which ignoring the source clock bias may provide a smaller localization MSE. Numerical examples were provided to validate the analysis and reveal that, in some cases, neglecting the clock bias can severely degrade the localization performance. As a result, a new efficient closed-form solution for joint synchronization and source localization was proposed. The new method can identify the source location and clock bias using a twostep approach. Theoretical performance analysis and simulations were conducted to show that it can achieve the CRLB accuracy for both source location and clock bias estimates under small Gaussian TOA measurement noise.

Some recent work such as [22] adopted a more completed model where both the time offset (clock bias) and skew are considered. However, the proposed least squares (LS) estimator in [22] cannot reach the CRLB accuracy. As a future topic, we would like to derive an efficient closed-form estimator for joint synchronization and source localization in the presence of time skew.

Acknowledgments

The authors contributed equally to this paper. L. Yang’s work is supported by the Startup Fund and the Youth Foundation of Jiangnan University (Contract no. JUSRP11234).

References

R. J. Kozick and B. M. Sadler, “Source localization with distributed sensor arrays and partial spatial coherence,” IEEE Transactions on Signal Processing, vol. 52, no. 3, pp. 601–616, 2004.

S. Gezici, Z. Tian, G. B. Giannakis et al., “Localization via ultra-wideband radios: a look at positioning aspects of future sensor networks,” IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 70–84, 2005.

T. Li, A. Ekpenyong, and Y. F. Huang, “Source localization and tracking using distributed asynchronous sensors,” IEEE Transactions on Signal Processing, vol. 54, pp. 3991–4003, 2006.

H. C. So, “Source localization: algorithms and analysis,” in Handbook of Position Location: Theory, Practice and Advances, S. A. Zekavat and R. M. Buehrer, Eds., chapter 2, pp. 25–66, Wiley-IEEE Press, Hoboken, NJ, USA, 2011.

M. Z. Win, A. Conti, S. Mazuelas et al., “Network localization and navigation via cooperation,” IEEE Communications Magazine, vol. 49, no. 5, pp. 56–62, 2011.

K. W. Cheung, H. C. So, W. K. Ma, and Y. T. Chan, “Least squares algorithms for time-of-arrival-based mobile location,” IEEE Transactions on Signal Processing, vol. 52, pp. 1121–1128, 2004.

P. Biswas, T.-C. Liang, K.-C. Toh, T.-C. Wang, and Y. Ye, “Semidefinite programming approaches for sensor network localization with noisy distance measurements,” IEEE Transactions on Automation Science and Engineering, vol. 3, pp. 360–371, 2006.

F. Chan, H. C. So, and W. K. Ma, “A novel subspace approach for cooperative localization in wireless sensor networks using range measurements,” IEEE Transactions on Signal Processing, vol. 57, no. 11, pp. 4548–4553, 2009.

H. Wymeersch, J. Lien, and M. Z. Win, “Cooperative localization in wireless networks,” Proceedings of IEEE, vol. 97, no. 2, pp. 427–450, 2009, Special issue on Ultra-Wide Bandwidth (UWB) technology and emerging applications

N. Patwari, A. O. Hero III, M. Perkins, N. S. Correal, and R. J. O’Dea, “Relative location estimation in wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 51, no. 8, pp. 2137–2148, 2003.

K. W. Cheung and H. C. So, “A multidimensional scaling framework for mobile location using time-of-arrival measurements,” IEEE Transactions on Signal Processing, vol. 53, no. 2, pp. 460–470, 2005.

A. Beck, P. Stoica, and J. Li, “Exact and approximate solutions of source localization problems,” IEEE Transactions on Signal Processing, vol. 56, no. 5, pp. 1770–1778, 2008.

M. Sun and K. C. Ho, “Successive and asymptotically efficient localization of sensor nodes in closed-form,” IEEE Transactions on Signal Processing, vol. 57, no. 11, pp. 4522–4537, 2009.

M. Sun, L. Yang, and K. C. Ho, “Accurate sequential self-localization of sensor nodes in closed-form,” Signal Processing, vol. 92, no. 12, pp. 2940–2951, 2012.

J. Elson, L. Girod, and D. Estrin, “Fine-grained network time synchronization using reference broadcasts,”,” in Proceedings of the 5th Symposium on Operating Systems Design and Implementation (OSDI '02), pp. 147–163, Boston, Mass, USA, December 2002.

S. Ganeriwal, R. Kumar, and M. B. Srivastava, “Timing-sync protocol for sensor networks,” in Proceedings of the 1st International Conference on Embedded Networked Sensor Systems (SenSys '03), pp. 138–149, November 2003.

M. Maroti, B. Kusy, G. Simon, and A. Ledeczi, “The flooding time synchronization protocol,” in Proceedings of the 8th International Conference on Embedded Networked Sensor Systems (SenSys '04), pp. 39–49, Baltimore, Md, USA, November 2004.

K. Römer and F. Mattern, “Towards a unified view on space and time in sensor networks,” Computer Communications, vol. 28, no. 13, pp. 1484–1497, 2005.

B. Denis, J. B. Pierrot, and C. Abou-Rjeily, “Joint distributed synchronization and positioning in UWB Ad Hoc networks using TOA,” IEEE Transactions on Microwave Theory and Techniques, vol. 54, no. 4, pp. 1896–1910, 2006.

H. Oliveira, E. Nakamura, and A. Loureiro, “Localization in time and space for sensor networks,” in Proceedings of the IEEE 21st International Conference on Advanced Information Networking and Applications (AINA '07), pp. 539–546, May 2007.

S. Zhu and Z. Ding, “Joint synchronization and localization using TOAs: a linearization based WLS solution,” IEEE Journal on Selected Areas in Communications, vol. 28, no. 7, pp. 1016–1025, 2010.

J. Zheng and Y. C. Wu, “Joint time synchronization and localization of an unknown node in wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 58, no. 3, pp. 1309–1320, 2010.

S. Bancroft, “An algebraic solution of the GPS equations,” IEEE Transactions on Aerospace and Electronic Systems, vol. 21, no. 1, pp. 56–59, 1985.

K. L. Noh, Q. M. Chaudhari, E. Serpedin, and B. W. Suter, “Novel clock phase offset and skew estimation using two-way timing message exchanges for wireless sensor networks,” IEEE Transactions on Communications, vol. 55, no. 4, pp. 766–777, 2007.

Z. Sahinoglu, “Improving range accuracy of ieee 802.15.4a radios in the presence of clock frequency offsets,” IEEE Communications Letters, vol. 15, no. 2, pp. 244–246, 2011.

Z. Sahinoglu and S. Gezici, “Enhanced position estimation via node cooperation,” in Proceedings of IEEE International Conference on Communications (ICC '10), pp. 23–27, May 2010.

M. R. Gholami, S. Gezici, M. Rydstrom, and E. G. Strom, “A distributed positioning algorithm for cooperative active and passive sensors,” in Proceedings of the 21st Annual IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC '10), pp. 1711–1716, September 2010.

M. R. Gholami, S. Gezici, E. G. Strom, and M. Rydstrom, “Hybrid TW-TOA/TDOA positioning algorithms for cooperative wireless networks,” in Proceedings of IEEE International Conference on Communications (ICC '11), June 2011.

B. Alavi and K. Pahlavan, “Modeling of the TOA-based distance measurement error using UWB indoor radio measurements,” IEEE Communications Letters, vol. 10, no. 4, pp. 275–277, 2006.

D. Dardari, A. Conti, U. Ferner, A. Giorgetti, and M. Z. Win, “Ranging with ultrawide bandwidth signals in multipath environments,” Proceedings of the IEEE, vol. 97, no. 2, pp. 404–425, 2009, Special issue on Ultra-Wide Bandwidth (UWB) technology and emerging applications

A. Conti, M. Guerra, D. Dardari, N. Decarli, and M. Z. Win, “Network experiment for cooperative localization,” IEEE Journal on Selected Areas in Communications, vol. 30, no. 2, pp. 467–475, 2012.

Y. Shen and M. Z. Win, “Fundamental limits of wideband localization—part I: a general framwork,” IEEE Transactions on Information Theory, vol. 56, no. 10, pp. 4956–4980, 2010.

L. L. Scharf, Statistical Signal Process, Detection, Estimation and Time Series Analysis, Addison-Wesley, Reading, Mass, USA, 1991.

Y. T. Chan and K. C. Ho, “A simple and efficient estimator for hyperbolic location,” IEEE Transactions on Signal Processing, vol. 42, no. 8, pp. 1905–1915, 1994.

S. M. Kay, Fundamentals of Statistical Signal Processing, Estimation Theory, Prentice Hall, Englewook Cliffs, NJ, USA, 1993.

K. C. Ho, X. Lu, and L. Kovavisaruch, “Source localization using TDOA and FDOA measurements in the presence of receiver location errors: analysis and solution,” IEEE Transactions on Signal Processing, vol. 55, no. 2, pp. 684–696, 2007.