- Research
- Open Access
- Published:

# Adaptive mobile tracking in unknown non-line-of-sight conditions with application to digital TV networks

*EURASIP Journal on Advances in Signal Processing*
**volume 2014**, Article number: 22 (2014)

## Abstract

This paper studies the problem of tracking a mobile device in mixed line-of-sight (LOS) and non-line-of-sight (NLOS) conditions. NLOS error is assumed to be Gaussian with unknown mean and variance. An adaptive Rao-Blackwellized particle filter (RBPF) is proposed for mobile tracking in such scenarios. An extended Kalman filter is used to approximately estimate the mobile state, and the particle filter is applied to estimate the posterior distribution of sight conditions and the unknown static parameters, the distribution of which is updated by sufficient statistics. To improve the efficiency of the particle filtering, we use the approximate optimal proposal distribution for particle inference.

Algorithm performance is investigated in the scenario of mobile tracking using signals of opportunity from digital TV (DTV) network. Simulation results show that the adaptive RBPF method is effective to infer the unknown NLOS parameter and can achieve good tracking accuracy using a small number of particles.

## 1 Introduction

Accurate and reliable positioning in non-line-of-sight (NLOS) conditions is a challenging task in many wireless positioning systems. In typical NLOS circumstances, e.g. indoors and urban canyons, the direct path from the transmitter is blocked by buildings and other obstacles and the propagation wave may actually travel excess path lengths due to reflection, refraction and scattering. In terms of range-based measurements such as time of arrival (TOA), time difference of arrival (TDOA) and received signal strength (RSS), these extra propagation distances impose positive biases on the true path and thereby cause large errors on the location estimations. The NLOS effect is very common, especially in dense urban scenarios. A field test in a Global System for Mobile Communications (GSM) positioning system has shown that the mean and standard deviation of NLOS range errors are on the order of 513 and 436 m, respectively [1].

Methods proposed to deal with NLOS can generally be divided into two categories: methods for static positioning and methods for mobile tracking. Reference [2] gives a thorough summary for NLOS detection and mitigation in static positioning. However, these methods are not effective for mobile tracking. Methods to mitigate NLOS errors in mobile tracking applications generally exploit the measurements in time series. The proposed methods include two-step Kalman filtering techniques for smoothing range measurements and mitigating NLOS errors [3], a Kalman filter-based interacting multiple model (IMM) smoother [4], grid-based Bayesian estimation [5], particle filter (PF) [6], a modified extended Kalman filter (EKF) bank [7], the improved Rao-Blackwellized particle filter (RBPF) [8], and the joint particle filtering and unscented Kalman filtering (UKF) method [9]. A posterior Cramér-Rao lower bound for this problem is further investigated in [10]. A limitation of these methods is their assumption of complete knowledge of the model parameters, especially the statistical parameters of the NLOS errors, which is not realistic in many practical situations because of the unpredictable characteristics of the wireless channels in harsh environments.

By assuming all the measurement errors including the LOS and NLOS errors as unknown, reference [11] introduces the Dirichlet Process Mixtures (DPM), a hierarchical mixture model to describe the unknown density; a Bayesian nonparametric method is further applied for density estimation. Although it offers flexibility for modeling, effectively tuning the hyperparameters of a DPM model is difficult, and nonparametric Bayesian inference on the infinite dimensional Dirichlet process is obviously complicated.

In this study, we consider the mobile tracking problem under the mixed LOS/NLOS conditions with a plausible assumption that the statistical properties of the LOS errors are known while the statistics of NLOS error are unknown. In essence, as assumed in [3–9], the LOS errors are treated as measurement noise in the receiver with zero mean Gaussian distribution and known variance. By considering a long period of observation, the distribution of NLOS range errors can be modeled by a positively biased Gaussian noise, while to be more reasonable, its mean and variance are assumed (static) unknown because the statistical parameters of the NLOS errors are highly dependent on the surroundings, and therefore can not be known beforehand [12].

Mobile tracking under such mixed LOS/NLOS conditions comes down to the problem of sequential state estimation with the inference of unknown static parameters. To tackle this kind of problem, Liu and West assumed an artificial dynamic evolution for the unknown parameter vector, which could be further included in the state vector [13]. In the literature of GPS (Global Positioning System), similar methods are applied in algorithms adapting the density estimation by mean jump [14] or variance jump [15] or both [16]. However, such treatment would enlarge the estimation (co)variance of the unknown parameter. Djurić et al. [17] suggested the use of density-assisted particle filters (DAPFs) as an alternative to jointly estimate the sequential state and the parameter without introducing an artificial dynamic model for the static parameters. However, in our situation, because of the Markov property of the sight condition, the derivation of the density update function for all the state variables, and the static parameter is not an easy task. Storvik [18] proposed to marginalize the static parameters out of the posterior distribution such that only the state vector needs to be considered. The method is applicable when the conditional posterior distribution of the parameters (given the observations and the states) can be compactly expressed in terms of a finite dimensional statistic that can be recursively updated. Because the analytical form of the marginal distribution is not available, we hereby perform an approximate (Monte Carlo) marginalization by sampling from the joint posterior distribution of the states and the parameters. An adaptive particle filtering method is therefore proposed, which uses an EKF to approximately estimate the mobile state and applies the particle filter to estimate the posterior density of sight conditions and the unknown static parameters, the distribution of which is updated by sufficient statistics.

The performance of the proposed algorithm is investigated by simulations of positioning with signals of opportunity from digital television (DTV) networks. Recently, research interest in positioning using DTV systems has grown rapidly after the DTV systems have been put into operation for massive use [19–28]. It has been recognized that, compared with GPS [29], the DTV signals have a range of potential advantages to achieve low cost and accurate positioning: a higher transmission power [30], larger signal transmission bandwidth [31], less Doppler effects and ionosphere disturbance [32], lower carrier frequency with better diffraction and better receiving quality for urban and indoor propagation [19], wide coverage of DTV transmitting stations to provide a substantially superior geometry [25] and the economic benefit of making use of existing facilities [27]. Rosum Company has shown that the location accuracy could reach meter-scale with the Advanced Television Systems Committee (ATSC) DTV signals [19]. However, ATSC DTV is a single-carrier modulation system, which is vulnerable to multipath fading. In order to mitigate the effect of multipath fading for making the technique suitable for mobile application [31], digital broadcasting signals based on orthogonal frequency division multiplex (OFDM) modulation have been investigated and tested. As the standards based on OFDM modulation, e.g. DVB-T/H/T2, T-DMB, etc, have been widely adopted in most countries, wireless position systems based on the multi-carrier OFDM may have a massive number of potential users in the future. Therefore, in this paper, we will focus on mobile tracking in the multi-carrier OFDM DTV networks.

The paper is organized as follows: Section 2 presents the system model of mobile tracking in the mixed LOS/NLOS conditions; Section 3 formulates the problem within the Bayesian framework; in Section 4, the RBPF based adaptive mobile tracking method is described in detail; numerical results and performance comparison are presented and discussed in Section 5; and Section 6 draws some conclusions.

## 2 System description

### 2.1 Motion model

Assume the mobile device of interest moves in a plane. The state at time instant *t*_{
k
} is defined as the vector ${\mathbf{x}}_{k}={\left[{x}_{k},{y}_{k},{\stackrel{\u0307}{x}}_{k},{\stackrel{\u0307}{y}}_{k}\right]}^{T},$ where [*x*_{
k
},*y*_{
k
}]^{T} are the east and north coordinates of the mobile position; ${\left[{\stackrel{\u0307}{x}}_{k},{\stackrel{\u0307}{y}}_{k}\right]}^{T}$ are the corresponding velocities. The mobile state with random acceleration can be modeled as:

where the transition matrix ${\Phi}_{k}=\left(\begin{array}{cc}{\mathbf{I}}_{2}& \Delta {t}_{k}\phantom{\rule{0.3em}{0ex}}{\mathbf{I}}_{2}\\ 0& {\mathbf{I}}_{2}\end{array}\right)$ models the state kinematics with **I**_{2} the 2×2 matrix and *Δ* *t*_{
k
}=*t*_{k+1}−*t*_{
k
}. The random process **w**_{
k
} is a white zero mean Gaussian noise with covariance matrix ${\mathbf{Q}}_{k}=\left(\begin{array}{cc}\frac{\Delta {t}_{k}^{4}}{4}\mathbf{Q}& \frac{\Delta {t}_{k}^{3}}{2}\mathbf{Q}\\ \frac{\Delta {t}_{k}^{3}}{2}\mathbf{Q}& \Delta {t}_{k}^{2}\mathbf{Q}\end{array}\right)$, where $\mathbf{Q}=\left(\begin{array}{cc}{\sigma}_{x}^{2}& 0\\ 0& {\sigma}_{y}^{2}\end{array}\right)$. For a derivation of this motion model see [33].

### 2.2 Measurement model

In DVB-T/H networks, different emitters are coordinated to GPS time and transmit the same DTV signals at the same time with the same frequency, which is called the single frequency network (SFN) transmission [31]. The feature of network synchrony makes the pseudorange measurements available in practice for mobile tracking. Moreover, the DTV signals are transmitted continuously, which allows the receiver to track the arrival of the signals in order to improve the measurement accuracy.

For the purpose of communications, the DTV receivers are required to extract the timing measurements and recover the frequency offset from the received signals, which result in a synchronization problem. For reliable communications, the OFDM systems have stringent requirements on the timing and frequency synchronization. However, for the purpose of high accuracy positioning and navigation, an even finer synchronization is required. To achieve a finer time delay estimation, methods to determine the time of arrival (TOA) of the OFDM DTV signals include the coarse timing acquisition with a sliding correlator, which detects the start of an OFDM symbol by using the property of OFDM cyclic prefix (CP) [34], and the fine timing tracking with a delay-locked loop (DLL) [29], which achieves the accuracy within fractional portion of chip duration.

Scatter pilots inserted in the OFDM symbols [22, 23] or full transmitted OFDM symbols [27] can be used for DLL tracking, which are the variants to the process inside conventional Global Navigation Satellite System (GNSS) receivers. Time delay can also be obtained using the time-domain synchronous OFDM signals [25] or the transmitter signature waveforms [26]. In DVB-T/H systems, the nominal signal bandwidth is designed as 6 to 8 MHz, which leads to a much higher chip/code rate than that in GPS L1 systems, thereby improving the precision of timing.

Based on the above discussion, we consider the pseudorange measurements for mobile tracking from different DTV emitters. The pseudorange is the product of the estimated time delay with the speed of the light.

Under possible NLOS propagation condition between the transceivers, the distance measured at time *t*_{
k
} is

where ${d}_{i,k}\triangleq {h}_{i,k}\left({\mathbf{x}}_{k}\right)=\sqrt{{({x}_{k}-{x}_{{\text{bs}}_{i}})}^{2}+{({y}_{k}-{y}_{{\text{bs}}_{i}})}^{2}}$ represents the true distance between the receiver’s position [*x*_{
k
},*y*_{
k
}]^{T} and the location of the *i* th DTV transmitter ${\left[{x}_{{\text{Tx}}_{i}},{y}_{{\text{Tx}}_{i}}\right]}^{T},i\in \{1,2,\dots ,M\}$, and *M* is the number of DTV transmitters. The Boolean variable *s*_{i,k}∈{0,1} represents LOS/NLOS condition between the receiver and Tx_{
i
}, with *s*_{i,k}=0 for LOS and *s*_{i,k}=1 for NLOS. In mobile tracking, the sight conditions undergo dynamical transitions, which can be further modeled as a time-homogeneous first-order Markov chain *s*_{i,k}∼MC(*π*_{
i
},**A**_{
i
}) with initial probability vector *π*_{
i
} and the transition probability matrix

where *p*_{0}=P(*s*_{i,k}=0|*s*_{i,k−1}=0) and *p*_{1}=P(*s*_{i,k}=1|*s*_{i,k−1}=1).

Assume that the measurement noise in the LOS condition has a zero mean Gaussian distribution $\mathrm{N}(0,{\sigma}_{n}^{2})$, while the NLOS error is modeled as a biased Gaussian distribution $\mathrm{N}({\mu}_{\text{NLOS}},{\sigma}_{\text{NLOS}}^{2})$[3, 4, 7]. Thus, *v*(*s*_{i,k})∼N(*m*(*s*_{i,k}),**R**(*s*_{i,k})) and

In this work, we assume that *σ*_{
n
} is known, while $\{{\mu}_{\text{NLOS}},{\sigma}_{\text{NLOS}}^{2}\}$ are static but unknown.

To sum up, the overall dynamic model of the mobile tracking can be represented as follows:

where *s*_{i,k} is the *i* th component of vector **s**_{
k
}, *i*∈{1,2,…,*M*}, $k\in \mathbb{N}$.

## 3 Bayesian inference and sequential Monte Carlo method

Denote the total observation sequence up to time *t*_{
k
} as **z**_{1:k}, where ${\mathbf{z}}_{k}\triangleq \phantom{\rule{0.3em}{0ex}}{[\phantom{\rule{0.3em}{0ex}}{z}_{1,k},{z}_{2,k},\dots ,{z}_{M,k}]}^{T}$. For brevity, let $\eta \triangleq {\sigma}_{n}^{2}+{\sigma}_{\text{NLOS}}^{2}$ and ** θ**={

*μ*

_{NLOS},

*η*}. The problem of mobile tracking in unknown NLOS conditions is to infer the current mobile state

**x**

_{ k }from the observation sequence

**z**

_{1:k}. Within the framework of Bayesian inference, the problem corresponds to computing the marginal posterior probability density function (pdf)

*p*(

**x**

_{ t }|

**z**

_{1:k}). The marginal posterior is the mixture

where the prior pdf *p*(**x**_{
k
},**s**_{
k
},** θ**|

**z**

_{1:k−1}) can be obtained via the Chapman-Kolmogorov equation:

In (5) and (6), the number of mixture components grows exponentially with time. Thus, the analytical solution to the posterior *p*(**x**_{
k
}|**z**_{1:k}) requires very high-dimensional integrals, which is prohibitive to compute in practice. The sequential Monte Carlo (SMC) method (also called particle filtering) has proven to be successful in tracking applications with nonlinear and non-Gaussian models [35]. Here, we resort to this kind of sample-based numerical approximate method.

Denote **y**_{
k
}={**x**_{
k
},**s**_{
k
}}, and suppose a set of *N* weighted samples ${\{{\mathbf{y}}_{k-1}^{j},{\mathit{\theta}}^{j},{w}_{k-1}^{j}\}}_{j=1}^{N}$ is used to approximate the posterior *p*(**y**_{k−1},** θ**|

**z**

_{1:k−1}) at time

*t*

_{k−1}with the following point-distribution:

where *δ*(·) is the Dirac delta measure, and the multiplication of Dirac deltas is the product measure of the one-dimensional Delta functions in each variable separately [36].

With the new reception of measurement **z**_{
k
}, the new samples at time *t*_{
k
} are generated from a suitably designed proposal distribution:

Accordingly, importance weights are updated as

In standard particle filtering (SPF), transition priors are used as the proposal distribution:

Thus, the weight update equation (7) can be simplified as:

In SPF, since {**x**_{
k
},**s**_{
k
},** θ**} constitutes a high-dimensional state estimation space, a large number of particles should be used to achieve good estimation results, which is computationally expensive. Additionally, using the transition prior as the proposal, which fails to consider the information of current measurements, would easily suffer from the ‘particle impoverishment’ problem. In order to overcome these two limitations, we re-formulate the inference problem using a RBPF.

## 4 Adaptive RBPF algorithm for mobile tracking

Factorize the posterior density of the hidden state *p*(**x**_{
k
},**s**_{
k
},** θ**|

**z**

_{1:k}) according to the chain rule:

If the marginal posterior density *p*(**s**_{
k
},** θ**|

**z**

_{1:k}) is represented by a set of weighted samples ${\{{\mathbf{s}}_{k}^{j},{\mathit{\theta}}^{j},{w}_{k}^{j}\}}_{j=1}^{N}$, i.e.,

then the marginal density *p*(**x**_{
k
}|**z**_{1:k}) can be approximately expressed by a mixture of densities:

Based on the system model in (4), the mixture component $p\left({\mathbf{x}}_{k}\right|{\mathbf{s}}_{k}^{j},{\mathit{\theta}}^{j},{\mathbf{z}}_{1:k})$ can be calculated by an EKF.

Thus, the decomposition (8) leads to developing a more efficient algorithm, with only the posterior density *p*(**s**_{
k
},** θ**|

**z**

_{1:k}) approximately calculated through a sampling method, while the mobile state

*p*(

**x**

_{ k }|

**s**

_{ k },

*θ*_{1:k}) is analytically computed. This method, motivated by the decomposition (8), is known as RBPF. In what follows, the method is illustrated in detail.

### 4.1 Inference for the mobile state

Let a set of weighted samples ${\{{\mathbf{s}}_{k}^{j},{\mathit{\theta}}^{j},{w}_{k}^{j}\}}_{j=1}^{N}$ represent the marginal posterior density of *p*(**s**_{
k
},** θ**|

**z**

_{1:k}). The sampling method and parameter update will be described in Sections 4.2 and 4.3. Then, the marginal posterior

*p*(

**x**

_{ k }|

**z**

_{1:k}) in (10) is approximated by a mixture of Gaussians, with the component $p\left({\mathbf{x}}_{k}\right|{\mathbf{s}}_{k}^{j},{\mathit{\theta}}^{j},{\mathbf{z}}_{1:k})$ approximately conforming to Gaussian distribution $\mathrm{N}({\widehat{\mathbf{x}}}_{k}^{j},{\widehat{\mathbf{P}}}_{k}^{j})$, where

${\widehat{\mathbf{x}}}_{k|k-1}^{j}$ is the predicted mean of ${\mathbf{x}}_{k-1}^{j}$:

and ${\widehat{\mathbf{P}}}_{k|k-1}^{j}$ is the corresponding predicted covariance:

The predicted mean of measurement ${\widehat{z}}_{i,k|k-1}^{j}$ is

The Kalman gain is

and ${\mathbf{H}}_{i,k}^{j}={\frac{\partial {h}_{i}}{\partial \mathbf{x}}|}_{\mathbf{x}={\widehat{\mathbf{x}}}_{k|k-1}^{j}}$. In LOS conditions, $m\left({s}_{i,k}^{j}\right)=0$ and $\mathbf{R}\left({s}_{i,k}^{j}\right)={\sigma}_{n}^{2}$, while in NLOS, $m\left({s}_{i,k}^{j}\right)={\mu}_{\text{NLOS}}^{j}$ and $\mathbf{R}\left({s}_{i,k}^{j}\right)={\eta}^{j}$.

### 4.2 Particle sampling and weights updating

The proposal trial distribution is said to be optimal if it minimizes the variance of the important weights [37]. To sample $\{{\mathbf{s}}_{k}^{j},{\mathit{\theta}}^{j}\}$, conditioned upon ${\mathbf{s}}_{k-1}^{j},{\mathbf{x}}_{k-1}^{j},{\mathit{\theta}}^{j},$ and **z**_{1:k}, the optimal proposal distribution is:

Accordingly, the weight is updated as:

Obviously, the importance weight update equation in (18) involves a high-dimension integral and it is difficult to get the closed-form solution. To approximate the integral in (18), we use the mean-value point approximation

where

where **E**(·) is the expectation operator.

Then, substituting (19) and (20a) into (18), the importance weight update is approximated by:

Applying the same mean-point approximation (20a) into (17a), the optimal proposal distribution can be approximated as

### 4.3 Update of the static parameters

In this work, the NLOS error is modeled as a Gaussian random variable with positive mean. To infer the statistical parameter ** θ**, we specify on them the Gaussian inverse chi-square prior, which is conjugate prior distribution and has computational convenience in Bayesian inference [38].

Suppose at time *t*_{k−1},

that is,

where the four hyperparameters $\{{\stackrel{\u0306}{\mu}}_{k-1}^{j},{\stackrel{\u0306}{\kappa}}_{k-1}^{j},{\stackrel{\u0306}{\nu}}_{k-1}^{j},{\stackrel{\u0306}{\eta}}_{k-1}^{j}\}$ can be identified as the location and scale of *μ*_{NLOS} and the degrees of freedom and scale of *η*.

At time *t*_{
k
}, the sampling density for ** θ** is updated according to

where the hyperparameters are updated as

where $\{{n}_{k}^{j},{\epsilon}_{k}^{j},{\stackrel{\u0304}{\epsilon}}_{k}^{j}\}$ are computed as

The ${n}_{k}^{j}$ is the total number of NLOS conditions in ${\mathbf{s}}_{k}^{j}$, and ${\stackrel{\u0304}{\epsilon}}_{k}^{j}$ is the mean of the innovation in NLOS conditions. Intuitively, the updated hyperparameters in (25) combine the prior information and the information contained in $\{{n}_{k}^{j},{\stackrel{\u0304}{\epsilon}}_{k}^{j}\}$. Obviously, for ${n}_{k}^{j}=0$, which means there is no NLOS condition inferred in ${\mathbf{s}}_{k}^{j}$, the hyperparameters remain unchanged from time *t*_{k−1} to *t*_{
k
} because there is no innovation to update. Thus, $\{{n}_{k}^{j},{\stackrel{\u0304}{\epsilon}}_{k}^{j}\}$ derived from $\{{\mathbf{s}}_{k}^{j},{\mathbf{x}}_{k}^{j},{\mathbf{z}}_{k}\}$ are the sufficient statistics to infer the parameter ** θ**.

To get the samples *θ*^{j} from the proposal distribution (22a), first sample ${\mathbf{s}}_{k}^{j}\sim \mathrm{P}\left({\mathbf{s}}_{k}\right|{\mathbf{s}}_{k-1}^{j})$. Then update the hyperparameters according to (25-26), and generate the particles *θ*^{j} from the Gaussian-Inverse-chi-square, which corresponds to ${\eta}^{j}\sim {\chi}^{-2}({\stackrel{\u0306}{\nu}}_{k}^{j},{\stackrel{\u0306}{\eta}}_{k}^{j})$ and ${\mu}_{\text{NLOS}}^{j}\sim \mathrm{N}({\stackrel{\u0306}{\mu}}_{k}^{j},{\eta}^{j}/{\stackrel{\u0306}{\kappa}}_{k}^{j})$. To get the weights at each time epoch in (21), the mean of *θ*^{j} is computed as:

### 4.4 Algorithm description

The importance weight ${w}_{k}^{j}$ in (18) only depends on the current measurement **z**_{
k
} and the particles of *t*_{k−1}, i.e. ${\{{\widehat{\mathbf{x}}}_{k-1}^{j},{\mathbf{s}}_{k-1}^{j}\}}_{j=1}^{N}$, while {**x**_{
k
},**s**_{
k
},** θ**} are all marginalized out. Thus, to improve the sample effectiveness, the particles of

*t*

_{k−1}could be first selected (resampled) based on the current measurement

**z**

_{ k }, and the fittest particles could be allowed to propagate. Then, ${\{{\mathit{\theta}}^{j},{\mathbf{s}}_{k}^{j}\}}_{j=1}^{N}$ is sampled from the (approximately) optimal distribution, with an EKF attached to every particle to further update the mobile state ${\mathbf{x}}_{k}^{j}$. Thus, the whole algorithm infers the sight conditions, the mobile state and adaptively learns the unknown statistical parameters of the biased NLOS errors. The proposed method is summarized in Algorithm 1.

#### Algorithm 1 **RBPF-based adaptive mobile tracking algorithm**

## 5 Simulation results

Computer simulations are performed in a DVB-T-based positioning system to evaluate the proposed adaptive mobile tracking algorithm in the mixed LOS/NLOS conditions.

It is assumed that at every epoch, the signals can be received from five transmitters. By considering the SFN coverage in 4K mode of DVB-T/H, which supports the application of the mobile reception [31], the DTV transmitters are set at [ −2 km, −1 km], [ −2 km, 6 km], [ 5 km, −1 km], [ 6 km, 5 km] and [ 1 km, −2 km], respectively. The receiver’s trajectories are generated according to the motion model described in Section 2.1, in which the initial position is set to [ −1.5 km, 1.5 km], and the initial velocity is set at [ 10 m/s,0 m/s]. The random acceleration variances ${\sigma}_{x}^{2},\phantom{\rule{0.3em}{0ex}}{\sigma}_{y}^{2}$ are both chosen to 0.5 (m/s^{2})^{2}. The simulated trajectory has *L*=1,000 time steps, and the time step size is *Δ* *t*=0.2 s. The measurement data are generated by adding the measurement noise and the NLOS noise to the true distance from receiver to each transmitter. The measurement noise is assumed to be a white Gaussian random variable with zero mean and standard deviation *σ*_{
n
}=15 m, which is in agreement with the theoretical analysis of the 4K mode in the DVB-T system in [23] and the simulation results in [27]. The NLOS measurement errors are usually much larger than LOS errors [1]. Thus, the parameters of the positive Gaussian distribution of the NLOS errors are assumed with mean *μ*_{NLOS}=50 m and standard deviation *σ*_{NLOS}=40 m. The mode transition probability is chosen to be *p*_{0}=*p*_{1}=0.8. The LOS or NLOS mode is generated by making a Markov chain transition every 10 steps, which simulates a highly dynamic environment during the tracking process.

In the algorithms, the initial estimation of sight condition is set to P(*s*_{i,0}=0)=P(*s*_{i,0}=1)=0.5, where *i*=1,⋯,5. The initial position is calculated by Chan’s algorithm [39] using the first five range measurements. For the lack of definitive prior information on the mobile state, the initial velocity is set as [ 0 m/s,0 m/s] and the covariance matrix ${\mathbf{C}}_{t0}=\left[\begin{array}{cc}1{5}^{2}\xb7{\mathbf{I}}_{2}& \mathbf{0}\\ \mathbf{0}& 1{0}^{2}\xb7{\mathbf{I}}_{2}\end{array}\right]$ corresponding to a standard deviation of 15 m for the position and 10 m/s (36 km/h) for the velocity of each coordinate. For vague prior information on NLOS parameter *θ*, the initial values of the hyperparameters are set as $\{{\stackrel{\u0306}{\mu}}_{0}=1,000,{\stackrel{\u0306}{\kappa}}_{0}=1,{\stackrel{\u0306}{\nu}}_{0}=1,{\stackrel{\u0306}{\eta}}_{0}={\left(5{\sigma}_{n}\right)}^{2}\}$ in the algorithm. The simulation results are obtained based on *n*_{MC}=20 Monte Carlo realizations with the same parameters.

We compare the performance of the adaptive RBPF method with the following three reference methods. The first is the RBPF method [8], where the NLOS parameters ** θ** is assumed known and only the mobile state

**x**

_{ k }and the sight condition

**s**

_{ k }have to be inferred. The second reference method assumes the sight conditions known for the whole trajectory and the adaptive RBPF is modified only to infer the mobile state

**x**

_{ k }and the static parameters

**(adaptive RBPF with**

*θ***s**

_{ k }known). The aim of the comparisons with these two methods is to show how the different parameters affect the final inference for the mobile tracking. In addition, to show the superiority of the RBPF method over the SPF method mentioned in Section 3, we also compare the results obtained by the SPF, which uses the transition prior as the proposal to sample the high dimensional state space {

**x**

_{ k },

**s**

_{ k },

**}, and applies the same parameter update method for**

*θ***in Section 4.3. In what follows, the third reference method is called adaptive SPF.**

*θ*The accuracies are compared in terms of root square error (RSE), position root-mean-square error (RMSE) and average RMSE. RSE is defined as $\text{RSE}\triangleq \sqrt{{({\widehat{x}}_{k}-{x}_{k})}^{2}+{({\u0177}_{k}-{y}_{k})}^{2}}$ and the position RMSE at time *t*_{
k
} as: ${\text{RMSE}}_{k}\triangleq \sqrt{\frac{1}{{n}_{\text{mc}}}\sum _{m=1}^{{n}_{\text{mc}}}[{({\widehat{x}}_{k,m}-{x}_{k})}^{2}+{({\u0177}_{k,m}-{y}_{k})}^{2}]}$, while $\text{average RMSE}\triangleq \frac{1}{L}{\sum}_{k=1}^{L}{\text{RMSE}}_{k}$.

### 5.1 Tracking accuracy of the proposed algorithm

In the algorithms of RBPF, adaptive RBPF with **s**_{
k
} known and adaptive RBPF, the number of the particles is set to 10. To achieve a comparable accuracy, 1,000 particles are used in the adaptive SPF. Figure 1 shows comparison of the CDF of RSE, in which the positioning accuracy of 66.7% and 99% errors are also marked and the comparison of the position RMSE versus time is in Figure 2. Figures 3 and 4 show one realization of updating the static parameters *μ*_{NLOS} and $\sqrt{\eta}$; 3- *σ* confidence intervals are also shown. In all these figures, the numbers in the brackets denote the number of particles used.

From Figures 1 and 2, the adaptive RBPF with **s**_{
k
} known achieves the smallest positioning errors with the 67% error 6.0 m and 95% error 10.0 m, while RBPF and adaptive RBPF has approximately the same accuracy. Although using 1000 particles, the adaptive SPF still has the largest error statistics. In Figures 3 and 4, the true values of the static parameters are within the 3- *σ* confidence interval, which suggests a good inference of the unknown parameters of all three algorithms. By further comparisons, the sequential estimations by the adaptive RBPF with **s**_{
k
} known are slightly better than the adaptive RBPF. Both of the above two algorithms achieve better estimations than the adaptive SPF method, which is obvious in Figure 3.

From Figures 1, 2, 3 and 4, among all the algorithms, the adaptive RBPF with **s**_{
k
} known achieves the best accuracy. The reason is that, in the proposed algorithm, the sufficient statistics for updating ** θ** and the mobile state inference for

**x**

_{ k }are largely dependent on the density estimation of

**s**

_{ k }. When the sight conditions are known during the whole trajectory, the algorithm could have more accurate estimation on NLOS parameter

*θ*, which further improves the estimation for the mobile state. Therefore, to improve the estimation accuracy, the accurate inference on the sight condition has the most importance. RBPF achieves slightly better performance over adaptive RBPF, which is reasonable, since in RBPF, the NLOS parameter is known and the parameter updating is not included. But the improvement is slight, as shown in Figures 1 and 2. Combined with the results of Figures 3 and 4, it is clear that the proposed adaptive RBPF can effectively estimate the unknown mean and variance. Even with 100 times more particles and having the most computation complexity, the adaptive PF still has the worst accuracy among all the algorithms, which suggests that using the prior transition as the proposal distribution is not effective to get the fittest particles in mobile tracking and the unknown parameter inference.

### 5.2 Complexity comparison

Table 1 compares the relative complexity and accuracy of the algorithms. It is clear that with the increase of the particle numbers, the computing time of the adaptive RBPF increases proportionally. Accordingly, the accuracy also increases. But the improvement is slight when the number is larger than 50. Also, the accuracy improvement can be omitted when considering the large increase of complexity between 100 and 1,000 particles. Thus, the adaptive RBPF (10) achieves a good tradeoff between complexity and accuracy.

In this section, we compared the performance of adaptive RBPF with the other three algorithms, i.e. the adaptive RBPF with **s**_{
k
} known, the RBPF which assumes the static parameters are known and the adaptive SPF method. Simulation results show that the accurate estimation of sight condition has an important effect on the ultimate accuracy of mobile tracking and parameter estimation. By applying the approximate optimal proposal distribution to sample the posterior distribution of the sight conditions, the adaptive RBPF method is effective to infer the unknown NLOS parameter and achieves a good tracking accuracy with small number of particles.

## 6 Conclusions

We have considered the problem of mobile tracking in the mixed LOS/NLOS conditions, where the statistical parameter of NLOS error is unknown. Under the sequential Monte Carlo Bayesian framework, an adaptive RBPF method is developed, which uses an analytical method to estimate the mobile state while applying the particle filter to estimate the posterior density of sight conditions and the unknown static parameters. The distribution of the static parameters is updated by sufficient statistics from the mobile state, the sight condition and the measurements at current step. To improve the efficiency of the particle filtering, the approximate optimal proposal distribution is used for particle inference.

Simulation tests of positioning with signals from DTV networks show that the adaptive RBPF method is effective to infer the mobile state and the unknown NLOS parameters simultaneously. With only 10 particles, the adaptive RBPF achieves a good accuracy of the mobile tracking while the complexity has not been increased much.

To achieve more accurate and robust mobile tracking, future work will investigate the impact of NLOS in the case of non-Gaussian distribution. We will also compare the unscented Kalman filter or the cubature Kalman filter with the current EKF method in the adaptive particle filter scheme.

## References

- 1.
Silventoinen MI, Rantalainen T: Mobile Station Emergency Locating in GSM. In

*Proceedings of IEEE International Conference on Personal Wireless Communications*. New Delhi; 19–21 Feb 1996:232-238. - 2.
Guvenc I, Chong C-C: A survey on TOA based wireless localization and NLOS mitigation techniques.

*IEEE Commun. Surv. Tutorials*2009, 11(3):107-124. - 3.
Le BL, Ahmed K, Tsuji H: Mobile location estimator with NLOS mitigation using Kalman filtering.

*IEEE Wireless Communications and Networking Conference, vol. 3. New Orleans, USA*18–20 March 2003, 1969-1973. - 4.
Liao J-F, Chen B-S: Robust mobile location estimator with NLOS mitigation using interacting multiple model algorithm.

*IEEE Trans. Wireless Commun*2006, 5(11):3002-3006. - 5.
Morelli C, Nicoli M, Rampa V, Spagnolini U: Hidden Markov models for radio localization in mixed LOS/NLOS, conditions.

*IEEE Trans. Signal Process*2007, 55(4):1525-1542. - 6.
Nicoli M, Morelli C, Rampa V: A jump Markov particle filter for localization of moving terminals in multipath indoor scenarios.

*IEEE Trans. Signal Process*2008, 56(8):3801-3809. - 7.
Chen L, Wu L: Mobile positioning in mixed LOS/NLOS conditions using modified EKF banks and data fusion method.

*IEICE Trans. Commun*2009, EB92(4):1318-1325. - 8.
Chen L, Ali-Löytty S, Piché R, Wu L: Mobile tracking in mixed line-of-sight/non-line-of-sight conditions: algorithm and theoretical lower bound.

*Wireless Personal Communications*2012, 65(4):753-771. 10.1007/s11277-011-0294-7 - 9.
Huerta JM, Vidal J, Giremus A, Tourneret J-Y: Joint particle filter and UKF position tracking in severe non-line-of-sight situations.

*IEEE J. Selected Topics Signal Process*2009, 3(5):874-888. - 10.
Chen L, Wu L, Piché R: Posterior Cramer-Rao bounds for mobile tracking in mixed LOS/NLOS conditions. In

*Proceedings of 17th European Signal Processing Conference*. Glasgow, Scotland; 2009:90-94. - 11.
Rabaoui A, Viandiery N, Maraisy J, Duflos E: DPMs for the density estimation in a dynamic nonlinear modeling Application to GPS positioning in urban canyons.

*IEEE Trans. Signal Process*2012, 60(4):1638-1655. - 12.
Chen L, Piché R: Mobile tracking and parameter learning in unknown non-line-of-sight conditions. In

*the Preceedings of the 13th Conference on Information Fusion (FUSION)*. Edinburgh, Scotland; 26–29 July 2010:1-6. - 13.
Liu J, West M: Combined parameter and state estimation in simulation-based filtering. In

*Sequential Monte Carlo Methods in Practice*. Edited by: Doucet A, Freitas JFG, Gordon NJ. New York: Springer; 2001:197-223. - 14.
Giremus A, Tourneret J-Y, Calmettes V: A particle filter approach for joint detection/estimation of multipath effects on GPS measurements.

*IEEE Trans. Signal Process*2007, 55(4):1275-1285. - 15.
Spangenberg M:

*Safe navigation for vehicles*. PhD thesis, Toulouse University; 2009. - 16.
Caron F, Davy M, Duflos E, Vanheeghe P: Particle filtering for multisensor data fusion with switching observation models: Application to land vehicle positioning.

*IEEE Trans. Signal Process*2007, 55(6):2703-2719. - 17.
Djurić PM, Bugallo MF, Míguez J: Density assisted particle filters for state and parameter estimation. In

*IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP ‘04)*. Mon0treal QC, Canada; 17–21 May 2004:701-704. - 18.
Storvik G: Particle filters for state-space models with the presence of unknown static parameters.

*IEEE Trans. Signal Process*2002, 50(2):281-289. 10.1109/78.978383 - 19.
Rabinowitz M, Spilker J: Positioning using the ATSC digital television signal. U.S. Patent 6,861,984 B2, 1 March 2005

- 20.
Wang X, Wu Y, Caron B: Transmitter identification using embedded pseudo random sequences.

*IEEE Trans. Broadcasting*2004, 50(3):244-252. 10.1109/TBC.2004.834027 - 21.
Wang X, Wu Y, Chouinard J: A new position location system using DTV transmitter identification watermark signals.

*EURASIP J. Appl. Signal Process*2006, 2006: 1-11. - 22.
Serant D, Thevenon P, Boucheret ML, Julien O, Macabiau C, Corazza S, Dervin M, Ries L: Development and Validation of an OFDM/DVB-T Sensor for Positioning. In

*Proceedings of the IEEE/ION Position Location and Navigation Symposium (PLANS10)*. IndianWells, Calif, USA; 4–6 May 2010:998-1001. - 23.
Serant D, Julien O, Macabiau C, Ries L, Thevenon P, Dervin M, Boucheret ML: Positioning Using OFDM-based Digital TV: New Algorithms and Tests with Real Signals. In

*Proceedings of the 24th International Technical Meeting of the Satellite Division of The Institute of Navigation*. Portland, OR, USA; 19–23 September 2011:3451-3460. - 24.
Thevenon P, Julien O, Macabiau C, Serant D, Ries L, Corazza S, Boucheret ML: Positioning principles with a mobile TV system using DVB-SH signals and a single frequency network. In

*Proceedings of the Digital Signal Processing, (DSP09)*. Santorini, Greece; 5–7 July 2009. - 25.
Dai L, Wang Z, Pan C, Chen S: Wireless positioning using TDS-OFDM signals in single-frequency networks.

*IEEE Trans. Broadcasting*2012, 58(2):236-246. - 26.
Yang J, Wang X, Rahman MJ, Park SI, Kim H-M, Wu Y: A new positioning system using DVB-T2 transmitter signature waveforms in single frequency networks.

*IEEE Trans. Broadcasting*2012, 58(3):347-359. - 27.
Chen L, Yang L-L, Chen R: Time delay tracking for positioning in DTV networks.

*Proceedings of the Ubiquitous Positioning, Indoor Navigation, and Location Based Service (UPINLBS), 2012*3–4 October 2012, 1-4. - 28.
Yan J, Wu L: A passive location system for single frequency networks using digital terrestrial TV signals.

*Eur. Trans. Telecommunications*2011, 22(8):487-499. 10.1002/ett.1498 - 29.
Kaplan ED(ed.),

*Understanding GPS: Principles and Applications*, (Artech House, Norwood, 1996) - 30.
Rabinowitz M, Spilker Jr JJ: A new positioning system using television synchronization signals.

*IEEE Trans. Broadcasting*2005, 51(1):51-61. 10.1109/TBC.2004.837876 - 31.
ETSI Standard, ETSI EN 302 744 V1.1.4 (2011)

- 32.
Ong JT, Yan H, Rao SV, Shanmugam G: Indoor DTV reception: measurement techniques.

*IEEE Trans. Broadcasting*2004, 50: 192-199. 10.1109/TBC.2004.828371 - 33.
Bar-Shalom Y, Li RX, Kirubarajan T:

*Estimation with Applications to Tracking and Navigation, Theory Algorithms and Software*. New York: John Wiley & Sons; 2001. - 34.
van de Beek J-J, Sandell M, Börjesson PO: ML estimation of time and frequency offset in OFDM systems.

*IEEE Trans. Signal Process*1997, 45: 1800-1805. 10.1109/78.599949 - 35.
Ristic B, Arulampalam S, Gordon N:

*Beyond the Kalman Filter, Particle Filters for Tracking Applications*. Artech House, Boston, London; 2004. - 36.
Haörmander L:

*The Analysis of, Linear Partial Differential Operators I, 2nd edn. Classics in Mathematics*. Berlin, Heidelberg, New York: Springer; 1990. - 37.
Arulampalam MS, Maskell S, Gordon N, Clapp T: A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking.

*IEEE Trans. Signal Process*2002, 50(2):174-188. 10.1109/78.978374 - 38.
Gelman AB, Carlin JS, Stern HS, Dubin DB:

*Bayesian Data Analysis*. London: Chapman & Hall; 2000. - 39.
Chan YT, Ho KC: A simple and efficient estimator for hyperbolic location.

*IEEE Trans. Signal Process*1994, 42(8):1905-1915. 10.1109/78.301830

## Acknowledgements

The authors wish to thank Professor Joaquín Míguez from the Department of Signal Theory & Communications in Universidad Carlos III de Madrid, Spain for the fruitful discussions and valuable comments. This research has been supported by the Academy of Finland Postdoctoral Research Funding (No.254232).

## Author information

### Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Chen, L., Piché, R., Kuusniemi, H. *et al.* Adaptive mobile tracking in unknown non-line-of-sight conditions with application to digital TV networks.
*EURASIP J. Adv. Signal Process.* **2014, **22 (2014). https://doi.org/10.1186/1687-6180-2014-22

Received:

Accepted:

Published:

### Keywords

- Parameter estimation
- Bayesian inference
- Signal of opportunity
- Digital TV
- Non-line-of-sight
- Rao-Blackwellized particle filter
- Wireless positioning
- Mobility tracking