Next Article in Journal
Influences of High-Sulphur Fly Ash on the Properties of Lightweight Cement-Treated Materials Subjected to Sulphate Corrosion
Next Article in Special Issue
UWB Indoor Localization Using Deep Learning LSTM Networks
Previous Article in Journal
A Novel Variable Step Size Incremental Conductance Method with an Adaptive Scaling Factor
Previous Article in Special Issue
A Fine Fast Acquisition Scheme for a Communication and Navigation Fusion System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gaussian Process Modeling of Specular Multipath Components

1
School of Electronics and Telecommunications, Hanoi University of Science and Technology, Hanoi 112400, Vietnam
2
Signal Processing and Speech Communication Lab, Graz University of Technology, 8010 Graz, Austria
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2020, 10(15), 5216; https://doi.org/10.3390/app10155216
Submission received: 30 May 2020 / Revised: 17 July 2020 / Accepted: 24 July 2020 / Published: 29 July 2020
(This article belongs to the Special Issue Indoor Localization Systems: Latest Advances and Prospects)

Abstract

:
The consideration of ultra-wideband (UWB) and mm-wave signals allows for a channel description decomposed into specular multipath components (SMCs) and dense/diffuse multipath. In this paper, the amplitude and phase of SMCs are studied. Gaussian Process regression (GPR) is used as a tool to analyze and predict the SMC amplitudes and phases based on a measured training data set. In this regard, the dependency of the amplitude (and phase) on the angle-of-arrival/angle-of-departure of a multipath component is analyzed, which accounts for the incident angle and incident position of the signal at a reflecting surface—and thus for the reflection characteristics of the building material—and for the antenna gain patterns. The GPR model describes the similarities between different data points. Based on its model parameters and the training data, the amplitudes of SMCs are predicted at receiver positions that have not been measured in the experiment. The method can be used to predict a UWB channel impulse response at an arbitrary position in the environment.

1. Introduction

Future 5G wireless communication technologies and the Internet of Things (IoT) paradigm will be characterized by supporting a variety of services with high quality requirements, addressing performance metrics such as reliability, latency, data throughput, and resource-efficient use of the infrastructure [1,2,3]. Spatial location information is expected to become an indispensable feature of these emerging wireless networks, considering that much more accurate position information will become available with those new generations of wireless networks [4,5,6]. These references elaborate, how such “location-awareness” can be exploited for improving the robustness and performance of wireless networks. Location-awareness is founded in the observation that many parameters of the propagation channel are directly linked with the geometry of the scenario and environment, starting from the distance-dependent path loss, spatial correlation of shadowing, to the arrival angle of individual multipath components. The availability of such geometric models thus allows for efficient and robust location-based routing and medium access schemes, or geometry-based beam steering, to name just a few popular examples.
The prerequisites for implementing location awareness are accurate position estimation and the mapping of the propagation conditions, both of which must go hand in hand. The well-known simultaneous localization and mapping (SLAM) algorithms can be used to build a geometric model from past received signals, while locating the mobile user [7,8,9]. Note that these references are based on the use of wireless signals, while SLAM was originally developed by the robotics community for efficiently fusing the data from optical and odometry sensors. By exploiting the position-related information acquired from the multipath channels, localization and mapping information is obtained. The obtained localization and mapping errors were continuously reduced from meters to centimeters, while moving from conventional wideband transceivers, through ultra-wideband (UWB), to mm-wave systems [7,8,9,10,11]. The robustness of the acquired location information strongly depends on the surrounding environment and hence on multipath propagation [12,13,14].
Strongly related to those works are papers on radio positioning systems that exploit map information. For instance, low-level descriptors of the channel, such as the RSS, the time-of-flight (ToF), the angle-of-arrival (AoA) are mapped to a location using fixed geometry models and triangulation/trilateration methods [15,16]. In Reference [17], delay information of specular multipath components (SMCs) is used to enhance the robustness and reduce the infrastructure demand of a UWB localization system, see also Reference [12]. The concept extends straigthforwardly towards mm-wave systems, which will also make use of angle information [18].
These works indicate a strong link between map and position information. Specifically, the arrival time, arrival angles, and departure angles of so-called “specular multipath components” are all directly related to simple geometric models. It is thus common to separate these SMCs and use geometry-based models to describe their parameters. A popular approach is the use of clusters of scatter points for the modeling, characterized by a delay spread and angle spread [19,20,21], allowing to quantify the second-order statistic of geometry-related multipath components, that is, SMCs. Vice-versa, SMCs allow for a prediction of these channel parameters, given that the geometric models are available.
In contrast to this delay and angle modeling by means of SMCs, the modeling of the amplitudes of the SMCs is much harder to achieve. This is because many factors need to be considered, beginning from free-space path loss, directional beampatterns of antennas, reflection characteristics of building materials, to the polarization state of the waves and the mutual coupling of the antennas with nearby objects. Ray tracing simulations make an attempt to solving this challenge [22,23,24]. However, ultimately, this approach needs exact models of the building materials which makes it infeasible for online use in channel prediction.
In this paper, we use a data-driven approach to solve this issue. A mirror-source model is used to describe the arrival times and arrival/departure angles of the SMCs, which can be obtained by using a SLAM algorithm as in References [8,9]. We then apply Gaussian process regression (GPR) as a tool to model the amplitudes of the SMCs. That is, GPR is used to learn the joint effect of the antenna patterns and the material properties of the environment. Our prime aim is on validating this concept based on experimental data. In order to simplify the involved mathematics as much as possible, we select a scenario with a constant polarization state, no antenna coupling, and we restrict the modeling to horizontal characteristics, that is, we simplify the approach to a 2D geometric model. This reduces the amplitude model to a model that depends only on the angle of arrival (or departure) of the SMCs in the horizontal plane.
An accurate amplitude model of SMCs leads to the possibility of predicting the amplitude of SMCs very accurately. The potential use could be in performance prediction, where an accurate amplitude model of SMCs may be related to the channel capacity at a certain receiver position [25]. In References [11,26], it was shown that SLAM using position-related multipath parameters (delay and AoA) can be significantly improved in robustness if the multipath components amplitude information is also used. The contributions of this paper are the following:
  • We derive a Gaussian process regression model for describing the spatial non-stationarity of SMC amplitudes in indoor environments.
  • We validate the prediction capability of the GP based real measurements acquired in an indoor environment.
  • We show that the GP is capable of modeling the angle-dependencies of the SMC amplitudes.
The remainder of the paper is structured as follows: In Section 2 we introduce the system and signal models. Section 3 presents the GP regression and the methods for quality evaluation. Section 4 reports numerical results using real measurements. Section 5 concludes the paper.

2. System and Signal Models

We consider indoor environments where fixed anchors communicate with a mobile agent by means of radio signals. Figure 1 illustrates the considered scenario, where two physical anchor positions are shown as blue crosses at positions a 1 ( 1 ) and a 1 ( 2 ) and the agent positions p are selected along a segmented trajectory. Furthermore, some exemplary virtual anchors (VAs) are shown. The VA positions are mirror images of the physical anchor positions that are induced by reflections at flat surfaces—typically walls—and thus depend on the surrounding environment (the floor plan) [27]. The position of the k-th VA of the jth physical anchor at position a 1 ( j ) is denoted as a k ( j ) . Note that in this work we consider horizontal propagation only. For brevity, we neglect the anchor index j from now on. Also, we denote L as the set of measurement points.

2.1. Signal Model

A baseband signal s ( t ) is transmitted at carrier frequency f c . The complex envelope of the received signal then reads
r ( t ; p ) = k = 1 K α k ( p ) s ( t τ k ( p ) ) + ( s ν ) ( t ; p ) + w ( t ) ,
where the first term describes the superposition of K specular multipath components (SMCs), the second term describes the convolution of s ( t ) and the diffuse/dense multipath component (DMC) ν ( τ ; p ) and the third term is the measurement noise w ( t ) modeled as AWGN with a two-sided power spectral density of N 0 . In order to minimize inter-symbol interference, root raised cosine (RRC) pulses are considered for the transmitted signal, with pulse width T s . The energy of s ( t ) is normalized to one, that is, | s ( t ) | 2 d t = 1 .
The SMCs: The kth SMC is characterized by delay τ k ( p ) and complex amplitude α k ( p ) . The delay is modeled deterministically as a function of the (varying) agent position p , that is, τ k ( p ) = 1 c a k p = 1 c d k ( p ) , where c is the propagation speed. The complex amplitude α k ( p ) C is given as
α k ( p ) = A ( ϕ k ( p ) ) Γ k ( ϕ k ( ( p ) ) d k ( p ) exp j 2 π f c 1 c d k ( p ) ,
where A ( ϕ k ( p ) ) comprises the antenna gains of the TX and RX antennas and the path loss at a reference distance d 0 = 1 m. The antenna gains are parameterized by the direction angle ϕ k ( p ) = ( p a k ) which is defined by the line connecting the agent position p and the (virtual) anchor position a k . We assume that the antenna gains are nearly constant over the observed frequency range and that the antennas are non-dispersive. Γ k ( ϕ k ( p ) ) describes the reflection coefficient of the (flat) surface corresponding to VA a k . For the line-of-sight (LOS) component, with no involved reflection, we set the coefficient Γ 1 ( ϕ 1 ( p ) ) = 1 . Equation (2) can be decomposed into a distance-dependent factor and an angle-dependent factor. The former consists of the path-loss of a (reflected) SMC in free space and a distance-dependent phase shift. The latter describes the angle-dependence of the agent and physical anchor antenna patterns (shown in A ( ϕ k ( p ) ) ) and the dependence of the reflection coefficient Γ k ( ϕ k ( p ) ) on the incident angle. Since we assume that the agent and physical anchor are at the same height, and the orientations of the agent and physical anchor antennas are fixed throughout the measurement, the incident angle (at the reflecting surface) and the angles at the agent and anchor (used to compute antenna gains) can be represented by one angle ϕ k ( p ) = ( p a k ) as shown in Figure 1. Thus, (2) can be rewritten as.
α k ( p ) = γ k ( ϕ k ( p ) ) d k ( p ) exp j 2 π f c 1 c d k ( p ) ,
where γ k ( ϕ k ( p ) ) = | γ k ( ϕ k ( p ) ) | exp j γ k ( ϕ k ( p ) denotes the combined angle-dependence of the antenna patterns in A ( ϕ k ( p ) ) and the reflection coefficient Γ k ( ϕ k ( p ) ) .
The DMC: For the diffuse component, we use a complex Gaussian process where we assume uncorrelated scattering [28] and hence describe ν ( τ ; p ) by the auto-correlation function
K ν ( τ , u ; p ) = E ν ( τ ; p ) ν ( u ; p ) = S ν ( τ ; p ) δ ( τ u ) ,
where E · is the expectation operator and S ν ( τ ; p ) is a power delay profile (PDP), which effectively models the diffuse multipath.

2.2. Amplitude Estimation

For UWB transmission, it is reasonable to assume that the SMC signals s k ( t ) = s ( t τ k ( p ) ) are well separated in delay domain and thus approximately orthogonal. The SMC amplitude can be estimated by projecting the received signal onto the corresponding shifted pulse [29], that is,
α ^ k ( p ) = 0 T r ( t ; p ) s k ( t ) d t = α k ( p ) + 0 T s ( t ) ν ( t ; p ) + w ( t ) s k ( t ) d t
= α k ( p ) + ν k ( p ) + w k .
The variance of the estimated amplitude as derived in Appendix A is given as
var [ α ^ k ( p ) ] = E | ν k ( p ) + w k | 2 = S ν ( τ k ; p ) T p + N 0 ,
where T p = | S ( f ) | 4 d f is the effective pulse duration (which is being discussed in Appendix A). The first term on the right-hand side is a Gaussian process corresponding to the DMC and the second term is measurement noise. The estimated amplitude is thus distributed according to
α ^ k ( p ) CN α k ( p ) , S ν ( τ k ; p ) T p + N 0 .
This equation indicates that the mean SMC amplitudes α k ( p ) depend deterministically on the agent position p . The same holds for the variance of the SMC amplitudes, being a result of the interfering diffuse multipath that is characterized by the PDP S ν ( τ k ; p ) , cf. Equation (4).

2.3. GP Modeling of the SMC Amplitudes

From (3) it is seen that α k ( p ) has a well-defined distance dependence, while the angle-dependence depends on many factors that are hard to express in analytic form, for example, antenna radiation patters and reflection coefficients of building materials. We hence propose to use a GP regression model to describe the (angle-dependent part of) α k ( p ) , that is, γ k ( ϕ k ( p ) ) . A Gaussian Process (GP) model is a regression tool to analyze and predict data assuming that the data have a normal distribution. This is a loose requirement because for a large enough database, the data can often be assumed to be Gaussian distributed. Furthermore, it is evident from (8) that our data obeys a (complex) Gaussian distribution itself.
However, we will not apply the GPR model to the complex amplitudes directly, because phase coherence of SMCs may be hard to achieve, in particular, when the training data has been recorded at uncertain positions. We therefore separate the complex amplitudes into real-valued absolute values and phase values. Note from the model (8) that the absolute values will be Rician distributed [29]. This distribution is characterized by the Rician K-factor, for which in the case of K 1 , it is possible to approximate the Rician distribution with a Gaussian distribution [30], and the same holds for the phase of complex Gaussian variable. The GPR model will hence be used to model the functions
ψ ( ϕ k ( p ) ) : = | γ k ( ϕ k ( p ) ) | + n ν , k abs = | α k ( p ) | d k ( p ) + n ν , k abs
ζ ( ϕ k ( p ) ) : = γ k ( ϕ k ( p ) ) + n ν , k ph = α k ( p ) + 2 π f c 1 c d k ( p ) + n ν , k ph ,
where n ν , k abs with variance σ ν , k abs 2 and n ν , k ph with variance σ ν , k ph 2 represent the DMC of the amplitude and of the phase, respectively.

3. Gaussian Process Regression

GP regression usually consists of two steps: analyzing data obtained from measurements and predicting data at other positions that have not been measured.

3.1. GP Model

Our goal is to model the angle-dependent SMC amplitude ψ ( ϕ k ( p ) ) and phase ζ ( ϕ k ( p ) ) using a Gaussian process model. For the amplitude we have
ψ ( ϕ k ( p ) ) GP μ GP ( ϕ k ( p ) ) , c GP ( ϕ k ( p ) , ϕ k ( p ) ) ,
where μ GP ( ϕ k ( p ) ) : R 2 R denotes the mean function and c GP ( ϕ k ( p ) , ϕ k ( p ) ) : R 2 × R 2 R the covariance function (also called kernel). Equation (11) is used equivalently for the phase ζ ( ϕ k ( p ) ) . Covariance functions are a key component of GPR, since they provide the information of similarities between observations at different data points. Covariance functions are usually parameterized and these parameters, the mean function, and the noise are usually lumped together into a set called ’hyperparameters’, which we will describe later in this section.
In this work, the mean value is modeled by the constant, that is,
μ GP ( ϕ k ( p ) ) = β k abs
(and β k ph for the phase). A squared exponential kernel is used for the covariance,
c GP ( ϕ k ( p ) , ϕ k ( p ) ) = σ k abs 2 exp ( ϕ k ( p ) ϕ k ( p ) ) 2 2 a k abs 2 + σ ν , k abs 2 δ ( p p ) ,
where σ k abs is the standard deviation of the correlation kernel, a k abs is the characteristic correlation angle, and σ ν , k abs is the standard deviation accounting the DMC as shown in (9) (and σ k ph , a k ph , and σ ν , k ph for the phase). Note the abuse of notation in the latter term, which is due to the mapping of p R 2 to ϕ k ( p ) R . The Dirac function δ ( p p ) indicates that the DMC are assumed to be spatially uncorrelated. Here, we have a pair of agent positions ( p , p ) mapped to one covariance value, that is, R 2 × R 2 R .
These models will be the basis for the GP modeling of the SMC amplitudes. Together with the data, the GP model reveals the angle-dependent functions ψ ( ϕ k ( p ) ) and ζ ( ϕ k ( p ) ) as well as a variance that relates to the DMC S ν ( τ k ; p ) . The full hyper-parameter vectors of the GP model which have to be learned from the amplitude and phase data of SMC k are given as
θ k abs = [ a k abs , σ k abs , σ ν , k abs , β k abs ] T
and
θ k ph = [ a k ph , σ k ph , σ ν , k ph , β k ph ] T ,
respectively.

3.2. Prediction

We next review the generic principle of GP regression (GPR). For each SMC k, we have a training database D k = { ϕ k , y k } with i = 1 N measurements y k , i = ψ ( ϕ k ( p i ) ) + ϵ (or y k , i = ζ ( ϕ k ( p i ) ) + ϵ ), recorded at training positions P k = { p 1 , p 2 , , p N } that are mapped to angles ϕ k = [ ϕ k ( p 1 ) , , ϕ k ( p N ) ] and ϵ N ( 0 , σ ϵ 2 ) represents the measurement noise with variance σ ϵ 2 (which is connected to the power spectral density of the AWGN as described in Appendix B). The measurements are stacked into the training data vector y k = [ y k , 1 , , y k , N ] T . After the hyperparameters are learned, the mean and variance can be predicted at test position p with angle ϕ k ( p ) . The mean and the variance of the angle-dependent amplitude ψ ( ϕ k ( p ) ) (or phase ζ ( ϕ k ( p ) ) ) conditioned on the data D k and θ k abs (or θ k ph ) can be expressed in closed form [31]
E [ ψ ( ϕ k ( p ) ) | D k , θ k abs ] = μ GP ( ϕ k ( p ) ) + k ( ϕ k ( p ) , ϕ k ) T K ( ϕ k ) 1 ( y k μ ( ϕ k ) ) ,
V [ ψ ( ϕ k ( p ) ) | D k , θ k abs ] = σ k 2 + σ ν , k abs 2 k ( ϕ k ( p ) , ϕ k ) T K ( ϕ k ) 1 k ( ϕ k ( p ) , ϕ k ) ,
where V [ · ] is the variance operator and with
k ( ϕ k ( p ) , ϕ k ) = [ c GP ( ϕ k ( p ) , ϕ k ( p 1 ) ) , , c GP ( ϕ k ( p ) , ϕ k ( p N ) ) ] T ,
μ ( ϕ k ) = [ μ GP ( ϕ k ( p 1 ) ) , , μ GP ( ϕ k ( p N ) ) ] T ,
and K ( ϕ k ) R N × N with [ K ( ϕ k ) ] i , j = c GP ( ϕ k ( p i ) , ϕ k ( p j ) ) + σ ϵ 2 δ i , j , where δ i , j is the Kronecker delta. In Appendix B, the relation between the predicted variance and the PDP of the DMC and the power spectral density of the AWGN is analyzed in detail.

3.3. Learning

The GP prediction in Section 3.2 assumes the hyper-parameter vector θ k abs given in (14) (or θ k ph given in (15)) is known. If the hyper-parameter vector θ k abs is not known, it needs to be estimated from the training database D k = { ϕ k , y k } for each SMC. The joint probability density function of the observed measurements y k conditioned on the training angles ϕ k is given by a Gaussian distribution with mean and covariance both dependent on the hyper-parameter vector θ k abs (see for (12) and (13)), which we make explicit by writing μ θ k abs and K θ k abs , respectively. The parameter vector θ k abs (or θ k ph ) is estimated by maximizing the log-likelihood function, that is,
θ ^ k abs = arg max θ k abs log p ( y k | ϕ k , θ k abs ) = arg max θ k abs 1 2 log | K θ k abs ( ϕ k ) | 1 2 ( y k μ θ k abs ( ϕ k ) ) T K θ k abs ( ϕ k ) 1 ( y k μ θ k abs ( ϕ k ) ) N 2 log 2 π .
In general, this maximum likelihood estimation cannot be done in closed-form. Therefore, numerical approximation needs to be used. Since the function may be highly non-convex, θ ^ k abs (or θ ^ k phi ) has to be found by a global search over the domains of a k abs , σ k abs , and σ ν , k abs (or a k ph , σ k ph , and σ ν , k ph ).

3.4. Evaluate the Quality of Prediction

Given a prediction method, we can evaluate the quality of prediction in several ways. Perhaps, the simplest is the squared error loss. However, this quantity is sensitive to the overall scale of the target values, so it makes sense to normalize by the variance of the targets of the test cases to obtain the standard mean squared error (SMSE) [31]. Additionally, for the GPR on the absolute values of SMC amplitudes, since we produce a predictive distribution at each test input, we can evaluate the negative log probability of the target under the model. This loss should be standardized and averaged to obtain the mean standard log loss (MSLL) [31].
In order to evaluate the quality, the data set will be divided into two sets, for training and testing. Besides the training data D k , the set of test data is given by D k = { ϕ k , y k } with N measurements y k , i at positions p i . The SMSE of the predicted mean E [ ψ ( ϕ k ( p i ) ) | D k , θ k abs ] in (16) of the angle-dependent amplitude ψ ( ϕ k ( p i ) ) (or phase ζ ( ϕ k ( p i ) ) ) based on the learned parameters θ k abs (or θ k ph ) in (16) is given by [31]
SMSE abs = i = 1 N y k , i E [ ψ ( ϕ k ( p i ) ) | D k , θ k abs ] 2 i = 1 N y k , i y ¯ k , i 2 ,
where y ¯ k , i = 1 / N i = 1 N y k , i is the mean value of the measurements (similarly for the phase we have SMSE ph ). In (21), the numerator is the sum of the squared error between the measured test data y k and the predicted mean. The denominator is the sum of squared error between the measured test data y k and the mean value of the measurements y ¯ k , i . Thus, an SMSE significantly smaller than one indicates a high prediction quality.
Another measure of quality is the MSLL, which also considers the predicted variance, that is, V [ ψ ( ϕ k ( p i ) ) | D k , θ k abs ] in (17). It is calculated by averaging over the standard log loss of the amplitude ψ ( ϕ k ( p i ) ) (or phase ζ ( ϕ k ( p i ) ) ) based on the learned parameters θ k abs (or θ k ph ) for each pair in D k , that is,
MSLL = 1 N i = 1 N log p ( y k , i | ϕ k ( p i ) , θ k abs ) = 1 N i = 1 N 1 2 log ( 2 π σ n , k , i 2 ) + y k E [ ψ ( ϕ k ( p i ) ) | D k , θ k abs ] 2 2 σ n , k , i 2 ,
where σ n , k , i 2 = V [ ϕ k ( p i ) | D k , θ k abs ] + σ ϵ 2 . Large negative values of the MSLL indicate a high prediction quality.

4. Results

4.1. Experimental Setup

Figure 1 shows the laboratory room at Graz University of Technology that was used for the experimental validation. The room consists of two plaster board walls and two reinforced concrete walls (shown as black outer lines), three glass windows at the north wall (shown as thick gray lines), one white board and one metal door at the south wall (indicated by A and C , respectively). We introduce the following labels to refer to the involved reflection surfaces:
  • EPB East plaster board.
  • SW South wall.
  • WW West wall.
  • NGW North glass wall.
To conduct the channel measurements, an Ilmsens Ultra-Wide band M-sequence device [32] was used, c.f. Figure 2a. The measurement principle is correlative channel sounding [33], that is, a binary code sequence with suitable autocorrelation properties is transmitted over the channel. At the receiver, the channel impulse response is recovered using a correlation with the known code sequence. The channel sounder has one transmitter port and two receiver ports. A 12-bit M-sequence has been employed, corresponding to a sequence length of 4095 samples. This allows for an unambiguous delay window of 589.2 ns at a clock rate of 6.95 GHz. The M-sequence is modulated onto a 6.95 GHz carrier, yielding a probing signal that covers a frequency band between approx. 3.5 and 10.5 GHz. The measurement data are available online in the “MeasureMINT” data base. Please refer to: https://www.spsc.tugraz.at/databases-and-tools/uwb-indoor-channel-experimental-data.html.
Each of the ports was connected to a dipole coin antenna as shown in Figure 2b. According to Reference [34], the coin antenna has a very wide bandwidth ranging from 3 to 9 GHz. It also has a nearly isotropic radiation pattern in the horizontal plane.
We used the two receiver ports as anchors and placed their antennas at fixed positions a 1 ( 1 ) and a 1 ( 2 ) . The transmitter port is connected to another antenna that was moved along a trajectory with 595 points p , as shown in the Figure 1, to obtain the same number of channel measurements. All antennas were mounted on tripods at the same height, therefore only the co-polarized, azimuth radiation pattern of the antenna has an impact on the data, c.f. Equation (2). The raw measurements at the receiver ports were filtered with an RRC pulse with center frequency 6.95 GHz, roll-off factor 0.5 and bandwidth 1 / T s = 2 GHz to obtain the received signals corresponding to the model in (1).

4.2. Measurement Pre-Processing

The agent trajectory is divided into 7 segments shown with different colors. The spacing between trajectory points was 2 cm for segments 1 and 2, and 4 cm elsewhere. During the experiment, the exact positions of the trajectory points were unknown. Hence, in a pre-processing step, we used the SLAM algorithm [8,10] to estimate the trajectory positions L = { p } and the locations of the VAs and physical anchors { a k } . The estimated trajectory and VAs are plotted in Figure 1. For Anchor 1, there are the VAs a 2 ( 1 ) , a 4 ( 1 ) and a 5 ( 1 ) —indicated here by additional superscript index ( 1 ) —, corresponding to single reflections on the EPB, WW and NGW respectively. For Anchor 2, there are the VAs a 2 ( 2 ) , a 3 ( 2 ) and a 4 ( 2 ) , indicated by index ( 2 ) , corresponding to single reflections on the EPB, SW and WW, respectively.
Given the estimated positions of transmitters and receivers, the expected delays, that is, τ k ( p ) = 1 c a k p are deduced. Then, the SMC amplitudes α ^ k ( p ) are estimated from the received signals, using these delays in (5). The use of (5) for estimating the SMC amplitudes relies on the assumption that individual SMCs do not overlap, reducing the interfering multipath to the DMC ν k ( p ) . We remove overlapping SMCs from our data set if two delays have | τ k ( p ) τ l ( p ) | T p 2 , for k l , yielding the subsets P k L of trajectory points. Finally, the distance dependence will be removed from the SMC amplitudes to get the normalized amplitude data ψ ( ϕ k ( p ) ) and phase data ζ ( ϕ k ( p ) ) according to Equations (9) and (10). Note that the trajectory positions and VA positions are mapped to the direction angles ϕ k ( p ) = ( p a k ) for the purpose of GPR.
Several datasets are hence deduced, which are denoted as D k = ϕ k ( p ) , y k ( p ) p P k , containing the direction angles ϕ k ( p ) and the data points y k ( p ) . The data points are either the absolute values ψ ( ϕ k ( p ) ) or the phases ζ ( ϕ k ( p ) ) . The index k relates here to the different SMCs and the two anchors. Specifically, for Anchor 1, we analyze data sets corresponding to the LOS, the EPB, and the NGW, for Anchor 2, we analyze the EPB and the SW.
On the basis of these data sets, the GP model is learned using the Matlab function fitrgp (Matlab version R2018a; Statistics and Machine Learning Toolbox). Prediction is done using the function predict. The regression result will be discussed to see if it is sensible to use GP regression to model the SMC amplitudes. We also evaluate the GPR performance numerically with the SMSE and the MSLL in Section 3.4.

4.3. GPR of SMC Amplitudes

Table 1 lists the learned hyperparameters θ k abs (or θ k ph ) obtained from the measurements. It is shown that the characteristic correlation angle of the LOS GP amplitude has the highest angle correlation, that is, the estimated LOS amplitude shows a small variation over a wide range of angles. This might come from two factors, (i) the influence of DMC on the LOS amplitudes is rather small shown by the DMC standard deviation σ ν , k abs in Table 1, and (ii) the estimated LOS amplitudes only depend on the antenna pattern A ( ϕ k ( p ) ) , while the reflection coefficient Γ 1 ( ϕ 1 ( p ) ) is defined to be one. The characteristic correlation angles of the SMC GP amplitudes are significantly smaller, that is, the estimated SMC amplitudes have a much larger variation over the range of angles. This is a strong indicator that the SMC amplitudes show an angle-dependent reflection coefficient Γ k ( ϕ k ( p ) ) and the influence of DMC is larger, which is in agreement with increased DMC standard deviation σ ν , k abs in Table 1. All these observations are supported by results shown in the figures below. In the following, all estimated and modeled SMC amplitudes are normalized as defined by Equation (9).
Figure 3 shows the regression results for the LOS amplitude associated with Anchor 1 and some SMCs associated with Anchors 1 and 2. Particularly, it shows the estimated SMC amplitudes (colored markers) for the individual segments, the predicted mean given in (16) (red solid line), and the predicted standard deviation given in (17) (black dashed and blue dash-dotted lines) as a function of the angle of departure ϕ 2 ( p ) . The deviation is shown as ± 2 σ , representing the upper and lower limits that contain 95% of the data points.
In Figure 3a, the LOS amplitudes are shown. It can be observed that the predicted mean value of the LOS GP amplitude only slowly changes with the angle ϕ 1 ( p ) according to the antenna pattern A ( ϕ 1 ( p ) ) . The predicted standard deviation is rather small since the learned DMC standard deviation and the characteristic correlation angle are large.
Figure 3b–e show clearly that the predicted mean values of the SMC GP amplitudes vary significantly with the angle ϕ k ( p ) . Those variations originate mostly from the angle-dependency of the reflection coefficient Γ k ( ϕ k ( p ) ) , because the angle-dependency of the antenna pattern A ( ϕ k ( p ) ) is slow over the entire observation angle. As expected, the predicted standard deviation is much larger than the one for the LOS amplitude. (A more detailed verification of the variance regression performance is found in Section 4.3.2.) For the SW (at a 3 ( 2 ) ) and the NGW (at a 5 ( 1 ) ), the effect of varying reflection coefficient due to different building materials can be seen; c.f. Figure 3d,e, respectively. In Figure 3d, sections A , C and B correspond to the whiteboard, the metal door and a small plasterboard section in between, c.f. Figure 1. It can be seen that, for angles between 260 to 230 degrees, the fluctuations are small because the SMC relates to the strong reflection by the metallic whiteboard. For angles above 230 degrees, the amplitude fluctuates severely, because the point-of-reflection moves to sections B and C comprising of other materials. Also note that the variance of the data points fluctuates. In general, data from positions which are quite far away from the VA a 3 ( 2 ) have larger variances. This is partly the result of a more significant impact of the DMC on those data but it is also the effect of multiplying the DMC with the longer travel distance to obtain the normalized SMC amplitude. In Figure 3e, the amplitude fluctuations due to varying reflection materials can be seen again. The various materials are noted accordingly, whereas sections A, C and E are glass windows. Even with the same material, separate angle ranges are covered where the normalized amplitude, and thus the reflection coefficient, is significantly different, for example, for lower angles, an increased coefficient can be seen.

4.3.1. Predictability

In order to test the ability to predict the SMC amplitudes using GPR, we divide the data set into two subsets, for training and testing. It is desirable for the training data to cover the largest possible range of angles. Therefore, we choose segments 4, 5, 6 for training, and 1, 2, 3, 7 for testing (c.f. Figure 1).
Figure 4a–d show the prediction results using the normalized amplitudes for EPB (at a 2 ( 1 ) ), LOS (at a 1 ( 1 ) ), NGW (at a 5 ( 1 ) ) and SW (at a 3 ( 2 ) ), respectively. Again, the red solid bold line depicts the predicted mean, the black and blue solid lines show the upper and lower ± 2 σ limits (comprising 95 % of the data points), plotted over the positions along the trajectory. It is shown that the predicted mean of the SMC amplitudes coincides very well with variations of the estimated SMC amplitudes. Moreover, the predicted upper and lower limits contain most of the SMC amplitudes obtained from measurements for the test data.
Table 2 shows the prediction quality evaluated with the SMSE and MSLL as defined in Section 3.4. The predictability for the LOS amplitude is clearly very reliable, indicated by the small SMSE value and large negative MSLL values. The predictions of the SMC amplitudes for EPB, NGW and SW are still reasonable good in comparison to the LOS amplitude given that the SMSE values are all smaller than 1 and MSLL values are all negative. However, the amplitude prediction for EPB outperforms the predictions for the NGW and SW. A possible explanation might be that the DMC standard deviation is smaller (see for Table 1) for the estimated SMC amplitudes of EPB and the material of the reflective surface is more homogeneous (only plasterboard instead of a combination of plasterboard, metal, and glass).

4.3.2. Variance Verification

In this section, we will validate the capability of the GPR to model the power of the DMC by means of the predicted variance, cf. Equation (17). Based on the estimated SMC of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), the empirical mean and the empirical standard deviation are calculated over a sliding window along the trajectory positions, that is,
μ meas , m = 1 N win i = m N win + m | α ^ k |
and
σ meas , m = 1 N win i = m N win + m ( | α ^ k | μ meas , m ) 2 ,
where m is the index of the agent position and N win is the sliding window size, which is chosen to be N win = 10 .
Figure 5a,b compare the empirical mean in (23) and the empirical standard deviation in (24) with the predicted mean (given by (16)) and the predicted standard deviation (square-root of (17)) of the GP along the agent trajectory. In Figure 5a, the blue solid line shows the estimated SMC amplitudes, the red solid line shows the empirical mean μ meas , m and the red dashed lines the empirical standard deviation ( ± 2 σ meas , m ) of the estimated amplitudes. In Figure 5b, the blue solid line shows the estimated SMC amplitudes, the green solid line shows the predicted mean and the green dashed lines the predicted standard deviation ( ± 2 σ GP ) of the SMC GP amplitude. Both, the mean and the standard deviation are in a good agreement. However, the predicted standard deviation (variance) in (17) has constant values for the DMC standard deviation (see Table 1), independent of the measurements used for prediction, and only changes deterministically with the angle via the first term of the correlation kernel in (13).

4.4. GPR of SMC Phases

To show the applicability of the GPR for phase prediction of the SMC phases, we exemplaryly analyze the estimated SMC amplitudes related to EPB and NGW.
Figure 6 shows the results, which indicates the phase coherence of the SMCs. The respective hyperparameters obtained from the GPR are also shown in Table 1. From this Table, it is shown that the phase of the NGW component has the largest angle dependence, given that its characteristic correlation angle is the smallest.

Predictability

Similar to Section 4.3.1, we separate the measurements into a training set containing data from Segments 4, 5, 6 and a test set containing data from Segments 1, 2, 3, and 7. Figure 7 shows the prediction results. It can be seen that the prediction of the mean and variance is reasonably good. Table 2 shows the prediction quality via SMSE and MSLL. The SMSE values for EPB a 2 ( 1 ) , EPB a 2 ( 2 ) and NGW a 5 ( 5 ) are all significantly smaller than 1, indicating that the GPR prediction is much better than the one of a trivial method. Similarly, the MSLLs are all negative, also showing good prediction quality. This evaluation indicates that a consideration of the carrier phase may become possible for multipath-based tracking filters, if GPR is used for modeling the phase shifts of the SMCs.

5. Conclusions

We demonstrate the applicability of Gaussian process regression (GPR) for the large-scale amplitude and phase modeling of specular multipath components (SMCs). UWB or mm-wave radio systems provide sufficient temporal and/or angular resolution to resolve such SMCs at the receiver side. The GPR model describes the amplitude and phase fluctuations of SMCs due to changing reflection characteristics of building materials and due to the angle-dependent gain of the TX and RX antennas. Our experimental results—performed with UWB channel measurements—show that the mean amplitude and phase information of SMCs can be predicted reliably, while some uncertainty results from the dense multipath component (DMC) that interferes with the SMCs. The variance predicted by the GPR model reflects the impact of this DMC.
Large-scale predictability of the radio channel is achieved through the combination of the GPR model (to describe the SMC amplitudes and phase shifts) and a geometric model of the SMCs to describe the arrival times and angles. That is, based on a set of training data and hyperparameters of the models learned from it, these core features of the channel impulse response can be predicted at any point in the surveyed environment. We obtain a new level of environment-awareness, useful to predict performance characteristics of wireless networks over large spatial scales.

Author Contributions

Conceptualization, K.W. and E.L.; Formal analysis, A.H.N., E.L. and M.R.; Investigation, A.H.N., M.R. and E.L.; Software, Validation, Visualization, A.H.N. and M.R.; Supervision, K.V.N. and K.W.; Writing—original draft, A.H.N.; Writing—review & editing, M.R., E.L. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge the financial support from the OEAD (Autrian Exchange Service) (for funding Anh Hong Nguyen to conduct research at Graz University of Technology), from Graz University of Technology (for funding Michael Rath within the Research Progam “Dependable Things”), and from the Austrian Science Fund (FWF) (for funding Erik Leitinger under the Schroedinger Fellowship program, grant number J4027).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Variance of νk

Here, the variance of DMC ν k ( p ) within one estimated SMC amplitude is derived. For the sake of readability, we omit the agent position p . The DMC of an estimated SMC amplitude is given as
ν k = s ( τ λ ) ν ( λ ) s ( τ τ k ) d τ d λ
where the energy of the signal s ( t ) is assumed to be one, that is, | s ( t ) | 2 d t = 1 . Using (A1), the variance of the DMC of an estimated SMC amplitude
E | ν k | 2 = E s ( τ λ ) ν ( λ ) s ( τ τ k ) s ( τ λ ) ν ( λ ) s ( τ τ k ) d τ d λ d τ d λ = s ( τ λ ) s ( τ λ ) s ( τ τ k ) s ( τ τ k ) S ν ( λ ) d τ d τ d λ = s ˜ ( τ k λ ) s ˜ ( τ k λ ) S ν ( λ ) d λ
where s ˜ ( t ) = s ( τ ) s ( t τ ) d τ . In the case a large bandwidth (UWB) is assumed, (A2) can be approximated by
E | ν k | 2 S ν ( τ k ) | s ˜ ( t ) | 2 d t = S ν ( τ k ) | S ( f ) | 4 d f
where S ( f ) is the Fourier transform of s ( t ) and T p = | S ( f ) | 4 d f is the effective pulse duration.
  • Assuming a signal s ( t ) with block-spectrum, the effective pulse duration is given as T p = | S ( f ) | 4 d f = T s , where T s is the Nyquist sampling time.
  • For root-raised-cosine, raised-cosine, or similar waveforms having a rectangular spectrum with tapered ends, T p is less than the Nyquist sampling time T s , under common symmetry conditions. For example, a raised-cosine pulse results in T p = | S ( f ) | 4 d f = 4 β 4 T s , where β is the roll-factor.

Appendix B. Predicted Variance

The predicted variance V [ y ( ϕ k ( p ) ) | D k , θ k ] at position p approximates the variance of the SMC amplitude values, which relates to the DMC term in (8) for position p . Assuming the Gaussian approximation to hold for | α ^ k ( p ) | , we have
V [ ψ ( ϕ k ( p ) ) | D k , θ k ] 1 2 S ν ( τ k ; p ) T p + N 0 d k ( p ) 2 .
While the PDP S ν ( τ k ; p ) decreases with delay τ k = 1 c d k ( p ) , the squared distance d k ( p ) 2 will be increasing, implying a counteracting effect. The left-hand side is independent of the distance. We can thus re-write the equation to illustrate the shape of the PDP as a function of τ k at some angle ϕ k ,
S ν ( τ k ; ϕ k ) T p 2 τ k 2 c 2 V [ ψ ( ϕ k ) | D k , θ k ] N 0 .
This result indicates a reciprocal squared decay in contrast to the usually assumed exponential decay.
On the other hand, from Appendix A, by neglecting measurement noise, we have:
E | ν k ( p ) + w k ( p ) | 2 E | ν k ( p ) | 2 = T p S ν ( τ k ; p )
So, the predicted variance obtained from GPR can be validated by comparing with the DMC power, c.f. Equations (A4) and (A5).

References

  1. Boccardi, F.; Heath, R.W.; Lozano, A.; Marzetta, T.L.; Popovski, P. Five disruptive technology directions for 5G. IEEE Commun. Mag. 2014, 52, 74–80. [Google Scholar] [CrossRef] [Green Version]
  2. Fettweis, G.P. The Tactile Internet: Applications and Challenges. IEEE Veh. Technol. Mag. 2014, 9, 64–70. [Google Scholar] [CrossRef]
  3. Simsek, M.; Aijaz, A.; Dohler, M.; Sachs, J.; Fettweis, G. 5G-Enabled Tactile Internet. IEEE J. Sel. Areas Commun. 2016, 34, 460–473. [Google Scholar] [CrossRef] [Green Version]
  4. Di Taranto, R.; Muppirisetty, S.; Raulefs, R.; Slock, D.; Svensson, T.; Wymeersch, H. Location-Aware Communications for 5G Networks: How location information can improve scalability, latency, and robustness of 5G. IEEE Signal Process. Mag. 2014, 31, 102–112. [Google Scholar] [CrossRef] [Green Version]
  5. Muppirisetty, L.; Svensson, T.; Wymeersch, H. Spatial Wireless Channel Prediction under Location Uncertainty. IEEE Trans. Wirel. Commun. 2016, 15, 1031–1044. [Google Scholar] [CrossRef] [Green Version]
  6. Fröhle, M.; Charalambous, T.; Nevat, I.; Wymeersch, H. Channel Prediction With Location Uncertainty for Ad Hoc Networks. IEEE Trans. Signal Inf. Process. Netw. 2018, 4, 349–361. [Google Scholar] [CrossRef] [Green Version]
  7. Gentner, C.; Jost, T.; Wang, W.; Zhang, S.; Dammann, A.; Fiebig, U.C. Multipath Assisted Positioning with Simultaneous Localization and Mapping. IEEE Trans. Wirel. Commun. 2016, 15, 6104–6117. [Google Scholar] [CrossRef] [Green Version]
  8. Leitinger, E.; Meyer, F.; Hlawatsch, F.; Witrisal, K.; Tufvesson, F.; Win, M.Z. A Belief Propagation Algorithm for Multipath-Based SLAM. IEEE Trans. Wirel. Commun. 2019, 18, 5613–5629. [Google Scholar] [CrossRef] [Green Version]
  9. Mendrzik, R.; Meyer, F.; Bauch, G.; Win, M.Z. Enabling Situational Awareness in Millimeter Wave Massive MIMO Systems. IEEE J. Sel. Top. Signal Process. 2019, 13, 1196–1211. [Google Scholar] [CrossRef]
  10. Leitinger, E.; Meyer, F.; Tufvesson, F.; Witrisal, K. Factor graph based simultaneous localization and mapping using multipath channel information. In Proceedings of the 2017 IEEE International Conference on Communications Workshops (ICC Workshops), Paris, France, 21–25 May 2017; pp. 652–658. [Google Scholar] [CrossRef]
  11. Leitinger, E.; Grebien, S.; Witrisal, K. Multipath-Based SLAM Exploiting AoA and Amplitude Information. In Proceedings of the 2019 IEEE International Conference on Communications Workshops (ICC Workshops), Shanghai, China, 20–24 May 2019; pp. 1–7. [Google Scholar] [CrossRef]
  12. Leitinger, E.; Meissner, P.; Rudisser, C.; Dumphart, G.; Witrisal, K. Evaluation of Position-Related Information in Multipath Components for Indoor Positioning. IEEE J. Sel. Areas Commun. 2015, 33, 2313–2328. [Google Scholar] [CrossRef] [Green Version]
  13. Shahmansoori, A.; Garcia, G.E.; Destino, G.; Seco-Granados, G.; Wymeersch, H. Position and Orientation Estimation Through Millimeter-Wave MIMO in 5G Systems. IEEE Trans. Wirel. Commun. 2018, 17, 1822–1835. [Google Scholar] [CrossRef] [Green Version]
  14. Mendrzik, R.; Wymeersch, H.; Bauch, G.; Abu-Shaban, Z. Harnessing NLOS Components for Position and Orientation Estimation in 5G Millimeter Wave MIMO. IEEE Trans. Wirel. Commun. 2019, 18, 93–107. [Google Scholar] [CrossRef]
  15. Gustafsson, F.; Gunnarsson, F. Mobile positioning using wireless networks: Possibilities and fundamental limitations based on available wireless network measurements. IEEE Signal Process. Mag. 2005, 22, 41–53. [Google Scholar] [CrossRef]
  16. Kumar, S.; Gil, S.; Katabi, D.; Rus, D. Accurate indoor localization with zero start-up cost. In Proceedings of the 20th Annual International Conference on Mobile Computing and Networking, Maui, HI, USA, 7–11 September 2014. [Google Scholar]
  17. Meissner, P.; Leitinger, E.; Witrisal, K. UWB for Robust Indoor Tracking: Weighting of Multipath Components for Efficient Estimation. Wirel. Commun. Lett. IEEE 2014, 3, 501–504. [Google Scholar] [CrossRef]
  18. Witrisal, K.; Meissner, P.; Leitinger, E.; Shen, Y.; Gustafson, C.; Tufvesson, F.; Haneda, K.; Dardari, D.; Molisch, A.F.; Conti, A.; et al. High-Accuracy Localization for Assisted Living: 5G systems will turn multipath channels from foe to friend. IEEE Signal Process. Mag. 2016, 33, 59–70. [Google Scholar] [CrossRef]
  19. Poutanen, J.; Salmi, J.; Haneda, K.; Kolmonen, V.; Vainikainen, P. Angular and Shadowing Characteristics of Dense Multipath Components in Indoor Radio Channels. IEEE Trans. Antennas Propag. 2011, 59, 245–253. [Google Scholar] [CrossRef]
  20. Karedal, J.; Tufvesson, F.; Czink, N.; Paier, A.; Dumard, C.; Zemen, T.; Mecklenbrauker, C.F.; Molisch, A.F. A geometry-based stochastic MIMO model for vehicle-to-vehicle communications. IEEE Trans. Wirel. Commun. 2009, 8, 3646–3657. [Google Scholar] [CrossRef] [Green Version]
  21. Virk, U.T.; Haneda, K.; Wagen, J. Dense multipath components add-on for COST 2100 channel model. In Proceedings of the 2015 9th European Conference on Antennas and Propagation (EuCAP), Lisbon, Portugal, 13–17 April 2015; pp. 1–5. [Google Scholar]
  22. Oestges, C.; Clerckx, B.; Raynaud, L.; Vanhoenacker-Janvier, D. Deterministic channel modeling and performance simulation of microcellular wide-band communication systems. IEEE Trans. Veh. Technol. 2002, 51, 1422–1430. [Google Scholar] [CrossRef]
  23. Athanasiadou, G.E.; Nix, A.R.; McGeehan, J.P. A microcellular ray-tracing propagation model and evaluation of its narrow-band and wide-band predictions. IEEE J. Sel. Areas Commun. 2000, 18, 322–335. [Google Scholar] [CrossRef] [Green Version]
  24. Kanatas, A.G.; Kountouris, I.D.; Kostaras, G.B.; Constantinou, P. A UTD propagation model in urban microcellular environments. IEEE Trans. Veh. Technol. 1997, 46, 185–193. [Google Scholar] [CrossRef]
  25. Nguyen, H.A.; Rath, M.; Leitinger, E.; Hinteregger, S.; Nguyen, V.K.; Witrisal, K. Channel Capacity Analysis of Indoor Environments for Location-Aware Communications. In Proceedings of the 2016 IEEE Globecom Workshops (GC Wkshps), Washington, DC, USA, 4–8 December 2016; pp. 1–6. [Google Scholar] [CrossRef]
  26. Leitinger, E.; Grebien, S.; Li, X.; Tufvesson, F.; Witrisal, K. On the Use of MPC Amplitude Information in Radio Signal Based SLAM. In Proceedings of the IEEE Statistical Signal Processing Workshop (SSP), Freiburg, Germany, 10–13 June 2018; pp. 633–637. [Google Scholar]
  27. Borish, J. Extension of the Image Model to arbitrary Polyhedra. JASA 1984, 75, 1827–1836. [Google Scholar] [CrossRef]
  28. Bello, P. Characterization of Randomly Time-Variant Linear Channels. IEEE Trans. Commun. Syst. 1963, 11, 360–393. [Google Scholar] [CrossRef] [Green Version]
  29. Meissner, P.; Witrisal, K. Analysis of position-related information in measured UWB indoor channels. In Proceedings of the 2012 6th European Conference on Antennas and Propagation (EUCAP), Prague, Czech Republic, 26–30 March 2012; pp. 6–10. [Google Scholar]
  30. Durgin, G. Space-time Wireless Channels, 1st ed.; Prentice Hall Press: Upper Saddle River, NJ, USA, 2002. [Google Scholar]
  31. Rasmussen, C.; Williams, C. Gaussian Processes for Machine Learning; Adaptive Computation and Machine Learning; MIT Press: Cambridge, MA, USA, 2006; p. 248. [Google Scholar]
  32. Sachs, J.; Herrmann, R.; Kmec, M.; Helbig, M.; Schilling, K. Recent Advances and Applications of M-Sequence based Ultra-Wideband Sensors. In Proceedings of the 2007 IEEE International Conference on Ultra-Wideband, Singapore, 24–26 September 2007; pp. 50–55. [Google Scholar] [CrossRef]
  33. Molisch, A.F. Wireless Communications, 2nd ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2010. [Google Scholar]
  34. Krall, C. Signal Processing for Ultra Wideband Transceivers. Ph.D. Thesis, Graz University of Technology, Graz, Austria, 2008. [Google Scholar]
Figure 1. Floor plan of the evaluation scenario. Bold black lines denote walls, thick gray lines represent glass windows, other lines illustrate other materials. Two blue crosses represent the physical anchors; orange circles denote virtual anchors (VAs) which were considered in the experimental validation. An agent moves along a trajectory segmented into seven parts indicated with distinct colors. Capital letters (with or without mark * or †) refer to sub-segments of different materials along each wall.
Figure 1. Floor plan of the evaluation scenario. Bold black lines denote walls, thick gray lines represent glass windows, other lines illustrate other materials. Two blue crosses represent the physical anchors; orange circles denote virtual anchors (VAs) which were considered in the experimental validation. An agent moves along a trajectory segmented into seven parts indicated with distinct colors. Capital letters (with or without mark * or †) refer to sub-segments of different materials along each wall.
Applsci 10 05216 g001
Figure 2. Photos of the measurement hardware showing (a) the Ilmsens ultra-wideband (UWB) M-Sequence channel sounder and (b) a coin UWB dipole antenna.
Figure 2. Photos of the measurement hardware showing (a) the Ilmsens ultra-wideband (UWB) M-Sequence channel sounder and (b) a coin UWB dipole antenna.
Applsci 10 05216 g002
Figure 3. GP regression of SMC amplitude for (a) LOS amplitude associated with anchor 1 at a 1 ( 1 ) , (b) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (c) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (d) the SMC amplitude of VA at a 3 ( 2 ) associated with anchor 2 (south wall), (e) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Figure 3. GP regression of SMC amplitude for (a) LOS amplitude associated with anchor 1 at a 1 ( 1 ) , (b) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (c) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (d) the SMC amplitude of VA at a 3 ( 2 ) associated with anchor 2 (south wall), (e) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Applsci 10 05216 g003
Figure 4. GP prediction of SMC amplitude for (a) the LOS amplitude associated with anchor 1 at a 1 ( 1 ) , (b) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window), (d) the SMC amplitude of VA a 3 ( 2 ) associated with anchor 2 (south wall).
Figure 4. GP prediction of SMC amplitude for (a) the LOS amplitude associated with anchor 1 at a 1 ( 1 ) , (b) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window), (d) the SMC amplitude of VA a 3 ( 2 ) associated with anchor 2 (south wall).
Applsci 10 05216 g004
Figure 5. Estimated alpha compared to (a) measurement and (b) GP statistics ( a 2 ( 1 ) ).
Figure 5. Estimated alpha compared to (a) measurement and (b) GP statistics ( a 2 ( 1 ) ).
Applsci 10 05216 g005
Figure 6. GP regression of the SMC phases for (a) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (b) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Figure 6. GP regression of the SMC phases for (a) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (b) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Applsci 10 05216 g006
Figure 7. GP prediction of the SMC phases for (a) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (b) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Figure 7. GP prediction of the SMC phases for (a) the SMC amplitude of VA at a 2 ( 1 ) associated with anchor 1 (east plaster board wall), (b) the SMC amplitude of VA at a 2 ( 2 ) associated with anchor 2 (east plaster board wall), (c) the SMC amplitude of VA at a 5 ( 1 ) associated with anchor 1 (north glass window).
Applsci 10 05216 g007
Table 1. Parameters of the Gaussian Process (GP) model for the absolute value of the specular multipath component (SMC) amplitudes.
Table 1. Parameters of the Gaussian Process (GP) model for the absolute value of the specular multipath component (SMC) amplitudes.
MPCConstant Mean β k abs DMC Stdv. σ ν , k abs Char. Corr. Angle a k abs GP Stdv. σ k abs
or β k ph [Linear or Rad]or σ ν , k ph [Linear or Rad]or a k ph [Deg.]or σ k ph [Linear or Rad]
abs. valueEPB—anchor 10.0620.0092.910.013
LOS—anchor 10.0990.00343.110.035
NGW—anchor 10.0750.0141.020.02
SW—anchor 20.090.0132.170.03
phaseEPB—anchor 1−0.9560.3879.831.03
EPB—anchor 2−12.740.77985.9915.57
NGW—anchor 1−3.790.8263.371.71
Table 2. Prediction quality.
Table 2. Prediction quality.
MPCSMSEMSLL
AmplitudeEPB—anchor 10.35−0.381
LOS—anchor 10.014−2.099
NGW—anchor 10.46−0.175
SW—anchor 20.762−0.242
PhaseEPB—anchor 10.143−0.834
EPB—anchor 20.058−1.25
NGW—anchor 10.294−0.54

Share and Cite

MDPI and ACS Style

Nguyen, A.H.; Rath, M.; Leitinger, E.; Nguyen, K.V.; Witrisal, K. Gaussian Process Modeling of Specular Multipath Components. Appl. Sci. 2020, 10, 5216. https://doi.org/10.3390/app10155216

AMA Style

Nguyen AH, Rath M, Leitinger E, Nguyen KV, Witrisal K. Gaussian Process Modeling of Specular Multipath Components. Applied Sciences. 2020; 10(15):5216. https://doi.org/10.3390/app10155216

Chicago/Turabian Style

Nguyen, Anh Hong, Michael Rath, Erik Leitinger, Khang Van Nguyen, and Klaus Witrisal. 2020. "Gaussian Process Modeling of Specular Multipath Components" Applied Sciences 10, no. 15: 5216. https://doi.org/10.3390/app10155216

APA Style

Nguyen, A. H., Rath, M., Leitinger, E., Nguyen, K. V., & Witrisal, K. (2020). Gaussian Process Modeling of Specular Multipath Components. Applied Sciences, 10(15), 5216. https://doi.org/10.3390/app10155216

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

Article Metrics

Back to TopTop