Radiation Noise Separation of Internal Combustion Engine Based on Gammatone-RobustICA Method

In the internal combustion engine noise source separation process, the combustion noise and the piston slap noise are found to be seriously aliased in time-frequency domain. It is difficult to accurately separate them. Therefore, the noise source separation method which is based on Gammatone filter bank and robust independent component analysis (RobustICA) is proposed. The 6-cylinder internal combustion engine vibration and noise test are carried out in a semianechoic chamber. The lead covering method is adopted to isolate the interference noise from numbers 1 to 5 cylinder parts, with only the number 6 cylinder parts left bare. Firstly, many mode components of the measured near-field radiated noise signals are extracted through the designed Gammatone filter bank. Then, the RobustICA algorithm is utilised to extract the independent components. Finally, the spectrum analysis, the continuouswavelet time-frequency analysis, the correlation functionmethod, and the drag test are employed to further identify the separation results. The research results show that the frequency of the combustion noise and the piston slap noise are, respectively, concentrated at 4025Hz and 1725Hz. Compared with the EWT-RobustICA method, the separation results obtained by the Gammatone-RobustICA method have very fewer interference components.


Introduction
The internal combustion engine as power heart is widely used in ships, vehicles, and other means of transportation.With the development of internal combustion engine to heavy load, high speed, and light direction, its vibration noise problem is becoming more and more serious, which has become the main vibration noise source of ships and vehicles.By reducing the level of vibration and noise of the internal combustion engine, the working performance of the internal combustion engine can be improved, and the people's living environment can be ameliorated [1].In order to draw up the vibration and noise reduction program of the internal combustion engine, the first step is to accurately separate and identify the noise sources of the internal combustion engine.Many parts of the internal combustion engine can produce noise.According to the source of the internal combustion engine, the noise can be divided into mechanical noise, combustion noise, and aerodynamic noise [2][3][4].Mechanical noise mainly includes piston slap noise, air valve knock noise, gear meshing noise, and fuel injection pump noise.Combustion noise is caused by the change of the cylinder pressure.Aerodynamic noise includes intake noise, exhaust noise, and fan noise.
Currently, there are many scholars who utilised the cyclic Wiener filtering method [5], blind source separation method [6][7][8], improved spectrofilter method [9], speed-varying filter [10], and other multichannel methods to separate and identify the noise sources of the internal combustion engine.
When the multichannel method is adopted to separate and identify the noise sources of the internal combustion engine, the multiple channel noise signals of the internal combustion engine need to be measured.However, in the actual engineering test, by the limit of cost and installation conditions, it is necessary to use a few sensors to achieve the same noise source separation effect.Therefore, scholars studied the single channel method to separate and identify the noise sources of internal combustion engine.For example, Zhang et al. [11] adopted the ensemble empirical mode decomposition, the coherent power spectrum analysis, and the improved analytic hierarchy process to separate and identify the noise sources of the diesel engine.Bi et al. [12] utilised the EEMD-RobustICA method to separate the combustion noise, piston slap noise, and exhaust noise of the gasoline engine.
When the above single channel method is employed to separate and identify the noise source of the internal combustion engine, the first step is to decompose the single channel noise signal by the EEMD algorithm.However, EEMDbased single channel method has many defects [13][14][15].The EMD algorithm lacks rigorous mathematical derivation.The EEMD algorithm is required to add Gaussian white noise before each step of the EMD algorithm.Therefore, the calculation cost of the EEMD algorithm is very large.Moreover, when the EEMD algorithm is applied to decompose the noise signal, the endpoint effect and modal aliasing problems will also exist.In recent years, empirical wavelet transform (EWT) [16], variational mode decomposition (VMD) [17], and the related improved methods are proved to be useful and better tools to decompose the signal.For example, Yao et al. [18] utilised variational mode decomposition and robust independent component analysis to separate the noise source of diesel engine.Moreover, these methods are widely used in bearing fault diagnosis [19][20][21].
In the real world, the human auditory system can distinguish the mixed speech signals in a noisy environment.At present, there are many scholars establishing the computational auditory scene analysis model and algorithm based on the human ear hearing system to separate and identify the mixed speech signals [22].The speech signals and the radiated noise signals of the internal combustion engine belong to the separable aliasing sound signals.Thus, the noise source of the internal combustion engine can be separated and identified by the human ear hearing model.Moreover, the Gammatone filter bank is a widely used human ear hearing model [23][24][25][26].Therefore, the Gammatone hearing filter bank can be introduced into the separation and identification of the internal combustion engine noise sources.
In this paper, the lead covering method is utilised to isolate the noise interference from 1 to 5 cylinders, with only the number 6 cylinder part left bare.The WP10-240 sixcylinder four-stroke high speed internal combustion engine vibration and noise test are carried out in a semianechoic chamber.The Gammatone-RobustICA method is proposed to separate and identify the combustion noise and the piston slap noise of the number 6 cylinder.
The paper is organized as follows.In Section 2, basic theory methods (such as the cochlear basement membrane model, the Gammatone filter bank, and the robust independent component analysis) are described.In Section 3, the simulation analysis is carried out.In Section 4, an internal combustion engine noise and vibration test is introduced.In Section 5, separation and identification methods are adopted to separate and identify the noise sources of an internal combustion engine.Finally, Section 6 presents conclusions.

Cochlear Basement Membrane
Model of the Auditory System.The human auditory system includes the auditory periphery and the auditory center [27,28].The auditory periphery consists of outer ear, middle ear, and inner ear.The inner ear consists of semicircular canal, vestibule, and cochlear implant.The speech signals are through the outer ear, middle ear, and ear bone to the cochlear.The cochlear transforms the mechanical energy of sound waves to neural coding signals and through the auditory nerve into the auditory center.Then the hearing will be produced.
In the auditory system, the cochlear as the human auditory receptor is one of the most important parts.The sound signal processing of the human ear is achieved by the frequency decomposition of cochlear basement membrane.It can not only convert different frequencies of the sound signals into different basement membrane positions, but also switch different sound intensity into different basement membrane vibration amplitude.Thus, the cochlear can complete the coding according to the sound frequency and the sound intensity.
In the auditory model, a set of overlapping band-pass filters' bank is usually used to simulate the cochlear basement membrane.Suppose there are  filters.The th filter is ℎ(, ).The response of the basement membrane to the signal () is defined as follows.
where (, ) is the output of the basement membrane.* represented the convolution.The characteristic of the Gammatone filter bank is highly consistent with the auditory characteristic which can well simulate the basement membrane frequency selectivity and spectral analysis characteristics.Thus, the Gammatone filter bank is widely adopted as the cochlear basement membrane model.

Gammatone Filter Bank.
The Gammatone filter bank is a standard cochlear auditory filter [29,30].It can simulate the band-pass filter characteristics of the basement membrane.The center frequency of each filter is logarithmically evenly distributed on the frequency axis.Suppose the center frequency of the th filter is   .The time domain expression of the Gammatone filter is where  is the filter order.  is the filter center frequency.  is the phase of the filter, due to the phase of the sound signals having small influence for hearing; thus the phase   is usually set as zero.() is the unit step function. is the attenuation factor of filter.It determines the decay rate of the impulse response which is related to the bandwidth of the filter.The bandwidth of the auditory band-pass filter depends on the center frequency.It usually adopted the equivalent rectangular bandwidth (ERB).The calculation formula is  = 1.019ERB (  ) = 1.019 (24.7 + 0.108  ) . ( The Gammatone filter bank is shown in Figure 1.The time domain waveform under different center frequencies   is shown in Figure 1(a).The corresponding amplitude frequency response curve is shown in Figure 1    x n (t) y n (t)

Robust Independent Component Analysis.
The blind source separation method is a powerful signal processing method [31].It can extract and recover the original signals which cannot be directly observed from several observed mixed signals.Due to the fact that the source signal is unknown and unobserved and the mixed system characteristics are unknown or only a small amount of prior knowledge is known in advance, it is a challenge to use the blind source separation method to separate the source signals.
The independent component analysis (ICA) is an important method to solve the blind source separation problems.The FastICA algorithm is a widely used blind source separation algorithm which is proposed by the Hyvärinen [32].The FastICA algorithm has high convergence speed and reliability.It is cubic convergence.However, it has some limitations.When the source signal spatial correlation is high, the FastICA algorithm will have not good separation results, and the kurtosis based on higher order cumulant has its own defects in the field value.To solve these problems, Zarzoso and Comon [33] proposed the RobustICA algorithm which has excellent blind source separation effect.The RobustICA algorithm utilised the kurtosis as the objective function.The kurtosis is an indicator to evaluate the non-Gaussian characteristic of the analysed signal.The kurtosis is defined as Suppose the observed signal is .The RobustICA algorithm does not need to carry out the whitening process for the observed signal.It only requires the mean value of the observed signal which is zero.By RobustICA algorithm, the unmixed signal is  =   .The specific calculation process of the RobustICA algorithm is shown as follows.
Step 1. Counter  = 1.Suppose the number  of source signals is the same with the number of observed signals.
Step 2.   is randomly assigned as the initial value   (0) and the norm is 1.The number of iteration steps  is set as 1.
Step 7. Make  =  + 1.Then go back to Step 2 for next iteration until  = .

Calculation Process of Gammatone-RobustICA Method.
For the measured single channel radiated noise signals, the Gammatone-RobustICA method is utilised to separate and identify the noise sources.First, the different mode components of the measured signals are extracted by the Gammatone filter bank.Because these different mode components are not always independent of each other, thus the blind source separation technique is still needed to further deal with the separated signal components to extract the independent components.Then the correlation coefficient between the mode components and the measured signal is calculated.
The mode components which have a higher correlation coefficient with the measured signal are retained.The retained mode components and the measured signal are combined together to form a new signal group.The RobustICA algorithm is used to extract the independent components.Finally, the obtained separation results are further identified through the spectral analysis, the continuous wavelet transform (CWT), the coherence function method, and the drag test of internal combustion engine.The calculation flow of the Gammatone-RobustICA method is shown in Figure 3.

Simulation Analysis
In order to illustrate the performance of the Gammatone-RobustICA method, some typical signals are selected to carry out the simulation analysis by MATLAB.The selected typical signals (1, 2, and 3) are shown as follows.
The sampling frequency of the selected signals is 3000 Hz, and the sampling number of the signals is 300.The time domain waveform of the selected source signals is shown in Figure 4.
For the mixed signal , the Gammatone-RobustICA method and the EWT-RobustICA method are, respectively, utilised to separate the selected source signals 1, 2, and 3.
When using the Gammatone-RobustICA method, the frequency range of the Gammatone filter bank is set to 45 Hz-205 Hz, and the channel number is set to 3. By the Gammatone filter bank, the three mode components of the mixed signals can be extracted.Then the RobustICA algorithm is adopted to further extract the independent components from the three mode components.The separation results of Gammatone-RobustICA method are shown in Figure 5(a).From Figures 4 and 5(a), it can be seen that the IC1, IC2, and IC3 are, respectively, corresponding to 1, 2, and 3.The results obtained by Gammatone-RobustICA method are similar to the time domain waveforms of the source signals.
When using EWT-RobustICA method, the separation results of EWT-RobustICA method are shown in Figure 5(b).From Figures 4 and 5(b), it can be seen that IC1, IC2, and IC3 are, respectively, corresponding to 1, 2, and 3.But  the IC2 has some difference with the 2.This can be seen from the red circle in Figure 5(b).Thus, it can be considered that Gammatone-RobustICA method has a better separation effect than EWT-RobustICA method.1.

Experimental Investigation
The internal combustion engine has six cylinders.It will produce many vibration and noise sources.It is very difficult to directly separate all the noise sources of the internal combustion engine.Therefore, the lead covering method is utilised to isolate the noise from the numbers 1-5 cylinders, and only the number 6 cylinder parts are exposed.The noise source of the specified number 6 cylinder is separated and identified.The lead cover of the internal combustion engine is shown in Figure 6.

Test Conditions.
The internal combustion engine is normally operated at a rated speed of 2200 rpm.Thus, it is necessary and meaningful to separate and identify the noise  Drag condition 2100 0 sources of the internal combustion engine at the rated speed.However, in the internal combustion engine test bench, the coupling equipment of the internal combustion engine has some problems.The internal combustion engine cannot reach the rated speed.In the test, the actual speed of the internal combustion engine is 2100 rpm.Moreover, when the internal combustion engine is in drag condition, the internal combustion engine will not produce the combustion noise and will only produce the mechanical noise.However, it is difficult to measure the independent piston slap noise because many moving parts of an internal combustion engine will produce noise.According to the relevant knowledge of internal combustion engine, when the piston impacts the inner wall of the cylinder, it will produce the vibration, and then the vibration will further produce piston slap noise.Thus the frequency of the piston slap vibration can be utilised to assess the accuracy of the separated piston slap noise.The test conditions of internal combustion engine are shown in Table 2.

Measuring System and Measuring Point Arrangement.
The vibration and noise measuring system of internal combustion engine is shown in Figure 7.The measuring system includes the NI 9234 and NI 9205 acquisition module that the highest sampling rate can be up to 51.2 kHz.The LC0158T accelerometers are adopted to measure the cylinder head vibration and the piston slap vibration that the sensitivity is 30 mV/g, the range is 166 g, and the frequency range is 0-15 kHz.The type of the cylinder pressure sensor is Kistler 7013C that the range is 25 MPa, with a single channel charge amplifier 5018A1000.The DGO9767CD electret microphone is applied to measure the noise signals that the sensitivity is 50 mV/Pa, and the frequency response range is 20 Hz-20 kHz.In fact, due to the limitation of the test conditions, it is difficult to measure the independent combustion noise.According to the related knowledge of internal combustion engine, the combustion noise is related to the drastic change of cylinder pressure.Drastic change of cylinder pressure can cause vibration of cylinder head and body surface, and then the vibration will produce the combustion noise.Thus, the correlation function of cylinder pressure and cylinder head vibration can be used to determine the frequency of the combustion noise and then further evaluate the accuracy of separated combustion noise.

Accelerometer for piston slap vibration
Due to the fact that the structure of internal combustion engine and movement trajectory of piston are known before, thus the place of the piston slap occurring can be determined according to the structure of internal combustion engine and movement trajectory of piston.The measured position of the piston slap vibration is corresponding to the position that the piston impacts the inner wall of the cylinder.The measured position of the piston slap vibration is shown in Figure 8.
The DGO9767CD electret microphone is arranged at 1 cm distance away from the number 6 cylinder body side.The LC0158T accelerometer is arranged at the piston slap place  to measure the piston slap vibration when the internal combustion engine is in drag condition.The specific arrangement of measuring point is shown in Figure 9.
In the test, the highest frequency of the internal combustion engine radiated noise is below 8000 Hz.According to the sampling theorem, the sampling frequency is greater than twice the highest frequency of the analysed signals.Thus, the sampling frequency can be set to 25.6 kHz.
The number 6 cylinder is set as the research object.When the internal combustion engine is in 2100 rpm and noload condition, the cylinder pressure signals ( 1 ), cylinder head vibration acceleration signals (), and the cylinder body side near-field radiated noise signals ( 2 ) are measured.It is shown in Figure 10.

Separation and Identification of the Near-Field Radiated Noise
In this part, the Gammatone-RobustICA method and the EWT-RobustICA method are, respectively, used to separate and identify the noise sources of the cylinder body side nearfield radiated noise.
In the process of the noise and vibration test, the measured noise signal may have random error components which will affect the subsequent calculation results.In order to reduce the random error components, the pretreatment such as eliminate trend items and slip average is carried out on the measured noise signal.The preprocessed noise signal is shown in Figure 11.5.1.Gammatone-RobustICA Method.When the Gammatone-RobustICA method is utilised to separate and identify the noise source of the internal combustion engine, the first step is to design an appropriate Gammatone filter bank to extract various mode components of the preprocessed noise signal.It is necessary to predefine the two important parameters of the Gammatone filter bank, the center frequency range, and the number of channels.For the center frequency range, on the one hand, the frequency range of human audible sounds is 20 Hz-20 kHz.On the other hand, the frequency range of the internal combustion engine radiated noise is usually below 8000 Hz.Thus, the center frequency of the Gammatone filter bank is set to 20 Hz-8000 Hz.For the number of channels of the Gammatone filter bank, on the one hand, considering the calculation accuracy and the computing costs, the higher the number of channels is, the higher the accuracy of the calculation is, and the higher the cost of the computation is; thus the number of channels of the Gammatone filter bank should not be too much.On the other hand, according to the related knowledge of internal combustion engine, the noise sources of internal combustion engine are combustion noise, piston slap noise, air valve knock noise, gear meshing noise, fuel injection pump noise, and so forth.After considering these factors, the number of channels of the Gammatone filter bank is set to 11.
By the Gammatone filter bank, the mode components from the preprocessed noise signal can be obtained.The decomposition results are shown in Figure 12.
In order to improve the efficiency of the calculation, the mode components with large correlation with the preprocessed noise signal need to be selected to carry out the next step calculation.
Suppose the correlation coefficient between the mode components and the preprocessed noise signal is   ( = 1, 2, . . ., );  is the number of the mode components.Correlation coefficient  is defined as follows.
where  represents the source signal. represents the separation component.cov(, ) represents the covariance of  and .() and () represent the variance of  and , respectively.
According to Pearson correlation coefficient theory [34], the correlation coefficient is in the range of −1 to +1.If the correlation coefficient is greater than zero, there is a positive correlation between the two variables.On the contrary, there is a negative correlation.In general, if the correlation coefficient is greater than 0.3, the two variables are correlated.In order to select the mode components, the threshold  is set to determine whether the mode component is retained or not.But the internal combustion engine has many noise sources, and there are many noise disturbances in the measurement process.Thus, if the threshold of the correlation coefficient is set to 0.3, the selected mode components will be less, and the separated independent noise source information will not be sufficient.In order to obtain more sufficient independent noise source information, in the actual study the threshold of the correlation coefficient needs to be set smaller.Here, the threshold of the correlation coefficient is defined as follows.

Shock and Vibration
where  is the ratio factor and  is set to 10.0.The calculation results of the correlation coefficient between the mode components and the preprocessed noise signal are shown in Table 3.
From Table 3, it can be seen that the maximum correlation coefficient is 0.1256.Thus, if the correlation coefficient is more than 0.01256, it should be retained.On the contrary, if the correlation coefficient is less than 0.01256, it should be taken out.Through analysing, it can be determined that m2, m3, m4, m5, m6, m7, m9, and m11 should be retained, and m1, m8, and m10 should be taken out.
Due to the fact that the retained mode components are not always independent of each other, the RobustICA algorithm needs to be used to extract the independent components.The retained mode components and the preprocessed  noise signal are combined together to form a new signal group.Then the RobustICA algorithm is employed to extract the independent components.The calculation results are shown in Figure 13.
According to the relevant knowledge of internal combustion engine, the fuel combustion in the cylinder can cause drastic pressure change in the cylinder, and it caused the combustion noise.Thus the combustion noise is high frequency signal.Through analysing the nine components, the IC6 is the high frequency signal.Thus it can be preliminarily judged that the IC6 component is the combustion noise.
The mechanical noise of internal combustion engine is usually below 3500 Hz [35].Through analysing the Robus-tICA calculation results, the frequency of IC9 is below 3500 Hz.Thus it can be preliminarily judged that the IC9 component is the piston slap noise.
Apart from IC6 and IC9, as for the rest 7 components, they may be the air valve knock noise, gear meshing noise, fuel injection pump noise, or other noise sources.Ideally, all the noise sources of the internal combustion engine should be separated and identified.However, since the noise sources are seriously mixed with each other, it is difficult to separate and identify all the noise sources, and they need further research.For example, the air valve knock noise is mainly caused by the impact of valve opening and closing, and timing of valve opening and closing can be determined according to the timing of inlet and exhaust of internal combustion engine.Moreover, through arranging the vibration sensor and the acoustic sensor near the valve, and combining the timing of valve opening and closing, the air valve knock noise can be further studied.The gear meshing noise is generated by the collision and friction between teeth during the gear meshing process.Every turn of the gear can produce the collision and friction; thus the gear meshing noise can be further studied.The fuel injection pump noise is mainly caused by the oil pressure change during the oil spraying process.The fuel injection pump noise can be studied through oil injection pressure and burning time of internal combustion engine.
Considering that combustion noise and piston slap noise are the main noise sources of internal combustion engine, thus the paper is mainly on the separation and identification of combustion noise and piston slap noise.
The FFT (Fast Fourier Transform) and CWT (Continuous Wavelet Transform) are combined to further identify the separation results.The calculation results are shown in Figures 14 and 15.The correlation function of cylinder pressure and cylinder head vibration is shown in Figure 16.
From Figure 14, it can be seen that the frequency of IC6 component focuses on 4025 Hz, and the time domain waveform of IC6 component has great change at 390 ∘ CA.From Figure 14(c), the frequency energy is large around the 390 ∘ CA, and it is mainly concentrated around 4000 Hz.According to the prior knowledge of the internal combustion engine, the ignition sequence of the internal combustion engine is 1-5-3-6-2-4.The ignition angle of the number 6 cylinder is around 390 ∘ CA.Due to the fact that the number 1-5 cylinders are covered by lead, the frequency energy of the number 6 cylinder is bigger than the number 1-5 cylinders.Moreover, from Figure 16, cylinder pressure and cylinder head vibration have a good correlation near 4000 Hz.The combustion noise is mainly caused by the change of cylinder pressure.Near 4000 Hz, the main information about combustion noise is included.Thus, it can be determined that the IC6 component is the combustion noise.
From Figure 15, it can be seen that the frequency of IC9 component focuses on 1725 Hz.Combined with Figure 16, the correlation function of cylinder pressure and cylinder head vibration is not good near 1725 Hz.From Figure 15(c corresponding with the ignition angle of the number 6 cylinder.In order to further confirm the noise component of the IC9, the internal combustion engine drag test is carried out.By the drag test, the piston slap vibration signal is measured separately.The FFT of IC9 and piston slap vibration is shown in Figure 17.From Figure 17, the IC9 spectrum is very consistent with the piston slap vibration spectrum.Combined with Figure 18, correlation function of IC9 can be seen and piston slap vibration is basically greater than 0.5, and therefore it can be considered that IC9 and piston slap vibration have a good correlation.Thus, it can be determined that the IC9 component is mainly the piston slap noise.

EWT-RobustICA Method.
In this part, the EWT-Robus-tICA method is used to separate the noise source of internal combustion engine.Firstly, the EWT algorithm is applied to decompose the preprocessed noise signal.The calculation results are shown in Figure 19.
After the preprocessed noise signal is decomposed by the EWT algorithm, empirical modal components can be obtained.But it may have some false empirical modal components.Thus, it is necessary to remove the false empirical modal components.Hence, the correlation coefficient between empirical modal components and preprocessed noise signal is utilised to determine whether the empirical modal component is removed or not.The correlation coefficient between empirical modal components and preprocessed noise signal is shown in Table 4.
From Table 4, it can be seen that the maximum value of the empirical modal components is 0.8532.Thus, if the empirical modal components are less than 0.08532, it should be removed.Therefore, F3, F4, F5, F6, F7, F8, and F9 should be removed.Then the remaining F1 and F2 and the preprocessed noise signal are combined together to form a new signal group.The RobustICA algorithm is utilised to extract the    20.By analysing, it can be preliminarily judged that the IC1 component and the IC2 component are the combustion noise and the piston slap noise.Then the FFT and CWT are utilised to further analyse the IC1 component and the IC2 component.The calculation results are shown in Figures 21  and 22.
From Figure 21, it can be seen that the amplitude of IC1 component has great change around 390 ∘ CA, and the ignition sequence of the number 6 cylinder is at 390 ∘ CA.The frequency range of the IC1 component is concentrated around 4000 Hz-6700 Hz.Combined with Figure 16, the correlation function of cylinder pressure and cylinder head vibration is good at 4000 Hz-5000 Hz, but the correlation function is not good at 5000 Hz-6700 Hz.Thus, it can be considered that the IC1 component is mainly the combustion noise.However, it   vibration is not very good.Thus, it can be determined that the IC2 component has some interference components except for the piston slap noise component.It can be seen in the yellow circle in Figure 22(c).

Comparison Discussion
. When adopting the Gammatone-RobustICA method to separate the combustion noise and piston slap noise, the calculation results are shown in Figures 14 and 15.The frequencies of combustion noise and piston slap noise are concentrated at 4025 Hz and 1725 Hz, respectively.When utilising the EWT-RobustICA method, the separated combustion noise and piston slap noise are shown in Figures 21 and 22. Comparing Figure 14 with Figure 21, it can be seen that there are many other frequency components except for the 4025 Hz components in Figure 21.Comparing Figure 15 with Figure 22, it can be seen that there are many other frequency components except for the 1725 Hz components in Figure 22.Moreover, from Figure 17, the IC9 spectrum is consistent with the piston slap vibration spectrum, while in Figure 23, the IC2 spectrum has some difference with the piston slap vibration spectrum.Therefore, compared with the EWT-RobustICA method, the Gammatone-RobustICA method can get more pure combustion noise and piston slap noise.The Gammatone-RobustICA method has a better separation and identification effect than the EWT-RobustICA method.

Conclusions
(1) In the aspect of the internal combustion engine test, the lead covering method is adopted to wrap the numbers 1-5 cylinders, and only the number 6 cylinder parts are exposed.It can effectively insulate the interference noise from the numbers 1-5 cylinders.Thus, it is conducive to separating and identifying the noise sources of the specified number 6 cylinder.
(2) Inspired by the human auditory system, the Gammatone-RobustICA method is proposed to separate and identify the combustion noise and the piston slap noise.Compared with the EWT-RobustICA method, the proposed Gammatone-RobustICA method has a better separation effect.The separation results obtained by the Gammatone-RobustICA method have very fewer other interference components. (b).
Amplitude frequency response curve of Gammatone function

Figure 2 :
Figure 2: A common blind source separation process.

Figure 4 :Figure 5 :
Figure 4: Time domain waveform of the selected source signals.

Figure 6 :
Figure 6: Lead cover of the internal combustion engine.

Figure 7 :
Figure 7: Vibration and noise measuring system of internal combustion engine.

Figure 9 :Figure 10 :
Figure 9: Specific arrangement of the measuring point.

Figure 12 :
Figure 12: Decomposition results by the Gammatone filter bank.

Figure 18 :
Figure 18: Correlation function of IC9 and piston slap vibration.

Table 1 :
Specific parameters of the WP10-240 type diesel engine.
Test bench contains the WP10-240 type high speed diesel engine, transmission shaft, Germany's Siemens 1PL6 AC motor, console, and other related accessories.The specific parameters of the WP10-240 type high speed diesel engine are shown in Table 4.1.Test Platform.The noise and vibration test of internal combustion engine are carried out in a fully enclosed semianechoic chamber.The size of the semianechoic chamber is length 7.04 m × width 6.79 m × height 5.95 m.The free sound field radius is not less than 2 meters, and the background noise is 18 dB.

Table 2 :
Internal combustion engine test conditions.

Table 3 :
Correlation coefficient between the mode components and the preprocessed noise signal.

Table 4 :
Correlation coefficient between empirical modal components and preprocessed noise signal.