Improvement of Fast Kurtogram Combined with PCA for Multiple Weak Fault Features Extraction

Demodulation plays an important role in fault feature extraction for rotating machinery. The fast kurtogram method was proved to be effective for rotating machinery demodulation. However, the demodulation effectiveness of fast kurtogram was poor for multiple fault features extraction under low signal-to-noise ratio. In this paper, an improved method of fast kurtogram, called P-kurtogram, is presented. The proposed method extracted the multiple weak fault features from multiple envelope signals-based principal component analysis. Compared with extracting features from one envelope signal of fast kurtogram, P-kurtogram showed a better demodulation performance for multiple faults. Combined with principal component analysis method, the proposed method also showed a good performance under low signal-to-noise ratio(SNR). By simulation analysis, the P-kurtogram method showed good performance for multiple modulation features extraction and robust performance in demodulation under low SNR. Then, the proposed method was demonstrated by applications of bearing faults detection and propeller detection. The results verified that the P-kurtogram has a better demodulation performance than fast kurtogram for multiple weak fault features extraction, especially under low signal-to-noise ratio. The proposed method provides a reliable basis for multiple weak fault features extraction of rotating machinery.


Introduction
Condition-based maintenance [1] is an effective maintenance policy in industrial enterprises. Robust fault diagnosis technique is a growing necessity in condition-based maintenance [2]. Fault features extraction methods are the key techniques for fault diagnosis of rotating machinery [3]. According to the mechanism of rotating machinery, amplitude modulation is the main modulation type. Demodulation is the best way to extract fault features in monitoring signals of rotating machinery. Various methods have been proposed for demodulation of rotating machinery, such as envelope demodulation [4], resonance demodulation [5], cyclostationary analysis [6,7], and spectrum kurtosis [8]. Among these demodulation techniques, spectral kurtosis is an effective and usual method for rotating machinery. Spectra kurtosis method was originally proposed by Dwyer to indicate the existence of transients in signals [9]. The core idea of spectral kurtosis is to select the narrow modulated band by calculating the kurtosis [10]. Antoni presented a formal definition of the spectral kurtosis in term of the Wold-Cramer representation, called kurtogram [11]. For the sake of on-line industrial applications, Antoni further developed a fast algorithm to enhance the computational efficiency, called fast kurtogram, which optimized the band-filter [12]. Spectral kurtosis has been widely applied to extract fault features of rotating machinery. Combining spectral kurtosis and squared envelope spectral analysis, Leite et al. presented an approach for bearing condition-based maintenance [13]. Mohamed et al. applied spectral kurtosis to wind turbines condition monitoring [14]. Lee et al. used spectral kurtosis to detect the signal arisen by the tip vortex cavitation [15]. Han et al. used the kurtosis of envelope spectrum to determine the cavitation inception speed quantitatively [16].
The early fault characteristics of rotating machinery carried by vibration signals are quite weak. Due to the low signal-to-noise ratio (SNR) and complex interferences, it is difficult to identify the weak fault features. Some researchers have proposed noise reduction methods to improve the SNR, such as Blind deconvolution [17], variational mode decomposition [18], Wavelet transform [19,20], sparse stacked denoising autoencoder [21,22], and higher order coherent spectra [23]. Therefore, some researchers have taken measures to enhancing the performance of kurtogram for weak fault feature extraction. Lei et al. proposed an improved kurtogram method, which took wavelet packet transform as filter instead of original filter [24]. Wang et al. proposed a meshing frequency modulation index, which is combined with kurtogram to identify the bearing fault features in significant interrupts [25]. Wang et al. presented an improved kurtogram for bearing fault feature extraction, which applied manifold learning method to signal denoise [26]. Maximum correlated kurtosis deconvolution was developed for enhancing the periodic impulses of fault signals, which was taken as a signal preprocessing method [27]. Xiang et al. applied probabilistic principal component analysis as denoise method for rolling element bearing fault signals [28]. Tomasz et al. presented a protrugram which had a better detection ability of modulating signals than fast kurtogram under low SNR, but a priori knowledge was required [29]. Antoni put forward an infogram method to capture the transient components in both domain, including time and frequency domain [30]. Miao proposed a gini index to improve the kurtogram for bearing faults detection [31]. Ali et al. presented an autogram which calculated the kurtosis of unbiased autocorrelation of the squared envelope of filtered signal instead of the envelope signal of filtered signal [32]. However, only one envelope signal of filtered signal is selected in spectral kurtosis, which is not conductive to enhance the performance for multiple weak fault features extraction.
In the present study, a novel demodulation method for multiple weak fault features extraction is proposed. The proposed method is an improvement of kurtogram, which combined kurtogram and principal component analysis (PCA). Comparing with fast kurtogram, the proposed method extracted features form multiple envelope signals, which shows a better performance for multiple weak fault features extraction under low SNR. The multiple weak fault signals are extracted by PCA, which is a classical method of multivariate statistical analysis for dimension reduction. PCA have been successfully applied in many applications [33,34]. This paper is organized as follows. In Section 2, an overview of fast kurtogram and the P-kurtogram method is presented. The proposed method is verified by simulations analysis in Section 3. Section 4 shows the application of the proposed method in bearing fault feature extraction and propeller features detection. In addition, the results were compared with the fast kurtogram to examine the performance of P-kurtogram method. Finally, the conclusions are drawn in Section 5.

The Overview of Fast Kurtogram
A brief overview of spectral kurtosis and fast kurtogram is presented in this subsection. Kurtosis is a statistical indicator which is used to detect faults impulsive component in signals [11]. For time series signals, kurtosis is defined as Equation (1).
where x(t) denotes time series signal, µ x is the mean value of x(t), t i and N denotes instant time and the length of x(t). Comparing with kurtosis for time series signal, spectrum kurtosis applied it to time-frequency domain. Therefore, K x ( f ) was used as an indicator to find the modulation component of filtered signal based on time-frequency analysis, which is defined as Equation (2) [12].
where < · > represents the time-frequency averaging operator, X(t, f ) is the time-frequency envelope of the band-pass filtered signal of x(t) around f.
Based on kurtosis values of filtered signals, the subband that includes impulsive component can be identified by kurtogram. The filter-bank structure is the core of kurtogram algorithm. A 1/3-binary tree filter bank strategy was proposed in fast kurtogram, which improved the computing efficiency. Using this technique, the whole band can be calculated efficiently. However, fast kurtogram has worse performance for multiple features extraction and low SNR condition, which restricts its application range.

P-Kurtogram
To improve the performance of fast kurtogram for multiple features extraction and low SNR condition extraction, an improvement method of fast kurtogram was proposed. The proposed method P-kurtogram combined fast kurtogram and PCA for multiple weak fault features extraction, which also could enhance the extraction performance under low SNR condition. The detail technique routines of fast kurtogram and P-kurtogram are shown in Figure 1. As shown in Figure 1, the P-kurtogram chooses multiple envelope signals of filtered signal based on kurtosis. Multiband selection based on fast kurtogram, PCA and feature extraction are three key steps for P-kurtogram. The number of data samples needs to be higher than the highest frequency of the resonant band. where ( ) denotes time series signal, is the mean value of ( ), and N denotes instant time and the length of ( ). Comparing with kurtosis for time series signal, spectrum kurtosis applied it to time-frequency domain. Therefore, ( )was used as an indicator to find the modulation component of filtered signal based on time-frequency analysis, which is defined as Equation (2) [12].
where <•> represents the time-frequency averaging operator, ( , )is the time-frequency envelope of the band-pass filtered signal of ( ) around f. Based on kurtosis values of filtered signals, the subband that includes impulsive component can be identified by kurtogram. The filter-bank structure is the core of kurtogram algorithm. A 1/3-binary tree filter bank strategy was proposed in fast kurtogram, which improved the computing efficiency. Using this technique, the whole band can be calculated efficiently. However, fast kurtogram has worse performance for multiple features extraction and low SNR condition, which restricts its application range.

P-Kurtogram
To improve the performance of fast kurtogram for multiple features extraction and low SNR condition extraction, an improvement method of fast kurtogram was proposed. The proposed method P-kurtogram combined fast kurtogram and PCA for multiple weak fault features extraction, which also could enhance the extraction performance under low SNR condition. The detail technique routines of fast kurtogram and P-kurtogram are shown in Figure 1. As shown in Figure 1, the Pkurtogram chooses multiple envelope signals of filtered signal based on kurtosis. Multiband selection based on fast kurtogram, PCA and feature extraction are three key steps for P-kurtogram. The number of data samples needs to be higher than the highest frequency of the resonant band. In this step, the reason for multiple filtered signals selection is put forward. Combined with fast kurtogram, the modulation components extraction process from multiple filtered signals is plotted in Figure 2.

Multiband Selection Based on Kurtogram
In this step, the reason for multiple filtered signals selection is put forward. Combined with fast kurtogram, the modulation components extraction process from multiple filtered signals is plotted in Figure 2. According to the mechanism of rotating machinery, signal modulation model can be expressed as Equation (3).
where, ( ) denotes raw signal, P denotes the total number of frequency bands, ( ) denotes modulation signal, denotes carrier wave, ( ) denotes the non-modulation signal and Q denotes its total number, ( ) denotes white noise.
According to the mechanism of rotating machinery, narrow band and broad band carrier signals exist all at once. The monitoring signal can be expressed as Equation (4), which are filtered into several bands equally. The modulation signal ( ) in frequency bands is correspond to modulation signal ( ) as expressed in Equation (5).
where, M denotes the total number of frequency bands, m denotes the index of frequency band, ( ) denotes modulation signal in frequency band m, which include the modulation waves, ( ) denotes carrier wave in frequency band m, ( ) denotes the non-modulation signal in frequency band m, ( )denotes white noise in frequency band m.
The filtered signal ( ) is expressed as Equation (6), which include modulation signal, nonmodulation signal and noise. For the filtered signal ( ), the envelop signal ( ) is expressed as Equation (7), which can be approximate as | ( )|.
From the signal modulation model of rotating machinery, weak fault modulated signal ( ) can be demodulated from multiple filtered signals including carrier waves ( ) . In order to increase the accuracy for multiple weak fault features extraction, first n filtered signals are both selected based on kurtosis value of kurtogram. In order to guarantee the performance of PCA, the total filtered signals' number n should be satisfied with: > 3. Based on multiple envelope signals selection, multiple components are extracted and weak fault signals are enhanced. Envelope signals matrix ( , ) is make up of multi envelope signals ( ) of filtered signal as shown in Equation (8). According to the mechanism of rotating machinery, signal modulation model can be expressed as Equation (3).
where, s(t) denotes raw signal, P denotes the total number of frequency bands, xm i (t) denotes modulation signal, xc i denotes carrier wave, xn j (t) denotes the non-modulation signal and Q denotes its total number, n(t) denotes white noise.
According to the mechanism of rotating machinery, narrow band and broad band carrier signals exist all at once. The monitoring signal can be expressed as Equation (4), which are filtered into several bands equally. The modulation signal sm m (t) in frequency bands is correspond to modulation signal xm i (t) as expressed in Equation (5).
where, M denotes the total number of frequency bands, m denotes the index of frequency band, sm m (t) denotes modulation signal in frequency band m, which include the modulation waves, sc m (t) denotes carrier wave in frequency band m, sn m (t) denotes the non-modulation signal in frequency band m, n m (t) denotes white noise in frequency band m.
The filtered signal s f m (t) is expressed as Equation (6), which include modulation signal, non-modulation signal and noise. For the filtered signal s f m (t), the envelop signal ES m (t) is expressed as Equation (7), which can be approximate as sm m (t) .
From the signal modulation model of rotating machinery, weak fault modulated signal sm m (t) can be demodulated from multiple filtered signals including carrier waves sc m (t). In order to increase the accuracy for multiple weak fault features extraction, first n filtered signals are both selected based on kurtosis value of kurtogram. In order to guarantee the performance of PCA, the total filtered signals' number n should be satisfied with: n > 3. Based on multiple envelope signals selection, multiple components are extracted and weak fault signals are enhanced. Envelope signals matrix ESM(i, j) is make up of multi envelope signals ES m (t) of filtered signal as shown in Equation (8).
where, i denotes ith envelope signal ES i (t), j denotes the instant time of ES n ( j).

Principal Component Analysis
In this step, the modulation components is extracted from envelope signals matrix ESM(i, j) by PCA as shown in Figure 3. PCA is statistical procedure, which is used as a dimensional reduction tool. For envelope signals matrix ESM(i, j), it is two-dimensional including sequence number i and instant time j. Based on Equation (8), the weak fault modulation signals can be seen as the principal components of ES i (t). We adopted PCA to reduce sequence number dimension of ESM(i, j) to multiple extract weak fault signals WFS(t) of envelope signal as shown in Figure 3.
where, i denotes ith envelope signal ( ), j denotes the instant time of ( ).

Principal Component Analysis
In this step, the modulation components is extracted from envelope signals matrix ( , ) by PCA as shown in Figure 3. PCA is statistical procedure, which is used as a dimensional reduction tool. For envelope signals matrix ( , ), it is two-dimensional including sequence number i and instant time j. Based on Equation (8), the weak fault modulation signals can be seen as the principal components of ( ). We adopted PCA to reduce sequence number dimension of ( , ) to multiple extract weak fault signals ( ) of envelope signal as shown in Figure 3. represents the covariance matrix with the corresponding column variances along the diagonal, which is expressed as Equation (9).
where (•) represents the covariance operator. Covariance matrix can be converted into a canonical form by eigendecomposition, whereby the matrix is represented in terms of its eigenvalues and eigenvectors as Equation (10).
[ , ] = eig( ) where eig(•) represents the eigendecomposition operator. V denotes diagonal matrix of eigenvalues as Equation (11). U denotes eigenvectors matrix whose columns are the corresponding eigenvectors as Equation (12). The criterion for choosing principal components is expressed as Equation (13).
= [ , , … , ] ∈ R × Where where k represents the total number of selected principal components. And the weak fault signals (t) can be reconstruct as Equation (14). ESMC cov represents the covariance matrix with the corresponding column variances along the diagonal, which is expressed as Equation (9).
where cov(·) represents the covariance operator. Covariance matrix ESMC cov can be converted into a canonical form by eigendecomposition, whereby the matrix is represented in terms of its eigenvalues and eigenvectors as Equation (10).
where eig(·) represents the eigendecomposition operator. V denotes diagonal matrix of eigenvalues λ i as Equation (11). U denotes eigenvectors matrix whose columns are the corresponding eigenvectors µ i as Equation (12). The criterion for choosing principal components is expressed as Equation (13). where where k represents the total number of selected principal components. And the weak fault signals WFS i (t) can be reconstruct as Equation (14).

Feature Extraction
From reconstructed principal component signals WFS i (t), characteristic frequency can be extracted by Fourier transform as in Equation (15). In order to achieve the accurate extraction of fault features, the spectral resolution needs to be greater than 0.2 Hz.
In conclusion, the proposed P-Kurtogram method is an improvement of Kurtogram. Based on PCA of envelope signals of filtered signals, the multiple modulation components are extracted and the weak fault features are enhanced under low SNR, which is verified by simulation analysis and applications in follow Sections.

Simulations
In this section, the effectiveness of proposed P-kurtogram method for multiple modulation features extraction was verified by simulations analysis. Due to line-spectrum, narrow band and broad band modulation signals exist simultaneously for rotating machinery, an amplitude modulation simulated signal is built as Equation (16) [35].
cos(2π f i t) (16) where i denotes the index of signal components, f m,1 , f m,2 = 15 Hz, 25 Hz denote the modulation wave frequencies, f c,1 , f c,2 = 900 Hz, 2300 Hz denote the carrier wave frequencies, f 1 , f 2 = 30 Hz, 500 Hz denote the frequencies of cosine signal, nc 1 (t) : 1100 Hz-1600 Hz and nc 2 (t) : 2400 Hz-2900 Hz denote the narrow band carrier waves. wc i (t) denotes the wide band carrier wave in form of white Gaussian noise. In order to study the performance of proposed P-kurtogram method for multiple modulation features extraction under low SNR condition, different white noise levels were added to signal x(t). The SNR levels of x(t) are −5 dB, −10 dB, −15 dB, respectively. The simulated signals were analyzed using fast kurtogram [11,12] and proposed method P-kurtogram simultaneously. The results are shown in Figures 4-6. As SNR = -5 dB, the kurtogram and proposed P-kurtogram method were preferable for feature detection. The raw signal and its amplitude spectrum is shown in the Figure 4a,b. Modulation wave frequencies 2 , , 2 , can be detected exactly in both method. The resonance band was detected by fast kurtogram as shown in Figure 4c, which is marked by solid line box. Amplitude spectrum of kurtogram is plotted in Figure 4d. For P-kurtogram, four envelope signals of filtered signals were selected for PCA as shown in Figure 4c. Two principal components were selected according to Equation (13)   Given SNR = −10 dB, the raw signal and its amplitude spectrum is shown in the Figure 5a,b. The resonance band is marked by solid line box as shown in Figure 5c. The fast kurtogram method was unable to extract the modulation wave frequency 2 , which was buried in the spectrum of kurtogram as shown in Figure 5d. By contrast, the proposed P-kurtogram had a better extraction performance. Three envelope signals of filtered signals were selected as shown in Figure 5c, which are marked by line box. Two principal components were selected as shown in Figure 5e. The modulation wave frequencies 2 , , 2 , can be clearly observed in Figure 5f.  Given SNR = −15 dB, the SNR is relatively low. The raw signal and its amplitude spectrum is shown in the Figure 6a,b. The fast kurtogram method was unable to extract the modulation wave frequency 2 , , 2 , which were buried in the spectrum of fast kurtogram as shown in Figure 6d. By contrast, the modulation wave frequencies 2 , , 2 , can be clearly observed in Figure 6f. Three envelope signals of filtered signals were selected as shown in Figure 6c. Two principal components were selected as shown in Figure 6e.
In conclusion, the proposed method P-kurtogram shows two advantages than fast kurtogram. The first advantages is that the proposed P-kurtogram has a good performance for multiple modulation components extraction. Secondly, the proposed P-kurtogram method has a better performance than fast kurtogram for weak feature extraction under low SNR.

Application of P-Kurtogram for Bearing Fault Signals
In this subsection, the proposed P-kurtogram was validated via bearing data. The bearing fault data in this study is provided by the Case Western Reserve University (CWRU) laboratory [36]. The vibration data were collected using accelerometers and sampling frequency is 12 KHz. Fan end bearing outer race fault and inner race fault data were adopted in this study. The main characteristic frequency of bearings consist of Ball Pass Frequency Outer , Ball Pass Frequency inner and As SNR = −5 dB, the kurtogram and proposed P-kurtogram method were preferable for feature detection. The raw signal and its amplitude spectrum is shown in the Figure 4a,b. Modulation wave frequencies 2 f m,1 , 2 f m,2 can be detected exactly in both method. The resonance band was detected by fast kurtogram as shown in Figure 4c, which is marked by solid line box. Amplitude spectrum of kurtogram is plotted in Figure 4d. For P-kurtogram, four envelope signals of filtered signals were selected for PCA as shown in Figure 4c. Two principal components were selected according to Equation (13) as shown in Figure 4e. P-kurtogram extracted the modulation frequency exactly as shown in Figure 4f.
Given SNR = −10 dB, the raw signal and its amplitude spectrum is shown in the Figure 5a,b. The resonance band is marked by solid line box as shown in Figure 5c. The fast kurtogram method was unable to extract the modulation wave frequency 2 f m,1 which was buried in the spectrum of kurtogram as shown in Figure 5d. By contrast, the proposed P-kurtogram had a better extraction performance. Three envelope signals of filtered signals were selected as shown in Figure 5c, which are marked by line box. Two principal components were selected as shown in Figure 5e. The modulation wave frequencies 2 f m,1 , 2 f m,2 can be clearly observed in Figure 5f.
Given SNR = −15 dB, the SNR is relatively low. The raw signal and its amplitude spectrum is shown in the Figure 6a,b. The fast kurtogram method was unable to extract the modulation wave frequency 2 f m,1 , 2 f m,2 which were buried in the spectrum of fast kurtogram as shown in Figure 6d. By contrast, the modulation wave frequencies 2 f m,1 , 2 f m,2 can be clearly observed in Figure 6f. Three envelope signals of filtered signals were selected as shown in Figure 6c. Two principal components were selected as shown in Figure 6e.
In conclusion, the proposed method P-kurtogram shows two advantages than fast kurtogram. The first advantages is that the proposed P-kurtogram has a good performance for multiple modulation components extraction. Secondly, the proposed P-kurtogram method has a better performance than fast kurtogram for weak feature extraction under low SNR.

Application of P-Kurtogram for Bearing Fault Signals
In this subsection, the proposed P-kurtogram was validated via bearing data. The bearing fault data in this study is provided by the Case Western Reserve University (CWRU) laboratory [36]. The vibration data were collected using accelerometers and sampling frequency is 12 KHz. Fan end bearing outer race fault and inner race fault data were adopted in this study. The main characteristic frequency of bearings consist of Ball Pass Frequency Outer f BPFO , Ball Pass Frequency inner f BPFI and Ball Spin Frequency f BSF , which can be calculated by Equation (17) [35]. Under bearing fault condition, these characteristic frequency can be detected by signal processing.
where, n denotes number of rollers, f denotes rotor speed, θ denotes contact angle, d p denotes pitch diameter and d b denotes the ball number.

Case 1: Fan End Bearing Inner Race Fault
Under fan end bearing inner race fault condition, the speed is 1772 rpm, fault diameter is 0.014 inches and motor load is 1 HP. The vibration signals of drive end were analyzed using fast kurtogram and proposed method P-kurtogram. The results are shown in Figure 7. Ball Spin Frequency , which can be calculated by Equation (17) [35]. Under bearing fault condition, these characteristic frequency can be detected by signal processing.
where, n denotes number of rollers, f denotes rotor speed, denotes contact angle, denotes pitch diameter and denotes the ball number.

Case 1: Fan End bearing inner Race Fault
Under fan end bearing inner race fault condition, the speed is 1772 rpm, fault diameter is 0.014 inches and motor load is 1 HP. The vibration signals of drive end were analyzed using fast kurtogram and proposed method P-kurtogram. The results are shown in Figure 7. The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 7a. The bearing inner race fault frequency was extracted, but other inferences frequency occur in the spectrum of the kurtogram as shown in Figure 7b. In this study, three filtered signals were selected for PCA and two principal components were selected as shown in Figure 7c. The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 7a. The bearing inner race fault frequency f BPFI was extracted, but other inferences frequency occur in the spectrum of the kurtogram as shown in Figure 7b. In this study, three filtered signals were selected for PCA and two principal components were selected as shown in Figure 7c. The bearing inner race fault frequency f BPFI was exactly extracted by P-kurtogram as shown in Figure 7d. The results indicate the proposed P-kurtogram method has a better performance than fast kurtogram for multiple features extraction.

Case 2: Fan End Bearing Outer Race Fault
Under fan end bearing outer race fault condition, the speed is 1750 rpm, fault diameter is 0.021 inches and motor load is 2 HP. The vibration signals of fan end were analyzed using fast kurtogram and proposed method P-kurtogram. The results are shown in Figure 8.  The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 8a, but the bearing outer race fault frequency is buried in the amplitude spectrum of the squared envelop as shown in Figure 8b. Three filtered signals were selected for PCA and two principal components were selected as shown in Figure 8c. The bearing outer race fault frequency was exactly extracted by P-kurtogram in Figure 8d. The results indicate the proposed P-kurtogram method has a better performance than fast kurtogram for weak feature extraction.

Application of P-Kurtogram for Acoustic Signals of Propeller
In this subsection, the proposed P-kurtogram was validated via acoustic signal of merchant ship propeller. The acoustic signals was collected by hydrophones. The hydrophone sensitivity is 39.7 mv/pa and its sampling frequency is set as 44,100 Hz. The shaft speed of propeller is 72 r/min and the propeller has six blades. Therefore, the characteristic shaft frequency is 1.2 Hz and blades frequency is 6 Hz. The acoustic signals of propeller was analyzed by fast kurtogram and proposed method P-kurtogram. The results are shown in Figure 9. The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 8a, but the bearing outer race fault frequency f BPFO is buried in the amplitude spectrum of the squared envelop as shown in Figure 8b. Three filtered signals were selected for PCA and two principal components were selected as shown in Figure 8c. The bearing outer race fault frequency f BPFO was exactly extracted by P-kurtogram in Figure 8d. The results indicate the proposed P-kurtogram method has a better performance than fast kurtogram for weak feature extraction.

Application of P-Kurtogram for Acoustic Signals of Propeller
In this subsection, the proposed P-kurtogram was validated via acoustic signal of merchant ship propeller. The acoustic signals was collected by hydrophones. The hydrophone sensitivity is 39.7 mv/pa and its sampling frequency is set as 44,100 Hz. The shaft speed of propeller is 72 r/min and the propeller has six blades. Therefore, the characteristic shaft frequency f shaft is 1.2 Hz and blades frequency f blade is 6 Hz. The acoustic signals of propeller was analyzed by fast kurtogram and proposed method P-kurtogram. The results are shown in Figure 9. The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 9a. The blade frequency is detected exactly, but the shaft frequency is buried in the amplitude spectrum of the squared envelop as shown in Figure 9b. Three filtered signals is selected for PCA and one principal component is selected as shown in Figure 9c. The blade frequency and shaft frequency are exactly extracted by P-kurtogram in Figure 9d. The results indicate the proposed P-kurtogram method has a better extraction performance than fast kurtogram for multiple components extraction, especially under weak feature extraction.

Conclusions
This paper proposed an improvement method of kurtogram for multiple modulation components extraction, called P-kurtogram. P-kurtogram extracted the impulsive and non-stationary signals from multi envelope signals of filtered signals based on principal component analysis, rather than one envelope signal of fast kurtogram. Combined with principal component analysis, the proposed method has a well ability for multiple fault features and weak fault signals extraction. By simulation analysis, the P-kurtogram method shows good performance for multiple modulation features extraction and robust performance in demodulation under low SNR. Then, the proposed method was demonstrated by applications of bearing faults detection and propeller detection. Results indicated that the proposed method provides a better multiple features extraction and weak characteristic fault frequencies identification than fast kurtogram. Therefore, the proposed method provides a reliable basis for multiple weak fault features extraction of rotating machinery.   The resonance frequency band using fast kurtogram method is marked by solid line box as shown in Figure 9a. The blade frequency f blade is detected exactly, but the shaft frequency f shaft is buried in the amplitude spectrum of the squared envelop as shown in Figure 9b. Three filtered signals is selected for PCA and one principal component is selected as shown in Figure 9c. The blade frequency f blade and shaft frequency f shaft are exactly extracted by P-kurtogram in Figure 9d. The results indicate the proposed P-kurtogram method has a better extraction performance than fast kurtogram for multiple components extraction, especially under weak feature extraction.

Conclusions
This paper proposed an improvement method of kurtogram for multiple modulation components extraction, called P-kurtogram. P-kurtogram extracted the impulsive and non-stationary signals from multi envelope signals of filtered signals based on principal component analysis, rather than one envelope signal of fast kurtogram. Combined with principal component analysis, the proposed method has a well ability for multiple fault features and weak fault signals extraction. By simulation analysis, the P-kurtogram method shows good performance for multiple modulation features extraction and robust performance in demodulation under low SNR. Then, the proposed method was demonstrated by applications of bearing faults detection and propeller detection. Results indicated that the proposed method provides a better multiple features extraction and weak characteristic fault frequencies identification than fast kurtogram. Therefore, the proposed method provides a reliable basis for multiple weak fault features extraction of rotating machinery.

Conflicts of Interest:
The authors declare no conflict of interest.