Academic Editor: A. G. Constantinides
Copyright © 2008 Hong Tang et al. 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
Wireless location becomes difficult due to contamination of measured time-of-arrival (TOA) caused by non-line-of-sight. In this letter, TOA measurements seen at base stations are adjusted by scale factors, and a modified deterministic model is built. An effective numerical solution is proposed to resolve the scale factors and mobile position. A simulation comparison of four algorithms indicates that the proposed algorithm outperforms the other three algorithms.
1. Introduction
For outdoor
wireless location, most time-of-arrival- (TOA-) based algorithms are derived
from deterministic model. The mobile position is commonly modeled as the
intersection of a set of circles defined by TOA ranges. Non-line-of-sight (NLOS)
propagation is one of the major concerns. The field trials in [1] performed by
Nokia showed that the positive range error due to NLOS is up to a thousand
meters in some channel environments. NLOS is a key problem to improve location
accuracy.
In recent years,
very important contributions to the TOA-based location have been made to
compute the mobile position using TOA measurements. And the least-square (LS) criterion
is often used to compute the estimate of mobile position, which is optimal if
the residual follows Gaussian distribution. For example, the “time of arrival
data fusion" based on LS criterion in [2] is a closed-form solution for mobile
position. A constrained nonlinear LS approach was proposed in [3] to minimize
the residual. In [4], the quadric programming technique is used to find the ML
estimation of the mobile position subjected to constraint from the inequalities
introduced due to NLOS propagation. In the range scaling algorithm (RSA) in [5],
the NLOS is mitigated using TOA adjustment by scaling NLOS-corrupted TOA
measurements using factors that are estimated from a constrained nonlinear
optimization process. But RSA is only valid for three base stations (BSs)
scenario, and the optimization process is very complex.
In this
letter, scale factors are employed to adjust TOA measurements seen at base
stations, and a modified deterministic model is proposed. An effective solution
is used to resolve the mobile position, which does not require prior knowledge
of NLOS distribution or time-based history measurements.
2. Deterministic Model for TOA-Based Location
It is assumed
that
BSs are involved in the location process. The
TOA-based range measurements are modeled as
(1) where
is the pseudorange between the
BS and mobile station (MS),
is the LOS distance, and
is the excess distance due to NLOS.
is the sum of other errors caused by system
delay, synchronization error, measurement noise, and so forth. Let the mobile
position be
and let the position of the
BS be
.
Define
as the residual, where
denotes the norm operator, and
represents the distance between vectors
and
. A generalized estimator normally used for TOA in determinate models is
(2) where
is the index set of BSs.
is an optimal function of the residual. The determinate model in [3] was chosen as
where
is the weight reflecting the
reliability of the signal received at BS
,
and
in [6] was selected as a robust function with
more tolerance for outliers. To locate mobile position in presence of NLOS
propagation, an alternative way is to shift view to
in a modified deterministic model.
3. A Modified Deterministic Model and Numerical Solution
The fact that
NLOS propagation causes positive range error means that the following
inequality holds:
(3) In order to change (3) into equality, let
be the scale factor on the range measurements,
which yields
(4) Each
must be less than or equal to one.
Furthermore, each
has a lower bound. For example, three BSs are involved
in the location process, as seen in Figure 1. The mobile must lie in the
overlapped region formed by
and
.
Let
be the lower bound of
.
We find that
.
In a similar way, we get
and
.
Thus,
must be constrained by
(5)
Figure 1: The geometry of three
circles shows the constraints on the scale factors.
Define
as the modified residual in this letter. This residual
will be zero or sufficient small if the scale factors are selected accurately. Thus,
the modified deterministic model is defined as
(6) Subject to
(7) Let the objective
function which is minimized for TOA location be
(8)
Clearly, (8)
will be exactly zero with the true scale factors and true position. The
sequential quadratic programming (SQP) technique is employed to study the
performance surface of the objective function. We find that it has many local
minimums in the overlapped region, and these minimums skirt the true mobile
position.
Because of the
difficulty in determining a closed-form solution for (8), we propose a numerical
solution to resolve the scale factors and mobile position by the following
steps.
(1)
Compute
the lower bound
for
.
(2)
Generate a random variable
over
with a predefined distribution. The predefined
distribution of
is different in different channel environments.
Typically, the channel environments can be classified as open area, mountains,
suburban, urban, and bad urban.
(3)
Let the
scale factors equal
for
and compute the mobile position as follows.Equation (8)
is quadratic about
and
.
Calculate the derivatives with respect to
and
,
and let the derivatives be zero. The estimate of the mobile position is
(9) An interesting result from (9) is that the mobile position estimate
is the linear combinations of the positions of all BSs. The tap weight of each
linear combination is
.
However, we must note that (9) is not the closed-form solution because
and
are still functions of
and
.
Following the similar steepest descendent technique as illustrated in [3], the estimate
of mobile position is
(10) where
and
are the step sizes, and
is the iteration number.
(4)
Substitute the scale factors
and mobile position estimate
into (8), and calculate the exact value of
.
(5)
Repeat
steps, (2)–(4)
times. We get a set of the scale factors
,
a set of mobile position
and a set of exact value of objective function
.
(6)
Find
the minimum
in
,
where
is the index. Then,
in
is the final estimate for mobile position, and
in
are the final estimates for scale factors.
The random variable
controls the step size of the scale factors
while we search them in solution space. Without any prior knowledge, we apply
the same
on each scale factor in step (3) to limit the
solution space. In step (6), both the scale factors and the mobile position that minimizes
the objective function are taken as the final estimates. This scheme to select
scale factors is reasonable because it follows the physical meaning of the modified
deterministic model. Clearly, the final estimates cannot make the objective
function exactly zero. But, they can guarantee it to be a relatively small
number. A notable advantage of the numerical process is that compared to the
complex nonlinear optimization processes as illustrated in RSA [5], the
complexity to search the scale factors is reduced. In addition, the proposed
algorithm can be applied to the scenarios of more than three BSs; however, RSA is difficult to be applied
in these scenarios.
4. Simulations
The
simulations are conducted in a cellular network composed of hexagonal cells
with radius of 1000 meters. The mobile station is assumed uniformly distributed
in the central cell. And the mobile position is computed based on TOA
measurements from the six nearest BSs. In scenario 1, the excess range
due to NLOS propagation is modeled as positive
random variables over
meters, generated according to CDSM [5, 7]. To
study the performance in a high NLOS environment, the “reverse” CDSM is used. Without
knowledge of the channel environment, the variable
in the proposed algorithm is uniformly generated
over
,
and the number of attempts is 10, that is,
.
The range error
is caused by system delay, synchronization
error, measurement noise and is typically smaller in magnitude than NLOS error.
As an example,
is modeled as a Gaussian variable with mean
100 meters and variance 30 meters. The performances are compared with the LS
algorithm, that is, “time of arrival data fusion" algorithm in [2], the
algorithm in [3], and the Huber-estimator mentioned in [6]. The root mean
square errors of the four algorithms with the two NLOS models are shown in Figure 2.
Figure 2: The root mean square errors
with CDSM and reverse CDSM.
In scenario 2, another model
considered is a range-dependent NLOS model [8]. The NLOS error is proportional
to the BS-MS range, that is, the NLOS error for the
th BS is taken to be
,
where
is a proportionality constant. The root mean
square errors with the proportional NLOS model are shown in Figure 3. It is
found that the performance of the other three algorithms degrade greatly with
increasing
and the proposed algorithm is robust with
.
Figure 3: The root mean square errors with
range-dependent model.
5. Conclusions
A modified
deterministic model with a new residual is built to compute the mobile position.
Due to the difficulty to find a closed-form solution, an effective numerical
solution is proposed to resolve the scale factors and the mobile position. A
simulation comparison among four algorithms is conducted to evaluate their
performance. The results indicate that the proposed algorithm can deal with
NLOS effectively and is simple to be implemented.
Acknowledgments
The authors are very grateful to the anonymous reviewers
for their useful comments. The authors wish to thank the Regional
Innovation Center (RIC) of Yeungnam University partially
funding this research.
References
- M. I. Silventoinen and T. Rantalainen, “Mobile station emergency locating in GSM,” in Proceedings of the IEEE International Conference on Personal Wireless Communications (ICPWC '96), p. 232, New Delhi, India, February 1996.
- A. H. Sayed, A. Tarighat, and N. Khajehnouri, “Network-based wireless location: challenges faced in developing techniques for accurate wireless location information,” IEEE Signal Processing Magazine, vol. 22, no. 4, 24 pages, 2005.
- J. Caffery, Jr. and G. L. Stuber, “Subscriber location in CDMA cellular networks,” IEEE Transactions on Vehicular Technology, vol. 47, no. 2, 406 pages, 1998.
- X. Wang, Z. Wang, and B. O'Dea, “A TOA-based location algorithm reducing the errors due to non-line-of-sight (NLOS) propagation,” IEEE Transactions on Vehicular Technology, vol. 52, no. 1, 112 pages, 2003.
- S. Venkatraman, J. Caffery, Jr., and H.-R. You, “A novel ToA location algorithm using LoS range estimation for NLoS environments,” IEEE Transactions on Vehicular Technology, vol. 53, no. 5, 1515 pages, 2004.
- G.-L. Sun and W. Guo, “Bootstrapping M-estimators for reducing errors due to non-line-of-sight (NLOS) propagation,” IEEE Communications Letters, vol. 8, no. 8, 509 pages, 2004.
- R. B. Ertel and J. H. Reed, “Angle and time of arrival statistics for circular and elliptical scattering models,” IEEE Journal on Selected Areas in Communications, vol. 17, no. 11, 1829 pages, 1999.
- S. Venkatraman, J. Caffery, Jr., and H.-R. You, “A novel ToA location algorithm using LoS range estimation for NLoS environments,” IEEE Transactions on Vehicular Technology, vol. 53, no. 5, 1515 pages, 2004.