Next Article in Journal
Position Tracking During Human Walking Using an Integrated Wearable Sensing System
Previous Article in Journal
An AST-ELM Method for Eliminating the Influence of Charging Phenomenon on ECT
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Signal Subspace Smoothing Technique for Time Delay Estimation Using MUSIC Algorithm

1
Information Engineering College, Shanghai Maritime University, Shanghai 201306, China
2
Institut d’Electronique et Télécommunications de Rennes, UMR CNRS 6164, Université de Nantes, 44306 Nantes CEDEX 3, France
3
Cerema, Project-Team ENSUM, 49136 Les Ponts de Cé, France
4
School of Electronic and Information Engineering, South China University of Technology, Wushan Road, Tianhe District, Guangzhou 510641, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(12), 2868; https://doi.org/10.3390/s17122868
Submission received: 11 November 2017 / Revised: 5 December 2017 / Accepted: 7 December 2017 / Published: 10 December 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
In civil engineering, Time Delay Estimation (TDE) is one of the most important tasks for the media structure and quality evaluation. In this paper, the MUSIC algorithm is applied to estimate the time delay. In practice, the backscattered echoes are highly correlated (even coherent). In order to apply the MUSIC algorithm, an adaptation of signal subspace smoothing is proposed to decorrelate the correlation between echoes. Unlike the conventional sub-band averaging techniques, we propose to directly use the signal subspace, which can take full advantage of the signal subspace and reduce the influence of noise. Moreover, the proposed method is adapted to deal with any radar pulse shape. The proposed method is tested on both numerical and experimental data. Both results show the effectiveness of the proposed method.

1. Introduction

Time Delay Estimation (TDE) has been a hot issue for many years. It has a great number of applications in radar, sonar, geophysics, medical imaging, communications, and so on. In [1], a unified approach to TDE from low-rate samples of the received signal is proposed for multipath channel estimation. A robust Capon beamformer based on time delay or time reversal is applied for ultrasound imaging in [2]. Moreover, the authors in [3] focus on TDE in both active and passive systems and take the envelope variation into account. In the field of civil engineering, time delays are important parameters for the quantitative interpretation of Ground-Penetrating Radar (GPR) data [4,5]. Within the range of centimeter wavelengths, GPR is often exploited for specific applications of stratified media, like roadways [6] or walls [7]. The structure (the layer thickness) of the media can be extracted from the time delays of the backscattered echoes associated with each interface and the dielectric constants of the various layers [6].
TDE is usually performed using the conventional FFT-based methods (inverse FFT or cross-correlation methods). However, the resolution of such methods is restricted by the frequency bandwidth of GPR. The case of small pavement thicknesses was studied in recent papers [8]. The main difficulty with data processing lies in the detection of close backscattered echoes. Some particular pavement materials are made up of thin layers (thickness 3 cm). The conventional methods are not able to distinguish close backscattered echoes (overlapped echoes). In this case, high resolution methods like MUSIC [9,10,11,12] and ESPRIT [13,14,15] are more suitable for TDE. Unlike the situations in [2,3], where the signals are supposed to be totally uncorrelated, in practice, the backscattered echoes are highly correlated (even coherent). Under this condition, the cross-correlation between the backscattered echoes may be too high to degrade the performance of the high resolution methods, such as the MUSIC algorithm, due to the rank loss of the data covariance matrix [16].
In order to apply high resolution methods, like in [1], for highly correlated echoes, sub-band averaging techniques are required to decorrelate the correlation between echoes. The well-known sub-band averaging technique, Spatial Smoothing Preprocessing (SSP), was firstly proposed in [17] and then developed in [18,19], which divides the whole frequency band into a series of overlapping sub-bands to obtain a new data covariance matrix with restored rank. Moreover, some improvements of SSP have been suggested in [14,15,20,21,22,23]. For example, [21] proposes an improved spatial smoothing technique, which takes full advantage of all the cross sub-band correlation matrices and auto-correlation matrices, while SSP only uses auto-correlation matrices.
However, the above preprocessing methods make use of the information from both the signal and noise subspaces of the data covariance matrix. In fact, it is not necessary to apply sub-band averaging techniques on the noise subspace. Furthermore, in this application, to apply the sub-band averaging technique, a whitening procedure is necessary. Due to the GPR pulse shape, the noise covariance matrix after the whitening procedure is no longer an identity matrix, which still contains the radar pulse. Therefore, in this paper, the MUSIC algorithm combined with an adaptation of the Signal Subspace Smoothing (SSS) technique [24] is proposed for TDE. The proposed method only makes use of the signal subspace and can be applied on any GPR pulse shape. The performance of the proposed method is tested on both numerical and experimental data. Both simulation and experimental results prove the effectiveness of the proposed method.
The rest of this paper is organized as follows: Section 2 gives the received radar data model. In Section 3, we present the proposed sub-band averaging technique. Simulation and experiment results are provided in Section 4. Finally, conclusions are drawn in Section 5.
Notations: ( . ) T , ( . ) * , ( . ) 1 and ( . ) H denote the transpose, conjugate, inverse and conjugate transpose operations, respectively. E ( : ) denotes the ensemble average. diag ( A ) denotes the diagonal elements of matrix A . Vectors and matrices appear in boldface lower case letters and boldface capital letters, respectively.

2. Signal Model

In the roadway survey, we focus on the first layers, which are low-loss media. For pavement materials, the conductivity typically ranges within the interval [ 10 3 ; 10 2 ] S/m, according to [25]. Thus, the media can be considered as a low-loss media. In addition, according to the work in [26], if the medium is slightly lossy, the dispersivity of the medium can be neglected. As a consequence, the echoes are simply time-shifted and attenuated copies of the transmitted signal [4,6,27,28]. Therefore, the received signal model can be written in the time domain as [8,14,22,29]:
r ( t ) = k = 1 d s k e ( t t k ) + n ( t )
For applying spectral analysis techniques to TDE, the received signal is usually formulated in the frequency domain. By using Fourier transform, the received signal model can be expressed as:
r ( f i ) = e ( f i ) k = 1 d s k exp ( j 2 π f i t k ) + n ( f i )
where d is the number of backscattered echoes, which can be composed of multiple reflection echoes. d is either assumed to be known or estimated with some detection criteria [30]; t k is the time delay of the k-th echo; e ( t ) and e ( f i ) are the radar pulse in the time and frequency domains, respectively; s k represents the amplitude of the k-th scattered echo; n ( f i ) is an additive white Gaussian noise, with zero mean and variance σ 2 ; with f i = f 1 + ( i 1 ) Δ f , i = 1 , 2 N the frequency, N the number of used frequencies, f 1 the lowest frequency of the studied frequency band and Δ f the frequency step. Although the model is established in this work, it can also be applied in many applications such as radar, sonar, telecommunication, etc., which require the estimation of time delay or frequency. The received signal model can be written in the following vector form:
r = Λ A s + n
with the following notational definitions:
  • r = r ( f 1 ) r ( f 2 ) r ( f N ) T is the ( N × 1 ) received signal vector, called the observation vector, which may represent either the Fourier transform of the measured GPR signal or the measurements by a step-frequency radar;
  • Λ = diag e ( f 1 ) , e ( f 2 ) , , e ( f N ) is a ( N × N ) diagonal matrix, whose diagonal elements are the Fourier transform of the radar pulse e ( t ) ;
  • A = a ( t 1 ) a ( t 2 ) a ( t d ) is the ( N × d ) mode matrix;
  • a ( t k ) = [ e 2 j π f 1 t k e 2 j π f 2 t k e 2 j π f N t k ] T is the mode vector;
  • s = s 1 s 2 s d T is a ( d × 1 ) vector;
  • n = n ( f 1 ) n ( f 2 ) n ( f N ) T is the ( N × 1 ) Gaussian noise vector with zero mean and covariance matrix σ 2 I .
According to (3) and assuming the noise to be independent of the echoes, the covariance matrix Y of r can be written as:
Y = E r r H = Λ A E s s H A H Λ H + E n n H = Λ A S A H Λ H + σ 2 I
where S is the ( d × d ) -dimensional covariance matrix of s ; I is the identity matrix. To apply the sub-band averaging techniques, the influence of radar pulse must be removed from the radar signal. Therefore, in the following, the data are divided by the pulse, then the new observation vector r can be written as r = Λ 1 r = A s + Λ 1 n = A s + b , with b the new noise vector. The new covariance matrix R can be written as:
R = E ( r r H ) = Λ 1 Y Λ H = A S A H + σ 2 Σ
with:
Σ = Λ 1 Λ H = diag ( 1 | e ( f 1 ) | 2 , 1 | e ( f 2 ) | 2 , , 1 | e ( f N ) | 2 )

3. Sub-Band Averaging Technique

In this section, we start with a short review of the conventional SSP and MSSP (Modified Spatial Smoothing Preprocessing) [22]. Then, the proposed SSS is presented.

3.1. Conventional SSP and MSSP

To mitigate the influence of the cross-correlation between the backscattered echoes, the conventional SSP/MSSP is applied on the data covariance matrix R . The whole frequency band with N sample points is partitioned into M overlapping sub-bands. Each sub-band is composed of L frequency points. Therefore, the maximum number of echoes that can be estimated is L 1 . N, M and L are related to each other by N = L + M 1 .
Let r l denote the ( L × 1 ) data vector on the l-th sub-band. It can be written as: r l = A 1 D l 1 s + b l , where b l is the ( L × 1 ) noise vector on the l-th sub-band; A 1 denotes the first ( L × d ) sub-matrix of A ( A 1 is independent of l); and D denotes the ( d × d ) diagonal matrix expressed as:
D = diag e 2 j π Δ f t 1 , , e 2 j π Δ f t d .
Therefore, the l-th sub-matrix of the covariance matrix R can be written as follows:
R l = A 1 D l 1 S ( D l 1 ) H A 1 H + σ 2 Σ l
where Σ l is the l-th sub-matrix of the noise matrix Σ .
According to SSP [22], the rank restored covariance matrix R SSP can be expressed as follows:
R SSP = 1 M l = 1 M R l .
Likewise, in MSSP [22], the modified covariance matrix R MSSP can be written as follows:
R MSSP = 1 2 M l = 1 M { R l + J R l * J }
where J is the ( L × L ) exchange matrix. From (9) and (10), the conventional sub-band averaging techniques are directly applied on the data covariance matrix.

3.2. Proposed SSS Technique

In this paper, the SSS technique is introduced, which is applied only on the signal subspace. It suffers the rank loss, but contains all the information of the backscattered echoes. Therefore, in the first step, we need to apply the EigenValue Decomposition (EVD) on the covariance matrix Y :
Y = U s Λ s U s H + U n Λ n U n H
where Λ s is a ( d × d ) diagonal matrix containing larger eigenvalues. For totally correlated echoes, the number of larger eigenvalues d is 1, d d . The corresponding eigenvectors are in the ( N × d ) matrix U s . Λ n is a ( ( N d ) × ( N d ) ) diagonal matrix, which contains the other smaller eigenvalues, and the associated eigenvectors are in the ( N × ( N d ) ) matrix U n . According to [9], there exists a ( d × d ) -dimensional full rank matrix T , with U s = Λ AT . The signal subspace can be expressed as:
Y s = U s Λ s U s H = Λ A T Λ s T H A H Λ H .
Dividing matrix Y s by the radar pulse, we have the new signal subspace matrix R s as:
R s = Λ 1 Y s Λ H = A T Λ s T H A H .
The l-th ( L × L ) sub-matrix of R s is written as R s l . Then, we apply the principle of the MSSP technique on R s . Accordingly, the rank restored covariance matrix R can be written as the average of M overlapping sub-matrices as follows:
R = 1 2 M l = 1 M { R s l + J R s l * J }
After decorrelation, subspace methods, like the MUSIC algorithm [31], can be applied for TDE. In addition, the amplitudes of echoes can be calculated by the least squares method [5], and then, the relative permittivity of each layer can be estimated. Finally, the layer thickness can be estimated using the estimated time delay and relative permittivity. Compared with conventional approaches, the proposed method makes full use of signal subspace and is more robust to the noise impact. Moreover, the proposed method is able to be applied on any shape of radar pulse.

4. Simulations and Experiment

4.1. Simulation Results

In this section, the performance of the proposed method is tested on simulated data obtained from (3). In the simulation, the media are composed of three homogeneous layers. The pavement configuration is shown in Figure 1. Assume the observation is performed at nadir and far-field. The parameters of the media are chosen as follows: the relative permittivities of the three layers are ε r 2 = 4 , ε r 3 = 8 and ε r 4 = 9 , respectively. The thicknesses of Layer 1 and Layer 2 are H 1 and H 2 , respectively. The thickness of Layer 3 is assumed as infinite. The frequency band used is [0.5–2.5] GHz, with a 0.05 -GHz frequency step (41 frequency samples). The number of sub-bands (M) is equal to 20. The data covariance matrix is estimated from 500 independent snapshots. The Signal-to-Noise Ratio (SNR) is defined as the ratio between the power of the last primary echo and the noise variance.
In the first simulation, the pseudo-spectrum of the proposed method is calculated and compared with that of MUSIC-MSSP. SNR is fixed at 10 dB. Two cases are studied with different thicknesses:
  • Case a. H 1 = 30 mm and H 2 = 60 mm: The received signal is made up of three primary echoes and one multiple echo (one multiple reflection inside the first layer is considered). The corresponding time delays ( t 1 , t 2 , t M 1 , t 3 ) are 1.00 ns, 1.40 ns, 1.80 ns and 2.53 ns, respectively. The first three echoes are overlapped.
  • Case b. H 1 = 20 mm and H 2 = 10 mm: The received signal is made up of three primary echoes. The corresponding time delays ( t 1 , t 2 , t 3 ) are 1.00 ns, 1.27 ns and 1.46 ns, respectively. The three echoes are overlapped.
The pseudo-spectrums of the proposed method and MUSIC-MSSP are presented in Figure 2 and Figure 3. The time delays of the backscattered echoes are well estimated in Cases a and b by the proposed method. The peaks of MUSIC-SSS correspond to the true time delays, thanks to its high resolution and robustness to the noise impact. However, in Case a, MUSIC-MSSP fails to estimate the weak multiple reflection echo; in Case b, MUSIC-MSSP cannot estimate the third primary echo, as the third echo is too close to the second echo.
In the second simulation, the performance of the proposed method versus SNR is evaluated by a Monte Carlo process with 200 independent runs. The Root-Mean-Squared Error (RMSE) of the studied parameter is defined as:
R M S E = 1 U u = 1 U z ^ u z 2
where z ^ u denotes the estimated parameter for the u-th run of the algorithm and z can be the true value of the k-th time delay ( t k ) or layer thickness ( H 1 or H 2 ); U is the number of independent runs. SNR varies from 0–30 dB. Only Case a is studied. The proposed method is compared with MUSIC-MSSP.
Figure 4, Figure 5, Figure 6 and Figure 7 show the RMSE on the estimated time delays by the proposed method and MUSIC-MSSP. The RMSE on the estimated time delays continuously decreases with increasing SNR for both compared methods. The proposed method has a more significant decrease of the RMSE than that of MUSIC-MSSP; MUSIC-MSSP fails to estimate the multiple echo at low SNR, as shown in Figure 6. Moreover, the RMSE also depends on the amplitude of the backscattered echo. The echo with the larger amplitude has smaller RMSE on TDE. Besides, the RMSE on t ^ 3 is smaller than that of the other three time delays because this echo is not overlapped with the others. It is clear that the RMSE obtained by the proposed method is smaller than that obtained by MUSIC-MSSP for all the considered values of SNR, which means that the proposed method has better accuracy and outperforms MUSIC-MSSP. Furthermore, Figure 8 shows the RMSE on the estimated thicknesses of the first and second layers by the proposed method. It estimates the layer thickness with small RMSE and gives relatively good performance in thickness estimation.
Moreover, the performance of the proposed method versus the number of snapshots is tested. The number of Monte Carlo processes is equal to 200. The parameter setting in Case a is applied here; SNR is fixed at 10 dB. Figure 9, Figure 10, Figure 11 and Figure 12 show the RMSE of the proposed method and MUSIC-MSSP as a function of the number of snapshots. Similar to the second simulation, the RMSE on the estimated time delays using the proposed method continuously decreases as the number of snapshots increases. In addition, it can be seen in Figure 11 that MUSIC-MSSP fails to estimate the multiple echo for each number of snapshots. Finally, it can be concluded that the proposed method has better accuracy than MUSIC-MSSP.

4.2. Experimental Results

In addition, the proposed method is tested on experimental data. A mono-static step frequency radar is used, which is composed of a Vector Network Analyzer (VNA) and a bistatic antenna device whose Transmitter (Tx) and Receiver (Rx) are close to each other. The antennas are about 70 cm above the tested media, which allows them to be in the far-field condition. As shown in Figure 13, the tested media are comprised of a PolyVinyl Chloride (PVC) slab set on a metal plane. The thickness of the PVC is about 4 cm, and its relative permittivity is ε r = 2.97 + 0.0015 j . The radar frequency bandwidth ranges from 1.0 GHz– 2.6 GHz, with a 0.02 -GHz frequency step (81 frequency samples). The radar pulse is measured with a metal plane [8]. In the experiment, the proposed method is tested with the number of sub-bands equal to nine. Figure 14 shows the experimental result of the proposed method. The estimated time delay between two echoes is 0.47 ns, and the estimated thickness of the PVC is approximately 4.09 cm (the relative error is 2.25 % ).

5. Conclusions

In this paper, we propose a MUSIC algorithm combined with the proposed SSS technique for TDE. Compared with conventional methods, the proposed method has three major merits: firstly, it takes full advantage of the signal subspace; secondly, it is robust to the influence of noise; thirdly, it can be applied to any radar pulse. The performance of the proposed method is compared with the conventional MSSP. Both numerical and experimental results show the stability and robustness of the proposed method in TDE.

Acknowledgments

This work is supported in part by China Scholarship Council (No. 201506040050), the National Natural Science Foundation of China (No. 61401157, 61501192) and Guangzhou Science & Technology Plan.

Author Contributions

The general idea was proposed by Meng Sun, Yide Wang and Yuehua Ding. Meng Sun designed and performed the simulations, analyzed the results and wrote the paper. All authors participated in amending the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TDETime Delay Estimation
GPRGround-Penetrating Radar
FFTFast Fourier Transform
MUSICMultiple Signal Classification
ESPRITEstimation of Signal Parameters via Rational Invariance Technique
SSPSpatial Smoothing Preprocessing
MSSPModified Spatial Smoothing Preprocessing
SSSSignal Subspace Smoothing
EVDEigenValue Decomposition
SNRSignal-to-Noise Ratio
RMSERoot Mean Square Error
VNAVector Network Analyzer
PVCPolyVinyl Chloride

References

  1. Gedalyahu, K.; Eldar, Y.C. Time-delay estimation from low-rate samples: A union of subspaces approach. IEEE Trans. Signal Process. 2010, 58, 3017–3031. [Google Scholar] [CrossRef]
  2. Wang, Z.; Li, J.; Wu, R. Time delay and time reversal based robust capon beamformers for ultrasound imaging. IEEE Trans. Med. Imaging 2005, 24, 1308–1322. [Google Scholar] [CrossRef] [PubMed]
  3. Ge, F.; Shen, D.; Peng, Y.; Li, V.O.K. Super-resolution time delay estimation in multipath environments. IEEE Trans. Circuits Syst. I Regul. Pap. 2007, 54, 1977–1986. [Google Scholar] [CrossRef]
  4. Spagnolini, U. Permittivity measurements of multilayered media with monostatic pulse radar. IEEE Trans. Geosci. Remote Sens. 1997, 35, 454–463. [Google Scholar] [CrossRef]
  5. Al-Qadi, I.; Lahouar, S. Measuring layer thicknesses with GPR–theory to practice. Constr. Build. Mater. 2005, 19, 763–772. [Google Scholar] [CrossRef]
  6. Spagnolini, U.; Rampa, V. Multitarget detection/tracking for monostatic ground penetrating radar: Application to pavement profiling. IEEE Trans. Geosci. Remote Sens. 1999, 37, 383–394. [Google Scholar] [CrossRef]
  7. Protiva, P.; Mrkvica, J.; Macháč, J. Estimation of wall parameters from time-delay-only through-wall radar measurements. IEEE Trans. Antennas Propag. 2011, 59, 4268–4278. [Google Scholar] [CrossRef]
  8. Le Bastard, C.; Baltazart, V.; Wang, Y.; Saillard, J. Thin-pavement thickness estimation using gpr with high-resolution and superresolution methods. IEEE Trans. Geosci. Remote Sens. 2007, 45, 2511–2519. [Google Scholar] [CrossRef]
  9. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  10. Devaney, A.J.; Marengo, E.A.; Gruber, F.K. Time-reversal-based imaging and inverse scattering of multiply scattering point targets. J. Acoust. Soc. Am. 2005, 118, 3129–3138. [Google Scholar] [CrossRef]
  11. Ciuonzo, D.; Romano, G.; Solimene, R. Performance analysis of time-reversal MUSIC. IEEE Trans. Signal Process. 2015, 63, 2650–2662. [Google Scholar] [CrossRef]
  12. Ciuonzo, D.; Salvo Rossi, P. Noncolocated time-reversal MUSIC: High-SNR distribution of null spectrum. IEEE Signal Process. Lett. 2017, 24, 397–401. [Google Scholar] [CrossRef]
  13. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  14. Qu, L.; Sun, Q.; Yang, T.; Zhang, L.; Sun, Y. Time-delay estimation for ground penetrating radar using ESPRIT with improved spatial smoothing technique. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1315–1319. [Google Scholar]
  15. Sun, M.; Le Bastard, C.; Wang, Y.; Pinel, N. Time-delay estimation using ESPRIT with extended improved spatial smoothing techniques for radar signals. IEEE Geosci. Remote Sens. Lett. 2016, 13, 73–77. [Google Scholar] [CrossRef]
  16. Pan, J.; Wang, Y.; Le Bastard, C.; Wang, T. DOA Finding with Support Vector Regression Based Forward-Backward Linear Prediction. Sensors 2017, 17, 1225. [Google Scholar] [CrossRef] [PubMed]
  17. Evans, J.E.; Sun, D.F.; Johnson, J.R. Application of Advanced Signal Processing Techniques to Angle of Arrival Estimation in ATC Navigation and Surveillance Systems; Massachusetts Institute of Technology Lexington Lincoln Laboratory: Lexington, MA, USA, 1982. [Google Scholar]
  18. Shan, T.J.; Wax, M.; Kailath, T. On spatial smoothing for direction-of-arrival estimation of coherent signals. IEEE Trans. Acoust. Speech Signal Process. 1985, 33, 806–811. [Google Scholar] [CrossRef]
  19. Pillai, S.U.; Kwon, B.H. Forward/backward spatial smoothing techniques for coherent signal identification. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 8–15. [Google Scholar] [CrossRef]
  20. Williams, R.T.; Prasad, S.; Mahalanabis, A.K.; Sibul, L.H. An improved spatial smoothing technique for bearing estimation in a multipath environment. IEEE Trans. Acoust. Speech Signal Process. 1988, 36, 425–432. [Google Scholar] [CrossRef]
  21. Du, W.; Kirlin, R.L. Improved spatial smoothing techniques for DOA estimation of coherent signals. IEEE Trans. Signal Process. 1991, 39, 1208–1210. [Google Scholar] [CrossRef]
  22. Yamada, H.; Ohmiya, M.; Ogawa, Y.; Itoh, K. Superresolution techniques for time-domain measurements with a network analyzer. IEEE Trans. Antennas Propag. 1991, 39, 177–183. [Google Scholar] [CrossRef]
  23. Weiss, A.J.; Friedlander, B. Performance analysis of spatial smoothing with interpolated arrays. IEEE Trans. Signal Process. 1993, 41, 1881–1892. [Google Scholar] [CrossRef]
  24. Grenier, D.; Bossé, E. A new spatial smoothing scheme for direction-of-arrivals of correlated sources. Signal Process. 1996, 54, 153–160. [Google Scholar] [CrossRef]
  25. Fauchard, C. Utilisation de Radars TréS Hautes Fréquences: Application á L’auscultation non Destructive des Chaussées. Ph.D. Thesis, University of Nantes, Nantes, France, 2001. [Google Scholar]
  26. Daniels, D.J. Ground Penetrating Radar, 1st ed.; IET: London, UK, 2004. [Google Scholar]
  27. Wu, R.; Li, X.; Li, J. Continuous pavement profiling with ground-penetrating radar. IEE Proc. Radar Sonar Navig. 2002, 149, 183–193. [Google Scholar] [CrossRef]
  28. Li, J.; Wu, R. An efficient algorithm for time delay estimation. IEEE Trans. Signal Process. 1998, 46, 2231–2235. [Google Scholar]
  29. Le Bastard, C.; Wang, Y.; Baltazart, V.; Derobert, X. Time delay and permittivity estimation by ground-penetrating radar with support vector regression. IEEE Geosci. Remote Sens. Lett. 2014, 11, 873–877. [Google Scholar] [CrossRef]
  30. Radoi, E.; Quinquis, A. A new method for estimating the number of harmonic components in noise with application in high resolution radar. EURASIP J. Appl. Signal Process. 2004, 2004, 1177–1188. [Google Scholar] [CrossRef]
  31. Stoica, P.; Moses, R.L. Introduction to Spectral Analysis, 1st ed.; Prentice Hall, Inc.: Upper Saddle River, NJ, USA, 1997. [Google Scholar]
Figure 1. Pavement configuration; t k represents the time delay of the k-th interface, and t M 1 represents the time delay of a multiple echo.
Figure 1. Pavement configuration; t k represents the time delay of the k-th interface, and t M 1 represents the time delay of a multiple echo.
Sensors 17 02868 g001
Figure 2. Case a, pseudo-spectrum of MUSIC-SSS and MUSIC-MSSP for TDE with SNR = 10 dB.
Figure 2. Case a, pseudo-spectrum of MUSIC-SSS and MUSIC-MSSP for TDE with SNR = 10 dB.
Sensors 17 02868 g002
Figure 3. Case b, pseudo-spectrum of MUSIC-SSS and MUSIC-MSSP for TDE with SNR = 10 dB.
Figure 3. Case b, pseudo-spectrum of MUSIC-SSS and MUSIC-MSSP for TDE with SNR = 10 dB.
Sensors 17 02868 g003
Figure 4. RMSE on the estimated time delay t ^ 1 versus SNR.
Figure 4. RMSE on the estimated time delay t ^ 1 versus SNR.
Sensors 17 02868 g004
Figure 5. RMSE on the estimated time delay t ^ 2 versus SNR.
Figure 5. RMSE on the estimated time delay t ^ 2 versus SNR.
Sensors 17 02868 g005
Figure 6. RMSE on the estimated time delay t ^ M 1 versus SNR.
Figure 6. RMSE on the estimated time delay t ^ M 1 versus SNR.
Sensors 17 02868 g006
Figure 7. RMSE on the estimated time delay t ^ 3 versus SNR.
Figure 7. RMSE on the estimated time delay t ^ 3 versus SNR.
Sensors 17 02868 g007
Figure 8. RMSE on the estimated thicknesses of the first and second layers: H 1 and H 2 versus SNR.
Figure 8. RMSE on the estimated thicknesses of the first and second layers: H 1 and H 2 versus SNR.
Sensors 17 02868 g008
Figure 9. RMSE on the estimated time delay t ^ 1 versus the number of snapshots.
Figure 9. RMSE on the estimated time delay t ^ 1 versus the number of snapshots.
Sensors 17 02868 g009
Figure 10. RMSE on the estimated time delay t ^ 2 versus the number of snapshots.
Figure 10. RMSE on the estimated time delay t ^ 2 versus the number of snapshots.
Sensors 17 02868 g010
Figure 11. RMSE on the estimated time delay t ^ M 1 versus the number of snapshots.
Figure 11. RMSE on the estimated time delay t ^ M 1 versus the number of snapshots.
Sensors 17 02868 g011
Figure 12. RMSE on the estimated time delay t ^ 3 versus the number of snapshots.
Figure 12. RMSE on the estimated time delay t ^ 3 versus the number of snapshots.
Sensors 17 02868 g012
Figure 13. Experimental setup.
Figure 13. Experimental setup.
Sensors 17 02868 g013
Figure 14. Pseudo-spectrum of MUSIC using the SSS technique with experimental data.
Figure 14. Pseudo-spectrum of MUSIC using the SSS technique with experimental data.
Sensors 17 02868 g014

Share and Cite

MDPI and ACS Style

Sun, M.; Wang, Y.; Le Bastard, C.; Pan, J.; Ding, Y. Signal Subspace Smoothing Technique for Time Delay Estimation Using MUSIC Algorithm. Sensors 2017, 17, 2868. https://doi.org/10.3390/s17122868

AMA Style

Sun M, Wang Y, Le Bastard C, Pan J, Ding Y. Signal Subspace Smoothing Technique for Time Delay Estimation Using MUSIC Algorithm. Sensors. 2017; 17(12):2868. https://doi.org/10.3390/s17122868

Chicago/Turabian Style

Sun, Meng, Yide Wang, Cédric Le Bastard, Jingjing Pan, and Yuehua Ding. 2017. "Signal Subspace Smoothing Technique for Time Delay Estimation Using MUSIC Algorithm" Sensors 17, no. 12: 2868. https://doi.org/10.3390/s17122868

APA Style

Sun, M., Wang, Y., Le Bastard, C., Pan, J., & Ding, Y. (2017). Signal Subspace Smoothing Technique for Time Delay Estimation Using MUSIC Algorithm. Sensors, 17(12), 2868. https://doi.org/10.3390/s17122868

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