Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
An Investigation of Fractional Bagley–Torvik Equation
Previous Article in Journal
Continuous Variables Graph States Shaped as Complex Networks: Optimization and Manipulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Fuzzy Logic Classifier Based on Multiscale Permutation Entropy and Its Application in Bearing Fault Diagnosis

1
College of Mechanical Engineering, North University of China, Taiyuan 030051, China
2
School of Mechanical and Electrical Engineering, North University of China, Taiyuan 030051, China
3
School of Mechanical Engineering, Guizhou University, Guiyang 550025, China
4
State Key Laboratory of Marine Resource Utilization in South China Sea, Hainan University, Haikou 570228, China
5
Collage of Information and Communication Engineering, Hainan University, Haikou 570228, China
6
School of Energy and Power Engineering, North University of China, Taiyuan 030051, China
7
School of Computer Science and Information Engineering, Chongqing Technology and Business University, Chongqing 400067, China
*
Authors to whom correspondence should be addressed.
Entropy 2020, 22(1), 27; https://doi.org/10.3390/e22010027
Submission received: 11 November 2019 / Revised: 12 December 2019 / Accepted: 19 December 2019 / Published: 24 December 2019

Abstract

:
The self-organizing fuzzy (SOF) logic classifier is an efficient and non-parametric classifier. Its classification process is divided into an offline training stage, an online training stage, and a testing stage. Representative samples of different categories are obtained through the first two stages, and these representative samples are called prototypes. However, in the testing stage, the classification of testing samples is completely dependent on the prototype with the maximum similarity, without considering the influence of other prototypes on the classification decision of testing samples. Aiming at the testing stage, this paper proposed a new SOF classifier based on the harmonic mean difference (HMDSOF). In the testing stage of HMDSOF, firstly, each prototype was sorted in descending order according to the similarity between each prototype in the same category and the testing sample. Secondly, multiple local mean vectors of the prototypes after sorting were calculated. Finally, the testing sample was classified into the category with the smallest harmonic mean difference. Based on the above new method, in this paper, the multiscale permutation entropy (MPE) was used to extract fault features, linear discriminant analysis (LDA) was used to reduce the dimension of fault features, and the proposed HMDSOF was further used to classify the features. At the end of this paper, the proposed fault diagnosis method was applied to the diagnosis examples of two groups of different rolling bearings. The results verify the superiority and generalization of the proposed fault diagnosis method.

1. Introduction

Rotating machinery has been widely used in various modern industries such as wind turbines, aero engines, water turbines, and gas turbines. As a key component of rotating machinery, rolling bearings play an important role in rotating machinery [1,2,3,4,5]. Due to the complicated structure and harsh operating environment, various faults of rolling bearings (inner ring fault, outer ring fault, ball fault) is inevitable; thus, it is of great significance to study the fault detection methods and diagnostic techniques of rolling bearings [6,7,8,9,10]. In order to find out the fault location, the current frequently used fault diagnosis method is the time-frequency analysis method, such as the variable mode decomposition, local mean mode decomposition, empirical mode decomposition and so on [11,12,13,14].
The fault type is found by the fixed frequency of the intrinsic mode function after time-frequency analysis. However, due to the presence of noise, the collected fault signal may be submerged. In order to enhance the signal, filters are often used, such as minimum entropy deconvolution (MED), maximum correlation kurtosis deconvolution (MCKD), and multipoint optimal minimum entropy deconvolution adjusted (MOMEDA) [15,16,17]. However, the above algorithm is not adaptive, and misdiagnosis may occur. The current intelligent diagnosis method is generally to combine fault extraction technology with the machine learning method, and the first step is to extract fault feature information from the vibration signal [18]. However, since the equipment is usually inevitably operated under friction, vibration, and shock conditions, the vibration signal will show nonlinear and non-stationary characteristics. Since linear analysis methods cannot extract fault features, nonlinear analysis methods are particularly important for fault diagnosis of bearings. In recent years, many nonlinear dynamic methods such as sample entropy (SE), fuzzy entropy (FE), permutation entropy (PE), multiscale sample entropy (MSE), multiscale fuzzy entropy (MFE), multiscale permutation entropy (MPE), and improved methods based on them are used to extract nonlinear fault features [19,20]. For example, Yan et al. [21] extracted bearing fault features with the improved multiscale discrete entropy (MDE) and input them into the extreme learning machine (ELM), obtaining satisfactory fault diagnosis results. Liu et al. [22] used local mean decomposition to denoise the vibration data and then used MSE to extract the fault characteristics from the denoised signal. A lot of research studies have found that MPE has a faster calculation speed and stronger robustness than MSE, MFE, and MDE, and can better extract fault feature information [21,22]. Therefore, this paper uses MPE to extract fault features.
It is well known that after multiscale entropy is used to extract multiscale feature sets, feature reduction is needed to eliminate redundant features and improve computational speed. At present, the commonly used feature dimension reduction methods are principal component analysis (PCA) and linear discriminant analysis (LDA) [23]. For example, Aouabdi et al. [24] used multiscale sample entropy to extract the fault features of gears and then used PCA to reduce the dimensionality of fault features. Chen et al. [25] applied PCA to the feature reduction of high-speed train fault diagnosis. After using empirical mode decomposition to decompose the bearing data into the intrinsic mode function, Su et al. [26] extracted the high-dimensional feature vector set from the intrinsic mode function and then reduced it with LDA. PCA ignores other components while retaining the principal components with larger variance, so critical fault information may be lost during data dimensionality reduction. As a commonly used data dimension reduction method, LDA has a simple principle and a short calculation time. The feature set with dimensionality reduction has more sensitive features and is easier to classify [24,26].
The next and most important step after dimensionality reduction is to input the dimensionality reduction feature set into the classifier. Classification is one of the hot issues in machine learning research. In recent years, various methods of machine learning have been used in the field of fault diagnosis, such as support vector machine (SVM), decision tree (DT), k-nearest neighbor (KNN), extreme learning machine (ELM), etc. [27,28,29,30,31]. The self-organizing fuzzy logic classifier (SOF) has not been used in the field of fault diagnosis since it was proposed in 2018. SOF has the advantages of fast calculation speed, high classification accuracy, and no parameters [32]. The classification process is divided into three stages: the offline training stage, online training stage, and testing stage. In the offline training and online training stage, fuzzy rules of different categories are constructed after various types of qualified prototypes are obtained through self-iterative updating of meta-parameters. However, in the testing stage, the testing samples are classified according to the maximum similarity between the testing samples and the prototypes in each category. This does not take into account the impact of other prototypes in the same category on the classification of testing samples, so classification accuracy may be affected. This paper have improved the testing stage of SOF from the classification decision of SOF. The harmonic mean difference (HMDSOF) proposed in this paper not only considers the influence of other prototypes on testing samples but also assigns different weights to different prototypes. In the experimental part, the influence of the parameter g on the classification result of HMDSOF is analyzed by the bearing fault data of Case Western Reserve University, and the default value of the parameter g is given. Then, by comparing the classification results of HMDSOF with SOF, SVM, DT, KNN, ELM, least squares support vector machine (LSSVM), and kernel extreme learning machine (KELM), the validity and rationality of the proposed HMDSOF are illustrated. Finally, the generalization of HMDSOF is verified by bearing testing data of coal washer.

2. Basic Theory

2.1. Multiscale Permutation Entropy

MPE can be defined as the set of permutation entropy values of time series at different scales, and its calculation can be described as:
(1) Assuming a one-dimensional time series { x ( i ) , i = 1 ~ N } of length N. Set the embedded dimension as m and set the delay time as τ , and then conduct phase space reconstruction to obtain the matrix in the following form:
[ x ( 1 ) x ( 1 + τ ) x ( 1 + ( m 1 ) τ ) x ( j ) x ( j + τ ) x ( j + ( m 1 ) τ ) x ( K ) x ( K + τ ) x ( K + ( m 1 ) τ ) ]
where K is the number of reconstruction vectors, K = N ( m 1 ) τ .
To explain Formula (1), let us give an example, assuming that x = ( 4 , 8 , 9 , 6 , 5 , 11 , 7 ) . When τ = 1 , m = 3 , five embedding vectors can be obtained as:
[ 4 8 9 8 9 6 9 6 5 6 5 11 5 11 7 ] .
(2) Arrange the reconstruction matrix of each row according to the increasing rule:
x ( i + ( j 1 1 ) τ ) x ( i + ( j 2 1 ) τ ) x ( i + ( j m 1 ) τ ) .
It is important to note that if two equal elements exist in the reconstructed vector, the two elements are arranged in the original order. That is to say, suppose that p and q are any two numbers between 1 and m , if x ( i + ( j p 1 ) τ ) = x ( i + ( j q 1 ) τ ) and p < q , the following formula can be obtained.
x ( i + ( j p 1 ) τ ) x ( i + ( j q 1 ) τ )
(3) The symbol sequence corresponding to k reconstruction vectors of one-dimensional time series, whose permutation entropy is expressed as:
P E P ( m ) = j = 1 k P j ln P j
where P j represents the probability of any time series occurring.
(4) After steps (1), (2), and (3), the permutation entropy of the first scale is calculated. When calculating multiscale permutation entropy, it is necessary to use Formula (6) to conduct multiscale coarse granulation treatment on time series.
y j ( s ) = 1 s i = ( j 1 ) s + 1 j s x i ,   j = 1 , 2 , , N s
where s = 1 , 2 , , is the scale factor. y j ( s ) represents the coarse granulation time series of length N s , it can be seen from Formula (6) that the coarse graining process is achieved by calculating the average value of the time series. y j ( 1 ) ( j = 1 , 2 , , N ) is the original time series.
(5) After coarse granulation of time series, according to steps (1), (2) and (3), permutation entropy of different scales is calculated.

2.2. Linear Discriminant Analysis

Theoretically, the extracted multiscale permutation entropy set can be used to identify fault categories. The high-dimensional feature contains a lot of redundant information, so it is necessary to use the dimensionality reduction algorithm to reduce the dimension of the initial, which can not only avoid the dimension disaster, but also improve the performance of fault diagnosis. The role of LDA is to project a high-dimensional matrix into a low-dimensional matrix with minimal intraclass dispersion and maximum interclass dispersion. Assume that the calculated multiscale permutation entropy set is Y = [ y 1 , y 2 , , y n ] d × n , where n is the total number of samples, d is the dimension, and d = s = 32 . LDA will supervise the learning of a linear transformation matrix W d × m ( m d ) by itself. After the calculation as below, the high-dimensional data set y d is mapped to the low-dimensional data set x m .
x = W T y
Y is classified as Y = [ ϖ 1 , ϖ 2 , , ϖ C ] , C represents the number of categories. ϖ i d × n i is the data set of category i , and n i is the number of data samples in the category i .
The optimal projection matrix W should satisfy the following formula:
S t 1 S w W = W Λ
where S w = i = 1 C y π i ( y y i ¯ ) ( y y i ¯ ) T is defined as an intraclass discrete matrix, y i ¯ = 1 n i y j π i y j . S t = i = 1 n ( y i y ¯ ) ( y i y ¯ ) T is defined as a discrete matrix of the whole class, y ¯ = 1 n i = 1 n y i , Λ is the eigenmatrix of S t 1 S w .

2.3. Self-Organizing Fuzzy Logic Classifier

SOF is a fuzzy rule classifier without parameters. The algorithm includes three stages: the offline training stage, online training stage, and testing stage. In the first two stages, the fuzzy rules of each category were constructed based on the prototype of each category after the meta-parameters were updated iteratively, and the test samples were classified in the testing stage. The specific process is as follows:

2.3.1. Offline Training Stage

The role of the offline training stage is to find prototypes from different categories and build fuzzy rules that belong to different categories. Suppose there are a total of K samples, (the sample here refers to the low-dimensional feature vector processed by LDA), and the sample set belonging to the category c is { x } K c c = { x 1 c , x 2 c , , x K c c } ( { x } K c c { x } K ) , where c = ( 1 , 2 , , C ) . Since the same sample may appear more than once, the unique sample set and the occurrence frequency of each sample in it are expressed as { u } U K c c = { u 1 c , u 2 c , , u U K c c } and { f } U K c c = { f 1 c , f 2 c , , f U K c c } , respectively. K c is the number of samples of { x } K c c . It can be concluded that c = 1 C K c = K and c = 1 C U K c = U K , and U K is the unique sample set for all categories. It should be noted that the selection of the prototypes is carried out under the premise of the same category, and there is no influence between the samples of different categories when selecting the prototypes. The specific process of the offline training stage is as follows:
(1) The multimodal density corresponding to each unique sample is calculated according to Formula (9), where d represents Euclidean distance.
D K c M M ( u i c ) = f i c l = 1 K c j = 1 K c d 2 ( x i c , x j c ) 2 K c j = 1 K c d 2 ( u i c , x j c ) ,   i = 1 , 2 , , U K c
(2) Sort the samples according to the calculated multimode density and mutual distance. The sorted sample set is { r } = { r 1 , r 2 , , r U K c } , where r 1 = arg max i = 1 , 2 , , U K c ( D K c M M ( u i c ) ) , and r 2 is the sample that has the smallest distance from r 1 , that is, r 2 = arg min ( d ( r 1 , u i c ) ) i = 1 , 2 , , U K 1 c . r 3 is a sample with the smallest distance from r 2 , and so on. The multimodal density set of the sorted sample set { r } is represented as { D K c M M ( r ) } . Then, select the initial prototype according to Formula (10).
i f ( D K c M M ( r i ) > D K c M M ( r i + 1 ) ) a n d ( D K c M M ( r i ) > D K c M M ( r i 1 ) )     t h e n ( r i { p } 0 )
where { p } 0 represents a collection of initial prototypes.
(3) In order to increase the number of initial prototypes, the initial prototypes selected by Equation (11) is used as the center to attract nearby samples to form a data cloud.
w i n n i n g   p r o t o t y p e = arg min ( d ( x i , p ) ) , p { p } 0 , x i { x } K c c
It is important to note here that as mentioned above, the sample x i may not be unique, so the data cloud may not consist of only two samples.
(4) Define the set { p } 0 of the initial prototypes obtained by the Formula (10) as { φ } 0 ; that is, define the set of the data cloud center as { φ } 0 . Recalculate the multimodal density according to Equation (12).
D K c M M ( φ i ) = S i l = 1 n j = 1 n d 2 ( φ l , φ j ) 2 K j = 1 n d 2 ( φ i , φ j )
where φ i { φ } 0 , S i is the number of samples in the i t h data cloud, and n is the number of elements in the set { φ } 0 .
(5) According to Formula (13), the set { φ } i n e i g b o r of adjacent centers of each data cloud center is composed.
i f ( d 2 ( φ i , φ j ) G K c c , L )   t h e n ( φ j { φ } i n e i g b o r )
G K c c , L is the average radius of the locally affected area around the data sample corresponding to the level of granularity L , with a default value of L = 12 . The calculation process is as shown in Equation (14):
G K c c , L = d 2 ( x , y ) Q K c c , L
where x , y { x } K c c , x y , d 2 ( x , y ) G K c c , L 1 , and G K c c , L 1 is the average radius of the granularity level L 1 . Q K c c , L is the number of times that the distance between any two samples in { x } K c c is less than G K c c , L 1 . Q K c c , 1 is the number of times that the distance between any two samples in { x } K c c is less than the average square distance d K c c ¯ .
d K c c ¯ = 1 ( K c ) 2 l = 1 K c j = 1 K c d 2 ( x l c , x j c )
(6) According to Formula (16), select the most representative prototype { p } c in the category c from the center of the data cloud.
i f ( D K c M M ( φ i ) > D K c M M ( φ ) )   t h e n ( φ i { p } c )
where φ { φ } i n e i g b o r .
(7) After determining the representative prototypes of category c , according to Formula (17), AnYa type fuzzy rules belonging to each category are constructed, where N c is the number of prototypes in { p } c .
i f ( x ~ p 1 c ) o r ( x ~ p 2 c ) o r o r ( x ~ p N c c )   t h e n ( x ( c a t e g o r y   c ) )
where x represents a training sample, and ~ represents similarity.

2.3.2. Online Training Stage

After the offline training stage, it is followed by the input of online training samples to continue training. The purpose of the online training stage is to continue to select prototypes, update the meta-parameters of the classifier, and improve the classification accuracy of the test samples. The online training process is based on the assumption that the samples are stream data that appear one by one. When the online training sample is input, it is assumed that the new sample of the category c is x K c + 1 c , and the sample set after increasing the sample is defined as { x } K c + 1 c . In order to improve the computational efficiency, the average radius of the locally affected area will be calculated according to the new formula:
G K c + 1 c , L = d K c + 1 c ¯ d K c c ¯ G K c c , L .
whether the sample x K c + 1 c is a prototype will be determined according to Formula (19)
i f ( l = 1 K c + 1 j = 1 K c + 1 d 2 ( x l c , x j c ) 2 ( K c + 1 ) j = 1 K c + 1 d 2 ( x K c + 1 c , x j c ) > max ( l = 1 m m j = 1 m m d 2 ( p l c , p j c ) 2 ( m m ) j = 1 m m d 2 ( p , p j c ) ) ) o r ( l = 1 K c + 1 j = 1 K c + 1 d 2 ( x l c , x j c ) 2 ( K c + 1 ) j = 1 K c + 1 d 2 ( x K c + 1 c , x j c ) < min ( l = 1 m m j = 1 m m d 2 ( p l c , p j c ) 2 ( m m ) j = 1 m m d 2 ( p , p j c ) ) ) t h e n ( x K c + 1 c { p } c )
where m m is the number of elements in { p } c , ( p , p l c , p j c ) { p } c .
If Formula (19) is not satisfied, we can continue to judge whether the sample x K c + 1 c is a prototype according to Formula (20).
i f ( min ( d 2 ( x K c + 1 c , p ) ) > G K c + 1 c , L )   t h e n ( x K c + 1 c { p } c )
If any of Formula (19) or Formula (20) is satisfied, the meta-parameter of the SOF is updated as follows:
N c + 1 N c , x K c + 1 c p N c c , 1 S N c c , { p } c + p N c c { p } c .
If neither Formulas (19) nor (20) is satisfied, the sample is assigned to the nearest prototype, that is, p n * c = arg min ( d ( x K c + 1 c , p ) ) . The corresponding meta-parameters are updated as follows:
S n * c S n * c + 1 p n * c + 1 S n * c + 1 x K c + 1 c p n * c , S n * c + 1 S n * c .
After that, Equation (17) will be updated accordingly. The SOF classifier is ready to process the next data sample or enter the testing stage.

2.3.3. Testing Stage

The role of the testing stage is to classify the input testing samples. Assuming that the testing sample set is { z 1 , z 2 , z v v } , in order to determine the category of a testing sample z i i , the classification process of SOF is as follows:
(1)
According to Formula (23), calculate the similarity between each prototype selected in the first two stages and the testing sample.
S i m i l a r i t y = e d 2 ( z i i , p ) , p { p } c , c = 1 , 2 , , C , i i = 1 , 2 , , v v
(2)
Classify the testing sample into the category of the prototype that has the greatest similarity to the testing sample.
l a b e l ( z i i ) = arg max c   S i m i l a r i t y ( z i i , p )

3. Proposed HMDSOF

After the offline training phase and the online training phase, SOF selects a number of representative prototypes from each category of samples, and when selecting prototypes, different categories of samples will not affect each other. However, in the testing stage, the classification of the testing sample is only related to the prototype with the greatest similarity to the testing sample. The effect of other prototypes in the same category on the testing sample classification decisions is not considered, which affects the classification accuracy. In order to improve the classification accuracy, in this paper, we propose a SOF classifier based on harmonic mean difference, which is called HMDSOF. The first two stages of HMDSOF are the same as those of SOF. In the testing stage of HMDSOF, the category of the testing sample is determined and assigned a label corresponding to the category by calculating the harmonic mean difference between the testing sample and each prototype. The content of innovation is mainly two points: (1) The influence of different prototypes in the same category on the classification of the testing sample is considered by calculating multiple local mean vectors in the samples of each category. (2) In order to distinguish the influence of different prototypes in the same category on the testing sample classification decision, the harmonic mean difference constructed by introducing the concept of the harmonic mean is used as the decision of the testing sample classification. In addition, prototypes that differ slightly from test samples have greater weight in their classification decisions. The main process is as follows.
(1) Calculate the similarity between each prototype in each category and the test sample using Equation (22), and arrange the results in descending order.
(2) In each category, the corresponding prototype is sorted according to the result of similarity ranking—that is to say, the prototype with greater similarity to the test sample z i i is sorted in front.
(3) According to Formula (25), calculate the local average vector of the prototype in category c after sorting.
a i c = 1 i h i p h R c , c = 1 , 2 , , C
where 1 i g , and g is a parameter set before the test stage. In addition, g cannot be bigger than the minimum of the number of prototypes in every category—that is g min ( N c ) . It is easy to conclude that g is also the number of local average vectors a i R c , and a 1 c = p 1 R c .
(4) Construct the harmonic mean difference by introducing the concept of the harmonic mean value. Suppose there is a sample set { y 1 , y 2 , , y g } with g elements, and its harmonic mean value is calculated as shown in Formula (26). The calculation of difference is shown in Formula (27). It can be seen that the value range of the difference is D i f f e r e n c e ( z i i , p ) 1 . The harmonic mean difference is the sum of the harmonic mean value of the difference between each prototype in the same category and the testing sample. This paper applies the proposed harmonic mean difference to the SOF-based classification decision. The harmonic mean difference is defined as H M D ( . ) , and the calculation process is as shown in Equation (28).
H M ( { y 1 , y 2 , , y g } ) = g i = 1 g 1 y i
D i f f e r e n c e = e d 2 ( z i i , p ) , p { p } c , c = 1 , 2 , , C , i i = 1 , 2 , , v v
H M D ( z i i , { a i c } i = 1 g ) = g i = 1 g 1 D i f f e r e n c e ( z i i , a i c )
To illustrate the role of harmonic mean differences in assigning different weights to different prototypes in the same category, Equation (29) is given. As can be seen from Equation (29), when setting parameter g , H M D ( z i i , { a i c } i = 1 g ) is a fixed value. As a result, the prototype that has a smaller difference from the testing sample will be given more weight.
  H M D ( z i i , { a i c } i = 1 g )   D i f f e r e n c e ( z i i , a i c ) = [ g i = 1 g 1 D i f f e r e n c e ( z i i , a i c ) ]   D i f f e r e n c e ( z i i , a i c ) = g × ( 1 / ( D i f f e r e n c e ( z i i , a i c ) × i = 1 g 1 D i f f e r e n c e ( z i i , a i c ) ) ) 2 = ( H M D ( z i i , { a i c } i = 1 g ) ) 2 g × ( D i f f e r e n c e ( z i i , a i c ) ) 2 = 1 g × ( H M D ( z i i , { a i c } i = 1 g ) ) 2 ( D i f f e r e n c e ( z i i , a i c ) ) 2
(5) Assign the testing sample to the category with the smallest harmonic mean difference. It should be noted that when g = 1 , Formula (30) can be converted into Formula (24); that is, when g = 1 , the HMDSOF degenerates into SOF. Equation (31) expresses the relationship between SOF and HMDSOF.
l a b e l ( z i i ) = arg min c H M D ( z i i , { a i c } i = 1 g )
H M D ( z i i , { a i c } i = 1 g ) = D i f f e r e n c e ( z i i , a 1 c ) = 1 S i m i l a r i t y ( z i i , a 1 c )

4. Proposed Fault Diagnosis Method

The fault diagnosis method proposed in this paper is shown in Figure 1: after the vibration signal is collected, the multiscale permutation entropy set is firstly extracted. The parameters of multiscale permutation entropy selected in this paper are an embedded dimension of m = 6 and delay time of τ = 1 . In order to obtain the signal features as much as possible, the scale factor is set as s = 32 [33,34]. It can be seen that such a feature set has many scales and the entropy values are crossed together, which is not conducive to the final classification. Therefore, in this paper, linear discriminant analysis (LDA) is used to conduct dimensionality reduction for the multiscale permutation entropy feature set and the dimension of the feature set after the dimension reduction is nine. Then, the reduced dimensional feature set is randomly divided into online training samples, offline training samples, and testing samples. Finally, the proposed HMDSOF classifier is used for classification. After the training parameters of the HMDSOF are updated in two training stages, the testing samples are classified. For the convenience of description, this fault diagnosis method is named MPE-LDA-HMDSOF.

5. Experiments

5.1. Experiment 1

In Experiment 1, the experimental data of rolling bearings provided by Case Western Reserve University (CWRU) is used to verify the effectiveness of the proposed method. The experimental equipment is shown in Figure 2. It consists mainly of a three-phase induction motor, a torque sensor, and a load motor. The testing bearings are 6205-2RS (SKF, Sweden) deep groove bearings. The vibration acceleration signal of the bearing is obtained from the driving end under the condition of a rotation speed of 1797 r/min and a sampling frequency of 12 kHz. The bearing vibration signals are first classified into four categories, namely ordinary rolling bearings (normal) and rolling bearings with ball failure (B), outer ring failure (OR), and inner ring failure (IR). The faulty bearing is formed on the normal bearing by using electro-discharge machining (EDM), and each fault condition is classified according to the fault size of 0.007, 0.014, and 0.021 inches (1 inch = 25.4 mm), so the bearing vibration signal is finally classified into 10 categories. The first 102,400 points under each category are divided into 50 non-overlapping data samples on average; that is, 2048 sampling points are taken as a sample, and 50 samples can be obtained for each category, for a total of 500 samples. A detailed description of the class label is given in Table 1. The time-domain waveforms of their typical vibration signals are shown in Figure 3. So, a multiscale permutation entropy feature set with the size of 500 × 32 is obtained. The results of the multiscale permutation entropy corresponding to the vibration signal of Figure 3 are shown in Figure 4. In this paper, 10 samples are randomly selected from each category to form the online training sample set. In the remaining samples, 10 samples are randomly selected in each category to form the offline training sample set, and then the remaining samples constitute the testing sample set. It is known that both the online training set and offline training set have 100 samples, and the test sample set has 300 samples.
The first line in Figure 3 is the time-domain signal corresponding to the normal bearing. The three time-domain waveforms on the left side below correspond to inner ring faults, ball faults, and outer ring faults, and their fault size is 0.007 inches. The three time-domain waveforms on the left side below correspond to inner ring faults, ball faults, and outer ring faults. In addition, their fault size is 0.014 inches. The faulty bearing of the three time-domain waveforms on the right has a fault size of 0.021 inches. In order to facilitate the classification of inner ring type faults of different scales, it is expressed as IRD = 0.007, IRD = 0.014, and IRD = 0.021. Correspondingly, the ball element faults of different sizes are expressed as BD = 0.007, BD = 0.014, and BD = 0.021. Outer ring faults of different sizes are expressed as ORD = 0.007, ORD = 0.014, ORD = 0.021.
Since the experiment in this paper is conducted under the condition of randomly select samples, in order to reduce the impact of contingency, the average value of 10 experiments is taken, and the maximum and minimum values of classification accuracy are given. In addition, the standard deviation of classification accuracy is given to analyze the stability of the classification method. In this paper, three different feature extraction methods (MPE, MPE-PCA, and MPE-LDA) are used to extract fault features and then used for the comparison between SOF and the proposed HMDSOF, and the comparison is listed in Table 2. All the methods are implemented on MATLAB R2016a version and tested on Intel Core CPU i5-6200U @2.30 GHz/4.00 GB RAM and a Win10 computer with a 64-bit operating system.
The contribution rate of each principal component of MPE after PCA treatment is listed in Table 3, and the first eight principal components of the cumulative contribution rate of 90% are selected to form a feature set. Since the minimum number of prototypes of the third category obtained after the end of training in this experiment is 4, the case of g 5 does not exist. It can be seen that the bigger the value of g is, the longer the classification time will be. When the fault feature extraction method is MPE (numbered 1–4), the average classification time of HMDSOF ( g = 3 ) consumes 1.0576 s more than that of SOF. The average classification accuracy of HMDSOF ( g = 3 ) is 0.7334% higher than that of SOF, and the standard deviation of the classification accuracy of HMDSOF ( g = 3 ) is 0.0005 lower than that of SOF. When the fault feature extraction method is MPE-PCA (numbered 5–8), compared with SOF, the average classification time of HMDSOF ( g = 3 ) is 1.0379 s longer, its average classification accuracy is 0.6% higher, and its standard deviation of classification accuracy is 0.0054 lower. When the fault feature extraction method is MPE-LDA (numbered 9–12), compared with SOF, the average classification time of HMDSOF ( g = 3 ) is 0.9893 s longer, and the classification accuracy standard deviation is reduced by 0.0022. In addition, the average accuracy of classification was only improved by 0.3667%, but the maximum accuracy of HMDSOF reached 100%, which was satisfactory. When the classification method is HMDSOF and different feature extraction methods are selected (for example, numbered 2, 6, 10, or numbered 3, 7, 11), the comparison of the five indicators shows the advantages of the proposed MPE-LDA-HMDSOF. In conclusion, three different fault extraction methods have shown a better classification effect than SOF after being used as an input of HMDSOF, which proves the effectiveness of the proposed HMDSOF. Under the premise of using the same classification method HMDSOF, the rationality of the proposed fault diagnosis method MPE-LDA-HMDSOF is proved by adopting different classifier inputs. In addition, as the value of g increases, the longer the classification takes, and when g = 3 , the classification efficiency of the HMDSOF classifier is optimal, so the default value of g is set to 3.
In order to make the proposed HMDSOF more convincing, this paper also compares it with other common classification methods, which are SVM, DT, KNN, ELM, least squares support vector machine (LSSVM), and kernel extreme learning machine (KELM), respectively. The input of each classification method is the features set processed by LDA after calculating multiscale permutation entropy. The training samples of the six classification methods as comparisons are the sum of the online training samples and offline training samples of the HMDSOF, and the test samples used by them are the same as those of HMDSOF. The penalty factor of a standard SVM is 100, and the kernel function is 0.01. The minimum number of father nodes of DT is 5. The nearest neighbor number of KNN is K = 5, and the number of hidden layer nodes of ELM is 100 [21,35]. The Gaussian kernel function of the LSSVM is 0.5. The kernel function of the KELM is RBF, and its regularization parameter is 10,000 [36,37,38,39]. The classification results are shown in Table 4.
It can be seen from Table 4 that the SVM has the lowest classification accuracy, and it can be seen from the standard deviation that the classification effect of this method on different testing samples is very different, and the classification algorithm is very unstable. The standard deviation of the classification accuracy of DT is 1.7525, the algorithm is very unstable, and the minimum classification accuracy is 9% lower than that of HMDSOF. The maximum classification accuracy of KNN is 99%, but the standard deviation of classification accuracy is 1.4915 higher than that of HMDSOF. The input of different samples has a great influence on KNN classification accuracy. The average classification accuracy of ELM is 2.0333% lower than HMDSOF, and the standard deviation of classification accuracy is 0.7316 higher than that of HMDSOF. Compared with SVM, the calculation speed and classification accuracy of LSSVM have been significantly improved. KELM has the fastest calculation speed, but its maximum and minimum classification accuracy are 1% lower than HMDSOF. In addition, from the standard deviation of classification accuracy, the KELM classification stability is not as good as the proposed HMDSOF. In a word, the classification accuracy of HMDSOF is the highest; thus, the classification result is the best.
In order to express the classification effects of various classification methods more intuitively, Figure 5 shows the classification results of various classification methods in the fifth experiment. SVM has the lowest classification accuracy. Seventy of the 300 samples do not match the real category. Among the 70 misclassified samples, 67 samples of different categories are classified into category 6, with an overall classification accuracy of 76.6667%. In the classification results of DT, 27 samples are misclassified, and the overall classification accuracy is 91%. In the classification results of KNN, six samples are misclassified, of which four samples in category 6 are classified as category 3, and one sample in category 6 is classified as category 9. In the nine categories, one sample is misclassified as category 5, and the overall classification accuracy of KNN reached 98%. A total of 10 samples in the classification result of ELM are misclassified, and its overall classification accuracy is 91%. In the classification results of SOF, four samples were misclassified, among which three samples in category 6 are classified as category 3, and one sample in category 9 is classified as category 5. The total classification accuracy of SOF is 98.6667%. There are 20 misclassified samples in the classification results of LSSVM, and its classification accuracy is 93.3333%. There are five misclassified samples in the classification results of KELM, and its classification accuracy is 98.3333%. In the classification results of proposed HMDSOF, there are no misclassified samples, and the classification accuracy is 100%.
In addition, in order to evaluate the results of this experiment from different perspectives, F s c o r e s was introduced [40]. Its calculation process is shown in Formulas (32)–(34).
F s c o r e s ( j ) = 2 × p r e c i s i o n ( j ) × Re c a l l ( j ) p r e c i s i o n ( j ) + Re c a l l ( j )   j = 1 , 2 , , 10
p r e c i s i o n ( j ) = C M j , j i = 1 10 C M j , i
Re c a l l ( j ) = C M j , j i = 1 10 C M i , j
where p r e c i s i o n ( j ) , Re c a l l ( j ) , and F s c o r e s ( j ) represent the precision, recall, and F-scores measures of the j-th predicted class; respectively [41]. The F s c o r e s of each category corresponding to the experimental results in Figure 5 is shown in Figure 6.

5.2. Experiment 2

The fan-end bearing of CWRU has proven to be a more complex database [35]. In Experiment 2, we use its data to verify the effectiveness of the proposed fault diagnosis method. All the parameters used in Experiment 2 are exactly the same as those in Experiment 1. The classification results of each classification method are shown in Table 5.

5.3. Experiment 3

This section uses the experimental data of the rolling bearing in the coal washer to verify the generalization of the proposed fault diagnosis method. The experimental device is shown in Figure 7a. The motor speed is 1500 r/min, and the sampling frequency is 10 KHz. There are two acceleration sensors used to measure the bearing signal, and the position of the measuring point is shown in Figure 7b. The two bearing models are NJ210 (NSK, Japan) and NJ405 (NSK, Japan), respectively. NJ210 has two states, normal and crack, and NJ405 also has two states, normal and peeling. Their fault status is shown in Figure 8. In order to distinguish the two bearings, NJ210 is defined as A and NJ405 as B, so the collected signals can be divided into four categories. Their classification is shown in Table 6, and the typical time-domain diagram corresponding to the four states is given in Figure 9.
After calculating the multiscale permutation entropy of the obtained experimental data, LDA is used for dimensionality reduction processing, and the feature set after dimensionality reduction is input into different classification methods for comparison. In this experiment, there are 200 samples for each state, among which 50 samples from each category were randomly selected as the online training samples of HMDSOF and SOF. Then, 50 samples from the remaining 150 samples were randomly selected for offline training, and the remaining 100 samples were used as the testing samples. There are 200 offline training samples, 200 online training samples, and 400 testing samples in HMDSOF and SOF. The training samples in SVM, DT, KNN, and ELM are the sum of the training samples and the offline training samples input to HMDSOF, and their testing samples are the testing samples used by HMDSOF. That is to say, among the four classification algorithms SVM, DT, KNN, and ELM, there are 400 training samples and 400 test samples. The other parameters used in Experiment 3 are the same as those used in Experiment 1, and the comparison results are shown in Table 7.
It can be concluded from Table 7 that among the six classification methods, SVM has the lowest classification accuracy and the worst classification effect. The classification result of KNN is the most unstable, the standard deviation of classification accuracy is the largest, and the classification time is the longest. From the four indicators of classification accuracy, the classification effect of SOF is better than that of SVM, DT, KNN, and ELM. Although the average classification time of HMDSOF is 0.717304 s more than SOF, its maximum classification accuracy is 1.25% higher than SOF, its minimum classification accuracy is 1.5% higher than SOF, its average classification accuracy is 1.375% higher than SOF, and its classification standard deviation is 0.108864 lower than SOF; such results are satisfactory. The standard deviation of the classification accuracy of LSSVM is very close to that of HMDSOF, but the average classification accuracy is 4.1% lower than that of HMDSOF. KELM has the fastest classification speed and the shortest classification time; However, its maximum classification accuracy is 0.5% lower than HMDSOF, and the average classification accuracy is 1.625% lower than HMDSOF.

6. Conclusions

In this paper, a new SOF classifier (HMDSOF) based on the harmonic mean difference is proposed. Based on this, a new bearing fault diagnosis method is proposed. The validity and generalization of the proposed fault diagnosis method are verified by the bearing experimental data of Case Western Reserve University and the bearing experimental data of coal washer. The following conclusions can be drawn in this paper.
(1)
As the parameter g increases, the classification time of HMDSOF increases. When g = 3 , the classification effect of HMDSOF is optimal.
(2)
Under the premise of the same input, the proposed classification effect of HMDSOF is always higher than that of SOF, and the classification effect is better. By comparing with SVM, DT, KNN, ELM, LSSVM, and KELM, the proposed HMDSOF has higher classification accuracy and can be better used for bearing fault diagnosis.
(3)
By changing the input of the classifier, it is proved that the proposed bearing fault diagnosis method MPE-LDA-HMDSOF has better classification performance, and the classification accuracy reaches 100%.

Author Contributions

Data curation, J.W., C.L. and L.W.; Formal analysis, J.W., J.Z. and Y.S.; Methodology, W.D., X.G. and Z.W.; Writing-original draft, X.G., H.G. and X.Y.; Writing-review and editing, M.Y., G.W. and H.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Hainan Province under Grant 617079; the National Natural Science Foundation of China under Grants 51905496, 51605061, 61865005, and 61601413; the Provincial Natural Science Foundation of China under Grants 201801D221237, 201801D121185, 201801D121186, 201801D121187, 201701D121061, 201701D221146, and cstc2017jcyjAX0183; the Science and technology innovation project of Shanxi Province University Grant 2019L0581, and in part by the Science Foundation of North University of China under Grant XJJ201802.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, X.X.; Shen, C.Q.; Shi, J.J.; Zhu, Z.K. Initial center frequency-guided VMD for fault diagnosis of rotating machines. J. Sound Vib. 2018, 435, 36–55. [Google Scholar] [CrossRef]
  2. Li, K.; Xiong, M.; Li, F.C. A novel fault diagnosis algorithm for rotating machinery based on a sparsity and neighborhood preserving deep extreme learning machine. Neurocomputing 2019, 350, 261–270. [Google Scholar] [CrossRef]
  3. Wang, Z.J.; Zhou, J.; Wang, J.Y. A novel Fault Diagnosis Method of Gearbox Based on Maximum Kurtosis Spectral Entropy Deconvolution. IEEE Access 2019, 7, 29520–29532. [Google Scholar] [CrossRef]
  4. Wang, X.L.; Zhou, F.C.; He, Y.L.; Wu, Y.J. Weak fault diagnosis of rolling bearing under variable speed condition using IEWT-based enhanced envelope order spectrum. Meas. Sci. Technol. 2019, 30. [Google Scholar] [CrossRef]
  5. Wang, Z.J.; Du, W.H.; Wang, J.Y. Research and application of improved adaptive MOMEDA fault diagnosis method. Measurement 2019, 140, 63–75. [Google Scholar] [CrossRef]
  6. Wang, X.L.; Yan, X.L.; He, Y.L. Weak Fault Feature Extraction and Enhancement of Wind Turbine Bearing Based on OCYCBD and SVDD. Appl. Sci. 2019, 9, 3706. [Google Scholar] [CrossRef] [Green Version]
  7. Gao, Y.D.; Villecco, F.; Li, M. Multi-Scale Permutation Entropy Based on Improved LMD and HMM for Rolling Bearing Diagnosis. Entropy 2017, 19, 176. [Google Scholar] [CrossRef] [Green Version]
  8. Song, W.Q.; Cattani, C.; Chi, C.H. Fractional Brownian Motion and Quantum-Behaved Particle Swarm Optimization for Short Term Power Load Forecasting: An Integrated Approach. Energy 2019, in press. [Google Scholar]
  9. Duan, J.; Shi, T.L.; Duan, J. A narrowband envelope spectra fusion method for fault diagnosis of rolling element bearings. Meas. Sci. Technol. 2018, 29, 125106. [Google Scholar] [CrossRef]
  10. Song, W.Q.; Cheng, X.X.; Cattani, C. Multi-Fractional Brownian Motion and Quantum-Behaved Partial Swarm Optimization for Bearing Degradation Forecasting. Complexity 2019, in press. [Google Scholar]
  11. Yang, Y.Z.; Yang, W.G.; Jiang, D.X. Simulation and experimental analysis of rolling element bearing fault in rotor-bearing-casing system. Eng. Fail. Anal. 2018, 92, 205–221. [Google Scholar] [CrossRef]
  12. Zhang, Y.S.; Gao, Q.W.; Lu, Y.X. A Novel Intelligent Method for Bearing Fault Diagnosis Based on Hermitian Scale-Energy Spectrum. IEEE Sens. J. 2018, 18, 6743–6755. [Google Scholar] [CrossRef]
  13. Wang, Z.J.; Wang, J.Y.; Cai, W.N.; Zhou, J.; Du, W.H. Application of an Improved Ensemble Local Mean Decomposition Method for Gearbox Composite Fault diagnosis. Complexity 2019, 2019. [Google Scholar] [CrossRef] [Green Version]
  14. Cai, Z.Y.; Xu, Y.B.; Duan, Z.S. An alternative demodulation method using envelope-derivative operator for bearing fault diagnosis of the vibrating screen. J. Vib. Control. 2018, 24, 3249–3261. [Google Scholar] [CrossRef]
  15. Ciabattoni, L.; Ferracuti, F.; Freddi, A.; Monteriu, A. Statistical Spectral Analysis for Fault Diagnosis of Rotating Machines. IEEE Trans. Ind. Electron. 2018, 65, 4301–4310. [Google Scholar] [CrossRef]
  16. Wang, Z.J.; He, G.F.; Du, W.H. Application of Parameter Optimized Variational Mode Decomposition Method in Fault Diagnosis of Gearbox. IEEE Access 2019, 7, 44871–44882. [Google Scholar] [CrossRef]
  17. Lv, Y.; Yuan, R.; Shi, W. Fault Diagnosis of Rotating Machinery Based on the Multiscale Local Projection Method and Diagonal Slice Spectrum. Appl. Sci. 2018, 8, 619. [Google Scholar] [CrossRef] [Green Version]
  18. Mo, Z.L.; Wang, J.Y.; Zhang, H.; Miao, Q. Weighted Cyclic Harmonic-to-Noise Ratio for Rolling Element Bearing Fault Diagnosis. IEEE Trans. Instrum. Meas. 2019, 1–11. [Google Scholar] [CrossRef]
  19. An, Z.H.; Li, S.M.; Xin, Y. An intelligent fault diagnosis framework dealing with arbitrary length inputs under different working conditions. Meas. Sci. Technol. 2019, 30, 125107. [Google Scholar] [CrossRef]
  20. Yuan, R.; Lv, Y.; Li, H.W.X.; Song, G.B. Robust Fault Diagnosis of Rolling Bearing Using Multivariate Intrinsic Multiscale Entropy Analysis and Neural Network Under Varying Operating Conditions. IEEE Access 2019, 7, 130804–130819. [Google Scholar] [CrossRef]
  21. Yan, X.A.; Jia, M.P. Intelligent fault diagnosis of rotating machinery using improved multiscale dispersion entropy and mRMR feature selection. Knowl. Based Syst. 2019, 163, 450–471. [Google Scholar] [CrossRef]
  22. Liu, H.H.; Han, M.H. A fault diagnosis method based on local mean decomposition and multi-scale entropy for roller bearings. Mech. Mach. Theory 2014, 75, 67–78. [Google Scholar] [CrossRef]
  23. Wang, Z.J.; Zheng, L.K.; Du, W.H. A novel method for intelligent fault diagnosis of bearing based on capsule neural network. Complexity 2019, 2019. [Google Scholar] [CrossRef] [Green Version]
  24. Aouabdi, S.; Taibi, M.; Bouras, S.; Boutasseta, N. Using multi-scale entropy and principal component analysis to monitor gears degradation via the motor current signature analysis. Mech. Syst. Signal Process. 2017, 90, 298–316. [Google Scholar] [CrossRef]
  25. Chen, H.T.; Jiang, B.; Lu, N.Y. A Newly Robust Fault Detection and Diagnosis Method for High-Speed Trains. IEEE Trans. Intell. Transp. Syst. 2019, 20, 2198–2208. [Google Scholar] [CrossRef]
  26. Su, Z.Q.; Tang, B.P.; Liu, Z.R.; Qin, Y. Multi-fault diagnosis for rotating machinery based on orthogonal supervised linear local tangent space alignment and least square support vector machine. Neurocomputing 2015, 157, 208–222. [Google Scholar] [CrossRef]
  27. Azami, N.; Arnold, S.E.; Sanei, S. Multiscale Fluctuation-Based Dispersion Entropy and Its Applications to Neurological Diseases. IEEE Access 2019, 7, 68718–68733. [Google Scholar] [CrossRef]
  28. Xia, M.; Li, T.; Xu, L. Fault Diagnosis for Rotating Machinery Using Multiple Sensors and Convolutional Neural Networks. IEEE/ASME Trans. Mechatron. 2018, 23, 101–110. [Google Scholar] [CrossRef]
  29. Udmale, S.S.; Singh, S.K. Application of Spectral Kurtosis and Improved Extreme Learning Machine for Bearing Fault Classification. IEEE Trans. Instrum. Meas. 2019, 68, 4222–4233. [Google Scholar] [CrossRef]
  30. Pan, Z.B.; Wang, Y.D.; Ku, W.P. A new k-harmonic nearest neighbor classifier based on the multi-local means. Expert Syst. Appl. 2017, 67, 115–125. [Google Scholar] [CrossRef]
  31. Wang, Z.J.; Zheng, L.K.; Wang, J.Y.; Du, W.H. Research of novel bearing fault diagnosis method based on improved krill herd algorithm and kernel Extreme Learning Machine. Complexity 2019, 2019. [Google Scholar] [CrossRef]
  32. Gu, X.W.; Angelov, P.P. Self-organising fuzzy logic classifier. Inf. Sci. 2018, 447, 36–51. [Google Scholar] [CrossRef] [Green Version]
  33. Dávalos, A.; Jabloun, M.; Ravier, P.; Buttelli, O. On the Statistical Properties of Multiscale Permutation Entropy: Characterization of the Estimator’s Variance. Entropy 2019, 21, 450. [Google Scholar] [CrossRef] [Green Version]
  34. Huo, Z.Q.; Zhang, Y.; Shu, L. A New Bearing Fault Diagnosis Method based on Fine-to-Coarse Multiscale Permutation Entropy, Laplacian Score and SVM. IEEE Access 2019, 7, 17050–17066. [Google Scholar] [CrossRef]
  35. Zheng, L.K.; Wang, Z.J.; Zhao, Z.Y.; Wang, J.Y.; Du, W.H. Research of Bearing Fault Diagnosis Method Based on Multi-Layer Extreme Learning Machine Optimized by Novel Ant Lion Algorithm. IEEE Access 2019, 7, 89845–89856. [Google Scholar] [CrossRef]
  36. Li, Y.J.; Zhang, W.H.; Xiong, Q. A rolling bearing fault diagnosis strategy based on improved multiscale permutation entropy and least squares SVM. J. Mech. Sci. Technol. 2017, 31, 2711–2722. [Google Scholar] [CrossRef]
  37. Rodriguez, N.; Alvarez, P.; Barba, L. Combining Multi-Scale Wavelet Entropy and Kernelized Classification for Bearing Multi-Fault Diagnosis. Entropy 2019, 21, 152. [Google Scholar] [CrossRef] [Green Version]
  38. Mao, W.T.; Feng, W.S.; Liang, X.H. A novel deep output kernel learning method for bearing fault structural diagnosis. Mech. Syst. Signal Process. 2019, 117, 293–318. [Google Scholar] [CrossRef]
  39. Rostaghi, M.; Ashory, M.R.; Azami, H. Application of dispersion entropy to status characterization of rotary machines. J. Sound Vib. 2019, 438, 291–308. [Google Scholar] [CrossRef]
  40. Rodriguez, N.; Cabrera, G.; Lagos, C. Stationary Wavelet Singular Entropy and Kernel Extreme Learning for Bearing Multi-Fault Diagnosis. Entropy 2017, 19, 541. [Google Scholar] [CrossRef] [Green Version]
  41. Ferri, C.; Hernández-Orallo, J.; Modroiu, R. An experimental comparison of performance measures for classification. Pattern Recognit. Lett. 2009, 30, 27–38. [Google Scholar] [CrossRef]
Figure 1. Fault diagnosis flow chart.
Figure 1. Fault diagnosis flow chart.
Entropy 22 00027 g001
Figure 2. Experimental equipment.
Figure 2. Experimental equipment.
Entropy 22 00027 g002
Figure 3. Time-domain diagram of a typical vibration signal of each state of bearing.
Figure 3. Time-domain diagram of a typical vibration signal of each state of bearing.
Entropy 22 00027 g003
Figure 4. Multiscale permutation entropy.
Figure 4. Multiscale permutation entropy.
Entropy 22 00027 g004
Figure 5. Classification results of each classification method in the fifth experiment. (a) Classification results of SVM; (b) Classification results of DT; (c) Classification results of KNN; (d) Classification results of ELM; (e) Classification results of LSSVM; (f) Classification results of KELM; (g) Classification results of SOF; (h) Classification results of HMDSOF.
Figure 5. Classification results of each classification method in the fifth experiment. (a) Classification results of SVM; (b) Classification results of DT; (c) Classification results of KNN; (d) Classification results of ELM; (e) Classification results of LSSVM; (f) Classification results of KELM; (g) Classification results of SOF; (h) Classification results of HMDSOF.
Entropy 22 00027 g005aEntropy 22 00027 g005b
Figure 6. F-scores for each method.
Figure 6. F-scores for each method.
Entropy 22 00027 g006
Figure 7. Experimental device and position of measuring points. (a) Experimental device; (b) position of measuring points.
Figure 7. Experimental device and position of measuring points. (a) Experimental device; (b) position of measuring points.
Entropy 22 00027 g007
Figure 8. Fault status of two bearings. (a) Bearing NJ210 with a crack; (b) Bearing NJ405 with a piece of peeling.
Figure 8. Fault status of two bearings. (a) Bearing NJ210 with a crack; (b) Bearing NJ405 with a piece of peeling.
Entropy 22 00027 g008
Figure 9. Time-domain diagram corresponding to each state.
Figure 9. Time-domain diagram corresponding to each state.
Entropy 22 00027 g009
Table 1. Labels for various states. B: Ball failure, IR: inner ring failure, OR: outer ring failure.
Table 1. Labels for various states. B: Ball failure, IR: inner ring failure, OR: outer ring failure.
StageFault Size (Inches)LabelStageFault Size (Inches)Label
Normal01B0.0146
IR0.0072OR0.0147
B0.0073IR0.0218
OR0.0074B0.0219
IR0.0145OR0.02110
Table 2. Classification results. HMDSOF: harmonic mean difference, LDA: linear discriminant analysis, MPE: multiscale permutation entropy, PCA: principal component analysis, SOF: self-organizing fuzzy.
Table 2. Classification results. HMDSOF: harmonic mean difference, LDA: linear discriminant analysis, MPE: multiscale permutation entropy, PCA: principal component analysis, SOF: self-organizing fuzzy.
Serial NumberMethodsClassification Accuracy (%)Time/s
MaximumMinimumAverageStd
1MPE-SOF98.333394.333396.73330.01390.4243
2MPE-HMDSOF ( g = 2 ) 98.666794.333397.23330.01391.3394
3MPE-HMDSOF ( g = 3 ) 98.666794.666797.46670.01341.4819
4MPE-HMDSOF ( g = 4 ) 98.666794.666797.46670.01341.6589
5MPE-PCA-SOF98.66679597.40.01280.5149
6MPE-PCA-HMDSOF ( g = 2 ) 98.66679597.83330.01141.3203
7MPE-PCA-HMDSOF ( g = 3 ) 9997980.00741.5528
8MPE-PCA-HMDSOF ( g = 4 ) 9997980.00741.7343
9MPE-LDA-SOF99.333396.666798.63330.00740.4532
10MPE-LDA-HMDSOF ( g = 2 ) 99.333397.333398.83330.00591.2799
11MPE-LDA-HMDSOF ( g = 3 ) 10098.3333990.00521.4425
12MPE-LDA-HMDSOF ( g = 4 ) 10098.3333990.00521.5943
Table 3. Results of PCA dimension reduction.
Table 3. Results of PCA dimension reduction.
Principal Component Eigenvalue (×10−4)Rate of Contribution %Cumulative Contribution Rate %Principal ComponentEigenvalue (×10−4)Rate of Contribution %Cumulative Contribution Rate %
110564.702764.7027170.72700.448095.1333
22314.173078.8757180.70930.437195.5704
37.31324.506583.3822190.67390.415395.9857
44.31342.658086.0402200.65570.404196.3898
52.36841.459487.4996210.64440.397196.7869
61.62891.003888.5034220.60920.375497.1623
71.30300.802989.3063230.57960.357297.5195
81.21430.748390.0546240.51470.317297.8367
91.14390.704990.7595250.49980.30898.1447
101.09410.674291.4337260.48870.301198.4458
110.99190.611292.0449270.48160.296898.7426
120.91730.565392.6102280.47050.289999.0325
130.88050.542693.1528290.43420.267699.3001
140.85760.528593.6813300.40590.250199.5502
150.83980.517594.1988310.38270.235899.786
160.78950.486594.6853320.34760.2142100.0002
Table 4. Classification results of various methods. DT: decision tree, ELM: extreme learning machine, KNN: k-nearest neighbor, LSSVM: least squares support vector machine, KELM: kernel extreme learning machine, SVM: support vector machine.
Table 4. Classification results of various methods. DT: decision tree, ELM: extreme learning machine, KNN: k-nearest neighbor, LSSVM: least squares support vector machine, KELM: kernel extreme learning machine, SVM: support vector machine.
Classification MethodClassification Accuracy (%)Time/s
MaximumMinimumAveragestd
SVM79.333373.333376.56671.86820.8304
DT95.333389.333392.53331.75251.556
KNN999498.53331.49670.7318
ELM97.66679596.96670.73710.1348
LSSVM94.666789.666792.61.41260.1149
KELM9997.333398.30.45830.0344
HMDSOF ( g = 3 ) 10098.3333990.00521.4425
Table 5. Classification results of various methods.
Table 5. Classification results of various methods.
Input of ClassifierClassification MethodsClassification Accuracy (%)Time/s
MaximumMinimumAverageStd
MPE-LDASVM8478.333381.96671.4566110.761977
DT92.33338589.43332.0494571.845077
KNN949192.36670.8749590.765174
ELM9187.333389.21.1469750.137345
SOF9490.333392.251.0562150.504949
LSSVM87.33338084.22.4594410.121701
KELM93.666791.333392.73330.7118210.025278
HMDSOF ( g = 3 ) 97.333392.666794.41.271911.791814
Table 6. Category labels for various states.
Table 6. Category labels for various states.
LabelCorresponding State
1A is normal, B is normal
2A has a crack, B is normal
3A is normal, B has a piece of peeling failure
4A has a crack, B has a piece of peeling failure
Table 7. Classification results of various methods.
Table 7. Classification results of various methods.
Input of ClassifierClassification MethodsClassification Accuracy (%)Time/s
MaximumMinimumAverageStd
MPE-LDASVM95.7591.2593.851.1467340.601842
DT95.7593.7594.7750.7198091.643364
KNN97.7591.596.051.6271142.029434
ELM96.7594.2595.90.8455770.17406
SOF9896.2597.050.5567760.426295
LSSVM9593.594.3250.4479110.095361
KELM98.7595.7596.81.150.029672
HMDSOF ( g = 3 ) 99.2597.7598.4250.4479121.143599

Share and Cite

MDPI and ACS Style

Du, W.; Guo, X.; Wang, Z.; Wang, J.; Yu, M.; Li, C.; Wang, G.; Wang, L.; Guo, H.; Zhou, J.; et al. A New Fuzzy Logic Classifier Based on Multiscale Permutation Entropy and Its Application in Bearing Fault Diagnosis. Entropy 2020, 22, 27. https://doi.org/10.3390/e22010027

AMA Style

Du W, Guo X, Wang Z, Wang J, Yu M, Li C, Wang G, Wang L, Guo H, Zhou J, et al. A New Fuzzy Logic Classifier Based on Multiscale Permutation Entropy and Its Application in Bearing Fault Diagnosis. Entropy. 2020; 22(1):27. https://doi.org/10.3390/e22010027

Chicago/Turabian Style

Du, Wenhua, Xiaoming Guo, Zhijian Wang, Junyuan Wang, Mingrang Yu, Chuanjiang Li, Guanjun Wang, Longjuan Wang, Huaichao Guo, Jinjie Zhou, and et al. 2020. "A New Fuzzy Logic Classifier Based on Multiscale Permutation Entropy and Its Application in Bearing Fault Diagnosis" Entropy 22, no. 1: 27. https://doi.org/10.3390/e22010027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop