MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2019/3089417 3089417 Research Article A Rolling Bearing Fault Diagnosis Method Based on EMD and Quantile Permutation Entropy https://orcid.org/0000-0002-0102-9820 Chen Qiang-qiang https://orcid.org/0000-0003-4841-3915 Dai Shao-wu Dai Hong-de Kratz Frederic Naval Aviation University Yantai 264000 China 2019 492019 2019 07 05 2019 04 07 2019 01 08 2019 492019 2019 Copyright © 2019 Qiang-qiang Chen et al. 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.

The vibration signals resulting from rolling bearings are nonlinear and nonstationary, and an approach for the fault diagnosis of rolling bearings using the quantile permutation entropy and EMD (empirical mode decomposition) is proposed. Firstly, the EMD is used to decompose the rolling bearings vibration signal, and several IMFs (intrinsic mode functions) spanning different scales are obtained. Secondly, aiming at the shortcomings of the permutation entropy algorithm, a new permutation entropy algorithm based on sample quantile is proposed, and the quantile permutation entropy of the first few IMFs, which contain the main fault information, is calculated. The quantile permutation entropies are accordingly seen as the characteristic vector and then input to the particle swarm optimization and support vector machine. Finally, the proposed method is applied to the experimental data. The analysis results show that the proposed approach can effectively achieve fault diagnosis of rolling bearings.

National Natural Science Foundation of China ZR2017MF036 Defense Science and Technology Project Foundation of China F062102009
1. Introduction

Rolling bearings are the most frequently used component in a rotary machines, which play an essential role in the modern industry . Rolling bearings failures could lead to uncontrollable losses for industrial equipment . Therefore, rolling bearings fault diagnosis has appealed a significant attention from the research and engineering communities. Rolling bearings fault diagnosis can be summarized in three steps: signal collection, feature extraction, and fault pattern classification . When the rolling bearings work with faults, their fault information always appears to be uncertain, and the vibration signals will also be nonlinear and nonstationary. Therefore, how to extract the operate state from the vibration signals is the key of rolling bearing fault diagnosis . With the application of the classical nonlinear algorithm, a lot of theory, such as chaos, on fractal dimension has been applied to the feature extraction field. Correlation dimension and Naïve Bayes are combined and used to bearing fault diagnosis by Zhang et al. . Chaotic self-synchronization was applied to ball bearing signals by Yau et al. . All of the classical nonlinear theory have made some accomplishments to fault diagnosis. However, these algorithms often require a lot of time, and it is vulnerable to noise interference in the application process.

The application of entropy theory provides a new idea for rolling bearing fault diagnosis. Pincus proposed the concept of approximate entropy (ApE) to measure the probability of new patterns in signals according to the characteristics of time series . However, ApE depends heavily on the record length, and its estimated value is uniformly lower than the expected one. In order to overcome the shortcomings of ApE, Richman and Moorman put forward a new form of entropy, sample entropy (SampE) . Nonetheless, the definition of SampE must include a template matching; otherwise, it is meaningless. A new kind of entropy named permutation entropy (PE) was put forward by Bandt and Pompe in 2002 . PE is a new method proposed for detecting the randomicity and dynamic changes of time series, which can be used in the field of fault diagnosis. PE and support vector machine (SVM) are combined and used to bearing fault diagnosis by Zhang et al. , and the superiority of PE is proved.

However, due to the complexity of mechanical systems, the randomicity and dynamic changes of the vibration signal behave on different scales, making it necessary to analyze the vibration signal with PE in a multiscale way [11, 12]. The commonly used signal decomposing method is empirical mode decomposition (EMD), and EMD is a kind of data-driven and adaptive nonlinear signal decomposition method and hence suitable for processing nonlinear signals of rolling bearing . In the calculation of PE, the selection of embedding dimension affects the performance of PE to a certain extend. In order to adaptively complete the embedding dimension selection problem and improve the feature extraction ability of PE, a new kind of entropy named sample quantile PE (SQPE) is defined by combining sample quantile and PE in this paper. Firstly, the EMD method is used to decompose the vibration signal, and several IMFs spanning different scales are obtained. Secondly, the SQPE of the first few IMF components, which contain the main fault information, is calculated. After extracting feature vectors by EMD-SQPE, the SVM which was optimized by particle swarm optimization (PSO) is used as a classifier. PSO enhances the adaptability and feasibility of SVM and improves the recognition ability of SVM.

Compared with the traditional feature extraction method, this paper solves the problem of embedding dimension selection in the PE calculation process by SQPE entropy calculation and improves the feature extraction ability of PE at the same scale. The proposed method is applied to the bearing experimental data. The results show that the method can effectively distinguish the fault type of rolling bearings and is an effective fault diagnosis method.

The remainder of this paper is organized as follows: Section 2 proposes a concept of SQPE based on PE. In Section 3, the proposed method based on EMD-SQPE-PSOSVM is introduced. In Section 4, several examples are presented to demonstrate the effectiveness of the proposed algorithm, and a conclusion is given in Section 5.

2. The Proposed Feature Extraction of EMD-SQPE

In order to realize the feature extraction of different vibration signals of rolling bearings, the multiscale decomposition method is first used to decompose the signals into different time scales. Then, using the SQPE proposed in this paper, the SQPE value of the IMF component containing the main fault information is calculated to deeply mine the rolling bearing fault information.

2.1. Empirical Mode Decomposition

EMD was put forward by Huang N. E. in 1998, which seemed much more promising and had appealed a lot of attention . The EMD algorithm can decompose a nonlinear signal self-adaptively into some IMFs and a residual series :(1)st=i=1nimfit+rnt,where st is the original time series, imfit is the ith IMF, and rnt is the residual series that represents the trend of the signal. The IMF component is a signal that has zero mean and a number of maxima and minima that are equal or differs by at most one.

The EMD algorithm can decompose a nonlinear and nonstationary signal self-adaptively into some IMFs. Each IMF contains fault information of different frequency bands and different scales of the original vibration signal .

2.2. Sample Quantile Permutation Entropy (SQPE)

The SQPE proposed in this paper is an improved algorithm based on PE, which mainly solves the problem of PE in embedding dimension selection and improves the feature extraction ability of PE for complex signals through the quantile principle. The definition of SQPE is based on PE, and PE is described as follows.

Consider a time series, xi,i=1,2,...,N, with length N, which can be reconstructed as (2)X1=x1,x1+τ,,x1+m1τ,Xj=xj,xj+τ,,xj+m1τ,Xk=xk,xk+τ,,xk+m1τ,where m is the embedding dimension and τ is the time delay. The element of Xj can be rearranged in an increasing order as follows:(3)xi+j11τxi+jm1τ,where if there exist two elements in Xj that have the same value, such as, xi+ji11τ=xi+ji21τ, then we order the quantities x according to the values of their corresponding j; namely, if ji1<ji2, then xi+ji11τxi+ji21τ is written. Accordingly, any vector Xj can be mapped onto a group of symbols as(4)Sl=j1,j2,,jm,where l=1,2,,k, km!. m! is the largest number of distinct symbols, and Sl is one of the m! permutations of m distinct symbols, which is mapped onto the m number symbols in m-dimensional embedding space. Then, the PE for the time series can be defined as the Shannon entropy for (5)Hpm=j=1kPjlnPj.

It is noticed that Hpm attains the maximum value, lnm!, when Pj=1/m!. For convenience, Hpm can be normalized by lnm! as(6)Hp=Hpmlnm!.

Obviously, the values of Hp is 0,1. A bigger value of Hp indicates that the signal is much more complex and the maximum value of Hp=1 means that the signal is very nonlinear and nonstationary. Therefore, PE is a very suitable tool for describing the local order structure and amplifying the dynamic changes of signal.

There are two parameters to be set in the calculation of PE: the embedding dimension m and time delay τ. The time delay τ has a few effect on the calculation of PE; that is, there are very small differences among the PEs between different time delay τ. So we set τ=1. The embedding dimension m has a big effect on the estimation of PE, and there are differences among the PEs between different embedding dimension m. There is no definite standard for the choice of embedding dimension m. Bandt recommended m=27, and this recommendation is not adaptive because different m values will effect the fault diagnosis results to some extent.

In order to improve the adaptability of the selection with embedding dimension and improve the feature extraction ability of PE. On the basis of PE, this paper introduces the concept of sample quantile and defines SQPE. Sample quantile divides the probability distribution range of a random variable into equal numerical points, without considering the time dimension characteristics of the time series, and directly and effectively characterizes the characteristics of the fluctuation of the signal time series data .

Given a time series X1,X2,,Xn, define the order statistic as X1X2Xn. The p quantile of the time series can be defined as follows :(7)mp=Xk,k=np+1,where represents integer operational symbols. p is often set to 0.25, 0.5, and 0.75. In this paper, we set p=0.75.

The sample quantile has the advantages of simple calculation and small calculation. And the sample quantile is a nonparametric statistic, which can accurately reflect the aggregation characteristics of the data at a certain quantile in the absence of the overall distribution prior information. The SQPE is defined as follows:

Given a time series xi, the PE is calculated under embedding dimension m=2,3,,7 and τ=1, and their PE values are denoted by PE=PEm=2,PEm=3,,PEm=7

The 0.75 sample quantile of PE=PEm=2,PEm=3,,PEm=7 is calculated, which is denoted by SQPE=quantile PEp=0.75

SQPE is used instead of PE as the feature information for fault diagnosis

The proposed SQPE algorithm can be described briefly as in Figure 1.

Flow chart of SQPE.

In the calculation of SQPE, the problem of embedding dimension selection is not considered, but the entropy value is solved by calculating the sample quantile of all embedding dimensions. This avoids the effect of embedding dimension on feature extraction to some extent.

3. The Proposed Fault Diagnosis Based on SQPE

In addition to EMD-SQPE, the proposed fault diagnosis also adopts the PSO-SVM algorithm. The SVM was originally a deterministic algorithm for finding the linear separating hyperplane of a binary labeled dataset. Similar with the neural network, it has machine learning features. The core algorithm of SVM is to change the original data space to a high-dimensional feature space through a nonlinear mapping . The ability of SVM classification depends largely on the kernel function parameter and penalty parameter. The PSO has the following advantages over other techniques. PSO has the characteristics of stable convergence, and it can get a high-quality solution within shorter calculation time. PSO was applied into adaptive selection of the best penalty parameter and kernel function parameter with K-fold cross validation error minimum as the optimization goal .

Since the introduction of the PSO algorithm in 1995, researchers have put much effort to improve the original version of PSO. Shi and Eberhart proposed the idea of inertia weight in order to balance the local and global search during the optimization process, in which the inertia weight is linearly decreasing over iterations . The performance of SVM can be improved by optimizing the parameters of SVM through PSO.

Based on the advantages of EMD, SQPE, and PSO-SVM, the proposed rolling bearing fault diagnosis algorithm is described as follows.

Decomposing the initial rolling bearing vibration signals into several IMF components using the EMD method.

SQPE of the first several IMF components which contain the main failure information is calculated and taken as the feature vector. According to , the main fault information is often contained in the first few IMFs. The cross correlation coefficients of several IMFs with the original signal are relatively large and retain more information of the original signal. Therefore, the IMFs having large cross correlation coefficients are used as the main components. The cross correlation can be defined as

(8)ρxy=Cxyσxσy,

where Cxy is the cross-covariance function. The cross correlation coefficients ρxy indicate the correlation between signals; the larger the value of ρxy, the stronger the correlation between two signals.

The feature vectors are seen as the inputs to the SVM classifier. Here, the feature information in (2) is divided into two parts: one is for training, and the other is for testing.

PSO is applied into adaptive selection of the best penalty parameter and kernel function parameter.

The output results of testing samples achieve to discriminate the fault categories automatically.

The proposed EMD-SQPE-PSOSVM fault diagnosis algorithm can be described briefly as in Figure 2.

Flow chart of the proposed method.

4. Simulation Results 4.1. Experimental Data

In order to validate the capability of the EMD-SQPE-PSOSVM method, experimental analyses on rolling bearing faults were conducted. The 6205-2RS JEM SKF deep groove ball bearing was used in the experimental, and single point faults were introduced to the test rolling bearings using electrodischarge machining with fault diameters of 0.3556 mm . The shaft rotating speeds of the bearing is 1772 r/min, and the sampling frequency is 12000 Hz. Bearing fault under consideration covers outer race fault (ORF), inner race fault (IRF), rolling element fault (REF), and the norm (NORM), and the length of the bearing vibration signal series is 2048. There are 80 vibration signals for each category and totally 320 chosen randomly from the data sets. In test 1∼3, among these 80 data, 40 samples are randomly selected as training data, and the rest 40 will be seen as testing data.

In the experiments, the feature vectors are seen as the inputs to the PSOSVM model, and the output of the PSOSVM model is the label value corresponding to each fault state, set to NORM (label1), IRF (label2), REF (label3), and ORF (label4).

4.2. Results 4.2.1. Test 1: PE(SQPE)-SVM and PE(SQPE)-PSOSVM

PE of the data sets is calculated and taken as the feature vector. The embedding dimension and the delay of PE were set to m=6 and τ=1. 40 samples of every category are used to train the SVM and PSOSVM classifiers. The classification results of the desired outputs and the outputs of the trained SVM are shown in Figure 3(a). And the results of PSOSVM are shown in Figure 3(b).

The output results and the desired outputs of testing sets with feature vector consisting of the original data’s PE: (a) SVM diagnosis results; (b) PSOSVM diagnosis results.

The accuracy of 160 testing sets with the feature vector consisting of PE by the SVM classifier is 84.375% (135/160), and the accuracy of the PSOSVM classifier is 87.5% (140/160). The best penalty parameter is 1.8232, and the best kernel function parameter is 645.2294 after PSO optimized.

In order to visually explain the fault diagnosis results of PE and SQPE, SQPE of the data sets is calculated and taken as the feature vector. The accuracy of 160 testing sets with the feature vector consisting of SQPE by the SVM classifier is 98.125% (157/160), and the accuracy of the PSOSVM classifier is 98.75% (158/160). The best penalty parameter is 0.1, and the best kernel function parameter is 2796.8326 after PSO optimized.

4.2.2. Test 2: EMD-PE-SVM and EMD-PE-PSOSVM

As shown in Figure 3, the fault categories cannot be classified from each other. So it is necessary to preprocess the original vibration signals by using EMD. The original signal is decomposed into several IMFs by using the EMD method, respectively. Since the main fault feature is often contained in the first few IMFs, the PE of first five IMFs are calculated and seen as the feature vector and used to train and test the SVM and PSOSVM models. The embedding dimension and the delay of PE were set to m=6 and τ=1. The classification results of the desired outputs and the outputs of the trained SVM are shown in Figure 4(a). And the results of PSOSVM are shown in Figure 4(b).

The output results and the desired outputs of testing sets with feature vector consisting of the EMD-PE: (a) SVM diagnosis results; (b) PSOSVM diagnosis results.

The accuracy of 160 testing sets with the feature vector consisting of EMD-PE by the SVM classifier is 95% (152/160), and the accuracy of the PSOSVM classifier is 96.875% (155/160). The best penalty parameter is 13.2931, and the best kernel function parameter is 1960.7666 after PSO optimized.

4.2.3. Test 3: EMD-SQPE-SVM and EMD-SQPE-PSOSVM

As shown in Figure 4, the accuracy of EMD-PE is obviously higher than PE, but it still cannot be fully classified. Using the EMD-SQPE feature extraction method presented in this paper, the original signal is decomposed into several IMFs by using the EMD method, and the SQPE of first five IMFs is calculated and seen as the feature vector and used to train and test the SVM and PSOSVM model. The delay of SQPE was set to τ=1, and the embedding dimension was set to m=2,3,7. The classification results of the desired outputs and the outputs of the trained SVM are shown in Figure 5(a). And the results of PSOSVM are shown in Figure 5(b).

The output results and the desired outputs of testing sets with the feature vector consisting of the EMD-SQPE: (a) SVM diagnosis results; (b) PSOSVM diagnosis results.

The accuracy of 160 testing sets with the feature vector consisting of EMD-SQPE by the SVM classifier is 98.75% (158/160), and the accuracy of the PSOSVM classifier is 100% (160/160). The best penalty parameter is 0.1, and the best kernel function parameter is 1307.5593 after PSO optimized.

In order to visually explain the fault diagnosis results on the vibration signal of the rolling bearing, the results obtained by each method are shown in Table 1. The numerical values in the table indicate the correct classification results obtained by the classification algorithm for the 40 test set samples in each category.

Accuracy rate comparison of three algorithms.

Fault diagnosis method NORM IRF REF ORF Accuracy (%)
PE-SVM 40 39 35 21 84.375
PE-PSOSVM 40 37 38 25 87.5
EMD-PE-SVM 40 38 40 34 95
EMD-PE-PSOSVM 40 39 40 36 96.875
EMD-SQPE-SVM 40 40 40 38 98.75
EMD-SQPE-PSOSVM 40 40 40 40 100

As shown in Table 1, optimized SVM by PSO can use adaptive selection of the best penalty parameter and kernel function parameter, which improved the accuracy rate of fault diagnosis. The fault diagnosis method based on EMD-PE can analyze the vibration signal with PE in a multiscale way, which improved the accuracy than PE. The fault diagnosis algorithm proposed in this paper has the best effect, and the recognition rate reaches 100%. It proves the superiority of the SQPE algorithm in fault feature extraction and also proves the effectiveness of the algorithm proposed in this paper.

4.2.4. Test 4: Different Percentages of the Samples Used for Training

To illustrate the influence of the number of training sets, the experiments were designed by different training set scales (10%, 20%, 30%, and 40% of total data sets), and the remaining sets are used for test; the accuracies of the EMD-SQPE-PSOSVM and EMD-PE-PSOSVM for different percentages of the samples used for training are presented in Table 2.

Accuracies with different training set scales.

Fault diagnosis method Percentage of the samples used for training
10% 20% 30% 40% 50%
EMD-PE-PSOSVM 97.9167% 96.875% 97.7679% 97.3958% 96.875%
EMD-SQPE-PSOSVM 100% 100% 100% 100% 100%

As shown in Table 2, the accuracy of EMD-PE-PSOSVM is not good enough to classify different bearing faults than EMD-SQPE-PSOSVM with different percentages of the samples used for training. Even with less training data sets, the EMD-SQPE-PSOSVM still has 100% accuracy.

5. Conclusion

In this paper, a new rolling bearing fault diagnosis method has been introduced based on the vibration signal analysis using EMD-SQPE and PSOSVM. The test data analysis conclusions illustrate that the proposed algorithm can not only analyze among different fault categories but also identify the level of fault severity. The experimental results show the following:

SQPE is a valid method to measure the complexity of nonlinear and nonstationary vibration signals. Compared with PE, SQPE can extract the features with high distinguishability.

SQPE can avoid the selection of embedding dimension, and it is not necessary to determine the specific embedding dimension in the solution process, which has better adaptability.

EMD has a good application prospect in analyzing the nonlinear signal. After the decomposition of the original signal by EMD, entropy is applied to measure the feature of the IMFs obviously.

PSO can select the best penalty parameter and kernel function parameter adaptively, which improved the performance of SVM.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the Shandong Natural Science Foundation of China (grant number ZR2017MF036) and Defense Science and Technology Project Foundation of China (grant number F062102009).

Wang H. Ke Y. Song L. Tang G. Chen P. A sparsity-promoted decomposition for compressed fault diagnosis of roller bearings Sensors 2016 16 9 1524 1532 10.3390/s160915242-s2.0-84988485591 Yasir M. Koh B.-H. Data decomposition techniques with multi-scale permutation entropy calculations for bearing fault diagnosis Sensors 2018 18 4 1278 10.3390/s180412782-s2.0-85045903597 Wu S.-D. Wu P.-H. Wu C.-W. Ding J.-J. Wang C.-C. Bearing fault diagnosis based on multiscale permutation entropy and support vector machine Entropy 2012 14 8 1343 1356 10.3390/e140813432-s2.0-84867606873 Wang F. Zhang Y. Zhang B. Su W. Application of wavelet packet sample entropy in the forecast of rolling element bearing fault trend Proceedings of the 2011 International Conference on Multimedia and Signal Processing May 2011 Guilin, China 10.1109/cmsp.2011.932-s2.0-80051867188 Zhang N. Wu L. Yang J. Guan Y. Naive Bayes bearing fault diagnosis based on enhanced independence of data Sensors 2018 18 2 463 10.3390/s180204632-s2.0-85041535454 Yau H. T. Wu S. Y. Chen C. L. Li Y. C. Fractional-order chaotic self-synchronization-based tracking faults diagnosis of ball bearing systems IEEE Transactions on Industrial Electronics 2016 63 6 3824 3833 10.1109/tie.2016.25229412-s2.0-84969940516 Pincus S. M. Approximate entropy as a measure of system complexity Proceedings of the National Academy of Sciences 1991 88 6 2297 2301 10.1073/pnas.88.6.22972-s2.0-0026015905 Richman J. S. Moorman J. R. Physiological time-series analysis using approximate entropy and sample entropy American Journal of Physiology-Heart & Circulatory Physiology 2000 278 6 H2039 H2049 10.1152/ajpheart.2000.278.6.h2039 Bandt C. Pompe B. Permutation entropy: a natural complexity measure for time series Physical Review Letters 2002 88 17 174102 10.1103/physrevlett.88.1741022-s2.0-84880568281 Zhang X. Liang Y. Zhou J. Zang Y. A novel bearing fault diagnosis model integrated permutation entropy, ensemble empirical mode decomposition and optimized SVM Measurement 2015 69 164 179 10.1016/j.measurement.2015.03.0172-s2.0-84926352537 Zheng J. Cheng J. Yang Y. Multiscale permutation entropy based rolling bearing fault diagnosis Shock and Vibration 2014 2014 8 154291 10.1155/2014/1542912-s2.0-84898034251 Li K. Su L. Wu J. Wang H. Chen P. A rolling bearing fault diagnosis method based on variational mode decomposition and an improved kernel extreme learning machine Applied Sciences 2017 7 10 1004 10.3390/app71010042-s2.0-85030651664 Huang N. E. Shen Z. Long S. R. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 1998 454 1971 903 995 10.1098/rspa.1998.01932-s2.0-5444236478 Li M. Wang H. Tang G. Yuan H. Yang Y. An improved method based on CEEMD for fault diagnosis of rolling bearing Advances in Mechanical Engineering 2014 6 676205 10.1155/2014/6762052-s2.0-84914694676 Rashid N. I. A. Shabri A. Samsudin R. New hybridization of empirical mode decomposition and least squares support vector machine model in forecasting Malaysia exchange rates Proceedings of the 2017 International Conference on Robotics, Automation and Sciences (ICORAS) November 2017 Melaka, Malaysia IEEE 1 5 10.1109/icoras.2017.83080692-s2.0-85050675152 Klionskiy D. M. Kaplun D. I. Geppener V. V. Empirical mode decomposition for signal preprocessing and classification of intrinsic mode functions Pattern Recognition and Image Analysis 2018 28 1 122 132 10.1134/s10546618180100912-s2.0-85044164162 Bandt C. Permutation entropy and order patterns in long time series Time Series Analysis and Forecasting 2016 Berlin, Germany Springer 61 73 10.1007/978-3-319-28725-6_5 Salazar-Alvarez M. I. Temblador-Pérez C. Conover W. J. Tercero-Gómez V. G. Cordero-Franco A. E. Beruvides M. G. Regressing sample quantiles to perform nonparametric capability analysis The International Journal of Advanced Manufacturing Technology 2016 86 5–8 1347 1356 10.1007/s00170-015-8285-62-s2.0-84953403222 Jentsch C. Leucht A. Bootstrapping sample quantiles of discrete data Annals of the Institute of Statistical Mathematics 2016 68 3 491 539 10.1007/s10463-015-0503-32-s2.0-84923261036 Moura M. D. C. Zio E. Lins I. D. Droguett E. Failure and reliability prediction by support vector machines regression of time series data Reliability Engineering & System Safety 2017 96 11 1527 1534 10.1016/j.ress.2011.06.0062-s2.0-80052449651 Lee C.-M. Ko C.-N. Time series prediction using RBF neural networks with a nonlinear time-varying evolution PSO algorithm Neurocomputing 2009 73 1–3 449 460 10.1016/j.neucom.2009.07.0052-s2.0-70449528756 García-Gonzalo E. Fernández-Martínez J. L. A brief historical review of particle swarm optimization (PSO) Journal of Bioinformatics and Intelligent Control 2012 1 1 3 16 10.1166/jbic.2012.1002 Wang H. Li R. Tang G. Yuan H. Zhao Q. Cao X. A compound fault diagnosis for rolling bearings method based on blind source separation and ensemble empirical mode decomposition PLoS One 2014 9 10 e109166 10.1371/journal.pone.01091662-s2.0-84907831000 Smith W. A. Randall R. B. Rolling element bearing diagnostics using the Case Western Reserve University data: a benchmark study Mechanical Systems and Signal Processing 2015 64-65 100 131 10.1016/j.ymssp.2015.04.0212-s2.0-84937975641