Bearing fault diagnosis based on feature extraction of empirical wavelet transform (EWT) and fuzzy logic system (FLS) under variable operating conditions
Fawzi Gougam^{1} , Chemseddine Rahmoune^{2} , Djamel Benazzouz^{3} , Boualem Merainani^{4}
^{1, 2, 3, 4}Solid Mechanics and Systems Laboratory (LMSS), University M’hamed Bougara Boumerdes, Algeria
^{1}Corresponding author
Journal of Vibroengineering, Vol. 21, Issue 6, 2019, p. 16361650.
https://doi.org/10.21595/jve.2019.20092
Received 19 July 2018; received in revised form 20 December 2018; accepted 7 February 2019; published 30 September 2019
JVE Conferences
Condition monitoring of rotating machines has become a more important strategy in structural health monitoring (SHM) research. For fault recognition, the analysis is categorized in two essential main parts: Feature extraction and classification; the first one is used for extracting the information from the signal and the other for decisionmaking based on these features. A higher accuracy is needed for sensitive places to avoid all kinds of damages that can lead to economic losses and it may affect the human safety as well. In this paper, we propose a new hybrid and automatic approach for bearing faults diagnosis. This method uses a combination between Empirical wavelet Transform (EWT) and Fuzzy logic System (FLS), in order to detect and localize the early degradation of bearing state under different working conditions. EWT build a wavelet filter bank to extract amplitude modulatedfrequency modulated component of signal. Modes presenting a high impulsiveness is then selected using the kurtosis indicator. Thereafter, time domain features (TDFs) are applied for the reconstructed signal to extract the fault features which are finally used as an inputs of FLS in order to identify and classify the bearing states. The experimental results shows that the proposed method can accurately extract and classify the bearing fault under variable conditions. Moreover, performance of EWT and empirical mode decomposition (EMD) are studied and shows the superiority of the proposed method.
 We proposed a new method of bearing fault diagnosis based on feature extraction and classification.
 EWT is an efficient decomposition method, it constructs a bank of adaptive filters on the acquired vibration signal.
 The original signal is reconstructed based on the decomposed modes which have more impulsivity of shock in them.
 Time domain features are introduced to extract the fault feature vector on the vibration signal.
 A samples part of the characteristic vector is used to learn the fuzzy system, while the rest of the samples are used for testing the approach.
 The proposed approach gives a satisfactory result with a higher efficiency rate.
Keywords: bearing faults, faults detection, empirical wavelet transform, fuzzy logic, features extraction, fault diagnosis.
1. Introduction
Rotating machines are very useful; they cover a large part in the industry and they are introduced in all mechanical equipment; so their diagnostic and maintenance strategies have an important role for factory economics and also for human safety, by avoiding unexpected accidents that can occur due to the presence and nondetection of defects in machines. Large industrial plants and areas which are sensitive to human safety require high precision, because even the slightest mistake can cause unwanted damages [1, 2]. Most researchers use vibration signals, for their simplicity since it requires just the placement of sensors for signal acquisition; without any complex intervention which is the ideal way to detect failure in any rotating machines caused by their cinematic faults. Machine vibration signatures can be used as an early warning to the operator to make a crucial decision before any serious problem or unwanted downtime. The amplitude of the vibration signature gives an indication of the severity of the problem, while the frequency can indicate the source of the fault [3]. For faults diagnosis, researchers have applied various techniques based on vibration analysis [4, 5] such as: Kurtogram, Cepstrum analysis, power spectrum, which can detect the nature of faults in the rotating machines. These techniques are introduced in frequency domain and all spectral techniques are based on FFT. However, to analyze and process signals in time domain, researchers focused on featurebased methods; the procedure strategy can be divided into two main parts: Feature extraction and classification to make a decision which is based on these features. So, various features time domain features (TDFs) are used such as: root mean square (RMS), standard deviation (SD), kurtosis (KU), variance (VAR), entropy (ENT) and etc. The features extraction actually describes the condition of the machine, they have an important role in classifying defects by introducing a set of statistical features as entries to the classification system for faults selfrecognition [6].
Mechanical faults generate nonstationary signal emerged with back ground of noise in case of an experimental signal acquisition, where the analysis of conventional methods like FFT technique are complex, so the main problem with the Fourier transform is its lack of information in temporal domain. This simply means that if you are able to fault any frequency that appears in the signal, you are unable to determine when it occurs. For more information, researchers applied other techniques such as time frequency analysis which can solve this problem and can detect the abnormality due to the presence of faults. So many methods have been applied for faults diagnosis such as Short Time Fourier Transform (STFT), Wavelet Transform (WT), [79], Wavelet Packet Transform (WPT), Wigner Ville Distribution (WVD) [10], Empirical Mode Decomposition (EMD) [1113], Empirical Wavelet Transform (EWT) [14] etc.
One of the most successful methods in machine faults diagnosis in timefrequency domain is wavelet transform (WT) for its distinct advantages [15]. It was used as a denoising technique for noisy working environment signal processing [16]. However, the high frequency band is not split where the modulation information of machine fault exists. This problem can be overcome by using a wavelet packet transform (WPT) which decomposes the signal in different levels due to their bandpass filters to extract the critical components. These decomposed signals are widely used for features extraction for machine faults diagnosis [17]. For lack of adaptability, WPT is sometimes unable to provide effective information about faults [18, 19]. EMD technique is an adaptive method that decompose a signal into a sum of Intrinsic Mode Function (IMF). However, EMD suffers from the mode mixing, distorted components and the end effects phenomena [20, 21]. By combining the adaptability of EMD and the advantages of WT, Gilles [22], has developed the EWT. Most of researchers introduce the EWT on diagnostic due to their efficiency and accuracy. It has been proved that is better than EMD method for rolling element bearing diagnosis [23, 24].
For automatic fault diagnosis, a classifier is employed to make a quick and reliable decision on the running condition of machinery. The recent works are introduced to give more accuracy to problems; when similar and confusing faults occur in complex systems, leading to difficult identification and localization of the faults. For that raison, a variety of artificial and intelligent techniques have been applied to classify various states. Once the input feature vectors are identified, the decisionmaking rules could be constructed. This could be done through the use of classification tools based on Artificial Intelligence such as: Artificial Neural Network (ANN) [25], Support Vector Machine (SVM) [26, 27] and Fuzzy Logic System (FLS) [28]. The major problem in pattern recognition based on neurons; for example ANN has a system learning like a black box. The exact number of layers is hard to define in order to achieve the best convergence of output. But in opposition, FLS provides a simple way to achieve a definite output (desired) based on vague and ambiguous information using rules for learning system [29].
In this paper, EWT is proposed to decompose the original signal into different modes, after that we calculate kurtosis values of each decomposed modes. As the kurtosis sensitive to impulsiveness of data in one hand and as the bearing fault generates periodic shocks which increases the kurtosis value of signal in the other, only modes having kurtosis value superior than 3 is kept. Thereafter, statistical features (standard deviation (SD) and the Upperbound (UPP)) are applied for the reconstructed signal to extract the information on machine condition. Finally, the obtained feature vectors are considered as inputs data of FLS to classify the various bearing condition states, where the classification is based on theory of fuzzy sets. The system output gave information of various faults types for specific size (0.007 inch) and different operating modes (0 to 3 horsepower). The accuracy of this approach is expressed by analyzing signals of bearing with the small fault size, precise identification and localization of the bearing faults under complex and variable working environment.
This paper is organized as follows; Section 1: Introduction, Section 4: EWT, FLS, Statistical features, experimental study and experiment results, Section 5: Conclusion.
2. Empirical wavelet transform
Combining the adaptability of EMD decomposition and the advantages of WT, Gilles [22] has constructed a wavelet filter bank adaptive to extract relevant frequencies of processed signal. In recent works most of researchers use this method for their performance and their reliability in signal processing and also is used for denoising signal, the Amplitude ModulatedFrequency Modulated (AMFM) components called modes which are the bandpass supports of decomposition signal. filter model mathematic of these modes is [22]:
where $f\left(t\right)$ is complex decomposed signal, ${f}_{k}\left(t\right)$ is frequencymodulated amplitudemodulated signal.
These filters depend on the spectrum of the input signal. First, we note that $f\left(t\right)$ is the input signal and $\widehat{f}\left(w\right)$ it’s Fourier spectrum ,we assume that the Fourier support is in $\left[0,\pi \right]$segmented in $N$, and we denote that ${w}_{i}$to be the limit between each segments ( for $i=$ 0: ${\omega}_{0}$, for $i=\pi $: ${\omega}_{n}$) , the Fig. 1. Corresponds to the segmentation axis of EWT wavelet construction, and we see the transitional phase around ${T}_{n}=\text{2}{\tau}_{n}$, and we define the complex filters for this support segment [0,$\pi $]:
Fig. 1. Fourier axis segmentation and EWT wavelets construction
The empirical wavelets are defined from these equations:
Note that the most used function that satisfies this property:
The detail coefficients ${W}_{f}^{\epsilon}\left(n,t\right)$ are obtained by the inner products of the input signal with the empirical wavelets:
The approximation coefficients with scaling function is given:
Finally, the original signal is decomposed into various empirical modes ${f}_{k}\left(t\right)$, which is given by:
The signal can be reconstructed as follows:
2.1. Fuzzy logic system
Instead of classic logic (Boolean logic) [30, 31], Lotfi Zadeh created a more precise method resting on the intermediary values between {0, 1}, this method is called Fuzzy Logic where this system is based on functions Membership and their rules. Mamdani and Assilian inspired of Mohamed Zadeh this logic [32], by proposing a more widely used model called Mamdani. Fig. 2 shows the structure of FLS and it can be summarized in three stages, executed in the order.
Fig. 2. Fuzzy system structure
1) Fuzzification: is the first step of fuzzy inference process. It is a domain transformation of crisp into fuzzy inputs to determine the truth’s degree for each rule premise.
2) Inference: is the second step in FIS where, for the premise, we calculated the truth value of each rule and applied it to the conclusion part of each rule. As a result, a blurred subset must be assigned to each output variable for each rule. Usually, the most useful inference rules method are minimum and product. In the case of a minimum, the output adhesion function is cut to a height corresponding to the calculated degree of truth of the rule. When selecting the product, the output adhesion function is scaled according to the calculated degree of truth of the rule.in this step a single fuzzy set is obtained for each output variable.
3) Defuzzification: the process of production results considering the fuzzy sets and degrees of membership, this step covert the crisp value which is obtained on fuzzy output to crisp values where it is performed with one of many defuzzification methods.
Instead of Neural Network learning, the use of rule has an essential part in Fuzzy Inference System which are introduced as follows:
Rule: if $x$ is $A$ and $y$ is $B$ than $z$ is $C$;
While: $x$ is $A$, $y$ is $B$: are the premise of our rule
$Z$ is $C$: Consequence rule;
And: Operator.
2.2. Basic elements of fuzzy logic (linguistic variables)
The essential notion of linguistic variable was introduced by Zadeh, it immediately suggests that the values of this variable are not numerical, but rather symbolic, in terms of words or expressions of natural language. A linguistic variable is characterized by a triplet ($x$, $T\left(x\right)$, $u$) where: $x$: the variable name; $T\left(x\right)$: variable term set, it is the set of names of the linguistic values of $x$ of which each value is a fuzzy subset defined in $U$ which is the universe of discourse.
2.3. Fuzzy sets and membership functions
Classical logic has two acceptable values 0 or 1, fuzzy logic intervenes to define intermediary of these two values based on the membership sets. Function of membership are defined as follows:
There are different forms of membership functions; trapezoidal and triangular are the most used ones:
2.4. Operators of fuzzy logic
2.5. Operator AND
The fuzzy AND operator is defined as the intersection of two fuzzy sets A and B:
2.6. Operator OR:
The fuzzy OR operator is defined as the union of two fuzzy sets A and B, respectively being the largest fuzzy set contained in A and B:
2.7. Gravity center defuzzification
The center of gravity (GC) method is based on the calculation of the abscissa gravity center of resulting surface (Eq. (14)) obtained on the fuzzy subset of the solution determined by the aggregation of fuzzy rules. Fig. 3. shows the procedure solution of fuzzy logic system using single input $y$ and the solution $x$ which is given by the method of the center of gravity:
Fig. 3. Fuzzy system structure
3. Statistical features
These statistical features are calculated to extract information of signal from various states of bearing, were these features are explained as follow [33]:
4. Experimental study
4.1. Experimental description
The signals vibration data are provided by Case Western Reserve University Bearing Data Center [34], the Fig. 4. shows the test rig of bearing in induction motor of 1.5 kw in the left, a dynamometer (in the right), a torque transducer/Encoder (center), the faults diameter is (7, 14, 21, 28) mils, (1 mil = 0.001 inch) the tested bearing is of type (62052RS JEM SKF) as it is presented in Fig. 4, knowing that we have two bearing in the motor FanEnd (FE) and the Driveend (DE).
Fig. 4. Test band
Single point defects outer race (ORF), Innerrace (IRF), ball (BF) was introduced separately to the test of bearing using electrodischarge machining (EDM) with different sizes. The vibration data was recorded for motor Horse power (1 HP = 1772 RPM) and for sampling frequency of 12 kHz. We know that the bearing faults produce a repetitive mechanical shockwaves over time, which will create instantaneous frequencies specific to each fault of these defects as shown in Table 1.
Table 1. Formula of expected bearing faults frequencies
Type of defect

Formula fault

Innerrace fault frequency

${f}_{i}=\frac{{N}_{b}}{2}{f}_{s}\left(1+\frac{{D}_{b}}{{D}_{c}}\mathrm{c}\mathrm{o}\mathrm{s}\alpha \right)$

Outerrace fault frequency

${f}_{o}=\frac{{N}_{b}}{2}{f}_{s}\left(1\frac{{D}_{b}}{{D}_{c}}\mathrm{c}\mathrm{o}\mathrm{s}\alpha \right)$

Rollingelement fault frequency

${f}_{b}=\frac{{D}_{c}}{{D}_{b}}{f}_{s}\left(1\frac{{{D}_{b}}^{2}}{{{D}_{c}}^{2}}{\mathrm{c}\mathrm{o}\mathrm{s}}^{2}\alpha \right)$

Cage fault frequency

${f}_{c}=\frac{{f}_{s}}{2}\left(1\frac{{D}_{b}}{{D}_{c}}\mathrm{c}\mathrm{o}\mathrm{s}\alpha \right)$

4.2. Proposed study of rolling bearing
The vibration signals presented in Fig. 5 are processed in order to classify the various bearing faults of 0.007 inch diameter, by using the procedure of mentioned flowchart (Fig. 6), the steps are described as follows:
1) Have an interesting database on the real signals of these two states (healthy and defective) in order to enrich the input data. The Empirical Wavelet Transform is used to decompose the original signal into several bandpass frequencies (modes).
Fig. 5. Vibration signals of various bearing fault
2) After decomposition, the signal is reconstructed from the significant modes by calculating the kurtosis value, this indicator has a direct relation with fault impulsivity because it causes repetitive shocks during signal acquisition and it increases the amplitude of vibration signature where the fault characteristic frequency depends on the faultlocalisation in bearing.
3) Features extraction and selection: the statistical features or time domain features (TDFs) are applied for the reconstructed signal to extract the health condition of bearing. After that, the selected features from their high performance of separability and stability are used as inputs; to avoid all kind of confusion and similiarity of fault states under various working conditions (0 to 3 HP). The Fig. 8. shows variation of feature samples from different states.
4) Features classification: The selected Features will be defined as the input vectors of fuzzy classifier (FLS) to classify the various faults of bearing under different working conditions.
Fig. 6. Flowchart of proposed method
Table 2. Description of bearing data set
Fault type

Data size vector

Diameter defect

Condition working

Healthy (H)

112

//

[0, 1, 2, 3] hp

Inner race (IR)

112

0.007 inch

[0, 1, 2, 3] hp

Outer race (OR)

112

0.007 inch

[0, 1, 2, 3] hp

Ball (B)

112

0.007 inch

[0, 1, 2, 3] hp

The corresponding modes shown in Fig. 7. contain the relevant information on bearing condition, we calculate time domain features (TDFs) on the reconstructed signal, knowing that the two indicators of standard deviation (SD) and upperbound (UPP) have been selected according to the separability of states and especially there is no interesting confusion between them and also the corresponding features values have a great stability; this avoids any kind of confusion states.
Fig. 7. Decomposed signal using EWT
Fig. 8. Feature extraction of 0.007 inch of various defects
a)
b)
c)
d)
For this study, a part of initial vector 35.71 % was used for learning system and the rest for testing. The minimum and the maximum has been extracted from the selected features as the Fig. 9 and Fig. 10. show, the interval minmax of features states are the inputs of FLS. The membership of input and output are shown in Fig. 11, Fig. 12 respectively.
When the fuzzification part is completed, the laws must be made for the defuzzification where this part is essential for classification. The result of the fuzzy output depends on the laws implemented rules for fuzzy learning system.
Rules are used for training fuzzy system that is based on the fuzzified memberships of inputs, the (SD, UPP) features have been introduced as two inputs for classifier to enhance and confirm the fuzzified entries (operator and) and eliminating the undefined values from a signal (operator OR) using the following rules:
• If: SD is H and UPP is H the output is H,
• If: SD is IR and UPP is IR the output is IR,
• If: SD is OR and UPP is OR the output is OR,
• If: SD is B and UPP is B the output is B,
• If: SD is H or UPP is H the output is H,
• If: SD is IR or UPP is IR the output is IR,
• If: SD is OR or UPP is OR the output is OR,
• If: SD is B or UPP is B the output is B.
Fig. 9. Minima and maxima of standard deviation feature
Fig. 10. Minima and maxima of upperbound feature
For testing this approach, 64.29 % of the input vector is introduced for the FLS classifier, in order to show the effectiveness of this method, we will test the fuzzy system with new samples which are not introduced in the learning system. The output correspond to forth classes of bearing states under variable working condition. The data contains various faults (IR ,OR,B) of 0.007 inch under different working conditions (0 to 3 HP) as it is represented in Table 4, where classification result is given in the Fig. 13.
Fig. 13 shows the classification of the proposed approach; all samples are classified with some fluctuations that are very closer to the desired values. Based on this result, we conclude that this classification result gives a good performance and satisfactory detection and identification for bearing fault.
To confirm the proposed diagnosis strategy, we will compare the given approach with other decomposition method using the same classifier, inputs and rules. First, we give the classification without EWT decomposition using SDUPP (standard deviationupperbound) features as inputs, the result is presented in the Fig. 14. From the last figure, we can see a correct classification of H samples and acceptable classification in B samples with some fluctuation in ORIR samples.
Table 3. Gravity center of output
Bearing states

Healthy

IRfault

ORfault

Ballfault

Working condition (HP)

0 1 2 3

0 1 2 3

0 1 2 3

0 1 2 3

Length of training data

40

40

40

40

Desired output

Gravity center of H

Gravity center of IR

Gravity center of OR

Gravity center of B

Value of desired classes

0.09808

0.3333

0.6642

0.9019

Fig. 11. Membership functions of inputs fuzzy logic
a) First input fuzzification
b) Second input fuzzification
Fig. 12. Membership functions of output fuzzy logic
Fig. 15 presents a classification of SDUPP features using EMD decomposition. We observed a large fluctuation of H state in 1797 rpm (HP0) which can lead to unsatisfactory classification. Generally, excepting Healthy samples, the classification is good and acceptable.
On the basis of compared results, the proposed approach gives a satisfactory classification to early fault detection and identification of rolling bearing under various working condition.
Fig. 13. Classification result of proposed method
Fig. 14. Classification result of proposed method without EWT
Fig. 15. Classification result of proposed method with EMD
5. Conclusions
Early detection of bearing faults is one of important research fields. It permits to reduce maintenance cost and avoiding economic losses and human risks. To monitor the condition of bearing for various condition modes and different fault types in noisy environment, which that constitutes a new challenge on rolling bearing diagnosis. For the complexity process and the reason of effectiveness, a significant research has been applied to health monitoring with emergence a variety of methods and techniques in signal processing. In this paper we present a Fuzzy Logic System decision using EWT and TDFs. The procedure is based on feature extraction and classification for obtaining the current health of rotating component. First, the vibration signal is decomposed into many modes to reconstruct new vibration signal. After that, SDUPP features are applied for the reconstructed signal and the obtained feature vectors will be defined as inputs of FLS; for an automatic and intelligent fault detection and identification. This proposed method was applied on real experimental data where the results show the effectiveness of the approach for bearing fault diagnosis under various working condition.
References
 Miao, Wang Q., Pecht D. A probabilistic description scheme for rotating machinery health evaluation. Journal of Mechanical Science and Technology, Vol. 24, Issue 12, 2010, p. 24212430. [Publisher]
 Yang Z., Chen X., Xie Y., et al. The hybrid multivariate analysis method for damage detection. Structural Control and Health Monitoring, Vol. 23, Issue 1, 2016, p. 123143. [Publisher]
 Goyal D., Pabla B. S. The vibration monitoring methods and signal processing techniques for structural health monitoring: a review. Archives of Computational Methods in Engineering, Vol. 23, Issue 4, 2016, p. 585594. [Publisher]
 Lei Y., He Z., et al. Application of an improved kurtogram method for fault diagnosis of rolling element bearings. Mechanical Systems and Signal Processing, Vol. 25, Issue 5, 2011, p. 17381749. [Publisher]
 Wang D., Guo W., Wang X. A joint sparse wavelet coefficient extraction and adaptive noise reduction method in recovery of weak bearing fault features from multicomponent signal mixture. Applied Soft Computin, Vol. 13, Issue 10, 2013, p. 40974104. [Publisher]
 Ali Jaouher Ben, Fnaiech Nader, Saidi Lotfi, et al. Application of empirical mode decomposition and artificial neural network for automatic bearing fault diagnosis based on vibration signals. Applied Acoustics, Vol. 89, 2015, p. 1627. [Publisher]
 Yan R., Gao R., Chen X. Wavelets for fault diagnosis of rotary machines: a review with applications. Signal Processing, Vol. 96, 2014, p. 115. [Publisher]
 Mallat S. A Wavelet Tour of Signal Processing. 3d Edition, The Sparse Way, 2009, p. 832. [CrossRef]
 Combastel C., Lesecq S., Petropol S., Gentil S. Modelbased and wavelet approaches to induction motor online fault detection. Control Engineering Practice, Vol. 10, Issue 5, 2002, p. 493509. [Publisher]
 Blodt M., Bonacci D., Regnier J., Chabert M. Online monitoring of mechanical faults in variablespeed induction motor drives using Wigner distribution. IEEE Transactions on Industrial Electronics, Vol. 55, Issue 2, 2008, p. 522533. [Publisher]
 Lei Y., Lin J., He Z., et al. A review on empirical mode decomposition in fault diagnosis of rotating machinery. Mechanical Systems and Signal Processing, Vol. 35, Issues 12, 2013, p. 108126. [Publisher]
 Gao Q., Duan C., Fan H. Rotating machine fault diagnosis using empirical mode decomposition. Mechanical Systems and Signal Processing, Vol. 22, Issue 5, 2008, p. 10721081. [Publisher]
 Li Y., Xu M., Wei Y., et al. An improvement EMD method based on the optimized rational Hermite interpolation approach and its application to gear fault diagnosis. Measurement, Vol. 63, Issue 2, 2015, p. 330345. [Publisher]
 Merainani B., Benazzouz D., Rahmoune C. Early detection of tooth crack damage in gearbox using empirical wavelet transform combined by Hilbert transform. Journal of Vibration and Control, Vol. 23, Issue 10, 2017, p. 16231634. [Publisher]
 Kankar P., Sharma S., And Harsha S. Fault diagnosis of ball bearings using continuous wavelet transform. Applied Soft Computing, Vol. 11, Issue 2, 2011, p. 23002312. [Publisher]
 Kulkarni G., Sahasrabudhe D. Application of wavelet transform for fault diagnosis of rolling element bearings. International Journal of Technology Enhancements and Emerging Engineering Research, Vol. 2, Issue 4, 2013, p. 138148. [CrossRef]
 Xu Q., Lia Z. Recognition of wear mode using multivariable synthesis approach based on wavelet packet and improved threeline method. Mechanical Systems and Signal Processing, Vol. 21, Issue 8, 2007, p. 31463166. [Publisher]
 Chen S. W., Chen Y. H. Hardware design and implementation of a wavelet denoising procedure for medical signal processing. Sensors, Vol. 15, 2015, p. 2639626414. [Publisher]
 Zhong J. H., Wong P. K., Yang Z. X. Simultaneousfault diagnosis of gearboxes using probabilistic committee machine. Sensors, Vol. 16, Issue 2, 2016, p. 185. [Publisher]
 Xiao Q., Li J., Bai Z., Sun J., Zhou N., Zeng Z. A small leak detection method based on VMD adaptive denoising and ambiguity correlation classification intended for natural gas pipelines. Sensors, Vol. 16, Issue 12, 2009, p. 2116. [Publisher]
 Wang T., Zhang M., Yu Q., Zhang H. Comparing the applications of EMD and EEMD on timefrequency analysis of seismic signal. Journal of Applied Geophysics, Vol. 83, 2012, p. 2934. [Publisher]
 Gilles J. Empirical wavelet transform. IEEE Transactions on Signal Processing, Vol. 61, Issue 16, 2013, p. 3999401. [Publisher]
 Maheshwari S., Pachori R. B., Acharya U. R. Automated diagnosis of glaucoma using empirical wavelet transform and correntropy features extracted from fundus images. EEE Journal of Biomedical and Health Informatics, Vol. 21, Issue 3, 2016, p. 803813. [Publisher]
 Jiang Y., Zhu H., Li Z. A new compound faults detection method for rolling bearings based on empirical wavelet transform and chaotic oscillator. Chaos Soliton Fractals, Vol. 89, 2016, p. 819. [Publisher]
 Wang C., Yuan K., Shen P., et al. Applications of fault diagnosis in rotating machinery by using time series analysis with neural network. Expert Systems with Applications, Vol. 37, Issue 2, 2010, p. 16961702. [Publisher]
 FernándezFrancos D., Martinez D., FontenlaRomero O., et al. Automatic bearing fault diagnosis based on oneclass vSVM. Computers and Industrial Engineering, Vol. 64, Issue 1, 2013, p. 357365. [Publisher]
 Wyk B., Wyk M., Qi G. Difference histograms: a new tool for time series analysis applied to bearing fault diagnosis. Pattern Recognition Letters, Vol. 30, 2013, p. 595599. [CrossRef]
 Bezdek J. C. Pattern Recognition with Fuzzy Objective Function Algorithms. Springer, 2008. [CrossRef]
 Krishnakumari A., Elayaperumal A., Saravanan M., Arvindan C. Fault diagnostics of spur gear using decision tree and fuzzy classifier. The International Journal of Advanced Manufacturing Technology, Vol. 89, Issues 912, 2017, p. 34873494. [Publisher]
 Hemba S., Islam N. Fuzzy logic: a review. International Journal of Computer Sciences and Engineering, Vol. 5, Issue 2, 2017, p. 6163. [CrossRef]
 Kasbe T., Pippal R. S. Dengue fever: stateoftheart symptoms and diagnosis. International Journal of Computer Sciences and Engineering, Vol. 4, Issue 6, 2016, p. 2630. [CrossRef]
 Dadios E. P. Fuzzy Logic Controls, Concepts, Theories and Applications. IntechOpen, 2012, p. 325350. [Publisher]
 Javed Kamran, Gouriveau Rafael, Zerhouni Noureddine, Nectoux Patrick A feature extraction procedure based on trigonometric functions and cumulative descriptors to enhance prognostics modeling. IEEE International Conference on Prognostics and Health Management, United States, 2013. [CrossRef]
 Case Western Reserve University Bearing Data Center, http://csegroups.case.edu/bearingdatacenter/home. [CrossRef]