Coal–Rock Cutting Sound Denoising Based on Complete Ensemble Empirical Mode Decomposition with Adaptive Noise and an improved Fruit Fly Optimization Algorithm

: The cutting sound signal of a coal mining shearer is an important signal source for identifying the coal–rock cutting mode and load state. However, the coal–rock cutting sound signal directly collected from the industrial ﬁeld always contains a large amount of background noise, which is not conducive to the subsequent feature extraction and recognition. Therefore, efﬁcient noise elimination for the original signal is required. An intelligent processing method based on an improved complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) denoising algorithm is constructed for the cutting sound signal in this paper. CEEMDAN ﬁrst decomposes the sound to generate a series of intrinsic modal functions (IMFs). Because the denoising threshold of each IMF is usually obtained by an experimental test or an empirical formula in the traditional CEEMDAN method, obtaining an optimal threshold set for each IMF is difﬁcult. The processing effect is often restricted. To overcome this problem, the fruit ﬂy optimization algorithm (FOA) was introduced for CEEMDAN threshold determination. Moreover, in the basic FOA, the scouting bee mutation operation and adaptive dynamic adjustment search strategy are applied to maintain the convergence speed and global search ability. The simulation result shows that the signal waveform processed by the improved CEEMDAN denoising algorithm is smoother than the other four typical eliminate noise signal algorithms. The output signal’s signal-to-noise ratio and mean square error are signiﬁcantly improved. Finally, an industrial application of a shearer in a coal mining working face is performed to demonstrate the practical effect.


Introduction
The shearer is one of the key pieces of equipment in the comprehensive mechanized coal mining face.The cutting sound signal generated by the shearer's colliding drum with coal and rock contains a lot of key information during its operation.This information can represent the shearer's current cutting state and load condition.By excavating the cutting sound signal, the load condition and fault information of the shearer could be recognized.In recent years, shearer cutting state identification based on the sound signal has gained the attention of scholars [1][2][3].However, there are many sources of sound in coal mining faces.The original sound signals collected on-site are mixed signals from multiple sound sources and cannot be used directly.How to remove the noise signal from the complex mixed sound signal to the maximum extent and retain the coal and rock cutting signal is extremely important for subsequent feature extraction and recognition.
Machines 2022, 10, 412 2 of 16 Preliminary research shows that the coal-rock cutting sound signal has the characteristics of being nonperiodic, nonstationary, and discontinuous [4].Traditional denoising algorithms based on fast Fourier transform, wavelet transform, and wavelet packet transform are difficult to achieve effectively, removing noise components when processing such signals [5,6].Empirical mode decomposition (EMD), which does not require selecting the basis function and decomposition layer, can adaptively decompose the signals according to different time scales when processing the signal and can extract the features of nonstationary signal variations [7].Torre et al. proposed the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) based on EMD [8].In CEEMDAN, Gaussian white noise is adaptively added to each order of intrinsic mode function (IMF) component, and then superimposed and canceled.In this way, the reconstruction error caused by adding noise will be eliminated in the decomposition iteration.It not only guarantees the accuracy of signal decomposition but also effectively avoids the modal blending in the traditional EMD decomposition process [9].Therefore, signal processing accuracy and completeness are improved.
The signal denoising method based on CEEMDAN is to decompose the signal into several components, and then shrink and transform the components in each order to eliminate the noise components in the mixed signals.Peng and others denoised the original signal through CEEMDAN and filtering algorithms and obtained certain results [10].A joint noise elimination method based on CEEMDAN and a nonlocal mean algorithm was proposed by Zhang and others.Finally, the feasibility and applicability of this method for signal noise reduction in practical engineering are verified [11].The lifting wavelet transform technique and the CEEMDAN method were combined in reference [12] to improve the computational efficiency and signal reconstruction accuracy and realize the denoising of fiber optic gyroscope signals.A denoising method based on CEEMDAN, least mean square difference criterion, and a least mean square adaptive filter were proposed by Li-y and others.By analyzing the simulation signal and the real underwater acoustic signal, the method was proved to have a better denoising effect [13].
The above studies investigated and applied the CEEMDAN-based denoising algorithm, achieving certain results.However, the determination of each order of IMF denoising thresholds is mostly obtained through comparative experiments or empirical formulas, which requires a significant amount of calculation.It is difficult to obtain the best denoising effect in the face of randomly varying noisy signals [14].In recent years, with the widespread use of artificial intelligence, adaptive threshold selection methods based on intelligent optimization algorithms have been gradually adopted.In [15], the CEEMDAN threshold was modified by a niche genetic algorithm, and an application example on ECG signals was shown to verify the effectiveness of the proposed work.A hybrid wind energy prediction model based on CEEMDAN, SVD, and PSO was presented by Zhang and others.The final prediction results showed that the proposed model can improve the effect of wind speed prediction, reduce the prediction error, and provide strong support for the stable operation of wind farms, and the grid connection of power plants [16].
The fruit fly optimization algorithm (FOA) is a new biological intelligent optimization algorithm proposed by Pan in 2012.FOA, as a meta-heuristic method, simulates the foraging behavior of fruit fly individuals and groups, flies in the direction of the highest food concentration, and conducts appropriate iterations to find the global optimal solution [17].The fruit fly is spatially and visually superior to other species.It can even sniff out food sources up to 40 km away and identify aggregations of other animals through its sensitive vision [18].Compared with the abovementioned optimization algorithms, FOA has the advantages of a simple structure, fewer setup parameters, direct applicability to practical applications, easy implementation, and fast convergence.Since FOA was proposed, it has been widely applied in time series prediction [19], neural network parameter optimization [20], task cooperative planning [21], and other fields.FOA has a strong convergence speed and the ability to deal with optimization problems.However, it is single in diversity and does not have a high probability of mutation.As a result, the search space will be limited.When the global optimal is shifted away from an overly converged swarm, it is difficult for FOA to jump away from the local extremum, and diversity is lost.To improve the performance of the basic FOA and avoid local optimization, certain modified approaches were proposed to strengthen the local and global optimal values.In [22], the back-propagation neural network based on an improved fruit FOA (IFOA-BP) was applied to load forecasting.The simulation results showed that the algorithm can improve the prediction accuracy of wind energy and load.Guo and others introduced two concepts, sensitivity, and pheromone, which improve the optimization strategy and position replacement of fruit flies and modify the global optimization characteristics of the algorithm.The experimental results showed that the model's minimum prediction error accuracy is only 0.55%, indicating a more robust prediction effect [23].Furthermore, an olfactory search and a visual search were adopted in FOA to effectively solve the inverse kinematics problem of redundant manipulators [24].Many improvements focusing on the FOA have been elaborated in recent years, but few researchers have been able to balance both local extremes and iterative rates.
To remove noise components from the original coal-rock cutting sound signal, an adaptive threshold noise elimination method based on improved CEEMDAN is proposed in this paper.To overcome the difficulty of determining the optimal denoising threshold for each order of IMF in the traditional CEEMDAN denoising, FOA iterative convergence has been introduced to perform global optimization for each order of the IMF threshold.At the same time, it is easy to fall into a local optimum in the iterative process to avoid the simple diversity of the FOA population.In an artificial bee colony (ABC), the mutation mode of scouting bees is combined with FOA to divide the fruit flies into the native and mutation groups separately.Furthermore, the adaptive dynamic adjustment search strategy is introduced to flexibly generate the range of variation population and improve the convergence accuracy and maintain the convergence speed.The rest of the paper is organized as follows: In Section 2, the basic principles of CEEMDAN and FOA are first introduced.Section 3 introduces the improvement of FOA and the improved CEEMDAN adaptive denoising method and optimization process.In Section 4, to verify the effectiveness and superiority of the proposed method, Gaussian white noise with a different SNR is added to the simulation signal, and different denoising algorithms are compared and verified using simulation.In Section 5, the improved CEEMDAN denoising algorithm is used to denoise the shearer coal-rock cutting sound signal, and a joint time-frequency domain analysis is performed to verify the practical application effect.Some conclusions and outlooks are summarized in Section 6.

CEEMDAN Denoising Algorithm
CEEMDAN adaptively added several pairs of positive and negative Gaussian white noise according to the signal characteristics, decomposed the first layer of IMF components, and immediately calculated the average to get the first IMF.The above operations continue for the residual signal R(t) until the decomposition ends when the residual signal R(t) is a monotone function.Each stage is selected by calculating a signal-to-noise control factor to add positive and negative Gaussian white noise with different signal-to-noise ratios [25].Each order of IMF must satisfy the following two conditions: in the whole time range, the difference between the number of zeros and extreme points in the function is not more than one; the average value of the upper and lower envelope, composed of local maxima and local minima, is zero at any point.CEEMDAN effectively solves the problem that the EMD decomposition process is prone to modal aliasing and improves the completeness, decomposition efficiency, and reconstruction accuracy of the decomposition process [26].Finally, the original signal x(t) can be expressed as: Machines 2022, 10, 412 CEEMDAN is a linear transformation.The coefficients of the decomposed noisy signal include the noise signal coefficient and the useful signal coefficient.CEEMDAN denoising is similar to the principle of wavelet denoising, and the threshold of each order IMF should be determined [27]: where T i denotes the threshold for the i-th IMF, C is a constant between [0, 1], N is the signal data length, and E i is the estimated energy of the first IMF.
where E i is the variance of the first IMF and α and β are constants determined by a large number of independent experiments or empirical formulas.Unlike wavelet threshold denoising, CEEMDAN uses intermittent threshold denoising.For each IMF with a random neighboring zero interval z i,j = [z i,j , z i,j+1 ], a hard threshold denoising function is defined within the interval: Additionally, the soft threshold denoising function is defined in the interval: where h i (r i,j ) is the maximum value in the interval and h i (r i,j ) is the data point in [z i,j , z i,j+1 ].
The final denoised signal is expressed as: where x g (t) is the denoised signal.
The CEEMDAN denoising process reveals that the selection of the soft and hard thresholding functions and the threshold values of each order are important factors influencing the denoising effect.The hard threshold function denoising signal is discontinuous, resulting in a vibration phenomenon [28].The signal is still a continuous signal with a high degree of signal reduction after the soft threshold function denoising.In practice, soft threshold function denoising is mostly used [29].For each order denoising threshold, if the threshold is chosen too small, it will lead to incomplete denoising; if the threshold is chosen too high, the detailed information of the signal will be easily ignored.Currently, the denoising thresholds for each order of IMF are primarily determined by simulation cross-comparison experiments or empirical formulas, which requires a significant amount of work and makes obtaining the optimal denoising threshold difficult.

Fruit Fly Optimization Algorithm
The fruit FOA is a new biological intelligent optimization algorithm proposed by Pan in 2012.By simulating the foraging behavior of fruit flies as individuals and groups, they fly in the direction of the highest food concentration and conduct appropriate iterations to find the global optimal solution [17].During the foraging process, information about the location of food is constantly exchanged between individuals so that the whole group can quickly follow the shortest path to find food.FOA has the advantages of simple structure, few setting parameters, direct application, easy implementation, and fast convergence.Since FOA was proposed, it has been widely applied in time series prediction, neural network parameter optimization, task cooperative planning, and other fields.The FOA optimization process is as follows: (1) Several key parameters of the fruit FOA are determined: population size P, maximum iteration N, search range LR, flight distance FR, and the initialization of the population locations.
where rand(•) is a random function that follows a uniform distribution.(2) Using the sense of smell, the direction and distance of each fruit fly in the population are determined.
where i = 1, 2, 3, . . .P, (X i , Y i ) is the position of the i-th fruit fly.(3) After the location of each fruit fly in the population is determined, the distance Dist i of each fruit fly to the origin is determined, and the reciprocal of Dist i represents the judgment value of taste concentration at this location.
(4) Smell i is solved by S i and Function.The optimal odor concentration for each fruit fly in the population is determined.
(5) The location of the best flavor concentration causes the rest of the population to fly toward that point.
where (X bestindex , Y bestindex ) is the individual position of a drosophila corresponding to its fitness value bestindex.

Improvement of FOA
The steps of the original FOA in searching for the optimal solution will reveal that it is single in diversity and does not have a high probability of mutation.Hence, the search

Improvement of FOA
The steps of the original FOA in searching for the optimal solution will reveal that it is single in diversity and does not have a high probability of mutation.Hence, the search space will be limited.When the global optimal is shifted away from a too converged swarm, it is difficult for FOA to jump away from the local extremum, and diversity is lost.
The ABC has the advantages of simple calculation, fast convergence speed, and strong global optimization ability [31].ABC's ability to jump out of the locally optimal solution is dependent on the scout bee's mutation role.During the search, if the nectar source has not been updated to a better nectar source after n iterations and reaches the optimal solution not updated the number of times L, the nectar source will be abandoned.Scouting bees will go out and randomly generate a new nectar source.To balance the ability of global and local search, some improvements to the original FOA are pointed out in this paper.Scouting bees' mutation method and FOA combined to increase the diversity of the population.In this way, while retaining the local search ability of the original FOA, the global optimization expansion ability of the population is improved, and the convergence ability of the population is maintained.In this case, the optimal solution is not updated.The number of times L is determined by the object applied by the algorithm.If the value is too large, the algorithm will have little improvement; if the value is too small, the local search ability will be weakened.The value of L needs to be determined according to the specific situation.Based on empirical values, L is set as the product of population size and problem dimension [32].
In addition, considering that the optimal fitness value may appear near the local optimal concentration value S i , it is also considered to further expand the search scope, jump out of the local optimal concentration value S i , and improve the convergence rate.Therefore, the adaptive dynamic adjustment search strategy is introduced, and the population position mutation range S i is set as [S i ± 2].The principle of the mutation is presented in Figure 2. The mutation population has the same number of individuals as the basic population.The flowchart of the IFOA can be shown in Figure 3.

Flow of the Proposed Denoising Method
To obtain the best denoising effect of CEEMDAN, this paper uses IFOA adaptive seeking to obtain the best denoising threshold for each order of IMF and proposes an adaptive denoising algorithm based on the improved CEEMDAN.The main process of the algorithm is as follows: Step 1: The parameters of the improved CEEMDAN are initialized, such as P, LR, N, FR, L, and the number of populations K. Positions of each fruit fly population, as well as the flight direction and distance of individuals in each fruit fly population are initialized.In the improved CEEMDAN denoising algorithm, each order IMF threshold is optimized by a fruit fly population.
Step 2: CEEMDAN decomposition is performed on the coal-rock cutting sound signal containing noise to obtain each order of IMF.The distance of each fruit fly individual to the origin Dist i and the flavor concentration judgment value S i are calculated.Each S i corresponds to a g i , which is a threshold feasible solution, so there are P N combinations of threshold feasible solutions.
A wider range of soft threshold functions are used for denoising in this paper, and each set of feasible solutions is used for threshold denoising of each order of IMF.Using the constructed fitness function f it, the fitness values under the feasible solutions corresponding to each IMF threshold are calculated.The maximum value Smellbest and the corresponding optimal threshold Smell g are recorded, updating each population position and the corresponding position of each fruit fly individual.The fitness evaluation function f it constructed by the signal independence theory is as follows: where r 22 is the autocorrelation of noise.Its rising value means that more of the original signal is attached to the noise, so the reconstructed signal will not recover well.hr 21 is the high order cross-correlation between useful signal and noise.Its decline means that the two signals have become more independent.Therefore, the original signal and noise gradually separate.r 11 is the autocorrelation of useful signals.An ascending value implies that its component is more than a component of noise.Hence, the restructured signal has a good recovery [33,34].r 22 , hr 21 , r 11 are defined as follows: where s 1 is the useful signal, s 2 is the noise component, cov[s i , is the mathematical expectation of s i and ϕ(s i ) = s 2 i + s 3 i .r22 and hr 21 tend to the minimum and r 11 tends to the maximum when the location of each fruit fly group is placed in the best.Then, g is the minimal and f it is the maximal [35].
Step 3: The fitness value under the feasible solution of each IMF threshold in each combination is calculated, and the maximum fitness value is recorded as bestsmell.The global maximum value Smellbest is compared with the current maximum value bestsmell.If the global maximum value is smaller than the current maximum value, bestsmell is used to replace Smellbest, and the current optimal threshold value best g is used to replace Smell g , and the location of each fruit fly population is updated.
Step 4: Judge whether the iteration number n reaches the maximum iteration number N or the error meets the preset accuracy.If not, proceed with Step 2-Step 3.If so, iteration is stopped.The CEEMDAN reconstruction is carried out, and a shearer coal-rock cutting sound signal corresponding to each IMF under optimal threshold denoising is output.The flowchart of the process can be shown in Figure 4. Step 4: Judge whether the iteration number n reaches the maximum iteratio number N or the error meets the preset accuracy.If not, proceed with Step 2-Step 3. so, iteration is stopped.The CEEMDAN reconstruction is carried out, and a sheare coal-rock cutting sound signal corresponding to each IMF under optimal threshold de noising is output.The flowchart of the process can be shown in Figure 4.

Experimental Data and Evaluation Indicators
The denoising algorithm was tested by combining simulation with engineering ex periments to verify its feasibility and effectiveness.The simulation environment of thi paper is the Matlab2020 software platform, built on the Windows 10 system.A piece o the pure train sound signal from the Matlab2020 sptool toolbox is extracted with a sam pling frequency of 8192 Hz and data length of 12,880.White Gaussian noise with inpu SNR (SNRin) of 0 dB, 10 dB, 20 dB, and 30 dB was added to the pure sound signal, re spectively.The original signal and the noisy signal were shown in Figure 5.It can b seen from the waveform of the original signal that the amplitude of the signal varie widely in the time domain, and there were two obvious ups and downs.The first am plitude rose and fell rapidly, while the second rose and fell relatively slowly.In add tion, with the decrease in SNRin, the waveform of the original signal was gradually sub merged.Taking SRNin = 10 dB as an example, the CEEMDAN decomposition of th noise-containing signal was then performed to obtain each order of IMF component and the residual signal as shown in Figure 6.The decomposed signal can be expresse as X = IMF1 + IMF2 + IMF3 + … + IMF14 + RES, which corresponds to Formula (1).It ca be seen that each IMF component had no mode aliasing phenomenon and the noise wa mainly distributed in the high-frequency inherent component.

Simulation 4.1. Experimental Data and Evaluation Indicators
The denoising algorithm was tested by combining simulation with engineering experiments to verify its feasibility and effectiveness.The simulation environment of this paper is the Matlab2020 software platform, built on the Windows 10 system.A piece of the pure train sound signal from the Matlab2020 sptool toolbox is extracted with a sampling frequency of 8192 Hz and data length of 12,880.White Gaussian noise with input SNR (SNR in ) of 0 dB, 10 dB, 20 dB, and 30 dB was added to the pure sound signal, respectively.The original signal and the noisy signal were shown in Figure 5.It can be seen from the waveform of the original signal that the amplitude of the signal varied widely in the time domain, and there were two obvious ups and downs.The first amplitude rose and fell rapidly, while the second rose and fell relatively slowly.In addition, with the decrease in SNR in , the waveform of the original signal was gradually submerged.Taking SRN in = 10 dB as an example, the CEEMDAN decomposition of the noise-containing signal was then performed to obtain each order of IMF components and the residual signal as shown in Figure 6.The decomposed signal can be expressed as X = IMF1 + IMF2 + IMF3 + . . .+ IMF14 + RES, which corresponds to Formula (1).It can be seen that each IMF component had no mode aliasing phenomenon and the noise was mainly distributed in the high-frequency inherent component.The signal-to-noise ratio of the output signal (SNRout) and mean square error (MSE) were introduced as evaluation indexes to compare the denoising effects of the five algorithms intuitively.The larger the SNR and the lower the MSE after denoising, the greater the denoising effect will be.The formulas of SNRout and MSE are as follows: where ( ) x t is the noisy sound signal, ( ) s t is the original pure signal,  ( ) x t is the re- constructed sound signal, and N is the signal length.The signal-to-noise ratio of the output signal (SNRout) and mean square error (MSE) were introduced as evaluation indexes to compare the denoising effects of the five algorithms intuitively.The larger the SNR and the lower the MSE after denoising, the greater the denoising effect will be.The formulas of SNRout and MSE are as follows:

Comparative Analysis
where ( ) x t is the noisy sound signal, ( ) s t is the original pure signal,  ( ) x t is the re- constructed sound signal, and N is the signal length.The signal-to-noise ratio of the output signal (SNR out ) and mean square error (MSE) were introduced as evaluation indexes to compare the denoising effects of the five algorithms intuitively.The larger the SNR and the lower the MSE after denoising, the greater the denoising effect will be.The formulas of SNR out and MSE are as follows:

Comparative Analysis
where x(t) is the noisy sound signal, s(t) is the original pure signal, x(t) is the reconstructed sound signal, and N is the signal length.In the case of different SNRin, SNRout, and MSE, they were shown in Table 1 after signal processing by five denoising algorithms.As can be seen from Table 1, with the decrease in the SNRin, SNRout gradually increases and MSE gradually decreases, indicating that the recovery degree of the denoising algorithm to the original signal was gradually enhanced.The improved CEEMDAN denoising algorithm showed that SNRout and MSE were better than the other four algorithms in the case of four different SNRin.In the case of different SNR in , SNR out , and MSE, they were shown in Table 1 after signal processing by five denoising algorithms.As can be seen from Table 1, with the decrease in the SNR in , SNR out gradually increases and MSE gradually decreases, indicating that the recovery degree of the denoising algorithm to the original signal was gradually enhanced.The improved CEEMDAN denoising algorithm showed that SNR out and MSE were better than the other four algorithms in the case of four different SNR in .Under different SNR in , compared with EMD denoising, CEEMDAN denoising, CEEMDAN-PSO denoising, and CEEMDAN-FOA denoising, the SNR out of the proposed algorithm improved by 8.92%, 5.05%, 3.11%, and 1.29% on average, respectively, and the maximum increase was 9.60% (SNR in = 10 dB), 5.42% (SNR in = 0 dB), 3.67% (SNR in = 20 dB), and 1.36% (SNR in = 0 dB).The MSE and SNR out had similar variation patterns.Compared with the four algorithms, MSE improved by 39.55%, 32.60%, 20.99%, and 12.19% on average, respectively, and the maximum improved by 42.71% (SNR in = 10 dB), 38.82% (SNR in = 0 dB), 37% (SNR in = 0 dB), and 21.68% (SNR in = 10 dB).This indicates that under different noise input conditions, the proposed algorithm can achieve better denoising effects compared with EMD denoising, CEEMDAN denoising, CEEMDAN-PSO denoising, and CEEMDAN-FOA denoising.

Industrial Application
The project team built an experimental platform for coal-rock cutting to denoise the sound signals collected on-site and compare the processing effects to further verify the effect of the algorithm proposed in this paper in practical engineering.The equipment layout and experimental site situation were shown in Figure 8.The sound sensor was installed on the rocker arm of the shearer facing the sample cutting and was fixed in the metal shell.The installed sound sensor belongs to the large diaphragm lateral acoustic-type external bias capacitive sensor.Shearer coal and rock cutting sound signals are similar to vibration, temperature, current, and other signals, all belonging to a one-dimensional signal sequence.The difference is that the characteristic frequencies of sound signals tend to be concentrated in the middle and high-frequency ranges (6.2016-8.9578kHz and 11.7141-17.2266kHz).According to the Nyquist sampling theorem, the sampling frequency of the sound sensor should be at least twice that of 20 kHz to avoid aliasing.To achieve better sampling, the sampling frequency of the sound sensor was set to 44.1 kHz, and the maximum frequency component in the signal was 22.05 kHz, which can realize an omnidirectional direction.During the experiment, the sound sensor was installed on the shearer, and the sound signal collected at the front end was transmitted to the signal processing device through the MIC-IN interface.The shearer cutting operation is not only identified by the sound signal of coal and rock cutting caused by the collision between the drum and rock.It also includes the low-frequency area sound generated by the shearer's walking mechanism, the sound generated by the hydraulic heightening system, the sound generated by the operation of other equipment on-site, and the white noise distributed throughout the frequency and the impact noise with a random frequency distribution.Thus, the collected sound signals contained a lot of background noise, and the SNR of the shearer coal-rock cutting sound signal is low.It has an impact on the subsequent identification of coal and rock distribution states and shearer cutting states.The original signal was a piece of sound signal collected on-site.The improved CEEMDAN denoising algorithm proposed in this paper was used to carry out adaptive denoising of the original signal to improve the SNR of the shearer coal-rock cutting sound signal.To verify the method's effect on the signals, the characteristics of the time domain and frequency domain before and after denoising were analyzed.The original signal and the denoised signal were shown in Figure 9.The FFT was performed on the original signal and the denoised signal, and the characteristics of the frequency domain were shown in Figure 10.In addition, for the particularity of coal-rock cutting sound signal, a further study was made.The characteristics of the middle and high-frequencys domain of the coal-rock cutting sound signal before and after denoising were analyzed and compared.It was shown in Figure 11.As can be seen from Figure 9, the original signal had obvious noise components and many burrs.However, the improved CEEMDAN denoising method significantly reduced the noise component of the reconstructed signal and retained the original coal-rock cutting sound signal to the maximum extent.As can be seen from Figure 10, the frequency components of the signal were distributed in the whole frequency range [0-22.05KHz], and the signal spectrum had strong divergence and extremely dense frequency.It has an influence on feature extraction and cutting state recognition of shearer.After denoising, the frequency of background noise was obviously reduced and some wave peaks appeared.At the same time, the signal amplitude decreased after denoising, which was due to the range mutation caused by the removal of noise components.As can be seen from Figure 11, the original signal frequency had strong divergence, frequency density, and chaotic distribution in the mid-to-high-frequency interval of the main distribution of the shearer coal-rock cutting sound signal.Compared with the original signal, the divergence of the signal after denoising was weakened, the frequency of the background noise was obviously reduced, and a more obvious wave peak appeared.The working state of the shearer can be identified according to the wave peak.It showed that the improved CEEMDAN denoising algorithm retained significant advantages in engineering signal processing.

Conclusions and Future Work
An adaptive denoising method based on improved CEEMDAN for shearer coal and rock cutting sound signals is proposed in this study.The scouting bee mutation method was combined with FOA, and an adaptive dynamic adjustment search strategy was introduced to avoid falling into the local extremum, improve the convergence accuracy, and maintain the convergence speed based on basic FOA.After CEEMDAN decomposition of the coalrock cutting sound signal containing noise, the highest SNR of the reconstructed signal was taken as the objective function, and IFOA was used to find the optimal threshold of each order of IMF adaptively through an iterative process to obtain the best denoising effect.The simulation results show that compared with EMD denoising, CEEMDAN denoising, CEEMDAN-PSO denoising, and CEEMDAN-FOA denoising, the SNR out and MSE of the proposed algorithm were improved in different input noise environments.It showed that the algorithm has obvious advantages in denoising under different input noise environments.In the denoising experiment of the shearer coal-rock cutting sound signal, it was found that the noise component of the signal in the time domain diagram and the noise frequency in the frequency spectrum diagram were significantly reduced.In addition, the frequency domain diagram shows a wave peak that determines the working state.The method's effectiveness in practical applications was proved.
However, this method also has some shortcomings and drawbacks, which can be listed as follows: On the one hand, although the proposed improved CEEMDAN denoising method is more effective than other denoising methods, the computational time is still a serious problem that cannot be ignored.On the other hand, the parameters in the optimization process are determined based on previous scholars and a large number of simulation experiments and lack a rigorous mathematical derivation.Several improvements to the proposed method are planned for the future.These may include an improved algorithm code with higher execution efficiency and an appropriate scheme for determining the optimization parameters.

( 6 ) 16 Figure 1 .
Figure 1.Process of foraging behavior for a fruit fly group.

Figure 1 .
Figure 1.Process of foraging behavior for a fruit fly group.
the location of each fruit fly population is updated.

Figure 4 .
Figure 4. Based on the improved CEEMDAN coal-rock cutting sound signal denoising process.

Figure 4 .
Figure 4. Based on the improved CEEMDAN coal-rock cutting sound signal denoising process.

Figure 6 .
Figure 6.IMF components of each order with residual signal (SNR in = 10 dB).

Figure 8 .
Figure 8.(a) Experimental site equipment arrangement; (b) Experimental process of coal-rock cutting.The sound sensor was installed on the rocker arm of the shearer facing the sample cutting and was fixed in the metal shell.The installed sound sensor belongs to the large diaphragm lateral acoustic-type external bias capacitive sensor.Shearer coal and rock cutting sound signals are similar to vibration, temperature, current, and other signals, all belonging to a one-dimensional signal sequence.The difference is that the characteristic frequencies of sound signals tend to be concentrated in the middle and high-frequency ranges (6.2016-8.9578kHz and 11.7141-17.2266kHz).According to the Nyquist sampling

Figure 8 .
Figure 8.(a) Experimental site equipment arrangement; (b) Experimental process of coal-rock cutting.

Table 1 .
SNR out and MSE of different denoising methods.

Table 1 .
SNRout and MSE of different denoising methods.