Purposeful Co-Design of OFDM Signals for Ranging and Communications

This paper analyzes the fundamental trade-offs that occur in the co-design of pilot resource allocations in orthogonal frequency-division multiplexing signals for both ranging (via time-of-arrival estimation) and communications. These trade-offs are quantified through the Shannon capacity bound, probability of outage, and the Ziv-Zakai bound on range estimation variance. Bounds are derived for signals experiencing frequency-selective Rayleigh block fading, accounting for the impact of limited channel knowledge and multi-antenna reception. Uncompensated carrier frequency offset and phase errors are also factored into the capacity bounds. Analysis based on the derived bounds demonstrates how Pareto-optimal design choices can be made to optimize the communication throughput, probability of outage, and ranging variance. Different pilot resource allocation strategies are then analyzed, showing how Pareto-optimal design choices change depending on the channel.

backdrop of these significant limitations and in view of the increasing importance of positioning in wireless networks, the current paper considers the fundamental theoretical trade-offs inherent in co-design of the pilot resource allocations in OFDM waveforms for both communications and positioning.It then proceeds to explore the pilot resource designs that optimize communications or positioning while satisfying a threshold requirement for the other.
Consider a scheme in which a UE attempts to determine its position through TOA estimates of signals arriving from multiple cellular base stations.The simplest approach the UE can take is to correlate its received samples against a local replica signal constructed from the known OFDM pilot resources.Through this correlation, the UE can extract a TOA estimate [6] relative to its local clock.This TOA estimate may then be used as a pseudorange measurement [3], a time-difference-of-arrival (TDOA) measurement if subtracted against measurements from other base stations [7], or as part of a round-trip-time (RTT) measurement if the base station measures the TOA in the uplink direction as well [7].Much like in traditional global navigation satellite system processing, the UE can then determine its location and time relative to the known locations of the base stations and to the network time.However, the requirement to collect TOA estimates from multiple base stations may not be easily satisfied in wireless networks, which are often designed to serve a given region with a single base station.In this setting, the signal-to-noise ratios (SNRs) of signals collected from more distant base stations may be so low that the associated TOA estimates experience a thresholding effect where their variance rises dramatically.This thresholding effect has been thoroughly studied in the analysis of fundamental bounds on parameter estimation, such as the Barankin bound [8,9] and the Ziv-Zakai bound [10], and also in the context of TOA estimation [11][12][13].The thresholding phenomenon occurs when sidelobes (secondary peaks) in the autocorrelation function become the dominant source of estimation error.Figure 1 depicts estimation errors that may occur along an autocorrelation function.At a sufficiently high SNR, the variance of the noise in the receiver's correlation output will be much smaller than the ratio between the peak power of the mainlobe and the peak powers of the sidelobes.As a result, the probability of estimation errors occurring on the sidelobes is negligibly low, and errors are instead concentrated on the mainlobe.In this "mainlobe-dominated regime, " TOA accuracy is determined by the shape of TOA estimate on the mainlobe TOA estimate on a sidelobe Fig. 1 A graphical depiction of an autocorrelation function, marking a time-of-arrival (TOA) estimate on the mainlobe and a TOA estimate on the sidelobe.The mainlobe estimate has a smaller difference in power relative to the peak than the sidelobe estimate.The mainlobe estimate also has a smaller error in delay the autocorrelation function near the mainlobe's peak, with a sharper peak yielding greater accuracy.This mainlobe-dominated regime is where the Cramer-Rao lower bound (CRLB) on the TOA error variance applies.However, at a sufficiently low SNR, the probability of estimates settling on the sidelobes becomes non-negligible, causing drastic increases in TOA error variance, as such TOA estimates are far from the mainlobe peak.The operating conditions under which this occurs will be referred to as the "sidelobe-dominated regime." As the SNR drops still further, TOA error variances continue to rise, eventually plateauing in an "ambiguous regime" where meaningful estimates cannot be obtained from the signal, and the best estimator is that which maximizes the prior belief about the TOA.Additional operating regimes may emerge if the UE uses more advanced techniques than simply correlating against known pilots.One such technique is decision-directed estimation, in which the UE first decodes the previously-unknown data resources, reconstructs the OFDM signal using the data resource estimates along with the known pilot resources, and correlates against this reconstructed signal to obtain an improved TOA estimate [14].Since data resources make up a significant portion of the energy in OFDM transmissions, such a decision-directed estimator enjoys noticeable gains in post-correlation SNR and thus a decrease in TOA estimation error when the bit error rate is sufficiently low.Decision-directed approaches have seen success in channel and Doppler estimation [15,16].A similar maximum-likelihood approach called non-dataaided (NDA) estimation has been used for time delay estimation [17].While decisiondirected estimators have been used for positioning [14], further research is needed to quantify their performance in modern cellular networks.Meanwhile, the current paper's focus will be on schemes where TOA is estimated only through correlation against known elements of an OFDM signal.
This paper explores how the design of OFDM pilot resources impacts both ranging precision and communication capacity in various propagation environments.Throughout this paper, "ranging" is used synonymously with TOA estimation at the receiver relative to its local clock.This TOA estimation precedes and is agnostic to the positioning method used (e.g., pseudorange multilateration, TDOA, RTT, etc.).To quantify ranging precision, it employs the Ziv-Zakai bound [10], which, unlike the CRLB, accurately captures estimator errors across all applicable regimes.To quantify communication capacity, Shannon capacity and probability of outage are computed, factoring in impairments due to multipath, block fading, carrier frequency offset (CFO) estimation error, and common phase error (CPE).The placement and power allocation of pilot resources within an OFDM block involves intricate trade-offs between ranging precision and communication throughput, especially when channel impairments are considered.This paper quantifies these trade-offs and proposes OFDM pilot resource designs that balance performance in ranging, capacity, and outages.

Contributions
The main contributions of this paper are as follows: • A derivation of the Shannon capacity bound and probability of outage for generic OFDM signals, accounting for channel estimation error, intercarrier interference, and common phase errors, together with a derivation of the Ziv-Zakai bound on range estimation variance for the same signals.These bounds also account for Rayleigh fading and multipath channels.• A method of co-designing the placement and power allocation of pilot resources in an OFDM signal to achieve both ranging and communication performance requirements through the use of Pareto curves plotting both the Shannon capacity and probability of outage against the Ziv-Zakai bound on ranging error variance.• An analysis of how different channel impairments, fading models, and multipath affect both communications and ranging performance, including their effect on Pareto-optimal pilot resource allocations.

Prior work
Prior work has studied how OFDM signals can be used for positioning, but much of this work operates only within existing protocols rather than proposing new signal designs.This is a broad field of work, covering several protocols of interest.TOA and ranging estimators for LTE signals have been analyzed in [18][19][20].The CRLB for TOA/range estimation is derived in both [19,20] to evaluate the performance of their estimators, but this bound is inapt for low SNR regimes.Outside of OFDM signals, [21] derives the CRLB for TOA estimation that exploits temporal correlation in fading channels but is similarly inapt at low SNRs.A comparison between OFDM and pseudonoise-based signals in [22] demonstrates that OFDM signals may provide improved time-based range estimation performance.Many publications have also focused on the field of opportunistic positioning and navigation, utilizing signals from LTE [3,4], FM OFDM [23], and mobile TV [24,25].While these studies provide valuable insights into the performance capabilities of such estimation and positioning algorithms, they do not address the design of the signals themselves, instead working within existing protocols.Some work has specifically addressed the design of OFDM signals for ranging.Driusso et al. considered signal design and study how the placement of positioning pilots within the LTE framework affects ranging performance by computing the Ziv-Zakai bound in additive white Gaussian noise (AWGN) [26].However, their bounds do not account for fading effects.Wang et al. provided ranging accuracy bounds for a generic OFDM signal model that includes multipath fading but only compute the CRLB, failing to address the SNR threshold effect [27].The study in [28] proposes a unique OFDM design strategy for selecting a sparse subset of dedicated bands such that TOA estimation can meet given requirements under multipath propagation environments.The multipath signal modeling is rigorous, and the estimation computationally complexity is significantly reduced using the proposed sparse design.However, the criterion used for optimization is the CRLB, which ignores sidelobes and thresholding effects at low SNR.Furthermore, the study does not directly address how such a ranging signal would coexist within an OFDM system also being used for communications.Another optimization technique is proposed in [29], in which pilots are allocated to optimize for both time-delay and channel estimation.Much like the previous paper, the CRLB is used in the optimization criterion, limiting this technique's applicability at low SNR, and communication capacity is not addressed.Karisan et al. take a similar approach [30] where the power allocation across pilots is designed to minimize the range estimation CRLB in the presence of interference.But, like others mentioned, this paper does not address low SNR thresholding effects or the trade-offs that such a design would have with a joint communication system.While not specifically addressing OFDM design for ranging, the impact of OFDM design parameters on sidelobe energy in the signal's autocorrelation function was analyzed in [31].The Ziv-Zakai bound has also been applied to TOA estimation in multipath fading channels from wideband pulse signals when the receiver has prior knowledge of the channel multipath [32,33].A multiband OFDM positioning testbed was developed in [34], demonstrating decimeter-level accuracy through real-world TOA-based positioning and a sparse selection of bands.
The communication capacity of OFDM systems has also been extensively studied in prior work.Goldsmith's textbook on wireless communications thoroughly covers the computation of channel capacity and outage in the presence of fading [35].Yoo and Goldsmith extend this analysis to multiple-input multiple-output (MIMO) channels with channel estimation error [36].Tang et al. analyzed the effect of channel estimation error in the presence of Rayleigh fading [37].Ohno and Giannakis provide analysis on the minimum mean square error (MMSE) channel estimation error in OFDM systems and its impact on channel capacity in block Rayleigh fading [38], using this work to propose optimal pilots to maximize capacity.The study in [39] analyzes the trade-off between ranging and communications using the Ziv-Zakai bound for millimeter wave signals as the power allocation is varied.However, the study does not analyze the impact of pilot placement or consider multipath effects, which are more common at sub-6-GHz frequencies.Zhu et al. [40] used the CRLB to optimize beamforming weights to balance communications rate and sensing, but did not analyze TOA estimation or pilot resource design.While the capacity of these systems alone has been thoroughly analyzed and some trade-offs with positioning have been provided, such analysis has not been combined with a ranging variance analysis under multipath fading environments to illuminate the trade-offs between capacity and ranging precision.The current paper is the first to jointly analyze the impact that OFDM pilot placements and power allocations have on both communications and ranging while considering multipath fading channels, degradation in communications caused by channel, CFO, and CPE estimation errors, and thresholding effects in TOA estimation quantified through the use of the Ziv-Zakai bound.
A previous conference publication [41] by the current paper's authors laid the groundwork for this current paper, providing simplified versions of the Shannon capacity, Ziv-Zakai bound, and analysis on the trade-offs experienced in different propagation environments.The current paper significantly extends the conference paper in the following ways.First, it models carrier frequency offset, intercarrier interference, and common phase errors, quantifies the impact that pilot resource allocations have on their estimation errors, and factors the impact of these errors into capacity.Second, it evaluates more realistic OFDM signals spanning multiple symbols, capturing the timedependent aspects of phase errors in the signal.Third, it analyzes the Ziv-Zakai bound with the optimal detector rather than the suboptimal detector in [41] when the receiver has no channel knowledge.Last, it provides a significantly expanded analysis of the results and further quantifies communications performance through outage probability.

Signal model
A user receives an OFDM signal that has propagated through a finite-impulse-response (FIR) channel in the presence of AWGN.The signal model is discrete time and baseband.The OFDM system operates with K subcarriers and a cyclic prefix length of L c .Let (s (m) [k]) K −1 k=0 be the complex modulation sequence for the mth OFDM symbol, where Assume that for the mth OFDM symbol, pilot resources are placed at subcarrier indices I (m)  pilot and data resources at subcarrier indices I (m)  data such that these sets are disjoint and exhaustive: These are modulated onto subcarriers using an inverse DFT, creating a complex time-domain signal with K K + L c samples per OFDM symbol.Suppose M OFDM symbols are transmitted.Then, the transmitted signal is This transmitted signal has a sampling rate of f s Hz , resulting in a sample period of T s = 1 f s s and a subcarrier spacing of 1 KT s Hz .With the cyclic prefix, each OFDM symbol has a period of T sym K T s s .The first sample of the mth OFDM symbol is w[m K ] .The time-domain signal propagates through a fading channel modeled as an order-L FIR filter with L L + 1 complex taps (h[l]) L l=0 , written in vector form as h .The signal is also subject to noise modeled as AWGN and expressed as the sequence Additionally, the signal experiences phase rotations due to CFO with frequency f cfo , and phase noise is modeled as a sampled Wiener process (φ[n]) M K −1 n=0 .Provisionally assuming zero time delay, the resulting received signal from the transmission of ( 1) is given by Without loss of generality, the sample indices in the CFO term are shifted by L c so that the first sample after the cyclic prefix in the first symbol experiences zero phase shift due to CFO.The cyclic prefix is assumed to be sufficiently long for both the communications and ranging functions, permitting the receiver to obtain the frequency-domain signal by discarding the cyclic prefix samples and taking the DFT of each symbol.For reference, (1) (2) the standard LTE cyclic prefix length of 4.69µhboxs has an equivalent range of 1406 m [3].Define the time-domain samples for the mth symbol after discarding the cyclic prefix as expressed in vector form as v (m) .Let F ∈ C K ×K be a DFT matrix, normalized such that ||F || F = √ K , and let F L represent the first L columns of F .The frequency-domain signal for symbol m is x(m) = F x (m) .The frequency-domain channel coefficients are denoted h = √ K F Lh , and the frequency-domain noise is denoted ṽ(m) = F v (m) .Assuming the process φ[n] is slowly varying and can be approximated as constant over the duration of one OFDM symbol, φ[n] is approximated as CPE through a Wiener process: . Since this process is constant over the duration of a symbol, the phase errors are common to all subcarriers.The intercarrier interference (ICI) coefficient from subcarrier i to k can be expressed as [42] and packaged in matrix form as I ICI ∈ C K ×K .The frequency-domain signal for symbol m after OFDM processing of (2) can then be expressed using ) for l = 0, 1, . . ., L.

Communications capacity
As a first step to exploring OFDM signal design trade-offs for ranging and communications, this section explores how the power, quantity, and placement of pilot resources affect CFO and CPE estimation accuracy, which, in turn, determines the extent to which residual ICI and phase error present in the received signal limit the communications rate.

Phase errors and intercarrier interference
OFDM is particularly susceptible to CFO errors, which cause subcarriers to no longer remain orthogonal.This loss of orthogonality causes ICI and results in a decreased signal-to-interference-and-noise ratio (SINR) in the receiver [43].Additionally, CPE arises due to oscillator errors [44].Both of these impairments can be estimated and corrected through the use of training sequences and pilots.This subsection derives (3) ( expressions for the CFO and CPE estimation errors, the residual ICI after correction, and the residual CPE after correction.
Assuming that the receiver experiences a CFO of f cfo Hz, and that the Wie- ner phase noise process φ[m] has independent Gaussian increments with , the received signal model in ( 4) can be rewritten to sep- arate the desired signal from the interference: and ṽ(m) To simplify expressions, it will be helpful to introduce a few new symbols.Define the average received signal power at subcarrier k during symbol m as , and the normalized fading power as . Define the ICI power on subcarrier k as Treating the signal on each subcarrier as an independent Gaussian sample, let s(m) ∼ CN 0, (m)   s , where (m) . In practice, the symbols s(m) will not be Gaussian and will instead adhere to some constellation.How- ever, ICI can be appropriately modeled as Gaussian by the central limit theorem [45], so only the 2nd moments of the symbols s(m) will be considered in this analysis.Finally, let With these preliminaries, conditioned on the channel coefficients and CFO, the ICI can be expressed in vector form as ṽ(m) ICI,k .Additionally, the gain caused by attenuation due to ICI is defined as P att = |I kk | 2 , which has no dependence on k.The instantaneous SINR can then be written as which takes a similar form to the SINR expressions in [43,46].
To isolate the CFO and CPE terms, the effect of the channel coefficients must be mitigated by exploiting the known pilot modulations.Under the block fading assumption, h[i] remains constant for all m.Therefore, the receiver can divide the received signal by the known pilot modulations for any resources at symbol m and subcarriers k ∈ I (m) pilot .Let the received signal in (6) divided by the known pilot modulation be ( 6) for k ∈ I (m) pilot .Following [47], the additive noise ṽ(m) total [k] can be approximated as phase noise 1 ṽ(m) φ [k] .In vector form, this is expressed as ṽφ derived in Appendix A. The noise vectors from all symbols can be stacked into ṽφ ∼ N 0, ṽφ , where ṽφ is a block-diagonal matrix consisting of submatrices (0) ṽφ , (1)  ṽφ , . . ., (M−1) ṽφ .
Additionally, the CPE process φ[m] can be written in vector form as φ ∼ CN 0, � φ , where φ nm = σ 2 φ min{n, m}.Some pilot resources will contain a follow-on pilot resource placed at the same subcarrier but in a subsequent symbol.For each m and each k ∈ I (m) pilot , let n mk be the nearest subsequent OFDM symbol index also containing a pilot resource at subcarrier k.Symbolically, pilot and n mk < ∞ , a phase difference measurement between a pilot resource and its follow-on pilot resource can be obtained from ( 8) and expressed as The number of symbols between each pilot resource and its follow-on pilot resource is defined as As a result, the phase difference measurement in (9) may be expressed as The phase difference measurements �∠ỹ (m) [k] can be stacked as �∠ỹ (m) and the sym- bol index differences d (m) [k] as d (m) .Each vector has dimension N (m) � × 1 , where � ×M , the measurement matrix for φ .Also stack the row vectors z T mk to create Z (m) ∈ R N (m) � ×MK , the measurement matrix for ṽφ .This results in the vectorized form of (10), Now that the phase differences are expressed in a linear form, linear minimum mean square error (LMMSE) estimates [49] of f cfo and φ can be obtained.Stacking the vector measurements in (11) for all m results in (8) which may be written compactly as �∠ỹ = tf cfo + U φ + Z ṽφ .
Define the stacked parameter vector β [f cfo , φT ] T .Recall that the prior covari- ance of φ is φ .Since the receiver has no prior for f cfo , a diffuse prior is assumed such that f cfo ∼ N (0, σ 2 cfo ) and σ cfo → ∞ .Define weight as the limit of the inverse of the prior covariance of β, is the frequency-domain noise after correction, distributed with the same variance σ 2 v as the noise in (2).Using these updated quantities, the corrected signal after the DFT is Now that the receiver has corrected for CFO and CPE, it can finally estimate the channel.This channel estimation step is also prone to errors which will depend upon the power, quantity, and placement of the pilots.(12)      �∠ỹ (0) �∠ỹ (1)  . . .
weight lim

Channel estimation and capacity
This paper's approach to computing the channel capacity is inspired by [38], which considers channel estimation error present in the LMMSE estimate of the channel coefficients h[l].This estimation error factors into the communication link's effective SINR.

Estimating channel coefficients
For simplicity of expression, the residual symbol-dependent phase errors will be dropped from the signal model in (16) for analyzing channel estimation error, 2 resulting in the simplified model where This signal can be written in vector form as x(m) corr = √ P att hcorr ⊙ s(m) + ṽ(m) total .Assume that the receiver knows L and has statistical knowledge of the channel coefficient distributions.Following [38], a LMMSE channel estimation error variance at subcarrier k is determined and denoted σ 2 which is identical to the expression for SINR in (7).The effective SINR in (18) factors in all of the modeled impairments and will be used when computing channel capacity.

Maximum-ratio combining
The signal modeling and analysis up to this point have only considered single-antenna receivers.Now consider a receiver with N RX receive antennas.Assume that the channel coefficients for each antenna are i.i.d. and that each antenna's signal experiences the same CFO and CPE.The channel coefficients are denoted h i for i ∈ [0, . . ., N RX − 1] .Under these assumptions, the receiver's CFO and CPE estimates can be improved by stacking the phase difference measurements in (12) from each antenna.After estimating and correcting for the CFO and CPE, the channels are independently estimated for each antenna. ( After channel estimation, the receiver can employ maximum-ratio combining (MRC) [35], improving the effective SINR.Denote the individual effective SINR for antenna i at subcarrier k and symbol m as ρ (m) ki for i = 0, 1, . . ., N RX − 1 , using (18).Then, the MRC SINR is ρ

Computing channel capacity
Now that an expression has been derived for the effective SINR that accounts for CFO, CPE, and channel estimation errors, channel capacity, and outage probability can be quantified.Given ρ (m)  k , the SINR on subcarrier k during symbol m, the instantaneous capacity [35] summed over all resources and normalized by the number of samples [38] is However, from ( 16), a residual, symbol-dependent phase rotation remains present in the signal.The effect of uncompensated phase errors on capacity is difficult to analyze without simulating symbol constellations, since its impact on symbol decoding error is fundamentally different from that of AWGN.While some work has studied this capacity from an information theoretic view [50,51], this paper opts for a simpler treatment.Residual CPE is tolerable so long as the phase rotation is small enough that symbols do not get rotated into incorrect decoding regions.If the phase rotation is large enough, however, the symbol error rate may increase rapidly.To remain agnostic to constellations and account for this behavior, a tolerable limit of phase rotation ǫ φ,max is assumed such that if |ǫ φ,total [m]| > ǫ φ,max , symbol m is unused, and no data are communicated.Define the indicator function of this condition as 1 ǫ φ,total [m] .Then, multiplying (20) by

total [m] results in a new expression for the instantaneous capacity
The ergodic channel capacity can then be defined as where the expectation is taken over the channel coefficients h , the residual CFO δ , and the residual phase error ǫ φ [m] .In addition to capacity, the probability of outage is quan- tified, i.e., the probability that the instantaneous capacity drops below a particular rate requirement.Given a rate requirement C min , this is If N RX > 1 , the receiver uses MRC.Then, the instantaneous capacity is (20 the ergodic capacity is and the probability of outage is Equations ( 21) through (26) are key to understanding the channel capacity and outage probability as a function of arbitrarily-placed pilot resources in OFDM-based communications.Such expressions, accounting fully for the signal-dependent estimation errors of CFO, CPE, and channel coefficients, are novel to the best of the authors' knowledge.

Ranging variance
The next step in quantifying the OFDM signal design trade-offs for ranging and communications is to consider the factors that affect the receiver's range estimation variance.Specifically, this section explores how the power, quantity, and placement of pilot resources impact the Ziv-Zakai bound on TOA estimate variance.The Ziv-Zakai bound will be used to bound the system's ranging variance by providing a lower bound on the variance of the receiver's TOA estimate.Denote the true TOA as τ in units of seconds, the TOA estimation error as ǫ τ , and the Ziv-Zakai bound as ZZB.
Then, E[ǫ 2  τ ] ≥ ZZB .To construct the ZZB, consider a binary detection problem where two hypotheses are equally likely: (1) The received signal experienced delay τ , and (2) the received signal experienced delay τ + τ z , where τ z = T s z , and z ∈ R is the additional delay in units of samples.The probability of error of the optimal detector between these two hypotheses is defined as P min (z) .If the receiver has a priori knowledge that the TOA is uniformly distributed in [0, T a ] , then the Ziv-Zakai bound can be written as [10,52] Derivations of P min (z) will now be explored for different types of fading and channel knowledge: (1) AWGN, (2) Rayleigh fading with perfect channel knowledge, (3) Rayleigh fading with perfect channel knowledge and MRC, and (4) Rayleigh fading with no channel knowledge.Throughout this section, it is assumed that the receiver's a priori TOA distribution falls within the duration of the cyclic prefix as visualized in Fig. 2, allowing the receiver to discard the cyclic prefix and perform circular correlation against each OFDM symbol as in [6].

Perfect channel knowledge
Under the first three channel types, the receiver has perfect knowledge of the channel coefficients, the CFO, and the CPE.As a result, no estimation errors are present, lowerbounding the ranging variance for each channel type.
Define the pilot vector . First consider the AWGN channel, where the received pilot signal is x (m) p = b (m) + v (m) .Consider a real- valued representation of the signal where the real and imaginary coefficients are stacked: (25) . Define a circular shift function f z (x) , which circularly shifts the vector x by z samples.Note that z may be a fractional num- ber of samples, in which case sinc interpolation is applied.Without loss of generality, assume τ = 0 .The probability density function (PDF) for the received signal with delay τ is and the PDF for the received signal with delay τ + τ z is The likelihood ratio between the delay τ and delay τ + τ z signals is then defined as and the detection test is log � > 0 .The detection statistic log can be simplified as follows: Without loss of generality, assuming the first hypothesis is true and x (m)  sep = µ (m) + η (m) , the test can be rewritten as Denote the real component of the circular autocorrelation function at delay l samples as a(l) = M−1 m=0 µ (m) T f l (µ (m) ) .Recognize in (32) that M−1 m=0 µ (m) T µ (m) = a(0) is the real component of the circular autocorrelation of b (m) evaluated at a delay of 0 samples, and ( 28) 2 A graphical depiction of the received signal and the signal after cyclic prefix removal.The receiver has a priori knowledge of the starting time of the signal within a duration T a shorter than the cyclic prefix is the real component of the same autocorrelation evalu- ated at a delay of z samples.Also recognize in (32) that M−1 m=0 η (m) T (µ (m) − f z (µ (m) )) is a linear combination of jointly-distributed Gaussian random variables and thus is itself distributed as a zero-mean Gaussian with variance Therefore, the probability of error for the detection test can be rewritten as Under Rayleigh fading with perfect channel knowledge, the receiver can distort its known reference signal to account for the channel prior to correlation.The distorted pilot signal can be described as b(m) = b (m) ⊛ h , where ⊛ is circular convolution.
Stack the real and imaginary coefficients as before, yielding ) be the real component of the circular correlation function at a delay of l samples.Conditioning on h , an identi- cal derivation as in the AWGN case results in After the expectation over channel realizations in (35) is taken, the minimum probability of error can be substituted into (27) to obtain the Ziv-Zakai bound.
Next, the probability of error will be quantified for the case of Rayleigh fading with perfect channel knowledge and MRC.Recall that this case assumes that the receiver has N RX receive antennas which experience independent multipath Rayleigh channels.The receiver distorts the pilot signal by each antenna's known channel and then performs correlation with a coherent combination across antennas.Let h i for i ∈ [0, . . ., N RX − 1] be the channel coefficients at each antenna and b(m) i = b (m) ⊛ h i be the distorted pilot signals at receiver i. Accordingly, the autocorrelation of the distorted pilot signal for receiver i at delay l is ǎi (l).
The received signal at receiver i follows as i and channels h i are independent between receive antennas, the same derivation can be followed by conditioning on all channels h i , resulting in (33)

No channel knowledge
So far, the derivations for the Ziv-Zakai bound have assumed that the receiver performing ranging has perfect channel knowledge.While important for understanding the fundamental limits of ranging precision, such information may not be available to the receiver.A new expression for P min (z) will be derived where the receiver has no chan- nel knowledge, but has statistical knowledge of the fading distribution.Without explicit channel knowledge, a practical receiver would also be susceptible to the errors caused by CFO and CPE.However, these errors will not be treated in this expression for P min (z) , simplifying the analysis and yielding a lower-bound that only captures the effects of unknown channel coefficients.
Let B (m) ∈ C K × L be a Toeplitz matrix whose lth column is equal to b (m) circularly shifted by l samples for l = 0, 1, . . ., L .Following [38], the received pilot signal can be expressed as x (m)  p = B (m) h + v (m) .Additionally, let f z (B) indicate a circular shift of the columns of B by z samples.Without any channel corrections applied, the received signal is zero-mean and has a covariance x , which is a block-matrix whose nth row partition and mth column partition are The real and imaginary coefficients are stacked, creating . The PDF for x sep with delay τ is and the PDF with delay τ + τ z is which uses a new covariance , where x z is a block-matrix whose nth row partition and mth column partition are The optimal detection test between these two hypotheses takes the form This form matches that of [53, eq. 3].The minimum probability of error can then be defined as (37) This quadratic form has a generalized Chi-squared distribution, which has no closedform expression for its cumulative distribution function (CDF).However, this CDF can be computed using Imhof 's method [54] or Monte Carlo methods.

Results
Simulated analysis of two example scenarios is now provided to illustrate the use of the ergodic capacity bounds in ( 22) and ( 25), the outage probabilities in ( 23) and ( 26), and the Ziv-Zakai ranging variance bound in (27).These results highlight the trade-offs between communication capacity and ranging variance as a function of pilot placements, pilot power allocations, and propagation environments.

Simulation setup
The simulation analysis adopts an OFDM signal structure having interspersed pilots in both time and frequency, with all other resource elements allocated for data transmission.The signal has nine symbols in one block with K = 72 subcarriers and a subcarrier spacing of 30 kHz.The carrier frequency is f c = 3.5 GHz , the cyclic prefix is L c = 18 samples, and the sampling rate is f s = 2.16 MHz .This results in a cyclic prefix length of 8.33µs and a total symbol duration of 41.66µs .The receiver has a priori knowledge of the TOA uniformly distributed over a duration equal to the cyclic prefix length T a = 8.33µs .Two variants of pilot placements are explored.In the first variant, visualized in Fig. 3a, pilots are equally spaced in both time and frequency.The spacing over time in units of symbols is p sym , and the spacing in frequency in units of subcarriers is p sc .In the second variant, visualized in Fig. 3b, pilots are instead placed at the upper and lower subcarrier limits, with equally-spaced-pilots only placed in the m = 0 symbol spaced two subcarriers apart to ensure that the channel can be estimated at all subcarriers.The number of pilots placed at each extremity is N p .The first variant will be referred to as the "equally-spaced-pilot signal, " while the second variant will be referred to as the "outermost-pilot signal." The equally-spaced-pilot signal is explored because it minimizes channel estimation error [38].The outer-most-pilot signal is explored because allocating power to the extremities of the band maximizes the mean-squared-bandwidth of the signal and minimizes the CRLB on TOA estimation [12].The equally-spaced-pilots in the first symbol of the outer-most-pilot signal are included to allow accurate channel estimation across all subcarriers.The expected values over channel coefficients (i.e., in (22), ( 25), (35), and ( 36)) and outage probabilities (i.e., in (23) and ( 26)) are computed through Monte Carlo methods.For the Ziv-Zakai bound, the integral in ( 27) is computed using Riemann integration sampled at step sizes of T s 64 s .Note that (22) gives the communication capacity bound in units of bit/s/Hz, and that (27) gives the TOA estimate variance in units of secondssquared.To provide a more intuitive understanding of the Ziv-Zakai bound, the results plot the root-mean-square error (RMSE) in units of meters: RMSE = c √ ZZB , where c is the speed of light in meters-per-second (m/s).
To explore the impact of power allocation between pilot and data resources, define a total pilot power P (m)  pilot and total data power P (m) data .Assuming power is equally distributed among the pilot resources, for k ∈ I (m) pilot .Similarly for the data resources, for k ∈ I (m) data .The total power is equivalent over all symbols, P (m) data + P (m) pilot = P total for m = 0, 1, . . ., M − 1 .Defining the set of symbols containing both pilot and data resources as M {m : |I (m)  pilot | � = 0, |I (m) data | � = 0} , a power allocation α is applied such that P (m) data = αP total and P (m) pilot = (1 − α)P total for m ∈ M .The average SNR used throughout the results is defined as Simulation results consider the three channel types: AWGN, frequency-flat Rayleigh fading, and frequency-selective Rayleigh fading.In both Rayleigh cases, signals are simulated with and without multiple receive antennas and MRC.Each antenna's channel power is scaled by 1  N RX to keep the total received signal power consistent, highlighting gains due to diversity rather than simply increased SINR.Recall that the receiver knows L and has statistical knowledge of the fading distributions.
First, the capacity, outage, and ranging variance bounds are evaluated for the equallyspaced-pilot signal under these fading distributions, demonstrating the impact that fading has on both communications and ranging performance.Ranging variance bounds are also shown for the outer-most-pilot signal for comparison.Second, the trade-offs between communications and ranging are analyzed by plotting capacity and probability of outage against ranging error.These plots show results with different pilot resource placements in frequency and power allocations, providing insight into Pareto-optimal signal design choices.Last, capacity is plotted for varying placements of pilot resources over time to demonstrate the impact of residual phase errors.

Bounds
Capacity and outage bounds are first shown for the equally-spaced-pilot signal, with p sym = 2 and p sc = 8 .Power is evenly distributed across subcarriers, resulting in the pilot resources receiving 1/8 of the total signal power in symbols 1, 3, 5, 7, and 9.The simulation uses an outage capacity of C min = 0.25 bit/s/Hz, a common phase noise vari- ance per sample of σ 2 φ = 0.018 rad 2 , and a phase rotation limit of ǫ φ,max = 15 deg .Fig- ure 4a shows the ergodic capacity under an AWGN channel, Rayleigh channels of orders L = 0 and L = 2 , and Rayleigh channels of the same orders with N RX = 4 and MRC.Fig- ure 4b shows the probability of outage for the same signal and channels.Predictably, the greatest capacity is achieved under the AWGN channel, while Rayleigh fading results in reduced throughput.While increasing the order of the channel only slightly decreases capacity in the Rayleigh channels, the order L = 2 channels exhibit reduced outages at higher SNRs.When MRC is used, the diversity gains result in an increased capacity and notably reduced probability of outage.
Figure 5a shows the Ziv-Zakai bound for the same equally-spaced-pilot signal.In this figure, a "known" channel refers to bounds derived using the error probability in (35) or, when MRC is applied, (36).The "unknown" channel refers to bounds derived using the error probability in (42).The AWGN channel achieves the minimum ranging RMSE.Below − 5-dB SNR, the signal enters the sidelobe-dominated regime where rang- ing errors are exacerbated by misdetections occurring on sidelobes in the signal's ambiguity function.Above − 5-dB SNR, the ranging RMSE flattens out.In the known channel case, Rayleigh fading notably degrades performance relative the AWGN channel.However, the ranging errors approach the AWGN curve as diversity is exploited through both increased channel orders and diversity gains with MRC.The unknown Rayleigh channels result in the greatest degradation in performance at low SNR.Interestingly, the order L = 2 unknown Rayleigh channel performs better than the order L = 0 unknown channel above − 8-dB SNR and better than the order L = 0 known channel above − 1-dB SNR. Figure 5a shows that multipath, diversity, and channel knowledge have a drastic impact on ranging RMSE.At 10-dB SNR, the AWGN channel has an RMSE of 2.52 m, the order L = 2 Rayleigh MRC channel has an RMSE of 2.69 m, the order L = 0 Rayleigh MRC channel has an RMSE of 2.91 m, and the order L = 2 Rayleigh known channel has an RMSE of 3.35 m.Significantly increasing from these values, the order L = 2 Rayleigh unknown channel has an RMSE of 11.92 m, the order L = 0 Rayleigh known channel has an RMSE of 18.12 m, and the order L = 0 Rayleigh unknown channel has an RMSE of 25.02 m.
In comparison, Fig. 5b shows the Ziv-Zakai bound for the outer-most-pilot signal, exhibiting structure similar to Fig. 5a.The mainlobe-dominated regime of this signal achieves smaller errors than the equally-spaced-pilot signal.For example, at 10-dB SNR, the AWGN RMSE is reduced from 2.52 to 1.58 m.The outer-most-pilot signal achieves a lower RMSE than the equally-spaced-pilot signal above 0 dB in the AWGN channel, above 1 dB in the order L = 2 Rayleigh MRC channel, above 4 dB in the order L = 0 Ray- leigh MRC channel, above 9 dB in the order L = 2 Rayleigh known channel, and above 6 dB in the order L = 2 Rayleigh unknown channel.For all SNR below these thresholds and for all other channels, the outer-most-pilot signal performs worse than the equallyspaced-pilot signal, due to the different shapes of each signal's autocorrelation function.The allocation of power to the extremities in the outer-most-pilot signal sharpens the mainlobe of the autocorrelation function at the expense of higher sidelobes, resulting in reduced ranging errors only when the SNR remains high enough to prevent detections on the sidelobes.Thus, this paper's bounds offer an immediate and valuable insight: Allocating pilot resources to the subcarrier extremities reduces ranging errors in high SNR channels when MRC and a priori multipath knowledge can be exploited, but an equally-spaced-pilot allocation achieves better performance at lower SNRs and when fading is frequency-flat.

Pareto curves
Now, the trade-off between ranging and communications performance is analyzed for both signal structures.Pareto curves are plotted between capacity and ranging RMSE as well as outage probability and ranging RMSE.Each simulation uses an SNR of 0 dB, an outage capacity of C min = 0.25 bit/s/Hz, a common phase noise variance per sample of σ 2 φ = 0.018 rad 2 , and a phase rotation limit of ǫ φ,max = 15 deg .The data power fraction α is varied from 0.1 to 0.9.
Figure 6a and b shows the capacity and outage Pareto curves for the equally-spacedpilot signal under an order L = 2 Rayleigh channel.For the ZZB, perfect channel knowl- edge is assumed.The pilot spacing across symbols is fixed at p sym = 2 , and each curve corresponds to a different pilot spacing across subcarriers, p sc .The maximum capacity of this system is achieved at p sc = 12 and α ≈ 0.65 .As more power is allocated to the pilot resources by decreasing α , ranging RMSE is reduced by moving downward along the p sc = 12 curve at the expense of data throughput.However, it becomes advanta- geous to use a pilot spacing of p sc = 8 to continue maximizing capacity if a reduced ranging RMSE is desired.In fact, the capacity-maximizing power allocation for p sc = 8 approximately halves the ranging RMSE compared to p sc = 12 while only sacrificing 0.003 bit/s/Hz.Interestingly, pilot spacings of p sc = 4 and p sc = 2 are never Pareto- optimal for 0.1 ≤ α ≤ 0.9 .Similar results are seen in the outage Pareto curve.One expla- nation for this behavior is that increasing p sc allows more resources to be allocated for communications, increasing capacity so long as power is allocated appropriately to handle channel and phase estimation errors.However, increasing this spacing also changes the shape and power of the sidelobes in the autocorrelation function, potentially increasing ranging errors.
Figure 7a and b shows the capacity and outage Pareto curves for the same simulation setup but using the outer-most-pilot signal.This signal structure results in a different shape for the Pareto-optimal design where N p = 8 is the Pareto-optimal choice when power is allocated appropriately.Along the N p = 8 curve below 0.33 bit/s/Hz, slightly lower ranging RMSEs are achieved compared to the equally-spaced-pilot signal.However, the capacity-maximizing power allocation in this signal structure loses approximately 0.09 bit/s/Hz compared to the maximum capacity of the equallyspaced-pilot signal.Similar results are seen in the probability of outage, where the Fig. 6 Pareto trade-off curves plotting capacity (a) and probability of outage (b) against ranging RMSE for the equally-spaced-pilot signal at 0-dB SNR.Translation upward along each curve corresponds to increasing α , the fraction of power allocated to the data resources when pilots are present.The black dots correspond to values of α that maximize the capacity.Curves are plotted for multiple pilot resource spacings across the subcarriers, p sc .The channel is order L = 2 Rayleigh, and channel coefficients are assumed known for ranging only Fig. 7 As Fig. 6 but for the outer-most-pilot signal outer-most-pilot signal cannot minimize outages as well as the equally-spaced-pilot signal but can achieve slightly reduced ranging RMSE if the increased outages are tolerable.This behavior may be caused by increased channel estimation errors since pilots are not distributed equally across the band in the outer-most-pilot signal.
This same analysis can be repeated with an order L = 2 Rayleigh channel with MRC.As seen in the previous analysis, MRC provides immense benefits to capacity, outage, and ranging RMSE. Figure 8a and b shows the capacity and outage Pareto curves for the equally-spaced-pilot signal with N RX = 4 and MRC.Compared to Fig. 6a, the Pareto-optimal design with MRC achieves greater capacity and reduced ranging RMSE for all power allocations.More notable, however, is the reduction in outages that MRC provides, reducing outage probability down to a minimum of 2.4 ×10 − 3 .MRC is able to achieve the results by decreasing the likelihood of experiencing deep fades and extremely low SNRs by exploiting antenna diversity.
The last Pareto plots in Fig. 9a and b show the capacity and outage for the equallyspaced-pilot signal in an order L = 2 unknown Rayleigh channel.Much like Fig. 6a and b, capacity is maximized on the p sc = 12 curve.However, significant improve- ments in ranging accuracy can be gained if the system instead opts for a pilot resource spacing of p sc = 8 or p sc = 4 , reducing the ranging RMSE at the capacity-maxi- mizing point from 79.0 m to 27.0 m at the expense of 0.003 bit/s/Hz, or to 19.0 m at the expense of 0.110 bit/s/Hz.These poor ranging results are caused by an inability of the receiver to exploit multipath diversity when the channel is unknown, making it susceptible to deep fades that drastically increase ranging errors.
Finally, the impact of pilot resource spacing across symbols is analyzed for the equally-spaced-pilot signal.The spacing p sym is varied between 1, 2, 4, and 8.The channel is order L = 2 Rayleigh and the pilot resource spacing in frequency is p sc = 8 .The relative capacity curves in Fig. 10 show that p sym = 1 maximizes capacity for SNRs below 0 dB, whereas an increased p sym reduces CFO, CPE, and channel estimation accuracy, degrading the throughput.Even at high SNRs where channel estimation errors are small, uncompensated phase errors increase when p sym is too large, reducing capacity because of the phase error requirement imposed in (21).However, p sym = 2 maximizes capacity from 0 dB to 16 dB SNR, and

Conclusion
This paper has demonstrated how OFDM pilot resource allocations can be analyzed and designed for both ranging purposes and communications.Bounds were derived to capture the impact that the placement and power allocation of pilot resources have on communications capacity, outages, and the Ziv-Zakai bound on ranging variance.These bounds highlighted how multipath and receive diversity can reduce outages and improve ranging accuracy.Furthermore, ranging accuracy is significantly impacted by the receiver's existing knowledge of the channel.Pareto-optimal pilot resource designs were explored, revealing that ranging errors can be improved with minimal degradation As users continue to demand precise positioning from communication networks, designers of next-generation protocols will need to rethink basic OFDM design.The results in this paper demonstrate that OFDM signals can be designed in a manner that balances the trade-offs between ranging and communications.

1 . 1 K K − 1 n=0
Defining A [t, U ] and �ṽ φ Z ṽφ Z T , the LMMSE estimate of β takes the form resulting in an estimation error ǫ β = β − β with covarianceǫ β = weight + A T −1 �ṽ φ A −Estimates fcfo off cfo and φ of φ are extracted from β , yielding estimation errors ǫ cfo = fcfo − f cfo and ǫ φ = φ − φ .CFO estimation error nor- malized by the subcarrier spacing is defined as δ KT s ǫ cfo and will be used in the ICI analysis.Element-wise, the CPE estimation errors are denoted ǫ φ [m] = φ[m] − φ[m].Having obtained CFO and CPE estimates, the receiver can correct for the estimated CFO in its time-domain signal in (2) by multiplying x[n] by exp −j2π fcfo T s n , perform DFT processing, and then correct for the estimated CPE in its frequencydomain signal by multiplying the mth symbol by exp −j φ[m] .Let ǫ φ,total [m] 2πδT sym m/(KT s ) + ǫ φ [m] .This correction results in updated ICI coeffi- cients I ik = exp j(2πδn/K ) exp j2π in/K exp −j2π kn/K and updated total interference-and-noise vectors where ṽ ǫ h [k] .Define the variance of the frequency-domain channel coefficients as σ 2 h = E[| hcorr [k]| 2 ] .Treating the channel estimation errors as additive noise, an effective SINR at subcarrier k and symbol m is defined as When σ 2 ǫ h [k] = 0 , (18) simplifies to Graff and Humphreys EURASIP Journal on Advances in Signal Processing (2024) 2024:20

Fig. 4 Fig. 5
Fig.4 The capacity (a) and probability of outage (b) of the equally-spaced-pilot signal as a function of SNR under different fading and channel types

Fig. 8
Fig. 8 Pareto trade-off curves plotting capacity (a) and outage (b) against ranging RMSE for the equally-spaced-pilot signal at 0-dB SNR.The channel is order L = 2 Rayleigh and MRC with N RX = 4 is used

Fig. 9 Fig. 10
Fig. 9 Pareto trade-off curves plotting capacity (a) and outage (b) against ranging RMSE for the equally-spaced-pilot signal at 0-dB SNR in an order L = 2 unknown Rayleigh channel Notation: Column vectors are denoted with lowercase bold, e.g., x .Matrices are denoted with uppercase bold, e.g., X .Scalars are denoted without bold, e.g., x.The ith entry of a vector x is denoted x[i].The Euclidean norm is denoted ||x|| , and the Frobenius norm is denoted ||X|| F .Real transpose is represented by the superscript T and conju- gate transpose by the superscript H. Element-wise vector multiplication is represented by ⊙ .Circular convolution is represented by ⊛ .The floor function is denoted as ⌊.⌋ .The Q-function is denoted as Q(.).Zero-based indexing is used throughout the paper.A superscript ˜ is used throughout this paper to indicate signals in the frequency domain after the receiver's OFDM discrete Fourier transform (DFT) processing.