 Research
 Open Access
 Published:
A novel time of arrival estimation algorithm using an energy detector receiver in MMW systems
EURASIP Journal on Advances in Signal Processing volume 2017, Article number: 83 (2017)
Abstract
This paper presents a new time of arrival (TOA) estimation technique using an improved energy detection (ED) receiver based on the empirical mode decomposition (EMD) in an impulse radio (IR) 60 GHz millimeter wave (MMW) system. A threshold is employed via analyzing the characteristics of the received energy values with an extreme learning machine (ELM). The effect of the channel and integration period on the TOA estimation is evaluated. Several wellknown EDbased TOA algorithms are used to compare with the proposed technique. It is shown that this ELMbased technique has lower TOA estimation error compared to other approaches and provides robust performance with the IEEE 802.15.3c channel models.
Introduction
Accurate estimation of the time of arrival (TOA) in wireless systems is a challenging problem due to intersymbol interference and multipath fading. Sixty gigahertz millimeter wave (MMW) is a promising technology for accurate TOA estimation due to its unique advantages which include high time resolution [1], high multipath resolution [2], and robustness to interference [3]. In the 60GHz frequency band, an energy detector (ED) receiver is preferred for TOA estimation because it can be implemented easily in hardware and does not require precise channel estimation and synchronization as with a matched filter (MF) [4]. As shown in Fig. 1, a conventional ED only consists of a bandpass filter (BPF), a squaring operator, an integrator, and a decision device.
The TOA can be estimated based on the first energy sample to exceed a threshold. Various EDbased TOA estimation methods have been proposed [5,6,7,8,9,10,11]. Guvenc first proposed a TOA estimation method via analyzing the kurtosis characteristics of the energy of the received pulses [6]. To account for signal variations, a normalized threshold was employed in [7], which is based on the maximum and minimum energy values. In [8], a fixed threshold based on the maximum energy value was used. The cell averaging constant false alarm rate (CACFAR) method is developed for the TOA estimate by employing the ED receiver in [9].A weighted ED receiver was proposed in [10] to reduce the effects of noise. Although an extreme learning machine was used to obtain the threshold [11], this algorithm performs poorly with a conventional ED in low signaltonoise ratio (SNR) environments, particularly when there are noiseonly energy values. We can see all of these algorithms perform poorly in low SNR environments because of the sensitivity of the threshold to noise. Using a conventional ED, these obtained energy samples are only noise which will have a significant effect on the TOA estimation \( {\widehat{\tau}}_{\mathrm{TOA}} \).
To reduce the effects of noise and improve the estimation accuracy, in this paper, (1) an improved ED receiver was presented for TOA estimation which employs impulse radio (IR) 60 GHz MMW signals. The received signals were processed based on empirical mode decomposition (EMD) [12]. This method has been used to analyze nonstationary and nonlinear signals [13]. EMD adaptively decomposes a signal into a series of intrinsic mode functions (IMFs) in descending order of frequency and a residual trend mode [13]. The IMFs and residual trend mode can be used to represent the noise which can then be deleted, resulting in an improved SNR [14]. (2) An EDbased TOA estimation technique was proposed which employs extreme learning machine (ELM). ELM was used to resolve a regression problem with the characteristics of the received ED values as inputs and the detector threshold as the output. Results are presented which show that this denoising approach provides superior performance with the IEEE 802.15.3c channel models compared to several wellknown techniques. And the simulation results were presented which show that ELM can provide robust and reliable TOA estimates. Further, the proposed approach can be used with other wireless systems.
The remainder of this paper is organized as follows: Section 2 presents the system model. An improved ED receiver was discussed using EMD in Section 3. Section 4 discusses the signal characteristics, and the proposed EDbased TOA estimation technique is introduced in Section 5. Some performance results are given in Section 6, and Section 7 provides some conclusions.
System model
The IEEE 802.15.3c Line of Sight (LOS) and Nonline of Sight (NLOS) channel models have been proposed by the TG3c group for 60 GHz indoor residential, office, and library environments [21]. Several types of signals have been considered for 60 GHz systems [18]. Pulse position modulation (PPM) with a truncated sinc pulse has been shown to better conform to the Federal Communications Commission (FCC) regulations and provide a lower error probability than other signals. Thus, PPM with a truncated sinc pulse is considered here.
The received 60GHz MMW signal can be expressed as:
where E _{ b } is the energy per bit, (*) denotes convolution, N _{ s } is the number of pulses transmitted per data symbol, j is the frame index, and T _{ s } is the frame duration. The timehopping codes are c _{ j } ∈ (0, 1, .…, N _{ h } − 1), where N _{ h } = T _{ s }/T _{ c } is the number of chip positions in a frame and the chip duration is T _{ c }, β is the PPM time shift employed when a _{ j } = 1 and there is no shift when a _{ j } = 0, n(t) is the additive white Gaussian noise (AWGN) with zero mean and twosided power spectral density N _{0}/2, and h(t, ϕ) is the channel realization from the triple SalehValenzuela model in the IEEE 802.15.3c standard.
The transmitted pulse sequence in (1) is:
where p(t) is the Gaussian pulse and f _{ c } is the carrier frequency.
The channel realization can be given by:
where δ(⋅) denotes the Dirac delta function, L denotes the cluster number, and K _{ l } denotes the rays number in thelth lth cluster. The scalar α _{ k, l } denotes the complex amplitude, α _{ k, l } and τ _{ k, l } denote the TOA, and ω _{ k, l } denotes the azimuth of the kth ray in the lth cluster. The average TOA and the average AOA of the lth cluster can be expressed as T _{ l } and Θ_{ l }, respectively.
The first term in (4) denotes the LOS component with
where
and PL _{ d } denotes the path loss, λ denotes the wavelength, A _{NLOS} denotes the attenuation caused due to the NLOS propagation, μ _{ d } denotes the average range, Γ_{0} denotes the reflection coefficient, and h _{1} and h _{2} denote the heights of the transmitted antenna and the received antenna, respectively. G _{ t1}, G _{ t2}, G _{ r1}, and G _{ r2} denote the gains of the transmitted antenna for the path 1 and path 2 and the gains of the received antenna for the path 1 and path 2, respectively.
The ED integrator output for a received signal r(t) can be expressed as:
where N _{ s } is the number of pulses per symbol, T _{ b } is the integration duration, 3 T _{ s }/2 is the integration interval, N _{ b } = 3 T _{ s }/2 T _{ b } is the number of samples, j is the frame index, and the frame duration is T _{ s }. The TOA estimate can be expressed as:
where η is the threshold which is based on a normalized threshold η _{norm} given by:
As min(z[n]) is typically the integration of only noise, this can have a detrimental effect on the receiver performance due to it is sensitivity to noise. In this paper, an extended threshold crossing (ETC) algorithm is proposed where the threshold η can be expressed as:
The problem is then how to obtain η _{norm}. Both curvefitting and fixed threshold (FT) have been employed, but these methods cannot provide precise TOA estimates because the first arriving signal component cannot be accurately determined. The simplest threshold selection was the maximum energy selection (MES) method, i.e., the maximum energy was considered as the threshold. As shown in Fig. 2, the actual TOA locates before the maximum energy, which is challenging to acquire due to the first arriving path which cannot be detected. As a result, the MES method cannot achieve precise TOA estimations as shown in Fig. 2.
Proposed ED receiver
Figure 3 shows the energy values for a 60GHz MMW signal (1) with the 3c CM1.1 model for an SNR of 10 dB and the signal without noise. These results indicate that many energy values are only noise which will have a significant effect on \( {\widehat{\tau}}_{\mathrm{TOA}} \). To reduce the effects of noise and improve the estimation accuracy, an improved ED receiver based on EMD will be presented in this section.
Figure 4 shows the proposed receiver with signal denoising using EMD. EMD adaptively decomposes the received signal into N IMFs, i.e., IMF_{1}, ..., IMF_{ N }, in descending frequency order and a residual trend mode. The residual trend mode is the difference between the signal and the IMFs. The noise can be estimated using the first few IMFs and the residual trend mode [17, 18].
To obtain the IMFs of a signal using EMD, let r(t) be the received signal, α(t) the residual component, m(t) the shifted component, IMF_{ v }(t) the vth IMF, and r _{ u }(t) (r _{ l }(t)) the upper (lower) envelope using a cubic spline function according to the local extrema.
The EMD algorithm is as follows [13]:

1)
Set v = 0 and α(t) = r(t);

2)
Let m(t) = α(t);

3)
Determine the local extrema (maximum and minimum) of m(t);

4)
Construct the upper and lower envelopes of the extrema using cubic spline functions;

5)
Determine the average envelope given by [14]:

6)
If the average envelope is zero or has only one extrema or zero crossing, go to step 7. Otherwise, calculate [15]:
and go to step 3.

7)
If α(t) is monotonic, stop. Otherwise, set v = v + 1, IMF_{ v }(t) = m(t), and α(t) = r(t) − IMF_{ v }(t) and go to step 2 [16].
A 60GHz MMW signal (1) with the 3c CM1.1 model and an SNR of 10 dB was decomposed using EMD into 11 IMFs and a residual trend mode. The first eight IMFs are shown in Fig. 5.
The boundary between the noise and signal components is determined as:
where
and N _{ v } is the number of values in IMF_{ v }, C _{ v } is the maximum value in IMF_{ v }, and O _{ v } is the number of local extrema in IMF_{ v }.
For a signal in AWGN, (14) is approximately a constant [17, 18]. As a result, IMF_{ v } can be regarded as noise when RP_{ v } < 1 [18]. It was determined by simulation that IMFs 1 to 3 can be considered as noise. The histograms of the first three IMFs of white Gaussian noise and the 60GHz signal (1) with the CM1.1 model and SNR = 10 dB are shown in Fig. 6. IMF_{1} for the 60GHz MMW and AWGN follow the bimodal distribution, and the other IMFs follow the Gaussian distribution. The corresponding histograms are presented in Fig. 7. These results indicate that the first three IMFs have similar distributions in both cases. Table 1 shows the means and standard deviations of the first four IMFs, which confirm that the first three can be considered as noise [18]. Note that when v > 3 the IMFs differ significantly.
The received signal can be approximated as:
And the noise estimated as:
Thus, the noise energy is estimated in part B in Fig. 4 as:
And the signal energy values obtained in part C of Fig. 4 as:
Figure 8 shows the energy values after denoising using EMD with the 3c CM1.1 model and an SNR of 10 dB. The similarity, i.e., the energy values with strong amplitude between these results in Fig. 8 and Fig. 3a, confirms the effectiveness of EMD denoising.
Energy signal characteristics
In this section, the curl of the energy values is considered. The other characteristics including the kurtosis, skewness, maximum slope, and the standard deviation have been discussed in [11] in our previous work. However, the proposed threshold selection method cannot achieve accurate TOA estimations in low SNR values.
Curl
In vector calculus, the curl is an operator that describes the infinitesimal rotation of a vector field. The maximum curl of the received energy values is obtained as follows. Let (U, V) be a 2D vector field with U = V and
The arrays X _{ U } and Y _{ V } are the coordinates of U and V, respectively, and X _{ U } is an N _{ b } × N _{ b } dimensional array given by:
Y _{ V } is also an N _{ b } × N _{ b } dimensional array with Y _{ V } = (X _{ U })^{T}. Let the γth difference of the matrix U in the Y _{ V }(:, γ) direction be ℜ(γ) and the γth difference of the matrix V in the X _{ U }(γ, :) direction be ℑ(γ). The curl is then given by:
The maximum curl is the maximum of C(γ).
Signal parameter characteristics
In this section, the characteristics of the parameters of the energy values are investigated in indoor residential LOS and NLOS scenarios for SNR within 4–30 dB. For each SNR, 1000 channels were generated and the received signals sampled at 10 GHz. The system parameters are T _{ s } = 200ns, T _{ b } = 4ns, T _{ c } = 1ns and \( \widehat{\tau}\sim U\left(0,{T}_f\right) \), where U(0, T _{ f }) denotes the uniform distribution.
Figure 9 shows the skewness and maximum curl of the energy values. This figure also shows the kurtosis, maximum slope, standard deviation, the ratios of the kurtosis to the skewness (K/S), the ratios of the maximum slope to the standard, and the products of the maximum slope and the standard deviation (MSD). These results show that the skewness and kurtosis increase as SNR increases, but the skewness changes more quickly. Conversely, the maximum curl, maximum slope, and standard deviation decrease as SNR increases, but the maximum curl changes more quickly. Moreover, the skewness changes little while the maximum curl changes quickly when SNR < 8 dB. On the other hand, when SNR > 8 dB, the skewness changes rapidly but the maximum curl has little change. As a result, these parameters cannot reflect a wide range of SNRs.
TOA estimation using machine learning
In this section, a new metric was proposed using the skewness and maximum curl of the energy samples. This metric can provide a better indication of the changes in SNR than the parameters considered in the previous section. Then, ELM is used to determine the relationship between this metric and the threshold values.
The new metric
The parameters presented in Section 3 are not sufficiently sensitive to changes over a wide range of SNRs. Thus, a new matric is presented in this section. The mean absolute error (MAE) of the TOA estimates is used to evaluate the TOA algorithms and is given by:
where τ _{ m } and \( {\widehat{\tau}}_m \) are the mth propagation time and TOA estimate, respectively, and N _{ m } is the number of TOA estimates.
The MAE of the TOA estimation presented in [7] is less than that in [6] because the skewness changes more rapidly than the kurtosis when SNR > 12 dB. However, neither of these approaches can provide good performance over a wide range of SNRs. Thus, a new metric is proposed which is given by:
where S is the skewness of the energy values and C is the maximum curl of these values. S _{max} and S _{min} are the maximum and minimum values of the skewness and C _{max} and C _{min} are the maximum and minimum values of the curl, respectively. Based on the extensive simulation results, Q and R should be less than 20 and satisfy Q ≥ 3 R. In this paper, Q = 20 and R = 2.
To test the sensitivity of J with respect to SNR, the average value of J was determined for the CM1.1 and CM1.2 channel models with T _{ b } = 1 and 4 ns. Figure 10 shows that J is a monotonic function of SNR and is more sensitive to SNR than the other parameters over the range of SNRs.
Threshold determination
To find the best normalized threshold η _{best}, the relationship between MAE and η _{norm} with respect to J, the channel, and the integration period are examined. Except for small values of J, the MAE decreases with increasing J and the minimum MAE decreases with increasing J. Figure 11 shows the relationship between the normalized threshold and the MAEs. Usually, the normalized threshold with the minimum MAE was referred as the best normalized threshold. The value of η _{norm} for the minimum MAE is defined as the best normalized threshold η _{best}. As shown in Fig. 12, the average η _{best} is defined as:
where
N _{ i } is the number of integration periods and i is the integration period in nanosecond.
The integration periods considered are 1 and 4 ns, so that
where
ELMassisted TOA estimation
Artificial neural networks (ANN) such as the support vector machine (SVM) and back propagation neural network (BPNN) has been widely employed to solve regression problems [18,19,20,22]. ELM is a new and efficient method to train a singlehiddenlayer feedforward neural networks (SLFNs) [23]. ELM has been applied in many fields such as signal and information processing due to the good general performance and fast learning with minimal human intervention [23,24,25,26,27,28,29,30,33]. ELM has been shown to provide better results than SVM and require less training time than BPNN. The most important property is that the parameters in the hidden layers can be automatically determined and a simple generalized inverse operation can be used to obtain the output weight. Further, the only parameter that needs to be specified is the number of nodes in the hidden layer [24]. Due to the complexity of the indoor channel environment, it is very difficult to estimate TOA accurately using traditional curvefitting methods. ELM can be used to resolve this problem by determining the relationship between the input and output parameters according to the training results.
In ANN and SVM, the parameters in different layers need to be tuned which can consume significant time. Further, the learning rate and learning epochs must be tuned through an iterative procedure which makes it difficult to obtain optimal values [23, 24]. Different from NN and SVM solutions, ELM is a tuningfree algorithm. The input weights and hidden bias can be generated randomly before training the SLFN based on the sample set.
In this paper, ELM is used to determine the relationship between J and η _{opt}. The number of nodes in the hidden layer is determined based on the training accuracy. The root mean square error (RMSE) was calculated for 1000 training trials for each number of nodes. The results obtained show that the RMSE is a decreasing function of the number of nodes, but when this number is 30, the RMSE is approximately constant. As a result, the number of node in the hidden layer is set to 30. The values J are rounded to the nearest integer or half integer, and these rounded values are used in training the ELM. As η _{opt} varies from 0 to 1, the sigmoid function is used as the activation function which is given by:
Results and discussion
In this section, the MAEs of several wellknown TC methods were examined for SNRs in the range 4 to 30 dB with the IEEE 802.15.3c channel models. Figure 13a shows the calculated MAEs of the obtained TOA estimations with the LOS propagation including the channel models CM 1.1 and CM 3.1, while the MAEs of the obtained TOA over the channel models CM 2.1 and CM 4.1 were shown in Fig. 13b.
Based on the obtained results, the following conclusions can be obtained:

(1)
The obtained MAEs of the new proposed ED receiver with respect to each given SNR within 4–30 dB indicate the excellent capability of removing the AWGN and show better performances on TOA estimations compared with several wellknown algorithms based on the ED receiver whether in LOS or NLOS scene for each given integration period as shown in Fig. 14.

(2)
The obtained MAEs of the improved ED received indicate that the ability to achieve TOA estimations in LOS models is better than that in NLOS models with regard to each given SNR within 4–30 dB. The difference is 15 ns at best.

(3)
Integration periods usually can cause bad effects on TOA accuracy; the TOA estimation accuracy increases with the increasing integration period regardless of the propagation environment as shown in Fig. 15. In the majority of cases, the MAEs with T _{ b } = 1 ns is less than that with T _{ b } = 4 ns and the difference is at most 2 ns.

(4)
It was turned out that the new ED receiver was independent of the propagation models and integration periods.
Here, ELM denotes the new ED method with the optimum thresholds acquired with the ELM algorithm, MES denotes the maximum energy selection method, K denotes the kurtosis based method, and FT denotes the fixed normalized threshold. In all cases, the average MAE of the improved ED receiver with the ELM algorithm is the lowest among the above algorithms.
Conclusions
In this paper, an improved energy detector (ED)based TOA estimation method for impulse radio (IR) 60 GHz wireless system was presented using the empirical mode decomposition (EMD). The proposed solution employs the curl, the kurtosis, and the skewness of the energy values with an extreme learning machine (ELM). The optimum thresholds were determined for the IEEE 802.15.3c channel models in both LOS and NLOS environments. Results were presented which show that the proposed ELM approach provides better performance than other EDbased TOA estimation methods.
References
 1.
N Leonor, R Caldeirinha, T Fernandes, et al., A simple model for average reradiation patterns of single trees based on weighted regression at 60 GHz, IEEE T. Antenn Propag 63(11), 5113–5118 (2015)
 2.
T Nishesh, R Thipparaju, A switched beam antenna array with butler matrix network using substrate integrated waveguide technology for 60 GHz wireless communications. Int J Electron Commun. 70, 850–856 (2016)
 3.
T Sakamoto, S Okumura, R Imanishi, et al., Remote heartbeat monitoring from human soles using 60 GHz ultrawideband radar. IEICE Electron Expr 12(21),16 (2015)
 4.
R Hazra, A Tyagi, A survey on various coherent and noncoherent IRUWB receivers. Wireless Pers Commun 79(3), 2339–2369 (2014)
 5.
X Liang, H Zhang, T Lu, et al., Energy detector based TOA estimation for MMW systems using machine learning. Telecommun Syst. 64(2), 417–427 (2017)
 6.
I Guvenc, Z Sahinoglu, Threshold selection for UWB TOA estimation based on kurtosis analysis. IEEE Commun Lett. 9(12), 1025–1027 (2005)
 7.
Guvenc, I., and Sahinoglu Z., Multiscale energy products for TOA estimation in IRUWB systems, in Proc. IEEE Global Telecommun. Conf., St. Louis, 209213, 2005.
 8.
Guvenc, I., and Sahinoglu Z., Thresholdbased TOA estimation for impulse radio UWB systems, in Proc. IEEE Int. Conf. on Ubiquitous Wireless Broadband, Zurich, 420425, 2005.
 9.
A Maali, A Mesloub, M Djeddou, H Mimoun, G Baudoin, A Ouldali, Adaptive CACFAR threshold for noncoherent IRUWB energy detector receivers. IEEE Commun. Lett. 13(12), 959–961 (2009)
 10.
W Feng, T Zhi, BM Sadler, Weighted energy detection for noncoherent ultrawideband receiver design. IEEE Trans Wirel Commun. 10(2), 710–720 (2011)
 11.
X Liang, H Zhang, T Lu, TA Gulliver, Extreme learning machine for 60 GHz millimeter wave positioning. IET Commun. 11(4), 483–489 (2017)
 12.
T Wang et al., An EMDbased filtering algorithm for the fiberoptic SPR sensor. IEEE Photon J 8(3), 18 (2016)
 13.
H Zhang, F Wang, D Jia, et al., Automatic interference term retrieval from spectral domain lowcoherence interferometry using the EEMDEMDbased method. IEEE Photon. J 8(3),19 (2016)
 14.
ZW Xu, T Liu, Vital sign sensing method based on EMD in terahertz band, EURASIP J. Adv Sig Pr 2014(1)(2014) doi:10.1186/16876180201475
 15.
YL Wu, BW Shen, An evaluation of the parallel ensemble empirical mode decomposition method in revealing the role of downscaling processes associated with African easterly waves in tropical cyclone genesis. J Atmos Ocean Tech. 33(8), 1611–1628 (2016)
 16.
Y Kopsinis, S McLaughlin, Development of EMDbased denoising methods inspired by wavelet thresholding. IEEE Trans. Signal Process. 57(4), 1351–1362 (2009)
 17.
Z Wu, NE Huang, A study of the characteristics of white noise using the empirical mode decomposition method. Proc R Soc A 460, 1597–1611 (2004)
 18.
R Chen, B Tang, J Ma, Adaptive denoising method based on ensemble empirical mode decomposition for vibration signal. J Vibration Shock 31(15), 83–86 (2012)
 19.
P Sunita, S Archana, PP Siba, A new training strategy for neural network using shuffled frogleaping algorithm and application to channel equalization. Int. J. Electron. Commun. 68, 1031–1036 (2014)
 20.
G Mahrokh, A Hamidreza, Performance analysis of neural network detectors in DS/CDMA systems. Int. J. Electron. Commun. 57(3), 220–236 (2003)
 21.
K Fabian, The experimental results of the bulkdriven quasifloatinggate MOS transistor. Int. J. Electron. Commun. 69, 462–466 (2015)
 22.
P Sotirios, S Katherine, Mobile radio propagation path loss prediction using artificial neural networks with optimal input information for urban environments. Int. J. Electron. Commun. 69, 1453–1463 (2015)
 23.
G Huang, Q Zhu, C Siew, Extreme learning machine: a new learning scheme of feedforward neural networks, Proc. IEEE Int Joint Conf Neural Netw, 2, 985990 (2005)
 24.
L Shi, B Lu, EEGbased vigilance estimation using extreme learning machines. Neurocomputing 102, 135–143 (2013)
 25.
Y Peng, WL Zheng, BL Lu, An unsupervised discriminative extreme learning machine and its applications to data clustering. Neurocomputing 174(A), 250–264 (2014)
 26.
J Chorowski, J Wang, JM Zurada, Review and performance comparison of SVMand ELMbased classifiers. Neurocomputing 128, 507–516 (2014)
 27.
H Zhong, C Miao, Z Shen, Y Feng, Comparing the learning effectiveness of BP, ELM, IELM, and SVM for corporate credit ratings. Neurocomputing 128, 285–295 (2014)
 28.
XR Lee, CL Chen, HC Chang, Y Lee, A 7.92 Gb/s 437.2 mW stochastic LDPC decoder chip for IEEE 802.15.3c applications. IEEE T CircuitsI 62(2), 507–516 (2015)
 29.
T Lu, Pulse waveforms for 60 GHZ Mary pulse position modulation communication systems. IET Commun. 7(2), 169–179 (2013)
 30.
M Xia, W Lu, J Yang, et al., A hybrid method based on extreme learning machine and knearest neighbor for cloud classification of groundbased visible cloud image. Neurocomputing 160, 238–249 (2015)
 31.
S Ding, H Zhao, Y Zhang, et al., Extreme learning machine: algorithm, theory and applications. Artif. Intell. Rev. 44(1), 103–115 (2015)
 32.
JP Nobrega, ALI Oliveira, Kalman filterbased method for online sequential extreme learning machine for regression problems. Eng. Appl. Artif. Intell. 44, 101–110 (2015)
 33.
W Li, D Wang, T Chai, Burning state recognition of rotary kiln using ELMs with heterogeneous features. Neurocomputing 102,144153 (2013)
Acknowledgements
This work was funded by the Nature Science Foundation of China (41527901), Major Program of China’s Second Generation Satellite Navigation System (GF**********03), Fundamental Research Funds for the Central Universities (201713018), National High Technology Research and Development Program of China (2012AA061403), National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (2014BAK12B00), National Natural Science Foundation of China (61501424), National Natural Science Foundation of China (61701462), Ao Shan Science and Technology Innovation Project of Qingdao National Laboratory for Marine Science and Technology (2017ASKJ01), and the Qingdao Science and Technology Plan (17117jch).
Author information
Affiliations
Contributions
XL conceived and designed the experiments, performed the experiments, analyzed the data, and wrote the paper. HX helped to carry out the experiment and analyzed the data. HZ, TL, and TAG helped to review and revise the whole paper. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Liang, X., Zhang, H., Lyu, T. et al. A novel time of arrival estimation algorithm using an energy detector receiver in MMW systems. EURASIP J. Adv. Signal Process. 2017, 83 (2017). https://doi.org/10.1186/s136340170520x
Received:
Accepted:
Published:
Keywords
 Empirical mode decomposition (EMD)
 Impulse radio (IR)
 Energy detection (ED)
 Millimeter wave (MMW)
 Time of arrival (TOA)