Next Article in Journal
Effect of Feed Rate on the Force and Energy in the Cutting Process Using Planar Technical Blades
Next Article in Special Issue
Hierarchical Model Updating Method for Vector Electric-Propulsion Satellites
Previous Article in Journal
Diffraction of the Field of a Grounded Cable on an Elongated Dielectric Spheroid in a Conducting Layer
Previous Article in Special Issue
AC Electric Powertrain without Power Electronics for Future Hybrid Electric Aircrafts: Architecture, Design and Stability Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cramér–Rao Lower Bounds on 3D Position and Orientation Estimation in Distributed Ranging Systems

Center for Wireless Information Systems and Computational Architectures (WISCA), Arizona State University, Tempe, AZ 85281, USA
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(3), 2008; https://doi.org/10.3390/app13032008
Submission received: 28 December 2022 / Revised: 26 January 2023 / Accepted: 31 January 2023 / Published: 3 February 2023
(This article belongs to the Special Issue Aerospace Vehicle Design under Uncertainties)

Abstract

:
As radio frequency (RF) hardware continues to improve, many technologies that were traditionally impractical have suddenly become viable alternatives to legacy systems. Two-way ranging (TWR) is often considered a poor positioning solution for airborne and other vehicular navigation systems due to its low precision, poor angular resolution, and precise timing requirements. With the advent of modern RF hardware and advanced processing techniques, however, modern studies have experimentally demonstrated TWR systems with an unprecedented, sub-centimeter ranging precision with low size, weight, power, and cost (SWaP-C) consumer-grade hardware. This technique enables a new class of positioning, navigation, and timing (PNT) capabilities for urban and commercial aircraft but also instigates new system design challenges such as antenna placement, installation of new electronics, and design of supporting infrastructure. To inform these aircraft design decisions, we derive 2D and 3D Cramér–Rao lower bounds (CRLBs) on position and orientation estimation in a multi-antenna TWR system. We specifically formulate these bounds as a function of the number of antennas, platform geometry, and geometric dilution of precision (GDoP) to inform aircraft design decisions under different mission requirements. We simulate the performance of several classic position and orientation estimators in this context to validate these bounds and to graphically depict the expected performance with respect to these design considerations. To improve the accessibility of these highly theoretical results, we also present a simplified discussion of how these bounds may be applied to common airborne applications and suggest best practices for using them to inform aircraft design decisions.

1. Introduction

As we continue to develop high-performance radio technologies, aerospace vehicles have become increasingly capable and versatile; unfortunately, the design of these platforms is a complex, multidisciplinary task rife with numerous design challenges [1]. As a result, many technologies that have demonstrated widespread success in terrestrial applications are not always suitable for airborne applications without careful consideration and redesign [2]. To successfully integrate a novel technology into modern aircraft, we must thoroughly understand both its own behavior and its interactions with the rest of the platform [3].
Positioning, navigation, and timing (PNT) are some of the most critical services for modern aircraft and numerous technologies provide these services through different mechanisms [4]. Two-way ranging (TWR) is traditionally considered a poor positioning solution for aerospace systems, but emerging technologies [5] have demonstrated sub-centimeter ranging performance with minimal size, weight, power, and cost (SWaP-C). These modern solutions not only enable a new class of distributed applications for urban and commercial aircraft but also introduce significant system design challenges such as antenna placement [6], installation of new electronics [7], and the design of supporting infrastructure [8].
To inform these aircraft design decisions, we derive novel 2D and 3D Cramér–Rao lower bounds (CRLBs) on position and orientation estimation in a multi-antenna TWR system. These bounds are both a closed-form and tractable and are specifically formulated as a function of the number of antennas, platform geometry, and geometric dilution of precision (GDoP) to inform design decisions under different mission requirements. These bounds directly relate the position and orientation estimation precision to the geometric distribution of antennas on the platform. This can be used to infer the expected performance of existing installations or to inform the required distribution for new installations to achieve an arbitrary precision target. Since these geometries are often summarized by a single GDoP value, we also trace the performance manifold between ranging precision, positioning precision, and GDoP to graphically summarize the bounds, which can be explored visually without performing any calculations. To improve the accessibility of these highly theoretical results, we also present simplified commentary on how to navigate the trade space between performance and practical design considerations. These results are summarized into simple “best practices” to provide intuitive, geometric interpretations of the bounds. An example multi-antenna TWR configuration to which these bounds apply is depicted in Figure 1.

1.1. Background

The proposed Cramér–Rao lower bounds (CRLBs) on position and orientation estimation apply to estimators that use distance estimates as the primary measurement mechanism. This work was originally motivated by emerging TWR technologies for urban air mobility (UAM) applications, but the results apply to any technique that generates multiple range estimates for a given target. These results apply to several different types of TWR systems and are primarily enabled by previous work in the field of time-of-arrival (ToA) estimation.

1.1.1. Two-Way Ranging

One-way ranging systems such as radars [9] typically measure unresponsive or otherwise-uncooperative targets. Numerous extensions have been developed to refine the performance of this class of systems [10], but the uncooperative nature of the target imposes fundamental limits on the achievable precision [11]. Unlike these systems, twoway ranging (TWR) refers to a class of systems in which both the observer and the target communicate and cooperate to produce ranging estimates [12]. When these technologies were first emerging, radio-frequency (RF) hardware was large, heavy, and expensive, so static, robust, and reliable infrastructure such as a radar installation was a popular solution to the ranging problem [13]. As this hardware becomes more powerful and continues to shrink, however, installing modern TWR electronics on aircraft is significantly more viable [14].
As low-cost RF hardware became more accessible in the early 2000s, double-sided, two-way ranging became a somewhat more popular research topic [15]. Early theoretical work [16] demonstrated several tractable approaches assuming sufficiently behaved oscillators, closely followed by initial experimental demonstrations [17]. Several algorithmic extensions have since been explored [18] as well as more robust experimental demonstrations [19].
A fair amount of this early work was focused on traditional RF frequency bands, but the falling cost of higher-frequency equipment has motivated ultra-wideband (UWB) and even optical extensions. Since the precision of these techniques is proportional to the inverse of the signal bandwidth [20], high-frequency extensions can enable significantly better performance. Some studies have developed direct UWB extensions of the previous results [21], while others have explored the implications of integrating these high-frequency techniques with existing passive infrastructure [22]. While these high-frequency extensions enable intrinsically better performance, they also incur additional calibration challenges, which were characterized in [23].
In parallel, optical two-way ranging techniques have also been explored as high-precision ranging solutions [24]. These high-frequency systems have even better-ranging resolution than RF and UWB systems, but they also tend to have a more limited range and field of view, so they are suitable for different types of applications [25]. Recent studies have explored sophisticated 2D [26] and 3D [27] array processing extensions to further improve this performance. A comprehensive theoretical analysis for this class of systems was presented in [28], and modern studies continue to improve the fidelity of these techniques in real-world demonstrations [29].

1.1.2. Time-of-Arrival Estimation

For most TWR systems, time-of-arrival (ToA) estimation is the primary measurement mechanism that drives the rest of the signal processing chain. Naturally, the precision of these measurements will dictate the overall ranging precision of the system and is the subject of numerous studies in different applications. In indoor environments, researchers have characterized the performance of both wideband [30] and ultra-wideband [31] ToA estimators in dense multipath channels. Similarly for outdoor environments, several ToA estimation techniques were explored in [32,33], whose corresponding lower bounds were investigated in [34]. Other researchers have also investigated the use of existing communications infrastructure to enable these measurements, including LTE [35], 5G [36], and Wi-Fi [37]. Neural networks have also been considered for large networks of ToA-based localization sensors [38].
ToA estimation is well characterized by different types of bounds. Ziv–Zakai and Weiss–Weinstein lower bounds on ToA estimation were derived and discussed in [39] and [40], respectively, while a Cramér–Rao lower bound on hybrid time-of-arrival/received signal strength (ToA/RSS) was derived in [41]. Notably, 2D Cramér–Rao lower bounds on localization were presented in [42], and optimal sensor placement for time-difference-ofarrival (TDoA) systems is discussed in [43]. The CRLBs proposed in this work extend this 2D positioning bound to three dimensions, explicitly includes the number and placement of antennas in the bound, and additionally considers orientation estimation.

1.2. Contributions

In this manuscript, we make the following contributions:
  • Derive a novel, closed-form, tractable CRLB on position estimation in TWR systems;
  • Derive a novel, closed-form, tractable CRLB on orientation estimation in TWR systems;
  • Implement the proposed CRLBs in a simple MATLAB simulation platform;
  • Benchmark several popular estimators against the proposed CRLBs;
  • Discuss how these results can directly inform aircraft design decisions.

2. Two-Way Ranging (TWR) Overview

In this section, we briefly define a timing and propagation model and outline some rudimentary time-of-arrival (ToA) and time-of-flight (ToF) estimation techniques.

2.1. Timing Model

Two-way ranging (as the name implies) involves a cooperative exchange between at least two users. Consider two users labeled A and B separated in space and operating with independent, imperfect clocks. At any given instant n, there is an offset between these clocks labeled T ( n ) and a distance between them d ( n ) . If these users are interacting using over-the-air electromagnetic waveforms, then each transmission takes some time τ ( n ) = d ( n ) / c to propagate between the platforms, where c is the speed of light. These interactions are depicted in Figure 2. The clock offsets T and propagation delays τ are the fundamental quantities of interest and can be estimated by measuring the time-of-arrival (ToA) of each reception and applying an appropriate time-of-flight (ToF) estimation algorithm.
User A transmits a waveform at time instant n, labeled t A , T x ( n ) . This waveform takes a small amount of time to propagate to user B, labeled τ ( n ) . User B receives the waveform at time t B , R x ( n ) . If the clocks were perfectly aligned, the receive timestamp would simply be t A , T x ( n ) + τ ( n ) . Because the clocks are misaligned, user B measures the event earlier with respect to clock B, thus the received timestamp, as perceived by user B, is
t B , R x ( n ) = t A , T x ( n ) + τ ( n ) T ( n ) .
At some later time instant n + 1 , user B transmits a waveform to user A. The received timestamp, as perceived by user A is
t A , R x ( n + 1 ) = t B , T x ( n + 1 ) + τ ( n + 1 ) + T ( n + 1 ) .

2.2. Propagation Model

A signal x is transmitted through a line-of-sight channel, during which it is distorted before the signal z is received. The traditional propagation model [20] for this channel is
z = a x ( t τ ) e j 2 π f c τ ,
where a is the complex channel attenuation, τ is the propagation delay, and f c is the carrier frequency. We assume a basic line-of-sight channel attenuation [44]
a 2 = λ 4 π d 2 G T x G R x ,
where G T x and G R x are the transmitter and receiver antenna gains, λ is the signal wavelength, and d is the distance between the two users.
Advanced ToA estimators require a much higher fidelity propagation model, which was previously explored in significant detail in [45]. The simplified model presented in Equation (3), however, is sufficient to motivate the following discussions and for the remainder of this manuscript.

2.3. Time-of-Arrival Estimation

Equations (1) and (2) form a system of two equations with two unknowns: T and τ . The transmit timestamps t T x are known quantities, and the receive timestamps t R x can be estimated using any number of ToA estimation techniques. ToA estimation is itself a rich field with both legacy and emerging solutions, so we briefly summarize the basic concept here and direct the reader to some relevant publications [40,45,46] for further reading.
Consider the correlation between received signal z ( t ) and known transmitted signal x ( t ) :
g ( τ ) = d t z ( t ) x ( t τ ) 2 ,
where τ is a time delay relative to some time reference (analogous but not equivalent to the τ described above). By inspection, this correlation is maximized when the signals are aligned in time, so by defining an arbitrary but fixed reference time t 0 , the local maximum likelihood (ML) ToA estimate is simply
t ^ R x = τ ^ t 0 ; τ ^ = arg max τ g ( τ ) .
This formulation may seem pedantic for the limited scope of this discussion, but it is important for building the more advanced estimators outlined in [40,45,46]. There are numerous extensions to this simple estimator, including computationally efficient hardware implementations [5], leveraging phase information to adjust the estimate [45], and various iterative refinement methods [40,46].

2.4. Time-of-Flight Estimation

Equations (1) and (2) form a system of 2 equations with 2 unknowns T and τ . If the known transmit timestamps and estimated receive timestamps are shared between users, we have enough information to trivially solve the system. In reality, however, the clocks driving these devices can drift and the platforms may be in motion, so T and τ may change between frames, thereby creating a system of two equations with four unknowns. For reasonable frame rates (≥10 Hz) and modern RF hardware, both T and τ are well approximated by simple first- or second-order Markov models. By injecting these models and collecting data over multiple frames, we can construct closed-form estimators [47] or even Kalman tracking filters [5] to estimate T and τ . Various design implications may inform your choice of algorithm, but fundamentally, they all estimate τ , from which we can infer the distance between platforms d. For multi-antenna systems, we can estimate the distance d i , j between each transmit antenna i and receive antenna j, creating spatial diversity that enables relative position and orientation estimation. This configuration is depicted in Figure 3.

3. Bounds on Position Estimation

In this section, we define the ranging model, derive 2D and 3D CRLBs on position estimation, and discuss geometric dilution of precision (GDoP).

3.1. Ranging Model

Consider N ground antenna labeled g i , which estimates the ToF τ i between a target antenna h, depicted in Figure 3. Assume that the ToF estimates τ i are unbiased estimators of the corresponding distance d i , such that
τ i N ( d i / c , σ τ i 2 ) ,
where c is the speed of light, N denotes the normal distribution, σ τ i 2 denotes the variance of the ToF estimator, and d i is the Euclidean distance between g i and h. The CRLBs for time-of-flight (ToF) estimation are well studied [45] and primarily depend on the signal-to-noise ratio (ISNR) and the signal bandwidth. ISNR is a function of numerous variables, but if the bulk distance between platforms G and H is large compared with the local distance between antennas g i , it is reasonable to assume that each link has comparable ISNR and, therefore, comparable variance, i.e., σ τ i 2 σ τ 2 i { 1 , 2 , , N } . This approximation is reasonable for most applications, but a notable counter-example is any near-field application where the antennas of one system surround the antennas of the other, for example, the Airbus DeckFinder system [48].
The variance of an unbiased estimator α ^ = [ α ^ 1 α ^ 2 α ^ q ] T is bounded by σ α i ^ [ I 1 ( α ) ] i i where I ( α ) is the Fisher information matrix (FIM) of size q × q and is defined [49] as
[ I ( α ) ] i j = E 2 α i α j ln p ( z ; α ) .
We use this definition of the CRLB for the remainder of this manuscript.

3.2. Two-Dimensional CRLB on Position

The probability density function of the ToF distribution summarized in Equation (7) is written explicitly as
p τ ( τ ; h ) = 1 2 π σ τ 2 exp 1 2 σ τ 2 i = 1 N τ d i c 2 ,
where d i = ( x x i ) 2 + ( y y i ) 2 is the Euclidean distance between antennae g i = [ x i , y i ] and h = [ x , y ] . The elements of the FIM I ( h ) are computed using Equation (8) as
[ I ( h ) ] 11 = 1 c 2 σ τ 2 i = 1 N ( x x i ) 2 d i 2 , I ( h ) ] 12 = [ I ( θ ) ] 21 = 1 c 2 σ τ 2 i = 1 N ( x x i ) ( y y i ) d i 2 , I ( h ) ] 22 = 1 c 2 σ τ 2 i = 1 N ( y y i ) 2 d i 2 .
By definition, the CRLB on estimating the 2-D position of h is σ h ^ 2 I 1 ( h ) = adj { I ( h ) } / | I ( h ) | , where the adjugate adj { . } and determinant | . | are written as
adj { I ( h ) } = 1 c 2 σ τ 2 i = 1 N ( y y i ) 2 d i 2 i = 1 N ( x x i ) ( y y i ) d i 2 i = 1 N ( x x i ) ( y y i ) d i 2 i = 1 N ( x x i ) 2 d i 2 ,
| I ( h ) | = 2 ( c 2 σ τ 2 ) 2 i = 1 N j = 1 N i j D 2 2 d i d j 2 ,
where D 2 is the determinant of the distance matrix
D 2 = x x i y y i x x j y y j .
To make this result more accessible to a design engineer, we can loosely visualize the CRLB as a scalar value by defining an overall “position variance” as σ p 2 = σ x 2 + σ y 2 tr { I 1 ( h ) } , which can be written as
σ p 2 c 2 σ τ 2 N 2 i = 1 N j = 1 N i j D 2 2 d i d j 2 1 .

3.3. Three-Dimensional CRLB on Position

Reconsider Equation (9) with the 3D definitions d i = ( x x i ) 2 + ( y y i ) 2 + ( z z i ) 2 , g i = [ x i , y i , z i ] , and h = [ x , y , z ] . The FIM now takes the form
I ( h ) = 1 c 2 σ τ 2 i = 1 N ( x x i ) 2 d i 2 i = 1 N ( x x i ) ( y y i ) d i 2 i = 1 N ( x x i ) ( z z i ) d i 2 i = 1 N ( x x i ) ( y y i ) d i 2 i = 1 N ( y y i ) 2 d i 2 i = 1 N ( y y i ) ( z z i ) d i 2 i = 1 N ( x x i ) ( z z i ) d i 2 i = 1 N ( y y i ) ( z z i ) d i 2 i = 1 N ( z z i ) 2 d i 2 .
The CRLB is again expressed as σ h ^ 2 I 1 ( h ) = adj { I ( h ) } / | I ( h ) | , where the adjugate now takes the form adj { I ( h ) } = 2 ( c 2 σ τ 2 ) 2 A , where
A = i = 1 N j = 1 N V i , j ( y , z ) 2 ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , z ) V i , j ( z , y ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , y ) V i , j ( y , z ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , z ) V i , j ( z , y ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , z ) 2 ( d i d j ) 2 i = 1 N j = 1 N V i , j ( y , x ) V i , j ( x , z ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , y ) V i , j ( y , z ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( y , x ) V i , j ( x , z ) ( d i d j ) 2 i = 1 N j = 1 N V i , j ( x , y ) 2 ( d i d j ) 2 ,
and the determinant takes the form
| I ( h ) | = 6 ( c 2 σ τ 2 ) 3 i = 1 N j = 1 N k = 1 N V i , j , k d i d j d k 2 ,
where V i , j , k is the volume of a tetrahedron formed by points h, g i , g j , and g k , as depicted in Figure 4. When represented as the Cayley–Menger determinant, this takes the form
V i , j , k = 1 6 x x i y y i z z i x x j y y j z z j x x k y y k z z k ,
Thus, the closed-form CRLB is written as
I ( h ) 1 = adj { I ( h ) } | I ( h ) | = c 2 σ τ 2 3 A i = 1 N j = 1 N k = 1 N V i , j , k d i d j d k 2 1
As in Equation (14), we write the “position variance” σ p 2 = σ x 2 + σ y 2 + σ z 2 tr { I 1 ( h ) } as
σ p 2 c 2 σ τ 2 3 i = 1 N j = 1 N V i , j ( x , y ) 2 + V i , j ( y , z ) 2 + V i , j ( x , z ) 2 ( d i d j ) 2 i = 1 N j = 1 N k = 1 N V i , j , k d i d j d k 2
where V i , j ( · , · ) is the area of the triangle in three-dimensional Cartesian space with vertices h, g i , and g j when projected onto the ( · , · ) plane, i.e.,
V i , j ( x , y ) = 1 2 x x i y y i x x j y y j , V i , j ( y , z ) = 1 2 y y i z z i y y j z z j , V i , j ( x , z ) = 1 2 x x i z z i x x j z z j
To the best of our knowledge, this is the first closed-form, three-dimensional CRLB on position estimation in a TWR system.

3.4. Geometric Interpretation of D

The determinant D 2 can be interpreted in several different ways:
  • The shortest (or perpendicular) distance between the target antenna h = [ x , y ] and the line joining ground node antennae g i = [ x i , y i ] and g j = [ x j , y j ] is d i , j = D 2 / d i , j , where d i , j = d ( g i , g j ) (refer to Figure 3). The lower bound then reduces [50] to
    σ p 2 c 2 σ τ 2 N 2 i = 1 N j = 1 N i j d i , j d i , j 2 d i d j 2 1 .
  • As a higher-order generalization, interpret D 2 / 2 ! as the volume of a 2-simplex, V i , j , formed by the antennae h, g i , and g j (refer to Figure 3). The bound then becomes
    σ p 2 c 2 σ τ 2 N 2 ! i = 1 N j = 1 N i j V i , j d i d j 2 1 .
  • The quantity ( d i , j d i , j / d i d j ) is called geometric conditioning A i , j and is a measure of the area of a parallelogram contained by vectors h g i and h g j , scaled by length of those vectors. In this case, the lower bound becomes
    σ p 2 c 2 σ τ 2 N 2 i = 1 N j = 1 N i j A i , j 2 2 1 ,
    where A i , j / 2 can be envisioned as the area of triangle Δ i , j enclosed by unit vectors in the same direction. It is interesting to note that this geometric conditioning is independent of the absolute distances between the ground and target antennae.

3.5. Geometric Dilution of Precision

By inspection, the CRLB derived in Equation (20) only depends on the ToF estimator performance ( σ τ ) and the geometry of the ground nodes g i . In related estimation problems, the geometric portion of this equation is often referred to as “geometric dilution of precision”, and effectively represents the penalty incurred by using a noisy estimate to further estimate other quantities, in this case using range to estimate position. We can rewrite Equation (20) in the form σ p 2 c 2 σ τ 2 GDoP , where geometric dilution of precision (GDoP) is a unit-less quantity defined by
GDoP = N 2 i = 1 N j = 1 N i j Δ i , j 2 = 2 N k = 1 M sin 2 ( γ k ) ,
where Δ i , j = 1 2 sin ( γ i , j ) is the area of a triangle with unit sides and angle γ i , j between them and M = N 2 is the number of unique triangles enclosed by h g i and h g j i , j { 1 , 2 , , N } . By inspection, observe that maximizing M = k = 1 M sin 2 ( γ k ) minimizes the CRLB for a given ranging precision σ τ . This is a convex problem within the range γ k ( 0 , π ] with unique unambiguous solutions { γ 1 , γ 2 , , γ M } that yield the best performance for a given number of ground antennae N. This is important because it directly informs the optimal antenna placement for a fixed number of antennas N.
In Figure 5, the positioning CRLB, from Equation (20) is plotted as a function of the ranging precision σ τ and the geometric dilution of precision from Equation (24). By reformulating the bound as a function of traditional metrics (ranging precision and GDoP) rather than abstract volumes, this formulation provides a visual reference of the trade space between the quantities that a design engineer actually cares about, namely, the ranging precision, the positioning precision, and the GDoP induced by the platform geometry.
The ranging precision σ τ is primarily governed by the choice of ranging technology. If the design constraints limit the available locations of the antennas on an aircraft, then this curve can help inform decisions about which ranging technology is most appropriate to achieve a target positioning precision. In contrast, if the ranging precision is fixed either by design or operating conditions, this curve can instead inform decisions about the optimal placement of antennas on an aircraft or ground installation. Additional insight in this regard is also provided in [43]. Naturally, all of this information is contained in the original formulation presented in Equation (20), this figure is simply included as a visual reference that does not require any calculations.
In Table 1, the best achievable GDoP value is enumerated as a function of the number of antennas N between three and eight. This particular positioning method requires at least three antennas to converge on a solution and experiences diminishing returns beyond eight. In Table 2, we include a qualitative description of GDoP values summarized from [51]. In many TWR and other positioning applications—particularly GPS—it is common to summarize the GDoP of an installation without expressly stating the actual geometry, so this information is included as a coarse reference.

4. Bounds on Orientation Estimation

In this section, we derive 2-dimensional (2D) and 3-dimensional (3D) CRLBs for an orientation estimator based on the position estimation described process in the previous section.

4.1. Model

In the previous section, we derived bounds on the estimation of the position of antennas h i with respect to reference antennas g i . For N target antennas, we form M = N 2 vectors labeled h i h j ¯ i , j { 1 , 2 , , N } and i j . We now use these estimated positions to estimate the relative orientation of platform H with respect to platform G by observing the rotation of the vectors h i h j ¯ relative to g i g j ¯ . Consider the normalized vectors
b ¯ k = h i h j ¯ / h i h j ¯ 1 / 2 , v ¯ k = g i g j ¯ / g i g j ¯ 1 / 2 , k { 1 , 2 , , M }
We assume [52] that the unit vectors b ¯ k follow a von Mises distribution with mean direction μ ¯ k and concentration κ k , i.e.,
b ¯ k M ( μ ¯ k , κ k ) ; μ ¯ k = R v ¯ k ,
where κ 0 , μ ¯ = 1 , and R is a proper rotation matrix. For κ 1 the von Mises distribution reduces to a uniform distribution in which b ¯ k are uniformly distributed on a circle irrespective of the rotation matrix R . For all other scenarios, if the variance in estimating position h i and h j is σ h i and σ h j , we consider the circular variance in b ¯ k as σ b ¯ k 2 = σ h i 2 + σ h j 2 ; hence, the mean resultant length ρ k of the wrapped distribution and the concentration factor κ k are given by
ρ k = e σ h i h j 2 / 2 ; κ k = A 1 ( ρ k )
where A ( · ) is a ratio of modified Bessel functions.

4.2. Two-Dimensional CRLB on Orientation

The probability density function [52] of 2D vectors b ¯ = [ b ¯ 1 , b ¯ 2 , , b ¯ M ] with mean μ k = R ( θ ) v ¯ k , relative orientation θ , proper 2D rotation matrix R ( θ ) , ground node vectors v ¯ = [ v ¯ 1 , v ¯ 2 , , v ¯ M ] , and concentration κ = [ κ 1 , κ 2 , , κ M ] is
p b ¯ ( b ¯ ; v ¯ , θ , κ ) = k = 1 M 1 2 π I 0 ( κ k ) exp k = 1 M κ k b ¯ k T R ( θ ) v ¯ k
where the normalizing constant I 0 ( · ) is a modified Bessel function of the first kind at order 0. The FIM I ( θ ) is computed as
I ( θ ) = E 2 θ 2 ln f ( ν ¯ ; v ¯ , θ , κ ) = k = 1 M κ k b ¯ k T R ( θ ) v ¯ k ,
since R ( θ ) = R ( θ ) . Thus, the CRLB on an unbiased orientation estimator θ ^ is
σ θ 2 I 1 ( θ ) = 1 k = 1 M κ k b ¯ k T R ( θ ) v ¯ k
It is important to note that the lower bound is independent of the rotation angle θ and only depends on the concentration factor κ and the number of observations M.

4.3. Three-Dimensional CRLB on Orientation

Assume the 3D observations b ¯ k follow a von Mises distribution [52] whose probability density function is given by
p b ¯ ( b ¯ ; v ¯ , θ , κ ) = k = 1 M C 3 ( κ k ) exp k = 1 M κ k b ¯ k T R ( θ ) v ¯ k
where concentrations κ k 0 , means μ ¯ k = 1 , μ ¯ k = R ( θ ) v ¯ k , and C 3 ( κ ) is a normalized constant
C 3 ( κ ) = κ 4 π sinh κ = κ 4 π ( e κ e κ )
The relative rotation between two 2D vectors can be sufficiently defined with a single angle, whereas at least three are required to do the same in a 3D space. For convenience, we chose Tait Bryan angles—roll ( α ), pitch ( β ) and yaw ( γ )—to represent the relative attitude of the target H with respect to the observing node G. Therefore, the rotation matrix R ( θ ) can be written as the product of three independent rotations by α , β , γ about x , y , z axis, i.e., R ( θ ) = R z ( γ ) R y ( β ) R x ( α ) where θ = [ α , β , γ ] , which becomes
R ( θ ) = cos ( β ) cos ( γ ) cos ( α ) sin ( γ ) + sin ( α ) sin ( β ) cos ( γ ) sin ( α ) sin ( γ ) + cos ( α ) sin ( β ) cos ( γ ) cos ( β ) sin ( γ ) cos ( α ) cos ( γ ) + sin ( α ) sin ( β ) sin ( γ ) sin ( α ) cos ( γ ) + cos ( α ) sin ( β ) sin ( γ ) sin ( β ) sin ( α ) cos ( β ) cos ( α ) cos ( β ) .
The FIM I ( θ ) = E 2 θ 2 ln p b ¯ ( b ¯ ; v ¯ , θ , κ ) reduces to
I ( θ ) = i = 1 N κ i b i T 2 α 2 R ( θ ) v ¯ i i = 1 N κ i b i T 2 α β R ( θ ) v ¯ i i = 1 N κ i b i T 2 α γ R ( θ ) v ¯ i i = 1 N κ i b i T 2 α β R ( θ ) v ¯ i i = 1 N κ i b i T 2 β 2 R ( θ ) v ¯ i i = 1 N κ i b i T 2 β γ R ( θ ) v ¯ i i = 1 N κ i b i T 2 α γ R ( θ ) v ¯ i i = 1 N κ i b i T 2 β γ R ( θ ) v ¯ i i = 1 N κ i b i T 2 γ 2 R ( θ ) v ¯ i
The CRLB is then computed by inverting the FIM,
I ( θ ) 1 = adj { I ( h ) } | I ( h ) | .
It is important to note that the structure of the rotation matrix provides simplifications that can be leveraged to derive the CRLB in closed form. The partial derivatives in the rotation matrix significantly complicate the computation of Equation (34), which can be mitigated by the following simplifications. Let S [ a , b , c ] represent a skew-symmetric matrix of the form
S [ a , b , c ] = 0 c b c 0 a b a 0 .
Useful properties of S include S [ 1 , 0 , 0 ] 2 = diag [ 0 , 1 , 1 ] , S [ 0 , 1 , 0 ] 2 = diag [ 1 , 0 , 1 ] , and S [ 0 , 0 , 1 ] 2 = diag [ 1 , 1 , 0 ] . The first-order partial derivatives in the rotation can be simplified using
α R x ( α ) = S [ 1 , 0 , 0 ] R x ( α ) ,
β R y ( β ) = S [ 0 , 1 , 0 ] R y ( β ) ,
γ R z ( γ ) = S [ 0 , 0 , 1 ] R z ( γ ) ,
Therefore,
α R ( θ ) = R z ( γ ) R y ( β ) S [ 1 , 0 , 0 ] R x ( α ) ,
β R ( θ ) = R z ( γ ) S [ 0 , 1 , 0 ] R y ( β ) R x ( α ) ,
γ R ( θ ) = S [ 0 , 0 , 1 ] R z ( γ ) R y ( β ) R x ( α ) ,
and the second-order derivatives of R ( θ ) can be simplified to
2 α 2 R ( θ ) = R z ( γ ) R y ( β ) S [ 1 , 0 , 0 ] 2 R x ( α ) ,
2 β 2 R ( θ ) = R z ( γ ) S [ 0 , 1 , 0 ] 2 R y ( β ) R x ( α ) ,
2 γ 2 R ( θ ) = S [ 0 , 0 , 1 ] 2 R z ( γ ) R y ( β ) R x ( α ) ,
2 α β R ( θ ) = R z ( γ ) S [ 0 , 1 , 0 ] R y ( β ) S [ 1 , 0 , 0 ] R x ( α ) ,
2 β γ R ( θ ) = S [ 0 , 0 , 1 ] R z ( γ ) S [ 0 , 1 , 0 ] R y ( β ) R x ( α ) ,
2 α γ R ( θ ) = S [ 0 , 0 , 1 ] R z ( γ ) R y ( β ) S [ 1 , 0 , 0 ] R x ( α ) .
To the best of our knowledge, this is also the first closed-form, three-dimensional CRLB on orientation estimation in a TWR system.

5. Simulation Results

In this section, we benchmark several common position and orientation estimators against the bounds derived in the previous sections on a MATLAB simulation platform. We also include a very brief discussion of position and orientation estimation for reference.

5.1. Estimation Preliminaries

For the given model, the expected log-likelihood function L ( θ ) = E ln p τ ( t ; h ) reduces to a least-square problem. We adopt the ordinary least squares (LS) formulation from [53] and present it here for the sake of completeness. Define a matrix X and elements of a column vector y for i { 1 , 2 , , N } as
X = 2 ( x 1 x ¯ ) 2 ( y 1 y ¯ ) 2 ( x 2 x ¯ ) 2 ( y 2 y ¯ ) 2 ( x N x ¯ ) 2 ( y N y ¯ ) ; y [ i ] = d ¯ i 2 d i 2 + X g ¯
where g i = [ x i , y i ] T is the location of N ground node antennae, g ¯ = [ x ¯ , y ¯ ] T is their mean location, and d ¯ i = g i g ¯ 1 / 2 is the distance between g i and g ¯ . The unbiased, linearized, least-square estimate of position h is simply
h ^ LS = ( X T X ) 1 X T Y
For the given model, we propose a maximum-likelihood orientation estimator. The expected log-likelihood function L ( θ ) = E ln p b ¯ ( b ¯ ; v ¯ , θ , κ ) is
L ( θ ) = k = 1 M κ k b ¯ k T R ( θ ) v ¯ k + constants .
Since this is a convex function, there exists a unique, unambiguous angle θ [ 0 , π ] for which L ( θ ) is maximized, which is computed as
d d θ L ( θ ) = k = 1 N κ k b ¯ k T S ( θ ) v ¯ k = 0 ,
where S ( θ ) = d d θ R ( θ ) . Therefore, the maximum-likelihood estimate of orientation is obtained by solving Equation (51) as
θ ^ M L = arctan k = 1 M κ k b ¯ k T R ( π 2 ) v ¯ k k = 1 M κ k b ¯ k T v ¯ k
It is well established that κ = A 1 ( ρ ) and, for a two-dimensional space, A ( · ) = I 1 ( · ) / I 0 ( · ) is a ratio of modified Bessel functions of the first kind and order one and zero; there are multiple approximations to A 1 ( · ) [54], and we adopt the one provided in [55]. We attain a comprehensive measure for concentration factor κ by assuming κ k = κ for k { 1 , 2 , , M } , σ h i = σ p for i { 1 , 2 , , N } , and approximating κ as
κ ^ = ρ ¯ ( 2 ρ ¯ 2 ) 1 ρ ¯ 2 ; ρ ¯ = e σ p 2 ,
where 0 ρ ¯ 1 is the mean resultant length. When the standard deviation in position estimation σ p > 1 meter, the concentration factor κ 1 . In those scenarios, the von Mises model assumption fails, and the vectors b ¯ k do not contain any information regarding rotations θ .

5.2. Position—3D CRLB

We implemented three position estimators in our MATLAB simulation platform: ordinary least squares (OLS); iteratively reweighted least squares (IRLS); and non-linear least squares (NLLS). We distributed N = 50 ground antennas randomly in a 10-meter sphere centered at the origin and placed the target antenna h at (0, 0, 50) meters. We swept the distance estimator standard deviation σ d between 10 6 and 10 m, running 1000 Monte Carlo trials at each test point. These results are depicted in Figure 6. Each of the three estimators closely approaches the CRLB in the region of interest, with the NLLS estimator achieving slightly better performance than the IRLS and OLS estimators. This result demonstrates that the proposed CRLB is tractable for a massive number of sensing elements ( 50 ) and that it agrees with existing well-known estimators.
We characterize the positioning performance as a function of GDoP by scaling the ground constellation to meet different integer GDoP values and plotting the CRLB at each value for the same sweep over σ d . These results are depicted in Figure 7. While this plot is not particularly interesting, it can directly inform system design decisions; if you have a fixed-ranging precision, you can determine the necessary antenna placement requirements (GDoP) to achieve a certain positioning performance. Likewise, if your antenna placements are fixed, then you can determine how precise your ranging estimates need to be to achieve a certain positioning performance.

5.3. Orientation—3D CRLB

We implemented four orientation estimators in our MATLAB simulation platform: Tri-Axial Attitude Determination (TRIAD) [56]; Davenport’s Q-method [57,58]; QUaternion ESTimator (QUEST) [59]; and the optimal linear attitude estimator (OLAE) [60]. We placed two ground nodes g 1 and g 2 at (0,0,1) and (0,1,0) meters and defined the normalized vector v 0 between nodes G and H. We then generated normalized antenna position vectors b 1 and b 2 by rotating v k using a direction cosine matrix (DCM) with angles α = β = γ = π 6 . We evaluated κ using Equation (53).
In Figure 8, we plot the performance of these estimators versus the 3D orientation CRLB as a function of the positioning precision σ p . There is some slight variation between the estimators in different regimes, but they are all consistent with the proposed bound. In Figure 9, we substitute κ for σ p following Equation (53) to instead plot the performance as a function of the concentration.

6. Conclusions

In this manuscript, we derived novel, closed-form, and tractable 2D and 3D CRLBs for position and orientation estimation in a multi-antenna, two-way ranging system. These bounds are specifically formulated as a function of the number of antennas, platform geometry, and geometric dilution of precision. They are also reformulated using different geometric interpretations and plotted as functions of ranging precision and geometric dilution of precision to make these theoretical results more accessible and intuitive to a system design engineer. The performance manifold in Figure 5 also depicts the trade space between ranging precision, positioning precision, and antenna placement, providing a visual reference that does not require any calculations. The proposed CRLBs were compared with several popular position and orientation estimators in a simple MATLAB simulation platform to demonstrate that the bounds are tight, achievable, and consistent with existing work.
This work offers a closed-form and tractable CRLB for position and orientation estimation, but only for the subset of multi-antenna, two-way ranging systems described in Figure 1. This is the natural conclusion to the original motivation, but possible future extensions include a characterization of different antenna geometries similar to the comparisons in [43] and an extension to an N-dimensional formulation that may enable space–time or other high-dimensional processing techniques.

Author Contributions

Conceptualization, S.S., A.H., and D.W.B.; data curation, S.S.; formal analysis, S.S.; funding acquisition, D.W.B.; investigation, S.S., S.W., and A.H.; methodology, S.S., A.H., and D.W.B.; project administration, D.W.B.; resources, A.H.; software, S.S., and S.W.; supervision, A.H., and D.W.B.; validation, S.S., and S.W.; visualization, S.S., S.W., and A.H.; writing—original draft preparation, S.S., S.W., 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.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

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:
2DTwo-dimensional
3DThree-dimensional
CRLBCramér–Rao lower bound
DCMDirection cosine matrix
FIMFisher information matrix
GDoPGeometric dilution of precision
IRLSIteratively reweighted least squares
ISNRIntegrated signal-to-noise ratio
LSLeast squares
MLMaximum likelihood
NLLSNon-linear least squares
OLAEOptimal linear attitude estimator
OLSOrdinary least squares
PNTPositioning, navigation, and timing
QUESTQUaternion ESTimator
RFRadio frequency
SWaP-C  Size, weight, power, and cost
TDoATime difference of arrival
ToATime of arrival
ToA/RSSTime of arrival/received signal strength
ToFTime of flight
TRIADTri-axial attitude determination
TWRTwo-way ranging
UAMUrban air mobility
UWBUltra-wideband

References

  1. Sobieszczanski-Sobieski, J.; Haftka, R.T. Multidisciplinary aerospace design optimization: Survey of recent developments. Struct. Optim. 1997, 14, 1–23. [Google Scholar] [CrossRef]
  2. Çetin, E.; Cano, A.; Deransy, R.; Tres, S.; Barrado, C. Implementing mitigations for improving societal acceptance of urban air mobility. Drones 2022, 6, 28. [Google Scholar] [CrossRef]
  3. Weisshaar, T.A. Morphing Aircraft Technology-New Shapes for Aircraft Design; Purdue Univ: Lafayette IN, USA, 2006. [Google Scholar]
  4. Han, S.; Gong, Z.; Meng, W.; Li, C.; Gu, X. Future alternative positioning, navigation, and timing techniques: A survey. IEEE Wirel. Commun. 2016, 23, 154–160. [Google Scholar] [CrossRef]
  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. Macnamara, T. Introduction to Antenna Placement and Installation; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  7. Carter, N.J. The past, present and future challenges of aircraft EMC. IEEE Electromagn. Compat. Mag. 2012, 1, 75–78. [Google Scholar] [CrossRef]
  8. Pahlavan, K.; Krishnamurthy, P.; Geng, Y. Localization challenges for the emergence of the smart world. IEEE Access 2015, 3, 3058–3067. [Google Scholar] [CrossRef]
  9. Skolnik, M.I. Introduction to radar. Radar Handb. 1962, 2, 21. [Google Scholar]
  10. Barton, D.K. Modern Radar System Analysis; Artech House: Norwood, MA, USA, 1988. [Google Scholar]
  11. Levanon, N. Radar Principles; Wiley-Interscience: New York, NY, USA, 1988. [Google Scholar]
  12. Jiang, Y.; Leung, V.C. An asymmetric double sided two-way ranging for crystal offset. In Proceedings of the 2007 International Symposium on Signals, Systems and Electronics, Montreal, QC, Canada, 30 July–2 August 2007; pp. 525–528. [Google Scholar]
  13. James, R.J. A history of radar. IEE Rev. 1989, 35, 343–349. [Google Scholar] [CrossRef]
  14. Maleki, L.; Byrd, J.; Savchenkov, A.; Liang, W.; Ilchenko, V.; Seidel, D.; Matsko, A. On the development of photonic RF oscillators and resonant electro-optic modulators for advanced RF front-end applications. In Proceedings of the 2011 IEEE MTT-S International Microwave Symposium, Baltimore, MD, USA, 5–10 June 2011; pp. 1–4. [Google Scholar]
  15. Kim, H. Double-sided two-way ranging algorithm to reduce ranging time. IEEE Commun. Lett. 2009, 13, 486–488. [Google Scholar] [CrossRef]
  16. Lee, J.X.; Lin, Z.; Shin, F.C.P. Symmetric double side two way ranging with unequal reply time. In Proceedings of the 2007 IEEE 66th Vehicular Technology Conference, Baltimore, MD, USA, 30 September–3 October 2007; pp. 1980–1983. [Google Scholar]
  17. 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]
  18. Neirynck, D.; Luk, E.; McLaughlin, M. An alternative double-sided two-way ranging method. In Proceedings of the 2016 13th Workshop on Positioning, Navigation and Communications (WPNC), Bremen, Germany, 19–20 October 2016; pp. 1–4. [Google Scholar]
  19. 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]
  20. Bliss, D.W.; Govindasamy, S. Adaptive Wireless Communications: MIMO Channels and Networks; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  21. Horváth, K.A.; Ill, G.; Milánkovich, Á. Passive extended double-sided two-way ranging algorithm for UWB positioning. In Proceedings of the 2017 Ninth International Conference on Ubiquitous and Future Networks (ICUFN), Milan, Italy, 4–7 July 2017; pp. 482–487. [Google Scholar]
  22. Laadung, T.; Ulp, S.; Alam, M.M.; Le Moullec, Y. Active-passive two-way ranging using UWB. In Proceedings of the 2020 14th International Conference on Signal Processing and Communication Systems (ICSPCS), Virtual, 14–16 December 2020; pp. 1–5. [Google Scholar]
  23. Shalaby, M.A.; Cossette, C.C.; Forbes, J.R.; Ny, J.L. Calibration and Uncertainty Characterization for Ultra-Wideband Two-Way-Ranging Measurements. arXiv 2022, arXiv:2210.05888. [Google Scholar]
  24. Lange, R.; Seitz, P.; Biber, A.; Schwarte, R. Time-of-flight range imaging with a custom solid state image sensor. In Proceedings of the Laser Metrology and Inspection, SPIE, Balatonfured, Hungary, 24–27 September 1999; Volume 3823, pp. 180–191. [Google Scholar]
  25. Massa, J.S.; Buller, G.S.; Walker, A.C.; Cova, S.; Umasuthan, M.; Wallace, A.M. Time-of-flight optical ranging system based on time-correlated single-photon counting. Appl. Opt. 1998, 37, 7298–7304. [Google Scholar] [CrossRef]
  26. Lange, R.; Seitz, P. Solid-state time-of-flight range camera. IEEE J. Quantum Electron. 2001, 37, 390–397. [Google Scholar] [CrossRef]
  27. Buttgen, B.; Seitz, P. Robust optical time-of-flight range imaging based on smart pixel structures. IEEE Trans. Circuits Syst. I Regul. Pap. 2008, 55, 1512–1525. [Google Scholar] [CrossRef]
  28. Frank, M.; Plaue, M.; Rapp, H.; Köthe, U.; Jähne, B.; Hamprecht, F.A. Theoretical and experimental error analysis of continuous-wave time-of-flight range cameras. Opt. Eng. 2009, 48, 013602. [Google Scholar]
  29. Hong, H.; Quan, R.; Xiang, X.; Xue, W.; Quan, H.; Zhao, W.; Liu, Y.; Cao, M.; Liu, T.; Zhang, S.; et al. Demonstration of 50 km Fiber-optic two-way quantum clock synchronization. J. Light. Technol. 2022, 40, 3723–3728. [Google Scholar] [CrossRef]
  30. Alsindi, N.; Li, X.; Pahlavan, K. Analysis of time of arrival estimation using wideband measurements of indoor radio propagations. IEEE Trans. Instrum. Meas. 2007, 56, 1537–1545. [Google Scholar] [CrossRef]
  31. Dardari, D.; Chong, C.C.; Win, M. Threshold-based time-of-arrival estimators in UWB dense multipath channels. IEEE Trans. Commun. 2008, 56, 1366–1378. [Google Scholar] [CrossRef]
  32. Falsi, C.; Dardari, D.; Mucchi, L.; Win, M.Z. Time of arrival estimation for UWB localizers in realistic environments. EURASIP J. Adv. Signal Process. 2006, 2006, 1–13. [Google Scholar] [CrossRef]
  33. Rabbachin, A.; Oppermann, I.; Denis, B. ML time-of-arrival estimation based on low complexity UWB energy detection. In Proceedings of the 2006 IEEE International Conference on Ultra-Wideband, Waltham, MA, USA, 24–27 September 2006; pp. 599–604. [Google Scholar]
  34. Dardari, D.; Chong, C.C.; Win, M.Z. Improved lower bounds on time-of-arrival estimation error in realistic UWB channels. In Proceedings of the 2006 IEEE International Conference on Ultra-Wideband, Waltham, MA, USA, 24–27 September 2006; pp. 531–537. [Google Scholar]
  35. Wang, P.; Morton, Y.J. Performance comparison of time-of-arrival estimation techniques for LTE signals in realistic multipath propagation channels. NAVIGATION J. Inst. Navig. 2020, 67, 691–712. [Google Scholar] [CrossRef]
  36. Shamaei, K.; Kassas, Z.M. Receiver design and time of arrival estimation for opportunistic localization with 5G signals. IEEE Trans. Wirel. Commun. 2021, 20, 4716–4731. [Google Scholar] [CrossRef]
  37. Voltz, P.J.; Hernandez, D. Maximum likelihood time of arrival estimation for real-time physical location tracking of 802.11 a/g mobile stations in indoor environments. In Proceedings of the PLANS 2004. Position Location and Navigation Symposium (IEEE Cat. No. 04CH37556), Monterey, CA, USA, 26–29 April 2004; pp. 585–591. [Google Scholar]
  38. Xu, C.; Liu, Q. An inertial neural network approach for robust time-of-arrival localization considering clock asynchronization. Neural Netw. 2022, 146, 98–106. [Google Scholar] [CrossRef]
  39. Dardari, D.; Win, M.Z. Ziv-Zakai bound on time-of-arrival estimation with statistical channel knowledge at the receiver. In Proceedings of the 2009 IEEE International Conference on Ultra-Wideband, Vancouver, BC, Canada, 9–11 September 2009; pp. 624–629. [Google Scholar]
  40. Bidigare, P.; Madhow, U.; Mudumbai, R.; Scherber, D. Attaining fundamental bounds on timing synchronization. In Proceedings of the Acoustics, Speech and Signal Processing (ICASSP), Kyoto, Japan, 25–30 March 2012; pp. 5229–5232. [Google Scholar]
  41. Catovic, A.; Sahinoglu, Z. The Cramer-Rao bounds of hybrid TOA/RSS and TDOA/RSS location estimation schemes. IEEE Commun. Lett. 2004, 8, 626–628. [Google Scholar] [CrossRef]
  42. Mailaender, L. On the geolocation bounds for round-trip time-of-arrival and all non-line-of-sight channels. EURASIP J. Adv. Signal Process. 2007, 2008, 1–10. [Google Scholar] [CrossRef]
  43. Yang, B.; Scheuing, J. Cramer-Rao bound and optimum sensor array for source localization from time differences of arrival. In Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, Philadelphia, PA, USA, 18–23 March 2005; Volume 4, pp. 961–964. [Google Scholar]
  44. Balanis, C.A. Antenna Theory: Analysis and Design; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  45. 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]
  46. Bidigare, P.; Oyarzyn, M.; Raeman, D.; Chang, D.; Cousins, D.; O’Donnell, R.; Obranovich, C.; Brown, D.R. Implementation and demonstration of receiver-coordinated distributed transmit beamforming across an ad hoc radio network. In Proceedings of the 2012 Conference Record of the Forty Sixth Asilomar Conference on Signals, Systems and Computers (ASILOMAR), Pacific Grove, CA, USA, 4–7 November 2012; pp. 222–226. [Google Scholar]
  47. 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]
  48. Deckfinder. 2021. Available online: https://www.airbus.com/en/products-services/defence/uas/deckfinder (accessed on 2 February 2023).
  49. Kay, S.M. Fundamentals of Statistical Signal Processing, Volume i: Estimation Theory (v. 1); PTR Prentice-Hall: Englewood Cliffs, NJ, USA, 1993. [Google Scholar]
  50. Patwari, N.; Hero, A.O.; Perkins, M.; Correal, N.S.; O’dea, R.J. Relative location estimation in wireless sensor networks. IEEE Trans. Signal Process. 2003, 51, 2137–2148. [Google Scholar] [CrossRef] [Green Version]
  51. Person, J. Writing Your Own GPS Applications: Part 2. 2008. Available online: https://www.codeproject.com/Articles/9115/Writing-Your-Own-GPS-Applications-Part-2 (accessed on 30 January 2023).
  52. Jupp, P.E.; Mardia, K.V. Directional Statistics; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  53. Navidi, W.; Murphy, W.S., Jr.; Hereman, W. Statistical methods in surveying by trilateration. Comput. Stat. Data Anal. 1998, 27, 209–228. [Google Scholar] [CrossRef]
  54. Sra, S. A short note on parameter approximation for von Mises-Fisher distributions: Furthermore, a fast implementation of I s (x). Comput. Stat. 2012, 27, 177–190. [Google Scholar] [CrossRef]
  55. Banerjee, A.; Merugu, S.; Dhillon, I.S.; Ghosh, J.; Lafferty, J. Clustering with Bregman divergences. J. Mach. Learn. Res. 2005, 6, 1705–1749. [Google Scholar]
  56. Black, H.D. A passive system for determining the attitude of a satellite. AIAA J. 1964, 2, 1350–1351. [Google Scholar] [CrossRef]
  57. Wahba, G. A least squares estimate of satellite attitude. SIAM Rev. 1965, 7, 409. [Google Scholar] [CrossRef]
  58. Davenport, P.B. A Vector Approach to the Algebra of Rotations with Applications; National Aeronautics and Space Administration: Washington, DC, USA, 1968.
  59. Shuster, M.D.; Oh, S.D. Three-axis attitude determination from vector observations. J. Guid. Control 1981, 4, 70–77. [Google Scholar] [CrossRef]
  60. Mortari, D.; Markley, F.L.; Singla, P. Optimal linear attitude estimator. J. Guid. Control Dyn. 2007, 30, 1619–1627. [Google Scholar] [CrossRef]
Figure 1. Example multi-antenna two-way ranging (TWR) system with a 4-antenna ground user A and 4-antenna aerial user B. By measuring the distance between each antenna pair (4 of which are depicted above using different colors to denote unique links), we can estimate the relative position and orientation of each user, enabling a wide variety of distributed airborne applications. We derive bounds on the performance of these estimators as a function of the number of antennas, platform geometry, and geometric dilution of precision (GDoP).
Figure 1. Example multi-antenna two-way ranging (TWR) system with a 4-antenna ground user A and 4-antenna aerial user B. By measuring the distance between each antenna pair (4 of which are depicted above using different colors to denote unique links), we can estimate the relative position and orientation of each user, enabling a wide variety of distributed airborne applications. We derive bounds on the performance of these estimators as a function of the number of antennas, platform geometry, and geometric dilution of precision (GDoP).
Applsci 13 02008 g001
Figure 2. Depiction of two interactions between radios A and B. Interactions are indexed by n. The time offset is labeled T ( · ) and the propagation delay is labeled τ ( · ) . Transmit events are labeled with up arrows, and corresponding receive events are labeled with down arrows. The clocks are misaligned, so the receive timestamps t ( · ) , R x ( · ) follow Equations (1) and (2).
Figure 2. Depiction of two interactions between radios A and B. Interactions are indexed by n. The time offset is labeled T ( · ) and the propagation delay is labeled τ ( · ) . Transmit events are labeled with up arrows, and corresponding receive events are labeled with down arrows. The clocks are misaligned, so the receive timestamps t ( · ) , R x ( · ) follow Equations (1) and (2).
Applsci 13 02008 g002
Figure 3. Example system configuration with a 4-antenna user G and a 4-antenna user H. The distances between each antenna pair are extracted from time-of-flight (ToF) estimates produced by a two-way ranging (TWR) system.
Figure 3. Example system configuration with a 4-antenna user G and a 4-antenna user H. The distances between each antenna pair are extracted from time-of-flight (ToF) estimates produced by a two-way ranging (TWR) system.
Applsci 13 02008 g003
Figure 4. Geometric interpretation of the volume of the tetrahedron V i , j , k ( d i , d j , d k ) , indicated in blue.
Figure 4. Geometric interpretation of the volume of the tetrahedron V i , j , k ( d i , d j , d k ) , indicated in blue.
Applsci 13 02008 g004
Figure 5. Performance manifold defined by Equation (20) as a function of ranging precision and GDoP for N = 4.
Figure 5. Performance manifold defined by Equation (20) as a function of ranging precision and GDoP for N = 4.
Applsci 13 02008 g005
Figure 6. Performance comparison of the OLS, IRLS, and NLLS position estimators vs. the proposed 3D CRLB for N = 50 sensing elements. This demonstrates that the proposed CRLB is tractable for a massive number of elements and closely aligns with well-known results.
Figure 6. Performance comparison of the OLS, IRLS, and NLLS position estimators vs. the proposed 3D CRLB for N = 50 sensing elements. This demonstrates that the proposed CRLB is tractable for a massive number of elements and closely aligns with well-known results.
Applsci 13 02008 g006
Figure 7. The 3D position CRLB vs. ranging precision for several GDoP values. These curves can help inform system design decisions such as antenna placement or required ranging precision. In general, a GDoP of 1 is considered “ideal”, between 2 and 10 is considered “good”, and anything beyond 20 is considered “poor”.
Figure 7. The 3D position CRLB vs. ranging precision for several GDoP values. These curves can help inform system design decisions such as antenna placement or required ranging precision. In general, a GDoP of 1 is considered “ideal”, between 2 and 10 is considered “good”, and anything beyond 20 is considered “poor”.
Applsci 13 02008 g007
Figure 8. The 3D orientation CRLB vs. positioning precision for several attitude estimators.
Figure 8. The 3D orientation CRLB vs. positioning precision for several attitude estimators.
Applsci 13 02008 g008
Figure 9. The 3D orientation CRLB vs. concentration for several attitude estimators.
Figure 9. The 3D orientation CRLB vs. concentration for several attitude estimators.
Applsci 13 02008 g009
Table 1. Best achievable GDoP vs. number of antennas.
Table 1. Best achievable GDoP vs. number of antennas.
N345678
Best GDoP1.1551.00.8940.8160.7560.707
Table 2. Qualitative rating of realistic GDoP values.
Table 2. Qualitative rating of realistic GDoP values.
GDoP12–34–67–89–2021–50
“Rating”“Ideal”ExcellentGoodModerateFairPoor
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.; Welker, S.; Herschfelt, A.; Bliss, D.W. Cramér–Rao Lower Bounds on 3D Position and Orientation Estimation in Distributed Ranging Systems. Appl. Sci. 2023, 13, 2008. https://doi.org/10.3390/app13032008

AMA Style

Srinivas S, Welker S, Herschfelt A, Bliss DW. Cramér–Rao Lower Bounds on 3D Position and Orientation Estimation in Distributed Ranging Systems. Applied Sciences. 2023; 13(3):2008. https://doi.org/10.3390/app13032008

Chicago/Turabian Style

Srinivas, Sharanya, Samuel Welker, Andrew Herschfelt, and Daniel W. Bliss. 2023. "Cramér–Rao Lower Bounds on 3D Position and Orientation Estimation in Distributed Ranging Systems" Applied Sciences 13, no. 3: 2008. https://doi.org/10.3390/app13032008

APA Style

Srinivas, S., Welker, S., Herschfelt, A., & Bliss, D. W. (2023). Cramér–Rao Lower Bounds on 3D Position and Orientation Estimation in Distributed Ranging Systems. Applied Sciences, 13(3), 2008. https://doi.org/10.3390/app13032008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop