^{1}

^{1}

Aiming at the problem that ICA can only be confined to the condition that the number of observed signals is larger than the number of source signals; a single channel blind source separation method combining EEMD, PCA and RobustICA is proposed. Through the eemd decomposition of the single-channel mechanical vibration observation signal the multidimensional IMF components are obtained, and the principal component analysis (PCA) is performed on the matrix of these IMF components. The number of principal components is determined and a new matrix is generated to satisfy the overdetermined blind source separation conditions, the new matrix input RobustICA, to achieve the separation of the source signal. Finally, the isolated signals are respectively analyzed by the envelope spectrum, the fault frequency is extracted, and the fault type is judged according to the prior knowledge. The experiment was carried out by using the simulation signal and the mechanical signal. The results show that the algorithm is effective and can accurately diagnose the location of mechanical fault.

In the traditional signal processing, it is generally necessary to know in advance some of the prior knowledge of the signal or the mathematical model of the signal mixing matrix, and then estimate the source signal by filtering or transforming. In practice, the prior knowledge of the signal is not easy to obtain, signal processing cannot solve the problem. The advantage of blind source separation is that it can use the least prior knowledge to obtain the greatest information, the problem originally comes from the cocktail effect, in the cocktail meeting, noisy voices and background music is great, but if people talk to a conversation Interested, still can hear this voices.

The classical blind separation algorithm can recover the source signal better when the number of observed signals is larger than the number of source signals, FastICA algorithm, RobustICA algorithm, second order blind identification algorithm, joint approximation diagonalization and so on. In case of the underdetermined condition, the above algorithm cannot solve the problem. For the underdetermined model, new algorithms are proposed. Based on the blind separation of time-frequency distribution, sparse signal analysis is the main method to solve this problem.

Single channel blind separation is an extreme condition of the underdetermined condition, that is, only through single channel observation signal to estimate the multichannel source signals, in real life, due to environmental or cost constraints, often encountered such extreme problems. In this case, some scholars decompose the signal with wavelet, the resulting signal component is subjected to ICA processing, finally get the source signals [

Based on the existing methods, EEMD, PCA and RobustICA are applied to the mechanical fault diagnosis. Firstly, the EEMD algorithm is used to increase the dimension of the single observation signal, and the number of the source signal is estimated by PCA algorithm. Then, the signal is processed by RobustICA. The experimental results show that the method can effectively isolate the mechanical fault of each part.

1998, Norden E. Huang, first proposed the intrinsic Mode Function (Instrinc Mode Function, the IMF) and its Decomposition method, the concept of Empirical Mode Decomposition (EMD) [

EMD is a new algorithm in the field of modern signal processing. With the popularity of the algorithm and a wide range of applications, its inherent defects are inevitably exposed. Modal aliasing is the most important and deadly defect in empirical mode decomposition. It refers to the fact that the signal is interrupted during the decomposition process due to the weak interference of the signal to be decomposed, so that the adjacent natural modal components are superimposed together, Masking the instantaneous nature of the source signal. In response to this problem, Huang et al. proposed a collective empirical mode decomposition (EEMD) method [

EEMD decomposition principle is as follows:

1) In the single channel mixed signal

2) Decomposes the signal

3) Repeat the above steps to continue adding the same Gaussian white noise (eg. K th decomposition) during the repetition process:

Perform EMD decomposition:

4) Repeat the N times to find the average value of the decomposed groups of intrinsic modal components and the signal margin:

EEMD decomposition results were obtained:

Since the Gaussian white noise with a mean value of zero is added to the single channel mixed signal that needs to be decomposed, after the EEMD decomposition, the intrinsic modal component are averaged, and the white noise in the result will eventually cancel each other, and the noise is eliminated at the same time, and the phenomenon of modal aliasing is avoided.

In order to achieve single channel blind separation, first estimate the number of source signals, that is, after EEMD decomposition with principal component analysis (PCA) to estimate the number of source signals. The purpose of principal component analysis is to find r (r less than n) new vectors, which are used to represent the main features of the entire n-dimensional vectors, thus reducing the dimensionality of the original vector and compressing the entire matrix. Each of the new variables is a linear combination of the original variables, that is, extracted “principal components”, each principal component is uncorrelated, and orthogonal, with practical physical meaning, can represent the original n-dimensional vector of the whole feature. By PCA processing, the n-dimen- sional vector is reduced to r dimension.

For PCA algorithm, how to solve the new vector number r is particularly important. Although the value of r in the algorithm is as small as possible, the smaller the r, the lower the dimension, making the result analysis simpler and less interfering, but this may make some of the key information lost. In order to solve the problem of data loss which may occur in the above algorithm, the solution is to analyze the contribution rate of any vector to information. Contribution rate refers to the proportion of the principal component used in all the data analysis. In determining the r principal component, unless there is a special requirement, the general requirements of the contribution rate to reach more than 85%, because the contribution rate to a certain extent, reflects the size of the reliability, that is, the greater the contribution of the principal component, the greater the reliability.

The blind source separation model is as follows: Assuming that

where

In the modern signal processing process, in most cases the noise can also be considered a class of source signals, or that through other methods to reduce the noise to a negligible level, so the content of this study does not take into account the noise problem, Equation (8) can be rewrite as:

The goal of the blind source separation algorithm is to solve the separation matrix

The key problem with blind source separation is to find the solution matrix W [

RobustICA uses the independent component analysis algorithm based on kurtosis and optimal step length to search the global optimal step by using the sentinel as the control function, find the solution matrix

The method of single channel blind separation based on empirical mode decomposition and principal component analysis is the process of recovering the source signal from the observed signal. Combined with EEMD, PCA and RobustICA, the problem of single channel blind separation and the problem of source number estimation are solved. The process is as follows:

1) The single channel observation signal x is decomposed by EEMD to obtain the IMFs components.

2) Using pca to reduce the IMFs component dimension, Select several elements whose contribution rate is 95% to constitute a new signal. The dimension of this multidimensional signal is the number of source signals estimated by the PCA.

3) The new multi-dimensional signal is processed by RobustICA, and the blind source is separated to obtain the separated source signal [

Here, the experimental signals are selected as sinusoidal signal

where

The source signal

PCA on this series of IMF components to reduce the dimension, select the main component whose contribution rate is 95% to constitute a new signal. The dimension of this multidimensional signal was selected 3 by MATLAB, and the reconstruction signal is shown in “

Regard the reconstructed signals as new observation signals, then use RobustICA to isolate the source signals, The results are shown in “

Calculated by MATLAB, the correlation coefficient between each signal in

“

The experimental data are from the electrical engineering laboratory of Case Western Reserve University. The bearing fault types include inner ring fault, outer ring fault and rolling element fault. The inner ring fault and outer ring fault are selected as the fault source signal from the bearing state. The source signal is linearly mixed with a 1 × 2 random matrix A to obtain a mixed fault signal, i.e. a single channel observation signal. In

The single-channel fault observation signal is processed by the EEMD-PCA- RobustICA method described in the previous chapter, and the resulting separation signal is shown in

the separation signal 1 can be diagnose as outer fault, and the separation signal 2 can be diagnose as inner fault, which verifies that the algorithm can accurately determine the type of mechanical fault [

In this paper, the EEMD-PCA-RobustICA method is proposed, and the simulation signal and the actual mechanical fault signal are used to experiment. The source signal is well separated from the single channel observation signal, which proves the effectiveness of the method. At the same time, the method used in mechanical fault diagnosis, can effectively determine the type of mechanical fault.

Xu, W. and Yan, X.Z. (2017) Application of Single Channel Blind Separation Algorithm Based on EEMD- PCA-RobustICA in Bearing Fault Diagnosis. Int. J. Communications, Network and System Sciences, 10, 138-147. https://doi.org/10.4236/ijcns.2017.108B015