# An extension of the RiMAX multipath estimation algorithm for ultra-wideband channel modeling

Brecht Hanssens, Emmeric Tanghe, Davy Gaillot, M. Lienard, Claude

Oestges, David Plets, Luc Martens, Wout Joseph

► To cite this version:

Brecht Hanssens, Emmeric Tanghe, Davy Gaillot, M. Lienard, Claude Oestges, et al.. An extension of the RiMAX multipath estimation algorithm for ultra-wideband channel modeling. EURASIP Journal on Wireless Communications and Networking, 2018, 2018 (1), pp.164. 10.1186/s13638-018-1177-3. hal-03185604

HAL Id: hal-03185604

https://hal.science/hal-03185604

Submitted on 24 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# An extension of the RiMAX multipath estimation algorithm for ultra-wideband channel modeling

Brecht Hanssens<fnref n="1" />* , Emmeric Tanghe<fnref n="1" />, Davy P. Gaillot², Martine Liénard², Claude Oestges³, David Plets<fnref n="1" />, Luc Martens<fnref n="1" /> and Wout Joseph<fnref n="1" />

**Abstract**

This work presents an extension of the high-resolution RiMAX multipath estimation algorithm, enabling the analysis of frequency-dependent propagation parameters for ultra-wideband (UWB) channel modeling. Since RiMAX is a narrowband algorithm, it does not account for the frequency-dependency of the radio channel or the environment. As such, the impact of certain materials in which these systems operate can no longer be considered constant with respect to frequency, preventing an accurate estimation of multipath parameters for UWB communication. In order to track both the specular and dense multipath components (SMC and DMC) over frequency, an extension to the RiMAX algorithm was developed that can process UWB measurement data. The advantage of our approach is that geometrical propagation parameters do not appear or disappear from one sub-band onto the next. The UWB-RiMAX algorithm makes it possible to re-evaluate common radio channel parameters for DMC in the wideband scenario, and to extend the well-known deterministic propagation model comprising of SMC alone, towards a more hybrid model containing the stochastic contributions from the DMC's distributed diffuse scattering as well.
Our algorithm was tested with synthetic radio channel models in an indoor environment, which show that our algorithm can match up to 99% of the SMC parameters according to the multipath component distance (MCD) metric and that the DMC reverberation time known from the theory of room electromagnetics can be estimated on average with an error margin of less than 2 ns throughout the UWB frequency band. We also present some preliminary results in an indoor environment, which indicate a strong presence of DMC and thus diffuse scattering. The DMC power represents up to 50% of the total measured power for the lower UWB frequencies and reduces to around 30% for the higher UWB frequencies.

**Keywords:** Ultra-wideband, Multipath propagation, Multipath estimation, RiMAX, Specular multipath components, Dense multipath components, Channel modeling, Channel sounding, Indoor environment

## 1 Introduction

In the last couple of years, the physical view of how the radio channel is composed has undergone certain changes. The radio channel used to be considered as a collection of specular multipath components (SMC) that have well-defined discrete locations in the different radio channel dimensions (such as the spatial-, frequency-, or time-delay domain). These SMC are the propagation paths which are considered to have a significant influence

in the total received power and are comprised of the strongest (specular) reflections. Presently, it is widely accepted that over these dimensions, a part of the radio channel is also continuous, originating mainly from distributed diffuse scattering on electrically small objects [1, 2], which are inherently more present at sub-20 GHz frequencies [3, 4]. These are the so-called dense multipath components (DMC). The ideology behind the DMC is to include all radio channel energy that cannot be associated with the SMC, disregarding whether they originate from reflections, diffractions, etc. This implies that certain (specular) multipath components, from a physical point of view, can be regarded as being part of the DMC

<figure><img src="image_3.png" /></figure>

<footnote>*Correspondence: brecht.hanssens@ugent.be</footnote>

<footnote n="1">¹INTEC - WAVES, Ghent University - imec, Ghent, Belgium</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
process because they cannot be reliably detected. This
could be, e.g., due to the limited accuracy of the radio
channel model used, a too low signal-to-noise ratio (SNR),
or because they lie very close to other multipath compo-
nents in the time-delay or angular domains [5, 6]. Since
the phase of the DMC is inaccessible for estimation, only
their average power can be modeled across the aforemen-
tioned domains. This opposes the SMC contribution to
the channel, where both the power and the phase of each
propagation path are accessible.

The objective of this work is to extend the high-resolution DMC-inclusive RiMAX algorithm [7] from its narrowband channel model to an ultra-wideband (UWB) one, facilitating the analysis of frequency-dependent propagation parameters for UWB channel modeling. UWB is currently standardized in wireless personal area network (WPAN) IEEE 802.15.4a, which thanks to its large bandwidth allows for data rates over 2 Gbit/s over a short distance [8]. UWB systems are characterized by their ability to transmit small pulses with a very low power density (limited to – 41.3 dBm/MHz) in a large frequency band (3.1–10.6 GHz) [9]. This enables such systems to harmlessly operate in frequency bands currently occupied by other applications. Combining UWB technology with a MIMO antenna configuration vastly increases the capacity of the UWB system, allowing for extremely high data rates [10], and the accurate localization of target nodes in wireless networks [11–16].

Several ray-tracing studies [17] and measurements
[18, 19] have already shown that SMC alone is not suf-
ficient to describe the complicated interactions to which
electromagnetic waves are exposed to in a real environ-
ment. This is especially true for high-resolution parameter
estimation in the field of channel sounding, where the
number of resolvable SMC is limited by the overall char-
acteristics of the measurement system, the obtained SNR,
and the design of the multipath estimation algorithm [20].
This led to the innovative introduction of DMC in the
RiMAX multipath estimation algorithm, to account for
parts of the radio channel that cannot be resolved as SMC.
However, in contrast to standardized UWB channel mod-
els such as [21], the inclusion of the DMC contributions in
it are missing, which this work aims to overcome.

Multidimensional channel sounding is a necessary pro-
cess to describe the geometric structure of the multiple-
input multiple-output (MIMO) radio channel in terms of
geometrical parameters such as the geometrical parame-
ters such as [22]. These parameters are needed to deduce
a geometry-based stochastic channel model (GBSCM) for
the evaluation of MIMO transmission systems [23] and
to investigate their exact propagation mechanisms. They
can also be used to develop localization algorithms for
location estimation or tracking. In particular, the double-
directional modeling of the radio channel has attracted

a lot of interest because it gives a better physical insight
into the wave propagation mechanisms in real environ-
ments. The effect of DMC has been investigated on the
angular properties of the radio channel in [5, 6, 20, 24],
the polarization characteristics in [23, 25, 26], the cluster-
ing of multipath components in [1, 27], and the channel
capacity in [28, 29] for outdoor and industrial scenarios,
respectively.

The novelty of our approach is that the newly devel-
oped UWB-RiMAX algorithm allows for the global
(i.e., frequency-wise) estimation of multipath parameters
throughout the UWB frequency band, whilst maintain-
ing the DMC-inclusive behavior of the radio channel. In
contrast with executing the RiMAX algorithm in multiple
UWB sub-bands, the advantage of our algorithm is that
the AoD, AoA, and ToA are kept constant over the entire
UWB bandwidth in the initialization procedure, ensuring
that geometrical propagation parameters do not appear or
disappear from one sub-band onto the next. This makes it
possible to estimate the physically most likely radio chan-
nel parameters for the SMC and DMC components in the
wideband scenario and to extend the well-known deter-
ministic propagation model comprising of SMC alone,
towards a hybrid model which also contains the stochas-
tic contributions due to the DMC's distributed diffuse
scattering.

The structure of this paper is as follows; Section 2 describes related work to our approach, whilst Section 3 explains the applied data model. Section 4 then covers the UWB-RiMAX multipath estimation algorithm, and Section 5 describes the applied evaluation metrics. Finally, Sections 6 and 7 highlight the simulation and measurement results of our algorithm, and Section 8 summarizes this paper with a conclusion and some ideas for future work.

## 2 Related work

It was shown previously that the contribution of the DMC
to the capacity of the channel in MIMO systems is quite
significant, and often larger than that of the SMC. In [29],
an experimental analysis of the DMC was conducted in an
industrial environment at 3 GHz, where the DMC covari-
ance structure of the RiMAX data model was validated.
The authors found that the DMC power accounted for 23
to 70% of the total channel power and found it was more
important than in office environments due to its highly
cluttered and metallic nature. Similar results are reported
in [28] for outdoor environments at 5.2 GHz, which found
DMC power contributions ranging from 10% of the total
channel power, even up to 90%. Currently, only a few
more studies have been conducted regarding the influ-
ence of DMC [5, 6, 26, 30, 31]. It was shown in [32, 33]
and [34] that the DMC increases the level of the recon-
structed eigenvalues, resulting in a better approximation
of the measured eigenvalue structure of the MIMO channel. This indicates that an accurate modeling of the DMC parameters is necessary to prevent the underestimation of the MIMO transmission performance [28].

Naturally, the introduction of DMC in the physical model of the radio channel means that common radio channel parameters have to be re-evaluated for diffuse scattering as well [26]. This includes parameters such as mean delay, delay spread, Ricean K-factor, shadowing, fading, cross polarization and ratios. Recent studies such as [26] found that the SMC and DMC power show, on average, a strong correlation of about 0.90 and 0.95 for line-of-sight (LoS) and obstructed-line-of-sight (OLOs) scenarios, respectively. This implies that the DMC can alternatively be interpreted as the non-coherent superposition of paths with weaker SNR, which still follow the specular power decay as a function of distance [5, 6]. The same study found that for OLOs scenarios, the cross-polar normalized DMC power even exceeds 60% on average, indicating that these channels could be modeled relatively accurately by only considering their DMC characteristics (as it is done in the original room electromagnetics model [35]). Hence, simple DMC models can be used to design more advanced channel models, as e.g. proposed in [36], where a distance-dependent model for the power delay profile (PDP) of in-room radio channels was developed. The PDP in this model assumes an early primary component and a DMC reverberant component responsible for the shaping of the tail in the PDP.

The physical reality of DMC raises the question of how well estimation algorithms which historically did not include DMC in their signal model, such as estimation of signal parameters via rotational invariance techniques (ESPRIT) [37] or space-alternating generalized expectation-maximization (SAGE) [38], can estimate the SMC part of the channel. This was investigated in [39, 40], in which the authors compared both ESPRIT and SAGE to the performance of the DMC-inclusive RiMAX algorithm. The results of this study demonstrated that SMC estimation in the presence of DMC is prone to large errors if the signal model is not accordingly modified to cope with DMC contributions, as is the case with ESPRIT or SAGE. This was also shown theoretically in [41]. Therefore, determining the DMC by simply subtracting the specular part (estimated by ESPRIT or SAGE) from the total channel response is flawed and must be avoided. For a reliable estimation of the SMC and/or DMC parameters, the use of DMC-inclusive algorithms such as RiMAX is highly recommended. Next to that, since the resolution and accuracy of classical signal processing algorithms is limited by the available measurement aperture in the space-frequency-time domain, parametric super-resolution algorithms such as ESPRIT, SAGE, and RiMAX are more suitable to enhance the time-delay resolution.

This is done by fitting an appropriate data model to the measured data, allowing the algorithm to overcome the Fourier limitation of the delay resolution.

In [42, 43], the UWB-SAGE algorithm was proposed, which is an extension to the SAGE channel estimation algorithm for UWB channel modeling. The UWB-SAGE algorithm estimates a certain number of individual propagation paths from the measured data and estimates the AoD, AoA, ToA, and the variation of the amplitude and phase for each path. UWB-SAGE is based on the assumption that the UWB channel can be expressed as the superposition of a certain number of sub-bands, in which the scattering loss and the antenna directivity is sufficiently constant. The log-likelihood of the whole UWB bandwidth is then defined as the sum of the log-likelihoods of its sub-bands. This process reduces the distortion effect of amplitude and phase caused by antennas when the parameters of the incident waves are estimated. However, at the same time, the resolution of time-delay decreases due to the sub-band processing, making it necessary to appropriately choose the bandwidth of the sub-bands in which the total UWB bandwidth is divided. However, this algorithm disregards the influence of DMC on the channel transfer function, and is thus not an appropriate algorithm for channel parameter estimation or channel modeling. In our work, the ideology of the UWB-SAGE algorithm will be incorporated in the DMC-inclusive RiMAX algorithm.

# 3 Channel model

## 3.1 Specular- and dense multipath components

To describe the geometric properties of the electromagnetic waves of the MIMO propagation channel in terms of AoD, AoA, and ToA, multidimensional frequency domain channel sounding must be performed. This can be done with (virtual) MIMO array systems, consisting of $M_T$ and $M_R$ antennas at transmitter (Tx) and receiver (Rx), sampled at $M_f$ frequency points. As such, the total amount of samples can be defined by $M$ as follows:

$$M = M_T M_R M_f. \qquad (1)$$

An observation of the frequency response of a MIMO radio channel $h$ can be modeled as the superposition of a deterministic part $x(\theta_{smc})$ (determined by the SMC parameter set $\theta_{smc}$) and a stochastic part $d(\theta_{dan})$ (diffuse scattering and noise; determined by the DMC and noise (DAN) parameter set $\theta_{dan}$). Both parameter sets will be described later in this section.

$$ \begin{aligned} h &= x(\theta_{smc}) + d(\theta_{dan}) \\ h &\in \mathbb{C}^{M \times 1}. \end{aligned} \qquad (2) $$

The deterministic part $x(\theta_{smc})$ of the data model acts as the first order statistics of the radio channel, so that it can be interpreted as the mean of $h$, whilst the stochastic part $d(\theta_{dan})$ describes the second-order statistics by means of
its covariance matrix $R(\theta_{\text{dan}}) \in \mathbb{C}^{M \times M}$, which will be discussed in Section 5.1.2. A realization of the radio channel $h$ can be considered as a random variable distributed according to a complex multivariate Gaussian distribution $h \sim N_c(x(\theta_{\text{smc}}), R(\theta_{\text{dan}}))$ [44] as follows:

$$p(h) = \frac{1}{\pi^M \det(R)} e^{-(h-x)^H R^{-1} (h-x)}. \quad (3)$$

An estimate of the most likely SMC and DMC and noise parameters can be found by maximizing the likelihood function of Eq. (3). Since this is not a trivial task, estimation frameworks such as the RiMAX algorithm will estimate $\hat{\theta}_{\text{smc}}$ and $\hat{\theta}_{\text{dan}}$ of the deterministic and stochastic arrays, such that they maximize the likelihood of observing the measured frequency response $h$ of the radio channel. The objective is thus to find the parameters $\hat{\theta}_{\text{smc}}$ and $\hat{\theta}_{\text{dan}}$ that maximize the correlation with the measurement data. A maximum likelihood (ML) estimator for the parameters $\theta_{\text{smc}}$ and $\theta_{\text{dan}}$ has been proposed in [7, 22], exploiting the fact that the parameters of the two components of the channel model are asymptotically independent. Therefore, one can decouple the estimation problem into two separate estimation problems. The resulting RiMAX algorithm is iterative and alternates between the maximization of the likelihood function with respect to the parameters $\theta_{\text{smc}}$ and $\theta_{\text{dan}}$. It has an approximately linear computational complexity in the number of propagation paths $P$ and in the number of data samples $M$ [7].

Based on the capability of the RiMAX algorithm to extract both parameter sets from the (virtual) array measurement data, the following structures for the deterministic and stochastic parameter sets ($\theta_{\text{smc}} \in \mathbb{C}^{P \times 4S}$ and $\theta_{\text{dan}} \in \mathbb{R}^{4 \times S}$, respectively) can be adopted:

$$\theta_{\text{smc}} = \begin{bmatrix} \varphi_D^T \\ \varphi_A^T \\ \tau_A^T \\ \gamma^T \end{bmatrix}^T \begin{array}{l} \leftarrow \text{SMC angle of departure [rad]} \\ \leftarrow \text{SMC angle of arrival [rad]} \\ \leftarrow \text{SMC time-delay of arrival [s]} \\ \leftarrow \text{SMC complex amplitude [s]}. \end{array} \quad (4)$$

In (4), $\varphi_D$, $\varphi_A$, $\tau_A$, and $\gamma$ are $P \times S$ matrices, where $P$ is the number of SMCs extracted from the measurement data, and $S$ is the number of sub-bands in which the total UWB bandwidth was partitioned (in analogy with the UWB-SAGE algorithm). As such, $\theta_{\text{smc}}$ is of size $4S \times P$. We will discuss the frequency-dependency of the geometrical parameters of the propagation paths in Section 4.2. Each row in the aforementioned matrices $\varphi_D$, $\varphi_A$, and $\tau_A$ contains the corresponding specular parameter for each of the $p \in \mathcal{P}$ propagation paths ($P = |\mathcal{P}|$ in total) and describes its frequency-dependency in each of the $s \in S$ sub-bands ($S = |S|$ in total). We note that the angular modeling was limited to that of the azimuthal plane, which is acceptable since most measurement campaigns are only

performed with a planar array, with which the estimation of elevation parameters is not possible. The extension of the data model to the elevation domain is straightforward, where the AoD and AoA of the SMC will now have an extra elevation component together with an azimuthal component. We also consider only a single snapshot of the channel, such that the covariance matrix $R(\theta_{\text{dan}})$ is only averaged over one observation of the channel. Whilst it would be more reliable to use several snapshots of the channel in a real measurement environment, we would then have to impose a parametric model to handle the time dependence of the SMC, or assume them to be independent and independent and identically distributed (i.i.d.) across snapshots in time. Since this is out of the scope for the purpose of this paper, we will leave this up for future work.

$$\theta_{\text{dan}} = \begin{bmatrix} \alpha_1 \\ \tau_d \\ \tau_r \\ \alpha_0 \end{bmatrix} \begin{array}{l} \leftarrow \text{DMC peak power [W]} \\ \leftarrow \text{DMC onset time [s]} \\ \leftarrow \text{DMC reverberation time [s]} \\ \leftarrow \text{Noise power [W]}. \end{array} \quad (5)$$

In (5), $\theta_{\text{dan}}$ contains the DMC parameters ($[\alpha_1 \in \mathbb{R}^{1 \times S}, \tau_d \in \mathbb{R}^{1 \times S}, \tau_r \in \mathbb{R}^{1 \times S}]$) and the noise parameters ($\alpha_0 \in \mathbb{R}^{1 \times S}$) for each sub-band $s \in S$. A discussion of the model for the DMC can be found in [7]. The model is based on the observation that the PDP $\psi(\tau)$ [W] of the DMC and noise, which describes how the power of a signal is distributed over the time-delay domain, has a base time-delay $\tau_d$ related to the distance between the transmitter and receiver, together with an exponential decay over time-delay (see Eq. (6)), corrupted by complex additive white Gaussian noise with power $\alpha_0$:

$$\psi(\tau) = \begin{cases} \alpha_1 e^{(-\frac{\tau-\tau_d}{\tau_r})} + \alpha_0, & \text{if } \tau > \tau_d \\ \alpha_0, & \text{otherwise.} \end{cases} \quad (6)$$

In Eq. (6), $\alpha_1$, $\tau_d$, $\tau_r$, and $\alpha_0$ are the four parameters which fully describe the DMC and noise characteristics of each sub-band and are gathered in the DMC and noise parameter vector $\theta_{\text{dan}}$. In this work, we will assume that the DMC is spatially white at the transmit and receive side of the measurement system, meaning that they have constant angular power densities. It should be noted that recent works will assume the DMC to be spatially correlated with the SMC. For example, [25, 45] report a correlation between the location of SMC and DMC in the angular domain. In [27], the DMC is modeled as local clusters around the SMC. In [20, 46], it is proposed that the Power Angular Profile (PAP) is to be modeled by a uni-modal Von-Mises distribution.

The covariance matrix $R(\theta_{\text{dan}})$ can be constructed by converting the model in Eq. (6) from the time-delay domain to the frequency domain. To do so, $\kappa(\theta_{\text{dan}}) \in \mathbb{C}^{N_F \times 1}$ [W] is first introduced, which denotes a band-limited sampled version of the Fourier transform of
Eq. (6), and can be defined for a bandwidth B and $M_f$ frequency samples as follows:

$$ \kappa(\theta_{\text{dan}}) = \alpha_0 \mathbf{e}_0 + \frac{\alpha_1}{M_f} \left[ \frac{1}{\beta_d}, \frac{e^{-j2\pi\tau_d}}{\beta_d + j2\pi \frac{1}{M_f}} \dots \frac{e^{-j2\pi\tau_d(M_f-1)}}{\beta_d + j2\pi \frac{(M_f-1)}{M_f}} \right], \quad (7) $$

in which $\mathbf{e}_0 \in \mathbb{N}^{M_f \times 1}$ is a unit vector defined as follows:

$$ \mathbf{e}_0 = [1, 0, \dots, 0]. \quad (8) $$

Furthermore, in Eq. (7), $\beta_d$ is the normalized coherence bandwidth (dimensionless), calculated as follows:

$$ \beta_d = \frac{B_d}{B} = \frac{1}{B \cdot \tau_r}, \quad (9) $$

in which $B_d$ is the coherence bandwidth (Hz), normalized by the bandwidth per sub-band B (Hz). Furthermore, the frequency sampling interval $\Delta_f$ (Hz) can be written as a function of B and $M_f$ as follows:

$$ \Delta_f = \frac{B}{M_f - 1}. \quad (10) $$

After calculating $\kappa(\theta_{\text{dan}})$, the covariance matrix $\mathbf{R}(\theta_{\text{dan}})$ of the DMC and noise can be calculated by applying the Toeplitz-operator [7] as follows:

$$ \mathbf{R}(\theta_{\text{dan}}) = \text{toep}(\kappa(\theta_{\text{dan}}), \kappa(\theta_{\text{dan}})^H). \quad (11) $$

The original data model of the RiMAX estimation algorithm follows the narrowband assumption, hence stating that the SMC and DMC are Kronecker-separable in the spatial and frequency domains in order to keep the algorithm computationally viable [7]. In our measurements section, prior to the processing of the measurement data, we will check for the uncorrelated scattering assumption, which needs to hold in order to apply the Kronecker model. For a large MIMO configuration in a given communication system, the dimensions of the covariance matrix $\mathbf{R}(\theta_{\text{dan}})$ becomes too large to allow for a reliable estimation (remember that it has a size of $[M \times M]$). This becomes a processing burden to accurately model the interaction between transmitter and receiver, so that both ends of a communication system need to be decoupled by applying the well-known Kronecker model. We refer to [47] for a discussion of its applicability and limitations.

In the Kronecker model, the covariance matrices at transmitter and receiver are assumed independent and separable, allowing the full covariance matrix of the channel to be expressed as the Kronecker-product of several smaller matrices. Thus, the covariance matrix $\mathbf{R}$ is assumed to have the following structure [7]:

$$ \mathbf{R} = \mathbf{R}_T \otimes \mathbf{R}_R \otimes \mathbf{R}_f + \alpha_0 \mathbf{I}_M, \quad (12) $$

in which $\mathbf{R}_T \in \mathbb{C}^{M_T \times M_T}$ and $\mathbf{R}_R \in \mathbb{C}^{M_R \times M_R}$ are the covariance matrices at the transmitter and receiver, respectively, describing the spatial distribution of the DMC at both ends, whereas $\mathbf{R}_f \in \mathbb{C}^{M_f \times M_f}$ is the covariance matrix in the frequency domain. The term $\alpha_0 \mathbf{I}_M$ denotes the amount of complex additive white Gaussian measurement noise, which contributes to the stochastic part of the data model. The matrix $\mathbf{I}_M$ is an identity matrix of size $[M \times M]$.

## 3.2 Modeling propagation paths

In this extension of the RiMAX algorithm, we will assume that the transmitter and receiver are equipped with a uniform circular array (UCA), given that its angular resolution capability is uniform since the effective aperture does not change with azimuth angle. The extension to other array configurations is straightforward since only the steering matrices in the following subsection have to be adjusted. It should be noted that other configurations for the antenna array can be used in this algorithm as well, without much modification to the hereinafter explained matrices. The broadband MIMO radio channel model can be expressed in matrix notation to map the time-delay and angles of a propagation path to its complex notation in the frequency domain. We will therefore define the matrix $\mathbf{B}_f(\tau_A) \in \mathbb{C}^{M_f \times P}$ (dimensionless), which maps the time-delays $\tau_A$ of each propagation path $p \in P$ to its related complex exponential $e^{-j2\pi \cdot m \Delta_f \cdot \tau_{A,p}}$ notation as follows:

$$ \mathbf{B}_f(\tau_A) = \begin{bmatrix} e^{-j2\pi \Delta_f \left(-\frac{M_f-1}{2}\right)} \tau_{A,1} & \dots & e^{-j2\pi \Delta_f \left(-\frac{M_f-1}{2}\right)} \tau_{A,P} \\ \vdots & \ddots & \vdots \\ e^{-j2\pi \Delta_f \left(\frac{M_f-1}{2}\right)} \tau_{A,1} & \dots & e^{-j2\pi \Delta_f \left(\frac{M_f-1}{2}\right)} \tau_{A,P} \end{bmatrix}, \quad (13) $$

Similarly, we can define the mapping of both the departing and the arriving angles $\varphi_D$ and $\varphi_A$ to the transmitting and receiving array responses $\mathbf{B}_T^s(\varphi_D) \in \mathbb{C}^{M_T \times P}$ (dimensionless) and $\mathbf{B}_R^s(\varphi_A) \in \mathbb{C}^{M_R \times P}$ (dimensionless) in each sub-band $s \in S$ as follows:

$$ \mathbf{B}_T^s(\varphi_D) = \begin{bmatrix} e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{D,1} - \rho_1)} & \dots & e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{D,P} - \rho_1)} \\ \vdots & \ddots & \vdots \\ e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{D,1} - \rho_{M_T})} & \dots & e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{D,P} - \rho_{M_T})} \end{bmatrix}, \quad (14) $$

and

$$ \mathbf{B}_R^s(\varphi_A) = \begin{bmatrix} e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{A,1} - \rho_1)} & \dots & e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{A,P} - \rho_1)} \\ \vdots & \ddots & \vdots \\ e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{A,1} - \rho_{M_R})} & \dots & e^{-j\frac{2\pi}{\lambda_s}r \cos(\varphi_{A,P} - \rho_{M_R})} \end{bmatrix}, \quad (15) $$
with $r$ being the radius of the UCA, $\lambda_s$ (m) the wavelength in each sub-band $s \in S$, and the vector function $\rho(M_{T/R}) \in \mathbb{R}^{M_{T/R} \times 1}$ (rad) mapping the angle between each antenna in the UCA and a chosen reference axis, defined as follows:

$$ \rho(M_{T/R}) = (0:1:M_{T/R}-1) \frac{2\pi}{M_{T/R}}. \qquad (16) $$

Both matrices $B_T(\varphi_D)$ and $B_R(\varphi_A)$ describe the complex frequency-dependent far-field beam pattern at each antenna array port at the transmitting and receiving array, respectively.

Now let us consider a measurement snapshot of the multidimensional UWB-MIMO channel $\mathbf{h}$. This snapshot can be defined as the instantaneous frequency domain transfer function of the channel between each MIMO antenna pair and is modeled as a superposition of $P$ discrete paths $\mathbf{x}$ plus a contribution of DMC and noise, written as follows:

$$ \mathbf{h} = \sum_{p=1}^{P} \mathbf{x}_p (\theta_{\text{smc},p}) + \mathbf{d}(\theta_{\text{dan}}). \qquad (17) $$

The model for a single propagation path $p$ in sub-band $s$ is given by:

$$ \mathbf{x}_p^s (\theta_{\text{smc},p}^s) = \text{vec} (\mathbf{B}_T^s (\varphi_{D,p}) \diamond \mathbf{B}_R^s (\varphi_{A,p}) \diamond \mathbf{B}_f (\tau_{A,p})) \gamma_p^s, \qquad (18) $$

where the operator $\diamond$ denotes Khatri-Rao (column-wise Kronecker) product and the variable $\gamma_p^s$ denotes the complex amplitude of path $p$ in sub-band $s$. The superposition of $P$ paths in sub-band $s$ can then be written as follows:

$$ \begin{aligned} \mathbf{x}^s (\theta_{\text{smc}}^s) &= \sum_{p=1}^{P} \mathbf{x}_p^s (\theta_{\text{smc},p}^s) \\ &= (\mathbf{B}_T^s (\varphi_D) \diamond \mathbf{B}_R^s (\varphi_A) \diamond \mathbf{B}_f (\tau_A)) \gamma^s \\ &= \mathbf{B}^s (\varphi_D, \varphi_A, \tau_A) \gamma^s \in \mathbb{C}^{M \times 1}, \end{aligned} \qquad (19) $$

in which $\mathbf{B}^s(\varphi_D, \varphi_A, \tau_A)$ represents the steering vector in space and frequency for the entire array, and for all propagation paths, in sub-band $s$.

# 4 Extension to UWB-RIMAX algorithm

## 4.1 Global overview of the algorithm

In a wideband scenario, the reflection coefficients of certain materials in the environment (e.g., furniture, cabinets) can no longer be considered constant with respect to frequency. In addition, the antenna array responses of transmitter and receiver will no longer be frequency-independent and can vary significantly over the total UWB bandwidth. To overcome this issue, the total UWB band was split up into $S$ sub-bands to assure frequency stationarity in each sub-band. This is in agreement with the multiband radio UWB principle [48], where the total

UWB band is split into multiple sub-bands that are separately processed by the receiver, to avoid problems with the restrictions on the analog RF circuit designs. Although the geometrical propagation characteristics (AoD, AoA, and ToA) are frequency-independent, the complex amplitude of each path will vary throughout the UWB frequency band. Hence, we can adopt the narrowband assumption in each sub-band, making it sufficient in terms of measurement accuracy to describe the directional characteristics of the antenna arrays at the center frequency of each sub-band.

In order to properly estimate propagation parameters over the UWB frequency band, an extension to the RiMAX algorithm was developed that can process UWB measurement data and will be referred to as UWB-RiMAX from now on. The global outline for this algorithm is broadly described in Fig. 1. In the next following subsections, we will go deeper into certain aspects of the algorithm.

## 4.2 Initialization (Fig. 1): modeling frequency-dependency

In this subsection, we will discuss how the frequency-dependency of the SMC was modeled in the initialization procedure. As outlined in Fig. 1, the geometrical parameters $\varphi_D, \varphi_A$, and $\tau_A$ are kept constant over frequency, ensuring that these could not appear or disappear from one sub-band onto the next. In a latter part of the algorithm (estimation and optimization of the SMC (see Fig. 1)), these geometrical parameters will be optimized in each sub-band separately, in order not to over-constrain the optimization procedure of the complex amplitudes $\gamma$. This also facilitates the optimization, as there is no strong dependency to be modeled between the sub-bands. Doing so, the geometrical parameters can thus vary slightly over frequency after the optimization procedure is carried out, but jointly estimating them in the initialization procedure ensures that they start their optimization from the same set of values.

### 4.2.1 Maximum likelihood expression

The parameter estimation for propagation paths is initialized by a successive interference cancelation type of grid search, where the detection of paths is based on a single path ML criterion [42]. Taking the logarithm of Eq. (3) and maximizing it with respect to the deterministic and stochastic part gives the following ML criterion:

$$ \begin{bmatrix} \hat{\theta}_{\text{smc}} \\ \hat{\theta}_{\text{dan}} \end{bmatrix} = \underset{\theta_{\text{smc}}, \theta_{\text{dan}}}{\arg \max} \left( -\ln(\det(\mathbf{R}(\theta_{\text{dan}}))) - (\mathbf{h} - \mathbf{x}(\theta_{\text{smc}}))^H \mathbf{R}^{-1}(\theta_{\text{dan}})(\mathbf{h} - \mathbf{x}(\theta_{\text{smc}})) \right). \qquad (20) $$

The parameter estimates $\hat{\theta}_{\text{smc}}$ and $\hat{\theta}_{\text{dan}}$ of the deterministic and stochastic arrays are defined as the maximizing
<figure><img src="image_4.png" /></figure>

arguments of the nonlinear objective function in Eq. (20). The estimation algorithm determines these parameters such that they maximize the likelihood of observing the measured frequency response $\mathbf{h}$ of the radio channel. Since the number of parameters that are nonlinear in this equation is quite large, the solution to this equation is far from easy since it is not practical to perform a multidimensional search for the joint maximization of this log-likelihood function. By splitting the problem into several smaller sub-problems, and choosing several parameter subsets, it becomes practically feasible to solve the joint maximization problem. We then have to maximize the objective function by alternating between the optimization procedures with respect to these subsets. It is evident that we choose $\theta_{smc}$ and $\theta_{dan}$ for both these parameter subsets. This approach exploits the fact that the parameters of the two components of the channel model are asymptotically independent, so that one can decouple the full estimation problem into two separate estimation problems. If we assume the covariance matrix $\mathbf{R}(\theta_{dan})$ of the stochastic process is known (thus the parameters $\theta_{dan}$ are known), the maximization problem of Eq. (20) reduces to:

$$ \hat{\theta}_{smc} = \arg\min_{\theta_{smc}} ((\mathbf{h} - \mathbf{x}(\theta_{dan}))^H \mathbf{R}^{-1}(\theta_{dan}) (\mathbf{h} - \mathbf{x}(\theta_{smc}))) . \quad (21) $$

The term $(\mathbf{h} - \mathbf{x})^H \mathbf{R}^{-1}(\mathbf{h} - \mathbf{x})$ in this equation is the so-called Mahalanobis norm [49]. The ML function in Eq. (21) can be regarded as a nonlinear weighted least squares problem since it is nonlinear in the structural parameters $\theta_{smc}$. More specifically, it is an optimally weighted least squares problem since there is no arbitrary weighting matrix used, but instead, we apply the inverse noise covariance matrix $\mathbf{R}^{-1}(\theta_{dan})$ as weighting matrix. Since the Mahalanobis norm is a non-convex function of the structural parameters $\mathbf{x}$, multiple solutions to this problem exist, and no closed-form solution is available. Therefore, an iterative procedure has to be followed in order to find the ML parameter sets $\hat{\theta}_{smc}$ and $\hat{\theta}_{dan}$ (denoted with a hat-operator) of the “true” values of $\theta_{smc}$ and $\theta_{dan}$, such that they maximize the likelihood of observing the measured frequency response $\mathbf{h}$ of the radio channel given these parameters, as will be discussed in Section 4.3.2.
If the stochastic part of the MIMO channel observation is a zero-mean circular Gaussian i.i.d. process with a covariance matrix equal to $\sigma^2 I$, then the minimization problem is reduced to a classical nonlinear least squares problem [50]. It then becomes a matter of searching for the value $\hat{\theta}_{smc}$ that minimizes the error $\mathbf{h} - \mathbf{x}(\theta_{smc})$. This can also be seen as minimizing the Euclidean norm (the so-called Frobenius norm):

$$ \hat{\theta}_{smc} = \arg \min_{\theta_{smc}} \|\mathbf{h} - \mathbf{x}(\theta_{smc})\|_{F}^{2}. \quad (22) $$

### 4.2.2 Global estimation of SMC parameters

Using the general structure of the data model as given by Eq. (18), which describes the contribution of the specular propagation paths to the channel in a certain sub-band $s \in S$, we can substitute this term in Eq. (19) as follows:

$$ \begin{bmatrix} \hat{\theta}_{\text{smc}}^s \\ \hat{\theta}_{\text{dan}}^s \end{bmatrix} = \arg \min_{\theta_{\text{smc}}, \theta_{\text{dan}}^s} \left( (\mathbf{h}^s - \mathbf{B}^s(\varphi_D, \varphi_A, \tau_A) \mathbf{\gamma}^s)^H \right. \quad (23) \\ \left. \mathbf{R}^{-1}(\theta_{\text{dan}}^s) (\mathbf{h}^s - \mathbf{B}^s(\varphi_D, \varphi_A, \tau_A) \mathbf{\gamma}^s) \right), $$

in which $\mathbf{h}^s$ denotes the frequency response of the MIMO radio channel $\mathbf{h}$ in sub-band $s$, which corresponds to the following range of frequency samples in $\mathbf{h}$: $[(s-1)\frac{M_f-1}{S} + 1 : s\frac{M_f-1}{S} + 1]$. Since the complex amplitudes $\mathbf{\gamma}^s$ are linear in Eq. (23), this minimization problem can be solved directly for $\hat{\mathbf{\gamma}}^s$ given a parameter set $\hat{\mathbf{B}}^s = \mathbf{B}^s(\varphi_D, \varphi_A, \tau_A)$. For any $\hat{\mathbf{B}}^s$, and $\mathbf{R}^s = \mathbf{R}(\theta_{\text{dan}}^s)$, the best linear unbiased estimate (BLUE) is given as follows:

$$ \hat{\mathbf{\gamma}}^s = ((\mathbf{B}^s)^H (\mathbf{R}^s)^{-1} \mathbf{B}^s)^{-1} (\mathbf{B}^s)^H (\mathbf{R}^s)^{-1} \mathbf{h}. \quad (24) $$

Inserting Eq. (24) in Eq. (23) yields the following expression for the ML criterion in sub-band $s$ [7]:

$$ C(\varphi_D, \varphi_A, \tau_A, \mathbf{h}^s) = \left| (\mathbf{B}_p^s)^H \mathbf{h}^s \right|^2, \quad (25) $$

in which $\mathbf{B}_p^s$ can be written as follows:

$$ \mathbf{B}_p^s = \text{vec} \left( \mathbf{B}_T^s (\varphi_{D,p}) \diamond \mathbf{B}_R^s (\varphi_{A,p}) \diamond \mathbf{B}_f (\tau_{A,p}) \right). \quad (26) $$

For the UWB-RIMAX approach, this extends to:

$$ C(\varphi_D, \varphi_A, \tau_A, \{\mathbf{h}^s\}_{s=1}^S) = \sum_{s=1}^{S} C(\varphi_D, \varphi_A, \tau_A, \mathbf{h}^s) \quad (27) $$

$$ = \sum_{s=1}^{S} \left| (\mathbf{B}_p^s)^H \mathbf{h}^s \right|^2 . \quad (28) $$

We can see that this ML criterion estimates the parameters of the electromagnetic waves which extract the largest power from the measured channel $\mathbf{h}$. As such, this approach can allow for a path occurring at a low frequency, but not surviving at a high frequency. In our results, this gives rise to that path having a significant power at the low frequency, but a negligible power at a

high frequency. Thus, we can rewrite this minimization procedure in the following three-dimensional simultaneous search:

$$ (\hat{\varphi}_D, \hat{\varphi}_A, \hat{\tau}_A) = \arg \max_{[\varphi_D, \varphi_A, \tau_A]} \left| C(\varphi_D, \varphi_A, \tau_A, \{\mathbf{h}^s\}_{s=1}^S) \right| \quad (29) $$

The contribution of the SMC can be expressed by the superposition of several individual specular propagation paths. As such, it becomes apparent to maximize the correlation function sequentially with respect to the parameters $\theta_{smc,p}$. This implies that we can minimize the objective function in Eq. (21) sequentially with respect to the different parameter subsets. This sequential optimization for all the geometrical propagation parameters can be written as follows:

$$ \hat{\tau}_{A,p} = \arg \max_{\tau_A} \left| C(\varphi_{D,p}, \varphi_{A,p}, \tau_A, \{\mathbf{h}^s\}_{s=1}^S) \right|, \quad (30) $$

$$ \hat{\varphi}_{D,p} = \arg \max_{\varphi_D} \left| C(\varphi_{D,p}, \varphi_{A,p}, \hat{\tau}_{A,p}, \{\mathbf{h}^s\}_{s=1}^S) \right|, \quad (31) $$

$$ \hat{\varphi}_{A,p} = \arg \max_{\varphi_A} \left| C(\hat{\varphi}_{D,p}, \varphi_{A,p}, \hat{\tau}_{A,p}, \{\mathbf{h}^s\}_{s=1}^S) \right|. \quad (32) $$

The corresponding specular power of this path can then be determined by inserting these ML parameter estimates $[\hat{\varphi}_{D,p}, \hat{\varphi}_{A,p}, \hat{\tau}_{A,p}]$ into Eq. (24). The total specular path can then be described as follows:

$$ \hat{\mathbf{x}}_p^s (\hat{\theta}_{\text{smc},p}^s) = \text{vec} (\mathbf{B}_T^s (\hat{\varphi}_{D,p}) \diamond \mathbf{B}_R^s (\hat{\varphi}_{A,p}) \diamond \mathbf{B}_f (\hat{\tau}_{A,p}) \hat{\gamma}_p^s). \quad (33) $$

The initialization procedure as described above, i.e., finding the global ML estimates for the AoD, AoA, and ToA over UWB frequencies of a propagation path $p$, will subsequently be used for the estimation and optimization of the DMC and SMC parameter sets (see Fig. 1).

## 4.3 Optimization of DMC and SMC

### 4.3.1 Estimation and optimization of DMC

After the initialization procedure (see Fig. 1), the UWB-RIMAX algorithm will estimate and optimize the DMC and noise parameter estimates in each sub-band $s \in S$ separately by using a Gauss-Newton algorithm. This algorithm relies on the calculation of the correlation matrix of the residual signal $(\mathbf{h} - \mathbf{x}(\theta_{smc}))$, which is used to maximize the log-likelihood function with respect to the DMC parameters (see Eq. (20)). It uses the Jacobian of this log-likelihood function to update the DMC parameter set in a next iteration, until convergence of the result is achieved. Since no adjustments were made to the estimation and optimization algorithm of the DMC and noise parameters, we refer to [7] for its full mathematical description.
### 4.3.2 Estimation and optimization of SMC

The optimization of the SMC was performed as presented in Fig. 2. As we can see from this figure, the globally estimated geometrical parameters [$\varphi_D, \varphi_A, \tau_A$] over UWB frequency of each specular propagation path $p$, together with its frequency-dependent amplitude $\gamma^s$ per sub-band $s \in S$, are used as initial values for the search of their optimal values [$\hat{\varphi}_D^s, \hat{\varphi}_A^s, \hat{\tau}_A^s, \hat{\gamma}^s$] in sub-band $s$. By using the globally estimated geometrical parameters as initial values for the Levenberg-Marquardt (LM) algorithm [7], we can ensure that the optimal values in each sub-band lie close to these globally estimated values. By also making these geometrical parameters frequency-dependent per sub-band, we allow for more freedom in the estimation of the optimal value $\gamma^s$ for the frequency-dependent amplitude. In the algorithm for the optimization of the SMC, the switch ensures that subsequent iterations will use the previously optimized (but not yet converged) SMC values for a better optimization. When eventually convergence is achieved, the next sub-band that is to be optimized will again use the globally estimated SMC as initial values for their optimization in sub-band $s + 1$.

## 4.4 Model order selection and the reliability of specular paths

### 4.4.1 Proposed method for model order selection

The total number of specular propagation paths $P$ that is to be extracted from the measurement data is an issue that should be treated with care. This number will naturally influence the ratio between the total power of the SMC components and those that are categorized as DMC. Algorithms like the Akaike information criterion [51] or

the minimum description length [52] can be used to tackle this problem. However, in this work, we will use the approach outlined in [7, 29] that is based on the estimated power of the extracted propagation paths. Because the RiMAX algorithm provides an estimate of the Fisher Information Matrix (FIM) as a by-product, the diagonal elements of the inverse of the FIM are estimates of the variance of the channel parameters in Eq. (5). For each estimated propagation path, it is possible to associate an SNR with it: a path $p$ with an estimate $\hat{\gamma}_p$ for its complex amplitude, has an accompanying SNR $\hat{\rho}_p$ (dimensionless), equal to the following:

$$ \hat{\rho}_p = \frac{|\hat{\gamma}_p|^2}{\mathrm{var}\{|\hat{\gamma}_p|\}} \qquad (34) $$

In Eq. (34), $\mathrm{var}(\cdot)$ denotes the variance of $\hat{\gamma}_p$. It can be proved that the $|\hat{\gamma}_p|$ estimator follows a half-normal distribution with variance $\mathrm{var}(|\hat{\gamma}_p|)$ [7]. The SNR $\hat{\rho}_p$ in Eq. (34) then follows a chi-squared distribution with $2^\circ$ of freedom ($\chi_2^2$). The propagation path $p$ is considered to be unreliable and removed from further analysis if its estimated SNR is smaller than the 90th percentile of $\chi_2^2$, equal to 6.63 dB [29].

### 4.4.2 Stop criterion for path detection

In the flowchart of the UWB-RiMAX algorithm (see also Fig. 1), we can see that an attempt is made to extract more propagation paths from the measured (or residual) channel. It should be noted that (UWB-) RiMAX is an iterative algorithm, implying that in each iteration, it tries to estimate a fixed number of new specular paths from the measured (or residual) channel response. The number of

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Flowchart of the estimation and optimization of the SMC. The first iteration of a new sub-band $s$ will use the globally estimated SMC from the initialization procedure as described in Section 4.2.2 as initial values for the search of its optimal values. The switch ensures that subsequent iterations will use the previously optimized (but not converged) values for a better optimization. When convergence is achieved, the next sub-band $s + 1$ will again use the globally estimated SMC as initial values for the optimization</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** APDP of the synthetic radio channel in the time-delay domain, together with the UWB-RiMAX retained SMC, the DMC, and noise</figcaption></figure>

new propagation paths per iteration was arbitrarily chosen as five, as originally proposed in [7]. It should be noted that this number can be chosen differently depending on the propagation environment. For example, indoor scenarios usually result in stronger multipath behavior than outdoor scenarios. If at least one of the five paths in an iteration succeeds the SNR threshold of 6.63 dB in all sub-bands, the algorithm keeps searching for new paths. Otherwise, the channel is considered to be exhausted of reliable specular paths, and the algorithm is stopped. Subsequently, the resulting SMC and DMC parameter sets are stored, after which they can be used for further analysis.

The model order selection based on Eq. (34) is more suited for this paper's topic than a selection based on information criteria such as the Akaike one. The latter approach calculates an optimal value for the size of the signal subspace as a whole without deciding on the reliability of individual specular paths. In contrast, the path SNR method checks each individual path for its reliability. This approach is more in agreement with the philosophy of DMC, stating that they can also comprise specular

paths which cannot be resolved reliably due to the limited apertures of the measurement equipment or the limited capabilities of the multipath estimation algorithm.

# 5 Evaluation

## 5.1 Synthetic radio channel modeling

In order to evaluate the UWB-RiMAX multipath estimation algorithm as described above, we have generated 200 synthetic channels with controlled parameters for the SMC and DMC parameters. In the following subsections, we first describe how these parameters were modeled in our evaluation procedure. Subsequently, we will describe how the estimated SMC and DMC are compared with their synthetically generated counterparts, in order to evaluate the performance of our algorithm.

### 5.1.1 Contribution of SMC

The generation of the SMC was performed by a ray-tracer in a virtual indoor environment of 15 m × 10 m × 3 m, where for each of the 200 channel realizations in total, random positions are chosen for the transmitter and

<figure><img src="image_7.png" /><figcaption>**Fig. 4** CDF of the MCD metrics between the input and output strongest 90, 95, and 99% propagation paths</figcaption></figure>
<table><caption>**Table 1** Percentiles of the MCD metric as a function of the number of strongest input propagation paths</caption><thead><tr><th rowspan="2">Number of strongest input paths (%)</th><th colspan="3">MCD metric [J] - percentiles</th></tr><tr><th>5th</th><th>50th</th><th>95th</th></tr></thead><tbody><tr><td>90</td><td>0.0227</td><td>0.0443</td><td>0.1066</td></tr><tr><td>95</td><td>0.0228</td><td>0.0477</td><td>0.1136</td></tr><tr><td>99</td><td>0.0235</td><td>0.0554</td><td>0.1404</td></tr></tbody></table>

receiver, with the restriction that they are at least 3 m separated from each other. The height of both transmitter and receiver was fixed at 1.5 m. Further details about the working of the ray-tracer are explained in [53]. For each channel realization, the ray-tracer launched several rays from the randomly positioned transmitter into the environment, each of which could undergo up to 7 reflections, until a total number of 40 paths were obtained at the receiver.

Both the transmitter and receiver consisted of a planar UCA with $M_T = M_R = 8$ antennas, with an inter-element spacing of 0.45 times the wavelength at the highest UWB frequency (10.6 GHz), resulting in a diameter $d$ (m) equal to the following:

$$d = \frac{0.45 \times \lambda_{10.6 \text{ GHz}}}{\sin(\pi/M_T/R)} \\ = 1.8 \text{ cm} \qquad (35)$$

The radio channel was then sampled at $M_f = 4501$ uniformly spaced frequency points ranging from 3.1 up to 10.6 GHz, resulting in a maximum measurable time-delay $\tau_{max}$ of 600 ns. An arbitrary sub-bandwidth of 250 MHz was assumed in the evaluation setup of the UWB-RIMAX algorithm, resulting in a total of 30 UWB frequency sub-bands in which the time-delay bin width $\Delta_{\tau}$ is 4 ns.

### 5.1.2 Contribution of DMC and noise

In order to generate a realization of the DMC and noise process $d(\theta_{dan})$ contributing to the radio channel, we first have to generate a circular Gaussian process with zero mean and a covariance matrix $R(\theta_{dan})$ as follows:

$$d(\theta_{dan}) \sim N_c(0, R(\theta_{dan})) \in \mathbb{C}^{M_f \times 1}. \qquad (36)$$

In order to do so, we will generate a multivariate i.i.d. circular Gaussian process $z$ as follows:

$$z \sim N_c(0, I) \in \mathbb{C}^{M_f \times 1}, \qquad (37)$$

and use a transformation matrix $L(\theta_{dan})$ satisfying:

$$R(\theta_{dan}) = L(\theta_{dan}) \cdot L^H(\theta_{dan}), \qquad (38)$$

so that we can finally compute $d(\theta_{dan})$ as follows:

$$d(\theta_{dan}) = L(\theta_{dan}) \cdot z. \qquad (39)$$

In the equations above, the transformation matrix $L(\theta_{dan})$ can be calculated by using the Cholesky decomposition or the singular value decomposition (SVD) of the covariance matrix $R(\theta_{dan})$. Based on Eq. (11), we can construct this matrix by controlling the four parameters in $\theta_{dan}$. The choice for $\alpha_0$ was chosen as the average noise-level of the SMC contribution to the channel, taking into account an SNR of 20 dB after calculating the average signal power of the SMC. Mathematically, $\alpha_0$ can be calculated as follows:

$$\alpha_0|_{\text{dBm}} = P_{\text{SMC}}|_{\text{dBm}} - \text{SNR}|_{\text{dB}}. \qquad (40)$$

The value for $\alpha_1$ was chosen as the peak value of the SMC contributions, minus 3 dB. The reverberation time $\tau_{rev}$ was calculated based on Eq. (9), in which the coherence bandwidth $B_d$ can be calculated from the root-mean-square (RMS) delay spread $\tau_{RMS}$ [s] as follows:

$$B_d = \frac{1}{\tau_{RMS}}, \qquad (41)$$

<figure><img src="image_8.png" /><figcaption>**Fig. 5** CDF of AoD differences between the input and output strongest 90, 95, and 99% propagation paths</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Fig. 6 CDF of AoA differences between the input and output strongest 90, 95, and 99% propagation paths</figcaption></figure>

and $\tau_{\text{RMS}}$ can be calculated from the SMC power as follows:

$$ \tau_{\text{RMS}} = \sqrt{\frac{\sum_{n} P_{\text{SMC}}(n) \tau_{n}^{2}}{\sum_{n} P_{\text{SMC}}(n)} - \left( \frac{\sum_{n} P_{\text{SMC}}(n) \tau_{n}}{\sum_{n} P_{\text{SMC}}(n)} \right)^{2}}, \quad (42) $$

in which $\tau_n$ (s) is the time-delay of the nth delay-bin, which is equal to $n \times \Delta_{\tau}$, and $P_{\text{SMC}}(n)$ is the PDP of the SMC, which can be calculated by performing an Inverse Discrete Fourier Transform (IDFT) operation on the signal $x(\theta_{\text{smc}})$. This DMC model was then superimposed onto the SMC contribution to the channel, which was generated with the help of the ray-tracer as described in Section 5.1.

## 5.2 Pairing of estimated multipath parameters

In order to evaluate how well the synthetically generated SMC are estimated by the algorithm, we will make use of the multipath component distance (MCD) [54, 55]. The MCD (dimensionless) can be seen as a metric to define the “closeness” between two parameter sets in multipath parameter distance space. It was previously shown that it outperforms the Euclidean distance, and is a suitable metric for combining parameters that have different units (as is the case here with both angles and delays). The pairing of each UWB-RIMAX estimated SMC with its synthetically generated counterpart is done by searching for the smallest MCD between each (output) estimated SMC and the (input) synthetically generated SMC. The MCD $\Omega_{\varphi_{i,j}}$ (dimensionless) between two distinct angles $\varphi_i$ and $\varphi_j$, i.e., the angular distance between both, can be calculated as follows:

$$ \Omega_{\varphi_{i,j}} = \frac{1}{2} \left| \begin{pmatrix} \cos(\varphi_i) \\ \sin(\varphi_i) \end{pmatrix} - \begin{pmatrix} \cos(\varphi_j) \\ \sin(\varphi_j) \end{pmatrix} \right| \quad (43) $$

The MCD $\Omega_{\tau_{i,j}}$ (dimensionless) between two distinct angles time-delay instances $\tau_i$ and $\tau_j$ can be written as follows:

$$ \Omega_{\tau_{i,j}} = \zeta \cdot \frac{|\tau_i - \tau_j|}{\Delta \tau_{\text{max}}} \cdot \frac{\tau_{\text{std}}}{\Delta \tau_{\text{max}}} \quad (44) $$

with $\tau_{\text{std}}$ being the standard deviation of all delays $\tau$, and $\Delta \tau_{\text{max}}$ calculated as follows:

$$ \Delta \tau_{\text{max}} = \max_{i,j} \{|\tau_i - \tau_j|\}. \quad (45) $$

The parameter $\zeta$ is a delay scaling factor to give the delay more weight in the MCD metric when necessary. This factor was chosen as 1, as originally proposed in [55]. Higher values can be chosen to give more weighting to the delay for the pairing of multipath parameters. The delay distance of the MCD was scaled with the normalized delay spread $\frac{\tau_{\text{std}}}{\Delta \tau_{\text{max}}}$. The resulting MCD metric $\Omega_{i,j}$ between two multipath parameter sets i and j can then be calculated as follows:

$$ \Omega_{i,j} = \sqrt{\Omega_{\varphi_{D_{i,j}}}^2 + \Omega_{\varphi_{A_{i,j}}}^2 + \Omega_{\tau_{i,j}}^2}, \quad (46) $$

which can be interpreted as the radius of a circle in the normalized multipath parameter distance space.

It should be noted that in our scenario, there is no one-to-one matching of the (output) estimated SMC to the (input) synthetically generated SMC. The fact that we allow multiple output SMCs to be matched to the same

<table><caption>**Table 2** Percentiles of the angular differences (AoD and AoA) as a function of the number of strongest input propagation paths</caption><thead><tr><th rowspan="2">Number of strongest input paths</th><th rowspan="2"></th><th colspan="3">Angular difference [°] - percentiles</th></tr><tr><th>5th</th><th>50th</th><th>95th</th></tr></thead><tbody><tr><td rowspan="2">90%</td><td>AoD</td><td>-30.47</td><td>0.96</td><td>35.64</td></tr><tr><td>AoA</td><td>-22.17</td><td>0.27</td><td>24.48</td></tr><tr><td rowspan="2">95%</td><td>AoD</td><td>-29.08</td><td>0.98</td><td>37.52</td></tr><tr><td>AoA</td><td>-22.75</td><td>0.27</td><td>25.13</td></tr><tr><td rowspan="2">99%</td><td>AoD</td><td>-29.24</td><td>1.01</td><td>44.12</td></tr><tr><td>AoA</td><td>-27.65</td><td>0.27</td><td>26.76</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Fig. 7** CDF of the ToA differences between the input and output strongest 90, 95, and 99% propagation paths</figcaption></figure>

input SMC is to overcome the case when, e.g., the remainder of the signal after subtraction of previously detected paths at a certain time-delay instance, still contains a sufficient amount of power at this time-delay instance. Newly detected paths then have the same time-delay value and directional parameters, so that we allow them to be matched to the same input SMC.

# 6 Results of simulations

Figure 3 shows the averaged power delay profile (APDP) in the time-delay domain for a randomly generated synthetic radio channel as described in the previous chapter, as well as the UWB-RIMAX estimated SMC, the residual DMC, and noise, and its exponential fit from Eq. (6). Looking at this figure, we can state that the SMC components are well estimated from the radio channel and that the DMC and noise exponential decay (linear in a dBm-scale) of the radio channel is clearly visible and perfectly estimated from the simulations.

## 6.1 Estimation of propagation paths

### 6.1.1 MCD metric

Figure 4 shows the cumulative distribution function (CDF) of the resulting MCDs between the 90, 95, and 99% strongest input paths, and their closest output path in multipath parameter space to the inputs. Table 1 shows the 5th, 50th (median), and 95th percentiles of the MCD metric as a function of the percentage of strongest input paths.

Figure 4 and Table 1 show that when taking into account the 90 or 95% strongest input paths, at least 95% of those are closely matched with the output paths when applying a 0.11 threshold for the MCD metric. When taking into account the 99% strongest input paths, at least 95% of those are closely matched with the output paths when applying a 0.14 threshold for the MCD metric. Since the MCD is a bounded dimensionless metric between 0 and

1, we can safely state that 99% of the strongest input paths are very well estimated by our algorithm.

### 6.1.2 Angular differences

Figures 5 and 6 show the CDF of the angular differences (AoD and AoA) between the 90, 95, and 99% strongest input paths, and their closest output path. Table 2 shows the 5th, 50th (median), and 95th percentiles of the angular differences as a function of the percentage of strongest input paths.

Figures 5 and 6 and Table 2 show that when taking into account the 99% strongest input paths, 90% of them have an AoD error somewhere between −29.24° and +44.12°, and an AoA error somewhere between −27.65° and +26.76°. Table 2 also shows that the median error (50th percentile) is smaller or equal to 1°, regardless of whether it is the AoD or AoA. The table also shows that the range between the 5 and 95th percentiles, so where 90% of the estimated values lie between, gets larger when taking into account more input paths. For example, when considering the 90% strongest input paths, the difference between the 5 and 95th percentiles AoA error is 46.65°, where it is 47.88° when considering the 95% strongest input paths, and 54.41° when considering the 99% strongest input paths.

### 6.1.3 ToA differences

Figure 7 shows the CDF of the ToA differences between the 90, 95, and 99% strongest input paths, and their closest

<table><caption>**Table 3** Percentiles of the ToA differences as a function of the number of strongest input propagation paths</caption><thead><tr><th rowspan="2">Number of strongest input paths (%)</th><th colspan="3">ToA difference (ns) - percentiles</th></tr><tr><th>5th</th><th>50th</th><th>95th</th></tr></thead><tbody><tr><td>90</td><td>-0.81</td><td>0.02</td><td>7.51</td></tr><tr><td>95</td><td>-0.81</td><td>0.02</td><td>7.49</td></tr><tr><td>99</td><td>-1.12</td><td>0.01</td><td>6.69</td></tr></tbody></table>
<figure><img src="image_11.png" /><figcaption>**Fig. 8** Box plot of the differences between the input and output reverberation time</figcaption></figure>

output path. Table 3 shows the 5th, 50th (median), and 95th percentiles of the ToA differences as a function of the percentage of strongest input paths.

Figure 7 and Table 3 show that when taking into account the 99% strongest input paths, 90% of them have a ToA error somewhere between –1.12 and +6.69 ns. Table 3 also shows that the median error (50th percentile) is around 0.02 ns, regardless of whether we take the 90, 95, or 99% strongest input paths.

## 6.2 Estimation of reverberation times

All synthetically generated propagation paths in the APDP were superimposed with a DMC model, as described in Section 5.1.2. Figure 8 shows the resulting differences between all generated and estimated reverberation times as a function of UWB sub-band. This data is represented by means of box plots per frequency band, which indicate the first, second (median), and third quartiles of these differences, and the 1.5 interquartile range from both the lower and upper quartiles.

Figure 8 shows that a median error of less than 4 ns can be obtained for the lower frequencies in the UWB band. For the higher frequencies in the UWB band, a median error of less than 2 ns can be obtained. Overall, more than 75% of our simulation results show an absolute difference between true and estimated reverberation times of less than 4.6 ns in the worst-case scenario, which are more than acceptable results given the relatively low SNR of 20 dB in our evaluation.

## 6.3 Estimation of signal powers

Figure 9 depicts the box plots of the difference between the input and output powers for the total channel, as well as the contributions of both the SMC and DMC power.

From Fig. 9, we can see that the total power in the channel is very well estimated by the UWB-RIMAX algorithm, both at the lower and higher frequencies. We can also observe median differences between the true and estimated SMC powers ranging between –0.3 dB at the lower UWB frequencies, increasing to about 1.5 dB at the higher

<figure><img src="image_12.png" /><figcaption>**Fig. 9** Box plot of the differences between the input and output powers for the total channel, as well as the contributions of both the SMC and DMC power</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 10** Measurement environment: schematic representation</figcaption></figure>

UWB frequencies. These are more than acceptable errors in the estimation of the SMC signal power, if we take into account that the error is still less than 1 dB up to 9 GHz. Next to that, we can observe median absolute differences between the true and estimated DMC powers of about –0.2 dB at the lower UWB frequencies, increasing to about –1 dB at the higher UWB frequencies. The fact that the power of both SMC and DMC is estimated slightly worse at the higher UWB frequencies might be resolved by softening the stop criterion at these frequencies, as it is currently influenced by the SNR of each path, which is thus lower at higher frequencies than at lower frequencies.

Overall, the values for the MCD metric show that our proposed algorithm is able to correctly estimate the most significant input propagation paths, and by analyzing the differences between the input and output SMC, DMC, and total powers, we can state that our proposed algorithm gives a fairly good agreement between the generated input values and the estimated output values.

# 7 Results of measurements

## 7.1 Measurement scenario

### 7.1.1 Measurement environment

In order to test our UWB-RIMAX algorithm, we have performed indoor measurements in a laboratory of Ghent University in Belgium. The schematic representation of the measurement environment is depicted in Fig. 10, with an indication of what we call the long and small side of the laboratory. The long side was approximately 16 m long and 5 m wide (see Fig. 11), and the small side adjacent to it was approximately 8.5 m long and 5 m wide (see Fig. 12). In this environment, most of the equipment consisted of metallic cabinets, tables, computers, and various other hardware. As can be seen from the pictures below, the environment was a very cluttered one.

In this laboratory environment, the indoor radio channel was measured between 15 spatially distinct receiver positions with respect to one single transmitter (see Fig. 10 for an outline of their positions). The exact positions of the receivers with respect to the transmitter were

<figure><img src="image_14.png" /><figcaption>**Fig. 11** Measurement environment: photograph of long side</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Fig. 12** Measurement environment: photograph of small side</figcaption></figure>

measured by means of a digital laser distance meter, which had an accuracy of 2 mm. In total, 8 of the 15 positions were considered LoS scenarios, where there is a direct (free space) path from the transmitter to the receiver. Next to that, 4 of the 15 positions were considered OLoS scenarios, where the free space path from transmitter to the receiver can only be reached via a reflection and/or a diffraction. Lastly, 3 of the 15 positions were considered non-line-of-sight (NLoS) scenarios, where the path from the transmitter to the receiver needs to undergo a transmission through a certain medium (in this case, a plasterboard wall).

### 7.1.2 Channel sounding procedure

At each of the 15 indoor positions, we have carried out channel sounding measurements with a vector network analyzer (VNA) of type Rohde & Schwarz ZNB8 to probe the radio channel ranging from 3.1 up to 10.6 GHz. A planar horizontal virtual [8 × 8] UCA was created at the transmit and receive side of the measurement system by an automatic positioning system (see Fig. 13). The UCA is preferred over other planar array configurations due to its circular configuration, through which it has a 360° full azimuthal field of view (in contrast to linear arrays). A virtual array ensures that the antennas do not suffer from mutual coupling. The inter-element spacing between two adjacent antennas on the virtual UCA was $0.45 \times \lambda_{10.6 \text{ GHz}} = 1.27$ cm. Both antennas were omni-directional UWB antennas in the azimuthal plane of type Electro-Metrics EM-6865 [56], placed 1.5 m above the ground.

In the 3.1 to 10.6 GHz UWB frequency band, the VNA sampled the complex gain between each pair of Tx and Rx antennas over $M_f = 7501$ uniformly spaced frequency point, with a resolution bandwidth of 10 kHz. This complex gain corresponds with the $S_{21}$-scattering parameter, which is the ratio of the output reflected power

wave divided by the input incident power wave, where all values are expressed as complex quantities. Our measurements resulted in a maximum measurable time-delay $\tau_{max}$ of 1000 ns. The cables connecting Tx and Rx antennas were included in the calibration of the VNA to exclude their influence from the measurement data. Measurements were conducted outside of working hours to get a static radio channel without any movement.

<figure><img src="image_16.png" /><figcaption>**Fig. 13** Automated positioning system and UWB antenna, to form a planar horizontal virtual UCA</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Fig. 14** Analysis of the stationarity bandwidth of the measurement data for different values of the MSR</figcaption></figure>

### 7.1.3 Frequency stationarity

Prior to the processing of the measurement data in the UWB band, we will first check the uncorrelated scattering (US) assumption which is often assumed for the modeling of wireless channels. The US assumption states that the second-order statistical description of the radio channel is independent of the absolute frequency, which implies that channel contributions at different time-delays are uncorrelated. In order to evaluate the US assumption in the UWB band, we will apply the procedure explained in [57], which proposes a test that defines a frequency stationarity region (or stationarity bandwidth) in which the US assumption holds. This test is based on the definition of a minimum stationarity region, which is a power spectral density in the frequency domain in which the US property locally holds. Subsequently, the test measures how many neighboring MSRs can be formed, of which their adjacent overlap in the frequency domain exceeds a certain threshold. The total number of neighboring MSRs which exceed this threshold thus form a frequency stationarity region, corresponding with

a certain number of frequency samples (and thus a certain bandwidth).

We have analyzed the stationarity bandwidths for the different measurement positions for MSR values of 50 MHz (corresponding with 50 frequency samples) and 100 MHz (corresponding with 100 frequency samples) with a threshold value of 0.9 to decide if neighboring MSRs can be considered stationary in the frequency domain. Choosing a smaller value for the MSR would deteriorate its resolution. The results of this analysis can be found in Fig. 14.

Figure 14 shows that it is reasonable to assume a sub-bandwidth of 250 MHz in the setup of the UWB-RIMAX algorithm, in which we can thus state that the US assumption holds. This analysis proves that it is even reasonable to assume sub-bandwidths of 300 MHz.

## 7.2 Results of SMC and DMC

Figure 15 shows the 10 strongest estimated geometrical propagation paths in the environment. Their length is an indication of their relative power.

<figure><img src="image_18.png" /><figcaption>**Fig. 15** Estimated geometrical propagation paths in the environment</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. 16** Mean values of the measured powers in the radio channel as a function of UWB frequencies: LoS scenario</figcaption></figure>

From Fig. 15, we can see that our algorithm is able to estimate the correct geometrical propagation paths in the environment. For receiver positions 11 and 12, it looks like the algorithm estimated the wrong angles, but it should be noted that there was a metallic cabinet in the long side of the environment on which the paths apparently scattered from transmitter to receiver.

Figures 16, 17, and 18 show the mean values of the measured SMC, DMC, and SMC+DMC (total reconstructed) powers in the radio channel as a function of UWB frequencies. The power of the SMC part of the radio channel in each sub-band $P_{SMC}^s$ can be calculated from Eq. (19), whilst the power of the DMC part of the radio channel in each sub-band $P_{DMC}^s$ can be calculated from Eq. (39). The total reconstructed power in each sub-band $P_{Tot}^s$ can then be calculated by summing over both $P_{SMC}^s$ and $P_{DMC}^s$. Figure 16 presents the results of the LoS scenarios, Fig. 17 presents the results of the OLoS scenarios, and Fig. 18 presents the results of the NLoS scenarios.

Figures 16, 17, and 18 show that the SMC+DMC (total reconstructed) powers match very well with the measured power in the channel across the UWB frequency band. For the LoS scenario, the difference between the

measured and the reconstructed power resulted in a minimum (underestimated) value of –1.61 dB and a maximum (overestimated) value of 1.20 dB. For the OLoS scenario, the minimum difference was –1.38 dB, and the maximum difference was 0.55 dB. For the NLoS scenario, the minimum difference was –1.73 dB, and the maximum difference was 0.52 dB. Figure 16 shows that the SMC power dominates the DMC power in a LoS scenario, especially for the higher UWB frequencies. This can be explained by the fact for higher frequencies, the wavelength of the transmitted ray is small, meaning that an incident ray on a surface will encounter little effect from the roughness of the material and will reflect on it specularly. In an OLoS scenario, Fig. 17 shows that both the SMC and DMC power result in comparable power levels. This can be explained by the fact that the inherent necessity of a reflection from transmitter to receiver will automatically generate more DMC in the channel, originating from the roughness of the surface on which the reflection of the path occurs. The more reflections a propagation path undergoes, the higher the chance that the incident wave at the receiver will have encountered diffuse scattering along the way. In an NLoS scenario, Fig. 18 shows that

<figure><img src="image_20.png" /><figcaption>**Fig. 17** Mean values of the measured powers in the radio channel as a function of UWB frequencies: OLoS scenario</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>Fig. 18 Mean values of the measured powers in the radio channel as a function of UWB frequencies: NLoS scenario</figcaption></figure>

the DMC power dominates the SMC power, especially for the lower UWB frequencies. This can be explained by the fact for lower frequencies, the wavelength of the transmitted ray is large, meaning that an incident ray on a surface with irregularities comparable in size to its wavelength will cause this ray to be scattered at many angles rather than just at one angle (as is the case with a specular reflection). This diffuse scattering typically occurs at lower frequencies, which explains why the contribution of the DMC is significantly higher than those of the SMC. The DMC power ratio will be discussed in the next section.

## 7.3 DMC power ratio

The DMC power ratio $p_r^s$ can be quantified as the relative power attributable to the DMC part of the measured radio channel and can be written as follows in each sub-band s:

$$p_r^s = \frac{\sum_{n=1}^{M_f} P_{\text{DMC}}^s(n)}{\sum_{n=1}^{M_f} P_{\text{Meas}}^s(n)} \quad (47)$$

Figure 19 shows that the relative power attributable to the DMC part of the measured radio channel is higher for NLoS scenarios than for OLoS and LoS scenarios. In a LoS scenario, the presence of a strong direct path (LoS component) will dominate in power over the reflected paths, such that the relative DMC power ratio is lower than for OLoS or NLoS scenarios. The DMC power represents up to 50% of the total measured power for the lower UWB frequencies down to 30% for the higher UWB frequencies. As in the previous section, we know that this is due to the fact that an incident ray on an electrically small surface will encounter more effect from the roughness of its material, causing it to reflect on it diffusely. In contrast, in an OLoS scenario, the path between transmitter and receiver has to undergo one or more reflections, giving rise to more diffuse scattering along the way. This effect can be especially seen around 3 to 4 GHz. Finally, in an NLoS scenario, the DMC power represents up to 60% of the total measured power for the lower UWB frequencies (3.1 to 7 GHz), whilst it still represents up to 50% of the total measured power for the higher UWB frequencies (7 to 10.6 GHz). The necessity of a transmission from

<figure><img src="image_22.png" /><figcaption>Fig. 19 Relative DMC power ratio $p_r$ as a function of UWB frequencies for LoS, OLoS, and NLoS scenarios</figcaption></figure>
Tx to Rx in these NLoS scenarios between different media gives rise to more diffuse behavior of the waves, caused by the inherent roughness of the surfaces which the wave has to pass through. Overall, we can see that the DMC power ratio is lower for these higher frequencies, due to the fact that the encountered surfaces from Tx to Rx are electrically larger, resulting in more specular reflections.

## 8 Conclusions

This work presented an extension of the RiMAX multipath estimation algorithm, facilitating the analysis of frequency-dependent propagation parameters for ultra-wideband (UWB) channel modeling. The newly proposed algorithm is capable of tracking both the specular and dense multipath components (SMC and DMC) over different UWB sub-bands by estimating the most likely geometrical propagation parameters such as angle of departure (AoD), angle of arrival (AoA), and time-delay of arrival (ToA) of the SMC occurring throughout these bands. In order to do so, a maximum likelihood (ML) criterion was proposed, allowing to estimate those propagation paths representing a significant amount of power in the measured radio channel over all sub-bands. The model for the DMC and noise was based on the observation that its power in the time-delay domain follows an exponential decay, and is related to the distance between transmitter and receiver. This algorithm was tested by generating synthetic radio channels in an indoor environment, with contributions to these channel consisting of SMC generated with a ray-tracer, and an associated DMC model. The geometric parameters were estimated from these synthetic models and matched with the generated parameter values with the help of the multipath component distance (MCD) metric. This metric was also used to define the closeness between the generated and estimated parameters in multipath parameter distance space. Next to that, the powers of the reconstructed SMC, DMC, and total channel were compared with their generated values.

Our simulation results show that the newly designed UWB-RiMAX algorithm can match up to 99% of the SMC parameters according to the MCD metric, and that the DMC reverberation time known from the theory of room electromagnetics can be estimated on average with an error margin of less than 2 ns throughout the UWB frequency band.

Our measurement results indicate the strong presence of DMC in an indoor environment, in which the DMC power represents up to 50% of the total measured power for the lower UWB frequencies, caused by the fact that the encountered surfaces from Tx to Rx are electrically smaller, resulting in more diffuse reflections. This DMC power ratio reduces to around 30% for the higher UWB frequencies since these surfaces will react more as specular reflectors.

Future work consists of performing a more extensive measurement campaign in line-of-sight (LoS), obstructed-LoS and non-LoS environments, and the estimation of SMC and DMC parameters in these different scenarios. The results from this measurement campaign will be used in a localization algorithm in order to estimate the location of a mobile receiver node with the help of a single transmitter node.

### Abbreviations

AoA: Angle of arrival; AoD: Angle of departure; APDP: Averaged power delay profile; BLUE: Best linear unbiased estimate; CDF: Cumulative distribution function; DMC: Dense multipath components; ESPRIT: Estimation of signal parameters via rotational invariance techniques; FIM: Fisher information matrix; GBSCM: Geometry-based stochastic channel model; IDFT: Inverse discrete fourier transform; i.i.d.: Independent and identically distributed; LM: Levenberg-Marquardt; LoS: Line-of-sight; MCD: Multipath component distance; MIMO: Multiple-input multiple-output; ML: Maximum likelihood; MSR: Minimum stationarity region; NLoS: Non-line-of-sight; OLoS: Obstructed-line-of-sight; PAP: Power angular profile; PDP: Power delay profile; RMS: Root-mean-square; Rx: Receiver; SAGE: Space-alternating generalized expectation-maximization; SMC: Specular multipath components; SNR: Signal-to-noise ratio; SVD: Singular value decomposition; ToA: Time-delay of arrival; Tx: Transmitter; UCA: Uniform circular array; US: Uncorrelated scattering; UWB: Ultra-wideband; VNA: Vector network analyzer; WPAN: Wireless personal area network

### Acknowledgements

The author would also like to express his gratitude to Dr. Rein Houthooft, M. Sc., for his many insights, fruitful discussions, and shedding brightening lights on this subject.

### Funding

Brecht Hanssens is funded by the Agency for Innovation and Entrepreneurship. Emmeric Tanghe is a Post-Doctoral fellow of the Research Foundation Flanders (FWO). This work was carried out in the framework of COST Action CA15104 IRACON, and the Excellence of Science (EOS) project MULTI-Service Wireless NETwork (MUSEWINET). This research was supported by the FWO project G027714N.

### Authors' contributions

BH was responsible for carrying out the measurement campaign, most of the data processing and writing this manuscript. ET, DPG, and ML were responsible for assisting and/or writing parts of the UWB-RiMAX algorithm. DP and LM were responsible for proofreading this work and providing meaningful feedback. CO and WJ were responsible to a greater extent for proofreading this work and providing meaningful feedback on the methodology of the algorithm. All authors read and approved the final manuscript.

### Availability of data and materials

The data can be made available for scientific purposes only by sending a request via e-mail to Brecht.Hanssens@UGent.be.

### Competing interests

The authors declare that they have no competing interests.

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

### Author details

¹INTEC - WAVES, Ghent University - imec, Ghent, Belgium. ²IEMN/TELICE, University of Lille 1, Lille, France. ³ICTEAM, Université Catholique de Louvain, Louvain-la-Neuve, Belgium.

Received: 13 March 2017 Accepted: 8 June 2018

Published online: 27 June 2018

### References

1. J Salmi, J Poutanen, K Haneda, A Richter, V-M Kolmonen, P Vainikainen, AF Molisch, in *European Conference on Antennas and Propagation (EuCAP)*.
Incorporating diffuse scattering in geometry-based stochastic MIMO channel models (IEEE, Barcelona, 2010), pp. 1–5

2. R Verdone, A Zanella, *Pervasive Mobile and Ambient Wireless Communications: COST Action 2100. Signals and Communication Technology*. (Springer, London, 2012). https://doi.org/10.1007/978-1-4471-2315-6

3. B Hanssens, M-T Martínez-Ingólés, A Bamba, E Tanghe, J-M Molina-García-Pardo, DP Gaillot, M Liénard, C Oestges, L Martens, W Joseph, in *European Conference on Antennas and Propagation (EuCAP)*. Estimation of room electromagnetics reverberation time for ultra-wideband indoor channels (IEEE, The Hague, 2014), pp. 527–528

4. A Bamba, M-T Martínez-Ingólés, DP Gaillot, E Tanghe, B Hanssens, J-M Molina-García-Pardo, M Lienard, L Martens, W Joseph, Experimental investigation of electromagnetic reverberation characteristics as a function of UWB frequencies. IEEE Antennas Wirel. Propag. Lett. 14, 859–862 (2015). https://doi.org/10.1109/LAWP.2014.2382873

5. J Poutanen, J Salmi, K Haneda, V-M Kolmonen, P Vainikainen, Angular and shadowing characteristics of dense multipath components in indoor radio channels. IEEE Trans. Antennas Propag. 59(1), 245–253 (2011). https://doi.org/10.1109/TAP.2010.2090474

6. F Mani, F Quitin, C Oestges, Directional spreads of dense multipath components in indoor environments: experimental validation of a ray-tracing approach. IEEE Trans. Antennas Propag. 60(7), 3389–3396 (2012). https://doi.org/10.1109/TAP.2012.2196942

7. A Richter, Estimation of radio channel parameters: models and algorithms, PhD thesis. Technische Universität Ilmenau, Fakultät für Elektrotechnik und Informationstechnik. Ilmenau, Germany (2005)

8. A Batra, Multi-band OFDM physical layer proposal for IEEE 802.15 task group 3a (2003). https://doi.org/IEEE802.15-03/268r2

9. (FCC) FCC, First report and order 02-48. Technical report (2002)

10. T Kaiser, F Zheng, E Dimitrov, An overview of ultra-wide-band systems with MIMO. Proc. IEEE. 97(2), 285–312 (2009). https://doi.org/10.1109/JPROC.2008.2008784

11. O Boric-lubecke, VM Lubeckel, A Host-madsen, D Samardzija, K Cheung, in *International Conference on Telecommunications in Modern Satellite, Cable and Broadcasting Services*. Doppler Radar Sensing of Multiple Subjects in Single and Multiple Antenna Systems, (2005), pp. 7–11. https://doi.org/10.1109/TELSKS.2005.1572052

12. S Gezici, Z Tian, GB Giannakis, H Kobayashi, AF Molisch, HV Poor, Z Sahinoglu, Localization via ultra-wideband radios: a look at positioning aspects for future sensor networks. IEEE Signal Proc. Mag. 22(4), 70–84 (2005). https://doi.org/10.1109/2010RS004561

13. F Viani, P Rocca, G Oliveri, D Trinchero, A Massa, Localization, tracking, and imaging of targets in wireless sensor networks: an invited review. Radio Sci. 46(5), 1–4 (2011). https://doi.org/10.1029/2010RS004561

14. L Zwirello, T Schipper, M Harter, T Zwick, UWB localization system for indoor applications: concept, realization and analysis. J. Electr. Comput. Eng. 1, 1–11 (2012). https://doi.org/10.1155/2012/849638

15. R Zetik, H Yan, E Malz, S Jovanoska, G Shen, RS Thomä, R Salman, T Schultze, R Tobera, H-I Willms, L Reichardt, M Janson, T Zwick, W Wiesbeck, T Deiler, J Thielecke, in *Ultra-Wideband Radio Technologies for Communications, Localization and Sensor Applications*. Chap. 9. Cooperative localization and object recognition in autonomous UWB sensor networks (IntechOpen, London, 2013). https://doi.org/10.5772/55077

16. B Hanssens, D Plets, E Tanghe, C Oestges, DP Gaillot, M Liénard, L Martens, W Joseph, in *IEEE Antennas and Propagation Society International Symposium (APS/URSI)*. An indoor localization technique based on ultra-wideband AoD/AoA/ToA estimation, (2016), pp. 1445–1446. https://doi.org/10.1109/APS.2016.7696429

17. V Degli-Esposti, D Guiducci, A De'Marsi, P Azzi, F Fuschini, An advanced field prediction model including diffuse scattering. IEEE Trans. Antennas Propag. 52(7), 1717–1728 (2004). https://doi.org/10.1109/TAP.2004.831299

18. A Richter, RS Thomä, in *COST 273 Technical Document TD(03) 198*. Parametric modelling and estimation of distributed diffuse scattering components of radio channels, (2003), pp. 24–26

19. A Richter, C Schneider, in *International Symposium on Wireless Personal Multimedia Communications (WPMC)*. Parameter estimation results of specular and dense multipath components in micro-cell scenarios (IEEE, Abano Terme (Padova), 2004), pp. 90–94

20. M Käske, RS Thomä, in *European Conference on Antennas and Propagation (EuCAP)*. Maximum-likelihood based estimation of angular parameters of dense-multipath-components (IEEE, Lisbon, 2015), pp. 1–6

21. AF Molisch, D Cassioli, CC Chong, S Emami, A Fort, B Kannan, J Karedal, J Kunisch, HG Schantz, K Siwiak, MZ Win, A comprehensive standardized model for ultrawideband propagation channels. IEEE Trans. Antennas Propag. 54(11), 3151–3166 (2006). https://doi.org/10.1109/TAP.2006.883983

22. Thomä, M Landmann, G Sommerkorn, A richter, in *IEEE Instrumentation and Measurement Technology Conference (IMTC)*. multidimensional high-resolution channel sounding in mobile radio, (2004), pp. 257–262. https://doi.org/10.1109/IMTC.2004.1351040

23. M Landmann, K Sivasondhivat, J-i Takada, I Ida, RS Thomä, Polarization behavior of discrete multipath and diffuse scattering in urban environments at 4.5 GHz. EURASIP J. Wirel. Commun. Netw. 2007(1), 60 (2007). https://doi.org/10.1155/2007/57980

24. M Käske, M Landmann, RS Thomä, in *European Conference on Antennas and Propagation (EuCAP)*. Modelling and synthesis of dense multipath propagation components in the angular domain (IEEE, Berlin, 2009), pp. 2641–2645

25. F Quitin, C Oestges, F Horlin, P De Doncker, in *European Conference on Antennas and Propagation (EuCAP)*. Diffuse multipath component characterization for indoor MIMO channels (IEEE, Barcelona, 2010), pp. 1–5

26. DP Gaillot, E Tanghe, W Joseph, P Laly, VC Tran, M Liénard, L Martens, Polarization properties of specular and dense multipath components in a large industrial hall. IEEE Trans. Antennas Propag. 63(7), 3219–3228 (2015). https://doi.org/10.1109/TAP.2015.2430374

27. F Quitin, C Oestges, F Bellens, S Van Roy, F Horlin, P De Doncker, in *IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC)*. Extracting specular-diffuse clusters from MIMO channel measurements, (2011), pp. 940–944. https://doi.org/10.1109/PIMRC.2011.6140107

28. A Richter, J Salmi, V Koivunen, in *European Conference on Antennas and Propagation (EuCAP)*. Distributed scattering in radio channels and its contribution to MIMO channel capacity (IEEE, Nice, 2006), pp. 1–7

29. E Tanghe, DP Gaillot, M Liénard, L Martens, W Joseph, Experimental analysis of dense multipath components in an industrial environment. IEEE Trans. Antennas Propag. 62(7), 3797–3805 (2014). https://doi.org/10.1109/TAP.2014.2321162

30. J Kunisch, J Pamp, in *IEEE Conference on Ultra Wideband Systems and Technologies (UWBST)*. Measurement results and modeling aspects for the UWB radio channel, (2002), pp. 19–24. https://doi.org/10.1109/UWBST.2002.1006310

31. A Richter, RS Thomä, in *IEEE Vehicular Technology Conference (VTC)*. Joint maximum likelihood estimation of specular paths and distributed diffuse scattering, (2005), pp. 11–15. https://doi.org/10.1109/VETECS.2005.1543239

32. UT Virk, K Haneda, J-F Wagen, in *European Conference on Antennas and Propagation (EuCAP)*. Dense multipath components add-on for COST 2100 channel model (IEEE, Lisbon, 2015), pp. 1–5

33. E Tanghe, DP Gaillot, W Joseph, M Liénard, W De Ketelaere, L Martens, Capacity simulation and analysis of an IEEE 802.11n System in a Residential House. Eur. Conf. Antennas Propag., 1671–1675 (2016). https://doi.org/10.1109/EuCAP.2016.7481559

34. K Saito, J-i Takada, M Kim, Dense multipath component characteristics in 11 GHz-band indoor environments. IEEE Trans. Antennas Propag. 65(9), 4780–4789 (2017). https://doi.org/10.1109/TAP.2017.2728087

35. JB Andersen, J Nielsen, GF Pedersen, G Bauch, M Herdin, Room electromagnetics. IEEE Antennas Propag. Mag. 49(2), 27–33 (2007). https://doi.org/10.1109/MAP.2007.376642

36. G Steinböck, T Pedersen, BH Fleury, W Wang, R Raulefs, Distance dependent model for the delay power spectrum of in-room radio channels. IEEE Trans. Antennas Propag. 61(8), 4327–4340 (2013). https://doi.org/10.1109/TAP.2013.2260513. arXiv:1011.1669v3

37. M Haardt, Efficient one-, two-, and multidimensional high-resolution array signal processing. PhD thesis, Technische Universität München, Shaker Verlag GmbH, Aachen, Germany (1996)

38. BH Fleury, M Tschudin, R Heddergott, D Dahlhaus, KI Pedersen, Channel parameter estimation in mobile radio environments using the SAGE algorithm. IEEE J Sel. Areas Commun. 17(3), 434–450 (1999). https://doi.org/10.1109/49.753729
39. DP Gaillot, E Tanghe, P Stefanut, W Joseph, M Liénard, P Degauque, L Martens, in *European Conference on Antennas and Propagation (EuCAP)*. Robustness of high-resolution channel parameter estimators in the presence of dense multipath components (IEEE, Rome, 2011), pp. 3619–3622

40. E Tanghe, DP Gaillot, W Joseph, M Lienard, P Degauque, L Martens, Robustness of high-resolution channel parameter estimators in presence of dense multipath components. Electron. Lett. **48**(2), 130 (2012). https://doi.org/10.1049/el.2011.3454

41. M Bengtsson, B Volcker, in *IEEE Vehicular Technology Conference (VTC)*. On the estimation of azimuth distributions and azimuth spectra, (2001), pp. 1612–1615. https://doi.org/10.1109/VTC.2001.956471

42. K Haneda, J-i Takada, in *Ultra-Wideband Systems and Technologies (UWBST)*. An application of SAGE algorithm for UWB propagation channel estimation, (2003), pp. 483–487. https://doi.org/10.1109/UWBST.2003.1267889

43. K Haneda, J-i Takada, T Kobayashi, in *Ultra-Wideband Systems and Technologies (UWBST)*. Experimental evaluation of a SAGE algorithm for ultra wideband channel sounding in an anechoic chamber, (2004), pp. 66–70. https://doi.org/10.1109/UWBST.2004.1320937

44. M Debbah, RR Müller, MIMO channel modeling and the principle of maximum entropy. IEEE Trans. Inf. Theory. **51**(5), 1667–1690 (2005). https://doi.org/10.1109/TIT.2005.846388

45. J Poutanen, J Salmi, K Haneda, V-M Kolmonen, F Tufvesson, P Vainikainen, Propagation characteristics of dense multipath components. IEEE Antennas Wirel. Propag. Lett. **9**, 791–794 (2010). https://doi.org/10.1109/LAWP.2010.2064751

46. CB Ribeiro, A Richter, V Koivunen, Joint angular- and delay-domain MIMO propagation parameter estimation using approximate ML method. IEEE Trans. Signal Process. **55**(10), 4775–4790 (2007). https://doi.org/10.1109/TSP.2007.896247

47. C Oestges, in *IEEE Vehicular Technology Conference (VTC)*. Validity of the Kronecker Model for MIMO correlated channels, (2006), pp. 2818–2822. https://doi.org/10.1109/VETECS.2006.1683382

48. WP Siriwongpairat, W Su, M Olfat, KJR Liu, Multiband-OFDM MIMO coding framework for UWB communication systems. IEEE Trans. Signal Process. **54**(1), 214–224 (2006). https://doi.org/10.1109/TSP.2005.861092

49. PC Mahalanobis, On the generalised distance in statistics. Proc. Natl Inst. Sci. India. **2**, 49–55 (1936)

50. LL Scharf, C Demeure, *Statistical Signal Processing: Detection, Estimation, and Time Series Analysis. Addison-Wesley Series in Electrical and Computer Engineering*. (Addison-Wesley Pub. Co., Boston, 1991)

51. H Akaike, A new look at the statistical model identification. IEEE Trans. Autom. Control. **19**(6), 716–723 (1974). https://doi.org/10.1109/TAC.1974.1100705. arXiv:1011.1669v3

52. Y Hua, AB Gershman, Q Cheng, *High-Resolution and Robust Signal Processing*. (Marcel Dekker, New York, 2004), p. 448

53. J Pascual-García, J-M Molina-García-Pardo, M-T Martínez-Inglés, J-V Rodríguez, L Juan-Llácer, Fast and accurate electric field estimation from a single ray tracing simulation. Appl. Comput. Electromagn. Soc. **J. 30**(6), 608–618 (2015)

54. M Steinbauer, H Ozcelik, H Hofstetter, CF Mecklenb Auker, E Bonek, How to quantify multipath separation. IEICE Trans. Electron. Spec. Issue Signals Syst. Electron. Technol. **85**(3), 552–557 (2002)

55. N Czink, P Cera, J Salo, E Bonek, J-P Nuutinen, J Ylitalo, Improving clustering performance using multipath component distance. Electron. Lett. **42**(1), 33 (2006). https://doi.org/10.1049/el:20063917

56. Electro-Metrics, Omni-directional EM-6865 wideband antenna 2 GHz–18 GHz. https://electro-metrics.com/product/antenna-omni-directional-wideband-em-6865-2-ghz-18-ghz/

57. L Bernadó, Non-stationarity in vehicular wireless channels. PhD thesis, Institut für Nachrichtentechnik und Hochfrequenztechnik, Technische Universität Wien, Vienna, Austria (2012)

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article