Research Article  Open Access
Haibin Wang, Junbo Long, "Applications of Fractional Lower Order Synchrosqueezing Transform Time Frequency Technology to Machine Fault Diagnosis", Mathematical Problems in Engineering, vol. 2020, Article ID 3983242, 19 pages, 2020. https://doi.org/10.1155/2020/3983242
Applications of Fractional Lower Order Synchrosqueezing Transform Time Frequency Technology to Machine Fault Diagnosis
Abstract
Synchrosqueezing transform (SST) is a high resolution time frequency representation technology for nonstationary signal analysis. The short time Fourier transformbased synchrosqueezing transform (FSST) and the S transformbased synchrosqueezing transform (SSST) time frequency methods are effective tools for bearing fault signal analysis. The fault signals belong to a nonGaussian and nonstationary alpha () stable distribution with and even the noises being also stable distribution. The conventional FSST and SSST methods degenerate and even fail under stable distribution noisy environment. Motivated by the fact that fractional low order STFT and fractional low order Stransform work better than the traditional STFT and Stransform methods under α stable distribution noise environment, we propose in this paper the fractional lower order FSST (FLOFSST) and the fractional lower order SSST (FLOSSST). In addition, we derive the corresponding inverse FLOSST and inverse FLOSSST. The simulation results show that both FLOFSST and FLOSSST perform better than the conventional FSSST and SSST under stable distribution noise in instantaneous frequency estimation and signal reconstruction. Finally, FLOFSST and FLOSSST are applied to analyze the time frequency distribution of the outer race fault signal. Our results show that FLOFSST and FLOSSST extract the fault features well under symmetric stable (SS) distribution noise.
1. Introduction
Synchrosqueezing transform is a new time frequency analysis technology for the nonstationary signals. Its principle is to calculate time frequency distribution of the signal, then squeeze the frequency of the signal in time frequency domain, and rearrange its time frequency energy, so as to improve time frequency resolution greatly. Synchrosqueezing transform mainly includes continuous wavelet transformbased synchrosqueezing transform [1], short time Fourier transformbased synchrosqueezing transform [2], and S transformbased synchrosqueezing transform [3]. Synchrosqueezing transform methods have been widely applied in seismic signal analysis [4], biomedical signal processing [5, 6], radar imaging [7], mechanical fault diagnosis, and other fields [8–11].
Daubechies et al. firstly gave synchrosqueezing transform concept based on the continuous wavelet transform and proposed a continuous wavelet transformbased synchrosqueezing transform (WSST) time frequency representation method and its inverse method. The method squeezes the time frequency energy of continuous wavelet transform in a certain frequency range to nearby instantaneous frequency of the signal, and the time frequency resolution was improved effectively [12], whereafter an adaptive wavelet transformbased synchrosqueezing transform based on WSST was brought up by Li et al. who applied a timevarying parameter to control the widths of the time frequency localization window according to the characteristics of signals [13]. The demodulated WSST and FSST methods have been proposed for instantaneous frequency estimation in [14, 15]. To improve the ability of processing the nonstationary signals with fast varying instantaneous frequency, a new demodulated high order synchrosqueezing transform method was presented in [11, 16], which can effectively show the time frequency distribution of the fault signal. Fourer et al. proposed a FSST method employing the synchrosqueezing transform and the Levenberg Marquardt reassignment in [17]; the idea of the method is very similar to the WSST method, which is reversible and adjustable. Yu et al. subsequently presented a synchroextracting short time Fourier transform, which is a postprocessing procedure of STFT [18]. Recently, they proposed an improved local maximum synchrosqueezing transform, which can discover more detailed features of the fault signals [19]. To compress and rearrange the S transform time frequency distribution of the signal, Huang et al. proposed a new S transformbased synchrosqueezing transform time frequency method employing synchrosqueezing transform and S transform, which can greatly improve time frequency resolution of S transform [20, 21]. Subsequently, they modified the calculation formula of the instantaneous frequency of the SSST time frequency method by using the second derivative of time frequency spectrum to time and frequency and proposed a new secondorder S transformbased synchrosqueezing transform method, which can obtain high time frequency resolution for the nonstationary signals whose instantaneous frequency varies nonlinearly with the time [22]. Recently, an adaptive short time Fourier transformbased synchrosqueezing transform method has been proposed with a timevarying parameter, and the corresponding 2ndorder adaptive FSST was also present in [23, 24].
Recently, it is verified that probability density function (PDF) of the mechanical bearing fault signals has an obvious trail, which is a nonstationary and nonGaussian distribution and belongs to stable distribution (); even the noises are also stable distribution [25–28]. The performance of the abovementioned methods degenerates under stable distribution environment, which even fail. Some of the ways they apply are the fractional low order time frequency representation methods to analyze the signals, such as fractional low order short time Fourier transform (FLOSTFT) [29, 30], fractional low order S transform (FLOST) [31, 32], and fractional low order WignerVille distribution [30]. However, the time frequency resolution of the methods is not very ideal and depends jointly on the geometry of the signal and the window function; even false spectral energies would be observed on the time frequency distribution at the locations where no spectral energies should exist. Hence, we propose the improved fractional low order short time Fourier transformbased synchrosqueezing transform (FLOFSST) and fractional low order S transformbased synchrosqueezing transform (FLOSSST) methods inspired by the FSST and SSST methods in this paper, and the derivation procedures of the inverse FLOFSST and inverse FLOSSST are introduced. The simulation results show that the performances of the FLOFSST and FLOSSST time frequency representation methods are superior to the existing ones under stable distribution noise environment; they have higher time frequency resolution than the existing FLOSTFT and FLOST methods and can better be suitable for the impulse noise environment than the FSST and SST methods. The IFLOFSST and IFLOSSST methods have smaller reconstruction MSEs than the IFSST and ISST methods under different and . Finally, we apply the FLOFSST and FLOSSST time frequency representation methods to analyze the bearing out race fault signal. The simulation results show that the FLOFSST and FLOSSST methods can work in Gaussian noise and stable distribution noise environment and extract the features of the outer race fault signal, which have some robustness; their performances are better than the existing FSST and SSST methods.
In this paper, the improved FLOFSST and FLOSSST time frequency representation technologies based on fractional lower order statistics and synchrosqueezing transform are presented for the bearing fault diagnosis under Gaussian and stable distribution environment. The paper is structured in the following manner. stable distribution and the bearing fault signals are introduced in Section 2. The improved fractional lower order synchrosqueezing transform methods and their inverse transforms are demonstrated in Section 3, and simulation comparisons with the existing time frequency representation methods based on secondorder statistics are performed to demonstrate superiority of the improved methods. Applications of the improved methods for the outer race fault signals diagnosis are demonstrated in Section 4. Finally, conclusions and future research are given in Section 5.
2. Stable Distribution and Bearing Fault Signals
2.1. Stable Distribution
Probability density function (PDF) of stable distribution is expressed aswhere , . stable distribution is a generalized Gaussian process, is its characteristic index, and its variance is infinite. When , which is Gaussian distribution, and when , it is low order stable distribution. is the location parameter and is the dispersion coefficient, respectively. is the symmetry parameter, when , which is called the symmetric stable distribution (). The waveforms of stable distribution are shown in Figure 1 under , 1.0, 1.5, and 2.0, and their PDFs are demonstrated in Figure 2.
(a)
(b)
(c)
(d)
(a)
(b)
2.2. Bearing Fault Signals
The real bearing fault signals data are obtained from the Case Western Reserve University (CWRU) bearing data center [33]. The experimental equipment adopts 62052RS JEM SKF type bearing, the outer race diameter is 20.472 inches, and the inner race and the ball diameter are 0.9843 inches and 0.3126 inches, respectively. The bearing outer race thickness is 0.5906 inches, motor load is 0 HP, and motor speed is 1797 rpm. The bearing faults of inner race, ball, and outer race are set, and the fault diameters are all 0.021 inches. The fault data are collected at 12,000 samples per second, and the outer race position relative to load zone centered at 6:00. The experiments are conducted with a 2 HP reliance electric motor, and the acceleration data are measured at the proximal and distal of the motor bearings; the points include the drive end accelerometer (DE), fan end accelerometer (FE), and base accelerometer (BA). The normal signal is given in Figure 3(a), and the fault signals of inner race, ball, and outer race are shown in Figures 3(b)–3(d), respectively. We can know that the waveforms of the fault signals have a certain impulse.
(a)
(b)
(c)
(d)
In order to further verify the pulse characteristics of bearing failure signals, we use stable distribution statistical model to estimate the parameters of inner race fault, ball fault, and outer race fault signals; the results are given in Table 1. As it can be seen, the characteristic index of the normal bearing signal is equal to 2, which is Gaussian distribution. However, the characteristic index of the bearing fault signals is greater than 1 but smaller than 2, and it belongs to nonGaussian stable distribution ().

PDFs of the inner race fault, ball fault, and outer race fault signals are shown in Figures 4(a)–4(f), respectively. Figures 4(b), 4(d), and 4(f) show that the PDFs of the normal have no tail, but the PDFs of the fault signals have heavy tail, and the PDFs of the fault signals in DE are especially serious. The parameters of the fault signals are approximately equal to zero in Table 1, and Figure 4 shows that PDFs of the fault signals are near symmetric. Hence, distribution is a more concise and accurate statistical model for the bearing fault signals.
(a)
(b)
(c)
(d)
(e)
(f)
3. Fractional Lower Order Synchrosqueezing Transform Methods
3.1. Fractional Lower Order STFT TransformBased Synchrosqueezing Transform Method
3.1.1. Principle
Short time Fourier transform (STFT) of the fault machinery vibration signal contaminated by distribution noise or Gaussian noise can be written asand its fractional low order short time Fourier transform (FLOSTFT) is given by [30]where is the frequency parameter; and are the time parameter. denotes the displacement parameter on the time axis. is the Gaussian window function related to the time. denotes order moment of the signal (). When is a real signal, , , and when is a complex signal, . is the characteristic exponent of distribution, and is the complex conjugate operation.
Letting , its complex conjugate function is ; according to Plancherel’s theorem and Fourier transform, we havewhere is fractional lower order Fourier transform (FLOFT) of and denotes frequency variant. is Fourier transform of , and is complex conjugate of . The right side of (4) is converted from the time domain to the frequency domain.
Letting , can be written as
Substituting (5) into (4), we have
Letting and , then , and its FLOFT can be expressed as
Substituting (7) into (6), then
Fourier transform of and makes cluster around and will be concentrated around . Substituting into (8), then
The instantaneous frequency (IF) formula of can be written as
After synchrosqueezing the frequency in (10), the discrete values can be obtained. Letting the frequency points in FLOSTFT time frequency spectrum, and . By centering and letting , the synchrosqueezing calculation is extended to the successive bins ; then, fractional lower order STFTbased synchrosqueezing transform can be defined as
The FLOFSST can “squeeze” a frequency interval to a frequency point in the time frequency domain; therefore, the process can greatly improve the time frequency resolution.
A multicomponent signal can be expressed aswhere . Then, the of th component can be expressed as
FLOSTFT is just as linear as STFT; then
The instantaneous frequency (IF) calculation method of th component may be written as
The corresponding instantaneous frequency calculation method of may be obtained by
By substituting (15) and (17) into (12), we can obtain the fractional low order STFTbased synchrosqueezing transform of the multicomponent signal.
According to the definition of inverse STFTbased synchrosqueezing transform in [18, 19], inverse fractional lower order STFT transformbased synchrosqueezing transform (IFLOFSST) of nd signal can be written byand the signal can be gotten employing .
3.1.2. The Steps of the FLOFSST Method
(i)Step 1: compute of each component for the signal employing (14).(ii)Step 2: compute instantaneous frequency of each component by substituting to(16).(iii)Step 3: solve by (17).(iv)Step 4: solve the discrete values employing .(v)Step 5: compute by substituting to (12).
3.1.3. Application Review
In this section, we design the following experiments to compare the proposed FLOFSST method with the existing STFT, FLOSTFT, and FSST methods. The simulation signal contaminated by the noise is defined aswhere is distribution noise or Gaussian noise. When the noise is distribution noise, generalized signal to noise ratio () can be used instead of , which is expressed aswhere is the dispersion coefficient of distribution noise. According to the given , the amplitude of the signal can be written as
Let , , and . The waveforms of and in time domain are shown in Figure 5. We apply the fractional lower order STFT transformbased synchrosqueezing transform method, the existing the STFT transformbased synchrosqueezing transform method, the fractional lower order STFT method, and traditional STFT method to estimate time frequency distribution of the signal under Gaussian distribution noise and stable distribution noise; the simulation results are shown in Figures 6 and 7.
(a)
(b)
(c)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
In order to compare the effectiveness of the IFSST and IFLOFSST methods, letting , where is the number of MonteCarlo experiment, is the reconstructed signal employing the IFSST method or the IFLOFSST method. Letting , the signal is reconstructed employing the IFSST and IFLOFSST methods under different ; their MSEs are shown in Figure 8(a). We apply the IFSST and IFLOFSST methods to reconstruct the signal when ; changes from to ; the simulations are shown in Figure 8(b).
(a)
(b)
3.1.4. Remarks
The STFT, FLOSTFT, FSST, and FLOFSST time frequency methods in Figure 6 all can estimate out the time frequency representation of the signal under Gaussian noise environment (), but the synchrosqueezing methods have better performance. The STFT method in Figure 7(a) and FSST time frequency method in Figure 7(c) fail under noise environment ( and ); the FLOSTFT method in Figure 7(b) can estimate out the time frequency representation of the signal , but its effect is not very ideal. The improved FLOFSST method in Figure 7(d) can better get the time frequency representation of the signal under noise environment, which has good toughness.
The STFT and FSST are unsuitable for noise environment, and the FLOSTFT method can work under noise environment, but has poor time frequency resolution and is controlled by the window function. The FSST method has better time frequency resolution, but cannot work in noise environment. The improved FLOFSST method can work under noise environment and has high time frequency resolution. As a result, the FSST time frequency method is only suitable to analyze the signals under Gaussian noise environment, but the improved FLOFSST can work under Gaussian and noise environment, which is robust.
Figure 8(a) is comparison under and different ; the experimental results show that of the IFSST method change from 2 dB to 230 dB when changes from 0.2 to 2, but of the IFLOFSST method are 2 dB. Hence, the IFLOFSST method has obvious advantage in reconstructing the signal under different ; particularly, the advantage of the IFLOFSST method is more obvious when .
Figure 8(b) is comparison under and different ; we can know that reconstruction of the IFSST method have a large variation, which changes from 14 dB to 78 dB; however, of the IFLOFSST method are stable in 2 dB. Hence, the improved IFLOFSST method has better stability in reconstructing the original signal.
3.2. Fractional Lower Order S TransformBased Synchrosqueezing Transform
3.2.1. Principle
The fault machinery vibration signal contaminated by the noise may be given bywhere is fault vibration signal and is distribution noise or Gaussian noise. Its S transform can be written asand its fractional lower order S transform was defined as [28]where is the frequency parameter and is the time parameter. denotes the displacement parameter on the time axis. is the Gaussian window function related to the frequency.
Equation (24) can be written as
Let , and its complex conjugate function is . Then, (26) changes as
The right side of (27) is converted from the time domain to the frequency domain based on Plancherel’s theorem and Fourier transform; then we obtainwhere is fractional lower order Fourier transform (FLOFT) of and denotes frequency variant. is Fourier transform of , and is complex conjugate of .
Letting and , then , and its FLOFT can be expressed as
By substituting (29) to (28), then
Fourier transform of and can assemble at around , and will be concentrated around . By substituting to (30), then
The instantaneous frequency of can be written as
By substituting (31) and (33) to (32), we can obtain the squeezed instantaneous frequency. Through synchrosqueezing the frequency with (32), the discrete values can be gotten. Letting the frequency points in FLOST time frequency spectrum, and . By centering and letting , extend the synchrosqueezing process to the successive bins ; then, fractional lower order S transformbased synchrosqueezing transform can be written as
For the calculation of IF and FLOSSST of a multicomponent signal, , we havewhere . Then, the FLOST of th component can be expressed as
FLOST is just as linear as ST; then
The IF calculation method of th component may be written as
The corresponding IF calculation method of may be obtained by
By substituting (37) and (39) into (34), fractional low order STFTbased synchrosqueezing transform time frequency representation of the multicomponent signal can be obtained.
3.2.2. Inverse Fractional Lower Order S TransformBased Synchrosqueezing Transform
Multiplying on both sides of (28) and taking the integral to , then
Let ; then
For the real signal , letting , we have
In the piecewise constant approximation corresponding to the binning in , we have
FLOST of the signal in (26) can be written aswhere is modulo operation of and is its phase position. Substituting (44) to (34), we obtain
Multiplying on both sides of (45) and letting , then
Substituting (46) into (43), we can deduce the following expression:where is real signal. According to , inverse fractional lower order S transformbased synchrosqueezing transform (IFLOSSST) of can be written as
We can reconstruct the signal in FLOSSST time frequency domain employing (48).
3.2.3. The Steps of the FLOSST Time Frequency Method
(i)Step 1: compute of each component for the signal employing (31).(ii)Step 2: solve by substituting of each component to (37).(iii)Step 3: compute instantaneous frequency of each component by substituting to (38).(iv)Step 4: solve of the signal by substituting to (39).(v)Step 5: compute the discrete values employing .(vi)Step 6: solve of the signal by substituting to (34).
3.2.4. Application Review
In this section, in (19) is used as the experiment signal. The proposed fractional lower order S transformbased synchrosqueezing transform method, the existing the S transformbased synchrosqueezing transform method, the fractional lower order S transform method, and traditional S transform method are used to display time frequency distribution of the signal under Gaussian distribution noise () and stable distribution noise (and ); the simulation results are shown in Figures 9 and 10.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Letting and , the ISSST and IFLOSSST methods are used to reconstruct the original signal; the results are shown in Figure 11. In order to further compare the effectiveness of the ISSST and IFLOSSST methods, letting , the signal is reconstructed employing the ISSST and IFLOSSST methods under different ; their are shown in Figure 12(a). We apply the ISSST and IFLOSSST methods to reconstruct the signal when ; changes from to ; the simulations are shown in Figure 12(b).
(a)
(b)
(c)
(a)
(b)
3.2.5. Remarks
Figure 9 is the time frequency representations of the signal under Gaussian noise environment () employing the ST, FLOST, SSST, and FLOSSST methods, respectively. We can see that all the methods can estimate out the time frequency distribution of the signal , but the synchrosqueezing transform methods have obvious advantages in time frequency resolution. The time frequency representations of the signal employing the ST, FLOST, SSST, and FLOSSST methods under noise environment ( and ) are shown in Figure 10. The results show that the ST method in Figure 10(a) and SSST method in Figure 10(c) fail; the FLOST method in Figure 10(b) can estimate out the time frequency distribution of the signal , but its effect is not very ideal. The improved FLOSSST method in Figure 10(d) can better get the time frequency representation of the signal , which has higher time frequency resolution.
The reconstructed signal employing the ISSST method is shown in Figure 11(b) under noise environment (and ); it can be seen that the signal is covered by noise; the ISSST method fails. Figure 11(b) is the reconstructed signal based on the IFLOSSST method under the same conditions; it shows that the reconstructed signal is very similar to the original signal . Figure 12(a) is reconstruction comparison under when changes from 0.2 to 2; the results show that the reconstruction of the IFSST method change from 1 dB to 290 dB, but the reconstruction of the IFLOFSST method have an obvious low level, which are stable at about 2 dB. Hence, the IFLOFSST method has obvious advantage in reconstructing the signal under different ; particularly, the advantage of the IFLOFSST method is more obvious when . Figure 12(b) is the reconstruction comparison under when changes from 14 dB to 78 dB; it shows that the reconstruction of the IFSST method have a large variation, but the reconstruction MSEs of the IFLOFSST method change from −2 dB to 8 dB. Hence, the improved IFLOFSST method has better stability in reconstructing the signal.
As a result, the SSST time frequency method and the ISSST signal reconstruction method are only suitable to analyze and reconstruct the signals under Gaussian noise environment, but the improved FLOSSST and IFLOSSST methods can work in Gaussian and stable distribution noise environment, which are robust.
4. Application Simulations
In this simulation, the experiment signal adopts the bearing outer race fault signal (DE) in Section 2. 0.2 seconds’ data is selected as the test signal, which is collected at 12,000 samples per second, and . The improved FLOFSST and FLOSSST methods are applied to analyze time frequency distribution of the outer race fault signal; the simulation results are shown in Figure 13.
(a)
(b)
Figures 13(a) and 13(b) are the time frequency representations of the outer race fault signal employing the FLOFSST and FLOSSST methods, respectively. It can be seen that two methods have a good lateral resolution, the lowfrequency shock pulse mainly includes 0 Hz to 4000 Hz, and the dominant frequency of the vibration components is approximately 600 Hz, 2800 Hz, and 3500 Hz. All the methods have a good vertical resolution; the gap between the impacts can be clearly seen, which regularly change. The time interval in A, B, C, D, E, and F is about 30 ms; the corresponding characteristic frequency of the outer race fault signal is about 33.333 Hz.
In order to further prove the advantages of the improved FLOFSST and FLOSSST methods, distribution noise ( and ) is added in the stable distribution outer race fault signal as the background noise of actual working environment. The improved methods and existing methods are applied to demonstrate time frequency representation of the outer race fault signal; the simulations are shown in Figure 14. The results show that the FSST method in Figure 14(a) and SST method in Figure 14(b) fail. However, the FLOFSST method in Figure 14(c) and FLOSSST method in Figure 14(d) can give out time frequency distribution of the fault signal under substandard conditions, which have certain ability in the horizontal and vertical time frequency representation, and we can know the dominant frequency and the time interval in A, B, C, D, E, and F, but the overall resolution is not high and needs to improve. Hence, the improved fractional lower order synchrosqueezing methods have better performance superiority than the existing synchrosqueezing methods, which are more suitable for fault analysis in complex environment and are robust.
(a)
(b)
(c)
(d)
5. Conclusions
stable distribution is a more appropriate statistical model for the bearing fault signals. STFT transformbased synchrosqueezing transform and S transformbased synchrosqueezing transform are two new time frequency representation methods for nonstationary signal processing; their time frequency resolution can be greatly improved by rearranging the time frequency energy of the signals. In order to make the FSST and SST methods applicable to Gaussian and stable distribution noise environment, the improved FLOFSST and FLOSSST time frequency representation methods are proposed by employing fractional low order statistics. The performances of the FLOFSST and FLOSSST methods are superior to the existing time frequency analysis methods; they have higher time frequency resolution than the existing FLOSTFT and FLOST methods because of the synchrosqueezing processing and can better suppress the impulse noise than the FSST and SST methods. The IFLOFSST and IFLOSSST methods have smaller reconstruction than the IFSST and ISSST methods under different and . We can apply the improved methods to analyze the stable distribution bearing fault signal; even stable distribution noise environment, the fault characteristic frequency, the dominant frequency, and the other fault frequency features of the fault signals can be clearly obtained. In the future, we can also further study time frequency filtering technology based on the proposed IFLOFSST and IFLOSSST methods, and the methods have a good application prospect in the field of the bearing fault analysis and detection.
Data Availability
The data used to support the findings of this study are provided in the Supplementary Materials.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was financially supported by the Natural Science Foundation of China (61962029), the Natural Science Foundation of Jiangxi Province, China (20192BAB207002), the Science and Technology Project of Provincial Education Department of Jiangxi (GJJ170954), and the Science and Technology Project of Jiujiang University, China (2014SKYB009).
Supplementary Materials
This section contains the original experimental data of this paper. (Supplementary Materials)
References
 X.L. Wang, C.L. Li, and X. Yan, “Nonstationary harmonic signal extraction from strong chaotic interference based on synchrosqueezed wavelet transform,” Signal, Image and Video Processing, vol. 13, no. 2, pp. 397–403, 2019. View at: Publisher Site  Google Scholar
 A. Rueda and S. Krishnan, “Augmenting dysphonia voice using Fourierbased synchrosqueezing transform for a CNN classifier,” in Proceedings of the IEEE International Conference on Acoustics, Speech and Signal ProcessingProceedings, pp. 6415–6419, Brighton, UK, May 2019. View at: Google Scholar
 G. Pan, S. Li, and Y. Zhu, “A timefrequency correlation analysis method of time series decomposition derived from synchrosqueezed S transform,” Applied Sciences, vol. 9, no. 4, p. 777, 2019. View at: Publisher Site  Google Scholar
 Q. Wang, J. Gao, N. Liu, and X. Jiang, “Highresolution seismic timefrequency analysis using the synchrosqueezing generalized Stransform,” IEEE Geoscience and Remote Sensing Letters, vol. 15, no. 3, pp. 374–378, 2018. View at: Publisher Site  Google Scholar
 H.T. Wu, Y.H. Chan, Y.T. Lin, and Y.H. Yeh, “Using synchrosqueezing transform to discover breathing dynamics from ECG signals,” Applied and Computational Harmonic Analysis, vol. 36, no. 2, pp. 354–359, 2014. View at: Publisher Site  Google Scholar
 P. Ozel, A. Akan, and B. Yilmaz, “Synchrosqueezing transform based feature extraction from EEG signals for emotional state prediction,” Biomedical Signal Processing and Control, vol. 52, pp. 152–161, 2019. View at: Publisher Site  Google Scholar
 B. Zang, M. Zhu, X. Zhou, and L. Zhong, “Application of Stransform in ISAR imaging,” Electronics, vol. 8, no. 6, p. 676, 2019. View at: Publisher Site  Google Scholar
 J. Li, M. Huang, P. Xie et al., “Synchrosqueezingcross wavelet transform and enhanced fault diagnosis of rolling bearing,” Acta Metrologica Sinica, vol. 39, no. 2, pp. 237–241, 2018. View at: Google Scholar
 S. Meignen, T. Oberlin, and D.H. Pham, “Synchrosqueezing transforms: from low to highfrequency modulations and perspectives,” Comptes Rendus Physique, vol. 20, no. 5, pp. 449–460, 2019. View at: Publisher Site  Google Scholar
 J. Hu, L. Chai, D. Xiong, and W. Wang, “A novel method of realizing stochastic chaotic secure communication by synchrosqueezed wavelet transform,” Digital Signal Processing, vol. 82, pp. 194–202, 2018. View at: Publisher Site  Google Scholar
 Y. Hu, X. Tu, and F. Li, “Highorder synchrosqueezing wavelet transform and application to planetary gearbox fault diagnosis,” Mechanical Systems and Signal Processing, vol. 131, pp. 126–151, 2019. View at: Publisher Site  Google Scholar
 I. Daubechies, J. Lu, and H.T. Wu, “Synchrosqueezed wavelet transforms: an empirical mode decompositionlike tool,” Applied and Computational Harmonic Analysis, vol. 30, no. 2, pp. 243–261, 2011. View at: Publisher Site  Google Scholar
 L. Li, H. Cai, and Q. Jiang, “Adaptive synchrosqueezing transform with a timevarying parameter for nonstationary signal separation,” Applied and Computational Harmonic Analysis, 2019, In press. View at: Publisher Site  Google Scholar
 S. Wang, X. Chen, G. Cai, B. Chen, X. Li, and Z. He, “Matching demodulation transform and synchrosqueezing in timefrequency analysis,” IEEE Transactions on Signal Processing, vol. 62, no. 1, pp. 69–84, 2014. View at: Publisher Site  Google Scholar
 Q. Jiang and B. W. Suter, “Instantaneous frequency estimation based on synchrosqueezing wavelet transform,” Signal Processing, vol. 138, pp. 167–181, 2017. View at: Publisher Site  Google Scholar
 X. Tu, Y. Hu, F. Li, S. Abbas, Z. Liu, and W. Bao, “Demodulated highorder synchrosqueezing transform with application to machine fault diagnosis,” IEEE Transactions on Industrial Electronics, vol. 66, no. 4, pp. 3071–3081, 2019. View at: Publisher Site  Google Scholar
 D. Fourer, F. Auger, and P. Flandrin, “Recursive versions of the LevenbergMarquardt reassigned spectrogram and of the synchrosqueezing STFT,” in Proceedings of the 41st IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, March 2016. View at: Publisher Site  Google Scholar
 G. Yu, M. Yu, and C. Xu, “Synchroextracting transform,” IEEE Transactions on Industrial Electronics, vol. 64, no. 10, pp. 8042–8054, 2017. View at: Publisher Site  Google Scholar
 G. Yu, Z. Wang, P. Zhao, and Z. Li, “Local maximum synchrosqueezing transform: an energyconcentrated timefrequency analysis tool,” Mechanical Systems and Signal Processing, vol. 117, pp. 537–552, 2019. View at: Publisher Site  Google Scholar
 Z.l. Huang, J. Zhang, T.h. Zhao, and Y. Sun, “Synchrosqueezing Stransform and its application in seismic spectral decomposition,” IEEE Transactions on Geoscience and Remote Sensing, vol. 54, no. 2, pp. 817–825, 2016. View at: Publisher Site  Google Scholar
 Z. Huang and J. Zhang, “Synchrosqueezing Stransform,” SCIENTIA SINICA Informationis, vol. 46, no. 5, pp. 643–650, 2016. View at: Publisher Site  Google Scholar
 Z. Huang, J. Zhang, and Z. Zhou, “A second order synchrosqueezing Stransform and its application in seismic spectral decomposition,” Chinese Journal of Geophysics, vol. 60, no. 7, pp. 2833–2844, 2017. View at: Publisher Site  Google Scholar
 L. Li, H. Y. Cai, H. X. Han, Q. T. Jiang, and H. B. Ji, “Adaptive shorttime Fourier transform and synchrosqueezing transform,” Signal Processing, vol. 166, Article ID 107231, 2020. View at: Publisher Site  Google Scholar
 H. Y. Cai, Q. T. Jiang, L. Li, and B. W. Suter, “Analysis of adaptive short time Fourier transformbased synchrosqueezing transform,” Analysis and Applications, 2020, In press. View at: Publisher Site  Google Scholar
 M. Shao and C. L. Nikias, “Signal processing with fractional lower order moments: stable processes and their applications,” Proceedings of the IEEE, vol. 81, no. 7, pp. 986–1010, 1993. View at: Publisher Site  Google Scholar
 G. Yu, C. Li, and J. Zhang, “A new statistical modeling and detection method for rolling element bearing faults based on alphastable distribution,” Mechanical Systems and Signal Processing, vol. 41, no. 12, pp. 155–175, 2013. View at: Publisher Site  Google Scholar
 J. Long, H. Wang, and Li Peng, “Applications of fractional lower order frequency spectrum technologies to bearing fault analysis,” Mathematical Problems in Engineering, vol. 2019, Article ID 7641383, 24 pages, 2019. View at: Publisher Site  Google Scholar
 Q. Xiong, W. Zhang, Y. Xu, Y. Peng, and P. Deng, “Alphastable distribution and multifractal detrended fluctuation analysisbased fault diagnosis method application for axle box bearings,” Shock and Vibration, vol. 2018, Article ID 1737219, 12 pages, 2018. View at: Publisher Site  Google Scholar
 Y. Yang, X. Sun, and Z. Zhong, “A parameter estimation algorithm for frequencyhopping signals with alpha stable noise,” in Proceedings of the 3rd IEEE Advanced Information Technology, Electronic and Automation Control Conference (IAEAC), Chongqing, China, October 2018. View at: Google Scholar
 J. Long, H. Wang, P. Li, and H. Fan, “Applications of fractional lower order timefrequency representation to machine bearing fault diagnosis,” IEEE/CAA Journal of Automatica Sinica, vol. 4, no. 4, pp. 734–750, 2017. View at: Publisher Site  Google Scholar
 T. Shu, X. Yu, D. Zha, and Q. TianShuang, “Timefrequency distribution of generalized Stransform under alpha stable distributed noise conditions,” Journal of Signal Processing, vol. 30, no. 6, pp. 634–641, 2014. View at: Google Scholar
 J. Long, H. Wang, D. Zha, P. Li, H. Xie, and L. Mao, “Applications of fractional lower order S transform time frequency filtering algorithm to machine fault diagnosis,” PLoS One, vol. 12, no. 4, Article ID e0175202, 2017. View at: Publisher Site  Google Scholar
 CWRU bearing data center, http://csegroups.case.edu/bearingdatacenter/pages/downloaddatafile.
Copyright
Copyright © 2020 Haibin Wang and Junbo Long. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.