Next Article in Journal / Special Issue
Vehicular Visible Light Communication for Intersection Management
Previous Article in Journal
Breast Density Transformations Using CycleGANs for Revealing Undetected Findings in Mammograms
Previous Article in Special Issue
A Survey on Optimal Channel Estimation Methods for RIS-Aided Communication Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extended Kalman Filter Design for Tracking Time-of-Flight and Clock Offsets in a Two-Way Ranging System

Center for Wireless Information Systems and Computational Architectures (WISCA), Arizona State University, Tempe, AZ 85281, USA
*
Author to whom correspondence should be addressed.
Signals 2023, 4(2), 439-456; https://doi.org/10.3390/signals4020023
Submission received: 30 December 2022 / Revised: 13 February 2023 / Accepted: 16 May 2023 / Published: 15 June 2023
(This article belongs to the Special Issue Advances in Wireless Sensor Network Signal Processing)

Abstract

:
As radio frequency (RF) hardware continues to improve, two-way ranging (TWR) has become a viable approach for high-precision ranging applications. The precision of a TWR system is fundamentally limited by estimates of the time offset T between two platforms and the time delay τ of a signal propagating between them. In previous work, we derived a family of optimal “one-shot” joint delay–offset estimators and demonstrated that they reduce to a system of linear equations under reasonable assumptions. These estimators are simple and computationally efficient but are also susceptible to channel impairments that obstruct one or more measurements. In this work, we formulate an extended Kalman filter (EKF) for this class of estimators that specifically addresses this limitation. Unlike a generic KF approach, the proposed solution specifically integrates the estimation process to minimize the computational complexity. We benchmark the proposed first- and second-order EKF solutions against the existing one-shot estimators in a MATLAB Monte Carlo simulation environment. We demonstrate that the proposed solution achieves comparable estimation performance and, in the case of the second-order solution, reduces the computation time by an order of magnitude.

1. Introduction

As flying, ad-hoc networks (FANETs) become increasingly accessible and affordable, they demand increasing access to the already-crowded frequency spectrum. Intelligent transport systems (ITSs) and internet-of-things (IoT) networks require numerous measurement modalities and are increasingly dissatisfied with the limited reliability and efficiency of legacy wireless sensor network (WSN) solutions. Modern solutions have transitioned towards cooperative and co-design techniques under the category of “spectral convergence” that promise significant performance enhancements and massively increased spectral capacity [1,2,3].
In a previous work [4], we demonstrated a novel two-way ranging (TWR) system that integrated these modern co-design techniques. Specifically motivated by the limitations of GPS and ZigBee, this system efficiently and simultaneously enabled localization, synchronization, and communication services for wireless unmanned aerial vehicle (UAV) networks. We experimentally demonstrated high-precision ranging (sub-cm) and synchronization (sub-ns) with limited spectral access (10 MHz) and transmit power (1 W) [5] and further addressed the growing vulnerabilities of GPS [6,7,8] using an integrated encryption and authentication protocol.
While these initial results were successful as a “proof of concept”, this early prototype employed a simplified processing chain that is susceptible to outages, cycle slips, shadowing, multipath fading, and numerous other realistic channel conditions. Since FANETs typically require more robust and reliable positioning services, this prototype was largely considered immature and insufficient for safety-critical applications. To address these limitations, we have developed a family of novel extended Kalman filter (EKF) extensions that enhance and make robust this system. These tracking filters improve the positioning performance by assuming more realistic first- and second-order Markov models and reduce the susceptibility to outages and cycle slips by estimating the parameters over multiple frames rather than the previous “one-shot” approach. We specifically evaluate the proposed methods against our previous results [9] using a simple MATLAB simulation platform and demonstrate notable improvements in performance and processing speed.

1.1. System Overview

The proposed two-way ranging system measures the time-of-flight (ToF) between multi-antenna radio platforms, as depicted in Figure 1. The radio platforms operate with independent clocks; so, the measured time-of-arrival (ToA) is affected by both the propagation time τ and the clock offset T between the platforms. Using each antenna on the platforms, we measure the ToF between each transmit–receive antenna pair, creating a spatial diversity that can then be used to estimate the relative position and orientation of the platforms. We can construct and solve (in closed form) a system of equations using τ and T and their first derivatives. This closed-form solution is attractive for its computational simplicity and tractability but is sensitive to environmental conditions such as shadowing and fading that might disrupt one or more individual measurements. To mitigate these issues, we propose a family of EKF solutions that achieve comparable or better performance without increasing the computation time; we even demonstrate that the proposed second-order EKF solution reduces the computation time by an order of magnitude compared to its one-shot equivalent.

1.2. Contributions

In this manuscript, we make the following contributions:
  • We derive a family of optimal extended Kalman filter (EKF) tracking algorithms that jointly estimate time-of-flight and clock offsets in a two-way ranging network.
  • We implement the proposed EKF solutions and the corresponding optimal one-shot estimators in a simple MATLAB simulation environment.
  • We demonstrate that the proposed solution achieves comparable estimation performance to the existing one-shot solutions and, in the case of the second-order solution, reduces the computation time by an order of magnitude.

1.3. Organization

The remainder of this manuscript is organized as follows: in Section 2, we briefly review the relevant literature to contextualize our problem and introduce the previous results upon which our proposed solution is built; in Section 3, we specify the problem statement and provide the necessary preliminary definitions; in Section 4, we derive the optimal one-shot estimators for this problem assuming first- and second-order Markov models; in Section 5, we formulate extended Kalman filter tracking algorithms assuming these same models; in Section 6, we benchmark these estimators and algorithms in a simple MATLAB simulation platform; and in Section 7, we provide concluding remarks.

2. Background

Robust and precise localization is a critical service for adaptive self-organizing ad hoc networks. Extensive reviews on positioning using the angle-of-arrival (AoA), time-of-arrival (ToA), time-difference-of-arrival (TDoA), and received-signal-strength (RSS) are available in [10,11,12,13,14,15]. ToA positioning methods [16] are increasingly popular and have received growing interest in the last decade. Multi-antenna systems may leverage ToA measurements to estimate the position and orientation, the precision of which is dictated by the quality of the original measurements [17,18,19]. For broadcast systems such as the Global Positioning System (GPS) that have extremely high-quality clocks [20,21], it is easy to ignore the timing aspect of this problem; however, for consumer-grade solutions, tracking the relative clock offsets is absolutely critical.

2.1. Clock Synchronization

Since TWR-derived positioning relies heavily on precise clock synchronization [22,23,24,25,26], distributed coherence has become a common requisite for these types of networks. Coarse approaches such as the Network Timing Protocol (NTP) are sufficient for communications applications [27,28,29] but do not support the precision we need for rapid localization. Several extensions have improved this performance [30,31,32,33,34,35,36], but many have imposed simplifying assumptions, are not computationally efficient, and lack generality to changes in the underlying models. In response, we proposed a family of novel one-shot estimators in [9] that specifically integrate clock synchronization in the localization algorithm. We demonstrated that these estimators were optimal for single measurements and reduced to a system of linear equations. In this work, we formulate an extended Kalman filter (EKF) extension [37,38,39] to these estimators to further reduce the computational complexity and improve the robustness to cycle slips and other outages.

2.2. Statistical Clock Models

Estimating clock errors involves solving stochastic differential equations whose coefficients are derived from the Allan variance [40] characteristics of the oscillators that drive them [41]. EKF-based clock synchronization is well-studied in this context [42,43]. We draw insights from [44,45] to actualize the joint distributed coherence and ToF tracking using an extended Kalman filter. We do not consider flicker noise while modeling the clocks, and we direct the reader to [46,47,48] for more information on the subject. Furthermore, we do not consider applications in which the line of sight (LoS) between platforms is obstructed, but this subject is discussed in greater detail in [16,49,50].

3. Problem Formulation

In this section, we define the problem statement, summarize the relevant terminology, and establish the underlying assumptions used in the remainder of the document.

3.1. Problem Setup

The TWR system presented in [5] enables localization, synchronization, and communications between a network of multi-antenna users. For simplicity, we limit the scope of this discussion to the interactions between just two users labeled A and B, as depicted in Figure 1. In this example, user A is a stationary ground station and node B is a UAV moving with some velocity v and acceleration a .
Each transmission takes some time τ to propagate between these users. We define the relative velocity and acceleration of one node with respect to the other in the direction of the line joining them as radial velocity ( τ ˙ ) and acceleration ( τ ¨ ). Users A and B are driven by independent clocks, which at any given time read t A and t B . The relative time offset T is the time difference between the two clocks, T = t A t B . By convention, a positive T denotes that clock B displays an earlier time than clock A. The relative frequency offset and drift between these two clocks are denoted T ˙ and T ¨ , respectively. These interactions and variable definitions are summarized in Figure 2.

3.2. Timing Exchange Protocol

Users A and B engage in a cooperative timing exchange in which they alternate transmitting and receiving a waveform that contains both communications feedback and navigation reference sequences. These exchanges are depicted in Figure 2. The transmit timestamps are deterministic and shared, while the receive timestamps are estimated using a variety of ToA estimation techniques. These timestamps are denoted by t ( . ) , ( . ) ( . ) ; the first subscript indicates the node at which the event occurs, the second subscript indicates whether it was a transmit or receive event, and the superscript indicates in which frame the event occurred. Two successive frames comprise a cycle that is L seconds long and is denoted { ( n 1 ) , ( n ) } where the successive frames ( n 1 ) and ( n ) are l seconds apart. Despite rigorous scheduling, the nodes disagree on time, so these cycle and frame lengths are variable and different at each node. These definitions are summarized in the Table 1.
In each cycle, user A transmits a message to B in the frame ( n 1 ) , B waits for an agreed frame separation l, and transmits a response to A during the next frame ( n ) . Each message contains the transmit and receive timestamps t ( . ) , T x , t ( . ) , R x from the previous frame. This exchange enables estimation of the clock offset T and relative ToF τ between the two users for each cycle.
For a transmission from A to B during frame ( n 1 ) , node B will receive the signal at time t B , R x ( n 1 ) , which is a function of the transmit time t A , T x ( n 1 ) , the time delay τ ( n 1 ) , and the time offset T ( n 1 ) , given by:
t B , R x ( n 1 ) = t A , T x ( n 1 ) + τ ( n 1 ) T ( n 1 ) .
Similarly, for a transmission from B to A during frame ( n ) , node A will receive the signal at time t A , R x ( n ) , given by:
t A , R x ( n ) = t B , T x ( n ) + τ ( n ) + T ( n ) .
For convenience, we define the transmit timestamp t B , T x ( n ) as perceived by node A as t ˜ A , T x ( n ) , which includes the clock offset depicted in Figure 2, as:
t ˜ A , T x ( n ) = t B , T x ( n ) + T ( n ) .
At user A, the frame length l is perceived as l A , and the cycle separation L is perceived as L A , which for the current cycle becomes:
l A ( n 1 ) = t ˜ A , T x ( n ) t A , T x ( n 1 ) ,
L A ( n 1 ) = t A , T x ( n 1 ) t A , T x ( n 3 ) .
For conciseness, we derive the estimators and tracking algorithms for T and τ at user A. The results are symmetrical and can be applied at user B with trivial substitutions but are excluded for brevity.

4. Optimal One-Shot Estimators

In this section, we extend the estimators presented in [19,35] by removing the simplifying assumptions and generalizing the results to higher-order models [9].

4.1. First-Order Models

Here, we assume that the relative propagation delay τ and time offset T between two nodes A and B follow a first-order Markov model. This implies that node B moves with constant radial velocity τ ˙ and the clocks on the two radios have constant relative frequency offset T ˙ . During the two-way cycle { ( n 1 ) , ( n ) } , the delay and offset in frame ( n ) are expressed as functions of the previous frame ( n 1 ) as
τ ( n ) = τ ( n 1 ) + τ ˙ ( n 1 ) l A ( n 1 ) ,
T ( n ) = T ( n 1 ) + T ˙ ( n 1 ) l A ( n 1 ) ,
where l A ( n 1 ) is the separation between the two frames. By substituting Equation (1) and some simplifying algebra, we estimate the first-order frame length l A 1 , ( n 1 ) as
l A 1 , ( n 1 ) = t B , T x ( n ) t A , T x ( n 1 ) + T ( n 1 ) 1 T ˙ ( n 1 ) .
Using Equations (6)–(8), we can simplify Equations (1) and (2) into a system of linear equations in τ ( n 1 ) and T ( n 1 ) .
τ ( n 1 ) T ( n 1 ) = δ ( n 1 ) ,
ε ( n 1 ) τ ( n 1 ) + ζ ( n 1 ) T ( n 1 ) = η ( n 1 ) ,
where δ , ϵ , ζ , and η are scalar coefficients defined as:
δ ( n 1 ) = t B , R x ( n 1 ) t A , T x ( n 1 ) ,
ε ( n 1 ) = L A ( n 1 ) + t B , T x ( n ) t A , T x ( n 1 ) + T ( n 3 ) ,
ζ ( n 1 ) = L A ( n 1 ) + t A , R x ( n ) t A , T x ( n 1 ) τ ( n 3 ) ,
η ( n 1 ) = t A , R x ( n ) ( T ( n 3 ) + L A ( n 1 ) ) + t B , T x ( n ) ( τ ( n 3 ) L A ( n 1 ) ) t A , T x ( n 1 ) ( τ ( n 3 ) + T ( n 3 ) ) .
Using Equations (9) and (), we can construct closed-form estimates of T and τ at node A for frame ( n 1 ) :
τ ( n 1 ) = η ( n 1 ) + ζ ( n 1 ) δ ( n 1 ) ε ( n 1 ) + ζ ( n 1 ) ,
T ( n 1 ) = η ( n 1 ) ε ( n 1 ) δ ( n 1 ) ε ( n 1 ) + ζ ( n 1 ) .
Since our two-way ranging configuration is symmetrical, the delay and offset estimators defined in Equations (15) and (16) can be implemented at both nodes A and B with a trivial change of variables or by following the same process outlined in this section with the appropriate substitutions. We therefore exclude the equivalent derivations for node B for brevity for the remainder of the manuscript.

4.2. Second-Order Models

The first-order solutions presented in Equations (15) and (16) are closed-form, computationally tractable, and relatively simple to implement. If the underlying behavior of the platform motion or the oscillator imperfections is not well modeled by a first-order Markov model, then this solution may underperform. In some cases, a second-order Markov model may better describe the behavior of these variables; so, we extend the derivation in the previous section to include an additional derivative on T and τ .
Here, we assume that the relative propagation delay τ and time offset T between two nodes A and B follow a second-order Markov model. During the two-way cycle { ( n 1 ) , ( n ) } , the delay and offset in frame ( n ) are expressed as functions of the previous frame ( n 1 ) as
τ ( n ) = τ ( n 1 ) + τ ˙ ( n 1 ) l A ( n 1 ) + 1 2 τ ¨ ( n 1 ) l A ( n 1 ) 2 ,
T ( n ) = T ( n 1 ) + T ˙ ( n 1 ) l A ( n 1 ) + 1 2 T ¨ ( n 1 ) l A ( n 1 ) 2 .
Using Equations (17) and (18), we can substitute and simplify Equations (1) and (2) into a system of equations given by:
1 2 α ( n 1 ) l A ( n 1 ) 2 β ( n 1 ) l A ( n 1 ) + γ ( n 1 ) = 0 ,
where the scalar coefficients α , β , and γ are given by:
1 1 α ( n 1 ) = T ¨ ( n 1 ) ,
β ( n 1 ) = 1 T ˙ ( n 1 ) ,
n 1 = T ( n 1 ) + t B , T x ( n ) t A , T x ( n 1 ) .
To solve this system, let
ν ( n 1 ) = β ^ ( n 1 ) 2 2 α ( n 1 ) γ ( n 1 ) .
A tractable solution to this quadratic equation, for a reasonable set of parameters, is
l A 2 , ( n 1 ) = β ( n 1 ) ν ( n 1 ) α ( n 1 ) ,
which exists if and only if ν ( n 1 ) 0 . By substituting Equation (24), Equations (17) and (18) reduce to a system of linear equations in τ ( n 1 ) and T ( n 1 ) . It is not convenient to explicitly write out these estimates, so we rely on real-time solvers to perform the estimation.
These first- and second-order estimators are referred to as “one-shot” estimators because they use a single set of measurements to generate a single set of estimates. These one-shot formulations are advantageous because they are relatively straightforward to understand and implement on an actual system, where accessing memory to perform iterative or memory-driven techniques is notoriously expensive. Unfortunately, these formulations are also sensitive to any channel impairments that may corrupt the measurements. This sensitivity can be mitigated by using any number of tracking filter approaches, which typically incur the aforementioned memory access penalties. In the following section, we derive an efficient extended Kalman filter that mitigates this penalty and still retains the near-optimal performance of the original one-shot estimators.

5. Extended Kalman Filter Tracking

In this section, we formulate an extended Kalman filter to track the propagation delay τ and clock offset T between two users A and B. The terminology introduced in this section is summarized in Table 2.

5.1. Tracking Preliminaries

In the previous section, we demonstrated that the timestamps from two successive frames { ( n 1 ) , ( n ) } contain enough information to estimate T and τ (see Equations (1) and (2)). Therefore, we employ tracking once every cycle, producing state space estimates at a resolution of 1 cycle = two frames. Kalman filters inherently assume Gaussian processes, which nonlinear problems fails to preserve. Extended Kalman Filters (EKF) [38] counter this issue by linearizing the nonlinear transformations, thus preserving Gaussianity. The algorithm can be visualized in two steps: (a) prediction and (b) correction.
The state space parameters x ^ at time instance ( n 1 ) are predicted using estimates from the previous cycle x ^ ( n 3 ) using the standard transition model given by
x ^ ( ) ( n 1 ) = F L A ( n 1 ) x ^ ( n 3 ) + w ( n 1 ) ,
where F ( · ) is the state transition, L A is the cycle length defined in Equation (5), and w ( · ) N ( 0 , Q ( . ) ) is the process noise, assumed to be drawn from a zero-mean multivariate normal distribution. The standard error covariance matrix P ^ in predicting state space parameters can be written as a function of the state transition F ( · ) and the process noise covariance matrix Q , given by:
P ^ ( ) ( n 1 ) = F ( n 1 ) P ^ ( n 3 ) F ( n 1 ) T + Q ( n 1 ) ,
where P ^ ( n 3 ) is the covariance matrix from the previous cycle. The standard measurement prediction z ^ is estimated using
z ^ ( ) ( n 1 ) = u ( n 1 ) + h x ^ ( ) ( n 1 ) + v ( n 1 ) ,
where u ^ represents the control parameters, h ( · ) is the nonlinear measurement transition function, and v ( · ) N ( 0 , R ( · ) ) is the measurement noise assumed to be drawn from a zero-mean multivariate normal distribution. We use the measurement Jacobian H , defined as
H ( n 1 ) = h ( x ^ ) x ^ x ^ ( ) ( n 1 ) ,
as the acting measurement transition function. The error covariance matrix in predicting measurements S is also typically written as a function of the acting measurement transition function H and the measurement noise covariance matrix R , which takes the same form as Equation (26):
S ( n 1 ) = H ( n 1 ) P ^ ( ) ( n 1 ) H ( n 1 ) T + R ( n 1 ) .
Using these predictions, we correct for the state space variables by evaluating a weighted sum of the state predictions and deviation in measurement predictions from observations. To do so, we first define a weighting factor K by normalizing the acting transition function H by the error covariance matrices in the predicting parameters and measurements, P ^ and S , respectively. This takes the form:
K ( n 1 ) = P ^ ( ) ( n 1 ) H ( n 1 ) T S ( n 1 ) 1 ,
with which we make the corrections
x ^ ( n 1 ) = x ^ ( ) ( n 1 ) + K ( n 1 ) ( z ( n 1 ) z ^ ( ) ( n 1 ) ) ,
P ^ ( n 1 ) = ( I | x | K ( n 1 ) H ( n 1 ) ) P ^ ( ) ( n 1 ) ,
where | x | is the number of state space parameters. To maintain the resolution of one estimate per frame, we extrapolate between cycles; the estimates x ^ ( n ) are extrapolated using
x ^ ( n ) = F l A ( n 1 ) x ^ ( n 1 ) .
Using this construction, we formulate the first- and second-order joint delay–offset tracking algorithm described in Algorithm 1.
Algorithm 1: Extended Kalman Filter Tracking Algorithm
1:
while flight duration do
2:
     x ^ ( ) ( n 1 ) = F ( L A ( n 1 ) ) x ^ ( n 3 )
3:
     P ^ ( ) ( n 1 ) = F ( n 1 ) P ^ ( n 3 ) F ( n 1 ) T + Q ( n 1 )
4:
     z ^ ( ) ( n 1 ) = u ( n 1 ) + h ( x ^ ( ) ( n 1 ) )
5:
     H ( n 1 ) = h ( x ^ ) x ^ x ^ ( ) ( n 1 )
6:
     S ( n 1 ) = H ( n 1 ) P ^ ( ) ( n 1 ) H ( n 1 ) T + R ( n 1 )
7:
     z ^ ( ) ( n 1 ) = u ( n 1 ) + h ( x ^ ( ) ( n 1 ) )
8:
     K ( n 1 ) = P ^ ( ) ( n 1 ) H ( n 1 ) T S ( n 1 ) 1
9:
     x ^ ( n 1 ) = x ^ ( ) ( n 1 ) + K ( n 1 ) ( z ( n 1 ) z ^ ( ) ( n 1 ) )
10:
     P ^ ( n 1 ) = ( I | x | K ( n 1 ) H ( n 1 ) ) P ^ ( ) ( n 1 )
end

5.2. First-Order Extended Kalman Filter

In this section, we assume that the relative propagation delay τ and time offset T between two nodes A and B follow a first-order Markov model. We model the progression of the state space parameters x using Equation (25), where the state space variables x and transition matrix F at any arbitrary frame are
x = τ τ ˙ T T ˙ , F l = 1 l 0 0 0 1 0 0 0 0 1 l 0 0 0 1 ,
and the process noise covariance matrix is simply modeled as
Q l = diag Q τ ( l ) , Q T ( l ) ,
where Q τ and Q T are the covariances the in delay and offset transition processes, repespectively (implying that the two are uncorrelated). The delay covariance Q τ = cov [ τ , τ ˙ ] under a Gaussian velocity perturbation model is derived in [44] and is given by:
Q τ l = 1 2 σ τ ˙ 2 l 2 σ τ ˙ 2 l σ τ ˙ 2 l σ τ ˙ 2 ,
where σ τ ˙ 2 is the variance of the relative velocity and relies on the perturbations induced by environmental attributes such as wind currents. The offset covariance Q T = cov [ T , T ˙ ] is the result of a two-state clock error model, which is derived in [45] and is given by:
Q T L A = σ T 2 l + 1 3 σ T ˙ 2 l 3 1 2 σ T ˙ 2 l 2 1 2 σ T ˙ 2 l 2 σ T ˙ 2 l ,
where σ T 2 and σ T ˙ 2 are the variances in estimating the relative clock and frequency offsets between the two radios, respectively. We previously defined the measurement transition in Equation (27) and now define and substitute the measurements z and controls u as:
z ( n 1 ) = t B , R x ( n 1 ) t A , R x ( n ) , u ( n 1 ) = t A , T x ( n 1 ) t B , T x ( n ) .
The transition function is formulated by simply substituting the transition model defined in Equations (6) and (7) (the first-order Markov model), which is written as:
h x ^ ( ) ( n 1 ) = τ ( n 1 ) T ( n 1 ) τ ( n 1 ) + τ ˙ ( n 1 ) l A ( n 1 ) + T ( n 1 ) + T ˙ ( n 1 ) l A ( n 1 ) ,
where the first-order frame length estimate is given by Equation (8), and the measurement covariance matrix can be computed using
R = diag σ t A , R x 2 , σ t B , R x 2 ,
where the variances σ t A , R x 2 and σ t B , R x 2 are the variances of the ToA estimators that generate the measurements [51]. Following the definition presented in Equation (28), the measurement Jacobian may be written explicitly as:
H ( n 1 ) = h ( x ^ ) τ h ( x ^ ) τ ˙ h ( x ^ ) T h ( x ^ ) T ˙ x ^ ( ) ( n 1 ) ,
where each partial derivative is included below for convenience:
h ( x ^ ) τ x ^ ( ) ( n 1 ) = 1 1 ,
h ( x ^ ) τ ˙ x ^ ( ) ( n 1 ) = 0 t B , T x ( n ) t A , T x ( n 1 ) + T ( n 1 ) 1 T ˙ ( n 1 ) ,
h ( x ^ ) T x ^ ( ) ( n 1 ) = 1 1 + τ ˙ ( n 1 ) 1 T ˙ ( n 1 ) ,
h ( x ^ ) T ˙ x ^ ( ) ( n 1 ) = 0 ( 1 + τ ˙ ( n 1 ) ) ( t B , T x ( n ) t A , T x ( n 1 ) + T ( n 1 ) ) ( 1 T ˙ ( n 1 ) ) 2 .

5.3. Second Order Models

In this section, we assume that the relative propagation delay τ and time offset T instead follow the second-order Markov model defined in Equations (17) and (18). The second-order transition matrix has the structure F ( L A ) = diag f ( L A ) , f ( L A ) , where the state space variables x and transition matrix f ( L A ) are defined as:
x = τ , τ ˙ , τ ¨ , T , T ˙ , T ¨ , f ( l ) = 1 l 1 2 l 2 0 1 l 0 0 1 .
The delay covariance Q τ = cov [ τ , τ ˙ , τ ¨ ] under a Gaussian acceleration perturbation model is derived in [44] and is given by:
Q τ l = 1 4 σ τ ¨ 2 l 4 1 2 σ τ ¨ 2 l 3 1 2 σ τ ¨ 2 l 2 1 2 σ τ ¨ 2 l 3 1 2 σ τ ¨ 2 l 2 σ τ ¨ 2 l 1 2 σ τ ¨ 2 l 2 σ τ ¨ 2 l σ τ ¨ 2 ,
where σ τ ¨ 2 is the variance in relative acceleration. The offset covariance Q T = cov [ T , T ˙ , T ¨ ] is a reflection of the clock errors expressed as functions of their Allan variance characteristics. This offset covariance is derived in [41] and is given by:
Q T l = σ T 2 q 1 l + σ T ˙ 2 q 2 l + σ T ¨ 2 q 3 l ,
q 1 l = l 0 0 0 0 0 0 0 0 , q 2 l = 1 3 l 3 1 2 l 2 0 1 2 l 2 l 0 0 0 0 , q 3 l = 1 20 l 5 1 8 l 4 1 6 l 3 1 8 l 4 1 3 l 3 1 2 l 2 1 6 l 3 1 2 l 2 l .
The model parameters σ T 2 , σ T ˙ 2 , and σ T ¨ 2 may be rigorously defined as functions of the Allan variances of the oscillators that drive these radios, but this examination is beyond the scope of this report. We refer the reader to [41] for more detailed information.
The measurement prediction model is given in Equation (27), where the measurements z and controls u remain unchanged from Equation (38), but the transition function is h ( · ) now written as
h x ^ ( ) ( n 1 ) = τ ( n 1 ) T ( n 1 ) τ ( n 1 ) + τ ˙ ( n 1 ) l A ( n 1 ) + τ ¨ ( n 1 ) l A ( n 1 ) 2 + T ( n 1 ) + T ˙ ( n 1 ) l A ( n 1 ) + T ¨ ( n 1 ) l A ( n 1 ) 2 ,
where the second-order frame length estimate [9] is given by Equation (24). Following the definition presented in Equation (41), the measurement Jacobian is now given by Equation (51). We compute and simplify these six partial derivatives using Equation (50) and summarize the results in Equations (52)–(57). These derivatives themselves include partial derivatives of l A with respect to T, T ˙ , and T ¨ , which we can compute using the definition provided in Equation (23). We include the simplified results in Equations (58)–(60) for convenience.
h ( x ^ ) x ^ x ^ ( ) ( n 1 ) = h ( x ^ ) τ h ( x ^ ) τ ˙ h ( x ^ ) τ ¨ h ( x ^ ) T h ( x ^ ) T ˙ h ( x ^ ) T ¨ x ^ ( ) ( n 1 ) ,
h ( x ^ ) τ x ^ ( ) ( n 1 ) = 1 1 ,
h ( x ^ ) τ ˙ x ^ ( ) ( n 1 ) = 0 l A ( n 1 ) ,
h ( x ^ ) τ ¨ x ^ ( ) ( n 1 ) = 0 1 2 l A ( n 1 ) 2 ,
h ( x ^ ) T x ^ ( ) ( n 1 ) = 1 1 + τ ˙ ( n 1 ) + T ˙ ( n 1 ) l A ( x ^ ) T x ^ ( ) ( n 1 ) + τ ¨ ( n 1 ) + T ¨ ( n 1 ) l A ( n 1 ) l A ( x ^ ) T x ^ ( ) ( n 1 ) ,
h ( x ^ ) T ˙ x ^ ( ) ( n 1 ) = 0 l A ( n 1 ) + τ ˙ ( n 1 ) + T ˙ ( n 1 ) l A ( x ^ ) T ˙ x ^ ( ) ( n 1 ) + τ ¨ ( n 1 ) + T ¨ ( n 1 ) l A ( n 1 ) l A ( x ^ ) T ˙ x ^ ( ) ( n 1 ) ,
h ( x ^ ) T ¨ x ^ ( ) ( n 1 ) = 0 1 2 l A ( n 1 ) 2 + τ ˙ ( n 1 ) + T ˙ ( n 1 ) l A ( x ^ ) T ¨ x ^ ( ) ( n 1 ) + τ ¨ ( n 1 ) + T ¨ ( n 1 ) l A ( n 1 ) l A ( x ^ ) T ¨ x ^ ( ) ( n 1 ) ,
l A ( x ^ ) T x ^ ( ) ( n 1 ) = 1 ν ( n 1 ) ( see Equation ( 23 ) ) ) ,
l A ( x ^ ) T ˙ x ^ ( ) ( n 1 ) = 1 T ¨ ( n 1 ) T ˙ ( n 1 ) 1 ν ( n 1 ) + 1 ,
l A ( x ^ ) T ¨ x ^ ( ) ( n 1 ) = T ( n 1 ) t A , T x ( n 1 ) + t B , T x ( n 1 ) T ¨ ( n 1 ) ν ( n 1 ) + ν ( n 1 ) + T ˙ ( n 1 ) 1 T ¨ ( n 1 ) 2 .

6. Simulation Results

In this section, we describe the creation of a simple MATLAB simulation platform and compare the simulated performance of the optimal one-shot estimators and the proposed extended Kalman filter solutions. We implemented the configuration depicted in Figure 1 and allowed user B to follow the 60 s flight trajectory depicted in Figure 3. The users interacted every 50 ms. The transmit timestamps were known and shared, while the receive timestamps were estimated to precision σ = 0.1 ns, which we have previously demonstrated experimentally [4]. We assumed that the noise covariance matrices Q and R were known. We considered the real-time adaptive estimation of these measures in an adjacent publication. We implemented the optimal one-shot estimators from [9] to benchmark the proposed extended Kalman filters. We measured the performance using the root-mean-square error (RMSE) on the propagation delay and clock offset. We compare the delay and offset estimation performance in Figure 4 and Figure 5, respectively.

6.1. Simulated Computation Time

We implemented the optimal first- and second-order one-shot estimators defined in [9] and the proposed first- and second-order EKF estimators in MATLAB. We coarsely profiled the average computation time of each of these four solutions in 1000 Monte Carlo simulations over the flight trajectory depicted in Figure 3. MATLAB is a notoriously time-inefficient processing architecture; so, we only used this information to make claims about the relative complexity of each solution, not the actual complexity in an integrated system. These results are summarized in Table 3, which demonstrate that the proposed first-order solution ran slightly slower than the one-shot equivalent, but the proposed second-order solution reduced the run time by over an order of magnitude compared to the second-order one-shot solution.

6.2. Simulated Time Delay Estimation Performance

In Figure 4, we plot the time delay estimates τ ^ for the duration of the simulated flight path. In Figure 4a, we plot the true time delay τ in red and the radial velocity τ ˙ in blue. In Figure 4b,c, we plot the estimates generated by the optimal first- and second-order one-shot estimators described in [9], labeled STP1 and STP2, respectively. Both of these estimators successfully estimated the time delay to a precision of about σ = 1 cm and the radial velocity to a precision of about σ = 15 cm/s. In Figure 4d,e, we plot the estimates produced by the proposed first- and second-order EKF tracking solutions, labeled EKF1 and EKF2, respectively. The first-order method, depicted in Figure 4d, achieved slightly better performance and converged to a solution after only a single interaction cycle. The second-order method depicted in Figure 4e, however, took several cycles to converge, creating a direct tradeoff between the computational complexity and convergence speed.
Figure 4. (a) True time of flight and radial velocity, (b,c) existing first- and second-order one-shot estimators [9], labeled STP1 and STP2, respectively, and (d,e) our proposed first- and second-order extended Kalman filter estimators labeled EKF1 and EKF2, respectively.
Figure 4. (a) True time of flight and radial velocity, (b,c) existing first- and second-order one-shot estimators [9], labeled STP1 and STP2, respectively, and (d,e) our proposed first- and second-order extended Kalman filter estimators labeled EKF1 and EKF2, respectively.
Signals 04 00023 g004

6.3. Simulated Time Offset Estimation Performance

In Figure 5, we plotted the time offset estimates T ^ for the duration of the simulated flight path. In Figure 5a, we plotted the true time offset T in red and the clock drift T ˙ in blue. In Figure 5b,c, we plotted the estimates generated by the optimal first- and second-order one-shot estimators described in [9], labeled STP1 and STP2, respectively. Both of these estimators successfully estimated the time offset to a precision of about σ = 0.05 ns and the clock drift to a precision of about σ = 1 ns/s. In Figure 5d,e, we plotted the estimates produced by the proposed first- and second-order EKF tracking solutions, labeled EKF1 and EKF2, respectively. The first-order method depicted in Figure 4d achieved comparable performance and converged to a solution after only a single interaction cycle, while the second-order method depicted in Figure 5e again took several cycles to converge, displaying the same tradeoff between the computational complexity and convergence speed as the previous results.
Figure 5. (a) True time offset and drift, (b,c) existing first- and second-order one-shot estimators [9] labeled STP1 and STP2, respectively, and (d,e) our proposed first- and second-order extended Kalman filter estimators labeled EKF1 and EKF2, respectively.
Figure 5. (a) True time offset and drift, (b,c) existing first- and second-order one-shot estimators [9] labeled STP1 and STP2, respectively, and (d,e) our proposed first- and second-order extended Kalman filter estimators labeled EKF1 and EKF2, respectively.
Signals 04 00023 g005

7. Conclusions

In this manuscript, we proposed an extended Kalman filter design for tracking time-of-flight (ToF) and clock offsets in a distributed two-way ranging (TWR) network. We proposed this solution to specifically address the sensitivity of the existing one-shot approaches to intermittent channel impairments. Tracking solutions often increase the computational complexity of the estimation processes, which we explicitly mitigate by integrating the estimation process into the EKF solution rather than filtering the outputs of another estimator. We implemented the proposed solutions in a simple MATLAB simulation environment and benchmarked their performance against the corresponding optimal one-shot estimators in the Monte Carlo simulations. We demonstrated that this solution achieves comparable performance to existing, optimal estimators and, in the case of the second-order formulation, actually reduces the computation run time by an order of magnitude.
In the future, we will extend the initialization of the noise covariance matrices to include a real-time adaptive estimation process instead of initializing them using fixed values based on either knowledge of the estimation problem or initial calibration. This is commonly referred to as adaptive extended Kalman filtering (AEKF) and makes the system more robust to wide changes in received signal strength or other channel impairments. In parallel, we are also currently investigating computationally efficient methods for extending these results to higher-order Markov models on the underlying processes.

Author Contributions

Conceptualization, S.S. and A.H.; Data curation, S.S.; Formal analysis, S.S. and A.H.; Funding acquisition, D.W.B.; Investigation, S.S. and A.H.; Methodology, S.S. and A.H.; Project administration, D.W.B.; Resources, D.W.B.; Software, S.S.; Supervision, D.W.B.; Validation, S.S. and A.H.; Visualization, S.S., A.H. and D.W.B.; Writing—original draft, S.S. and A.H.; Writing—review and editing, S.S., A.H. and D.W.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the policies of Arizona State University regarding pending patents and publications.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AoAangle of arrival
EKFextended Kalman filter
FANETflying ad hoc network
GPSGlobal Positioning System
IoTInternet of Things
ITSintelligent transport system
LoSline of sight
NTPNetwork Timing Protocol
RMSEroot mean square error
RSSreceived signal strength
TDoAtime difference of arrival
ToAtime of arrival
ToFtime of flight
TWRtwo-way ranging
UAVunmanned aerial vehicle
WSNwireless sensor network

References

  1. Paul, B.; Chiriyath, A.R.; Bliss, D.W. Survey of RF communications and sensing convergence research. IEEE Access 2017, 5, 252–270. [Google Scholar] [CrossRef]
  2. Chiriyath, A.R.; Paul, B.; Jacyna, G.M.; Bliss, D.W. Inner bounds on performance of radar and communications co-existence. IEEE Trans. Signal Process. 2015, 64, 464–474. [Google Scholar] [CrossRef]
  3. Chiriyath, A.R.; Paul, B.; Bliss, D.W. Radar-Communications Convergence: Coexistence, Cooperation, and Co-Design. IEEE Trans. Cogn. Commun. Netw. 2017, 3, 1–12. [Google Scholar] [CrossRef]
  4. Herschfelt, A.; Yu, H.; Wu, S.; Lee, H.; Bliss, D.W. Joint Positioning-Communications System Design: Leveraging Phase-Accurate Time-of-Flight Estimation and Distributed Coherence. In Proceedings of the 2018 52nd Asilomar Conference on Signals, Systems, and Computers, Pacific Grove, CA, USA, 28–31 October 2018; pp. 433–437. [Google Scholar]
  5. Herschfelt, A.; Yu, H.; Wu, S.; Srinivas, S.; Li, Y.; Sciammetta, N.; Smith, L.; Rueger, K.; Lee, H.; Chakrabarti, C.; et al. Joint positioning-communications system design and experimental demonstration. In Proceedings of the 2019 IEEE/AIAA 38th Digital Avionics Systems Conference (DASC), San Diego, CA, USA, 8–12 September 2019; pp. 1–6. [Google Scholar]
  6. Shepard, D.P.; Bhatti, J.A.; Humphreys, T.E.; Fansler, A.A. Evaluation of smart grid and civilian UAV vulnerability to GPS spoofing attacks. In Proceedings of the 25th International Technical Meeting of the Satellite Division of The Institute of Navigation, Nashville, TN, USA, 17–21 September 2012. [Google Scholar]
  7. Humphreys, T. Statement on the Vulnerability of Civil Unmanned Aerial Vehicles and Other Systems to Civil GPS Spoofing; University of Texas at Austin: Austin, TX, USA, 2012; pp. 1–16. [Google Scholar]
  8. Kerns, A.J.; Shepard, D.P.; Bhatti, J.A.; Humphreys, T.E. Unmanned aircraft capture and control via GPS spoofing. J. Field Robot. 2014, 31, 617–636. [Google Scholar] [CrossRef]
  9. Srinivas, S.; Herschfelt, A.; Bliss, D.W. Joint positioning-communications system: Optimal distributed coherence and positioning estimators. In Proceedings of the 2019 53rd Asilomar Conference on Signals, Systems, and Computers, Pacific Grove, CA, USA, 3–6 November 2019; pp. 317–321. [Google Scholar]
  10. Sayed, A.H.; Tarighat, A.; Khajehnouri, N. Network-based wireless location: Challenges faced in developing techniques for accurate wireless location information. IEEE Signal Process. Mag. 2005, 22, 24–40. [Google Scholar] [CrossRef]
  11. Gustafsson, F.; Gunnarsson, F. Mobile positioning using wireless networks: Possibilities and fundamental limitations based on available wireless network measurements. IEEE Signal Process. Mag. 2005, 22, 41–53. [Google Scholar] [CrossRef]
  12. Gezici, S. A survey on wireless position estimation. Wirel. Pers. Commun. 2008, 44, 263–282. [Google Scholar] [CrossRef]
  13. Liu, H.; Darabi, H.; Banerjee, P.; Liu, J. Survey of wireless indoor positioning techniques and systems. IEEE Trans. Syst. Man, Cybern. Part C (Appl. Rev.) 2007, 37, 1067–1080. [Google Scholar] [CrossRef]
  14. Patwari, N.; Ash, J.N.; Kyperountas, S.; Hero, A.O.; Moses, R.L.; Correal, N.S. Locating the nodes: Cooperative localization in wireless sensor networks. IEEE Signal Process. Mag. 2005, 22, 54–69. [Google Scholar] [CrossRef]
  15. Mesmoudi, A.; Feham, M.; Labraoui, N. Wireless sensor networks localization algorithms: A comprehensive survey. arXiv 2013, arXiv:1312.4082. [Google Scholar] [CrossRef]
  16. Guvenc, I.; Chong, C.C. A survey on TOA based wireless localization and NLOS mitigation techniques. IEEE Commun. Surv. Tutorials 2009, 11, 107–124. [Google Scholar] [CrossRef]
  17. Lanzisera, S.; Zats, D.; Pister, K.S. Radio frequency time-of-flight distance measurement for low-cost wireless sensor localization. IEEE Sensors J. 2011, 11, 837–845. [Google Scholar] [CrossRef]
  18. Weiss, A.; Weinstein, E. Fundamental limitations in passive time delay estimation–Part I: Narrow-band systems. IEEE Trans. Acoust. Speech Signal Process. 1983, 31, 472–486. [Google Scholar] [CrossRef]
  19. Bidigare, P.; Madhow, U.; Mudumbai, R.; Scherber, D. Attaining fundamental bounds on timing synchronization. In Proceedings of the 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Kyoto, Japan, 25–30 March 2012; pp. 5229–5232. [Google Scholar]
  20. Lewandowski, W.; Thomas, C. GPS time transfer. Proc. IEEE 1991, 79, 991–1000. [Google Scholar] [CrossRef]
  21. Lewandowski, W.; Azoubib, J.; Klepczynski, W.J. GPS: Primary tool for time transfer. Proc. IEEE 1999, 87, 163–172. [Google Scholar] [CrossRef]
  22. Sivrikaya, F.; Yener, B. Time synchronization in sensor networks: A survey. IEEE Netw. 2004, 18, 45–50. [Google Scholar] [CrossRef]
  23. Ganeriwal, S.; Kumar, R.; Srivastava, M.B. Timing-sync protocol for sensor networks. In Proceedings of the 1st International Conference on Embedded Networked Sensor Systems, Los Angeles, CA, USA, 5–7 November 2003; pp. 138–149. [Google Scholar]
  24. Li, Q.; Rus, D. Global clock synchronization in sensor networks. IEEE Trans. Comput. 2006, 55, 214–226. [Google Scholar]
  25. Schenato, L.; Gamba, G. A distributed consensus protocol for clock synchronization in wireless sensor network. In Proceedings of the 2007 46th IEEE Conference on Decision and Control, New Orleans, LA, USA, 12–14 December 2007; pp. 2289–2294. [Google Scholar]
  26. Maggs, M.K.; O’keefe, S.G.; Thiel, D.V. Consensus clock synchronization for wireless sensor networks. IEEE Sensors J. 2012, 12, 2269–2277. [Google Scholar] [CrossRef]
  27. Mills, D.L. Internet time synchronization: The network time protocol. IEEE Trans. Commun. 1991, 39, 1482–1493. [Google Scholar] [CrossRef] [Green Version]
  28. Mills, D. Network Time Protocol (version 3) Specification, Implementation and Analysis. Technical Report, Internet Society. 1992. Available online: https://www.rfc-editor.org/rfc/rfc1305 (accessed on 15 May 2023).
  29. Eidson, J.; Lee, K. IEEE 1588 standard for a precision clock synchronization protocol for networked measurement and control systems. In Proceedings of the Sensors for Industry Conference, 2002. 2nd ISA/IEEE, Reston, VA, USA, 3–5 December 2002; Volume 10. [Google Scholar]
  30. Sichitiu, M.L.; Veerarittiphan, C. Simple, accurate time synchronization for wireless sensor networks. In Proceedings of the 2003 IEEE Wireless Communications and Networking, 2003. WCNC 2003, New Orleans, LA, USA, 16–20 March 2003; Volume 2, pp. 1266–1273. [Google Scholar]
  31. Van Greunen, J.; Rabaey, J. Lightweight time synchronization for sensor networks. In Proceedings of the 2nd ACM international conference on Wireless sensor networks and applications, San Diego, CA, USA, 19 September 2003; pp. 11–19. [Google Scholar]
  32. Brown III, D.R.; Klein, A.G.; Wang, R. Non-Hierarchical Clock Synchronization for Wireless Sensor Networks. arXiv 2012, arXiv:1212.1216. [Google Scholar]
  33. Preuss, R.D.; Brown, D.R. Retrodirective distributed transmit beamforming with two-way source synchronization. In Proceedings of the 2010 44th Annual Conference on Information Sciences and Systems (CISS), Princeton, NJ, USA, 17–19 March 2010; pp. 1–6. [Google Scholar]
  34. Brown, D.R.; Klein, A.G. Precise timestamp-free network synchronization. In Proceedings of the 2013 47th Annual Conference on Information Sciences and Systems (CISS), Baltimore, MD, USA, 20–22 March 2013; pp. 1–6. [Google Scholar]
  35. Bidigare, P.; Pruessing, S.; Raeman, D.; Scherber, D.; Madhow, U.; Mudumbai, R. Initial over-the-air performance assessment of ranging and clock synchronization using radio frequency signal exchange. In Proceedings of the Statistical Signal Processing Workshop (SSP), Ann Arbor, MI, USA, 5–8 August 2012; pp. 273–276. [Google Scholar]
  36. Denis, B.; Pierrot, J.B.; Abou-Rjeily, C. Joint distributed synchronization and positioning in UWB ad hoc networks using TOA. IEEE Trans. Microw. Theory Tech. 2006, 54, 1896–1911. [Google Scholar] [CrossRef]
  37. Gelb, A. Applied Optimal Estimation; MIT Press: Cambridge, MA, USA, 1974. [Google Scholar]
  38. Anderson, B.D.; Moore, J.B. Optimal Filtering; Courier Corporation: Chelmsford, MA, USA, 2012. [Google Scholar]
  39. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Fluids Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  40. Galleani, L.; Tavella, P. The characterization of clock behavior with the dynamic Allan variance. In Proceedings of the IEEE International Frequency Control Symposium and PDA Exhibition Jointly with the 17th European Frequency and Time Forum, Tampa, FL, USA, 4–8 May 2003; pp. 239–244. [Google Scholar]
  41. Zucca, C.; Tavella, P. The clock model and its relationship with the Allan and related variances. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 289–296. [Google Scholar] [CrossRef] [PubMed]
  42. Greenhall, C.A. Forming stable timescales from the Jones–Tryon Kalman filter. Metrologia 2003, 40, S335. [Google Scholar] [CrossRef]
  43. Giorgi, G.; Narduzzi, C. Performance analysis of Kalman-filter-based clock synchronization in IEEE 1588 networks. IEEE Trans. Instrum. Meas. 2011, 60, 2902–2909. [Google Scholar] [CrossRef]
  44. Li, X.R.; Jilkov, V.P. Survey of maneuvering target tracking. Part I. Dynamic models. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1333–1364. [Google Scholar]
  45. Galleani, L. A tutorial on the two-state model of the atomic clock noise. Metrologia 2008, 45, S175. [Google Scholar] [CrossRef]
  46. Davis, J.A.; Greenhall, C.; Stacey, P. A Kalman filter clock algorithm for use in the presence of flicker frequency modulation noise. Metrologia 2005, 42, 1. [Google Scholar] [CrossRef] [Green Version]
  47. Rutman, J. Characterization of phase and frequency instabilities in precision frequency sources: Fifteen years of progress. Proc. IEEE 1978, 66, 1048–1075. [Google Scholar] [CrossRef]
  48. Percival, D.B. Stochastic models and statistical analysis for clock noise. Metrologia 2003, 40, S289. [Google Scholar] [CrossRef] [Green Version]
  49. Chan, Y.T.; Tsui, W.Y.; So, H.C.; Ching, P.C. Time-of-arrival based localization under NLOS conditions. IEEE Trans. Veh. Technol. 2006, 55, 17–24. [Google Scholar] [CrossRef]
  50. Seco, F.; Jiménez, A.R.; Prieto, C.; Roa, J.; Koutsou, K. A survey of mathematical methods for indoor localization. In Proceedings of the 2009 IEEE International Symposium on Intelligent Signal Processing, Budapest, Hungary, 26–28 August 2009; pp. 9–14. [Google Scholar]
  51. Herschfelt, A. Simultaneous Positioning and Communications: Hybrid Radio Architecture, Estimation Techniques, and Experimental Validation. Ph.D. Thesis, Arizona State University, Tempe, AZ, USA, 2019. [Google Scholar]
Figure 1. Example system configuration of the previously developed two-way ranging (TWR) system that efficiently enables localization, synchronization, and communications in a distributed radio network. We propose a Kalman filter extension to this system that tracks the time offset T, time delay τ , and the first and second derivatives thereof between the distributed users. These extensions make the system more robust to cycle slips and reduce the computational complexity by an order of magnitude.
Figure 1. Example system configuration of the previously developed two-way ranging (TWR) system that efficiently enables localization, synchronization, and communications in a distributed radio network. We propose a Kalman filter extension to this system that tracks the time offset T, time delay τ , and the first and second derivatives thereof between the distributed users. These extensions make the system more robust to cycle slips and reduce the computational complexity by an order of magnitude.
Signals 04 00023 g001
Figure 2. Interactions between two users A and B. Successive transmissions are separated by l A , and successive cycles are separated by L A . These users exchange timing information each frame to drive the localization algorithm.
Figure 2. Interactions between two users A and B. Successive transmissions are separated by l A , and successive cycles are separated by L A . These users exchange timing information each frame to drive the localization algorithm.
Signals 04 00023 g002
Figure 3. Simulated flight trajectory of user B for 60 s. User A is stationary at the origin.
Figure 3. Simulated flight trajectory of user B for 60 s. User A is stationary at the origin.
Signals 04 00023 g003
Table 1. Timing Model Notation.
Table 1. Timing Model Notation.
t A , T x Transmit event timestamps at node A
t A , R x Receive event timestamps at node A
t B , T x Transmit event timestamps at node B
t B , R x Receive event timestamps at node B
t ˜ A , T x Transmit event time-stamp t B , T x ( . ) at node A
l A Frame length at node A
L A Cycle length at node A
τ Relative ToF
τ ˙ Relative velocity
τ ¨ Relative acceleration
TRelative time offset
T ˙ Relative frequency offset
T ¨ Relative frequency drift
Table 2. Extended Kalman Filter Notation.
Table 2. Extended Kalman Filter Notation.
x ^ ( ) , x ^ Predicted and estimated state parameters
| x | Cardinality of state space
F State transition matrix
z ^ ( ) , z ^ Predicted and observed measurements
u Control parameters
h , H Measurement transition and Jacobian
w , v State and measurement noise
Q , R State and measurement noise covariance matrix
P ^ ( ) , P ^ Predicted and estimated state covariance matrix
S Measurement covariance matrix
K Kalman gain
Table 3. Average Estimator Computation Time Across 1000 Trials.
Table 3. Average Estimator Computation Time Across 1000 Trials.
Optimal One-Shot Estimator—First-Order0.93 ms
Optimal One-Shot Estimator—Second-Order49.0 ms
Extended Kalman Filter Tracking—First-Order1.30 ms
Extended Kalman Filter Tracking—Second-Order2.60 ms
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Srinivas, S.; Herschfelt, A.; Bliss, D.W. Extended Kalman Filter Design for Tracking Time-of-Flight and Clock Offsets in a Two-Way Ranging System. Signals 2023, 4, 439-456. https://doi.org/10.3390/signals4020023

AMA Style

Srinivas S, Herschfelt A, Bliss DW. Extended Kalman Filter Design for Tracking Time-of-Flight and Clock Offsets in a Two-Way Ranging System. Signals. 2023; 4(2):439-456. https://doi.org/10.3390/signals4020023

Chicago/Turabian Style

Srinivas, Sharanya, Andrew Herschfelt, and Daniel W. Bliss. 2023. "Extended Kalman Filter Design for Tracking Time-of-Flight and Clock Offsets in a Two-Way Ranging System" Signals 4, no. 2: 439-456. https://doi.org/10.3390/signals4020023

APA Style

Srinivas, S., Herschfelt, A., & Bliss, D. W. (2023). Extended Kalman Filter Design for Tracking Time-of-Flight and Clock Offsets in a Two-Way Ranging System. Signals, 4(2), 439-456. https://doi.org/10.3390/signals4020023

Article Metrics

Back to TopTop