Accuracy Bounds for Array-Based Positioning in Dense Multipath Channels †

The accuracy of radio-based positioning systems will be limited by multipath interference in realistic application scenarios. This paper derives closed-form expressions for the Cramér–Rao lower bound (CRLB) on the achievable time-of-arrival (ToA) and angle-of-arrival (AoA) estimation-error variances, considering the presence of multipath radio channels, and extends these results to position estimation. The derivations are based on channel models comprising deterministic, specular multipath components as well as stochastic, diffuse/dense multipath. The derived CRLBs thus allow an evaluation of the influence of channel parameters, the geometric configuration of the environment, and system parameters such as signal bandwidth and array geometry. Our results quantify how the ToA and AoA accuracies decrease when the signal bandwidth is reduced, because more multipath will then interfere with the useful LoS component. Antenna arrays can (partly) compensate this performance loss, exploiting diversity among the multipath interference. For example, the AoA accuracy with a 16-element linear array at 1MHz bandwidth is similar to a two-element array at 1GHz, in the magnitude order of one degree. The ToA accuracy, on the other hand, still scales by a factor of 100 from the cm-regime to the m-regime because of the dominating influence of the signal bandwidth. The position error bound shows the relationship between the range and angle information under realistic indoor channel conditions and their different scaling behaviors as a function of the anchor–agent placement. Specular multipath components have a maximum detrimental influence near the walls. It is shown for an L-shaped room that a fairly even distribution of the position error bound can be achieved throughout the environment, using two anchors equipped with 2×2-array antennas. The accuracy limit due to multipath increases from the 1–10-cm-range at 1GHz bandwidth to the 0.5–1-m-range at 100MHz.


Introduction
High-accuracy positioning will be a key enabler for a wide range of novel applications in different sectors of industry, including manufacturing, logistics, retail, and transportation. Outdoors, global navigation satellite systems have become the backbone of such location-aware systems [1], exploiting the time-of-flight of radio signals, while indoors, a wide range of wireless systems has been considered [2]. The focus on radio technologies is due its advantages, e.g., small size, low cost, and signals in terms CRLB expressions for hybrid ToA-AoA positioning in [23]. The CRLB for positioning including agent movement and showing the contribution of the Doppler shift to the obtainable AoA information was examined in [24]. In [8,25], the influence of unknown array orientation on the resulting positioning accuracy is examined, allowing the inclusion of prior information about position, orientation and multipath parameters. In [26], a polarimetric signal model is used to derive general expressions for the waveform-based CRLB for channel parameter estimation including dense multipath components, but the results were not analyzed with respect to the position estimation problem. Most recently, the performance of massive-MIMO systems equipped with mm-wave antenna arrays at the agent and anchor sides has been studied [27][28][29], evaluating aspects such as the orientation estimation accuracy, the difference between uplink and downlink channels, the potential use of MPCs, and different transceiver frontend configurations.
In this paper, we extend our previous work on performance limits for high-accuracy localization in multipath channels [19,20] towards antenna arrays and joint ToA and AoA estimation by deriving the waveform-based CRLB. We further examine the resulting achievable positioning accuracy to be expected in dense multipath channels, exploiting these ToA and AoA estimates.
We derive a canonical expression for the Fisher information matrix (FIM), which allows a deeper insight into the dependence on the signal and system parameters, e.g., bandwidth, SNR, or number of array elements.
A geometric description of the generic antenna array as employed throughout this work is shown in Figure 1. We stress that arbitrary array geometries can be used for the derivations, as long as the assumptions introduced later on remain valid. Furthermore, the analyzed signal model is very generic and can be adapted to many signaling schemes, as long as the transmitted pulse shape is known. Thus, the results are not limited to UWB signals, but are also applicable to other modulating schemes, including direct-sequence-spread-spectrum and OFDM.
A possible application of these accuracy bounds is to analyze and quantify the positioning performance limits for specific system setups, using a wide range of wireless technologies, e.g., the IEEE 802. 15 positioning including agent movement and showing the contribution of the Doppler shift to the obtainable AoA information. In [25] and [8], the influence of unknown array orientation on the resulting positioning accuracy is examined, allowing an inclusion of prior information about position, orientation and multipath parameters. In [26], a polarimetric signal model has been used to derive Thus, we address the following research questions: • How does dense multipath influence the position estimation based on ToA and AoA estimation?
• How do the system and signal parameters influence the ToA, AoA and position estimation in presence of dense multipath? • How does "path-overlap" by specular components influence the position estimation (based on line-of-sight ToAs and AoAs) in presence of dense multipath?
The main novel aspects we present are: • We formulate the CRLB for array-based AoA estimation in presence of dense multipath (DM) and compare the CRLB for ToA estimation in such channels, providing insight on the impact of system parameters (e.g., bandwidth, antenna configuration, and carrier frequency) and radio channel parameters. (This novel aspect was previously presented at the ICL-GNSS in 2018 [30].) • We analyze and evaluate the position error bound for a multi-anchor scenario based on joint ToA and AoA estimation, showing the trade-off between these two measurement parameters and their scaling behaviors with respect to bandwidth and anchor-agent geometry. • We formulate and analyze "path-overlap", the interference of the useful LoS component by specular MPCs, and illustrate its impact with respect to the environment geometry.
The paper is structured as follows: In Section 2, we introduce the signal model, including the array processing for the LoS and the DM process. In Section 3, we derive the CRLB for ranging and angulation. Next, we numerically evaluate these two bounds, the ranging error bound (REB) and the angulation error bound (AEB) and discuss the influence of the DM process and system parameters on these bounds. In Section 4 we derive the CRLB for positioning called position error bound (PEB). Furthermore, we extend the signal model to include multiple specular components and DM and derive the PEB again to include the effect of path-overlap. Then, we analyze and discuss the PEB for different positioning schemes, including AoA-only, ToA-only and a combined AoA-ToA version. Finally, Section 5 concludes the paper by summing up the findings.

Signal Model
We consider the case of a single agent at an unknown position p = x, y T in Cartesian coordinates which transmits a signal to an anchor at known position q ( ) T . (The theoretical framework described in this paper can be straightforwardly modified to allow for different Tx-Rx configurations.) Similar to the signal models in [19,24], the received signals are modeled as complex baseband equivalent signals, allowing coherent processing of the signals at each antenna element due to the explicit inclusion of the carrier phase. A known wideband signal s(t) ∈ R transmitted from an agent to the anchor yields the received signal at array element m ∈ 1, . . . , M, with M denoting the total number of antennas [20] r ( ) where the first term models the deterministic line-of-sight (LoS) component with complex-valued amplitude a ( ) m ∈ C and time delay τ ( ) m at the antenna elements, and c is the speed of light. The array geometry is depicted in Figure 1 and described in detail in Section 2.1.
whereᾱ ( ) m ∈ R + is the absolute value of the amplitude and the phase is separated into delay-dependent and random phase parts. The latter represents the unknown initial phase of the agent and any influences on the phase due to the antenna responses or other effects. We use j = √ −1 as the imaginary unit. We assume that the array aperture is small enough that the amplitude differences at each antenna element m are negligibly small (Assumption S1), i.e.,ᾱ m ≈ᾱ described by the approximation with α ( ) ∈ C. The second term in Equation (1) accounts for all occurring reflections, representing the dense multipath [19,20]. It is modeled by the convolution of the transmit signal s(t) with a random process ν ( ) m (t). To characterize the process, we use the assumption of uncorrelated scattering (Assumption S2) in the delay domain [31], through which the auto-correlation function of ν ( ) m (t) can be written as In this equation, δ(·) is the Dirac delta, δ[·] is the Kronecker delta, and S ν,m (t) describes the power delay profile (PDP) of the DM at array element m. For simplicity, it is assumed that the DM is uncorrelated between antenna elements m and m (Assumption S3), which is a valid assumption under conditions of λ/2-spaced linear arrays and uniformly distributed scatterers in 3D [32]. It is furthermore assumed that the DM is quasi-stationary in the spatial domain, i.e., the PDP does not change in the vicinity of q ( ) the DM only depends on the agent position p and anchor reference point q ( ) , resulting in an identical PDP shape at each array element (Assumption S4). Similar to Witrisal et al. [19], we model the DM process ν ( ) m (t) by a zero-mean Gaussian process. The last term in Equation (1), w m (t), describes additive white Gaussian noise with double-sided power spectral density of N 0 /2.
In this section, we introduced Assumptions (S1)-(S4). When we invoke one of these assumptions in the following sections, we reference to the specific assumptions.

Relation to Array Geometry
An overview of the array geometry is shown in Figure 1, which is similar to [24]. For simplicity-and in accordance with many practical scenarios-we assume that anchor and agent are located on a plane, removing the need for estimating the elevation angle. The extension to three dimensions in terms of environment as well as arrays is straightforward. For some given array, the time delay, corresponding to the received signals r ( ) m (t) at element m, is related to the array geometry by As shown below, the most practical choice of the reference point is the center of gravity of the antenna element positions (see also [23,25]).
For the sake of a more concise notation, we omit the anchor index in the following, proceeding with a single-anchor scenario, and reintroduce it when needed.

Cramér-Rao Lower Bound for AoA and ToA
A popular measure for evaluating estimator performance is the Cramér-Rao lower bound (CRLB) which gives a lower bound on the achievable variance of any unbiased estimator. Before we examine the position error bound (PEB), we perform a separate evaluation of the ranging error bound (REB) and the angulation error bound (AEB).
The general form of the Fisher information matrix (FIM) for a parameter vector θ that parameterizes the probability density p(r; θ) for an observation r becomes [7,33] from which the CRLB for an estimatorθ of said parameter vector is defined as The operation A B on two matrices indicates that the matrix A − B is positive semidefinite.
The parameter vector for the setup as shown in Figure 1 is θ = φ, τ, Rα, Iα T ∈ R 4 , where Rα and Iα are the real and imaginary parts of the LoS amplitude α ∈ C. Under the Gaussian model, the likelihood function for the received signal vector r m at array element m is where C m = σ 2 n I N + C c,m ∈ R N×N is the overall covariance matrix of the sampled noise processes at element m, I N×N is an N × N identity matrix and σ 2 n = N 0 T s is the noise variance. The covariance of the DM component can be written in matrix notation as C c,m = SS ν,m S H , where S = s 0 ,s 1 , · · · ,s N−1 T is the convolution matrix of the sampled pulsess i = s(iT s ), s((i − 1)T s ), · · · , s((i − N + 1)T s ) T ∈ R N and S ν,m is the correlation matrix of the DM at element m, which is a diagonal matrix under the uncorrelated scattering (US) Assumption (S2). Similar to Witrisal et al. [19], we use an eigendecomposition of the covariance matrix where u i,m and λ i,m ∈ R are the ith eigenvector and corresponding eigenvalue of C c,m which make up the ith column of U m and diagonal element of Λ m . We further introduce a weighted inner product defined as and the induced norm as x 2 H m = x, x H m for a Hilbert space H m characterized by a specific covariance matrix C m (see Appendix in [19]). The weighting can be interpreted as a whitening operation.
A straightforward result from Equations (3) and (4) is that the transmitted signals have equal induced norms, i.e., s m The zero elements of the block matrices are due to the chosen array reference point and the resulting contribution of each array element's FIM [23], which is shortly demonstrated in Appendix A.
To gain insight into the parameters of interest, we use the equivalent Fisher information matrix (EFIM), which is well established in the literature [8,19,24]. As we focus on AoA and ToA estimation, we obtain the corresponding EFIM [J θ ] −1 2×2 for the truncated parameter vector θ 1,2 = φ, τ T by applying the Schur complement on Equation (14) [ The connection between EFIM and CRLB is then The structure of the FIM in Equation (14) indicates that the estimation of the ToA τ and AoA φ of the LoS component in DM are decoupled. The structure of B further shows that the estimation of the LoS amplitude α only influences the estimation of τ. Further insights are examined in the following sections in terms of the achievable bounds on ToA and AoA, as well as the resulting positioning accuracy in terms of the PEB.

Ranging Error Bound (REB)
By evaluating Equation (18), we obtain the ranging error bound (REB) R = J −1 τ , which is the square root of the lower bound var{τ} ≥ J −1 τ of the variance for estimating the time delay of the LoS component.
The result for the REB is very similar to the one shown in [19], the only difference being the scaling by the number of antenna elements M. The EFIM for the range estimation problem is where β 2 is the mean-squared bandwidth [34] defined as for a normalized pulse s 2 T s = 1 andṡ being the sampled derivative of the pulse with respect to τ. The parameter γ τ is a "whitening gain" for the delay estimation, representing the increased ranging information due to the whitening operation introduced in Equations (11) and (12) (cf., [19]), defined as γ τ = β 2 w /β 2 , where β 2 w is the corresponding mean squared bandwidth of the "whitened" pulse, defined as β 2 The factor sin 2 (ξ) ∈ [0, 1] is a loss factor attributable to the estimation of the nuisance parameter α, arising from the distortion of s through the whitening operation (see Appendix A and [19]), where ξ is the angle between s andṡ in the corresponding Hilbert space H. The factors SINR = |α| 2 s 2 H T s /N 0 and SINR τ = γ τ SINR sin 2 (ξ) in Equations (20) and (21) describe the signal-to-interference-plus-noise-ratio (SINR) and the effective SINR for the delay estimation, respectively. For negligible DM, which defaults to the AWGN case, Equation (20) simplifies to the known CRLB for delay estimation [7] and signal-to-noise-ratio SNR = |α| 2 s 2 T s /N 0 . Introducing the DM, it holds that SINR ≤ SNR (see Equation (12) and Section 3.4).

Angulation Error Bound (AEB)
In a similar fashion as for the REB, we obtain the angulation error bound (AEB) A = J −1 φ from Equation (18), which denotes the lower bound of the achievable variance var{φ} ≥ J −1 φ for estimating the AoA of the LoS component. The EFIM element for the angle estimation problem is In Equation (23), c +β 2 is the whitening gain for the AEB, and D m (φ) is a factor determined by the array geometry. In Equation (24), the SINR and the whitening gain for the AEB are combined into the effective SINR for the angle estimation SINR φ . The EFIM for the negligible DM case is given by [28] In contrast to the REB, the AEB can be simplified to obtain a more intuitive expression. The bandwidth term β 2 is negligible due to the usually much higher carrier frequency, i.e., f 2 c β 2 , resulting in γ φ ≈ 1 and further in SINR φ ≈ SINR such that we find the AEB to be proportional to By setting the carrier frequency in relation to the array geometry, we can define a normalized squared array aperture D 2 λ (φ) which captures all the system parameters that influence the AEB. It is defined as Geometric interpretation of the AEB: As becomes obvious when comparing the REB (Equation (20)) with the AEB (Equations (23) and (29)), the latter depends on the direction φ. For a uniform linear array (ULA), the CRLB is larger from end-fire direction, i.e., for signals impinging parallel to the array axis, and lower for sources from broadside direction, i.e., perpendicular to the array axis.
Special case M-ary ULA: For an M-ary ULA with inter-element spacing of λ/2 and resulting d m = (2m − 1)λ/4, oriented parallel to the y-axis, i.e., ψ m = ±π/2, Equation (30) is found to be According to Equation (33), sources from end-fire direction have infinite AEB. As angle estimation is bounded, e.g., for ULAs to A ULA ∈ [− π 2 , π 2 ], the angulation error approaches a uniform distribution on the defined range as the AoA φ approaches end-fire direction.

Simulation Environment
The theoretical bounds for AEB and REB derived in the previous sections were evaluated numerically and validated using simulations. The evaluation was performed for different pulse bandwidths 1/T p of the transmitted signals, where T p is the duration of the transmit pulse. We assumed that the PDP exhibits a double-exponential shape [19,35] with parameters γ rise = 5 ns, γ decay = 20 ns, χ = 1, K LoS = |α| 2 /Ω 1 = 0 dB and SNR = 30 dB. For simplicity, we used an M-ary ULA with inter-element spacing of λ/2 in the simulations, although the theoretical bounds are valid for arbitrary arrays as long as Assumptions (S1)-(S4) are valid. As we linked the inter-element spacing to the carrier frequency f c , we kept f c = 7 GHz fixed.
As transmit pulse s(t), a root-raised-cosine (RRC) pulse was used, as it is often encountered in wireless communications. We set the roll-off factor to β roll−off = 0.6 and used the corresponding minimum sampling rate f s = (1 + β roll−off )/T p needed to avoid aliasing throughout the simulations.

Simulation Results
In Figure 2, the EFIM parameters are depicted as a function of the bandwidth according to Equations (20) and (23), separated into parameters for the REB shown in Figure 2a and for the AEB in Figure 2b. In both subfigures, the SINR tends towards the Rician K-factor of the channel at narrow and to the SNR at high bandwidth (BW), which relates to the fading of the LoS signal component in a DM-channel: At high BW, the LoS is separated from the DM such that no fading occurs, while at low BW the complete DM interferes with the LoS leading to a flat fading channel. This implies that the "detectability" of the LoS component is linked to the SINR. The loss factor sin 2 (ξ) observable in the REB and shown in Figure 2a quantifies the SINR that is lost due to estimating the nuisance parameter α in the presence of the non-white, non-stationary DM process ν m (t) [19].  . REB and AEB parameters over different pulse bandwidths 1/T p : effective SINR for the delay estimation SINR τ , effective SINR for the angle estimation SINR φ , SINR, whitening gain for the delay estimation γ τ , whitening gain for the angle estimation γ φ , and loss factor sin 2 (ξ).
that the PDP exhibits a double-exponential shape [19,35] with parameters γ rise = 5 ns, γ decay = 20 ns, 276 χ = 1, K LoS = |α| 2 /Ω 1 = 0 dB and SNR = 30 dB. For simplicity, we use an M-ary ULA with 277 inter-element spacing of λ/2 in the simulations, although the theoretical bounds are valid for arbitrary 278 arrays as long as assumptions (S1) to (S4) are valid. As we linked the inter-element spacing to the 279 carrier frequency f c , we keep f c = 7 GHz fixed.

280
As transmit pulse s(t) a root-raised-cosine (RRC) pulse was used, as it is often encountered in 281 wireless communications. We set the roll-off factor to β roll−off = 0.6 and used the corresponding 282 minimum sampling rate f s = (1 + β roll−off )/T p needed to avoid aliasing throughout the simulations.

284
In Fig. 2, the EFIM parameters are depicted as a function of the bandwidth according to (20)  the non-white, non-stationary DM process ν m (t) [19]. 293 The whitening gain γ τ quantifies the SINR gain due to the whitening operation, where the known 294 DM statistics are critical to whiten the observed signal, suppressing the DM. Note that it behaves 295 complementary to the SINR curve (see Fig. 2a). The effective SINR for the delay estimation, SINR τ , 296 summarizes these parameters. It shows the same behavior as the SINR for high BW, while for small 297 BW it increases towards the SNR again due to the whitening gain. The effective SINR τ is thus tied to 298 the distortion effect of the DM. At high and low BW, no pulse distortion occurs as either the DM is 299 Figure 2. REB and AEB parameters over different pulse bandwidths 1/T p : effective SINR for the delay estimation SINR τ , effective SINR for the angle estimation SINR φ , SINR, whitening gain for the delay estimation γ τ , whitening gain for the angle estimation γ φ , and loss factor sin 2 (ξ).
The whitening gain γ τ quantifies the SINR gain due to the whitening operation, where the known DM statistics are critical to whiten the observed signal, suppressing the DM. Note that it behaves complementary to the SINR curve (see Figure 2a). The effective SINR for the delay estimation, SINR τ , summarizes these parameters. It shows the same behavior as the SINR for high BW, while, for small BW, it increases towards the SNR, again due to the whitening gain. The effective SINR τ is thus tied to the distortion effect of the DM. At high and low BW, no pulse distortion occurs as either the DM is separated from the LoS component or it interferes with it completely. However, at intermediate BW, the distortion effect reduces the effective SINR for the delay estimation [19].
For the angle estimation problem (see Figure 2b), no significant whitening gain is possible, as the factor γ φ = f 2 c +γ τ β 2 f 2 c +β 2 does not substantially exceed 1. Remember that β/ f c 1 at typical carrier frequencies of radio signals. Furthermore, at typical carrier frequencies, if the fractional BW is high, also the absolute BW is high, meaning that the factor γ τ is small. Thus, the whitening factor γ φ is approximately 1, which in turn ties the effective SINR φ for the angle estimation to the SINR, as assumed in Equation (29).
The results for the REB and the AEB are shown in Figure 3 for

312
In accordance to the results in [19], the MFE starts to deviate from the theoretical bound below In case of lower BW, the MFE again converges to the bound as the MFE starts to make use of the 320 DM and LoS power. It even outperforms the REB, which is somewhat pessimistic at low bandwidth 321 due to the neglected trace-term in the derivation of the FIM (see Appendix A). For a numerical 322 evaluation of this effect, see [19]. The whitening operation of the MLE on the other hand reduces the 323 SNR [19] at low BW, as already discussed.

324
The effect of the DM on the AEB is different to that on the REB: the DM causes the AEB to deviate 325 from the AWGN bound below 1 GHz, but the whitening operation does not influence the AEB. The In accordance with the results in [19], the MFE starts to deviate from the theoretical bound below a BW of 400 MHz, whereas the MLE follows the bound down to 40 MHz for M = 2 and down to 1 MHz for M = 16 (see Figure 3a). The MLE follows the bound longer due to the whitening operation which assumes known AWGN and DM statistics. For an increasing number of antenna elements M, the MLE follows the REB much longer due to the diversity gain, i.e., due to the M-fold increase of the SINR by performing multiple measurements with M antenna array elements. The REB improves correspondingly with 1/ √ M.
In the case of lower BW, the MFE again converges to the bound as the MFE starts to make use of the DM and LoS power. It even outperforms the REB, which is somewhat pessimistic at low bandwidth due to the neglected trace-term in the derivation of the FIM (see Appendix A). For a numerical evaluation of this effect, see [19]. The whitening operation of the MLE on the other hand reduces the SNR [19] at low BW, as already discussed.
The effect of the DM on the AEB is different to that on the REB: the DM causes the AEB to deviate from the AWGN bound below 1 GHz, but the whitening operation does not influence the AEB. The estimators perform similarly (see Figure 3b): The MFE leaves the bound at higher BW and returns to it at lower BW. The MLE follows the AEB down to 40 MHz and 1 MHz corresponding to the detectability of the LoS component within the noise, after the whitening. Again, the diversity gain is evident by the improved convergence for M = 16 in comparison with M = 2. The improvement of AoA accuracy, when comparing the results for M = 2 and M = 16, is attributed to the normalized squared array aperture D 2 λ (φ) as shown in Equation (33), which directly depends on the total number M of λ/2-spaced array elements.
An important observation from the theoretical results is the fact that the performance of both (ToA and AoA) estimators is limited by the influence of DM. Increasing the SNR of the radio signal will not improve the performance, while antenna diversity will help. Furthermore, we would like to emphasize the influence of the DM on the AEB, especially at smaller bandwidths. While the AWGN AEB does not depend on the bandwidth, we show that, by considering the influence of the DM, the AEB depends on the bandwidth (see Figure 3b).

Cramér-Rao Lower Bound for Positioning
Building on the REB and AEB derived in the previous sections, we join the two bounds to gain insight into the position error bound (PEB) which gives the achievable accuracy for estimating the agent position p = x, y T considering one or multiple fixed anchors at positions q ( ) equipped with antenna arrays. We therefore reintroduce the anchor index at the proper quantities. For simplicity of notation, we assume that each anchor is equipped with identical arrays and performs its measurements with the agent in a similar fashion, e.g., with the same bandwidth β and carrier frequency f c , which is a reasonable assumption for most communication systems.

Position Error Bound (PEB)
From the FIM for the parameter vector ψ ( ) = x, y, Rα ( ) , Iα ( ) T , containing the agent position and the LoS amplitudes at each anchor , we acquire the corresponding single-anchor CRLB as The FIM for the parameter vector ψ ( ) can be computed from the FIM for the AoA and ToA as with T ( ) denoting the Jacobian matrix containing the partial derivatives of the parameter vector ψ ( ) with respect to the parameters in θ ( ) [20,24,25] with Each anchor -having a line-of-sight connection to the agent-provides position information about the agents position p by performing independent measurements. Thus, we define the PEB via the EFIM for the position utilizing all anchors from the set of LoS anchors ∈ L 0 as where the EFIM for the position J p is obtained by The matrices A ( ) , B ( ) and D ( ) are found using Equations (15)-(17) for each anchor . By further inserting Equation (18) and the corresponding EFIMs for the AoAs and ToAs from Equations (21) and (29), we obtain the multiple anchor EFIM for the position as The matrices J p,φ and J p,τ are the position EFIMs obtained from AoAs and ToAs, respectively. They also show the available information for either AoA-only or ToA-only positioning. The matrix J r (φ) indicates the direction from which the information arrives at the agent and is usually called ranging direction matrix [8]. It arises from the Jacobian transforming the FIM (or EFIM) for ToA and AoA to the position and has the form where e(φ ( ) ) = cos φ ( ) , sin φ ( ) T is a unit vector pointing in direction of the agent φ ( ) , which is also the sole eigenvector of the ranging direction matrix. This formulation clearly shows that both range-only and angle-only positioning are ill-posed problems for the case of a single LoS anchor, as the EFIM for the position will be rank deficient and thus not invertible. Geometric interpretation of the PEB: The EFIM in Equation (41) further supports the well established result that ranging provides position information in direction of the anchor, whereas angulation is responsible for information perpendicular to the anchor direction (see e.g., [24]). A graphical interpretation illustrating Equation (41) is shown in Figure 4 where we showcase the most important parameters that influence the achievable positioning accuracy. As discussed in Section 3.4, the angulation accuracy is inversely proportional to the SINR, which in turn depends on the signal bandwidth β. The achievable ranging accuracy depends on the bandwidth as well, thus the PEB in a DM channel scales with the bandwidth both in direction and perpendicular to the agent.  Figure 4. PEB as stated in (42) for two anchors q (1) and q (2) equipped with antenna arrays, visualizing the influence of the different parts attributable to ranging (in direction of the agent) and angulation (perpendicular to the agent direction) that contribute to the positioning accuracy of an agent p.
The matrices A ( ) , B ( ) and D ( ) are found using (15), (16) and (17) for each anchor . By fur erting (18) and the corresponding EFIMs for the AoAs and ToAs from (21) and (29) we obtain ltiple anchor EFIM for the position as The matrices J p,φ and J p,τ are the position EFIMs obtained form AoAs and the ToAs respectiv y also show the available information for either AoA-only or ToA-only positioning. The ma ) indicates the direction from which the information arrives at the agent and is usually ca ging direction matrix [8]. It arises from the Jacobian transforming the FIM (or EFIM) for ToA A to the position and has the form

Position Error Bound (PEB) Considering Path-Overlap
The previously described signal model included only a single specular MPC, namely the LoS. Especially for signals of high bandwidth or in indoor scenarios, the assumption of a single specular component is no longer viable. Thus, we adapt the signal model to include multiple specular components k ∈ {1, . . . , K} with ToAs τ k , AoAs φ k and complex-amplitudes α k . This enables us to examine the influence of path-overlap (see, e.g., [18]) between the LoS component and later arriving specular MPCs in dense multipath channels.
For each anchor (we again omit the index in the following equations for readability), we get the extended signal model as where we again assume that Equation (3) and (S1) hold. In comparison to the LOS only signal model, see Equation (4), we redefine the auto-correlation function of the DM process as where the onset of the DM process is the ToA of the LoS signal τ m1 . Without any information about the room geometry, we can still only use the LoS signal for positioning. The parameters of all other specular components are regarded as nuisance parameters in the estimation process; the extended parameter vector is thus θ = φ 1 , . . . , φ K , τ 1 , . . . , τ K , Rα 1 , . . . , Rα K , Iα 1 , . . . , Iα K T . The previous definition of the overall covariance matrix C m of the sampled noise processes at each element m is still valid, thus the likelihood function for r m is with S m = s m1 , . . . , s mk , . . . , s mK , s mk = s(−τ mk ), s(T s − τ mk ), · · · , s((N − 1)T s − τ mk ) T and a m = a m1 , . . . , a mk , . . . , a mK . Noise and DM are still assumed uncorrelated at each element m (S3), therefore the likelihood function p(r; θ) for the stacked observation vector r is again factorized as p(r; θ) = ∏ M m=1 p(r m ; θ). Using the definition of the FIM (Equation (7)), we find J θ with the matrix blocks The elements of the matrix blocks in Equations (47)-(49) are described in detail in Appendix B. By examining the elements of the matrix blocks in Equation (46), we find that path-overlap occurs if the Fourier-weighted inner product between two specular MPCs k and κ is non-zero, i.e., the terms s mk , s mκ H m = 0 in J θ .
To find the FIM for the parameter vector ψ = x, y, Rα 1 , . . . , Rα K , Iα 1 , . . . , Iα K T , we introduce the Jacobian matrix T as with the sub-block Finally, to derive the EFIM for positioning using multiple anchors, we re-introduce the anchor index at all relevant quantities. With the assumption of independent measurements obtained between each anchor and the agent, we get the EFIM as from which we obtain the PEB using Equation (38).
Comparison to single-component PEB: For decreasing path-overlap, the matrices B φ in Equation (48) as well as J τφ and J φτ in Equation (47) will approach zero due to s mk , s mκ H m ≈ 0 ∀ m for k = κ, and in combination with the chosen reference point (e.g., see Appendix A). Furthermore, the matrices J RαRα and J IαIα in Equation (49) and J τRα and J τI α in Equation (48) will become diagonal due to the same mechanism and we will experience subtractive terms leading to the factor sin 2 (ξ), which allow formulation of the effective SINR SINR τ for the LoS component. As we can only use the LoS signal component for positioning, the PEB in Equation (52) will become Equation (41). In contrast, if path-overlap occurs, the subtractive term becomes none-zero, corresponding to a loss of useful position information. The resulting PEB can still be evaluated numerically as demonstrated in the following section. ] T and q (2) = [2, 1] T of the two anchors used for simulations are labeled as A1 ( ) and anchor is equipped with a 2 × 2 uniform rectangular array (URA) oriented parallel to the s with identical element spacings of λ/2 in x-and y-direction. We examine the achievable ccuracy over the whole room. The transmitting agent, equipped with a single antenna, is id-points within the room spaced by 2 cm in x-and y-directions. At each grid-point, the iants of the PEB are evaluated. dance to the examination of REB and AEB, we assume that the statistical properties of the stimated from a number of snapshots taken at the same position and are thus assumed . The DM parameters used are identical to the ones used in Section 3.4.1, i.e. we set γ decay = 20 ns and χ = 1. These are further kept constant over the room to allow a valuation of the PEB. 3 To model the amplitudes of the specular components we assume oefficient of η = 1.6 characterizing an indoor LoS environment [36] and assume that on at a wall results in a further reduction of 3 dB. We compute the amplitude of each e to the LoS amplitude at a distance of 1 m, which we set to α 1 m = 1 without loss of he Rician K-factor K LoS of the LoS power w.r.t. the mean power of the DM was again set to | 2 /Ω 1 = 0 dB, also defined at a distance of 1 m. The SNR is kept constant over the whole = 29.5 dB. The used transmit pulse shape is the same as in  In accordance with the examination of REB and AEB, we assumed that the statistical properties of the DM can be estimated from a number of snapshots taken at the same position and are thus assumed to be known. The DM parameters used are identical to the ones used in Section 3.4.1, i.e., we set γ rise = 5 ns, γ decay = 20 ns and χ = 1. These were kept constant over the room to allow a qualitative evaluation of the PEB. (The authors are aware that this is by no means a realistic assumption as the position of the agent within the environment might influence the DM statistics. It is considered beyond the scope of this work to formulate a realistic simulation model for the behavior of the DM across a room.) To model the amplitudes of the specular components, we assumed a pathloss coefficient of η = 1.6 characterizing an indoor LoS environment [36] and assumed that each reflection at a wall results in a further reduction of 3 dB. We computed the amplitude of each MPC relative to the LoS amplitude at a distance of 1 m, which we set to α 1 m = 1 without loss of generality. The Rician K-factor K LoS of the LoS power with respect to the mean power of the DM was again set to K LoS = |α 1 m | 2 /Ω 1 = 0 dB, also defined at a distance of 1 m. The SNR was kept constant over the whole room at SNR = 29.5 dB. The used transmit pulse shape is the same as in Section 3.4.1, i.e., a RRC-pulse with β roll-off = 0.6 and a sampling rate f s = (1 + β roll−off )/T p .

Simulation Results
We started by evaluating the PEB over the room for different channels: (i) AWGN-only (The AWGN-only channel, meaning that no DM is simulated, is used as a baseline as it has been evaluated in [22][23][24].); (ii) AWGN-plus-DM (cf. Section 4.1); and (iii) AWGN-plus-DM with path-overlap (cf. Section 4.2). Figure 6 shows the PEB over the room alongside the PEB ellipses for the different types of positioning, namely range-only and angle-only as well as the contributions of each anchor individually (Equation (40)) and the resulting full PEB (Equations (41) and (52)).
To examine the bandwidth dependence of the PEB, we performed the simulations at bandwidths of 1 GHz and 100 MHz shown in the left and right columns of Figure 6, respectively. Starting with the case of AWGN-only in Figure 6a,b, we observed a rather homogeneous distribution of the PEB in regions where both anchors are visible. When the agent and the anchors become increasingly co-linear, the accuracy degrades as expected. Furthermore, the well known distance dependent degradation is seen. This is explained by the lower accuracy obtained from the AoA, as its accuracy is inversely proportional to the distance, see, e.g., Equations (41) and (52) and Figure 4. Comparing the PEB for different bandwidths of 1 GHz (see Figure 6a) and 100 MHz (Figure 6b), we see that the range induced accuracy depends on the bandwidth and not on the distance, while the contribution of the AoA is independent of the bandwidth but in turn decreases with distance.
Comparing the cases of AWGN-only (Figure 6a,b) and AWGN-plus-DM (Figure 6c,d), the negative effect of the DM is clearly visible as SINR < SNR (see Figure 2). The PEB of the AWGN-plus-DM case will thus be larger than for the AWGN-only case. The SINR is bandwidth dependent, hence, the influence of the DM is stronger at smaller bandwidth. This was already shown in Figure 3 for the AEB and the REB, which are in turn tightly related to the PEB, e.g., Equation (41).
It is interesting to note that, in contrast to the AWGN-only case, angle estimation becomes bandwidth dependent in DM. This can be seen clearly by comparing the angle-only PEB for the different bandwidths between AWGN-only and AWGN-plus-DM. For the AWGN-only case, the AoA-only PEB is worse than the REB-only PEB, because it is dependent on the LoS signal component by direction and distance. In the AWGN-plus-DM case, this effect is even more severe, because the angle does not profit from the whitening operation, thus, SINR φ = SINR ≤ SINR τ .
The influence of specular components overlapping with the LoS signal component is shown in Figure 6e,f, which evaluates the PEB from Equation (52). At the high bandwidth of 1 GHz, the influence is negligible over most of the room, due to the short pulse duration of 1 ns ≈ 0.3 m. It is only slightly visible in room corners. The effect of path-overlap is much more severe when we reduced the bandwidth to 100 MHz. Due to the now much longer pulses, overlap occurs in regions near the walls, resulting in higher PEB.
These result show that especially in the lower bandwidth region below 1 GHz the AWGN-only assumption results in an overly optimistic PEB, as it neglects both the influences of the DM as well as the possibility of path-overlap caused by specular components. As we usually see specular components apart from the LoS in indoor scenarios, these results again highlight the difficulties due to DM and path-overlap.  Figure 7 examines the contribution of a single anchor in more detail by comparing the PEB for a single-anchor setup using only A1 (see Figure 7b) and comparing it to the full PEB when using both anchors (Figure 7a). The case of AWGN-plus-DM with path-overlap was considered for a bandwidth of 100 MHz. The trivial result, when using only a single anchor, is that, at non-LoS locations, no position can be obtained. From the orientation of the PEB ellipses, we again observed the decreasing accuracy perpendicular to the anchor direction due to the angulation accuracy being inversely proportional to the distance. It should further be noted that both angle and range estimation are influenced by path-overlap, resulting in general in a location-dependent deterioration of the PEB, most severely in room corners. The distance deterioration is rather well compensated by adding a second anchor. performance has been analyzed, considering these two measurement types, multiple anchors, and also 508 the influence of path-overlap by specular multipath components (MPCs).

509
The multipath signal model makes our numerical results significantly more realistic than previous  Comparing the error ellipses in both plots shows that range and angle information are perpendicular to one another. We again observe that multiangulation suffers from the increasing distance and path-overlap, whereas multilateration is independent of the distance between anchor and agent but still subject to path-overlap. Even thoguh the position related information contained in angle and range measurements are perpendicular, both cases suffer greatly from co-linearity of anchors and agent. Joint positioning using angles and ranges can nevertheless obtain high accuracy as shown in performance has been analyzed, considering these two measurement types, multiple anchors, and also 508 the influence of path-overlap by specular multipath components (MPCs).

509
The multipath signal model makes our numerical results significantly more realistic than previous

Conclusions
The ranging and angulation error bounds (REB and AEB, respectively) were evaluated for time-of-arrival (ToA) and angle-of-arrival (AoA) estimation in dense multipath channels, examining the influence of the number and geometry of array elements as well as the effect of the bandwidth of the transmitted signal. Simulation results are presented for a matched filter estimator (MFE) and a maximum likelihood estimator (MLE), validating these results. Furthermore, the positioning performance was analyzed, considering these two measurement types, multiple anchors, and the influence of path-overlap by specular multipath components (MPCs). The multipath signal model makes our numerical results significantly more realistic than previous analytical results. Most notably, it is shown that, due to the impact of dense multipath, the achievable accuracy for AoA estimation becomes bandwidth dependent. This is a novel finding and stands in contrast to the known results for the AWGN-only channel, where the AoA accuracy appears to be widely independent of the employed signal bandwidth. Specifically, the AEB scales with the squared array aperture, the number of array-antenna elements, and the signal-to-interference-plus-noise ratio (SINR). The REB scales with the squared signal bandwidth, the number of antennas, and the SINR. The SINR itself quantifies the influence of the dense multipath. It also scales with the bandwidth, because the interfering multipath can be better resolved at a higher signal bandwidth. The relevance of our work lies in the exact numerical quantification of these different influencing factors, which allows for the evaluation of trade-offs between various system configurations under realistic channel conditions. The analysis of the position error bound (PEB) illustrates the relationship between the ToA and AoA information in multi-anchor positioning scenarios and the influence of the anchor-agent placement. The different measurement types can complement one another. The AoA will be most useful at close range-allowing for accurate single-anchor positioning at close range-while, at far range, only the ranging information remains accurate. Interfering specular multipath components have a maximum detrimental influence near the walls.
Future work may couple the presented theoretical framework with a detailed, parameterized channel model of the application environment to study the expected, site-specific positioning performance as a function of a wide range of system parameters.
Author Contributions: T.W. derived the formal analysis (theoretic results), developed the experiments for validation, and performed the visualization. S.G. helped with the formal analysis, the derivations, and the validation. U.M. provided scientific support and supervision. K.W. helped with the methodology, formal analysis, and derivations, and provided scientific support and supervision. The manuscript was written by T.W. and revised by S.G., and K.W.
Funding: This research was partly funded by the Austrian Research Promotion Agency (FFG) within the project UBSmart (project number: 859475). The financial support by the Christian Doppler Research Association, the Austrian Federal Ministry for Digital and Economic Affairs and the National Foundation for Research, Technology and Development is gratefully acknowledged.

Conflicts of Interest:
The authors declare no conflict of interest.