Novel complete ensemble EMD with adaptive noisebased hybrid filtering for rolling bearing fault diagnosis
Xiaojun Song^{1} , Hongwei Sun^{2} , Liwei Zhan^{3}
^{1, 2}Power Engineering, Harbin Power Vocational Technology College, Harbin, 150030, People’s Republic of China
^{3}Aero Engine Corporation of China Harbin bearing Co., Ltd, Harbin, 150500, People’s Republic of China
^{3}Corresponding author
Journal of Vibroengineering, Vol. 21, Issue 7, 2019, p. 18451858.
https://doi.org/10.21595/jve.2019.20100
Received 24 July 2018; received in revised form 20 February 2019; accepted 7 April 2019; published 15 November 2019
JVE Conferences
A feature extraction of fault bearing has attracted considerable attention in recent years. However, weak fault feature is difficult to extract under heavy background noise. To solve this problem, a novel multilayer filtering method is proposed to filter out noise and extract weak fault feature. The first layer introduces a metric based on detrended fluctuation analysis (DFA) to identify intrinsic mode function (IMF) that reflect period impulsive information for vibration signal adaptively. The second layer uses nonlocal mean (NLM) method as a prefilter of the third layer to realize extraction of singular value decomposition (SVD) which reflect the most information of IMFs. The last layer introduces a relative energy difference criterion of a singular value to extract important feature of Hankel matrix of IMFs. The filtered signal is obtained by reconstructed signal from identified singular value of SVD. Experiment results on simulation and real vibration signals indicate that the hybrid filtering method removes heavy noise successfully and extract weak fault feature of rolling bearing effectively.
 It proposes criterion of the smallest Hurst exponent to identify the periodic impulse adaptively.
 It introduces the nonlocal meaning filter to filter out noise in the intrinsic mode functions.
 The relative energy difference spectrum is proposed to remove the weak noise component.
Keywords: hybrid filtering, detrended fluctuation analysis, relative energy difference, fault feature.
1. Introduction
Bearing is a critical component in rotary machine. The reliability directly affects machine health status. Thus, an effective condition monitoring (CM) and fault diagnosis (FD) method is required to detect and recognize potential abnormalities which reduce risk of unexpected severity and other disastrous events [14].
Vibration analysis technology is a key topic in CM and FD research domain. This technology may extract fault feature of rolling bearing effectively. Conventional methods, such as fast Fourier transform (FFT), are used to address stationary and linear signals [5] widely. However, FFT is unsuitable to extract fault feature on account of characteristic of nonlinear and nonstationary for vibration signal. Wavelet transform is a typical application in bearing fault diagnosis [68]. This technique has attracted considerable attention of researchers in the past two decades. However, wavelet transform is difficult to select wavelet base function. Recently, empirical mode decomposition (EMD) has been proposed to process complex nonlinear and nonstationary signal [9]. EMD is an adaptive decomposition method that decomposes complicated signals into a set of multicomponent with intrinsic mode function (IMF). Every IMF indicates a monocomponent amplitude and frequencymodulated signal. EMD has been applied to bearing [10], gearbox [11], and rubbing [12] fault diagnoses. However, EMD have several limitations, such as mode mixing and end effect. To solve the problem, Huang [13] proposed an ensemble EMD (EEMD) to alleviate the limitations of EMD. This method reflected the frequency of every IMF accurately and effectively. Moreover, EEMD has attracted significant attention in fault diagnoses domain. For example, Wang [14] proposed a method based on multiscale principal component analysis to extract fault frequency of slewing bearing. Chen [15] improved EEMD by adding white noises and adaptively determining an ensemble number to extract fault features. However, EEMD has caused spurious mode and residual noise because of interaction of noise and signal. Subsequently, Torres [16] introduced a complete EEMD with adaptive noise (CEEMDAN) to overcome this drawback. CEEMDAN has been applied to fault diagnoses [17, 18] successfully.
For vibration signal being immersed by heavy noise, it is important to develop a denoised method to extract the early detection of weak fault. However, filtering method based on EMD or its improved versions (EEMD and CEEMDAN) do not filter out noise in vibration signal fully. The reason is that IMF do not effectively reflect frequency component of fault because of noise interference. Recently, the hybrid filtering method based on singularvalue decomposition (SVD) and EMD or its improved versions have been proposed to extract weak fault feature. Cheng [19] used EMD to decompose impulse signal, and used SVD and support vector machine to extract fault feature and classify fault pattern of gearbox and bearing. Yang [20] adopted EEMD with sample entropy and SVD to diagnose fault by acoustic emission signal. It is found that SVD is an effective method to filter out noise and extract fault feature. However, important fault information losses because of spurious mode and heavy noise interference in IMF.
To solve this problem, this study introduces a multilayer filtering (hybrid filtering) method to extract weak fault feature of bearing. In the first layer, vibration signal is decomposed to obtain IMFs by CEEMDAN. Then, a detrended fluctuation analysis (DFA) is utilized to identify effective IMF and spurious mode adaptively. In second layer, a nonlocal mean (NLM) denoised method is used to filter out noise in each IMF. In last layer, SVD is used to remove residual noise in vibration signal. The hybrid filtering method is applied to simulated and measured vibration signal from rolling bearing. The result indicates that this method extracts weak fault effectively and operate well for bearing fault diagnose.
The remainder of this paper is organized as follows: Section 2 introduces the CEEMDAN algorithm. Section 3 describes the hybrid filtering algorithm. Section 4 verifies the effectiveness of the proposed algorithm for fault extraction. Section 5 presents the conclusion.
2. CEEMDAN algorithm
Traditional EEMD have high computational cost and residue noise because of interaction of signal and added noise. To solve this problem, CEEMDAN is introduced in reference [16]. Principle of CEEMDAN is as follows:
(1) Generate the noisy signal ${x}^{i}\left(t\right)$, ${x}^{i}\left(t\right)=x\left(t\right)+{\omega}_{0}{\xi}^{i}\left(t\right)$ ($i=$1,…, $N$), where ${\xi}^{i}\left(t\right)$ is white noise with unit variance, and ${\omega}_{0}$ is coefficient of added white noise.
(2) Obtain the first IMF (${c}_{1}^{i}\left(t\right)$) of each noisy signal ${x}^{i}\left(t\right)$ by EMD;
(3) Obtain the first decomposed IMF (${c}_{1}\left(t\right)$) by taking average of each ${c}_{1}^{i}\left(t\right)$:
(4) Compute the first residue by Eq. (2);
(5) Apply EMD to decompose ${r}_{1}\left(t\right)+{w}_{1}{E}_{1}\left({\epsilon}^{i}\left(t\right)\right)$, and extract the first IMF to obtain the decomposed second IMF (${c}_{2}\left(t\right)$), where ${E}_{m}\left(\cdot \right)$ signifies that EMD is used to decompose ${x}^{i}\left(t\right)$ to get the $m$th IMF:
(6) Compute the $m$th residue ($m=$ 2,…, $K$), and extract the first decomposed IMF to obtain $\left(m+1\right)$th decomposed IMF by Eq. (4):
(7) Repeat above steps until the residue $R\left(t\right)$ contain less than two extrema:
The signal $x\left(t\right)$ can be written as:
3. Hybrid filtering method
3.1. Selection of spurious IMFs
When CEEMDAN is used to decompose vibration signal, the obtained IMF contains periodic impulsive component (IPIC) and nonIPIC (NIPIC). The IMF with impulsive component can reflect physical meaning of vibration signal, and the residue is a noneffective IMF, which is called “spurious mode” or “spurious IMF”. It is necessary to develop a method to distinguish IPIC and NIPIC effectively. It is found that there is obvious fluctuation difference between IPIC and NIPIC. Recently, DFA [2123] has been regarded as analysis tool of scale and is suitable to analyze fluctuation of nonlinear and nonstationary time series. Thus, DFA is proposed to evaluate the fluctuation of IPIC and NIPIC. The following is the principle of DFA:
(1) Generate the integrated time series, where $\stackrel{}{x}$ indicates the mean of time series $x\left(i\right)$:
(2) Divide $y\left(k\right)$ into $n$ length sections;
(3) Use least square method to determine the local trend ${y}_{n}\left(k\right)$;
(4) Get fluctuation function $F\left(n\right)$ by ${y}_{n}\left(k\right)$:
(5) Obtain different $F\left(n\right)$ using different length segments;
(6) Calculate slope between $\mathrm{l}\mathrm{o}\mathrm{g}\left(F\left(n\right)\right)$ and $\mathrm{l}\mathrm{o}\mathrm{g}k$, where the slope is called Hurst exponent $\alpha $:
Reference [24] indicated that the larger value of Hurst exponent ($\alpha $) is, the smoother a time series is. Numerous experiments found that NIPIC is smoother than IPIC. It means that Hurst exponent ($\alpha $) decreases with the increase of IMF index. When the Hurst exponent ($\alpha $) of the $\left(k+\text{1}\right)$th IMF become large suddenly, and the IMF that corresponds to the index is a spurious mode. Then, effective vibration signals with impulsive components are adaptively constructed by subtracting first $k$th IMFs. The identified steps are summarized as follows:
(1) Decompose signal $x\left(t\right)$ into a set of IMFs by CEEMDAN:
(2) Calculate Hurst exponent (${\alpha}_{i}$) of every IMF by DFA ($i=$ 1,…, $N$);
(3) Identify index of the relevant $k$th from ${\alpha}_{i}$:
(4) Extract all IMFs up to index $k$th.
To verify the effectiveness of proposed method, the following simulated vibration signal [25] is used to identify IPIC and NIPIC:
where, $y\left(t\right)$ is noisy signal, $x\left(t\right)$ indicates the real signal and $n\left(t\right)$ is the noise. Eq. (13) and (14) is expression of $x\left(t\right)$:
here, $a$ is the modulated factor, and $T$ is the sample period. $k$ is the sample points ($k=$ 1000), ${f}_{m}$ is fault frequency of bearing, ${f}_{r}$ is the rotation frequency, ${f}_{c1}$ and ${f}_{c2}$ are the carrier frequencies from system resonance. Table 1 are the parameters of simulated vibration signal. Where the sample rate is 1.6 kHz, and the input signaltonoise ratio $(SN{R}_{in}$) is −10 dB.
Table 1. Parameters of the simulated vibration signal
$a$

$T$ (s)

${f}_{m}$ (Hz)

${f}_{r}$ (Hz)

${f}_{c1}$ (Hz)

${f}_{c2}$ (Hz)

800

1/16000

160

30

2000

7000

Fig. 1(a) and (b) are the simulated noisy signal and the corresponding fast Fourier transform (FFT), respectively. It is found that FFT cannot identify characteristic frequency of simulated fault signal, and extract only the carrier frequencies (${f}_{c1}$ and ${f}_{c2}$).
To extract the fault frequency, CEEMDAN is first used to decompose noisy vibration signals. Fig. 2 are the decomposed IMFs. It is found that there are nine IMFs and residue. The first five IMFs contain PICs obviously, and periodic information of fifth IMF is shown red box in the fifth IMF. However, it is obvious that the sixth IMF mainly contains NPIC. To identify the PIC, DFA is used to calculate Hurst exponent ($\alpha $) of each IMF to distinguish PIC and NPIC. Fig. 3 demonstrates the DFA of IMF1, and Fig. 4 exhibits the Hurst exponent ($\alpha $) for all IMFs. It is found that Hurst exponent ($\alpha $) decreases with increase of IMF index. The Hurst exponent ($\alpha $) of IMF6 rapidly become large when Hurst exponent ($\alpha $) of IMF5 is smallest. It indicates that the spurious modes are from IMF6 to IMF10 (expected residue). The identified result is in accordance with the analysis result in Fig. 2. It means that the proposed DFA method is suitable to identify IPIC and NIPIC. Fig. 5 presents the combination of signals with first five IMFs. It is found that the signal still contains heavy noise.
Fig. 1. Noisy signal and corresponding FFT
a) Noisy signal
b) Noisy signal FFT
Fig. 2. CEEMDAN decomposed IMFs for simulation vibration signal
a)
b)
Fig. 3. DFA of IMF1
Fig. 4. Hurst exponent of each IMF
Fig. 5. Summed signal of the first five IMFs
3.2. Nonlocal mean (NLM) denoising method
In this section, NLM method is used to filter out noise in each IMF. Furthermore, the method also improves effectiveness of third layer filtering to enhance feature extraction of fault signal. NLM [2628] is first used in image processing to process 2D signal rather than 1D signal. In this paper, NLM is introduced to be considered as prefilter of third layer filtering. The following is the principle of NLM:
where $x\left(t\right)$ is the real signal, and $n\left(t\right)$ is the noisy signal. Recovered signal $\stackrel{~}{x}\left(t\right)$ is the weighted sum of value at point $t$ being located in “search neighborhood” $N\left(s\right)$ and is expressed as:
where, $s$ is sample point. The expression of $Z\left(s\right)$ is the following:
The weights $w\left(s,t\right)$ are as follows:
where, $\lambda $ is the bandwidth, $\mathrm{\Delta}$ is the local patch of surrounding $s$ and contains ${L}_{\mathrm{\Delta}}$ length samples.
The second filtering layer principle is expressed as follows:
(1) The NLM method is used to filter out noise of identified PIC;
(2) The filtered PIC are combined to obtain the second layer denoised signal.
Fig. 6 depicts filtered result. It is found that periodic impulse component is more obvious than signal in Fig. 5 (parameter in NLM is determined on the basis of SURE criterion [29]).
3.3. Singular value decomposition denoising
3.3.1. Singular value decomposition principle
The final layer filtering utilizes SVD to remove extra noise. SVD demonstrates an excellent performance in noise removal and is extensively used in fault extraction of rolling bearing [3032]. It is assumed that noisy series $x\left(k\right)$ is expressed as:
Fig. 6. Second filtering layer
Signal $x\left(k\right)$ is consisted of real signal and noise. It is expressed as follows:
where, $s\left(k\right)$ is real signal, and $n\left(k\right)$ is noise. The following is Hankel matrix $H\left(H\in {R}^{p\times q}\right)$ of signal $x\left(k\right)$:
where, $p+q1=N$, $p\ge q$, and $H\left(i,j\right)=x\left(i+j1\right)$. Eq. (21) is furtherly expressed as:
The Eq. (22) is simplified as:
here, ${H}_{s}\left(i,j\right)=s\left(i+j1\right)$ and ${H}_{n}\left(i,j\right)=n\left(i+j1\right)$ are Hankel matrices of the real signal and noisy signal, respectively. Thus, the main task is to found ${H}_{s}$ to obtain denoised signal $\widehat{s}\left(k\right)$. Here, $H$ is first decomposed into the following expression:
here, $U\left(U\in {R}^{p\times q}\right)$ and $V\left(V\in {R}^{p\times q}\right)$ are two orthogonal matrices, $D$ is a diagonal matrix and defined as $D=\left[\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({\sigma}_{1},{\sigma}_{2},\cdots ,{\sigma}_{q}\right),O\right]$, $O$ is a zero matrix, and ${\sigma}_{1}>{\sigma}_{2},\cdots ,{\sigma}_{i},\cdots ,>{\sigma}_{q}>0$. Here, ${\sigma}_{i}$ is assumed to be last singular value of Hankel matrix ${H}_{s}$. Other values are set to zero if the previous $i$ singular values are preserved. Then, reverse reconstruction can obtain the denoised signal. Thus, parameter $i$ must be determined. The next part uses a relative energy difference spectrum of a singular value to identify important parameter $i$.
3.3.2. Relative energy difference spectrum
In principle, real signal $s\left(k\right)$ and noise $n\left(k\right)$ are noncorrelated. The energy of the real signal is concentrated. However, the noise is relatively decentralized. It indicates that the singular values of real signal components are rounded to certain numerical value, and singular value of noise components are around other numerical value. It means that the energy of first $i$ singular value (singular value of real signal component) is different from energy of the others ($ri$) singular value (noise component). If the energy difference of adjacent singular value is done, it is found that the energy difference of singular value are almost same as energy for real signal component or noise component, respectively. However, the energy difference of singular value is biggest in demarcation point of real signal component and noise component. It proves that the relative energy difference spectrum of SVD changes at the demarcation point of singular value of the real and noise signal components abruptly. This paper introduces the following relative energy difference spectrum to identify singular value ${\sigma}_{i}$:
when the singular value of real signal components are extracted by identifying demarcation point by proposed method, the denoised signal can been obtained by Eq. (24). It means that the signal can been reconstructed by finding the maximum ${\sigma}_{i}$. Here, the size of Hankel matrix is determined on the basis of Reference [33].
Now, the proposed denoising method is used to filter out residual noise of vibration signal in Fig. 6. Fig. 7 displays singular value of Hankel matrix. Figs. 8(a) and (b) depict the local large of singular value (first 50 singular values are displayed) and the corresponding relative energy difference spectrum, respectively. It is found that the 10th singular value is the maximum in all different spectra. Thus, the first 10 singular values are selected to reconstruct the denoised signal. Fig. 9 are comparison result of secondlayer and thirdlayer filtered signal. It is found that periodical information of the third layer filtered signal is more obvious than that of the second layer.
Fig. 7. Singular value
Fig. 8. Singular value and corresponding difference spectrum
a) Singular value
b) Difference spectrum of singular value
Fig. 9. Comparing of filtered signal with second the third layer filtering method
Analysis technique of envelope spectrum is used to extract fault feature. Fig. 10 presents the envelope spectrum of the filtered signal. It is found that fault characteristics $\text{1}X$ and $\text{2}X$ (159.7 and 319.4 Hz, respectively) and rotational frequency $\text{2}X$ (63.87 Hz) are extracted in the envelope spectrum. It indicates that the proposed hybrid filtering method is effective in extracting fault frequency of bearing.
Fig. 10. Envelope spectrum of hybrid filtered signal
Finally, the threelayer denoising method is presented as follows:
(1) Decompose vibrated signal $x\left(t\right)$ into a set of IMFs by CEEMDAN;
(2) Select the real and spurious IMF component by DFA adaptively;
(3) Filter out noise in each real IMF by NLM;
(4) Combine the filtered IMFs to obtain effective signal ${x}^{\mathrm{\text{'}}}\left(t\right)$;
(5) Use SVD to filter out noise to obtain signal ${x}^{\mathrm{\text{'}}\mathrm{\text{'}}}\left(t\right)$;
(6) Calculate the envelope spectrum of signal ${x}^{\mathrm{\text{'}}\mathrm{\text{'}}}\left(t\right)$;
(7) Extract the weak feature of vibration signal ${x}^{\mathrm{\text{'}}\mathrm{\text{'}}}\left(t\right)$;
Fig. 11 is flow chart of the hybrid filtering method.
Fig. 11. Flow chart of hybrid filtering method
4. Results and discussions
In this section, the fault data from the website of Case Western Reserve University [34] is used to evaluate the performance of proposed method. Fig. 12 exhibits the experiment rig which consists of a 2 HP motor, a torque transducer, a dynamometer, and control electronics. The test bearing is 62052RS SKF deep groove ball bearing. Table 2 is the bearing parameters.
Table 2. Bearing parameters for 62052RS deep groove ball bearing
Bearing type

Ball number $n$

Pitch diameter $D$ (mm)

Ball diameter $d$ (mm)

Contact angle $\alpha $ (°)

62052RS

9

52

8

0

Fig. 12. Bearing experiment rig: (1) – motor, (2) – torque transducer, (3) – dynamometer
The theoretical characteristic frequency of inner race ${f}_{in}$ can be calculated using the following expression:
where, ${f}_{r}$ is the rotation frequency.
In this experiment, single point fault is introduced to the test bearing using electrodischarge machining with fault diameters of 0.018, 0.036, 0.054, and 0.072 mm. Fault data from fault diameter of 0.018 mm is used to identify fault recognition as early as possible. The sample rate is 12 kHz. The selected data length is 1000. Figs. 13(a) and (b) display the time domain signal and corresponding FFT, respectively. It is found that FFT of vibration signal do not extract frequency feature of vibration signal, and the vibration signal contain heavy noise. Here, the proposed method is first used to remove noise by flow chart in Fig. 11, then analysis technology of envelope spectrum of vibration signal is employed to extract fault feature. Moreover, typical filtering methods, namely, EMDSVD [25], SVD, and NLM, are also used to extract fault feature for the same vibration signal. The filtered result and the corresponding envelope spectrum are demonstrated in Figs. 14(b), (c), and (d). It is found that more noise is observed using EMDSVD, SVD, and NLM methods than proposed method (CEEMANDFANLMSVD) for a time domain waveform. Moreover, for envelope spectrum, it also has similar result and contains a noisy frequency component (highfrequency component or interference frequency as shown elliptical annotation in red box in Fig. 14) among EMDSVD, SVD, and NLM methods (envelope spectrum is only displayed in the first 2 kHz). However, the envelope spectrum of the proposed method does not reflect a high frequency and only reflect fault frequency of bearing (lowfrequency information in Fig. 14). It indicates that the proposed method more effectively remove noise than other denoising methods. The envelope spectrum of the CEEMANDFANLMSVD presents the maximum frequency (155.7 Hz) and multiplication frequency (311.4 Hz). The identified frequency is nearly consistent with the theoretical frequency (151.44 Hz) by Eq. (26). For the EMDSVD, SVD method, although the fault frequency and corresponding multiplication frequency can be reflected in Fig. 14, the multiplication frequency is difficult to identify because of interference frequency. For NLM method, it is not found fault frequency and multiplication frequency in the envelope spectrum, and it also produce the interference frequency in red box in Fig. 14. In addition, the amplitude of envelope spectrum of proposed method is bigger than the other methods. It makes the fault frequency be more identified. Based on above analysis, the proposed method demonstrates an outstanding performance than the other filtering methods.
Fig. 13. Noisy signal and corresponding FFT
a) Noisy signal
b) Noisy signal FFT
Fig. 14. Comparing of filtered results and feature extraction
5. Conclusions
This study proposes a novel hybrid filtering method based on the CEEMDANDFANLMSVD to extract weak fault features under heavy background noise. On the basis of the abovementioned research results, the conclusions are as follows:
1) The first layer filtering introduces criterion of the smallest Hurst exponent to identify the periodic impulse among the IMFs by CEEMDAN adaptively.
2) The second filtering layer uses NLM to filter out noise in each IMF that contains a PIC.
3) In the final layer filtering, the relative energy difference spectrum, which identifies the singular value of Hankel matrix of real signal, is applied to adaptively remove noise for the secondlayer filtered signal.
Finally, the proposed method is used to effectively filter out noise and extract weak fault features of the simulation and measured signal. Furthermore, the method is also compared with typical filtering methods, respectively. The result indicates that the proposed method demonstrates the optimum performance among the filtering methods for extracting fault feature of bearing.
References
 Zhang L., Xiong G., Liu H., et al. Bearing fault diagnosis using multiscale entropy and adaptive neurofuzzy inference. Expert Systems with Applications, Vol. 37, Issue 8, 2010, p. 60776085. [Publisher]
 Rai V. K., Mohanty A. R. Bearing fault diagnosis using FFT of intrinsic mode functions in HilbertHuang transform. Mechanical Systems and Signal Processing, Vol. 21, Issue 6, 2007, p. 26072615. [Publisher]
 Khadersab A., Shivakumar S. Vibration analysis techniques for rotating machinery and its effect on bearing faults. Procedia Manufacturing, Vol. 20, 2018, p. 247252. [Publisher]
 Deng F., Yang S., Tang G., et al. Self adaptive multiscale morphology AVGHat filter and its application to fault feature extraction for wheel bearing. Measurement Science and Technology, Vol. 28, Issue 4, 2017, p. 045011. [Publisher]
 Bluestein L. I. A linear filtering approach to the computation of discrete Fourier transform. IEEE Transactions on Audio and Electroacoustics, Vol. 18, Issue 4, 1970, p. 451455. [Publisher]
 Aminian M., Aminian F. Neuralnetwork based analogcircuit fault diagnosis using wavelet transform as preprocessor. IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing, Vol. 47, Issue 2, 2000, p. 151156. [Publisher]
 Lou X., Loparo K. A. Bearing fault diagnosis based on wavelet transform and fuzzy inference. Mechanical Systems and Signal Processing, Vol. 18, Issue 5, 2004, p. 10771095. [Publisher]
 Wu J. D., Chen J. C. Continuous wavelet transform technique for fault signal diagnosis of internal combustion engines. NDT & E International, Vol. 39, Issue 4, 2006, p. 304311. [Publisher]
 Huang N. E., Shen Z., Long S. R., et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proceedings Mathematical Physical and Engineering Sciences, Vol. 454, Issue 1971, 1998, p. 903995. [Publisher]
 Yu D., Cheng J., Yang Y. Application of EMD method and Hilbert spectrum to the fault diagnosis of roller bearings. Mechanical Systems and Signal Processing, Vol. 19, Issue 2, 2005, p. 259270. [Publisher]
 Cheng J., Yu D., Tang J., et al. Application of frequency family separation method based upon EMD and local Hilbert energy spectrum method to gear fault diagnosis. Mechanism and Machine Theory, Vol. 43, Issue 6, 2008, p. 712723. [Publisher]
 Qi K., He Z., Zi Y. Cosine windowbased boundary processing method for EMD and its application in rubbing fault diagnosis. Mechanical Systems and Signal Processing, Vol. 21, Issue 7, 2007, p. 27502760. [Publisher]
 Wu Zhaohua, Huang Norden E. Ensemble empirical mode decomposition: a noiseassisted data analysis method. Advances in Adaptive Data Analysis, Vol. 1, Issue 1, 2009, p. 141. [Publisher]
 Yang J., Chen J., Hong R., et al. Multiscale fault frequency extraction method based on EEMD for slewing bearing fault diagnosis. Lecture Notes in Electrical Engineering, Vol. 334, 2015, p. 363370. [Publisher]
 Chen L., Tang G., Zi Y., et al. Improved EEMD applied to rotating machinery fault diagnosis. Applied Mechanics and Materials, Vols. 128129, 2011, p. 154159. [Publisher]
 Torres M. E., Colominas M. A., Schlotthauer G., et al. A complete ensemble empirical mode decomposition with adaptive noise. IEEE International Conference on Acoustics, Speech and Signal Processing, 2011, p. 41444147. [CrossRef]
 Kuai M., Cheng G., Pang Y., et al. Research of planetary gear fault diagnosis based on permutation entropy of CEEMDAN and ANFIS. Sensors, Vol. 18, Issue 3, 2018, p. 782. [Publisher]
 Guo T., Deng Z. An improved EMD method based on the multiobjective optimization and its application to fault feature extraction of rolling bearing. Applied Acoustics, Vol. 127, 2017, p. 4662. [Publisher]
 Cheng J., Yu D., Tang J., et al. Application of SVM and SVD technique based on EMD to the fault diagnosis of the rotating machinery. Shock and Vibration, Vol. 16, Issue 1, 2013, p. 8998. [Publisher]
 Yang Z. X., Zhong J. H. A hybrid EEMDbased SampEn and SVD for acoustic signal processing and fault diagnosis. Entropy, Vol. 18, Issue 4, 2016, p. 112. [Publisher]
 Chen Z., Hu K., Carpena P., Bernaola Galvan P., Stanley H. E., Ivanov P. C. Effect of nonlinear filters on detrended fluctuation analysis. Physical Review E, Vol. 71, 2005, p. 011104. [CrossRef]
 Chen Z., Ivanov P. C., Hu K., Stanley H. E. Effect of nonstationarities on detrended fluctuation analysis. Physical Review E, Vol. 65, 2002, p. 041107. [Publisher]
 Horvatic D., Stanley H. E., Podobnik B. Detrended crosscorrelation analysis for nonstationary time series with periodic trends. Europhysics Letters, Vol. 94, 2011, p. 18007. [Publisher]
 Leistedt S., Dumont M., Lanquart J. P., et al. Characterization of the sleep EEG in acutely depressed men using detrended fluctuation analysis. Clinical Neurophysiology, Vol. 118, Issue 4, 2007, p. 940950. [Publisher]
 Han T., Jiang D., Wang N. The fault feature extraction of rolling bearing based on EMD and difference spectrum of singular value. Shock and Vibration, Vol. 2016, 2016, p. 5957179. [Publisher]
 Buades A., Coll B., Morel J. M. A Review of image denoising algorithms, with a new one. Siam Journal on Multiscale Modeling and Simulation, Vol. 4, Issue 2, 2005, p. 490530. [Publisher]
 Yang J., Fan J., Ai D., et al. Local statistics and nonlocal mean filter for speckle noise reduction in medical ultrasound image. Neurocomputing, Vol. 195, 2016, p. 8895. [Publisher]
 Said A. B., Hadjidj R., Melkemi K. E., et al. Multispectral image denoising with optimized vector nonlocal mean filter. Digital Signal Processing, Vol. 58, 2016, p. 115126. [Publisher]
 Ville D. V. D., Kocher M. SUREbased nonlocal means. IEEE Signal Processing Letters, Vol. 16, Issue 11, 2009, p. 973976. [Publisher]
 Qiao Z., Pan Z. SVD principle analysis and fault diagnosis for bearings based on the correlation coefficient. Measurement Science and Technology, Vol. 26, Issue 8, 2015, p. 085014. [Publisher]
 Jiang H., Chen J., Dong G., et al. Study on Hankel matrixbased SVD and its application in rolling element bearing fault diagnosis. Mechanical Systems and Signal Processing, Vols. 5253, Issue 1, 2015, p. 338359. [Publisher]
 Golafshan R., Sanliturk K. Y. SVD and Hankel matrix based denoising approach for ball bearing fault detection and its assessment using artificial faults. Mechanical Systems and Signal Processing, Vols. 7071, 2016, p. 3650. [Publisher]
 Zhao X., Ye B. Selection of effective singular values using difference spectrum and its application to fault diagnosis of headstock. Mechanical Systems and Signal Processing, Vol. 25, Issue 5, 2011, p. 16171631. [Publisher]
 Loparo K. A. Bearing data Center. Case Western Reserve University, 2000. http://csegroups.case.edu/bearingdatacenter/home. [CrossRef]