Local maximum synchrosqueezes form scaling-basis chirplet transform

Abstract

In recent years, time-frequency analysis (TFA) methods have received widespread attention and undergone rapid development. However, traditional TFA methods cannot achieve the desired effect when dealing with nonstationary signals. Therefore, this study proposes a new TFA method called the local maximum synchrosqueezing scaling-basis chirplet transform (LMSBCT), which is a further improvement of the scaling-basis chirplet transform (SBCT) with energy rearrangement in frequency and can be viewed as a good combination of SBCT and local maximum synchrosqueezing transform. A better concentration in terms of the time-frequency energy and a more accurate instantaneous frequency trajectory can be achieved using LMSBCT. The time-frequency distribution of strong frequency-modulated signals and multicomponent signals can be handled well, even for signals with close signal frequencies and low signal-to-noise ratios. Numerical simulations and real experiments were conducted to prove the superiority of the proposed method over traditional methods.

Abstract

In recent years, time-frequency analysis (TFA) methods have received widespread attention and undergone rapid development. However, traditional TFA methods cannot achieve the desired effect when dealing with nonstationary signals. Therefore, this study proposes a new TFA method called the local maximum synchrosqueezing scaling-basis chirplet transform (LMSBCT), which is a further improvement of the scaling-basis chirplet transform (SBCT) with energy rearrangement in frequency and can be viewed as a good combination of SBCT and local maximum synchrosqueezing transform. A better concentration in terms of the time-frequency energy and a more accurate instantaneous frequency trajectory can be achieved using LMSBCT. The time-frequency distribution of strong frequency-modulated signals and multicomponent signals can be handled well, even for signals with close signal frequencies and low signal-to-noise ratios. Numerical simulations and real experiments were conducted to prove the superiority of the proposed method over traditional methods.

Citation: Hou Y, Wang L, Luo X, Han X (2022) Local maximum synchrosqueezes form scaling-basis chirplet transform. PLoS ONE 17(11): e0278223. https://doi.org/10.1371/journal.pone.0278223

Editor: Bashar Ibrahim, Gulf University for Science and Technology, GERMANY

Received: August 12, 2022; Accepted: November 12, 2022; Published: November 29, 2022

Copyright: © 2022 Hou et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The bat signal used in this paper is publicly available from the Rice University dataset (https://web.archive.org/web/20160403234536/ http://dsp.rice.edu/software/bat-echolocation-chirp). The vibration signals used in this paper are publicly available from the Case Western Reserve University dataset (https://engineering.case.edu/bearingdatacenter/apparatus-and-procedures).

Funding: The manuscript is funded by:(1)Science and Technology Innovation Project of Colleges and Universities in Shanxi Province(China),the award number is 2020L0301; (2) Fundamental Research Program of Shanxi Province(China), the award number is 20210302124545.

Competing interests: The authors have declared that no competing interests exist.

1. Introduction

In reality, many signals are nonstationary [1, 2], and the most common feature of these signals is that their frequency is a one-dimensional non-constant function with respect to a time variable. The frequency used in a Fourier transform, however, is an average description of the signal’s frequency change over time. Nonstationary signals can no longer be analyzed using the conventional frequency approach [3–5]. To represent this crucial time-frequency characteristic of the change in the frequency of nonstationary signals with time, the concept of instantaneous frequency has been developed [6, 7]. An effective method for processing nonstationary signals is time-frequency analysis (TFA) [8–10]. The result is a time-frequency representation (TFR) in the form of a time-frequency-density function [11], which contains rich signal information, including the distribution pattern of signal energy in the time-frequency plane, instantaneous frequency characteristics, and instantaneous bandwidth. Using only conventional signal analysis techniques makes it challenging to perform accurate analyses because mechanical equipment have numerous excitation sources, complicated generation mechanisms [12], sound transmission paths, and has non-smooth, non-linear characteristics [13]. Consequently, fault detection ideas that were previously guided by signal analyses are now being extended [14].

There are various TFA methods available at present. The concept of TFA methods originated from the Gabor expansion theory proposed by a Hungarian physicist, Gabor, in 1946. The famous linear time-frequency transform method, the short-time Fourier transform (STFT), was developed based on it. Classical TFA methods also include the Wigner-Ville distribution (WVD), wavelet transform (WT), and S-transform, which have allowed multi-resolution analyses of signals to be performed [15, 16]. Owing to the constraints of the Heisenberg-Gaber uncertainty principle, linear TFA methods should strike a balance between time and frequency resolutions. Bilinear time-frequency distribution (TFD) methods are also limited in applications dealing with multi-component analysis signals owing to the interference of the cross terms [17, 18]. The chirplet transform (CT) has been proposed to improve the energy concentration of the time-frequency map. The CT is a new TFA method that can be considered as a generalization of STFT and WT. However, when the frequency of a signal shows a nonlinear variation with time, the resolution of CT is low, and the accuracy of its analysis cannot be guaranteed. Therefore, polynomial chirplet transform (PCT) [19], modified spline-kernelled chirplet transform (MSCT) [20], velocity synchronous linear chirplet transform (VSLCT) [21], general linear chirplet transform (GLCT) [22], and scaling-basis chirplet transform (SBCT) [23] have been developed based on CT. However, the energy concentration of these TFA methods is not satisfactory, and they often exhibit poor noise resistance [9].

To overcome the aforementioned limitations, three researchers, Kodera, Gendrin, and Villedary, used the phase information of an analyzed signal to collect the scattered time-frequency energy in the time-frequency plane, which marked the development of the first time-frequency post-processing techniques. The reassignment method (RM) algorithm with a solid theoretical foundation was proposed as a post-processing TFA method. This method is mainly used to enhance the effect of time-frequency representation; however, it does not support signal reconstruction. Daubechies et al. then proposed the synchrosqueezing wavelet transform (SWT) in 2011, SWT rearranges the time-frequency coefficients using the synchrosqueezing operator, shifts the TFD of the signal at any point in the time-frequency plane to the center of gravity of the energy, and enhances the energy of the instantaneous frequency [24, 25]. This can solve the time-frequency ambiguity problem of traditional TFA methods. The synchroextracting transform (SET) and local maximum synchrosqueezing Transform (LMSST) were also proposed by Yu et al. [26, 27]. Unlike the classical synchrosqueezing transform theory, SET is only concerned with the instantaneous frequencies corresponding to the characteristic components of a signal. The divergent time-frequency energy coefficients are removed by the simultaneous extraction operator, and only the time-frequency ridge coefficients are retained. Tu proposed a horizontal synchrosqueezing transform (HST) that solves the problems of traditional SST by applying a local estimation of the group delay [28]. Zhu proposed a synchroextracting transform based on CT (SECT). SECT shows a better performance than certain advanced TFA methods [29].

With these applications of SST, we realize that it can be considered as a postprocessing technique based on the traditional TFA methods (STFT, WT, and CT). Fundamentally, the effectiveness of these postprocessing methods and the superiority of the processed signals also depend on the original TFA methods. The more accurate the instantaneous frequency obtained by the traditional TFA method, the clearer the separation of multi-component signals; therefore, postprocessing can play an important role in multicomponent signal analysis. However, once these traditional TFA methods have obtained the instantaneous frequency deviation coupled with the interference of noise, the post-processing results may become misleading. Therefore, it is necessary to first guarantee that the instantaneous frequency after TFA matches the actual instantaneous frequency of the signal. The SBCT method reconstructs a new chirplet, and this transform can match each instantaneous frequency slope in a multicomponent signal within the same window length. Although higher energy concentrations can be achieved, the frequency resolution is not clear enough. To obtain a higher frequency resolution, this study extends the synchrosqueezing transform to SBCT and proposes a new TFA method called the local maximum synchrosqueezing scaling-basis chirplet transform (LMSBCT), which can suppress the interference of noise more effectively and has a higher time-frequency aggregation compared with the existing TFA methods, to obtain the TFR.

The remainder of this paper is organized as follows: In Section 2, the SBCT and LMSST methods are introduced. In Section 3, the theory of the LMSBCT method proposed in this study and its algorithm implementation are introduced in detail. In Sections 4 and 5, the superiority of the proposed algorithm proposed is discussed and demonstrated through simulation experiments and real cases. Finally, Section 6 concludes the study.

2. Theoretical principles

2.1 SBCT theory

The expression of CT is given by(1)where, s(u) denotes the Hilbert transform of the signal x(u), h(u) denotes the real even Gaussian window function, and φ is a phase function defined as . The second-order derivative of the phase function yields C: It follows that the rotation angle θ has a constant value.

(2)(3)

At one window length, when the instantaneous frequency trajectory of the signal changes with time, different C values are required to achieve a higher energy resolution, and when the signal has multiple components, different C values are also required to match the frequency trajectory of the signal simultaneously. To overcome the limitations of the traditional CT method for determining the chirp rate, SBCT reconstructs a new phase function.

(4)

The second-order derivative of this function gives the equation concerning θ.

(5)(6)

When the signal is at the moment u = t_c, i.e. , θ can be expressed as:(7)that is, when the signal is a multicomponent signal. Simultaneously, different signal components have different θ values corresponding to the component time-frequency spine. Thus, the time-frequency resolution will be greatly improved, and the energy concentration will be higher.

when the signal is at moment t_c+Δu, i.e. ,(8)that is, when the signal is a strong time-varying signal, different moments will have various θ values to correspond with the time-frequency ridge.

Substituting the new phase function in Eq (4) into Eq (1), the SBCT can be expressed as(9)

2.2 LMSST theory

SST is a TFA postprocessing technique that builds on the obtained TFD and uses the local behavior (phase information) near the time-frequency point to perform frequency rearrangement of the TFD. Its significant contribution is to increase the time-frequency aggregation and time-frequency ridges in more detailed.

For a signal to be measured, x(u) should satisfy f∈L²(R). |G(t,ω)| denotes the spectrogram of the STFT, and SST is calculated as(10)where, δ denotes the Dirac function and ω₀(t,ω) can be obtained using the following equation:(11)

G^h′(t,ω) denotes the spectrum obtained after deriving the window function. This postprocessing operation results in a higher energy aggregation of the instantaneous frequencies and better frequency resolution.

LMSST is an improved algorithm based on SST that redefines a new frequency operator using the following definition rules:(12)where, Δ denotes the discrete frequency interval. When the two signal components reach a frequency of , at which point the window function reaches a maximum value of zero, the frequency operator is again given a new rule:(13)

Therefore, LMSST can be expressed as follows(14)

3. LMSBCT

3.1 Theory

Inspired by the LMSST, this study reassigns new time-frequency coefficients in the frequency direction by further processing the SBCT analysis results. According to Eqs (11) and (12), the instantaneous frequency in Eq (11) should be calculated twice through STFT. One is obtained through the conventional STFT, and the other is obtained through STFT, by deriving the window function. To reduce the computational effort, this study used the frequency operator of the local maximum, which only performs the TFA method once.

Therefore, this study proposes a new TFA method, which is expressed as follows:(15)

The signal formula is(16)where, A(u) represents the instantaneous amplitude and ∫v(u)du represents the instantaneous phase. At this point, the ideal instantaneous frequency is derived from the instantaneous phase as v(u). The Taylor expansion of v(u) can be written as(17)

Substituting Eqs (16) and (17) into Eq (9) yields(18)

The coefficients a₁ and a₂ can be derived based on the following assumptions:(19)(20)

M and N must be set in advance, and different sizes means that the chirp rate is divided into different segments from to . Larger M and N values indicate a higher frequency resolution. As(21)(22)the presetting of m₀ and n₀ reduces the computational load. Thus, choosing β₁ and β₂ is another problem. According to kurtosis theory, a larger kurtosis indicates a higher energy concentration. Kurtosis can be expressed as follows [23, 30].

(23)

Therefore, the selection rules of optimal β₁ and β₂ are as follows.

(24)

SBCT is defined as(25)

Therefore, the new frequency operator can be calculated using the following equation:(26)

3.2 Algorithm implementation

LMSBCT Algorithm

1. Initialization: Input L; M; N; m₀; n₀.

2. Calculate SBCT:

for i = 1 to M

for j = 1 to N

sub-SBCTs(: , : , i, j) ← SBCT(i, j);

end for

end for

Find Output SBTC(f,t_c)

3. Local maximum synchrosequeezing

Calculate ω_m(f,t_c)

Find ω_m(f,t_c) =

for t = 1 to T

for f = 1 to F

η←ω_m(f,t_c)

end for

end for

Output T_LMSBCT(t_c,η)

4. Simulation analysis

In this section, three sets of simulated signals are used to demonstrate the superiority of the proposed TFA method with a good time-frequency aggregation and high time-frequency resolution. The selected comparison algorithms were STFT, SBCT, GLCT, VSLCT, SST, SET, and RM.

4.1 Monocomponent signal

To construct a strongly time-varying monocomponent signal, the simulated signal model is considered as follows:(27)

The sampling frequency was set to 1024 Hz and the sampling time was 1 s. Fig 1 shows the ideal instantaneous frequency of the signal. The frequency of the signal varies with time, and the quantity of instantaneous frequency change decreases as time increases. The signals were processed using various TFA algorithms. Fig 2(A)–2(C) shows the TFDs of the STFT, SBCT, and GLCT algorithms, respectively. The results of these algorithms are energy-dispersive and have a poor frequency resolution, which is not sufficient to satisfy the signal-characteristics analysis. Fig 2(D) and 2(E) shows the processing results of SST and RM, which have higher energy concentrations than the three aforementioned algorithms and can describe the instantaneous frequency of the signal. The algorithm proposed in this paper is a postprocessing procedure for SBCT. The results are shown in Fig 2(F), which effectively improves the concentration of the time-frequency representation and accurately obtains the instantaneous frequency curve.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

TFA results obtained through (a) STFT, (b) SBCT, (c) GLCT, (d) SST, (e) RM, and (f) LMSBCT.

To be able to compare different post-processing results more clearly, Fig 2(D)–2(F) are partially enlarged, and the time between 0.4–0.45 s is selected with ideal frequencies. The red line in Fig 3 represents the ideal instantaneous frequency. The edge of the SST is partially blurred, and the energy aggregation is bad. The result of the RM is slightly blurred. However, a detailed analysis could not be conducted. The improved algorithm proposed in this paper not only has a higher time-frequency resolution but also a better energy focus on the spectrum. It is clearer and more accurate for describing the transient characteristics of the signal.

[Figure omitted. See PDF.]

(a) Magnified plots of the SST, (b) RM, and (c) LMSBCT results.

To objectively evaluate the different TFA methods, this study introduces the concept of Rényi entropy. Entropy is a method of quantitatively evaluating information uncertainty. The more random the signal gets, the greater the uncertainty, and the higher the corresponding entropy value, and vice versa. In the field of TFA, the more concentrated the time-frequency energy distribution, the smaller the uncertainty and the smaller the corresponding entropy value. Therefore, the entropy value can be used to judge the degree of concentration of the time-frequency spectrum energy and evaluate the superiority of the results of different TFA methods [31]. As shown in Table 1, LMSBCT has the smallest Rényi entropy among all methods, which implies the highest time-frequency aggregation. This represents a breakthrough in the field of TFA.

[Figure omitted. See PDF.]

4.2 Multicomponent noise-added signal

To demonstrate the generalizability of the method proposed in this study, the second set of numerically simulated signals consists of multicomponent strong frequency-modulated (FM) signals and noise interference. Under the strong noise interference, the ideal TFA method can effectively identify time-varying features, that is, correctly extract the time-frequency ridges of different components, thus demonstrating its good applicability. The numerical analog signal is specifically represented as(28)(29)(30)(31)

The ideal instantaneous frequency corresponding to S₁, S₂, and S₃ can be expressed as(32)(33)(34)

Fig 4 shows the ideal time-frequency ridge and TFD of the multicomponent signal. To test noise robustness, Gaussian white noise was added to the signal, and the calculated signal-to-noise ratio (SNR) was 3 dB. The setting time was 4 s, and the sampling frequency was 120 Hz. Several mainstream TFA methods were compared with the method proposed in this study, and the comparison results are plotted in Fig 5. In terms of the coarseness of the time-frequency ridges, the time-frequency energy spread of STFT and SBCT is large, and cannot describe the TFD of the signal well. Although the results of the three methods (SST, SET, and RM processing) are improved, the existence of noise interferes with the identification of the characteristics of the signal. As shown in Fig 5(F), the processing result of LMSBCT has a clean background and no noise interference, which indicates that the method proposed in this study can substantially improve the energy aggregation and has good noise robustness. This is highly consistent with the ideal time-frequency spectrum of the multicomponent simulated signal shown in Fig 4(A) and achieves the optimum time and frequency resolutions. Fig 6 illustrates the instantaneous amplitude spectra of different time-frequency analysis methods at the time of 0.5 s. The spectral bandwidths of STFT and SBCT are large, and there is no boundary between the frequency components. In contrast, the energy of the frequency spectrum of RM and LMSBCT is concentrated in a narrow bandwidth. The frequency component IF3 was locally amplified (represented in red). Comparing Fig 6(C) and 6(D), LMSBCT has clear boundaries between each frequency component, more concentrated energy, and high noise immunity.

[Figure omitted. See PDF.]

(a) Ideal instantaneous frequency and (b) ideal TFD.

[Figure omitted. See PDF.]

TFA results obtained through (a) STFT, (b) SBCT, (c) RM, (d) SET, (e) SST, and (f) LMSBCT.

[Figure omitted. See PDF.]

TF slices obtained through (a) STFT, (b) SBCT, (c) RM, and (d) LMSBCT at time t = 0.5 s.

To further compare the noise immunity performance of various methods, in this study, the Rényi entropy of different TFA methods was calculated for an input SNR of 0–30 dB. As shown in Fig 7, the Rényi entropy of each method gradually decreases as the SNR gradually increases. The Rényi entropy of SST, SET, and RM are smaller than that of STFT at any SNR, indicating that postprocessing enhances normal time-frequency methods to achieve better resolution. This proves the superiority of the proposed method. In Fig 7, that the Rényi entropy is not significantly affected by the SNR of the input signal, that is, the method in this study has a stronger noise robustness.

[Figure omitted. See PDF.]

4.3 Signals with close instantaneous frequency trajectories

The above two simulation experiments prove that the algorithm can obtain the results of TFD with good time-frequency aggregation in strong time-varying signals and multicomponent signals. The signals in this section are defined by the following equations:(35)(36)(37)

The sampling frequency was 20Hz and the sampling time was 70s. The results of GLCT in Fig 8(A) show a large amount of background noise. The instantaneous frequency ridges of the signal components are completely submerged, and it is impossible to see the number of components present. Fig 8(B) shows the results of STFT, in which the energy dispersion is serious, and the instantaneous frequency trajectory is indistinguishable. Owing to the poor STFT results, the results of the subsequent processing algorithms, SET and SST, also failed to correspond with the expected results. Fig 8(E) and 8(F) presents the analysis results of VSLCT and SBCT, respectively. In the plots, there is no crossover of transient frequencies, but the energy concentration is insufficient. The transient frequency traces in Fig 9 are clear, and there is no cross-mixing. This indicates that even if the signal components are close to each other, the LMSBCT transient frequency division is relatively accurate. Therefore, the energy divergence and overlapping phenomena are well-resolved. The improved algorithm proposed in this study has a high time-frequency concentration and is valid for signals with close components, outperforming other typical traditional algorithms. Table 2 shows the Rényi entropy values of the different algorithms, and the proposed algorithm has the smallest Rényi entropy.

[Figure omitted. See PDF.]

TFA results obtained through(a) GLCT,(b) STFT,(c) SET,(d) SST,(e) VSLCT, and (f) SBCT.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

5. Experimental analysis

5.1 Bat signal

The classical bat signal was used as the standard library to validate the method. This signal was first used by Rice University to validate a new method proposed by other researchers [27]. Because the bat signal contains an FM signal, a full down signal, and an echo delay signal, it can effectively verify the results in complex environments. Digitized echolocation pulse emitted by large brown bat Eptesicus fuscus. The signal had a sampling point count of 400 and sampling frequency of 140 kHz. It is difficult to accurately understand the nonlinear behavior of bat echolocation in the time-domain signals. Moreover, it is difficult to grasp the time-varying characteristics of the signal using only a one-dimensional time-domain or frequency-domain analysis. Extending the one-dimensional time-frequency domain to a two-dimensional time-frequency domain can produce more necessary information. In Fig 10, the time-frequency spectra obtained based on the STFT, GLCT, and SBCT algorithms have a heavy energy divergence, and the resolution is coarse. Noise interferes with the SET method when the frequency is above 60kHz. Fig 11 shows that the LMSBCT substantially improved the effect of the time-frequency spectrum resolution. The energy dispersion phenomenon was resolved effectively. The energy is gathered at the real instantaneous frequency of the signal. In this study, the results of the five TFA methods were also compared using the Rényi entropy as an evaluation index, as shown in Table 3. The Rényi entropy obtained from the LMSBCT method was the smallest, indicating its optimal performance.

[Figure omitted. See PDF.]

TFA results obtained by (a) STFT,(b) GLCT,(c) SBCT, and (d) SET.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

5.2 Vibration signals from the CWRU dataset

Without loss of generality, the bearing dataset provided by Case Western Reserve University was selected in this study to verify the effectiveness of the proposed algorithm [32]. Experiments were carried out using a 2-horsepower Reliance Electric motor with acceleration data measured near and away from the motor bearings. As shown in Fig 12, the test stand included a 2 hp motor, torque transducer/encoder, dynamometer, and control electronics. Vibration data were collected using accelerometers attached to housing with magnetic bases. It can be observed from Fig 13 that the signal consists of two frequency components. All five TFA algorithms could still roughly identify the instantaneous frequency of the nonstationary signal under the interference of noise, but their noise robustness was different. In Fig 13(A), there is a certain amount of mixing between the two frequency components, which deviates from the true instantaneous frequency of the multicomponent signal. In the interference of noise, a large number of interwoven noise textures were generated in the background of Fig 13(B) and 13(C). Their identified time-frequency ridges are inlaid in the noise textures. Fig 13(D) shows the approximate instantaneous frequency of the vibration signal at each moment. However, there are some interwoven textures in the background and the originally smooth and continuous time-frequency curve is broken and partially distorted. With the appearance of noise, Fig 13(E) clearly and accurately shows the instantaneous frequency of the signal at each moment. It does not produce excessive noise in the background. Based on the above analysis of the experimental results, the following conclusions can be drawn: Among the five TFA algorithms, the algorithm proposed in this paper is superior.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

TFA results obtained by(a)SBCT, (b)SST, (C)RM, (d)SET, (e) LMSBCT.

6. Conclusion

Inspired by LMSST, this study proposed a new TFA method, LMSBCT, based on SBCT, which redistributes the new instantaneous frequency operator by extracting the local maxima of the spectrogram in the frequency direction. This method overcomes the shortcomings of traditional TFA methods, improves the aggregation of signals, and achieves high-precision analysis of instantaneous signal frequencies. Three sets of simulation experiments demonstrate the three advantages of this algorithm. (1) The frequency changes of strong time-varying signals can be analyzed effectively. Compared with other TFA methods, the Rényi entropy of LMSBCT can be reduced to 9.6438. (2) The Rényi entropies of the LMSBCT algorithm were always lower than those of the other methods when the SNR was reduced from 30 dB to 1 dB. This implies that the multicomponent signals can be effectively separated, even at low SNRs. (3) This method can also obtain an elaborate TFR when the instantaneous frequencies of the signals are close to each other. Even if the frequency interval of the signal is less than 1Hz, the Rényi entropy of the LMSBCT is the smallest compared to the other methods, which is 8.0927. In this study, bat and vibration signals from the CWRU were selected to demonstrate the effectiveness of the algorithm. Compared with other methods, the algorithm in this study can describe and characterize the time-varying features accurately and in detail, which is beneficial for the subsequent extraction and analysis of signal features. Both the simulation signal and the actual vibration signal application results show that the method has the advantages of high time-frequency resolution and good energy aggregation.

Citation: Hou Y, Wang L, Luo X, Han X (2022) Local maximum synchrosqueezes form scaling-basis chirplet transform. PLoS ONE 17(11): e0278223. https://doi.org/10.1371/journal.pone.0278223

About the Authors:

Yating Hou

Roles: Methodology, Writing – original draft

Affiliations State Key Laboratory of Dynamic Testing Technology, North University of China, Taiyuan, China, School of Information and Communication Engineering, North University of China, Taiyuan, China

https://orcid.org/0000-0002-5894-0033

Liming Wang

Roles: Validation

E-mail: [email protected]

Affiliations State Key Laboratory of Dynamic Testing Technology, North University of China, Taiyuan, China, School of Information and Communication Engineering, North University of China, Taiyuan, China

https://orcid.org/0000-0002-3128-2747

Xiuli Luo

Roles: Visualization

Affiliations State Key Laboratory of Dynamic Testing Technology, North University of China, Taiyuan, China, School of Information and Communication Engineering, North University of China, Taiyuan, China

Xingcheng Han

Roles: Investigation

Affiliations State Key Laboratory of Dynamic Testing Technology, North University of China, Taiyuan, China, School of Information and Communication Engineering, North University of China, Taiyuan, China

References

1. Chen P, Wang K, Zuo MJ, Wei D. An ameliorated synchroextracting transform based on upgraded local instantaneous frequency approximation. Measurement. 2019;148: 106953.

2. Jablonski A, Dziedziech K. Intelligent spectrogram-A tool for analysis of complex non-stationary signals[J]. Mechanical Systems and Signal Processing, 2022;167: 108554.

3. Zheng J, Pan H, Yang S, Cheng J. Adaptive parameterless empirical wavelet transform based time-frequency analysis method and its application to rotor rubbing fault diagnosis. Signal Process. 2017;130: 305–314.

4. Wang Y, Yang W, Li D, Zhang JQ. A novel time-frequency model, analysis and parameter estimation approach: Towards multiple close and crossed chirp modes. Signal Process. 2022;201: 108692.

5. Stanković L, Brajović M, Daković M, Mandic D. On the decomposition of multichannel nonstationary multicomponent signals. Signal Process. 2020;167: 107261.

6. Boashash B, Khan NA, Ben-Jabeur T. Time-frequency features for pattern recognition using high-resolution TFDs: A tutorial review. Digit Signal Process. 2015;40: 1–30.

7. Boashash B. Estimating and interpreting the instantaneous frequency of a signal. I. Fundamentals. Proc IEEE. 1992;80: 520–538.

8. Boashash B, Ouelha S. Designing high-resolution time-frequency and time-scale distributions for the analysis and classification of non-stationary signals: a tutorial review with a comparison of features performance. Digit Signal Process. 2018;77: 120–152.

9. Ding C, Zhao M, Lin J, Liang K, Jiao J. Kernel ridge regression-based chirplet transform for non-stationary signal analysis and its application in machine fault detection under varying speed conditions. Measurement. 2022;192: 110871.

10. Guo Y, Yang L D. LFM signal optimization time-fractional-frequency analysis:Principles, method and application[J]. Digital Signal Processing, 2022;126: 103505.

11. Feng Z, Yu X, Zhang D, Liang M. Generalized adaptive mode decomposition for nonstationary signal analysis of rotating machinery: Principle and applications. Mech Syst Signal Process. 2020;136: 106530.

12. Schmidt S, Zimroz R, Heyns PS. Enhancing gearbox vibration signals under time-varying operating conditions by combining a whitening procedure and a synchronous processing method. Mech Syst Signal Process. 2021;156: 107668.

13. Feng Z, Liang M, Chu F. Recent advances in time-frequency analysis methods for machinery fault diagnosis: A review with application examples. Mech Syst Signal Process. 2013;38: 165–205.

14. Lucà F, Manzoni S, Cigada A, Frate L. A vibration-based approach for health monitoring of tie-rods under uncertain environmental conditions. Mech Syst Signal Process. 2022;167: 108547.

15. Qian Shie, Chen Dapang. Joint time-frequency analysis. IEEE Signal Process Mag. 1999;16: 52–67.

16. Stankovic L, Stankovic S, Dakovic M. From the STFT to the Wigner Distribution [Lecture Notes]. IEEE Signal Process Mag. 2014;31: 163–174.

17. Al-Sa’d M, Boashash B, Gabbouj M. Design of an Optimal Piece-Wise Spline Wigner-Ville Distribution for TFD Performance Evaluation and Comparison. IEEE Trans Signal Process. 2021;69: 3963–3976.

18. Zhu X, Zhang Z, Li Z, Gao J, Huang X, Wen G. Multiple squeezes from adaptive chirplet transform. Signal Process. 2019;163: 26–40.

19. Peng ZK, Meng G, Chu FL, Lang ZQ, Zhang WM, Yang Y. Polynomial Chirplet Transform With Application to Instantaneous Frequency Estimation. IEEE Trans Instrum Meas. 2011;60: 3222–3229.

20. Ma Y, Lv Y, Yuan R, Ge M. Synchro spline-kernelled chirplet extracting transform: A useful tool for characterizing time-varying features under noisy environments and applications to bearing fault diagnosis. Measurement. 2021;181: 109574.

21. Guan Y, Liang M, Necsulescu D-S. Velocity Synchronous Linear Chirplet Transform. IEEE Trans Ind Electron. 2019;66: 6270–6280.

22. Yu G, Zhou Y. General linear chirplet transform. Mech Syst Signal Process. 2016;70–71: 958–973.

23. Li M, Wang T, Chu F, Han Q, Qin Z, Zuo MJ. Scaling-Basis Chirplet Transform. IEEE Trans Ind Electron. 2021;68: 8777–8788.

24. Li Z, Sun F, Gao J, Liu N, Wang Z. Multi-Synchrosqueezing Wavelet Transform for Time-Frequency Localization of Reservoir Characterization in Seismic Data. IEEE Geosci Remote Sens Lett. 2022;19: 1–5.

25. Daubechies I, Lu J, Wu H-T. Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool. Appl Comput Harmon Anal. 2011;30: 243–261.

26. Yu G, Wang Z, Zhao P, Li Z. Local maximum synchrosqueezing transform: An energy-concentrated time-frequency analysis tool. Mech Syst Signal Process. 2019;117: 537–552.

27. Yu G, Yu M, Xu C. Synchroextracting Transform. IEEE Trans Ind Electron. 2017;64: 8042–8054.

28. Tu X, He Z, Hu Y, Abbas S, Li F. Horizontal Synchrosqueezing Transform: Algorithm and Applications. IEEE Sens J. 2020;20: 4353–4360.

29. Zhu X, Zhang Z, Gao J, Li B, Li Z, Huang X, et al. Synchroextracting chirplet transform for accurate IF estimate and perfect signal reconstruction. Digit Signal Process. 2019;93: 172–186.

30. Hua Z, Shi J, Zhu Z. Matching Linear Chirplet Strategy-Based Synchroextracting Transform and Its Application to Rotating Machinery Fault Diagnosis. IEEE Access. 2020;8: 185725–185737.

31. He Y, Hu M, Jiang Z, Feng K, Ming X. Local maximum synchrosqueezes from entropy matching chirplet transform. Mech Syst Signal Process. 2022;181: 109476.

32. Yu G, Wang Z, Zhao P. Multisynchrosqueezing Transform. IEEE Trans Ind Electron. 2019;66: 5441–5455.