Next Article in Journal
Field Demonstration of a Distributed Microsensor Network for Chemical Detection
Next Article in Special Issue
Integration of Microstrip Slot Array Antenna with Dye-Sensitized Solar Cells
Previous Article in Journal
MODIS Sensor Capability to Burned Area Mapping—Assessment of Performance and Improvements Provided by the Latest Standard Products in Boreal Regions
Previous Article in Special Issue
Small Antennas for Wearable Sensor Networks: Impact of the Electromagnetic Properties of the Textiles on Antenna Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Letter

Timestamp Estimation in P802.15.4z Amendment †

Communication Department, Technical University of Cluj-Napoca, 400027 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper: Domuta, I.; Palade, T.P.; Puschita, E.; Pastrav, A. Localization in 802. 15.4z Standard. In Proceedings of the 2020 International Workshop on Antenna Technology (iWAT), Bucharest, Romania, 25–28 February 2020; pp. 1–4, doi:10.1109/iWAT48004.2020.1570615511.
Sensors 2020, 20(18), 5422; https://doi.org/10.3390/s20185422
Submission received: 9 July 2020 / Revised: 18 September 2020 / Accepted: 20 September 2020 / Published: 22 September 2020
(This article belongs to the Special Issue Antennas and Propagation)

Abstract

:
Due to the known issue that the ranging in the 802.15.4™-2015 standard is prone to external attacks, the enhanced impulse radio (EiR), a new amendment still under development, advances the secure ranging protocol by encryption of physical layer (PHY) timestamp sequence using the AES-128 encryption algorithm. This new amendment brings many changes and enhancements which affect the impulse-radio ultra-wideband (IR-UWB) ranging procedures. The timestamp detection is the base factor in the accuracy of range estimation and inherently in the localization precision. This paper analyses the key parts of PHY which have a great contribution in timestamp estimation precision, particularly: UWB pulse, channel sounding and timestamp estimation using ciphered sequence and frequency selective fading. Unlike EiR, where the UWB pulse is defined in the time domain, in this article, the UWB pulse is synthesized from the power spectral density mask, and it is shown that the use of the entire allocated spectrum results in a decrease in risetime, an increase in pulse amplitude, and an attenuation of lateral lobes. The paper proposes a random spreading of the scrambled timestamp sequence (STS), resulting in an improvement in timestamp estimation by the attenuation lateral lobes of the correlation. The timestamp estimation in the noisy channels with non-line-of-sight and multipath propagation is achieved by cross-correlation of the received STS with the locally generated replica of STS. The propagation in the UWB channel with frequency selective fading results in small errors in the timestamp detection.

1. Introduction

The wireless localization is a key part of many emerging technologies: internet of things (IoT), intelligent transportation systems (ITS), autonomous robots, or unmanned aerial vehicles. For many critical applications, localization accuracy is a basic requirement of localization systems. Due to its large bandwidth, the impulse radio ultra-wideband (IR-UWB) technology provides the best precision in range measurement by time-of-flight (ToF) estimation.
The basic feature in the accuracy of estimating the ToF is the shape of the pulse, more precisely the speed of increase of the pulse front. Nowadays the IEEE Task Group 4z, (TG4z), [1] is working on the enhanced impulse radio (EiR) project focused on localization safety improvement. This new amendment proposes a new UWB reference pulse and a time domain mask. In this paper, the UWB pulse is synthesized from power spectral density (psd) specified by the 802.15.4-2015 standard [2]. Two shapes of pulses were synthesized and compared, the first pulse being synthesized using only the central lobe of the psd mask and the second one being synthesized from the entire allocated spectrum.
The EiR amendment proposes that the timestamp estimation is validated by the cross-correlation of locally generated STS replica with received STS sequence. To avoid the interferences, the pulses are spread out on a symbol. This article proposes a supplementary spreading by a bit position modulation with a randomly generated sequence and shows, by simulated experiments, that lateral lobes of cross-correlation are mitigated by this modulation.
The behavior of the proposed methods is analyzed in a noisy radio channel with non line-of-sight (NLOS) and multipath propagation. The channel impulse response is estimated and subsequently used for the generation of a local replica of STS.
The main contributions of this article are that it:
-
demonstrates that the inclusion of lateral lobes with very low power spectral density (−51.3 dBm and −59.3 dBm) in pulse synthesis leads to a pulse with a tighter shape and a steeper rising edge than the pulse recommended by EiR;
-
shows that, by a random spreading of STS, that it results in easier extraction of the main lobe by mitigating lateral lobes of cross-correlation.
This paper is organized as follows. Section 2 presents state of the art research in the field. In Section 3, the UWB pulses are synthesized. Section 4 presents the random spreading of STS sequence and timestamp estimation in noisy channels and NLOS propagation. All sections incorporate simulated experiments, and because the simulation results from a subsection are used in the subsequent ones, they will be presented along with the theoretical aspects in the corresponding subsection.

2. Literature Review

The UWB radio holds a large bandwidth, but the harmonized standards [3] impose upper limits for power spectral density (−41.3 dBm/MHz, which is under the noise floor), resulting in great difficulty in the extraction of signal from noise. The research in [4,5] shows the presence of intra-symbol interference (IASI), inter-symbol interference (ISI), and multipath interference (MUI). In order to minimize the interferences, the pulse should have a small duration of the leading lobe and a high attenuation of the side lobes. It must be noted that this small duration of the main lobe can lead to an excess in bandwidth. The pulse shape has to be a compromise between the regulatory compliance, the need of low voltage and low power supply, low duration for maximization of data rate, ranging accuracy, and minimalization of interferences. In most cases, the UWB pulse is synthesized from Gaussian impulse [6], its derivatives [7], or a linear combination of Gaussian pulses [8]. Keshavarz et al. [9] infer the weights of the derivatives in the impulse structure by particle swarm optimization (PSO) algorithm, and PSO algorithm is used for optimization of the architecture of a UWB transmitter [10]. A linear combination of Gaussian monocycles with weight optimization by semidefinite programming is used for pulse synthesis [11]. Baranauskas and Zelenin present a direct waveform synthesis of UWB pulse by high speed DAC [12].
The EiR amendment specifies two pulses [13] as boundaries for the UWB pulse and a time domain mask [14] as a constraint for pulse shape. The UWB pulse synthesized from the entire allocated spectrum falls into the time domain mask, has high energy, and can be used as a reference in UWB pulse design.
In the 802.15.4 standard [2], the UWB PHY is specified in detail and, the transceivers manufactured in this technology are widely used in localization. However, several researches show that the range measurement in the current technology is prone to external attacks. Francillon et al. [15] present a relay attack, Taponecco et al. [16] show a delay attack and Singh et al. [17] propose a modulation scheme that secures the distance measurement against relay attack. The EiR project [1] brings a lot of improvements, including UWB reference pulse shape, preamble symbols revision, addition of scrambled timestamp sequence for secure ranging, an increase in data rate and PHY payload length, and the modification and addition of a new MAC primitive for key management. An overview of the EiR standard is presented in the work of Sedlacek et al. [18].
This article is limited to estimating the timestamp, without going into detail regarding ranging or location methods. Alarifi et al. perform a deep analysis of ultrawide band indoor positioning Technologies [19]. Several works deal with the wireless localization in internet of things (IoT) [20,21,22] and many research depict the localization in vehicular technologies [23,24,25].

3. UWB Pulse Synthesis

The pulse shape plays an essential role in the ranging accuracy and in reaching the maximum distance, while maintaining regulatory compliance. The 802.15.4-2015 standard, hereinafter called ‘old standard’ has been defined as a root raised cosine reference pulse. As this pulse has a precursor, it can mask the attenuated first path signal. The TG4 proposes [13] that the transmitted pulse shape p(t) to be constrained by the time domain mask, specified by the standard. The EiR specifies that the pulse risetime, 10–90%, for 500 MHz channels has to be maximum 2 ns.
This paper proposes the synthesis of UWB pulse from the compliant power spectral density mask (psd) [2] by Kolmogorov factorization (detailed in Appendix A) [26], because this method provides a minimum phase pulse, as it is specified by the EiR. The old standard specifies the psd mask as it is depicted in Figure 1a, namely the trace 802.15.4a mask. The power is expressed in Watts, in order to get the pulse amplitude in volts (1 Ω load). For pulse synthesis, a new spectral mask is designed, trace 1.5 ns risetime impulse, using a raised cosine profile
H = { p s d H | f f c | < ( 1 β ) F c   ( p s d H + p s d L 1 ) ( 1 + cos [ π 2 β ( f f c F c 1 + β ) ] ) + p s d L ( 1 β ) F c | f f c | ( 1 β ) F c 0 o t h e r w i s e
where central frequency f c = 0 , cutoff frequency F c = 315   MHz , roll-off factor β = 0.25 and p s d H   a n d   p s d L are the high and low value of psd mask. The sampling frequency is f s = 10   GHz and the window length is N = 10 4 samples for 1   MHz frequency resolution.
Different from standard [3], where the pulse energy is averaged on a 1   ms interval, in the paper, the mean power is computed on a 1   μ s interval, considering that the pulse repetition frequency P R F = 1   MHz . Accordingly, the pulse amplitude is inferred based on this PRF. The EiR defines many mean PRFs and, in order to comply with the regulations, the determined amplitude has to be scaled with P R F   [ MHz ] . Furthermore, the regulation imposes the pulse peak power to a value that shall not exceed 0   dBm   on   50   MHz bandwidth, and the impulse has to respect this restriction too.
Figure 1b shows the synthesized pulse, trace synthesized 1.5 ns impulse, compared to EiR compliant pulses for the 499.2 MHz bandwidth (i.e., 1.7 ns risetime impulse and 1.2 ns risetime impulse traces). The synthesized pulse is situated between recommended pulse limits, so it respects the specifications. The traces 1.7 ns risetime and 1.2 ns risetime in Figure 1a show that the pulses suggested by the EiR do not fit exactly in the standard psd mask, the first exceeding the maximum psd and the second exceeding the bandwidth.
Recently, TG4z has defined a time domain mask for UWB impulse [12], as illustrated in Figure 2b. Based on this mask, it is appropriate to search for a new UWB pulse shape which falls in this mask for impulse energy maximization and risetime reduction. The Cramer–Rao lower bound in ToF estimation is inversely proportional to effective bandwidth [27], so it is convenient to use the lateral lobes of low power for pulse synthesis. In order to do this, a new psd mask is designed using a sum of raised cosine profiles, Figure 2, having the following parameters: ( F c 1 = 315   MHz ,   β 1 = 0.05 ) ;   ( F c 2 = 400   MHz ,   β 2 = 0.05 ) ;   ( F c 3 = 500   MHz ,   β 3 = 0.05 ) .
Figure 2b shows that the pulse synthesized from complete mask has a smaller risetime and a stronger attenuation of lateral lobes compared to the pulse synthesized from partial psd mask. As such, the pulse synthesized from complete mask allows for more precise timestamp estimation and less interferences. The pulse synthesized from complete mask has the amplitude 0.187   V / 1   Ω / 1   μ s , higher energy, low risetime T 10 % 90 % = 1.2   ns and smaller lateral lobes than the pulse synthesized from partial psd mask which has an amplitude of 0.164   V / 1   Ω / 1   μ s and a risetime of T 10 % 90 % = 1.5   ns . The pulse risetime plays a key role in timestamp estimation precision. Therefore, it is preferable to choose the pulse synthesized from complete mask as UWB reference pulse.
The degree of spectrum usage is evaluated by normalized effective signal power (NESP) [6]
NESP = | S ( f ) | 2 d f M ( f ) d f
where | S ( f ) | 2 is spectral density of the impulse and M ( f ) is allocated spectral mask.
Because the signal does not have a uniform distribution on the entire spectrum, the effective bandwidth β is defined
β = f 2 | S ( f ) | 2 d f | S ( f ) | 2 d f
Usually, the quality of the UWB pulse synthesis is defined in terms of occupied bandwidth and pulse duration. Table 1 compares the quality parameters of proposed pulse with P802.15.4z reference pulse.
The minimum uncertainty, σ , for range estimation in the time of arrival method is quantified by Cramér–Rao lower bound (CRLB) [28]
σ = c β 8 π 2 S N R
where β is effective bandwidth, c is speed of the light, and S N R is signal to noise ratio.
Relation (100) clearly shows that for minimizing the uncertainty the effective bandwidth has to be maximized, but within the limit of regulations.

4. Timestamp Estimation in the 802.15.4z Standard

4.1. Timestamp Estimation by Random Spreading of STS

The new standard brings in a new physical protocol data unit (PPDU) structure by incorporating the STS for secure ranging. The STS is encrypted using the AES-128 algorithm, the time of arrival estimation is achieved on STS, and the range measurement is validated only if the received STS cross correlated with the locally generated reference exceeds the “match level” [29] threshold.
The default PHY frame format proposed by EiR [29] is depicted in Figure 3, where SHR is the synchronization header (preamble), STS is the scrambled timestamp sequence, and PHR is the PHY header.
To avoid the inter-pulse interferences, [30] stipulates that every component B k of Ipatov ternary symbol (ITS), or STS, is spread out on a symbol of length δ L by n = 0 δ L 1 δ ( n )   p ( t n T c h ) ; where δ ( n ) is Kronecker delta, p ( t ) is UWB pulse and T c h = 2   n s is chip duration. To mitigate the side lobes of correlation, this paper proposes a supplementary spreading by a randomly generated sequence S k . The sequence of length N is:
s ( t ) = k = 0 N 1 B k n = 0 δ L 1 ( δ ( n ) · p ( t ( k · δ L + S k n s + n ) T c h ) ) .
The STS sequence of length N = 128 , s(t), is generated by taking B k = 1 2 A k , A being the result of AES encryption. The peak pulse repetition frequency (PRF) is 499.2 MHz, the mean PRF is 62.4 MHz, resulting δ L = 8   chips . Figure 4a depicts the STS sequences: STS standard spread out is generated according to the standard specifications, S k = 0 , STS reference signal is the signal needed to achieve correlation and for STS randomly spread out, the S k is generated based on a linear feedback shift register with the characteristic polynomial x 3 + x + 1 . The cross-correlation of the reference signal and STSs are shown in Figure 4b. The ratio between the main lobe and maximum side lobe is η 10 for cross-correlation randomly spread out and η 3.33 for cross-correlation standard spread out, which proves the efficiency of random spreading.
Figure 4a shows that the random spreading increases the risk of IASI. Therefore, it is necessary to perform the analysis of timestamp estimation for propagation in noisy channel and for multipath propagation.

4.2. Sounding the Channel with Multipath Propagation

The channel sounding is the estimation of the channel impulse response (CIR) using preamble sequence, in order to remove the noise and design the channel equalizer.
The frequency-dependent path gain G ( d , f ) is modeled considering an isotropic radiation pattern, with the “antenna attenuation factor” [31] of ½:
G ( d , f ) = P R x ( d , f ) P T x ( f ) = 1 2   G 0 η T x η R x ( f / f c ) 2 ( κ + 1 ) ( d / d 0 ) n   ,
where P R x ( d , f ) , P T x ( f ) are received and transmitted power, η T x , η R x are the transmission and reception antenna gains, G 0 is the path gain at reference distance d 0 , d is the distance between transmitter and receiver, n is the path gain exponent, f c is the carrier frequency, f is the frequency and κ is the frequency decaying factor.
The frequency decaying factor follows the Friis equation and, in Equation (6), it has the value κ = 0 . The path loss varies from n = 1.2 in industrial LOS, n = 1.76 in outdoor LOS, to n = 4.58 in residential NLOS. It is noted that multipath propagation in the industrial environment leads to an increase in the path gain.
Using the Saleh-Valenzuela statistical model, the propagation paths are designed as the sum of clusters, every cluster having multiple rays. The impulse response h m ( t ) is
h m ( t ) = l = 0 L r = 0 K a r , l e x p ( j r , l ) δ ( t T l τ r , l ) ,
where a r , l is the tap weight of r ray in cluster l , r , l is the ray phase, T l is the delay of l cluster and τ r , l is the delay of r ray relative to cluster l front. The intervals between the time arrivals of the clusters is modeled as Poisson process with the arrival rate Λ , and the ray delays inside the cluster are modeled as a mixture of two Poisson processes with arrival rates λ 1 ,   λ 2 and mixing weight β . The mean power of arriving clusters follows an exponential decay with time constant Γ , having a normal distribution around the mean value σ c l u s t e r , and the cluster shape also bears to an exponential decay with time constant γ .
In the above CIR model, the first path has the highest energy. In non-line-of-sight propagation, there are cases when the first path is strongly attenuated. For such situations, [31] proposes a new modeling for the first path.
E { a k , 0 } ( 1 χ e x p ( τ k , l γ r i s e ) ) e x p ( τ k , l γ 1 ) ,
where χ describes the attenuation of the first path, γ r i s e determines how fast the power delay profile (PDP) increases, and γ 1 determines the profile decay. By joining Equations (7) and (8), the CIR, h(t) can be found.
The preamble sequence consists of a string of 32 or 64 Ipatov ternary symbols, every symbol having 91 elements with 81 non-zero elements [30]. An Ipatov symbol has “perfect” periodic autocorrelation, i.e., all side lobes of autocorrelation are zero, and using Wiener-Hopf equation, by cross-correlation of the received signal y ( t ) with the input Ipatov sequence I ( t ) results immediately the CIR, h ( t )
y ( t ) I ( t ) = h ( t ) * I ( t ) I ( t ) = h ( t ) * R I , I ( t ) ,
where autocorrelation of ITS is: R I , I ( 0 ) = N ; R I , I ( t ) = 0   for   t 0 .
For outdoor, NLOS and multipath propagation environment, the channel PDP, h ( t ) , is modeled based on Equations (7) and (8), with the parameters retrieved from [31], and synthesized in Table 2.
For channel sounding, the EiR standard specifies a mandatory preamble sequence of 32 or 64 ITSs, the Ipatov symbol having a number of 91 elements with 10 zero elements. In this paper, the CIR is estimated only using one ITS with a total of 57 elements including 8 zero elements [32]. The transmitted sequence of pulses, i ( t ) , is generated based on Equation (5), with B k = I k , where I k is the kth element of ITS. The received signal y i ( t ) = h ( t ) * i ( t ) is the convolution of PDP with the emitted sequence. The CIR estimation, h ( t ) ^ , is achieved by the cross-correlation of the received signal y ( t ) with the Ipatov symbol I . Figure 5 shows that the CIR estimation, estimated PDP, is close enough to the true PDP. By successive transmission of ITSs contained in the preamble, the estimations are accumulated and averaged, resulting in SNR reduction and CIR estimation improvement.
This CIR estimate will be used in the following paragraphs, for analysis of timestamp estimation in channels with multipath propagation. This is a classical channel model, more detailed models are being published in the recent research [33,34].

4.3. Timestamp Estimation in Channel with NLOS and Multipath Propagation

The EiR specifies that STS consists of 32, 64, 128 AES-128 sequences, successively transmitted, but in this section only one AES sequence is considered for analysis the impact of multipath propagation.
The STS is generated based on Equation (5) with random spread out, and for propagation simulation it is convoluted with h ( t ) .
In the usual way, the received signal is passed through an equalizer filter and cross-corelated with the STS reference s r ( t ) . In this article, an easier way for timestamp detection is proposed, that is, to generate a virtual propagated STS reference by convolution of the STS reference with the estimated CIR, h ( t ) ^ , and cross-correlation of received sequence y s ( t ) with locally generated replicas y r ( t ) ^ as depicted in Algorithm 1.
Algorithm 1: Timestamp estimation in multipath propagation
Inputs: STS reference (AES-128 sequence), s r ( t ) ; Received STS, y s ( t )
 1. Generate locally replica of STS:
y r ( t ) ^ = h ( t ) ^ * s r ( t )
 2. Compute the received sequence:
y s ( t ) = h ( t ) * s ( t )
 3. Timestamp estimation:
r r s ( t ) = y r ( t ) ^ y s ( t )
Output: return r r s ( t )
Figure 6 shows that by cross-correlating the received STS, y s ( t ) , with the STS reference, s r ( t ) , the timestamp is not detectable, (see the correlation with STS reference trace), and that the cross-correlation of received signal with the locally generated replica, y r ( t ) ^ , the side lobes are strongly attenuated (see the correlation with STS reference convoluted with PDP trace).

4.4. Timestamp Estimation in Noisy Channel with NLOS and Multipath Propagation

The EiR standard specifies that the typical range of radio is 100 m. For timestamp estimation in noisy channel we consider that the transmitter is situated at 22 m and the noise source is situated at reference distance of 1 m.
The transmitter emits with maximum compliant power (Figure 2a) and the noise source emits with floor noise level (i.e., −41 dBm/MHz) [35]. In order to mitigate the effect of noise, consider that the receiver has a 1 GHz passband filter on the input. The standard deviation of noise is σ = 0.0228   V and the signal to noise ratio on the receiver is S N R = 38   dB . In this situation, the timestamp is undetectable from only a single AES-128 sequence (Figure 7b, trace 1 STS), so STS will consist of multiple AES sequences as the EiR standard specifies.
The timestamp detection is detailed in Algorithm 2 and the results are shown in Figure 7.
Algorithm 2: Timestamp estimation in noisy channel
Inputs: STS reference (AES-128 sequence), s r ( t ) ; Received STSs, y s i ( t ) ; Number of STS, N ;
 1. Generate locally replica of STS:
y r ( t ) ^ = h ( t ) ^ * s r ( t )
 2. Compute the received sequences:
f o r   i [ 0 , N 1 ]     d o :
      y s i ( t ) = h ( t ) * s ( t ) + n i ( t )
 3. Average the received sequences:
y s ( t ) = 1 N i = 0 N 1 y s i ( t )
 4. Timestamp estimation:
r r s ( t ) = y r ( t ) ^ y s ( t )
 Where n i ( t ) is the noise with standard deviation σ
Output: return r r s ( t )

4.5. UWB Channel with Frequency Selective Fading

Depending on antenna design, or antenna position relative to external objects, it is possible to encounter frequency selective fading [35,36]. To analyze the impact of frequency fading in timestamp detection, consider UWB channel 9 with the spectral mask profile power spectral density, illustrated in Figure 8, having a selective fading of 12   dBm at carrier frequency and frequency decaying factor of κ = 0.13 . The channel frequency response, H ( f ) , is achieved by Kolmogorov factorization (Appendix A).
To estimate the CIR, H ( f ) ^ , the Ipatov symbol, i ( t ) , is shifted on to the carrier, f c , by complex modulation with exp ( j 2 π f c t ). The cross-correlation is performed in the frequency domain, and to avoid the circular correlation, the series is padded with zeros. The propagation of STS in the faded channel is simulated by the shifting of the STS on to the carrier and by the convolution with H ( f ) . The timestamp estimation is detailed in Algorithm 3.
Algorithm 3: Timestamp estimation in channel with frequency selective fading
Inputs: Ipatov symbol, i ( t ) ; STS, s ( t ) ;
   CIR estimation:
 1. Lift up the Ipatov symbol on the carrier:
i c ( t ) = e j 2 π f c t · i ( t )
 2. Complete the series with zero and perform FFT.
I c ( f ) = F F T ( i c ( t ) )
 3. Convolve the Ipatov symbol with CIR:
Y I ( f ) = H ( f ) · I c ( f )
 4. CIR estimation by cross-correlation:
H ( f ) ^ = Y I ( f ) · I c ( f ) ¯ I c ( f ) · I c ( f ) ¯
   Timstamp estimation from STS:
 5. Shift STS to carrier, complete with zero and perform FFT:
S c ( f ) = F F T ( e j 2 π f c t · s ( t ) )
 6. Generate locally replica of STS:
S r ( f ) ^ = H ( f ) ^ · S c ( f )
 7. Simulate the received STS by convolution:
Y S ( f ) = H ( f ) · S c ( f )
 8. Cross-correlation in the frequency domain:
R r s ( f ) = S r ( f ) ^ · Y S ( f ) ¯
 9. Correlation in time domain:
r r s ( t ) = I F F T ( R r s ( f ) )
Output: return r r s ( t )
Figure 9a, psd at the emission, shows that Ipatov symbol has a uniform distribution over the entire bandwidth, that it follows the compliant spectral mask (Figure 2a), and that after propagation, psd after fading, it borrows CIR spectral mask.
The propagation in a channel with frequency fading leads to a small error in timestamp estimation, as shown in Figure 9b, correlation with STS reference, and this error is cancelled if the received STS is correlated with the locally generated replica, correlation with STS convolved with H.
It should be highlighted that the spectral leakage in the FFT leads to the introduction of nonuniformity in the estimated CIR spectrum, Figure 10a, resulting in incorrect results regarding timestamp estimation, Figure 10b.

5. Results

The NESP shows the usage efficiency of the allocated bandwidth. Table 3 displays the occupied bandwidth and NESP for the pulse synthesized from entire spectral mask which is compared with previously published related works.
The uncertainty in range estimation for noisy channel with multipath propagation is computed based on Equation (4). The effective bandwidth is computed based on Equation (3). The leading lobe in Figure 7b is extracted, the samples series is completed with zero and then transformed in frequency domain. The effective bandwidth is nearly constant β = 234   MHz .
The propagation loss is determined considering outdoor LOS channel with path loss exponent n = 1.76 ,
L = 10 · n · log ( d d o )
where d 0 = 1   m is reference distance.
The noise has psd at floor level P n 41   dBm and 7   GHz bandwidth.
The SNR and minimum uncertainty are depicted in Table 4.

6. Discussion

The synthesis of UWB impulse from a complete psd mask (Figure 2a) results in a risetime reduction, amplitude increasing, and the attenuation of side lobes relative to the impulse proposed by EiR (Table 1). Furthermore, the pulse psd fully complies with the standard spectral mask. Therefore, this impulse is advisable as a UWB reference pulse.
The random spreading of STS leads to an attenuation of lateral lobes of cross-correlation, but also leads to an increase in the probability of intra-symbol interferences. Thus, more investigations should be performed for testing these interferences in noisy channels with multipath propagation.
The cross-correlation in a noisy channel with multipath propagation displays a wide main lobe, which leads to a decrease in the accuracy of the timestamp estimation. This result seems to be due to the fact that the AES-120 sequence does not have uniform spectral distribution in the entire frequency band. The whitening of the AES sequence or of the entire STS by additional randomization could lead to the main lobe narrowing, and this is a subject of interest in our research.

Author Contributions

Conceptualization, I.D. and T.P.P.; methodology, E.P.; software, I.D.; investigation and writing, I.D., E.P. and A.P.; visualization, A.P.; supervision, T.P.P.; funding acquisition, T.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This study was accomplished within the Centre of competence for wireless intrasatellite technologies (IntraSAT-Tech), Technical University of Cluj-Napoca.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Kolmogorov Factorization
Algorithm A: Kolmogorov Factorization
Inputs: correlation { r t } , psd S ( ω )
 1. Inputs test:
{ r t } ? Continue: Goto 3
 2. Compute psd:
S ( ω ) = F F T ( { r t } )
 3. Compute cepstrum coefficients:
R e H ( ω ) = 1 2 ln S ( ω )
 4. Determine imaginary part of log filter by Hilbert transform:
R E H ( ω ) = F F T ( R e H ( ω ) )
H l i = { 0                                             i f     i 0 j   R E H i               i f   i > 0 j   R E H i                     i f   i < 0
I m H ( ω ) = I F F T ( H l )
 5. Transfer function in frequency domain:
H ( ω ) = e R e H + j   I m H
 6. Transfer function in time domain of moving average type:
h = I F F T ( H ( ω ) )
 7. Autoregressive function:
a = h 0 h
Outputs: return H ( ω ) ,   h ,   a

References

  1. IEEE 802.15 WPAN™ Task Group 4z Enhanced Impulse Radio. Available online: http://www.ieee802.org/15/pub/TG4z.html (accessed on 19 October 2019).
  2. IEEE Computer Society. IEEE Standard for Low-Rate Wireless Networks Amendment 2: Ultra-Low Power Physical Layer; IEEE Computer Society: Washington, DC, USA, 2016; Volume April. [Google Scholar]
  3. ETSI. Electromagnetic Compatibility and Radio Spectrum Matters (ERM); Short Range Devices (SRD) using Ultra Wide Band (UWB); Transmission Characteristics; Part 2: UWB Mitigation Techniques; ETSI: Sophia Technology Park, France, 2014. [Google Scholar]
  4. Luo, C.; Wu, X.; Cao, Y. Multipath interference analysis of IR-UWB systems in indoor office LOS environment. In Proceedings of the 2011 6th International ICST Conference on Communications and Networking in China, CHINACOM 2011, Harbin, China, 17–19 August 2011; pp. 846–850. [Google Scholar]
  5. Stanciu, M.I.; Azou, S.; Rədoi, E.; Serbənescu, A. A statistical analysis of multipath interference for impulse radio UWB systems. J. Frankl. Inst. 2015, 352, 5952–5967. [Google Scholar] [CrossRef]
  6. Li, B.; Zhou, Z.; Zou, W.; Li, D.; Zhao, C. Optimal waveforms design for ultra-wideband impulse radio sensors. Sensors 2010, 10, 11038–11063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Sharma, A.; Sharma, S.K. Spectral efficient pulse shape design for UWB communication with reduced ringing effect and performance evaluation for IEEE 802.15.4a channel. Wirel. Netw. 2019, 25, 2723–2740. [Google Scholar] [CrossRef]
  8. Liu, X.; Premkumar, A.B.; Madhukumar, A.S. Pulse shaping functions for UWB systems. IEEE Trans. Wirel. Commun. 2008, 7, 1512–1516. [Google Scholar] [CrossRef]
  9. Keshavarz, S.N.; Hamidi, M.; Khoshbin, H. A PSO-Based UWB Pulse Waveform Design Method. In Proceedings of the 2010 Second International Conference on Computer and Network Technology, Bangkok, Thailand, 23–25 April 2010; pp. 249–253. [Google Scholar] [CrossRef]
  10. Ben Issa, D.; Samet, M. Design and Optimization of Dual-Band Energy-Efficient OOK UWB Transmitter Via PSO Algorithm. J. Circuits Syst. Comput. 2019. [CrossRef] [Green Version]
  11. Wu, X.; Tian, Z.; Davidson, T.N.; Giannakis, G.B. Optimal waveform design for UWB radios. IEEE Trans. Signal Process. 2006, 54, 2009–2021. [Google Scholar] [CrossRef] [Green Version]
  12. Baranauskas, D.; Zelenin, D. A 0.36W 6b up to 20GS/s DAC for UWB Wave Formation. In Proceedings of the 2006 IEEE International Solid State Circuits Conference—Digest of Technical Papers, San Francisco, CA, USA, 8–10 February 2006; pp. 2380–2389. [Google Scholar] [CrossRef]
  13. McLaughin, M.; Niewczas, J.; Verso, B. IEEE P802.15-19-0443-01-004z, Text to Address Comment id r1-0820. 2018. Available online: http://www.ieee802.org/15/pub/TG4z.html (accessed on 19 October 2019).
  14. McLaughin, M.; Hammerschmidt, J.; Ibrahim, B.; Verso, B. IEEE P802.15-20-0089-01-004z, Pulse Shape Text Changes for HRP UWB PHY. 2020. Available online: http://www.ieee802.org/15/pub/TG4z.html (accessed on 1 June 2020).
  15. Francillon, A.; Danev, B.; Capkun, S. Relay Attacks on Passive Keyless Entry and Start Systems in Modern Cars. Netw. Distrib. Syst. Secur. Symp. 2011, 431–439, 8529521. [Google Scholar]
  16. Taponecco, L.; Perazzo, P.; D’Amico, A.A.; Dini, G. On the feasibility of overshadow enlargement attack on IEEE 802.15.4a distance bounding. IEEE Commun. Lett. 2014, 18, 257–260. [Google Scholar] [CrossRef]
  17. Singh, M.; Leu, P.; Capkun, S. UWB with Pulse Reordering: Securing Ranging against Relay and Physical-Layer Attacks. NDSS 2019. [Google Scholar] [CrossRef]
  18. Sedlacek, P.; Slanina, M.; Masek, P. An overview of the IEEE 802.15.4z standard its comparison and to the existing UWB standards. In Proceedings of the 2019 29th International Conference Radioelektronika, RADIOELEKTRONIKA 2019—Microwave and Radio Electronics Week, MAREW, Pardubice, Czech Republic, 16–18 April 2019; pp. 1–6. [Google Scholar] [CrossRef]
  19. Alarifi, A.; Al-Salman, A.; Alsaleh, M.; Alnafessah, A.; Al-Hadhrami, S.; Al-Ammar, M.A.; Al-Khalifa, H.S. Ultra wideband indoor positioning technologies: Analysis and recent advances. Sensors (Switzerland) 2016, 16, 707. [Google Scholar] [CrossRef]
  20. Khelifi, F.; Bradai, A.; Benslimane, A.; Rawat, P.; Atri, M. A Survey of Localization Systems in Internet of Things. Mob. Netw. Appl. 2019, 24, 761–785. [Google Scholar] [CrossRef]
  21. Shit, R.C.; Sharma, S.; Puthal, D.; Zomaya, A.Y. Location of Things (LoT): A review and taxonomy of sensors localization in IoT infrastructure. IEEE Commun. Surv. Tutor. 2018, 20, 2028–2061. [Google Scholar] [CrossRef]
  22. E Silva, P.F.; Kaseva, V.; Lohan, E.S. Wireless positioning in IoT: A look at current and future trends. Sensors (Switzerland) 2018, 18, 2470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kuutti, S.; Fallah, S.; Katsaros, K.; Dianati, M.; Mccullough, F.; Mouzakitis, A. A Survey of the State-of-the-Art Localization Techniques and Their Potentials for Autonomous Vehicle Applications. IEEE Internet Things J. 2018, 5, 829–846. [Google Scholar] [CrossRef]
  24. Balico, L.N.; Loureiro, A.A.F.; Nakamura, E.F.; Barreto, R.S.; Pazzi, R.W.; Oliveira, H.A.B.F. Localization prediction in vehicular Ad hoc networks. IEEE Commun. Surv. Tutor. 2018, 20, 2784–2803. [Google Scholar] [CrossRef]
  25. Del Peral-Rosado, J.A.; Seco-Granados, G.; Kim, S.; López-Salcedo, J.A. Network Design for Accurate Vehicle Localization. IEEE Trans. Veh. Technol. 2019, 68, 4316–4327. [Google Scholar] [CrossRef]
  26. Scharf, L. Statistical Signal Processing: Detection, Estimation and Time Series Analysis, 1st ed.; Addison-Wesley Publishing Company, Inc.: Boston, MA, USA, 1991; pp. 424–429. ISBN 0-201-19038-9. [Google Scholar]
  27. Qi, Y.; Kobayashi, H.; Suda, H. On time-of-arrival positioning in a multipath environment. IEEE Trans. Veh. Technol. 2006, 55, 1516–1526. [Google Scholar] [CrossRef]
  28. Álvarez, R.; Díez-González, J.; Alonso, E.; Fernández-Robles, L.; Castejón-Limas, M.; Perez, H. Accuracy Analysis in Sensor Networks for Asynchronous Positioning Methods. Sensors 2019, 19, 3024. [Google Scholar]
  29. Niewczas, J.; Verso, B.; Fagan, T.; Leong, F.; Hammerschmidt, J.; Ibrahim, B.; Sasoglu, E.; Knobloch, D.; Reisinger, T. IEEE P802.15-19-0134-00-004z, Security vs. Sequence Length Considerations. 2019. Available online: http://www.ieee802.org/15/pub/TG4z.html (accessed on 21 October 2019).
  30. Verso, B.; Leong, F.; Hammerschmidt, J.; Niewczas, J.; Ibrahim, B.; Shah, T.; Reisinger, T.; Daniel Knobloch, D. IEEE P802.15-18-0477-00-004z, HRP UWB PHY Enhanced Mode Converged Consensus. 2018. Available online: http://www.ieee802.org/15/pub/TG4z.html (accessed on 19 October 2019).
  31. Molisch, A.F.; Cassioli, D.; Chong, C.C.; Emami, S.; Fort, A.; Kannan, B.; Karedal, J.; Kunisch, J.; Schantz, H.G.; Siwiak, K.; et al. A comprehensive standardized model for ultrawideband propagation channels. IEEE Trans. Antennas Propag. 2006, 54, 3151–3166. [Google Scholar] [CrossRef] [Green Version]
  32. Li, X.; Fan, P.; Mow, W.H. Existence of ternary perfect sequences with a few zero elements. In Proceedings of the 5th International Workshop on Signal Design and Its Applications in Communications, IWSDA’11, Guilin, China, 10–14 October 2011; pp. 88–91. [Google Scholar]
  33. Kram, S.; Stahlke, M.; Feigl, T.; Seitz, J.; Thielecke, J. UWB channel impulse responses for positioning in complex environments: A detailed feature analysis. Sensors (Switzerland) 2019, 19, 5547. [Google Scholar] [CrossRef] [Green Version]
  34. Nwadiugwu, W.P.; Kim, D.S. Ultrawideband Network Channel Models for Next-Generation Wireless Avionic System. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 113–129. [Google Scholar] [CrossRef]
  35. International Telecommunication Union (ITU). Presentation: Advanced wireless technologies and spectrum management. In Radio Spectrum Management for a Converging World; ITU: Geneva, Switzerland, 2004; Available online: http://handle.itu.int/11.1002/pub/800c7e3d-07f64a9d-en (accessed on 25 October 2019).
  36. Otim, T.; Díez, L.E.; Bahillo, A.; Lopez-Iturri, P.; Falcone, F. Effects of the body wearable sensor position on the UWB localization accuracy. Electronics 2019, 8, 1351. [Google Scholar] [CrossRef] [Green Version]
Figure 1. UWB impulse synthesis: (a) Spectral masks; (b) The pulses in time domain.
Figure 1. UWB impulse synthesis: (a) Spectral masks; (b) The pulses in time domain.
Sensors 20 05422 g001
Figure 2. UWB impulse synthesis from complete spectral mask: (a) Spectral mask from multiple raised cosine profiles; (b) Time domain mask and synthesized pulses.
Figure 2. UWB impulse synthesis from complete spectral mask: (a) Spectral mask from multiple raised cosine profiles; (b) Time domain mask and synthesized pulses.
Sensors 20 05422 g002
Figure 3. PHY frame format in the 802.15.4z amendment.
Figure 3. PHY frame format in the 802.15.4z amendment.
Sensors 20 05422 g003
Figure 4. Timestamp estimation by STS: (a) UWB pulse position in STS sequence; (b) cross-correlation of STS sequences with reference signal.
Figure 4. Timestamp estimation by STS: (a) UWB pulse position in STS sequence; (b) cross-correlation of STS sequences with reference signal.
Sensors 20 05422 g004
Figure 5. CIR estimation for channel with multipath propagation.
Figure 5. CIR estimation for channel with multipath propagation.
Sensors 20 05422 g005
Figure 6. Timestamp estimation in channel with NLOS and multipath propagation.
Figure 6. Timestamp estimation in channel with NLOS and multipath propagation.
Sensors 20 05422 g006
Figure 7. Timestamp estimation in noisy channel with multipath propagation.
Figure 7. Timestamp estimation in noisy channel with multipath propagation.
Sensors 20 05422 g007
Figure 8. Psd mask and channel response.
Figure 8. Psd mask and channel response.
Sensors 20 05422 g008
Figure 9. Timestamp estimation in channel with selective fading.
Figure 9. Timestamp estimation in channel with selective fading.
Sensors 20 05422 g009
Figure 10. Spectral leakage in timstamp estimation: (a) Nonuniformity in CIR estimation; (b) Erroneous estimation of timestamp.
Figure 10. Spectral leakage in timstamp estimation: (a) Nonuniformity in CIR estimation; (b) Erroneous estimation of timestamp.
Sensors 20 05422 g010
Table 1. The synthesis quality parameters.
Table 1. The synthesis quality parameters.
Synthesis MethodBandwidth β [MHz]NESP [%]Risetime
[ns]
Amplitude
[V]
802.15.4z reference pulse [11]53394.51.50.164
Pulse synthesized from entire spectral mask (Figure 2a)53598.51.20.187
Table 2. The UWB channel parameters for outdoor NLOS propagation.
Table 2. The UWB channel parameters for outdoor NLOS propagation.
Name [unit]SymbolValue
Path gain:
Reference distance [m] d 0 1
Gain at reference distance [dB] G 0 −73
Path gain exponent n 2.5
Frequency decaying factor κ 0.13
Power delay profile:
Expected number of clusters L ¯ 10.5
Clusters arrival rate [1/ns] Λ 0.0243
Clusters decay time constant [ns] Γ 104.7
Rays arrival rate [1/ns] λ 0.223
Intracluster decay time constant [ns] γ 9.3
First path attenuation χ 0.65
Path increase time constant [ns] γ r i s e 5
Table 3. UWB pulse parameters.
Table 3. UWB pulse parameters.
Synthesis MethodBandwidth β [GHz]NESP [%]
Pulse synthesized from entire spectral mask (Figure 2a) 0.53798.5
Gaussian with NESP maximization [6]7.16098.7
Linear combination of Gaussian pulses [8] 7.20091
VC oscillator with PSO optimization [10] 1.332 NA
Gaussian with semidefinite programing [11]1285.5
Table 4. CRLB for timestamp estimation by STS in noisy channels and multipath propagation.
Table 4. CRLB for timestamp estimation by STS in noisy channels and multipath propagation.
Distance [m]110204080160
SNR [dB]−10.2−27.8−33.1−38.4−43.7−49
CRLB [m]0.1665 1.25 2.314.257.7814.33

Share and Cite

MDPI and ACS Style

Domuta, I.; Palade, T.P.; Puschita, E.; Pastrav, A. Timestamp Estimation in P802.15.4z Amendment. Sensors 2020, 20, 5422. https://doi.org/10.3390/s20185422

AMA Style

Domuta I, Palade TP, Puschita E, Pastrav A. Timestamp Estimation in P802.15.4z Amendment. Sensors. 2020; 20(18):5422. https://doi.org/10.3390/s20185422

Chicago/Turabian Style

Domuta, Ioan, Tudor Petru Palade, Emanuel Puschita, and Andra Pastrav. 2020. "Timestamp Estimation in P802.15.4z Amendment" Sensors 20, no. 18: 5422. https://doi.org/10.3390/s20185422

APA Style

Domuta, I., Palade, T. P., Puschita, E., & Pastrav, A. (2020). Timestamp Estimation in P802.15.4z Amendment. Sensors, 20(18), 5422. https://doi.org/10.3390/s20185422

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