A Novel Denoising Method for an Acoustic-Based System through Empirical Mode Decomposition and an Improved Fruit Fly Optimization Algorithm

Generally, the sound signal produced by transmission unit or cutting unit contains abundant information about the working state of a machine. The acoustic-based diagnosis system presents some distinct advantages in some severe conditions particularly due to its unique non-contact measurement and unlimited use at the installation site. However, the original acoustic signal collected from manufacture process is always polluted by various background noises. In order to eliminate noise components from machinery sound effectively, an empirical mode decomposition (EMD) threshold denoising method optimized by an improved fruit fly optimization algorithm (IFOA) is launched in this paper. The acoustic signal was first decomposed by the adaptive EMD to obtain a series of intrinsic mode functions (IMFs). Then, the soft threshold function was applied to shrink the IMF coefficients. While the threshold of each IMF was determined by statistical estimation and empirical value for traditional EMD denoising, the denoising effect was often not desired and time-consuming. To solve these disadvantages, fruit fly optimization algorithm (FOA) was introduced to search global optimal threshold of each IMF. Moreover, to enhance the group diversity during production of the next generation of fruit flies and balance the local and global searching ability, a variation coefficient and a disturbance coefficient was introduced to the basic FOA. Then, a piece of simulated acoustic signal produced by the train was applied to validate the proposed EMD and IFOA threshold denoising (EMD-IFOA). The simulation results, which decreased 35.40% and 18.92% in mean squared error (MSE) and percent root mean square difference (PRD) respectively, and increased 40.36% in signal-to-noise ratio improvement (SNRimp) compared with basic EMD denoising scheme at SNR = 5 dB, illustrated the effectiveness and superiority of the proposed approach. Finally, the proposed EMD-IFOA was conducted on an actual acoustic-based diagnosis system for cutting state recognition of the coal mining shearer to demonstrate the practical effect.


Introduction
Over the last years, customers who have wanted to monitor the real-time running condition of the machine's working process have studied and developed various monitoring methods.An important difference in these methodologies can be summarized as the different measuring system.Generally, the machine monitoring system can be classified as utilizing either contact or non-contact measurement.
Contact approaches are mainly aimed at vibration and strain signals to collect dynamic information, however, carrying out of such methods bears the disadvantages of limited detecting positions, which are strict with respect to the measuring directions, and have difficult-to-maintain detectors in some severe situations [1].Non-contact methods such as acoustic-based, image-based and laser vibrometry-based systems, however, can be used conveniently and do not affect the machine, thereby increasing efficiency and allowing online supervisory control [2][3][4].The acoustic signal has attracted much attention in industrial monitoring due to its strong robustness to illumination intensity, air quality and barrier objects compared with the image [5].In addition, it involves a lower cost and simpler structure than the laser vibrometry measuring system.Unfortunately, one of the most important unfavorable factors for the acoustic-based monitoring system is the substantial noise component in the initial sound signal [6].Thus, the denoising operation is essential and the performance of noise elimination directly influences the effect of subsequent processing.
Typical denoising methods can be divided into three types: time domain analysis, frequency domain analysis and joint time-frequency analysis [7].Many researchers pointed out that methods considering both scales perform better than those focusing on only one.The denoising process can be translated into elimination of noise elements and keeping useful components as much as possible in the two domains.With respect to the development of signal denoising, the most influential time-frequency joint analysis methods were short-time Fourier transform (STFT)-based schemes [8], wavelet transform (WT)-based schemes [9] and empirical mode decomposition (EMD)-based schemes [10].Due to the fixed window function during STFT, there exists an obvious contradiction between the time domain resolution and the frequency domain, which is inappropriate for signaling with complicated components [11].On the basis of STFT, the concepts of multi-scale and multi-resolution were introduced into WT to overcome the problem.Then, the soft and hard threshold denoising methods were proposed by Donoho based on WT and became the fundamental theory of modern denoising approaches [12,13].However, the selection of wavelet basis and decomposition level depended on previous experience and trails, and strict mechanisms were scarce.Since it was proposed by Huang et al. in 1998 [14], EMD has been widely applied for biomedical [15], fault diagnosis [16,17], and audio recognition [18] applications, among others.EMD exhibits obvious advantages in terms of non-linear, non-stationary and intermittent signals compared to STFT and WT due to its adaptive basis, which usually offers a physically meaningful representation of the underlying processes.As a powerful adaptive decomposition tool, EMD decomposes the original signal into a series of intrinsic mode functions (IMFs) and is also well suited to estimating the noise or frequency in measurement domains.The EMD-based denoising system was widely used in Lidar detection [19], the inertial navigation system [20], and acoustic signal [21], among others.The first attempt at applying EMD as a noise elimination tool emerged from the need to determine whether a specific IMF primarily contains useful information or noise [22].The EMD denoising scheme can be regarded as partially reconstructing the signal, retaining only the IMFs that include useful information and discarding those contain primarily noise.To cancel industrial noise from speech signal, EMD was applied to generate the components of signal and the selection of IMFs for signal regeneration is optimized by using objective function of the intelligent particle swarm optimization algorithm (PSO) [23].Then, the EMD denoising procedure inspired by wavelet threshold was proposed and achieved some preliminary outcomes.However, the direct application of wavelet threshold into EMD was inappropriate in principle and could result in incontinuity of the reconstructed signal [24].Immediately afterwards, the concept of interval threshold denoising was introduced to overcome the problem and was widely applied in the noise cancellation process.In [21], Baskar et al. compared Fourier Bessel (FB) expansion-based and EMD-FB-based algorithms to denoise underwater acoustic signals as the sea water usually contained unwanted background noise.In [25], an adaptive thresholding denoising technique by employing ensemble empirical mode decomposition (EEMD) and the intelligent genetic algorithm (GA) was applied for electrocardiogram signals.The experimental results indicated that the proposed denoising method outperformed other denoising methodologies in terms of signal-to-noise ratio, mean square error and percent root mean square difference.
The fruit fly optimization algorithm (FOA) was proposed by Pan in 2012 [26].FOA is a meta-heuristic optimization method simulating the intelligent foraging process of the fruit fly swarm.The natural fruit fly has superior ability in the senses of olfactory and vision system compared with other species.It can even smell food odor from as far as 40 km away and locate other flocks by its sensitive vision.Compared with other existing meta-heuristics, FOA involves an easily understandable principle, a simple structure and fewer parameters.Since the FOA was proposed, it has been widely applied in parameter optimization [27], machine scheduling [28], and pattern recognition [29], etc.The FOA searches the optimal solution only in the neighborhood of previous group location, which indicates that the location of each fruit fly during the iteration process is determined according to previous swarm coordinate and the preset fly distance range.The optimization is therefore simplified furthest to decrease the convergence time and enhance the local searching speed.Excessive pursuit of local searching ability has also reduced the diversity of the fruit fly swarm.Many researchers point out that the basic FOA can easily to fall into local extreme due to poor global searching capacity [30].In [29], a genetic proportion coefficient was introduced into the basic FOA to enhance the global searching ability.A simulation example, with the cutting pattern clustering accuracy of 95%, and some comparisons, verified the effectiveness and superiority of the proposed method and an industrial test proved its practical effect.In [31], an ensemble parameter FOA (EPFOA) was proposed to improve the precision and robustness of complex parts quality inspection.The proposed algorithm was superior to other algorithms in terms of accuracy and robustness by some comprehensive comparisons with other registration methods and some variants of FOA.Shan et al. adopted a new linear generation mechanism of candidate solution to enhance the searching efficiency and improve the searching quality of basic FOA [30].
Motivated by the above-mentioned observations, an adaptive EMD threshold denoising method for acoustic signals is proposed in this paper.The FOA is introduced into the threshold determination process of the denoising approach.In order to balance the local and global searching ability, a variation coefficient and a disturbance coefficient are added into the basic FOA.The rest of this paper is organized as follows.In Section 2, the basic EMD noise elimination method and optimization process of FOA are presented.In Section 3, the improved FOA (IFOA), introducing the variation coefficient and disturbance coefficient and the denoising process based on the EMD and IFOA (EMD-IFOA), is elaborated.In Section 4, the noisy train sound signal with different signal-to-noise ratios (SNRs) is analyzed to verify the superiority and effectiveness of the proposed noise elimination approach, and a piece of sound signal collected from a practical acoustic-based diagnosis system for the long wall coal mining shear is applied to validate the practical effect.Some conclusions and outlooks are summarized in Section 5.

Empirical Mode Decomposition-Based Denoising
EMD-based denoising consists of two main parts: EMD and threshold denoising.The EMD was pioneered by Huang et al. as an automatic tool for completely data-driven and adaptive multi-component signal processing methods in several IMFs so that the sum of them is equal to the original signal.Each IMF is required to satisfy two basic conditions: the number of the extremum points and the number of zero crossings must be equal or must differ by one at most; the mean of the upper and lower envelopes determined by the local maxima and minima should be zero at any point [16].
Based on the principles of IMF, the process of extracting IMFs can be summarized in the following steps [14]: Step 1.1:For an arbitrary signal x(t), all extrema are searched at first.The upper and the lower envelopes are developed by connecting all the maximal and the minimal points through cubic splines respectively.The average value of the two envelopes is recorded as m 1 (t) and the difference between x(t) and m 1 (t) is defined as h 1 (t), which can be expressed as follows: If h 1 (t) satisfies the two requirements of IMF, then it becomes the first IMF of the original signal, named as imf 1 .Else, it is treated as the x(t) and the above operation is repeated.
Step 1.2:The reminder r 1 (t) is obtained by separating imf 1 from x(t): Then r 1 (t) is regarded as the original signal and the above-mentioned steps are executed again until the k-th remainder r k becomes a monotonic function.r k can be presented as follows: Step 1.3: x(t) can be decomposed by the sum of k IMFs and a residual, which can be shown as follows: where res is the residual.Generally, lower-order IMFs have higher frequency components while higher-order IMFs contain lower frequency units.In addition, res represents the average trend of x(t).If x(t) contains some noise, the IMFs and residual will also be polluted by the noise components.In order to eliminate the useless signal, researchers attempt a series of EMD-based denoising methods inspired by wavelet thresholding.Among these methods, the EMD interval thresholding scheme (EMD-IT) is mostly applied, which can be expressed as follows [20]: Step 2.1: The denoising thresholds of each IMF are first calculated.In the EMD-IT study, multiples of the IMF dependent universal thresholds are used as follows: where T i is the threshold of i-th IMF, i = 1, 2, 3, . . ., k, C is a constant, N is the length of signal and E i is the estimated energy of the corresponding IMF, which can be calculated as follows: where E 1 is the variance of the first IMF, i = 2, 3, 4, . . ., k, α and β are parameters that have be estimated by a large number of independent noise realizations and their IMFs.
Step 2.2: Unlike the direct application of the hard or soft threshold function in wavelet thresholding, interval thresholding scheme is used in EMD-IT based on the mode cell filter which can reduce the discontinuity.For the arbitrary two adjacent zero-crossings z i,j = [z i,j , z i,j+1 ] in i-th IMF, the denoising process can be translated into: for hard thresholding and for soft thresholding, where h i (r i,j ) donates the extremum of the corresponding zero-crossing interval and h i (z i,j ) represents the points from z i,j to z i,j+1 in i-th IMF.
Step 2.3: The denoised signal can be presented as follows: where the introduction of parameters M 1 and M 2 gives flexibility on the elimination of the noisy low-order IMFs and on the optional thresholding of the high order ones which contain little noise energy.

Fruit Fly Optimization Algorithm
FOA is an interactive evolutionary optimization method based on the principle of selecting the superior and eliminating the inferior evolved.As the fruit fly has obvious superiority over other species in olfactory and vision system, the fly group can find food quickly.FOA is evolved from the food finding behavior of the fruit fly and applied to the adaptive optimization process.The detailed iterative process of FOA can be concluded as follows [26]: Step 3.1: The maximal iteration number is set as N, the number of fruit flies is P, the initial fruit fly group coordinate is (X fg , Y fg ), the fruit fly location range LR and the random fly distance range is FR.
Step 3.2: The flight direction and distance of each fly individual are achieved randomly according to its osphresis: where i = 1, 2, 3, . . ., P and i-th fly location donates (X i , Y i ).
Step 3.3: The distance of the fly individual to the origin is named as D i and the reciprocal of the distance is named smell concentration judgment value S i , which can be presented as follows: The smell concentration or fitness smell i is calculated subsequently according to fitness function: where f is the fitness function.The maximal smell concentration and the corresponding coordinate of fly are recoded as bestsmell and bestindex respectively, where [bestsmell, bestindex] = max(smell).
Step 3.4: If the present bestsmell is bigger than the former global smellbest, then smellbest is replaced by the bestsmell and the location of fruit fly group is also changed: Step 3.5: If the smellbest reaches the preset conditions or the iteration number comes to the maximal iteration N, then the circulation stops.Otherwise, repeat steps 3.2 to 3.4.

Improvement on FOA
It can be seen that basic FOA search optimal solution is in the neighborhood of previous group location only.In other words, the location of each fruit fly during the iteration process is determined according to previous swarm coordinate and the preset fly distance range.The optimization is therefore simplified further to shorten the program code and decrease the convergence time.Unfortunately, an excessive pursuit of local searching ability has also reduced the diversity of the fruit fly swarm.The basic FOA can easily to fall into a local extreme due to poor global searching capacity.In order to balance the optimization process, a variation coefficient C v and a disturbance coefficient C d are introduced during the process of producing new fruit flies to enhance the group diversity.The principle of the two coefficients is presented in Figure 1.The smell concentration of each individual is calculated to search the maximal fitness in Step 3.3.Then, the concentrations are ranked as (smell1, smell2, smell3, …, smelli, …, smellP), where smell1 is the biggest and smellP is the worst.For any smelli > Cv, the location of the corresponding fruit fly is updated according to Equation (11).That is to say, the fly moves towards the best one.If smelli < Cv, the variation operation is conducted and the fly obtains its location in a new way, which can be expressed as follows: where LR' = max(LR, Xfg, Yfg), the reason lies in the coordinate range of the fruit fly swarm may exceed the original scope with the optimization process.The location of the variation fly is updated in a bigger range randomly to search appropriate solution.Moreover, the disturbance coefficient Cd is introduced when the position of the group remained unchanged for several iterations.Assuming that the smellbest after m-th iteration could be presented as smellbestm and smellbestn for n-th, if smellbestm = smellbestn and n − m> Cd, then the location of each fruit fly can be obtained as follows: The smell concentration of each individual is calculated to search the maximal fitness in Step 3.3.Then, the concentrations are ranked as (smell 1 , smell 2 , smell 3 , . . ., smell i , . . ., smell P ), where smell 1 is the biggest and smell P is the worst.For any smell i > C v , the location of the corresponding fruit fly is updated according to Equation (11).That is to say, the fly moves towards the best one.If smell i < C v , the variation operation is conducted and the fly obtains its location in a new way, which can be expressed as follows: where LR = max(LR, X fg , Y fg ), the reason lies in the coordinate range of the fruit fly swarm may exceed the original scope with the optimization process.The location of the variation fly is updated in a bigger range randomly to search appropriate solution.
Moreover, the disturbance coefficient C d is introduced when the position of the group remained unchanged for several iterations.Assuming that the smellbest after m-th iteration could be presented as smellbest m and smellbest n for n-th, if smellbest m = smellbest n and n − m > C d , then the location of each fruit fly can be obtained as follows: where i = 1, 2, 3, . . ., P. The aim of the disturbance operation is to expand the searching area when the preset fly distance range does not work in finding a better smell concentration for the fly swarm.
The flowchart of the proposed IFOA is shown in Figure 2.

Flow of the Proposed Denoising Method
In order to eliminate noise component from the acoustic signal adaptively, an EMD-based denoising solution optimized by an improved FOA is proposed in this paper.As some key denoising parameters in Equations ( 5) and ( 6), such as C, α and β, are estimated according to statistic and experiment analysis, the determination of denoising threshold is always tough work and consumes much time.Moreover, an appropriate threshold gives significant efficiency of the EMD-based noise elimination process.In order to solve the above problems, the denoising

Flow of the Proposed Denoising Method
In order to eliminate noise component from the acoustic signal adaptively, an EMD-based denoising solution optimized by an improved FOA is proposed in this paper.As some key denoising parameters in Equations ( 5) and ( 6), such as C, α and β, are estimated according to statistic and experiment analysis, the determination of denoising threshold is always tough work and consumes much time.Moreover, an appropriate threshold gives significant efficiency of the EMD-based noise elimination process.In order to solve the above problems, the denoising thresholds in Equation ( 7) or (8) are determined automatically through the improved FOA instead of Equations ( 5) and (6).The flow of the proposed EMD-IFOA can be summarized as follows: Step 4.1: Parameter initialization.The parameters in EMD-IFOA such as the lower-order coefficient M 1 and high-order coefficient M 2 are given firstly, so the number of fruit fly group is The maximal iteration number is N, the number of fruit flies in each swarm is P, the terminal fitness η, the fruit fly location range LR , the random fly distance range is FR, the variation coefficient C v and disturbance coefficient C d .The initial locations of the fruit fly groups are obtained by Equation (10).
Step 4.3: Threshold optimization.The location and smell concentration judgement of each fly individual can be obtained by the improved FOA subsequently.It should be pointed out that any M 2 − M 1 + 1 smell concentration from different fruit fly swarm can be treated as a potential optimal combination of threshold set.For example, S ij donates i-th smell concentration judgement in j-th swarm, where i = 1, 2, 3, . . ., P, j = 1, 2, 3, . . ., M 2 − M 1 + 1.Then, the useful and noise signal are separated by hard threshold function in Equation (7) or soft function in Equation (8).The estimated signal is reconstructed as the final series according to Equation (9).In order to evaluate the denoising performance of different thresholds, a fitness function or smell concentration function in IFOA based on the independent measurement principle is established as follows: where r 11 and r 22 are second-order correlation coefficient of useful signal and noise signal, respectively, and hr 21 is high-order correlation coefficient between the two series.An increasing value of r 11 implies that the useful component is more than the component of noise.If r 22 and hr 21 are decreasing, then it means that the two signals become more independent to each other.Thus, the original signal and noise are separated gradually.In general, the noise elimination process performances are better with larger r 11 and smaller r 22 and hr 21 .
Step 4.4: Iteration termination.If the iteration number reaches N or the fitness reaches the preset value η, the iteration process is terminated; otherwise, the optimization is continued.The flow of the proposed denoising method based on EMD-IFOA can be presented as Figure 3.
decreasing, then it means that the two signals become more independent to each other.Thus, the original signal and noise are separated gradually.In general, the noise elimination process performances are better with larger r11 and smaller r22 and hr21.
Step 4.4: Iteration termination.If the iteration number reaches N or the fitness reaches the preset value η, the iteration process is terminated; otherwise, the optimization is continued.The flow of the proposed denoising method based on EMD-IFOA can be presented as Figure 3.

Simulation and Analysis
In order to validate the effectiveness and superiority of the proposed method, a piece of simulated train acoustic signal taken from Matlab 2015b Signal Processing Toolbox (8.6, MathWorks Inc., Natick, MD, USA) was extracted as the original pure signal.Then, Gaussian white noise with different signal-to-noise ratio (SNR) was added into the pure signal.Three evaluation parameters were introduced to make a comprehensive comparison amongst wavelet threshold denoising, traditional EMD denoising, EMD-PSO denoising, EMD-FOA denoising and the proposed EMD-IFOA.

Signal Preprocessing
To test the performance of the denoising methods, a piece of pure sound signal series was necessary.As machinery sound collected directly from industrial field was always polluted by the environment background noise, the most convenient method to obtain a clean signal was by synthetizing using simulation.In this paper, a piece of simulated train acoustic signal with sampling frequency of 8192 Hz and length of 12,880 taken from Matlab 2015b Signal Processing Toolbox was extracted as the original pure signal y.Then Gaussion white noise n was added into the original signal.The SNR was introduced to measure the content of noise in the noisy signal x, which could be presented as follows: where P y and P n were effective power of pure signal and noise respectively.The original train sound signal and the noise-added signal with SNR = 5 dB were presented in Figure 4. where Py and Pn were effective power of pure signal and noise respectively.The original train sound signal and the noise-added signal with SNR = 5 dB were presented in Figure 4.

Signal Decomposition and Denoising
Five different denoising solutions with SNR = 5 dB were conducted as follows: (1) Wavelet threshold denoising.The noisy signal was first decomposed by db2 wavelet at five levels [32], and the decomposition result was presented in Figure 5.Then, the soft threshold function was adopted to shrink the wavelet coefficients.Inverse wavelet transform was conducted consequently to reconstruct the signal.

Signal Decomposition and Denoising
Five different denoising solutions with SNR = 5 dB were conducted as follows: (1) Wavelet threshold denoising.The noisy signal was first decomposed by db2 wavelet at five levels [32], and the decomposition result was presented in Figure 5.Then, the soft threshold function was adopted to shrink the wavelet coefficients.Inverse wavelet transform was conducted consequently to reconstruct the signal.(2) EMD denoising.The noisy signal was first decomposed by EMD adaptively, which could be shown in Figure 6.The noisy signal could be decomposed as Then, the shrinkage threshold of each IMF was determined according to Equations ( 5) and ( 6).Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters in Equations ( 5), ( 6) and ( 9) were set as follows [10]: C = 0.7, α = 0.719, β = 2.01, M1 = 2 and M2 = 4. (2) EMD denoising.The noisy signal was first decomposed by EMD adaptively, which could be shown in Figure 6.The noisy signal could be decomposed as x = imf 1 + imf 2 + imf 3 + . . .+ imf 12 + res.
Then, the shrinkage threshold of each IMF was determined according to Equations ( 5) and ( 6).Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters in Equations ( 5), ( 6) and ( 9) were set as follows [10]: C = 0.7, α = 0.719, β = 2.01, M1 = 2 and M2 = 4. (3) EMD-PSO denoising.The noisy signal was first decomposed by EMD adaptively.Then, the threshold of each IMF was determined according to PSO intelligently.Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters during the EMD-PSO were set as recommended in [23]: M1 = 2, M2 = 4, the dimensionality of (3) EMD-PSO denoising.The noisy signal was first decomposed by EMD adaptively.Then, the threshold of each IMF was determined according to PSO intelligently.Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters during the EMD-PSO were set as recommended in [23]: M 1 = 2, M 2 = 4, the dimensionality of each particle was M 2 − M 1 + 1 = 3, the learning factor c 1 and c 2 were set as 2, the iteration number was 100 and the particle number of the population was 10. (4) EMD-FOA denoising.The noisy signal was first decomposed by EMD adaptively.Then, the threshold of each IMF was determined according to FOA automatically.Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters during the EMD-FOA were set as follows: M 1 = 2, M 2 = 4, the number of fruit fly population was M 2 − M 1 + 1 = 3, the iteration number was 100, the fly number in each population was 10, the fruit fly location range and the random fly distance range was 1. (5) EMD-IFOA denoising.The noisy signal was first decomposed by EMD adaptively.Then, the threshold of each IMF was determined according to IFOA automatically.Soft threshold function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters during the EMD-IFOA were set as follows: M 1 = 2, M 2 = 4, the number of fruit fly population was M 2 − M 1 + 1 = 3, the iteration number was 100, the fly number in each population was 10, the fruit fly location range and the random fly distance range was 1, the variation coefficient C v = 0.2 and disturbance coefficient C d = 5.The denoised signal of the five noise elimination schemes are presented in Figure 7.
function was applied and the denoised signal was reconstructed as Equation ( 9).Some parameters during the EMD-IFOA were set as follows: M1 = 2, M2 = 4, the number of fruit fly population was M2 − M1 + 1 = 3, the iteration number was 100, the fly number in each population was 10, the fruit fly location range and the random fly distance range was 1, the variation coefficient Cv = 0.2 and disturbance coefficient Cd = 5.The denoised signal of the five noise elimination schemes are presented in Figure 7.In order to evaluate the noise reduction performance of the five methods, three appraisal parameters such as the mean squared error (MSE), SNR improvement (SNRimp) and percent root mean square difference (PRD) were adopted and defined as follows: In order to evaluate the noise reduction performance of the five methods, three appraisal parameters such as the mean squared error (MSE), SNR improvement (SNR imp ) and percent root mean square difference (PRD) were adopted and defined as follows: 20) where y was the original signal, x i was the denoised signal and x was the noise-added signal.The results were shown in Table 1.It can be seen from the Figure 7 and Table 1 that adaptive noise cancellation methods based on intelligent searching technology performed better than wavelet and traditional EMD threshold denoising.The selection of wavelet basis and decomposition level might be inappropriate.
The recommended threshold for wavelet and EMD denoising shrank the decomposed coefficients excessively, which resulted in an amplitude of denoised signal that was much smaller compared with the original signal.The EMD-FOA exhibited slightly better than the EMD-PSO.Moreover, the EMD-IFOA had obvious preponderance of the MSE, signal-to-noise ratio improvement (SNR imr ) and PRD compared with the previous four schemes, which indicated that basic FOA and PSO fell into the local extreme during searching optimal threshold and the population diversity of the EMD-IFOA was enhanced to produce more threshold candidates.The MSE of EMD-IFOA decreased 5.77% and 4.99% compared with EMD-PSO and EMD-FOA.The PRD reduced 2.34% and 2.03% respectively.The SNR imr of the proposed improved 8.43% and 7.41% compared to the other two intelligent denoising approaches.
In order to research the denoising performance at different SNRs, a further comparison was made and the average value of ten simulation results was discussed.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB, 20 dB noisy conditions in terms of MSE is shown in Figure 8.It can be seen from the figure that the MSE did not have an obvious difference when the SNR was greater than 10 dB for wavelet and EMD threshold denoising.The other three denoising approaches based on intelligent optimization solutions outperformed both in higher and lower noise component.As expected, the proposed EMD-IFOA showed the best performance in all noisy conditions and decreased 8.24% and 6.73% compared with EMD-PSO and EMD-FOA even in the worst case scenario, which proved the proposed method had stronger global searching ability.Another comparison illustrated the performance in term of SNR imr at different SNR value and was presented in Figure 9.It can be seen from the figure that the SNR imr decreased almost linearly with the SNR.The proposed method obtained the biggest value in all SNRs than other approaches.Moreover, the PRD of the five methods was presented in Figure 10 and exhibited similar changing tendency with MSE.The PRD of the EMD-IFOA deceased 6.77%, 5.17%, 4.18%, and 3.39% respectively, compared with wavelet, EMD, EMD-PSO and EMD-FOA denoising schemes.In summary, the proposed method exhibited obvious superiority in the denoising effect, which revealed its terrific searching ability.It can be seen from the Figure 7 and Table 1 that adaptive noise cancellation methods based on intelligent searching technology performed better than wavelet and traditional EMD threshold denoising.The selection of wavelet basis and decomposition level might be inappropriate.The recommended threshold for wavelet and EMD denoising shrank the decomposed coefficients excessively, which resulted in an amplitude of denoised signal that was much smaller compared with the original signal.The EMD-FOA exhibited slightly better than the EMD-PSO.Moreover, the EMD-IFOA had obvious preponderance of the MSE, signal-to-noise ratio improvement (SNRimr) and PRD compared with the previous four schemes, which indicated that basic FOA and PSO fell into the local extreme during searching optimal threshold and the population diversity of the EMD-IFOA was enhanced to produce more threshold candidates.The MSE of EMD-IFOA decreased 5.77% and 4.99% compared with EMD-PSO and EMD-FOA.The PRD reduced 2.34% and 2.03% respectively.The SNRimr of the proposed method improved 8.43% and 7.41% compared to the other two intelligent denoising approaches.
In order to research the denoising performance at different SNRs, a further comparison was made and the average value of ten simulation results was discussed.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB, 20 dB noisy conditions in terms of MSE is shown in Figure 8.It can be seen from the figure that the MSE did not have an obvious difference when the SNR was greater than 10 dB for wavelet and EMD threshold denoising.The other three denoising approaches based on intelligent optimization solutions outperformed both in higher and lower noise component.As expected, the proposed EMD-IFOA showed the best performance in all noisy conditions and decreased 8.24% and 6.73% compared with EMD-PSO and EMD-FOA even in the worst case scenario, which proved the proposed method had stronger global searching ability.Another comparison illustrated the performance in term of SNRimr at different SNR value and was presented in Figure 9.It can be seen from the figure that the SNRimr decreased almost linearly with the SNR.The proposed method obtained the biggest value in all SNRs than other approaches.Moreover, the PRD of the five methods was presented in Figure 10 and exhibited similar changing tendency with MSE.The PRD of the EMD-IFOA deceased 6.77%, 5.17%, 4.18%, and 3.39% respectively, compared with wavelet, EMD, EMD-PSO and EMD-FOA denoising schemes.In summary, the proposed method exhibited obvious superiority in the denoising effect, which revealed its terrific searching ability.

Application
In order to test practical performance of the proposed automatic noise cancellation method for the machinery acoustic signal based on EMD-IFOA, an acoustic-based diagnosis system was constructed in the fully-mechanized coal mining working face.The shearer is a kind of essential equipment in modern coal mining, and the operating condition diagnosis for the shearer is of great importance.Conventional diagnosis systems are mainly based on the vibration signal [33] and cutting stress signal [34].Working lives of the vibration sensors and strain gauges are very short due to the harsh working conditions and contact measurement.The coal output is always restricted due to the frequent maintenance.In January 2016, a working state diagnosis system for coal mining shearers was set up in National Coal Mining Equipment Research and Experiment Center at the China Coal Zhangjiakou Coal Mining Machinery Co., Ltd.(Zhangjiakou, China) to identify whether the shearer was cutting coal seam, rock seam or coal seam-gripping rock.The arrangement of the shearer and the cutting sound collection system were presented in Figure 11.An industrial microphone was installed to collect the sound signal.In order to assure the normal use of the collector under the condition of excessive dust and gas, dust-proof and explosion-proof processing satisfying the coal mining requirements was carried out for the microphone.
As there existed a large number of noise signals among the cutting signal, eliminating noise from the field collected sound signal was a crucial step for improving the recognition accuracy.In order to validate the effectiveness and superiority of the proposed denoising method, a piece of acoustic signal with the sampling frequency of 44.1 kHz and the time duration of 1 s was collected and analyzed.As the frequency spectrum was often applied to further state diagnosis, the frequency distribution of the original field signal and denoised signal were also calculated and presented in Figure 12.It could be seen from the figure that the original signal saltated sharply in amplitude due to the random noise element and the denoised signal turned out to be smooth.Moreover, the frequency component of the original sound was presented disorderly in the whole domain, which resulted in difficulty determining the cutting state of the shearer.While it became regular for the

Application
In order to test practical performance of the proposed automatic noise cancellation method for the machinery acoustic signal based on EMD-IFOA, an acoustic-based diagnosis system was constructed in the fully-mechanized coal mining working face.The shearer is a kind of essential equipment in modern coal mining, and the operating condition diagnosis for the shearer is of great importance.Conventional diagnosis systems are mainly based on the vibration signal [33] and cutting stress signal [34].Working lives of the vibration sensors and strain gauges are very short due to the harsh working conditions and contact measurement.The coal output is always restricted due to the frequent maintenance.In January 2016, a working state diagnosis system for coal mining shearers was set up in National Coal Mining Equipment Research and Experiment Center at the China Coal Zhangjiakou Coal Mining Machinery Co., Ltd.(Zhangjiakou, China) to identify whether the shearer was cutting coal seam, rock seam or coal seam-gripping rock.The arrangement of the shearer and the cutting sound collection system were presented in Figure 11.An industrial microphone was installed to collect the sound signal.In order to assure the normal use of the collector under the condition of excessive dust and gas, dust-proof and explosion-proof processing satisfying the coal mining requirements was carried out for the microphone.
As there existed a large number of noise signals among the cutting signal, eliminating noise from the field collected sound signal was a crucial step for improving the recognition accuracy.In order to validate the effectiveness and superiority of the proposed denoising method, a piece of acoustic signal with the sampling frequency of 44.1 kHz and the time duration of 1 s was collected and analyzed.As the frequency spectrum was often applied to further state diagnosis, the frequency distribution of the original field signal and denoised signal were also calculated and presented in Figure 12.It could be seen from the figure that the original signal saltated sharply in amplitude due to the random noise element and the denoised signal turned out to be smooth.Moreover, the frequency component of the original sound was presented disorderly in the whole domain, which resulted in difficulty determining the cutting state of the shearer.While it became regular for the processed signal, some wave peaks which reflected the working state of the shearer appeared in the spectrogram, and other areas were relatively stable.Furthermore, a zoom on the range of [0, 2000 Hz] was shown in Figure 13.It can be seen from the figure that the amplitude of the denoised signal fluctuates at 50 Hz, 260 Hz and 480 Hz.The working state of the shearer can be diagnosed according to the corresponding frequency of the fluctuation and the value of the spectrum.The denoised signal was more suitable for further processing.
processed signal, some wave peaks which reflected the working state of the shearer appeared in the spectrogram, and other areas were relatively stable.Furthermore, a zoom on the range of [0, 2000 Hz] was shown in Figure 13.It can be seen from the figure that the amplitude of the denoised signal fluctuates at 50 Hz, 260 Hz and 480 Hz.The working state of the shearer can be diagnosed according to the corresponding frequency of the fluctuation and the value of the spectrum.The denoised signal was more suitable for further processing.

Conclusions
In order to eliminate noise elements from the sound signal for the acoustic-based diagnosis system, an adaptive threshold denoising method-based EMD and an improved FOA was designed in this paper.Improved strategies on the basis of introducing a variation coefficient and a disturbance coefficient into the basic FOA were applied in the noise cancellation process.The modified FOA is introduced into the threshold determination process of the traditional EMD denoising approach to replace the threshold selection based on statistics estimation and numerous trials.The proposed method can represent an important pretreatment on field sound signals for an acoustic-based diagnosis system.To demonstrate the effectiveness and superiority of the launched WTD-IFOA, a series of simulations were conducted and some comparisons were presented.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB, processed signal, some wave peaks which reflected the working state of the shearer appeared in the spectrogram, and other areas were relatively stable.Furthermore, a zoom on the range of [0, 2000 Hz] was shown in Figure 13.It can be seen from the figure that the amplitude of the denoised signal fluctuates at 50 Hz, 260 Hz and 480 Hz.The working state of the shearer can be diagnosed according to the corresponding frequency of the fluctuation and the value of the spectrum.The denoised signal was more suitable for further processing.

Conclusions
In order to eliminate noise elements from the sound signal for the acoustic-based diagnosis system, an adaptive threshold denoising method-based EMD and an improved FOA was designed in this paper.Improved strategies on the basis of introducing a variation coefficient and a disturbance coefficient into the basic FOA were applied in the noise cancellation process.The modified FOA is introduced into the threshold determination process of the traditional EMD denoising approach to replace the threshold selection based on statistics estimation and numerous trials.The proposed method can represent an important pretreatment on field sound signals for an acoustic-based diagnosis system.To demonstrate the effectiveness and superiority of the launched WTD-IFOA, a series of simulations were conducted and some comparisons were presented.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB, processed signal, some wave peaks which reflected the working state of the shearer appeared in the spectrogram, and other areas were relatively stable.Furthermore, a zoom on the range of [0, 2000 Hz] was shown in Figure 13.It can be seen from the figure that the amplitude of the denoised signal fluctuates at 50 Hz, 260 Hz and 480 Hz.The working state of the shearer can be diagnosed according to the corresponding frequency of the fluctuation and the value of the spectrum.The denoised signal was more suitable for further processing.

Conclusions
In order to eliminate noise elements from the sound signal for the acoustic-based diagnosis system, an adaptive threshold denoising method-based EMD and an improved FOA was designed in this paper.Improved strategies on the basis of introducing a variation coefficient and a disturbance coefficient into the basic FOA were applied in the noise cancellation process.The modified FOA is introduced into the threshold determination process of the traditional EMD denoising approach to replace the threshold selection based on statistics estimation and numerous trials.The proposed method can represent an important pretreatment on field sound signals for an acoustic-based diagnosis system.To demonstrate the effectiveness and superiority of the launched WTD-IFOA, a series of simulations were conducted and some comparisons were presented.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB,

Conclusions
In order to eliminate noise elements from the sound signal for the acoustic-based diagnosis system, an adaptive threshold denoising method-based EMD and an improved FOA was designed in this paper.Improved strategies on the basis of introducing a variation coefficient and a disturbance coefficient into the basic FOA were applied in the noise cancellation process.The modified FOA is introduced into the threshold determination process of the traditional EMD denoising approach to replace the threshold selection based on statistics estimation and numerous trials.The proposed method can represent an important pretreatment on field sound signals for an acoustic-based diagnosis system.To demonstrate the effectiveness and superiority of the launched WTD-IFOA, a series of simulations were conducted and some comparisons were presented.The performance of the proposed and other four denoising methods in 0 dB, 5 dB, 10 dB, 15 dB, and 20 dB noisy conditions were presented.Overall, the EMD-IFOA had the best results in all noisy conditions.The MSE and PRD of the denoised signal decreased 10.01% and 5.17% respectively even in the worst conditions, and the SNR imp increased from 0.0151 dB to 0.4767 dB compared with traditional EMD solution.The simulations and comparison results proved that the automatic denoising method could effectively eliminate noise and outperform the others.Finally, an industrial application for the shearer in a fully-mechanized coal mining face was performed to validate field effect.In future studies, the authors plan to investigate some improvements in determining the variation coefficient and disturbance coefficient automatically with the proposed approach.
swarm coordinate and the preset fly distance range.The optimization is therefore simplified further to shorten the program code and decrease the convergence time.Unfortunately, an excessive pursuit of local searching ability has also reduced the diversity of the fruit fly swarm.The basic FOA can easily to fall into a local extreme due to poor global searching capacity.In order to balance the optimization process, a variation coefficient Cv and a disturbance coefficient Cd are introduced during the process of producing new fruit flies to enhance the group diversity.The principle of the two coefficients is presented in Figure1.

Figure 1 .
Figure 1.(a) Principle of the variation operation; (b) principle of the disturbance operation.

Figure 1 .
Figure 1.(a) Principle of the variation operation; (b) principle of the disturbance operation.

Figure 4 .
Figure 4.The original signal and white noise-added signal.SNR: signal-to-noise ratio.

Figure 4 .
Figure 4.The original signal and white noise-added signal.SNR: signal-to-noise ratio.

Figure 5 .
Figure 5. Wavelet decomposition result of the noisy signal.

Figure 5 .
Figure 5. Wavelet decomposition result of the noisy signal.

Figure 6 .
Figure 6.Empirical mode decomposition (EMD) result of the noisy signal.IMF: intrinsic mode function.

Figure 6 .
Figure 6.Empirical mode decomposition (EMD) result of the noisy signal.IMF: intrinsic mode function.

Figure 7 .
Figure 7.The denoised signal of the five noise elimination schemes.PSO: particle swarm optimization algorithm; WT: wavelet transform; FOA: fruit fly optimization algorithm.

Figure 7 .
Figure 7.The denoised signal of the five noise elimination schemes.PSO: particle swarm optimization algorithm; WT: wavelet transform; FOA: fruit fly optimization algorithm.

Figure 8 .
Figure 8.The mean squared error of the five methods in different SNRs.MSE: mean squared error.Figure 8.The mean squared error of the five methods in different SNRs.MSE: mean squared error.

Figure 8 .
Figure 8.The mean squared error of the five methods in different SNRs.MSE: mean squared error.Figure 8.The mean squared error of the five methods in different SNRs.MSE: mean squared error.

Figure 9 .
Figure 9.The SNR improvement of the five methods in different SNRs.SNRimr: signal-to-noise ratio improvement.

Figure 9 .
Figure 9.The SNR improvement of the five methods in different SNRs.SNR imr : signal-to-noise ratio improvement.

Figure 9 .
Figure 9.The SNR improvement of the five methods in different SNRs.SNRimr: signal-to-noise ratio improvement.

Figure 10 .
Figure 10.The percent root mean square difference of the five methods in different SNRs.PRD: percent root mean square difference.

Figure 10 .
Figure 10.The percent root mean square difference of the five methods in different SNRs.PRD: percent root mean square difference.

Figure 11 .
Figure 11.(a) The arrangement of the shearer; (b) the arrangement of the cutting sound collection system.

Figure 12 .
Figure 12.(a) The amplitude of the original and the denoised signal; (b) the frequency spectrum of the original and the denoised signal.

Figure 13 .
Figure 13.The zoom on the range of [0, 2000 Hz] of the frequency spectrum.

Figure 11 .
Figure 11.(a) The arrangement of the shearer; (b) the arrangement of the cutting sound collection system.

Figure 11 .
Figure 11.(a) The arrangement of the shearer; (b) the arrangement of the cutting sound collection system.

Figure 12 .
Figure 12.(a) The amplitude of the original and the denoised signal; (b) the frequency spectrum of the original and the denoised signal.

Figure 13 .
Figure 13.The zoom on the range of [0, 2000 Hz] of the frequency spectrum.

Figure 12 .
Figure 12.(a) The amplitude of the original and the denoised signal; (b) the frequency spectrum of the original and the denoised signal.

Figure 11 .
Figure 11.(a) The arrangement of the shearer; (b) the arrangement of the cutting sound collection system.

Figure 12 .
Figure 12.(a) The amplitude of the original and the denoised signal; (b) the frequency spectrum of the original and the denoised signal.

Figure 13 .
Figure 13.The zoom on the range of [0, 2000 Hz] of the frequency spectrum.

Figure 13 .
Figure 13.The zoom on the range of [0, 2000 Hz] of the frequency spectrum.

Table 1 .
The denoising performance of the five methods (with SNR = 5 dB).SNR: signal-to-noise ratio; MSE: mean squared error; PRD: percent root mean square difference; SNR imp : signal-to-noise ratio improvement; EMD: empirical mode decomposition; PSO: particle swarm optimization algorithm; FOA: fruit fly optimization algorithm; IFOA: improved fruit fly optimization algorithm.