Improving the Robustness of Time Difference of Arrival Estimation Based on the Energy Center of Gravity Rearrangement

An accurate estimation of the time difference of arrival (TDOA) is crucial in localization, communication, and navigation. However, a low signal-to-noise ratio (SNR) can decrease the reliability of the TDOA estimation result. Therefore, this study aims to improve the performance of the TDOA estimation of dual-channel sensors for single-sound sources in low-SNR environments. This study introduces the theory of time rearrangement synchrosqueezing transform (TRST) into the time difference of arrival estimation. While the background noise TF points show random time delays, the signal time-frequency (TF) points originating from uniform directions that exhibit identical lags are considered in this study. In addition, the time difference rearrangement synchrosqueezing transform (TDST) algorithm is developed to separate the signal from the background noise by exploiting its distinct time delay characteristics. The implementation process of the proposed algorithm includes four main steps. First, a rough estimation of the time delay is performed by calculating the partial derivative of the short-time cross-power spectrum. Second, a rearrangement operation is conducted to separate the TF points of the signal and noise. Third, the TF points on both sides of the time-delay energy ridge are extracted. Finally, a refined TDOA estimation is realized by applying the inverse Fourier transformation on the extracted TF points. Furthermore, a second-order-based time difference reassigned synchrosqueezing transform algorithm is proposed to improve the robustness of the TDOA estimation by enhancing the TF energy aggregation. The proposed algorithms are verified by simulations and experiments. The results show that the proposed algorithms are more robust and accurate than the existing algorithms.


Introduction
Time delay estimation (TDE) is crucial in many fields, including radar, navigation, sonar, seismic monitoring, acoustic source localization, environmental monitoring, and space exploration applications [1][2][3][4][5][6].Accurate TDE is crucial for source localization systems, such as wireless communication and instantaneous finding systems, which track the relative location of a goal [7][8][9][10].There are two types of TDE distinguished by their active and passive detection methods, namely the time of arrival (TOA) estimation and the time difference of arrival (TDOA) estimation.The former estimates the lag between a signal's emission and its echo's reception.The TOA estimation has mainly been used in active systems, such as active sonars and wireless navigation systems.However, the TDOA estimation determines the variation in the wavefront travel time between distant receiving sensors in an array.Typically, it has been used in passive systems, such as passive radars and microphone arrays.Still, there are numerous challenges in TDOA estimation that must be addressed to meet the requirements for high accuracy and robustness of different applications [11][12][13][14].In recent years, many TDOA algorithms have been proposed, each with strengths and weaknesses.The cross-correlation (CC) algorithm is a prominent TDOA technique, assuming that signals satisfy the stationarity criteria.This approach has been highly favored due to its efficient computational capabilities.The CC algorithm for the TDOA estimation requires statistical knowledge about the independence of the signal and noise.Despite having a low SNR, the CC algorithm still provides accurate estimates.However, suppose the instantaneous correlation energy of the effective signal is lower than that of the noise.In that case, the CC approach will incorrectly estimate the delay since it will attribute the delay to the noise instead of the effective signal [15].
To enhance the TDOA estimation efficiency in low-SNR environments, recent research has proposed the multi-channel cross-correlation coefficient (MCCC) technique.This technique uses the spatial redundancy of various sensors to decrease the noise's influence on the TDE result.Compared to the CC technique, the MCCC technique is much more resilient to noise and can achieve significantly higher estimation accuracy under the same SNR conditions.However, in a reverberant environment, the MCCC technique's performance can substantially decline [16].In Ref. [17], multi-channel spatiotemporal sparse prediction (MCSTSP) and multi-channel spatiotemporal group sparse prediction (MCSTGSP) were proposed for the TED.By exploiting the sparsity of the spoken signal prediction parameter matrix, these algorithms can develop an optimization approach that incorporates the MCCC and MCSTP to improve the TDE performance.Compared to the MCCC, the MCSPSP and MCSTGSP are much more effective in handling reverberation.In Ref. [18], a frequencysliding generalized cross-correlation (FS-GCC) method was proposed.This method extracts the phases of different frequency bands of the inter-correlation spectrum to construct the CC matrix.In addition, the singular value decomposition (SVD) and weight SVD methods are used to obtain robust TDEs from the cross-correlation matrix.In Ref. [19], an improved algorithm for multiple-signal classification (IP-MUSIC) was proposed.This algorithm uses a smoothing technique to express the CC result in the form of a covariance matrix, and the TDE is obtained by the MUSIC algorithm.Further, a technique for estimating time delays below the Rayleigh distance using the eigenvalue decomposition and a function called the exponential kernel correlation (EKC) was proposed.The non-Gaussian impulsive noise is considered to have less of an influence due to the EKC function.An advanced algorithm that combines the EKC function and the MUSIC algorithm was proposed to improve estimation accuracy in noisy environments.This algorithm projects the signal subspace to the noise subspace.However, it requires knowledge of the number of sources and may have reduced accuracy under the conditions of a low SNR or limited data [20].An improved solution called the Dual-Path Transfer Function (DPTF) model was introduced to measure time delay in reflection environments, and a rapid DPTF algorithm that uses gradient-based techniques to obtain accurate values was developed.However, further improvements are still necessary to enhance the model's robustness, especially when operating under the low-SNR conditions [21].
The aforementioned studies on the TDOA estimation have achieved good estimation results [22][23][24].However, the TDE performance can be further improved in a low-SNR environment.In view of that, this study proposes an innovative TDOA algorithm aiming at non-stationary signals.In addition, the analysis in this paper is founded upon two hypotheses: (1) After rearranging the short-time cross-power spectrum (STCPS) based on the delay characteristics, the TF points with stronger SNR are closer to the actual time-delay ridges; (2) In each narrow frequency band, the noise power changes slowly over time but remains relatively constant overall.
The main contributions of this work are as follows: (1) The theory of time rearrangement synchrosqueezing transform (TRST) is extended to the time difference of arrival estimation, and a time difference rearrangement synchrosqueezing transform algorithm is proposed to enhance the performance of the TDOA estimation by separating the signal from the background noise by exploiting the distinct time delay characteristics; (2) A second-order-based time difference reassigned synchrosqueezing transform algorithm is developed to improve the robustness of the TDOA estimation by enhancing the TF energy aggregation.
The rest of this paper is structured as follows.Section 2 proposes two time difference reassignment transform algorithms.Section 3 introduces the evaluation of time difference reassigned synchrosqueezing transformation in the TDOA estimation.Finally, Section 4 presents the results and highlights the main achievements compared to the existing algorithms.

Proposed Algorithm
This section introduces two innovative TDOA estimation algorithms based on the rearrangement method.Specific details are provided in Sections 2.1 and 2.2.
This paper presents the research on the time delay estimation of a shock wave.This research belongs to the study of firing point localization.The shock wave is an impulsive signal produced by a bullet squeezing the air as it travels through the air at supersonic speeds.The shock wave expands as a cone behind the bullet, having the wavefront propagating outward at the speed of sound.The TDE of the shock waves plays a crucial role in fire position localization systems.The application of the TDOA in the passive detection of a fire position through gunfire is shown in Figure 1.

UAV Shock Wave
Bullet Trajectory Fire position In the TF analysis, the TF points are rearranged from the position (ω, u) computed by the short-time Fourier transform to the position of the center of energy (ω, û(ω, u)) according to the TRST theory [25][26][27].This rearrangement operation improves the timefrequency aggregation and accurately represents the non-stationary signal characteristics.In the far-field source localization model, distances from the sound source to different array elements of the array are different.Assuming that the first element of the array is a reference element, the time difference between the arrival of a signal at the reference element and its arrival at the other elements of the array remains constant if the direction of the signal source is unchanged.However, the time difference between the background noise received by different array elements of the array is random.According to the different nature of the time difference between the signal and noise received by the array elements, the TRST rearrangement theory is used to separate the signal TF points from the noise TF points in the short-time cross-power spectrum (STCPS).Each TF point of the STCPS is reassigned from its computed position (ω, u) to the position of the center of the gradient of the TDOA (ω, τ(ω, u)).After the rearrangement, the TF points of the signal are clustered, whereas the TF points of the noise are still randomly distributed in the time-delay domain.
The signal and noise can be separated by extracting the TF points at the center of gravity of the TDOA.
The workflow of the proposed algorithms is presented in Figure 2, where it can be seen that it consists of three steps.First, TDOA can be estimated only by data from multiple spatially separated sensors.Second, the TF points with a stronger SNR are extracted by the rearrangement methods.Finally, the maximum value of the cross-correlation function is determined to realize a more robust estimate of the TDOA.
The flowchart of the proposed method.

Time Difference Reassigned Synchrosqueezing Transform Algorithm
In this study, the rearrangement method based on the TDOA is the time difference rearrangement synchrosqueezing transform (TDST) algorithm.The implementation process of the TDST algorithm is illustrated in Figure 3. Figure 3 illustrates the process of separating the signal TF points from the noise TF points.Particularly, Figure 3c represents the STCPS, which is obtained from the dot product of the STFT of x 1 (t) and x 2 (t).Further, Figure 3d represents the time delay at each TF point in the STCPS, obtained by the derivative of the STCPS with respect to the frequency variable.The time delays of partial TF points are deviated due to the noise interference, as shown in red.Each TF point in Figure 3c is rearranged along the time-delay axis according to the time-delay in Figure 3d.As shown in Figure 3e, the TF points of the signal are aligned on the same time delay (as shown in blue), and the TF points of the noise are randomly distributed on both sides of the time delay energy ridges.

Noise-Free TDST Model
Assume that in a noise-free environment, the transmitted signal is unknown, and the TDOA can be estimated only by data from multiple spatially separated sensors.It should be noted that the signal received by the first sensor is usually set as a reference signal.
where x 1 (t) and x 2 (t) represent the output signals of sensors #1 and #2, respectively; a 1 and a 2 denote the amplitude attenuation coefficients of the signals received by sensors #1 and #2, respectively; τ 0 is the time difference of the signal arrival at sensor #2 relative to sensor #1.
To analyze the frequency features of a non-stationary signal, the STFT is computed by the method of a sliding window in the frequency domain, which is expressed as follows: where X 1 (ω, u) denotes the STFT of x 1 (t); A 1 (ω, u) and φ 1 (ω, u) are the magnitude and phase of X 1 (ω, u), respectively; the amplitude and phase of X 2 (ω, u) can be obtained from A 2 (ω, u) and φ 2 (ω, u), respectively; (3) and ( 4) can be used to derive the STCPS of the signal.
where (.) * represents the conjugate operation, ξ denotes the frequency shift variable, and G(ω) is the frequency domain form of g(t).
To obtain time information from the phase data, the partial derivative of ( 6) with respect to the frequency variable ω can be used, which is expressed as follows [28] The STFT with the window function dG(ω) dω is denoted as X G 1 (ω, u).The time-delay operator can be represented by where τ(ω, u) is the time delay operator that contains the time delay of all TF points of a signal; since the TF points of the signal have the same time delay, the time delay operator also represents the TDOA estimation.
In a noisy environment, the time delay operator contains the time delay data on the TF points of both the signal and the noise, so the time delay operator cannot represent the TDOA estimation of the signal.However, it is possible to separate the signal's TF points from the noise's TF points by using the time-delay operator.A precise TDOA estimation can be obtained by extracting the signal's TF points.

Noise-Optimized TDST Model
Under the noisy settings, the signals received by the two sensors can be modeled as follows: x x where s(t) is the unknown source signal, and w 1 (t) and w 2 (t) are the additive noise terms assumed to be uncorrelated with s(t).
The STFT can be calculated by where X 1 (ω, u) is the STFT of x 1 (t); A 1 (ω, u) and φ 1 (ω, u) are the amplitude and phase of X 1 (ω, u), respectively; W 1 (ω, u) is the STFT of the n 1 (t).
Next, the STCPS can be obtained by where P stn (ω, u) represents the STCPS, which includes the two sensor signals' STCPS and noise STCPS.
The partial derivative of P stn (ω, u) with respect to the frequency variable ω is expressed as follows The time delay operator, including the noise, can be calculated by where τw (ω, u) denotes the time delay operator, including the noise interference.
Thus, an accurate TDOA estimation cannot be achieved directly from the time delay operator.
The signal's TF points originating from the same directions exhibit identical lags, while the noise's TF points show random time delays.The rearrangement method can separate the signal from the noise and enhance the SNR value.
where Q e (ω, τe (ω, u)) denotes a TF point whose TF energy is greater than the energy threshold at time delays; )dω represents the energies at different time delays; ε is the energy threshold.How to choose the value of ε may be an important issue.The value of ε may depend on the real situation or physical environment.If the value of ε is chosen improperly, numerical simulations may yield inaccurate or unreasonable results.
Further, we extract the TF points on each side of the time-delay energy ridge.A onedimensional frequency distribution is obtained by accumulating the TF points at different time delays along the time-delay axis as follows where Q e (ω) denotes the frequency on both sides of the time-delay energy ridge.The cross-correlation (CC) function constructed from Q e (ω) can be derived as follows where R(τ) represents the CC function.
The time delay corresponding to the peak of the CC function indicates the TDOA estimation, which is given by τ0 = arg max τ R(τ) (19) where τ0 denotes the TDOA estimation.
The accuracy and robustness of the delay estimation can be improved by performing the rearrangement and extraction operations.The discrete form of the algorithm is given in Algorithm 1.

Algorithm 1:
The time difference reassigned synchrosqueezing transform algorithm.

Second Order-Based Time-Delay Reassignment Synchrosqueezing Transform Model
Based on the idea of improving the TF aggregation result to enhance the performance of the TDOA estimation, an S-TDST algorithm is proposed.In the S-TDST method, each value of the TF points is shifted from its original position at (ω, u) to the new position of (ω, τs (ω, u)).The first-order Taylor expansion of phase and second-order time-delay operator are combined to achieve this.The second-order time-delay operator is given by [33,34] where τs (ω, u) denotes the second-order time-delay operator; τ(ω, u) is the time delay operator; D(ω, u) represents the gradient operator [35,36], and it is defined by where ∂ ω ξ(ω, u) represents the frequency operator's partial derivative concerning ω.Further, ∂ ω τ(ω, u) refers to the partial derivative of the time delay operator in relation to ω.
, others According to (24), the second-order time-delay operator τs (ω, u) can accurately describe the delay time estimation result.To separate the signal's TF points from the noise's TF points, a rearranging process can be adopted as follows where Q s (ω, τs n (ω, u)) denotes the rearranged TF spectrum.Equation (25) shows that the TF points are rearranged from position (ω, u) to the position (ω, τs n (ω, u)).To improve the SNR, the TF points on each side of the signal time-delay energy ridge are extracted.
Q s e (ω, τs e (ω, u)) = arg max where Q s e (ω, τs e (ω, u)) denotes the TF point where the TF energy is larger than the energy threshold.
Next, the TF points on each side of the time-delay energy ridge are extracted.A onedimensional frequency distribution is obtained by accumulating the TF points at different time delays along the time-delay axis.
where Q s e (ω) is the cross-power spectrum.The CC function constructed from Q s e (ω) can be derived as follows where Rs (τ) represents the CC function.
The time delay corresponding to the peak of the CC function denotes the TDOA estimation result, which is given by τs 0 = arg max τ Rs (τ) (29) where τs 0 denotes the TDOA estimation result.The TF points of the signal are separated from the TF points of the noise based on their different time-delay properties.Since the TF points of the signal have the same time delay while the time delay of the noise is random, the time-delay energy ridges of the signal are formed during the rearrangement process.The SNR is improved by extracting the TF points on the energy ridge and both of its sides, which enhances the accuracy and robustness of delay estimation.Algorithm 2 outlines the specific steps of the procedure.

Simulation and Experimental Results
The effectiveness of the proposed algorithms was first verified by using the simulation data, and then their robustness and accuracy were verified by using real data.Furthermore, the proposed algorithms were compared with the related algorithms, including the generalized cross-correlation (GCC) transform, singular value decomposition frequency-sliding generalized cross-correlation (GCCsvd), and weight singular value decomposition frequency-sliding generalized cross-correlation (GCCwsvd) [18].The proposed algorithms' evaluation parameters included the Gaussian window with a window length of σ = 1.2 × 10 −4 and a window overlap of 75%, using 512 points for Fast Fourier Transform.The FS-GCC algorithm used a Hanning window with a length of 32 points, a 75% overlap, and a Fast Fourier Transform of 512 points for optimal performance.An accurate estimation of the time delay was defined as follows where τi is the estimated time delay of the ith Monte Carlo simulation, and M is the number of simulation runs, and it was set to 500 in this study.Moreover, the root mean square error (RMSE) was used to evaluate the performance of the time delay estimation of different algorithms, and it was calculated as follows where τ 0 denotes the actual time delay.

Simulation Results and Analysis
In the simulations, it was assumed that the transmitting signal was an impulsive signal.The sampling frequency was set to Fs = 32 kHz (i.e., the sampling interval was T s = 1/Fs), the signal length was 2 ms, and the time delay was τ 0 = −3T s .The attenuation coefficient of the signal reaching the second sensor was 85% of the signal reaching the first sensor.The source signal s(t) was defined in the frequency domain as follows where, S 3 (ω) = e −j(0.6ω−2e 0.2−0.03ω .* sin(0.2ω))( Experiment 1.This experiment was conducted to verify the effectiveness of the proposed algorithms in a noise-free environment.Figure 4 displays the TDE process implemented by the proposed algorithms.
Figure 4 demonstrates the effectiveness of the proposed algorithms.As presented on the left side of Figure 4b,c, a straight line is approximately perpendicular to the time-delay axis, which indicates that the proposed algorithms could rearrange the TF points according to the time-delay property.This straight line is called the time-delay energy ridge.The peaks of the red and blue curves coincide with the black dashed line, as displayed on the right side of Figure 4b,c, which demonstrates that the proposed algorithms could accurately estimate time delay based on the time-delay energy ridges.
Experiment 2. This experiment was performed to verify the effectiveness of the proposed algorithms in a noisy environment.The performance of the proposed methods in a noisy environment is shown in Figure 5.
Figure 5 illustrates the effectiveness and accuracy of the proposed approaches in estimating the time delay in noisy environments.As shown on the right side of Figure 5b,c, the peaks of the red and blue curves coincide with the black dashed line, which indicates that the proposed algorithms could accurately estimate the time delay in noisy environments.Figure 5a displays the impulsive signals' waveforms of the dual channels and STCPS.Figure 5b shows the time-delay energy ridges extracted by the TDST and their TDE spectra.The time-delay energy ridges and TDE spectrum obtained by the S-TDST method are depicted in Figure 5c.Comparing the results presented on the right side of Figures 5a and 4a, it can be concluded that many TF points were disturbed by the noise, but some TF points were still unaffected.This demonstrated why the algorithm can accurately estimate the time delay in a noisy environment.Comparing the results displayed on the left side of Figures 5b and 4b, it can be concluded that the TF points unaffected by noise were observed to be distributed along the time-delay energy ridge.In contrast, the TF points interfered with by the noise were randomly distributed on both sides of the time-delay energy ridge.Figure 5 illustrates that extracting the TF points on the time-delay energy ridges that were not disturbed by the noise could enable TDOA estimation.

Experiments under Real-Life Settings
The experiment was conducted to acquire the shock wave signal, as illustrated in Figure 6.
Figure 6 shows the test equipment for acquiring shock waves in an outdoor far-field environment.In Figure 6, the second black square to the right above the tripod represents an embedded quadrature planar microphone array system used to capture and measure the shock waves.In Figure 6, directly below the tray is the instrumentation, which was used to calibrate the distance and bearing of the firing point.The embedded quadrature planar microphone array system was fixed horizontally on a tripod, and the target azimuth was equivalently changed from small to large by rotating the array system horizontally clockwise, with the center of the array as the coordinate origin.The firing points and targets were at the same height as the microphone array.In these experiments, the target azimuth was increased by a 15-degree step from zero, and 20 independent trials were conducted at each target azimuth value.Further, there was no obstruction between the microphone array and the firing point location, and the distance between them was set to 1000 m.Two array elements were randomly selected from the array, and the shock wave data that were collected by both elements were acquired.The shock wave data were digitized to 16 bits with a sampling rate of Fs = 32 kHz.The proposed algorithm's evaluation parameters include a Gaussian window with a window length of σ = 1.2 × 10 −4 and a window overlap of 75%, using 512 points for Fast Fourier Transform.The FS-GCC algorithm used a Hanning window with a length of 32 points, a 75% overlap, and a Fast Fourier Transform of 512 points for optimal performance.7a and 5a, it can be concluded that the STCPS of the muzzle blast signals was similar to the STCPS of the ideal impulsive signal in a noisy environment.The TF points of the STCPS of muzzle blast were aligned on the time-delay energy ridge, as shown on the left side of Figure 7b and the right side of Figure 7a.
Comparing the results presented on the left side of Figures 7b and 5b, it can be concluded that the time-delay energy ridge of the shock waves was similar to that of the impulsive signals in the noise environment.This proved the effectiveness of the proposed algorithm in the TDE of the shock wave signals.Figure 8 illustrates that the TDE performances of all the algorithms were enhanced under increased SNR value.Figure 8a,b show that all algorithms, except the GCC algorithm, could accurately estimate the time delay.The main lobes and sidelobes in the GCC algorithm were almost of the same height, which caused this algorithm to perform poorly.In Figure 8d, it can be seen that only the S-TDST algorithm could achieve valid delay estimates.Although the TDST algorithm could obtain accurate time delay, its sidelobe was considerably high.Further, the comparison algorithm yielded incorrect delay estimates.Figure 8 demonstrates that the performance of the proposed algorithms outperformed those of the comparison algorithms in the low-SNR environments.
Experiment 5.This experiment was conducted to verify the robustness of the proposed algorithms in TDE with various signal sizes.Figure 9 displays the performance of all the algorithms obtained under varying signal lengths.
Figure 9 illustrates that both the proposed and the comparison algorithms had good adaptability to the data size variations.Figure 9a-d show that all algorithms could accurately estimate the time delay.By comparing the results presented in Figure 9a,d, it can be observed that the proposed S-TDST algorithm had the narrowest main lobe and the lowest sidelobe among all algorithms.This indicated that the robustness of the S-TDST outperformed that of the comparison algorithms.Experiment 6.This experiment was performed to verify the effectiveness of the proposed algorithms under the multi-path conditions.The SRN level was set to 20 dB. Figure 10 shows the average delay estimation spectrum versus different multi-path numbers.Figure 10 shows that the proposed algorithms had better robustness than the comparison algorithms in the multi-path environment.Figure 10a reveals that all algorithms could accurately estimate the time delay in the multi-path-free environment.However, contrastive algorithms exhibit higher sidelobes, which reduce their robustness.Further, Figure 10b shows that the proposed algorithms could accurately estimate the time delay of both direct and multi-path signals, whereas the contrastive algorithms failed to do so.
Experiment 7.This experiment was performed to verify the accuracy and robustness of proposed algorithms in the varying-SNR environments.Figure 11 shows the robustness and accuracy of the proposed and comparison algorithms under varying-SNR conditions.Figure 11 shows that the robustness and accuracy of the proposed algorithms outperformed the comparison algorithms in the varying-SNR environment.In Figure 11a, the TDE robustness values of all the algorithms are displayed for varying SNR levels.The results indicated that the contrastive approach showed certain instability in low-SNR environments despite achieving better estimation stability under high-SNR settings.As SNR decreased from 5 dB to −10 dB, the red and blue curves demonstrated smooth changes, indicating excellent robustness for the proposed algorithms.However, the green curve increased steeply, indicating the poor robustness of the GCC algorithm in low-SNR environments.Figure 11b displays the accuracy performance of all algorithms with the variation in the SNR value.The results demonstrated that the proposed algorithms had good TDE accuracy, particularly when the SNR value was higher than 10 dB, with an accuracy close to 100%.Compared to the proposed algorithm, the GCCwsvd and GCCsvd algorithms had the same estimation accuracy and poorer estimation accuracy in low-SNR environments and the same accuracy in high-SNR environments.Comparing the results presented in Figure 11a,b, it can be observed that the proposed algorithms outperformed the comparison algorithms in terms of robustness and accuracy.
Experiment 8.This experiment was performed to verify the computational efficiency of the proposed algorithms and compare it with those of the comparison algorithms.Table 1 shows that the running time varied for different algorithms.

Figure 1 .
Figure 1.Illustration of the localization process of a shock wave source.

Figure 3 .
Figure 3.The schematic diagram of the TDST algorithm: (a) STFT of x 1 (t); (b) STFT of x 2 (t); (c) shorttime cross-power; (d) the time delay of the TF points; (e) short-time cross-power after rearrangement.

Figure 4 .
Figure 4.The time delay estimation process of the proposed algorithms; the red curve indicates the TDST, and the blue curve denotes the S-TDST.

Figure 5 .
Figure 5.The TDE performance of the proposed algorithms in a noisy environment at SNR = 10 dB.

Figure 6 .
Figure 6.Experimental site.Experiment 3.This experiment was performed to verify the effectiveness of the proposed algorithms in the TDE of shock wave signals.The performance of the proposed methods is presented in Figure7.Figure7demonstrates the effectiveness of the proposed algorithm in the TDE of the shock wave signals.As shown in Figure7b,c, the black dashed line aligns with the highest points of the red and blue curves, which indicates that the proposed algorithms could accurately estimate the time delay.Figure7adisplays the double channel's shock waveform and the STCPS; Figure7bdepicts the STCPS of the rearrangement process and the TDE spectrum obtained by the TDST; Figure7cdisplays the STCPS of the rearrangement process and the TDE spectrum obtained by the S-TDST.Comparing the results presented on the right side of Figures7a and 5a, it can be concluded that the STCPS of the muzzle blast signals was similar to the STCPS of the ideal impulsive signal in a noisy environment.The TF points of the STCPS of muzzle blast were aligned on the time-delay energy ridge, as shown on the left side of Figure7band the right side of Figure7a.Comparing the results presented on the left side of Figures7b and 5b, it can be concluded that the time-delay energy ridge of the shock waves was similar to that of the impulsive signals in the noise environment.This proved the effectiveness of the proposed algorithm in the TDE of the shock wave signals.

Figure 7
demonstrates the effectiveness of the proposed algorithm in the TDE of the shock wave signals.As shown in Figure7b,c, the black dashed line aligns with the highest points of the red and blue curves, which indicates that the proposed algorithms could accurately estimate the time delay.
Figure 7a displays the double channel's shock waveform and the STCPS; Figure 7b depicts the STCPS of the rearrangement process and the TDE spectrum obtained by the TDST; Figure 7c displays the STCPS of the rearrangement process and the TDE spectrum obtained by the S-TDST.Comparing the results presented on the right side of Figures

Figure 7 .
Figure 7.The performances of the proposed algorithms in the TDE of the shock wave signals at SNR = 10 dB.

Experiment 4 .
This experiment was performed to verify the proposed algorithms' robustness in the varying SNR environments.The SRN level was set to 5 dB, 0, −5 dB, and −10 dB, in turn.

Figure 9 .
Figure 9. Performance results of the algorithms obtained for different signal lengths: (a) signal length was 64 sampling points; (b) signal length was 54 sampling points; (c) signal length was 32 sampling points; (d) signal length was 16 sampling points; the black dashed line denotes the real delays; SNR = 0 dB.

Figure 10 .
Figure 10.Performance comparison in the multi-path environment.

Figure 11 .
Figure 11.Performance comparison under different SNR values.