Next Article in Journal
Comparing YOLOv3, YOLOv4 and YOLOv5 for Autonomous Landing Spot Detection in Faulty UAVs
Next Article in Special Issue
MineBL: A Battery-Free Localization Scheme with Binocular Camera for Coal Mine
Previous Article in Journal
A Microfluidic Dielectric Spectroscopy System for Characterization of Biological Cells in Physiological Media
Previous Article in Special Issue
Human Activity Recognition for Indoor Localization Using Smartphone Inertial Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Amplitude Modeling of Specular Multipath Components for Robust Indoor Localization

1
School of Electronics and Telecommunications, Hanoi University of Science and Technology, Hanoi 11615, Vietnam
2
Signal Processing and Speech Communication Lab, Graz University of Technology, 8010 Graz, Austria
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(2), 462; https://doi.org/10.3390/s22020462
Submission received: 8 November 2021 / Revised: 24 December 2021 / Accepted: 30 December 2021 / Published: 8 January 2022

Abstract

:
Ultra-Wide Bandwidth (UWB) and mm-wave radio systems can resolve specular multipath components (SMCs) from estimated channel impulse response measurements. A geometric model can describe the delays, angles-of-arrival, and angles-of-departure of these SMCs, allowing for a prediction of these channel features. For the modeling of the amplitudes of the SMCs, a data-driven approach has been proposed recently, using Gaussian Process Regression (GPR) to map and predict the SMC amplitudes. In this paper, the applicability of the proposed multipath-resolved, GPR-based channel model is analyzed by studying features of the propagation channel from a set of channel measurements. The features analyzed include the energy capture of the modeled SMCs, the number of resolvable SMCs, and the ranging information that could be extracted from the SMCs. The second contribution of the paper concerns the potential applicability of the channel model for a multipath-resolved, single-anchor positioning system. The predicted channel knowledge is used to evaluate the measurement likelihood function at candidate positions throughout the environment. It is shown that the environmental awareness created by the multipath-resolved, GPR-based channel model yields higher robustness against position estimation outliers.

1. Introduction

The vision of an Internet-of-Things (IoT) has been discussed extensively, in recent years, covering end-user applications in a wide range of application domains, including industrial [1], health-care and assisted living [2], smart cities, and smart mobility [3]. Similarly, 5G and beyond-5G wireless networks are currently under development, aiming for instance at an end-to-end latency in the range of 1 ms, high user densities, and unprecedented data rates per volume for a heterogeneous set of end-user scenarios [4,5,6]. Location information is a critical component of many of these applications [7], creating context-awareness with different service quality requirements [8]. The relevant performance criteria extend far beyond the required level of accuracy. For assisted and autonomous driving, for example, location information is of mission-critical importance, making it a prime example of a high-end location-based application with a need for 100% availability [9,10,11,12]. Robustness, security, and privacy issues in IoT applications have recently been discussed in [13]. High-accuracy location information is seen as a key enabler for location-based services in the IoT and in 5G networks [12,14].
Much work has been done to address the various issues involved in accurate and robust position estimation, demonstrating the heterogeneity of scenarios and applications as discussed before. RSS-based techniques are very suitable for positioning based on signals-of-opportunity, because it does not require dedicated time stamping and synchronization. RSS measurements can be used to compute the position from path-loss-based range estimates [15,16], from ray-tracing based propagation predictions [17,18], from fingerprinting/machine-learning-based methods [18,19,20], and from tomography-inspired, device-free measurements [21]. The main disadvantage of this approach is the low accuracy and low robustness compared with other measurement principles [22], in particular, when multipath fading and shadowing are taken into account.
Time-of-flight (ToF) measurements and multilateration are usually the basis for robust, high-accuracy positioning systems, including global navigation satellite systems (GNSS) [23]. It is well known that the ToF ranging error variance scales reciprocally with the squared bandwidth and the signal-to-noise ratio (SNR) [24], which motivates the advantage of having a high signal bandwidth available. The SNR depends only weakly on the range (ensured e.g., due to power control mechanisms in the network), however, multipath will have a severe influence [25,26,27,28]. These basic insights have been used to study the positioning performance of wireless networks [26,29], showing the influence of the geometry (and number) of the nodes that are being involved in the position computation.
The availability of sufficiently accurate location information is a facilitator for exploiting “location awareness” to enhance the performance of the radio access network [30]. Location-awareness leverages geometric information to optimize algorithms and protocols of a wireless communication system at different layers of the protocol stack, to satisfy the needs for higher throughput, lower latency, higher robustness, and other performance indicators [30]. At higher layers of the stack, location information has been employed directly for routing algorithms [31,32]. At lower layers, a channel quality metric (CQM) is often introduced, which maps between a network performance indicator, e.g., the channel capacity or the signal-to-interference-plus-noise ratio (SINR) and the position, yielding a so-called “radio environment map” (REM) [33]. The REM is obtained from a database, which is acquired from measurements conducted beforehand. Crowd-sensing has been proposed for the data collection, where data covering a wide spatial area is being collected by a large number of agents navigating through an environment [34]. A regression tool is needed to handle the large database and perform predictions. Gaussian process regression (GPR) is of particular interest because it enables the development of inference algorithms whose complexity does not grow with the number of observations [34,35].
In [36,37], situational awareness has been defined, a concept where the location information of radio nodes is collected together with information about the propagation environment. The environment model is comprised of specular multipath components (SMCs), their angle-of-arrivals (AOAs), angle-of-departures (AODs), the corresponding delays as well as points of reflections. This information can be deduced from channel measurements by means of a Simultaneous Localization And Mapping (SLAM) algorithm [38,39], where SMCs are modeled by virtual anchors (VAs), i.e., mirror images of the anchor w.r.t. different reflecting surfaces. Using these VA positions together with a data association algorithm, multipath-assisted, single-anchor localization and tracking algorithms can be formulated. Such algorithms exploit the VA-based environment model, achieving centimeter accuracy when UWB signals are used [14,38,40,41]. Crowd-sourcing of such map information has been studied in [42]. Situational awareness is of specific interest for the beam-prediction in mmWave radio systems [37].
In this paper, we combine the VA-based geometric environment model and a GPR model. The former is used to predict delay, AOA, and AOD of SMCs, while the latter is used to predict the amplitudes of SMCs, as described in [43]. Extending over our previous paper [43], the following main contributions are made:
  • We formulate a REM that predicts the SMC amplitudes, the SINR of SMCs, and the position error bound (PEB) at any position throughout the whole floor plan.
  • We analyze a database of UWB radio channel measurements to show to what extent the REM represents the multipath components in this set of measured channel impulse response (CIR) data. To this end, we study the number of resolvable SMCs and the energy capture (EC) of these SMCs, which is the fraction of energy carried by the SMCs in relation to the total energy of the CIR.
  • We show the usefulness of the REM in improving the robustness of a multipath-based single-anchor localization algorithm.
The structure of the paper is as follows. The UWB indoor signal model and geometry-based channel model are described in Section 2. The parameters represented by the REM, the SMC amplitudes, the SINRs, and the PEB, are also formulated. The experimental data, as well as data pre-processing steps, are described in Section 3. Section 4 discusses the Gaussian process model (GPM) based channel predictions for the above REM. Section 5 shows how the REM is exploited for a multipath-based single-anchor localization algorithm. Section 6 shows the localization results. Section 7 concludes the paper.

2. System and Environment Model

2.1. System Model

We consider an indoor wireless communication environment with a fixed physical anchor (PA) at position a 1 , as illustrated in Figure 1, and a mobile agent at position p . The agent is moved along a segmented trajectory. Some VAs a k up to order two are also shown, i.e., mirror images of the PA at the corresponding flat surfaces. The VAs establish a geometric model of SMCs originating from these surfaces. The VA-model can be used as a deterministic description of the delay, AOA, and AOD of these SMCs, given the position of the agent.
We describe the radio channel from the PA located at position a 1 and the agent node located at position p via a geometry-based stochastic channel model by:
h ( τ ; p ) = k = 1 K α k ( p ) δ ( τ τ k ( p ) ) + ν ( τ ; p )
The first term on the RHS of (1) is the sum of K SMCs with amplitude α k and delay τ k described as functions of the mobile agent position p . The delay τ k ( p ) corresponds to the distance between the agent and the PA (for k = 1 ) or VA (for k = 2 , , K ). That is τ k ( p ) = p a k / c , where a k is the position of the respective PA or VA, and c is the speed of light. In Figure 2, the path lengths d k ( p ) = p a k w.r.t various VAs at positions a k are plotted on top of the measured CIRs. The second term denotes the diffuse multipath (DM) ν ( τ ; p ) , modelled as a zero-mean Gaussian random process. We assume uncorrelated scattering along the delay axis, hence the autocorrelation function of ν ( τ ; p ) is given by:
K ν ( τ , u ; p ) = E ν ( τ ; p ) ν * ( u ; p ) = S ν ( τ ; p ) δ ( τ u ) ,
where S ν ( τ ; p ) is the power delay profile (PDP) of the DM at agent position p . The DM process is assumed to be quasi-stationary in the spatial domain, i.e., S ν ( τ ; p ) does not change in the vicinity of p [44].
A base band signal s ( t ) is transmitted at carrier frequency f c . The complex envelope of the received signal then read:
r ( t ; p ) = k = 1 K α k ( p ) s ( t τ k ( p ) ) + ( s ν ) ( t ; p ) + w ( t )
The third term of (3) is the measurement noise w ( t ) modelled as additive white Gaussian noise (AWGN) with double-sided power spectral density N 0 / 2 . To use (3) for estimating parameters of SMCs, it is useful to minimize autocorrelation side lobes of s ( t ) . A root-raised-cosine (RRC) pulse is thus considered for the transmitted signal with pulse duration T s . The energy of s ( t ) is normalized to one, i.e., | s ( t ) | 2 d t = 1 .

2.2. Position Error Bound

The channel model in (1) has been used in [27,38,47,48] to analyze multipath-assisted indoor positioning systems and to develop algorithms for it. In [27,48], the PEB has been derived for the signal model. It quantifies how the delay information of each SMC contributes useful information that can be exploited by a positioning algorithm. In this paper, the PEB is modeled as a spatial field for the current application environment, providing a measure for the achievable positioning accuracy in the given scenario. To develop this “environment model”, first, the PEB is explained in detail.
The PEB for the multipath-assisted positioning problem is the Cramer-Rao lower bound on the position error, the square-root of the trace of the upper left 2-by-2 submatrix of the inverse Fisher information matrix (FIM) J θ [48]
P p = tr J θ 1 2 × 2 = tr J p 1
where θ = p T R α T I α T T is a parameter vector containing the unknown location of the mobile agent, and the real and imaginary parts of the SMC complex amplitudes. J p is the equivalent FIM (EFIM) [26]. According to [48], we have
J p = 8 π 2 β 2 c 2 k = 1 K SINR k ( p ) J r ϕ k .
in which β 2 = f 2 | S ( f ) | 2 d f is the effective bandwidth ( S ( f ) is the Fourier transform of s ( t ) ) and c is the speed of light. The summation shows that the information from individual SMCs is added up. Each SMC contributes information proportional to the SINR k and the ranging direction matrix J r ϕ k [26]. It is obvious that the more SMC information is available, the lower the PEB. The SINR k ( p ) is defined as
SINR k ( p ) = | α k ( p ) | 2 N 0 + T p S ν ( τ k ; p )
with T p being the effective pulse duration [43]. It quantifies the interference of the DM, expressed by the PDP S ν ( τ k ; p ) , onto the k-the SMC with energy | α k ( p ) | 2 . The ranging direction matrix describes the geometry of the setup. The delay information of the k-th SMC is directed towards the k-th VA, expressed by
J r ( ϕ k ) = cos 2 ( ϕ k ) cos ϕ k sin ϕ k cos ϕ k sin ϕ k sin 2 ( ϕ k )
where ϕ k is the angle from the corresponding VA to the mobile agent at position p .
From (6) it can be seen that the range information contributed by the k-th SMC is quantified by SINR k ( p ) . Knowledge of SINR k ( p ) , therefore, allows to properly select and weigh the SMCs to be used by positioning algorithms. The idea of this paper is to model the SINR k ( p ) of the SMCs as a spatial field in order to obtain a map of the expected PEB for the current application environment. Such an “environment map” will provide the side information needed to make a positioning algorithm more robust.

2.3. Radio Environment Map Using Gaussian Process Regression

According to [34], a spatial field is the distribution of a physical quantity such as RF interference, pollution, temperature, humidity, light intensity, etc., which, if known, will foster the development of a new context-aware application and algorithms. In the context of radio systems, these spatial distributions are referred to as radio environment maps (REMs).
In this paper, we are interested in investigating the characteristics of individual SMCs, i.e., the direct link and the ones that resulted from reflections on flat surfaces. Spatial fields are obtained from mapping the agent positions with these SMC characteristics. In the UWB communication system, thanks to its very high time resolution, these SMCs can be resolved and predicted as studied in [43].
The SMC characteristics that are modeled in this paper are the SMC amplitudes and variances, thus the SINRs and the PEB can be predicted using (5) and (6). In Section 5, this information will be exploited to enhance the robustness of a multipath-assisted positioning algorithm. The SMC amplitudes and variances can also be used to predict the channel capacity and bit-error rate. The individual SMC metrics can then be useful for efficient beam selection and beam finding in mmWave wireless systems, exploiting prior information on potentially useful SMCs.
To construct the REM for the PEB, we need a spatial field model for the SMC amplitudes and the interference power of DM, c.f. (6). A Gaussian Process model (GPM) is used to learn such an amplitude model from observed data, as described in [43]. In a training phase, the SMC amplitudes α ^ k ( p ) are estimated from the CIR. Given the known positions of the VAs and the agent, the expected delays τ k ( p ) = 1 c a k p are deduced, yielding
α ^ k ( p ) = s ( t τ k ( p ) ) , r ( t ; p ) .
Next, the distance dependence is removed from the SMC amplitudes, yielding the normalized amplitude data γ k ( p ) = α ^ k ( p ) a k p exp j 2 π f c c a k p , which is assumed to depend on the direction angle ϕ k ( p ) = ( a k p ) , only, describing the antenna pattern and the reflectivity of the reflecting surface [43]. This direction dependence is modeled by means of GPR. Several dataset are obtained, i.e., D k = ϕ k ( p ) , γ k ( p ) p P , where P is the set of measured positions. Please note that whenever there are overlapping SMCs detected within a CIR, the corresponding amplitudes are neglected/discarded.
In this paper, the dataset D k abs = ϕ k ( p ) , | γ k ( p ) | p P is used, which deals with only the SMC amplitude’s abscissa. According to [43], the hyper-parameters of the GPM are θ k abs = a k abs , σ k abs , σ ν , k abs , β k abs T , in which a k abs is the coherence angle, σ k abs is the standard deviation of the correlation kernel, σ ν , k abs is the standard deviation accounting for the DM, and β k abs is the mean of the abscissa of the k-th SMC amplitude.
From these data sets and hyper-parameters θ k abs of the GPM, one can estimate the mean E [ ψ ( ϕ k ( p * ) ) | D k , θ k abs ] and variance V [ ψ ( ϕ k ( p * ) ) | D k , θ k abs ] of the normalized SMC amplitudes at any agent position of interest p * [43], and hence the SIN R k ( p * ) and the PEB P ( p * ) , which are the REMs investigated in this paper. Specifically, we get
SIN R k ( p * ) E 2 [ ψ ( ϕ k ( p * ) ) | D k , θ k abs ] 2 V [ ψ ( ϕ k ( p * ) ) | D k , θ k abs ] ,
c.f. Appendix B in [43]. Note that the predicted SIN R k ( p * ) is independent of the anchor-agent distance a k p * .
Section 4 shows how the spatial fields are learnt. Section 5 and Section 6 show how this information is utilized for the single-anchor multipath-based localization algorithm.

3. Experiment

3.1. Experiment Setup

An experiment was conducted in a laboratory room at the Graz University of Technology, as illustrated in Figure 1. The room consists of two plasterboard walls and two reinforced concrete walls (shown as black outer lines), three glass windows at the north wall (shown as thick gray lines), one whiteboard, and one wooden door at the south wall. In the following sections, we label these walls by east plasterboard (EPB), south wall (SW), west wall (WW), and north glass window (NGW), respectively.
A measurement data set was acquired with an Ilmsens Ultra-Wideband M-sequence channel sounder [49], c.f. Figure 3a, which provides a transmitted signal with 6.95 GHz carrier frequency to acquire a frequency band between approximately 3.5 and 10.5 GHz. The measurement principle is correlative channel sounding [50], in which a binary-code sequence with suitable correlation properties are transmitted and correlation with the known-code sequence is performed at the receiver to obtain the CIR.
Dipole coin antennas, shown in Figure 3b, have been used for the measurements. According to [51], the coin antenna has a very wide bandwidth ranging from 3 to 9 GHz and a nearly isotropic radiation pattern in the horizontal plane. As shown in Figure 1, a receive antenna was placed at a fixed position a 1 , whereas a transmit antenna (at position p ) was moved along a trajectory segmented into 7 parts for easier interpretation. Both 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. The raw measurements at the receiver port were filtered with an RRC pulse with center frequency 6.95 GHz, roll-off factor 0.5, and bandwidth 1 / T p = 2 GHz to obtain the received signals corresponding to the model in (3).

3.2. Experiment Pre-Processing

For the experiment, the exact positions of the mobile agent and the VAs were unknown. Thus, in a pre-processing step, we used the SLAM algorithm [38,45] and CIRs from the measurement to estimate the trajectory positions { p } and the locations of the VAs { a k } and PAs [43].
Different VAs can be distinguished from different traces of the path lengths, as illustrated in Figure 2.
Compared with [43], more VA positions are explored in this paper. The estimated trajectory and VAs are plotted in Figure 1. VAs 2, 3, 4, and 6 represent single reflections on the EPB, WW, SW, and NGW respectively. VAs 7, 8, 10, 11, 12, 13, 16, 17, and 19 represent double reflections on EPB then WW, EPB then SW, EPB then NGW, WW then EPB, WW then SW, SW then NGW, SW then WW, NGW then SW, and NGW then WW respectively.
Given the estimated positions of VAs { a k } , and the receiver positions { p } , the expected delays τ k ( p ) are computed. Next, the SMC amplitudes α ^ k ( p ) are estimated from the CIR using (8). Finally, the distance dependence will be removed from the SMC amplitudes, yielding the normalized amplitude data ψ ( ϕ k ( p ) ) , where ϕ k ( p ) = ( a k p ) is the direction angle. Several dataset D k = ϕ k ( p ) , ψ ( ϕ k ( p ) ) are then obtained.
On the basis of these data sets, the GPM is learned, using the built-in Matlab function fitrgp (Matlab version R2018a), yielding the hyper-parameters θ k abs = a k abs , σ k abs , σ ν , k abs , β k abs T . From these and the data, prediction of the SMC amplitudes is performed using the function predict. Other spatial field indicators, such as SINR, PEB are then deduced, c.f. (4) and (6). The resulting REM will be discussed in Section 4, and exploited in Section 5 and Section 6.

4. Analysis of the GPR-Based Channel Model

In this section, various REM parameters, i.e., SINR, PEB, etc will be analyzed, using the measurement data. We first analyze the energy “captured” by the SMCs in relation to the total energy of the observed CIRs, denoted as energy capture (EC) [52] and the number of visible SMCs along the trajectory. Secondly, the GPM of the SMC amplitudes and the resulting REMs for SMC amplitudes, SINR-values, and the PEB are illustrated.

4.1. Energy Capture

We analyze the EC of the SMCs, which is defined as the energy ratio between one individual or a few SMCs over the total energy of the CIR [52]. Figure 4 plots the EC of individual SMCs, while Figure 5 shows the combined ECs of different sets of SMCs. The LOS component captures between 20% to 80% of the total received signal power; the closer the agent to the anchor, the higher the EC. Individual first-order SMCs capture up to 20% of the total received power, while the accumulated EC of all first-order SMCs is between 5% to 35%. For second-order SMCs, the ECs are between 0.1% to 5%, except for VA 8 and 10 where the maximum EC is 25% and 20% respectively. Their accumulated EC is mostly between 2% to 20%.
Figure 5 also plots the number of visible VAs at each trajectory point. For first-order and especially second-order SMCs, the accumulated EC correlates with the number of visible VAs, as shown in Figure 5a,b. It could be because for second-order VAs, the energy are relatively the same, as shown in Figure 4b. For first-order VAs, the energies for different VAs are quite different.
Figure 5c plots the EC of the LOS signal versus the total number of visible VAs. It is shown that the closer the PA to the anchor, the higher the EC and the fewer the number of visible VAs, which are both understandable.
Figure 5d illustrates the EC of all SMCs versus the number of visible VAs. It is seen that, for our data set, between 50% and 90% of the received power can be modeled by means of SMCs. Correspondingly, between 50% and 90% of the received power can be predicted accurately by means of a GPR amplitude model. A beam selection and beam searching algorithm can select strong first and second-order SMCs as candidate beams, in case the LOS gets blocked (e.g., by the user itself).

4.2. SMC Amplitudes and SINR

Figure 6 shows the GPR for various normalized SMC amplitudes | α ^ k ( p ) d k ( p ) | as a function of the direction angle ϕ k ( p ) of the propagation path. The estimated amplitudes from the data set are shown along the mean and the ± 2 σ interval, as well as the predicted SINR k ( ϕ k ( p * ) ) given by (9). It is evident that there is a correlation between the SMC amplitudes and the direction angle, as described in Section 2.3, which is exploited by the GPR. The normalization does not allow for a direct comparison of the SMC amplitudes. But the comparison of the SINR-values indicates how much ranging information each of the SMCs can yield, c.f. (5).
It is shown that the SINR for reflections on the EPB is high, around 40, see Figure 6b which is much higher than the SINR for other SMCs, see Figure 6c–n, except the LOS link. Please note that on the east side of the room, there are fewer obstacles which makes the SMCs more reliable than others.
For the reflections on the WW, c.f. Figure 6c, the SINR is very low, indicating low ranging information, except at a direction approximately perpendicular to the wall, where the SINR increases steadily. The result is reasonable since, near the west wall, there are quite a few obstacles, e.g., computers, measurement equipment, etc., that could produce scattering effects.
For the SMC that corresponds to VA 4, shown in Figure 6d, the reflection on the whiteboard (visible at angles < 1.5 rad) has very high SINR, much higher than reflections on the other materials, e.g., the wooden door and the plasterboard wall. From VA 6, c.f. Figure 6e, it is shown that reflections on the NGW with metal-coated glass windows yield very significant ranging information, reaching up to an SINR of around 20.
As shown in Figure 6f,g,i,k, for VA 7, 8, 11 and 13, the SINR is not as good as the first-order SMCs, but better than the other second-order SMCs shown in Figure 6h,j,l,m,n.
Implementing GPR on measurement data allows us to predict the SINR as well as the amplitude of the SMCs for each point p * in an environment at which an SMC is visible. Figure 7 shows this prediction result in terms of the SINR k ( p * ) , while Figure 8 shows the predicted mean amplitude E [ ψ ( ϕ k ( p * ) ) | D k , θ k abs ] / a k p * . It is evident that the LOS path has very high SINR compared to other SMCs, thus it yields very high range information. However, in the case that the LOS is not available, choosing a good NLOS SMC is necessary. Two parameters could be used to choose between SMCs, the SINR and the SMC amplitude. The SINR is an appropriate measure to describe the ranging information provided by some SMC, which is an important metric for (multipath-assisted) positioning algorithms (see Section 5), while the amplitude translates to communication-performance metrics like channel capacity [53] and bit error rate.
From Figure 8, it is seen that the LOS amplitude is dependent on the distance as well as the antenna radiating pattern, while for other SMCs, the amplitudes are also affected by the reflection coefficients of the plane surfaces.

4.3. Position Error Bound

Figure 9 shows the predicted PEB throughout the floor plan, computed from (5), when using information from various SMCs. The information provided by each SMC is quantified by SINR k ( p * ) . The sub-figures on the left-hand side include the LOS component from the PA, while the LOS is neglected in all sub-figures on the right. In the top row, all VAs are considered, while first-order VAs only are considered in the middle, and at the bottom, only VA 2, VA 4, VA 6, VA 7, and VA 8 are used, which contribute most to the PEB.
It is observed that the second-order VAs play an important role in any localization algorithm, despite their low powers. With information from them, the PEB improves significantly even without the LOS component. With only first-order VAs, the PEB is worse.

5. Exploiting the GPR Model for (Multipath-Assisted) Positioning

To formulate a signal model for position estimation, Equation (3) can be rewritten as
r = S ( p ) α + w CN S ( p ) α , C
where r = , r ( n T s ) , C N is a sampled version of the received signal, S ( p ) = , s ( τ k ( p ) ) , is the parameterized signal matrix with s ( τ ) = , s ( n T s τ ) , T C N , denoting the SMCs with delays τ k ( p ) that are determined by p , α = , α k , T C K are the SMC amplitudes, and w is a noise vector with covariance C = E w w H . The noise vector w may account for AWGN only or for the DM as well. A maximum likelihood (ML) estimation p ^ = arg max p * L r | p * considers α as nuisance parameters and finds the position by maximizing the concentrated log-likelihood function (LLHF)
L r | p * = max α L r | p * , α .
The LLHF evaluated at candidate position p * follows from (10),
L r | p * , α = log det ( C ) r S ( p ) α H C 1 r S ( p ) α .
The maximization in (11) is solved by a weighted least square solution for r S ( p ) α 2 which yields α ^ = S H ( p ) C 1 S ( p ) 1 S H ( p ) C 1 r . This ML positioning algorithm is our reference algorithm, denoted as “SALMA-light” and described in [54,55] if AWGN is assumed, i.e., C = σ n 2 I .
The GPM model provides prior knowledge about the absolute value of α but not on its phase. We thus expand the amplitude vector as α = Φ x , where diagonal matrix Φ = diag ( ϕ ) and [ ϕ ] k = exp ( j α k ) = exp ( j ζ k ) . From the GPM, we know that x | p N μ ( p ) , 1 2 Λ ( p ) with [ μ ( p ) ] k = E ψ ( ϕ k ( p ) ) | D k , θ k abs a k p , and the diagonal matrix has elements [ Λ ( p ) ] k , k = 2 V ψ ( ϕ k ( p ) ) | D k , θ k abs a k p 2 . The signal model then becomes
r = S ( p ) Φ μ ( p ) + S ( p ) x ˜ + w
with x ˜ | p CN 0 , Λ ( p ) is introduced to account for the variance of the SMC amplitudes and the uncertainty of the phases in Φ .
The concentrated LLHF is now obtained by maximizing the LLHF for signal model (13) by maximizing for the nuisance parameter vector ϕ and σ n 2 (maximizing the AWGN model for C ), given candidate position p * ,
F p * = max σ n , ϕ L r | σ n 2 , ϕ , p *
with
L r | σ n , ϕ , p * = log det ( C ˜ ) r S ( p ) M ( p ) ϕ H C ˜ 1 r S ( p ) M ( p ) ϕ
where M ( p ) = diag ( μ ( p ) ) and
C ˜ = σ n 2 I + S ( p ) Λ ( p ) S ( p ) H
for the noise covariance.
First, we will find the phase ϕ , and noise variance σ n 2 that maximize (15), i.e.,
ϕ ^ = arg max ϕ L r | σ n , ϕ , p *
σ ^ n = arg max σ n L r | σ n 2 , ϕ ^ , p *
Let’s denote the phase of the k-th SMC amplitude as ζ k . To solve for the phase, we have:
r S ( p ) M ( p ) ϕ H C ˜ 1 r S ( p ) M ( p ) ϕ = r k k s ( τ k ) Λ k e j ζ k s ( τ k ) μ k e j ζ k H C ˜ 1 r k k s ( τ k ) μ k e j ζ k s ( τ k ) μ k e j ζ k = r k s ( τ k ) μ k e j ζ k H C ˜ 1 r k s ( τ k ) μ k e j ζ k = f ( ζ k )
Vector r k is the residual of the CIR measurement after subtracting all SMCs except the k-th.
In order to find ζ k , let’s denote e j ζ k = x k + j y k . Thus:
x ^ k , y ^ k = arg min x k , y k f ( x k , y k ) s . t . g ( x k , y k ) = x k 2 + y k 2 1 = 0
Using Lagrange multiplier, we have
x ^ k , y ^ k , λ ^ k = arg min x k , y k , λ k f ( x k , y k ) + λ k g ( x k , y k )
and obtain
L x k = 2 μ k s ( τ k ) H C ˜ 1 Re r k + 2 x k μ k 2 s ( τ k ) H C ˜ 1 s ( τ k ) + 2 x k λ k
L y k = 2 μ k s ( τ k ) H C ˜ 1 Im r k + 2 y k μ k 2 s ( τ k ) H C ˜ 1 s ( τ k ) + 2 y k λ k
L λ k = x k 2 + y k 2 1
Equating Equations (22)–(24) to 0, we have:
λ k = | s H ( τ k ) C ˜ 1 μ k r k | μ k 2 s ( τ k ) H C ˜ 1 s ( τ k )
and
x k = s ( τ k ) H C ˜ 1 μ k Re { r k } μ k 2 s ( τ k ) H C ˜ 1 s ( τ k ) + λ k
Similarly for y k . Or, in short we have
ζ k = s ( τ k ) H C ˜ 1 μ k r k
and if the k-th SMC does not overlap with any other SMC, we have equivalently
ζ ^ k = s ( τ k ) C ˜ 1 r
which is a weighted projection of the measurement on the k-th SMC template s ( τ k ) . The weighting (by C ˜ 1 ) takes into account the variance of the SMC from the GPR model, which is essential for reducing secondary maxima in the LHF. Unfortunately, the weighting has a negative influence on the reliability of this phase estimation step. We thus replace for this estimator the covariance by an identity matrix, obtaining the approximate (but more robust) phase estimator
ζ ^ k s ( τ k ) r .
From (15), we can write:
L r | σ n , ϕ , p * = log det C ˜ trace C ˜ 1 C ˜ ^
where C ˜ ^ = r S ( p ) M ( p ) ϕ r S ( p ) M ( p ) ϕ H is a realization of C ˜ . According to [56], we have:
σ n ^ 2 = 1 N K trace Π S ( τ ) C ˜ ^
in which
Π S ( τ ) = I N S ( τ ) S ( τ ) H S ( τ ) 1 S ( τ ) H .
Thus σ ^ n can be computed. The localization problem reduces to
p ^ = arg max p * L r | σ ^ n , ϕ ^ , p * .

6. Result

To implement the GPR-based positioning algorithm, the data described in Section 3 is divided into two groups. The first group is used for GP training, in which both the agent positions and the corresponding CIRs are known. The GPM hyper-parameters are estimated. The second group contains the remaining data, where only the CIRs are used. The corresponding agent positions are unknown and will be estimated by exploiting the REM, i.e., the SMC amplitudes and SINRs predicted using the GPM hyper-parameters. Errors are computed and their CDFs are plotted in Figure 10. SALMA-light [54,55], the reference single-anchor localization method described in Section 5, is used for comparison.
It is shown that, with situational awareness through the REM information, the proposed algorithm outperforms SALMA-light in terms of the outlier robustness, while the accuracy remains comparable.
To explain this performance improvement, Figure 11 compares the LLHF for the two algorithms. The LLHF has been evaluated at candidate positions lying on a radius that corresponds to the estimated line-of-sight distance [55]. It is shown that the REM-assisted algorithm reduces the number of local maxima in the LLHF and emphasizes the correct maximum, yielding the performance improvement over SALMA-light.
Figure 12 shows the error CDF for a modified experiment, where a visibility test and overlap checking are conducted to reduce the number of SMCs that are used for position estimation to the most useful ones. This is done for both, SALMA-light and the proposed algorithm. It is shown that again the proposed algorithm outperforms SALMA-light in reducing outliers. The difference in results between Figure 10 and Figure 12 is because, after visibility test and overlap checking, the number of VAs used is reduced, leading to reduced robustness for SALMA light while the use of the REM avoids such loss of robustness. Note that also the accuracy is improved in this variant.
These results suggest that the proposed algorithm might perform better in obstructed-LOS or in the presence of humans or obstacles, even though the algorithm has not been tailored to specifically address this problem. The recent publications [57,58] describe algorithms that are designed to be robust to blocked LOS (in the former) and blocked SMCs (in the latter), in particular when the blocking is due to the influence of the human body on the visibility of SMCs. The belief propagation SLAM algorithm from [38] provides visibility information of SMCs in addition to the amplitude and delay information. This visibility information can be used to pre-select the data points for the construction of the GPR-based REM such that path blocking does not influence the quality of the REM. It is expected that the “environment part” of the propagation channel, which is expressed by the GPR-based REM, and the “device part” can be separated very well, once a consistent REM has been learned from a large number of independent measurements obtained from different user equipment. The rigorous formulation of such methods is left for future work, however.

7. Conclusions

In this paper, radio environment maps (REMs) have been formulated, taking into account amplitude information of specular multipath components (SMCs) that are learned by means of Gaussian process regression (GPR). The regression method allows predicting the mean and variance of the SMC amplitudes. From these parameters, the position-error-bound (PEB) has been computed as a performance metric for positioning algorithms that exploit range information contained in the SMCs, so-called multipath-based localization schemes.
The amplitude information can also be used to formulate REMs of communication performance metrics, e.g., the channel capacity or the bit-error rate, in order to realize more efficient and robust, location-aware protocols for medium access control. Using a set of training measurements from the environment, the REMs can be computed for any point within the environment. The paper evaluates to what extent the REM represents the multipath components in the analyzed data. It has been found that 50% to 90% of the received power can be attributed to the resolvable SMCs contained in our model, which corresponds to the fraction of power that can be predicted by the REM.
The REM is applied in this paper to improve the robustness of a single-anchor multipath-based localization method which sometimes suffers from outliers resulting in large errors. Using the REM as a channel prior knowledge helps to avoid these outliers.

Author Contributions

Conceptualization, H.A.N. and K.W.; methodology, H.A.N. and K.W.; software, H.A.N.; validation, H.A.N.; formal analysis, H.A.N.; investigation, H.A.N.; resources, K.W.; writing—original draft preparation, H.A.N.; writing—review and editing, H.A.N. and K.W.; visualization, H.A.N.; supervision, V.K.N. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

The research of K.W. was partly funded by TU Graz under the Lead Project “Dependable Things”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The measurement data are available on request from the authors, c.f. https://www.spsc.tugraz.at/databases-and-tools/uwb-indoor-channel-experimental-data.html.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, L.D.; He, W.; Li, S. Internet of Things in Industries: A Survey. IEEE Trans. Ind. Inform. 2014, 10, 2233–2243. [Google Scholar] [CrossRef]
  2. Acampora, G.; Cook, D.; Rashidi, P.; Vasilakos, A. A Survey on Ambient Intelligence in Healthcare. Proc. IEEE 2013, 101, 2470–2494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Borgia, E. The Internet of Things vision: Key features, applications and open issues. Comput. Commun. 2014, 54, 1–31. [Google Scholar] [CrossRef]
  4. Andrews, J.G.; Buzzi, S.; Choi, W.; Hanly, S.V.; Lozano, A.; Soong, A.C.K.; Zhang, J.C. What Will 5G Be? IEEE J. Sel. Areas Commun. 2014, 32, 1065–1082. [Google Scholar] [CrossRef]
  5. 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]
  6. Jungnickel, V.; Manolakis, K.; Zirwas, W.; Panzner, B.; Braun, V.; Lossow, M.; Sternad, M.; Apelfrojd, R.; Svensson, T. The role of small cells, coordinated multipoint, and massive MIMO in 5G. IEEE Commun. Mag. 2014, 52, 44–51. [Google Scholar] [CrossRef]
  7. del Peral-Rosado, J.A.; Raulefs, R.; López-Salcedo, J.A.; Seco-Granados, G. Survey of Cellular Mobile Radio Localization Methods: From 1G to 5G. IEEE Commun. Surv. Tutor. 2017, 20, 1124–1148. [Google Scholar] [CrossRef]
  8. Perera, C.; Zaslavsky, A.; Christen, P.; Georgakopoulos, D. Context Aware Computing for The Internet of Things: A Survey. IEEE Commun. Surv. Tutor. 2014, 16, 414–454. [Google Scholar] [CrossRef] [Green Version]
  9. Wymeersch, H.; Seco-Granados, G.; Destino, G.; Dardari, D.; Tufvesson, F. 5G mm-Wave Positioning for Vehicular Networks. IEEE Wireless Commun. Mag. 2017, 24, 80–86. [Google Scholar] [CrossRef] [Green Version]
  10. Hult, R.; Campos, G.R.; Steinmetz, E.; Hammarstrand, L.; Falcone, P.; Wymeersch, H. Coordination of Cooperative Autonomous Vehicles: Toward safer and more efficient road transportation. IEEE Signal Process. Mag. 2016, 33, 74–84. [Google Scholar] [CrossRef]
  11. Karlsson, R.; Gustafsson, F. The Future of Automotive Localization Algorithms: Available, reliable, and scalable localization: Anywhere and anytime. IEEE Signal Process. Mag. 2017, 34, 60–69. [Google Scholar] [CrossRef] [Green Version]
  12. Koivisto, M.; Hakkarainen, A.; Costa, M.; Kela, P.; Leppänen, K.; Valkama, M. High-Efficiency Device Positioning and Location-Aware Communications in Dense 5G Networks. IEEE Commun. Mag. 2017, 55, 188–195. [Google Scholar] [CrossRef] [Green Version]
  13. Chen, L.; Thombre, S.; Jarvinen, K.; Lohan, E.S.; Alen-Savikko, A.K.; Leppakoski, H.; Bhuiyan, M.Z.H.; Bu-Pasha, S.; Ferrara, G.N.; Honkala, S.; et al. Robustness, Security and Privacy in Location-Based Services for Future IoT: A Survey. IEEE Access 2017, 5, 2169–3536. [Google Scholar] [CrossRef]
  14. 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]
  15. Mazuelas, S.; Bahillo, A.; Lorenzo, R.; Fernandez, P.; Lago, F.; Garcia, E.; Blas, J.; Abril, E. Robust Indoor Positioning Provided by Real-Time RSSI Values in Unmodified WLAN Networks. IEEE J. Sel. Top. Signal Process. 2009, 3, 821–831. [Google Scholar] [CrossRef]
  16. Fadzilla, M.A.; Harun, A.; Shahriman, A.B. Localization Assessment for Asset Tracking Deployment by Comparing an Indoor Localization System with a Possible Outdoor Localization System. In Proceedings of the 2018 International Conference on Computational Approach in Smart Systems Design and Applications (ICASSDA), Kuching, Malaysia, 15–17 August 2018; pp. 1–6. [Google Scholar] [CrossRef]
  17. Raspopoulos, M.; Laoudias, C.; Kanaris, L.; Kokkinis, A.; Panayiotou, C.G.; Stavrou, S. 3D Ray Tracing for device-independent fingerprint-based positioning in WLANs. In Proceedings of the 2012 9th Workshop on Positioning, Navigation and Communication, Dresden, Germany, 15–16 March 2012; pp. 109–113. [Google Scholar] [CrossRef]
  18. Hatami, A.; Pahlavan, K. Comparative statistical analysis of indoor positioning using empirical data and indoor radio channel models. In Proceedings of the CCNC 2006, 2006 3rd IEEE Consumer Communications and Networking Conference, Las Vegas, NV, USA, 8–10 January 2006; Volume 2, pp. 1018–1022. [Google Scholar] [CrossRef]
  19. Chai, X.; Yang, Q. Reducing the Calibration Effort for Probabilistic Indoor Location Estimation. IEEE Trans. Mob. Comput. 2007, 6, 649–662. [Google Scholar] [CrossRef] [Green Version]
  20. Youssef, M.; Agrawala, A. The Horus WLAN Location Determination System. In Proceedings of the 3rd International Conference on Mobile Systems, Applications, and Services, MobiSys ’05, Seattle, WA, USA, 6–8 June 2005; ACM: New York, NY, USA, 2005; pp. 205–218. [Google Scholar] [CrossRef]
  21. Cassarà, P.; Potortì, F.; Barsocchi, P.; Girolami, M. Choosing an RSS device-free localization algorithm for Ambient Assisted Living. In Proceedings of the 2015 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Banff, AB, Canada, 13–16 October 2015; pp. 1–8. [Google Scholar] [CrossRef]
  22. 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]
  23. Piccinni, G.; Avitabile, G.; Coviello, G.; Talarico, C. Real-Time Distance Evaluation System for Wireless Localization. IEEE Trans. Circuits Syst. I Regul. Pap. 2020, 67, 3320–3330. [Google Scholar] [CrossRef]
  24. van Trees, H.L. Detection, Estimation, and Modulation Theory, Part I: Detection, Estimation, and Linear Modulation Theory; John Wiley & Sons: Hoboken, NJ, USA, 2001. [Google Scholar]
  25. Aditya, S.; Molisch, A.F.; Behairy, H.M. A Survey on the Impact of Multipath on Wideband Time-of-Arrival-Based Localization. Proc. IEEE 2018, 106, 1183–1203. [Google Scholar] [CrossRef]
  26. Shen, Y.; Win, M. Fundamental Limits of Wideband Localization—Part I: A General Framework. IEEE Trans. Inform. Theory 2010, 56, 4956–4980. [Google Scholar] [CrossRef]
  27. Witrisal, K.; Meissner, P. Performance Bounds for Multipath-assisted Indoor Navigation and Tracking (MINT). In Proceedings of the 2012 IEEE International Conference on Communications (ICC), Ottawa, ON, Canada, 10–15 June 2012. [Google Scholar]
  28. Witrisal, K.; Leitinger, E.; Hinteregger, S.; Meissner, P. Bandwidth Scaling and Diversity Gain for Ranging and Positioning in Dense Multipath Channels. IEEE Wirel. Commun. Lett. 2016, 5, 396–399. [Google Scholar] [CrossRef]
  29. Han, Y.; Shen, Y.; Zhang, X.P.; Win, M.Z.; Meng, H. Performance limits and geometric properties of array localization. IEEE Trans. Inform. Theory 2016, 62, 1054–1075. [Google Scholar] [CrossRef] [Green Version]
  30. 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]
  31. Marmasse, N.; Schmandt, C. Location-Aware Information Delivery withComMotion. In Handheld and Ubiquitous Computing; Thomas, P., Gellersen, H.W., Eds.; Springer: Berlin/Heidelberg, Germany, 2000; pp. 157–171. [Google Scholar]
  32. Cadger, F.; Curran, K.; Santos, J.; Moffett, S. A Survey of Geographical Routing in Wireless Ad-Hoc Networks. IEEE Commun. Surv. Tutor. 2013, 15, 621–653. [Google Scholar] [CrossRef]
  33. Yilmaz, H.B.; Tugcu, T.; Alagöz, F.; Bayhan, S. Radio environment map as enabler for practical cognitive radio networks. IEEE Commun. Mag. 2013, 51, 162–169. [Google Scholar] [CrossRef]
  34. Arias-de Reyna, E.; Closas, P.; Dardari, D.; Djuric, P.M. Crowd-Based Learning of Spatial Fields for the Internet of Things: From Harvesting of Data to Inference. IEEE Signal Process. Mag. 2018, 35, 130–139. [Google Scholar] [CrossRef]
  35. Muppirisetty, L.S.; Svensson, T.; Wymeersch, H. Spatial Wireless Channel Prediction under Location Uncertainty. IEEE Trans. Wirel. Commun. 2016, 15, 1031–1044. [Google Scholar] [CrossRef] [Green Version]
  36. 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]
  37. Wang, Y.; Narasimha, M.; Heath, R.W. MmWave Beam Prediction with Situational Awareness: A Machine Learning Approach. In Proceedings of the 2018 IEEE 19th International Workshop on Signal Processing Advances in Wireless Communications (SPAWC), Kalamata, Greece, 25–28 June 2018; pp. 1–5. [Google Scholar] [CrossRef] [Green Version]
  38. 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]
  39. 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]
  40. Meissner, P.; Leitinger, E.; Lafer, M.; Witrisal, K. Real-time demonstration of multipath-assisted indoor navigation and tracking (MINT). In Proceedings of the 2014 IEEE International Conference on Communications Workshops (ICC), Sydney, Australia, 10–14 June 2014; pp. 144–149. [Google Scholar] [CrossRef]
  41. Meissner, P. Multipath-Assisted Indoor Positioning. Ph.D. Thesis, Graz University of Technology, Graz, Austria, 2014. [Google Scholar]
  42. Ulmschneider, M. Cooperative Multipath Assisted Positioning. Ph.D. Thesis, Hamburg University of Technology, Hamburg, Germany, 2021. [Google Scholar] [CrossRef]
  43. Nguyen, A.H.; Rath, M.; Leitinger, E.; Nguyen, K.V.; Witrisal, K. Gaussian Process Modeling of Specular Multipath Components. Appl. Sci. 2020, 10, 5216. [Google Scholar] [CrossRef]
  44. Molisch, A.F. Ultra-Wide-Band Propagation Channels. Proc. IEEE 2009, 97, 353–371. [Google Scholar] [CrossRef]
  45. 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]
  46. Rath, M. Signal Process for Localization and Environment Mapping. Ph.D. Thesis, Graz University of Technology, Graz, Austria, 2021. [Google Scholar]
  47. Meissner, P.; Leitinger, E.; Witrisal, K. UWB for Robust Indoor Tracking: Weighting of Multipath Components for Efficient Estimation. IEEE Wirel. Commun. Lett. 2014, 3, 501–504. [Google Scholar] [CrossRef]
  48. Leitinger, E.; Meissner, P.; Rüdisser, 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]
  49. 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]
  50. Molisch, A.F. Wireless Communications, 2nd ed.; Wiley-IEEE Press: Hoboken, NJ, USA, 2010. [Google Scholar]
  51. Krall, C. Signal Processing for Ultra Wideband Transceivers. Ph.D. Thesis, Graz University of Technology, Graz, Austria, 2008. [Google Scholar]
  52. Win, M.; Scholtz, R. On the robustness of ultra-wide bandwidth signals in dense multipath environments. IEEE Commun. Lett. 1998, 2, 51–53. [Google Scholar] [CrossRef]
  53. 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]
  54. Rath, M.; Kulmer, J.; Leitinger, E.; Witrisal, K. Single-Anchor Positioning: Multipath Processing with Non-Coherent Directional Measurements. IEEE Access 2020, 8, 88115–88132. [Google Scholar] [CrossRef]
  55. Großwindhager, B.; Rath, M.; Kulmer, J.; Bakr, M.S.; Boano, C.A.; Witrisal, K.; Römer, K. SALMA: UWB-based Single-Anchor Localization System using Multipath Assistance. In Proceedings of the 16th ACM Conference on Embedded Networked Sensor Systems, Shenzhen, China, 4–7 November 2018. [Google Scholar]
  56. Stoica, P.; Nehorai, A. On the concentrated stochastic likelihood function in array signal processing. Circuits Syst. Signal Process. 1995, 14, 669–674. [Google Scholar] [CrossRef]
  57. Venus, A.; Leitinger, E.; Tertinek, S.; Witrisal, K. A Message Passing based Adaptive PDA Algorithm for Robust Radio-based Localization and Tracking. In Proceedings of the 2021 IEEE Radar Conference (RadarConf21), Atlanta, GA, USA, 7–14 May 2021. [Google Scholar] [CrossRef]
  58. Wilding, T.; Leitinger, E.; Witrisal, K. Multipath-based Localization and Tracking considering Off-Body Channel Effects. arXiv 2021, arXiv:2110.09932. [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. One blue cross represents the PA; orange circles denote VAs which were considered in the experimental evaluation. An agent moves along a trajectory segmented into seven parts indicated with distinct colors.
Figure 1. Floor plan of the evaluation scenario. Bold black lines denote walls, thick gray lines represent glass windows, other lines illustrate other materials. One blue cross represents the PA; orange circles denote VAs which were considered in the experimental evaluation. An agent moves along a trajectory segmented into seven parts indicated with distinct colors.
Sensors 22 00462 g001
Figure 2. Measured CIRs over the trajectory at a bandwidth of 2 GHz. The blue dash lines indicate expected paths of the SMCs corresponding to VAs up to order two, estimated using the SLAM algorithm from [38,45] and refined as described in [46]. The locations of VAs are shown in Figure 1.
Figure 2. Measured CIRs over the trajectory at a bandwidth of 2 GHz. The blue dash lines indicate expected paths of the SMCs corresponding to VAs up to order two, estimated using the SLAM algorithm from [38,45] and refined as described in [46]. The locations of VAs are shown in Figure 1.
Sensors 22 00462 g002
Figure 3. Equipments used for the experiment: (a) Ilmsens UWB M-sequence sounder and Ilmsens power supply, (b) coin antenna.
Figure 3. Equipments used for the experiment: (a) Ilmsens UWB M-sequence sounder and Ilmsens power supply, (b) coin antenna.
Sensors 22 00462 g003
Figure 4. EC of individual SMC: (a) VA 1 to VA 6, (b) VA 7 to VA 19.
Figure 4. EC of individual SMC: (a) VA 1 to VA 6, (b) VA 7 to VA 19.
Sensors 22 00462 g004
Figure 5. EC of certain SMCs vs. number of VAs visible: (a) first-order VAs only, (b) second-order VAs only, (c) LOS (PA) only, (d) all VAs and LOS (PA). The brown solid line shows the number of VAs that are visible from the corresponding trajectory point. The dashed brown line shows the number of visible VAs that are within the set of (a) first-order VAs or (b) second-order VAs.
Figure 5. EC of certain SMCs vs. number of VAs visible: (a) first-order VAs only, (b) second-order VAs only, (c) LOS (PA) only, (d) all VAs and LOS (PA). The brown solid line shows the number of VAs that are visible from the corresponding trajectory point. The dashed brown line shows the number of visible VAs that are within the set of (a) first-order VAs or (b) second-order VAs.
Sensors 22 00462 g005
Figure 6. GPR on SMC amplitudes corresponding to: (a) LOS link, (b) reflection on EPB, (c) reflection on WW, (d) reflection on the south part of the room containing various materials, e.g., white board, wall, metal door, (e) reflection on the north part of the room, containing glass windows and walls, (f) double reflection on EPB then WW, (g) double reflection on EPB then SW, (h) double reflection on EPB then NGW, (i) double reflection on WW then EPB, (j) double reflection on WW then SW, (k) double reflection on SW then NGW, (k) double reflection on SW then WW, (l) double reflection on SW then WW, (m) double reflection on NGW then SW, and (n) double reflection on NGW then WW. The regressed amplitude is | α k ( p ) d k ( p ) | . SINR computed from GPR is plotted in purple with axis on the right-hand side.
Figure 6. GPR on SMC amplitudes corresponding to: (a) LOS link, (b) reflection on EPB, (c) reflection on WW, (d) reflection on the south part of the room containing various materials, e.g., white board, wall, metal door, (e) reflection on the north part of the room, containing glass windows and walls, (f) double reflection on EPB then WW, (g) double reflection on EPB then SW, (h) double reflection on EPB then NGW, (i) double reflection on WW then EPB, (j) double reflection on WW then SW, (k) double reflection on SW then NGW, (k) double reflection on SW then WW, (l) double reflection on SW then WW, (m) double reflection on NGW then SW, and (n) double reflection on NGW then WW. The regressed amplitude is | α k ( p ) d k ( p ) | . SINR computed from GPR is plotted in purple with axis on the right-hand side.
Sensors 22 00462 g006aSensors 22 00462 g006b
Figure 7. Predicted SINR for the whole floor plan for various links: (a) LOS, (b) reflection on EPB, (c) reflection on WW, (d) reflection on on south part of the room, (e) reflection on the north part of the room, (f) double reflection on EPB then WW, (g) double reflection on WW then SW, (h) double reflection on EPB then NGW, (i) double reflection on WW then EPB, (j) double reflection on WW then SW, (k) double reflection on SW then GW, (l) double reflection on SW then WW, (m) double reflection on NGW then SW, (n) double reflection on NGW then WW. SINR is obtained by applying GPR on the CIRs obtained from measurement at 595 points along the trajectory, as shown in Figure 1. The same dB scale is used in all sub figures.
Figure 7. Predicted SINR for the whole floor plan for various links: (a) LOS, (b) reflection on EPB, (c) reflection on WW, (d) reflection on on south part of the room, (e) reflection on the north part of the room, (f) double reflection on EPB then WW, (g) double reflection on WW then SW, (h) double reflection on EPB then NGW, (i) double reflection on WW then EPB, (j) double reflection on WW then SW, (k) double reflection on SW then GW, (l) double reflection on SW then WW, (m) double reflection on NGW then SW, (n) double reflection on NGW then WW. SINR is obtained by applying GPR on the CIRs obtained from measurement at 595 points along the trajectory, as shown in Figure 1. The same dB scale is used in all sub figures.
Sensors 22 00462 g007aSensors 22 00462 g007b
Figure 8. Predicted SMC amplitude α k ( p ) for the whole floor plan for various SMCs: (a) LOS, (b) reflection on EPB, (c) reflection on WW, (d) reflection on on south part of the room, (e) reflection on the north part of the room, (f) double reflection on EPB then WW, (g) double reflection on WW then EPB, (h) double reflection on the SW then NGW. GPR was applied the set of measurement along the trajectory. The same dB scale is used in all sub-figures.
Figure 8. Predicted SMC amplitude α k ( p ) for the whole floor plan for various SMCs: (a) LOS, (b) reflection on EPB, (c) reflection on WW, (d) reflection on on south part of the room, (e) reflection on the north part of the room, (f) double reflection on EPB then WW, (g) double reflection on WW then EPB, (h) double reflection on the SW then NGW. GPR was applied the set of measurement along the trajectory. The same dB scale is used in all sub-figures.
Sensors 22 00462 g008
Figure 9. GPR-predicted PEB considering information from SMCs corresponding to: (a) all VAs and the PA, (b) all 1st and 2nd-order VAs, (c) 1st-order VAs and the PA, (d) all 1st-order VAs, (e) selected VAs and the LOS, (f) selected VAs.
Figure 9. GPR-predicted PEB considering information from SMCs corresponding to: (a) all VAs and the PA, (b) all 1st and 2nd-order VAs, (c) 1st-order VAs and the PA, (d) all 1st-order VAs, (e) selected VAs and the LOS, (f) selected VAs.
Sensors 22 00462 g009
Figure 10. CDF of the localization error using (a) training segments 1, 2, 3, 7, (b) training segments 4, 5, 6. The blue solid line plots the error CDF obtained when SALMA-light is used. The orange solid line plots the error CDF obtained when optimizing the LLHF with prior knowledge from GPR. The phase ζ k is approximated by using the phase from CIR. Overlap checking is performed first before implementing the algorithms.
Figure 10. CDF of the localization error using (a) training segments 1, 2, 3, 7, (b) training segments 4, 5, 6. The blue solid line plots the error CDF obtained when SALMA-light is used. The orange solid line plots the error CDF obtained when optimizing the LLHF with prior knowledge from GPR. The phase ζ k is approximated by using the phase from CIR. Overlap checking is performed first before implementing the algorithms.
Sensors 22 00462 g010
Figure 11. LLHF values obtained using (a) SALMA-light, (b) LHF+GPR with exact sigma. LLHF values at evaluation points are plotted with color code. The red circle denotes the estimation result, while the filled black circle denotes the true position.
Figure 11. LLHF values obtained using (a) SALMA-light, (b) LHF+GPR with exact sigma. LLHF values at evaluation points are plotted with color code. The red circle denotes the estimation result, while the filled black circle denotes the true position.
Sensors 22 00462 g011
Figure 12. CDF of the localization error using (a) training segments 1, 2, 3, 7, (b) training segments 4, 5, 6. The blue solid line plots the error CDF obtained when SALMA-light is used. The orange solid line plots the error CDF obtained when optimizing the LLHF with prior knowledge from GPR. The phase ζ k is approximated by using the phase from CIR. Visibility test is carried out, followed by overlap checking, before implementing the algorithms.
Figure 12. CDF of the localization error using (a) training segments 1, 2, 3, 7, (b) training segments 4, 5, 6. The blue solid line plots the error CDF obtained when SALMA-light is used. The orange solid line plots the error CDF obtained when optimizing the LLHF with prior knowledge from GPR. The phase ζ k is approximated by using the phase from CIR. Visibility test is carried out, followed by overlap checking, before implementing the algorithms.
Sensors 22 00462 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nguyen, H.A.; Nguyen, V.K.; Witrisal, K. Amplitude Modeling of Specular Multipath Components for Robust Indoor Localization. Sensors 2022, 22, 462. https://doi.org/10.3390/s22020462

AMA Style

Nguyen HA, Nguyen VK, Witrisal K. Amplitude Modeling of Specular Multipath Components for Robust Indoor Localization. Sensors. 2022; 22(2):462. https://doi.org/10.3390/s22020462

Chicago/Turabian Style

Nguyen, Hong Anh, Van Khang Nguyen, and Klaus Witrisal. 2022. "Amplitude Modeling of Specular Multipath Components for Robust Indoor Localization" Sensors 22, no. 2: 462. https://doi.org/10.3390/s22020462

APA Style

Nguyen, H. A., Nguyen, V. K., & Witrisal, K. (2022). Amplitude Modeling of Specular Multipath Components for Robust Indoor Localization. Sensors, 22(2), 462. https://doi.org/10.3390/s22020462

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