Acoustical signals from mechanical systems reveal the operating conditions of mechanical components and thus benefit for machinery condition monitoring and fault diagnosis. However, the acoustical signals directly measured by the sensors in essential are the mixed signals of all the sources, and normally it is very difficult to be used for source identification or operating feature extraction. Therefore, this paper studies the acoustical source tracing problem using independent component analysis (ICA) and identifies the sources using correlation analysis: the measured acoustical signals are separated into independent components by independent component analysis method, and thus all the independent information of all the sources is obtained; these independent components are identified based on the prior information of the sources and correlation analysis. Therefore, all the source information contained in the measured acoustical signals can be independently separated and traced, which can provide more purer source information for condition monitoring and fault diagnosis.
1. Introduction
Vibration and acoustical signals caused by the collision and friction of mechanical components normally can reduce the operational precision and even shorten the service life of the machinery. However, the vibration and acoustical signals provide important information of the operating conditions, and thus machinery condition monitoring and fault diagnosis can be carried out based on a nondestructive measurement. As the measured acoustical signals are normally mixed signals of all the sources and noises, they have great significance to trace the source information from the measured signals and thus provide pure and independent information of mechanical components for an effective condition monitoring or fault diagnosis.
Generally, the acoustical signals directly measured by the sensors are complicated and rough information of the mechanical systems because they are mixed signals of all the sources and the transmission effects of the transmission paths. Recently, many researchers have devoted their efforts to transmission effects of vibration and acoustical signals. Xin et al. [1] investigated the vibroacoustic performance of a rectangular double-panel partition with enclosed air cavity and simply mounted on an infinite acoustic rigid baffle. Huang and Nutt [2] presented an analytical study of sound transmission through unbounded panels of functionally graded materials. Bravo et al. [3] analyzed the sound absorption and transmission properties of microperforated panels backed by an air cavity and a thin plate. Qian et al. [4] described an approach on the prediction of sound transmission loss for a finite sandwich panel with honeycomb core. Posson et al. [5] developed an analytical model for the sound transmission through an annular stator row in a configuration without mean flow. Lee and Kim [6] investigated the effects of structure on sound absorption and transmission loss of a composite sheet. Chandra et al. [7] analyzed the vibroacoustic and transmission loss characteristics of functionally graded material plates. All these articles studied the transmission characteristics of structures and can help to carry out a passive vibration and sound monitoring and control. However, normally it is a challenging task and costs plenty of time to build a precise model of acoustical transmission for complicated mechanical systems, and thus the engineering applications are limited.
To reveal the operating conditions of the mechanical components, signal processing methods are developed to extract the features hidden in the noisy response signals and thus benefit for system analysis, condition monitoring, or fault diagnosis. To clearly reveal the operating conditions of each mechanical component, blind source separation (BSS) and independent component analysis (ICA) are developed to separate the mixed signals into uncorrelated components without the knowledge of the sources and their mixing modes. Jutten and Herault [8] solved blind source separation problem with an adaptive algorithm based on neuromimetic architecture. Belouchrani et al. [9] proposed a technique to blindly separate the sources using high-order statistics. Comon [10] reexamined the concept of independent component analysis (ICA) and proposed an efficient ICA algorithm. Hyvärinen [11, 12] proposed a very famous ICA algorithm called fast-fixed algorithm for independent component analysis, which is very effective and efficient for linear systems. Since then, independent component analysis has attracted a worldwide attention and has been widely used in many engineering applications, such as image processing [13], EEG/MEG analysis [14], communications [15], radar [16], and modal parameter estimation [17]. In recent years, independent component analysis has been applied to mechanical signal processing, and some practical algorithms have been constructed. Antoni [18] addressed the issues of blind separation of vibration components. Cheng et al. studied the source number estimation [19], source separation [20], and source contribution evaluation [21] methods for mechanical systems based on an effective source separation. Henriquez et al. [22] provided a review of automatic fault diagnosis using independent component analysis and vibration signals. In this paper, independent component analysis is applied to acoustical source tracing (acoustical source separation and identification), and the separation performances of independent component analysis for acoustical signals are quantitatively evaluated by a numerical case study and an experimental study on a mechanical system with shell structures.
The remainder of this paper is organized as follows. In Section 2, basic theory and key principals of independent component analysis are introduced. In Section 3, the separation performance of the independent component analysis algorithm is tested by a numerical case study. In Section 4, the independent component analysis algorithm is applied to separate the acoustical signals, and then acoustical source tracing is carried out based on the source separation and correlation analysis. In Section 5, the conclusions of this study are summarized.
2. Theory of Independent Component Analysis2.1. Principal Theory of Blind Source Separation
Blind source separation (BSS) including independent component analysis (ICA) is a powerful tool to separate a set of source signals S(t)=s1(t),s2(t),…,sn(t)T from the mixed signals X(t)=x1(t),x2(t),…,xm(t)T without any information of the sources and their mixing approaches [8, 23, 24], which can be described as follows: (1)xit=∑j=1naijsjt,i=1,2,…,m,j=1,2,…,n,X=ASMatrixform,where m is the number of the mixed signals, n is the number of the sources, xi is the ith mixed signal collected at location i, sj is the jth source signal, aij is a mixing coefficient, and A={aij} is a mixing matrix.
The key process of BSS or ICA is to find a separating matrix W=ω1,ω2,…,ωnT and recover the estimated sources Y(t)=y1t,y2t,…,yntT by(2)Y=WX=WAS.
2.2. Measures of Independence
To effectively separate the mixed signals into independent components, some quantitative measures of the independence for a random variable should be developed. A well-known theory of independence is “non-Gaussian is independent,” and a classical measure of non-Gaussianity is kurtosis or the fourth-order cumulant. The kurtosis of y (zero-mean and unit-variance) is classically defined by(3)kurty=Ey4-3Ey22.
The main reason of kurtosis as a widely used measure of independence is its computational and theoretical simplicity. Computationally, kurtosis can be estimated simply by using the fourth moment of the sample data. Theoretically, kurtosis has good performances for linear transform: (4)kurty1+y2=kurty1+kurty2,kurtαy=α4kurty,where α is a scalar.
Another and very important measure of non-Gaussianity is given by negentropy, which is based on the information theoretical quantity of entropy. The more “random” the variable is, the larger its entropy will be. Entropy is normally defined as (5)Hy=-∫fylogfydy,where f(y) is the density of a random vector y.
Negentropy is used to obtain an effective measure of non-Gaussianity, and can be defined as follows: (6)Jy=Hygauss-Hy,where ygauss is a Gaussian variable of the same covariance matrix as y.
The advantage of using negentropy as an effective measure of non-Gaussianity is that it is well justified by statistical theory. However, the estimation of negentropy is difficult, and therefore some approximations have to be used to simplify the calculation. The classical method of approximating negentropy is using higher-order moments: (7)Jy≈112Ey32+148kurty2.
2.3. Efficient ICA Algorithm2.3.1. Centering
The most basic and necessary process before an effective source separation is to center x:(8)x=x-meanx.
2.3.2. Whitening
One popular method for whitening is to use the eigenvalue decomposition (EVD) of the covariance matrix E{xxT}=EDET, where E is the orthogonal matrix of eigenvectors of E{xxT} and D is the diagonal matrix of its eigenvalues, D=diag(d1,…dn): (9)x~=ED-1/2ETx,where D-1/2=diag(d1-1/2,…,dn-1/2).
2.3.3. Fixed Point Iteration Scheme for One Component
Choose an initial weight vector w.
Let w+=E{xg(wTx)}-E{g′wTx}w, where g is the derivatives of the nonquadratic function and normally has (10)g1u=tanhau,g2u=uexp-u22,
where 1≤a≤2 is a suitable constant.
Let w=w+/w+.
If not converged, go back to (2).
2.3.4. Fixed Point Iteration Scheme for Several Components
Let wp+1=wp+1-∑j=1pwp+1Twjwj.
Let wp+1=wp+1/wp+1Twp+1.
Comparing the separating vector w for one component, the separating matrix W=w1,…,wnT for several components can be obtained as follows: (11)W=WWWT,W=32W-12WWTW.
2.4. Source Identification
For mechanical systems, some information of the sources normally can be obtained by theory study or instructions. Therefore, waveform correlation can be used to identify the source information by a correlation analysis between the sources and the separated signals. For discrete signals s and y, waveform correlation coefficient ρW is defined as(12)ρW=maxRsy(k)maxRss(k)×maxRyy(k),0≤k≤T,where Rsy(k)=∑i=0T-k-1s(i+k)y(i), Rss(k)=∑i=0T-k-1s(i+k)s(i), Ryyk=∑i=0T-k-1yi+kyi, ∀k≥0, and T is the data length.
3. Numerical Case Study
In this section, typical vibration and acoustical signals of mechanical systems are artificially generated to test the separation performance of the ICA algorithm based on negentropy. Since acoustical signals transmit from the sources to the measured points through air, the mixing mode of all the acoustical source signals tends to be a linear superposition. Furthermore, the generating mechanism of different sources is far from each other. Therefore, all the sources are considered to be linearly mixed and independent of each other.
The given source signals are typical vibration and acoustical signals of mechanical systems: signal s1(t) is a sinusoidal wave that simulates a vibration signal of rotational equipment; signal s2(t) is a periodic wave of oscillating attenuation that simulates mechanical shocks; signal s3(t) is a periodic wave that simulates amplitude modulation; signal s4(t) is a periodic wave that simulates frequency modulation; and signal s5(t) is a white noise that simulates noises produced by environment and structural transmission. The generating functions of sources are as follows:(13)St=s1(t)s2(t)s3(t)s4(t)s5(t)=sin15πt∑sin1000t-0.1n-0.02exp-150t-0.1n-0.022ut-0.1n-0.02sin10πtcos150πtsin6πcos6πtnt,where u(t) is a step function.
Generally, the measured acoustical signals are mixed signals of all the sources with noises. Furthermore, as the acoustical signals transmit from the sources to the measuring points through air, linear superposition is applied to artificially produce the mixed signals. The mixing matrix A is randomly generated as the following matrix:(14)A=0.23-0.580.73-0.810.950.28-0.890.49-0.750.970.560.15-0.720.580.99-0.38-0.520.450.610.99-0.480.380.84-0.510.980.88-0.610.54-0.200.95.
The waveforms of the source signals and the mixed signals are shown in Figures 1 and 2. Obviously, all the source signals have typical waveform features. However, it is very difficult to directly identify the features of the mixed signals as the complicated waveforms. Therefore, generally signal processing methods are required to help to extract the desired features clearly. In this numerical case study, the ICA algorithm based on negentropy is applied to separate the mixed signals into 5 independent components, and then these separated components are identified by the correlation analysis.
Waveforms of the source signals.
Waveforms of the mixed signals.
With the ICA algorithm, 6 mixed signals are separated into 5 independent components, and the waveforms of these independent components are shown in Figure 3. Comparing Figure 3 with Figure 1, obviously the waveform features of each source signal are well separated. Waveform correlation analysis is used to quantitatively evaluate the separation performances of the ICA algorithm and identify the sources, and the waveform correlation coefficients ρW are shown in the following matrix:(15)ΩW=0.990.010.040.000.040.010.990.030.170.120.030.030.990.010.090.000.150.010.990.030.050.120.090.040.99.
Separated components by DSS.
Each correlation coefficient ρij in ΩW={ρij} indicates the similarity between the separated component i and the source j. From the correlation matrix ΩW={ρij}, the elements in the diagonal show the correlation coefficients between the separated components and related source signals: the correlation coefficient of s1 (s2, s3, s4, and s5) and related separated component 1 (2, 3, 4, and 5) is 0.99, which indicates that all the source information of s1 (s2, s3, s4, and s5) has been completely separated from the mixed signals by the ICA algorithm. The high correlation coefficients indicate that the ICA algorithm can effectively separate the mixed signals into independent components for the given sources and linearly mixed mode. However, the correlation coefficients between the independent components and the unrelated sources are all less than 0.17, which indicates that all the separated components have good independent properties. Therefore, the ICA algorithm based on negentropy is applied to extract the source information from the measured acoustical signals from real mechanical systems, and the separated information is used to validate the effectiveness of ICA algorithm and trace the source information hidden in the measured signals.
4. Experimental Study4.1. Introductions of the Test Bed
A test bed of a mechanical system with shell structures is constructed to test the separation performances of the ICA algorithm based on negentropy, which is composed of four components: an end cover, a shell, two clapboards, and supports. Rubber air springs are used to support the whole test bed and reduce the effects of the ground vibrations and environmental noises. Three acoustical sources are designed in the test bed: two loudspeakers controlled by the signal generators and one motor controlled by the frequency converter. The structure of the test bed is shown in Figure 4.
The structure of the test bed.
Six sound pressure sensors are used to measure the acoustical signals, and they are located in six directions of the test bed with a distance of 500 millimeters. HBM Gen2i data acquisition system is applied to collect the acoustical data from these six sensors. The framework of the measuring system is shown in Figure 5, and the test parameters are shown in Table 1.
The test parameters of the measuring system.
Parameters
Values and units
Sound pressure sensors
6
HBM Gen2i data acquisition system
1
Test bed with shell structure
1
Sampling frequency
10240 Hz
Sampling length
10 seconds
Rotational speed of motor
1800 r/min (20 Hz)
Frequency of Loudspeaker I
1750 Hz
Frequencies of Loudspeaker II
500, 1700, and 2000 Hz
The framework of the measuring system.
4.2. Acoustical Signals of the Test Bed
The acoustical signals measured by the six sound pressure sensors are the mixed signals of all the sources as they are working together. In the experimental study, the acoustical signals from sensors 1, 4, and 6 are used as the mixed signals so as to satisfy the assumptions of ICA: the number of the mixed signals should be no less than the number of the source signals. Furthermore, the directions of these sensors represent a diversified mixing mode of all the sources.
The waveforms and spectrums of the mixed signals are shown in Figures 6 and 7, respectively. From waveforms in Figure 6, it is difficult to identify the waveform features of the source signals except some periodic waves, which indicates that the waveforms of the mixed signals are complicated, and normally signal processing method is required to extract the desired features. From spectrums in Figure 7, some major components of 20, 1100, 1700, 1750, and 2000 Hz are clear, which represent the characteristic features of the sources as the experimental settings. Generally, the independent information of the source signals cannot be directly identified from the measured mixed signals as each mixed signal contains all the major components of the source signals.
Waveforms of the mixed signals.
Spectrums of the mixed signals.
4.3. Acoustical Source Separation
The ICA algorithm based on negentropy is applied to separate the mixed signals into independent components, and 3 independent components are extracted from the given mixed signals. The waveforms and spectrums of the independent components are shown in Figures 8 and 9, respectively. Figure 8 clearly shows that the waveform of the separated component y1 has typically periodic and amplitude-modulation features, and the basic components are sine waves, which are normally caused by the eccentric vibration of mechanical systems. The spectrums in Figure 9 also show that the separated component y1 has major components of 20 and 2000 Hz. The waveform of the separated component y2 has typical features of sine waves, and its spectrum also clearly shows the major component of 1750 Hz. The waveform of the separated component y3 has an obvious feature of sine waves with a periodic amplitude-modulation, and its spectrum contains three major components of 500, 1700, and 2000 Hz.
Waveforms of separated components.
Spectrums of separated components.
Comparing the spectrums of the separated components with the parameters of the experimental settings, it can be speculated that the separated component y1 represents the typical feature of the source 1 from a motor, while the separated components y2 and y3 represent the typical features of the sources 2 and 3 from Loudspeakers I and II, respectively. However, this is just based on the parameters of the experimental settings and it is still not a convincing and reliable source tracing method.
4.4. Acoustical Source Tracing and Validation
To intelligently trace the sources and validate the effectiveness of the ICA algorithm in real mechanical systems, acoustical source signals are measured independently by the closest sensors in the condition that only one source is working with the given parameters. The independent source waveforms from sensors 1, 4, and 6 in the condition that only the motor, Loudspeaker I, or Loudspeaker II is working with the given experimental settings are shown in Figure 10, and their spectrums are shown in Figure 11.
Waveforms of the source signals.
Spectrums of the source signals.
Comparing the waveforms and spectrums of the source signals with those of the independent components, the waveforms of the independent components are very similar to those of the related source signals: the independent component y1 and the source signal s1 have a typical feature of sine waves with an uncertainty amplitude modulation, which can be caused by the eccentric vibration of the motor, and the uncertainty amplitude modulation can be caused by the rubbing; the independent component y2 and the source signal s2 have the same major component of 1750 Hz, which agrees with the experimental settings; the independent component y3 and the source signal s3 have three same major components of 500, 1700, and 2000 Hz as the experimental settings. Therefore, all the major components of the source signals can be effectively traced by the ICA algorithm based on negentropy.
To quantitatively and intelligently trace the sources, the waveform correlation analysis is used to evaluate the similarity between independent components and the sources. All the independent components are made a correlation analysis with all the sources, and the correlation coefficients are listed in the correlation matrix ΩW: (16)ΩW=0.900.070.060.050.760.080.100.160.85.
The correlation matrix ΩW shows that the correlation coefficients between the independent components and the related sources are 0.90, 0.76, and 0.85, which indicate high correlation coefficients and high similarity between the independent components and the related sources (Liu et al. [25] obtained waveform correlation coefficients of 0.77±0.03 for ECG signals with noises, and Farina et al. [26] obtained correlation coefficients of 0.70±0.09 for nonstationary surface myoelectric signals), while the correlated coefficients between the independent components and the unrelated sources are less than 0.16, which indicates that all the independent components and all the sources have good independence property. Therefore, a threshold γ can be set as γ∈(0.16,0.76) (in practice γ∈(0.50,0.70)) to intelligently identify and trace the acoustical sources, and the high correlation coefficients between the independent components and the related sources validate the effectiveness of the ICA algorithm based on negentropy in source separation and acoustical source tracing.
5. Conclusions
This paper presents the fundamental theory of blind source separation and key principals of independent component analysis and validates the effectiveness of the ICA algorithm based on negentropy according to a numerical case study and an experimental study on a mechanical system with shell structures. The experimental study indicates that the acoustical sources can be effectively separated and intelligently traced.
In the numerical case study, five typical acoustical source signals of mechanical systems are effectively separated from six linearly mixed signals, and the correlation coefficients between the independent components and the related source signals are all more than 0.99, which indicates a complete and highly effective source separation of the ICA algorithm for the given mixed signals while in the experimental study on a mechanical system with shell structures, the correlation coefficients between the independent components and the related source signals are all more than 0.76, which also reveals an effective acoustical source separation. If artificially giving a threshold γ∈(0.50,0.70) for the correlation coefficients, all the acoustical sources can be intelligently identified and traced.
This work can provide pure source information for machinery condition monitoring and fault diagnosis, and the independent source information can also be of benefit for noise identification, reduction, and control.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the projects of National Nature Science Foundation of China (no. 51305329, 51035007), the China Postdoctoral Science Foundation (nos. 2013M532032 and 2014T70911), the Doctoral Foundation of Education Ministry of China (no. 20130201120040), and the Shaanxi Postdoctoral Scientific Research Project.
XinF. X.LuT. J.ChenC. Q.Sound transmission through simply supported finite double-panel partitions with enclosed air cavity2010132110.1115/1.40004662-s2.0-77955265522HuangC. Z.NuttS.An analytical study of sound transmission through unbounded panels of functionally graded materials201133061153116510.1016/j.jsv.2010.09.0202-s2.0-78650714523BravoT.MauryC.PinhedeC.Sound absorption and transmission through flexible micro-perforated panels backed by an air layer and a thin plate2012131538533863QianZ. C.ChangD. Q.LiuB. L.LiuK.Prediction of sound transmission loss for finite sandwich panels based on a test procedure on beam elements2013135610.1115/1.40238422-s2.0-84889065706PossonH.BériotH.MoreauS.On the use of an analytical cascade response function to predict sound transmission through an annular cascade2013332153706373910.1016/j.jsv.2013.02.0132-s2.0-84876929919LeeB. C.KimS. R.Effect of structure on sound absorption and sound transmission loss of composite sheet201423431932510.1080/09243046.2013.8687122-s2.0-84898001284ChandraN.RajaS.GopalK. V. N.Vibro-acoustic response and sound transmission loss analysis of functionally graded plates2014333225786580210.1016/j.jsv.2014.06.031JuttenC.HeraultJ.Blind separation of sources, part I: an adaptive algorithm based on neuromimetic architecture199124111010.1016/0165-1684(91)90079-x2-s2.0-0026191274BelouchraniA.Abed-MeraimK.CardosoJ.-F.MoulinesE.A blind source separation technique using second-order statistics199745243444410.1109/78.5543072-s2.0-0031078854ComonP.Independent component analysis, a new concept?199436328731410.1016/0165-1684(94)90029-92-s2.0-0028416938HyvärinenA.Fast and robust fixed-point algorithms for independent component analysis199910362663410.1109/72.7617222-s2.0-0032629347HyvärinenA.OjaE.Independent component analysis: algorithms and applications2000134-541143010.1016/s0893-6080(00)00026-52-s2.0-0042826822TsaiD.-M.LaiS.-C.Independent component analysis-based background subtraction for indoor surveillance200918115816710.1109/TIP.2008.2007558MR27417412-s2.0-58049211805HyvärinenA.RamkumarP.ParkkonenL.HariR.Independent component analysis of short-time Fourier transforms for spontaneous EEG/MEG analysis201049125727110.1016/j.neuroimage.2009.08.0282-s2.0-70350070225GaoJ.ZhuX.NandiA. K.Independent component analysis for multiple-input multiple-output wireless communication systems201191460762310.1016/j.sigpro.2010.08.0082-s2.0-78650857272JiJ.Robust approach to independent component analysis for SAR image analysis20126328429110.1049/iet-ipr.2009.0084MR29624242-s2.0-84859031753Rodríguez-RochaR.Rivero-AngelesF. J.Gómez-RamírezE.Application of the baseline stiffness method for damage detection of an RC building without baseline modal parameters utilizing the independent component analysis for modal extraction2013133125007410.1142/s02194554125007452-s2.0-84875690512AntoniJ.Blind separation of vibration components: principles and demonstrations20051961166118010.1016/j.ymssp.2005.08.0082-s2.0-27644480038ChengW.LeeS.ZhangZ. S.HeZ. J.Independent component analysis based source number estimation and its comparison for mechanical systems2012331235153516710.1016/j.jsv.2012.06.0212-s2.0-84864971737ChengW.HeZ. J.ZhangZ. S.A comprehensive study of vibration signals for a thin shell structure using enhanced independent component analysis and experimental validation2014136404101110.1115/1.4027545ChengW.ZhangZ. S.LeeS.HeZ. J.Source contribution evaluation of mechanical vibration signals via enhanced independent component analysis2012134202101410.1115/1.40058062-s2.0-84859903676HenriquezP.AlonsoJ. B.FerrerM. A.TraviesoC. M.Review of automatic fault diagnosis systems using audio and vibration signals201444564265210.1109/tsmcc.2013.22577522-s2.0-84899744851ComonP.JuttenC.HeraultJ.Blind separation of sources. Part II: problems statement1991241112010.1016/0165-1684(91)90080-32-s2.0-0026188054SorouchyariE.Blind separation of sources. Part III. Stability analysis1991241212910.1016/0165-1684(91)90081-s2-s2.0-0026188053LiuH.ChangC. Q.LukK. D. K.HuY.Comparison of blind source separation methods in fast somatosensory-evoked potential detection201128217017710.1097/wnp.0b013e31821213bd2-s2.0-79955054423FarinaD.FévotteC.DoncarliC.MerlettiR.Blind separation of linear instantaneous mixtures of nonstationary surface myoelectric signals20045191555156710.1109/tbme.2004.8280482-s2.0-1642415267