Coal and Rock Hardness Identification Based on EEMD and Multi-Scale Permutation Entropy

This study offers an efficient hardness identification approach to address the problem of poor real-time performance and accuracy in coal and rock hardness detection. To begin, Ensemble Empirical Mode Decomposition (EEMD) was performed on the current signal of the cutting motor to obtain a number of Intrinsic Mode Functions (IMFs). Further, the target signal was selected among the IMFs to reconstruct the current signal according to the energy density and correlation coefficient criteria. After that, the Multi-scale Permutation Entropy (MPE) of the reconstructed signal was trained by the Adaboost improved Back Propagation (BP) neural network, in order to establish the hardness recognition model. Finally, the cutting arm’s swing speed and the cutting head’s rotation speed were adjusted based on the coal and rock hardness. The simulation results indicated that using the energy density and correlation criterion to reconstruct the signal can successfully filter out noise interference. Compared to the BP model, the relative root-mean-square error of the Adaboost-BP model decreased by 0.0633, and the prediction results were more accurate. Additionally, the speed control strategy based on coal and rock hardness can ensure the efficient cutting of the roadheader.


Introduction
In coal mining operations, the dentification of coal rock hardness is a prerequisite for high-efficiency cutting [1,2]. Due to the complex and variable physical properties of coal and rock, density and hardness are constantly changing and exceedingly unpredictable. Because of the fluctuating loads on the cutting head, the cutting motor's output power is unstable, and it often performs in an overload or underload state. Therefore, accurate coal hardness detection can assure high coal mining efficiency [3,4]. In recent years, recognition methods have mostly focused on coal-rock image analysis [5,6], multi-sensor information fusion [7,8], acoustic signal analysis [9,10], etc. The idea of the image recognition of coal hardness and type is to compare image features of the working face with those in the database [11]. The identification method selects various features, such as gray scale, texture, shape information, and their combination [12], to represent coal and rock traits. Research has been conducted to determine the type of coal and rock by comparing the gray distribution, average gray value [13], and texture characteristics [14] of the images. However, this method does not always extract the rich feature information in the coal image properly, resulting in an unsatisfying final output. To tackle this problem, Si et al. [15] investigated a new method of coal rock recognition based on the deep convolutional neural network (CNN). Multi-source data fusion is the main way to boost the accuracy of coal and rock recognition. The new features of coal and rock recognition are formed by the fusion of diverse signals such as current, vibration, acoustic emission, and infrared in the way of weight [16,17]. Liu et al. [18] developed a vibration feature extraction method based on Hilbert spectrum information entropy. Through the empirical mode decomposition of the vibration signal of the tail arm bracket of the excavator, it was found that the • It was theoretically analyzed that the fluctuation in load torque not only produces an amplitude-frequency modulation of the stator current but also produces phase modulation, which provides a theoretical basis for stator current as a condition for coal hardness identification.

•
The EEMD algorithm was used to decompose and reconstruct the stator current signal for highlighting the power frequency characteristics of the stator current. In addition, multi-scale permutation entropy was used to describe the weak changes of current from multiple scales. Finally, the multi-scale permutation entropy was identified by the Adaboost-BP network to improve the accuracy of coal rock hardness identification. • An adaptive speed control strategy of the roadheader based on coal hardness was proposed to ensure efficient cutting of the roadheader by adjusting the rotational speed of the cutting head and swing speed of the cutting arm.
The rest of the paper is organized in the following manner. Section 2 investigates the relationship between the cutting stator current and coal rock load. In Section 3, EEMD is applied to decompose and reconstruct the current signal according to the energy density and correlation coefficient criteria. In Section 4, the Adaboost-BP network is utilized to train the MPE eigenvalues derived from the reconstructed signal in order to establish a coal and rock hardness recognition model. In Section 5, to realize expeditious identification of coal rock hardness via the above methods, the Adaboost-BP neural network is used to identify the real-time collected current characteristics, and to control the optimized cutting parameters under a certain hardness by a servo valve and converter to achieve efficient and adaptive cutting of the roadheader.

Relationship between Load and Current
The variation in the load will influence the current value of the motor. The current is proportional to the load torque. The dynamic formula of the motor drive system is [28]: The relationship between the output torque and the load torque of the cutting head is: where T e and T L are the electromagnetic torque and output torque, respectively; J is the moment of inertia of the motor; p n is the pole number of the motor; ω m is the motor speed; Z i is the cutting resistance of the i-th pick in the working area, which is related to the rock hardness P k ; r i is the cutting radius of the i-th pick in the working area; µ is the reduction ratio of the reducer; and η is the transmission efficiency of the reducer. It can be seen from Equation (1) that T L and ω m are constant when the motor is stable. If the load torque M of the cutting head induces a periodic fluctuation T s in the output torque T L , the final steady state of the motor is in dynamic equilibrium. Corresponding fluctuations also occur in electromagnetic torque and motor speed [29]. Therefore, Fourier series expansion of each fluctuation component is carried out: where T L0 , T e0 , and ω m0 are the average values of output torque, electromagnetic torque, and motor speed, respectively; f i is the frequency corresponding to the periodic fluctuation component; T Li , T ei , and ω mi correspond to the amplitude of output torque, electromagnetic torque, and motor speed fluctuation. ϕ Li , ϕ ei , and ϕ mi correspond to the phase of the output torque, electromagnetic torque, and motor speed fluctuation. When the motor reaches a steady state, the direct current components and cosine component, respectively, match the following: where T Ls , T es , and ω ms are periodic fluctuations in each variable.
The electromagnetic torque equation of the motor with id = 0 vector control mode is: According to Equations (3) and (5), i q is the direct current components plus the cosine component of the corresponding frequency: where i q0 and i qi denote the average value of the q-axis current and the amplitude of the fluctuation component, respectively; ϕ qi is the phase. Transforming the current in the d-q coordinate system to the a-b-c three-phase coordinate system, the following is obtained: The rotor position information of the motor is: where θ 0 is the initial position of the rotor, f e0 is the fundamental frequency of the motor current, and θ tz is the modulation component corresponding to the angular velocity fluctuation. The expression of a-phase current can be obtained by Equations (6)- (8): the load torque impacts not only the current value, but also the motor current signal in the frequency domain.

EEMD Decomposition of Cutting Motor Current Signal
The IMF decomposed by EMD is clearer than the original signal. In addition, the frequency of IMFs decreases with the order, which is beneficial to the subsequent signal processing [30]. However, the EMD method has the mode mixing phenomenon, so this paper used the EEMD algorithm to decompose the stator current. The EEMD method is a step forward from the EMD method, which can decompose any nonlinear and nonstationary signal into the sum of IMFs and a residual value (r) [31,32]. The decomposition steps are as follows: Step 1: Add white noise n(t) with a mean value of 0 and a constant standard deviation to the original signal y(t) in order to obtain a noisy signal y (t): Step 2: Decompose the noise signal y (t) by EMD to obtain a set of denoised IMFs as c k (t) and a residual component r(t): where c k (t) is the k-th IMF after EMD decomposition of the noise-containing signal, and n is the number of IMFs.
Step 3: Repeat steps 1 and 2 N times. The total average of the corresponding IMF is derived on the principle that the statistical mean of the unrelated random series is 0. The IMF after EEMD decomposition is C k (t), which eliminates the influence of adding white noise to the real IMF, where k = 1, 2, 3, . . . , n.
where c kj (t) is the IMF obtained by EMD by adding white noise to the original signal y(t) for the j-th time.
In the first section of this paper, a simulation model of the cutting motor was constructed to obtain the stator current signal of the cutting motor, as shown in Figure 1. where ckj(t) is the IMF obtained by EMD by adding white noise to the original signal y(t) for the j-th time.
In the first section of this paper, a simulation model of the cutting motor was constructed to obtain the stator current signal of the cutting motor, as shown in Figure 1. The current signal of the cutting motor was decomposed by the EEMD method. The current signal was decomposed into 10 IMFs with high-to low-frequency distributions, as shown in Figure 2. IMF1 is the original signal and IMF10 is the residual component r. The motor stator current signal contains abundant characteristic information of the motor drive system. The amplitude variation of the stator current mainly occurred in the power frequency region. Combined with the 50 Hz power frequency characteristics of the current signal, it can be considered that the IMF component that has a strong correlation with the original current signal contained a strong 50 Hz signal. The decomposed IMF components were reconstructed by the correlation coefficient and energy density to further highlight the power frequency characteristics.  The current signal of the cutting motor was decomposed by the EEMD method. The current signal was decomposed into 10 IMFs with high-to low-frequency distributions, as shown in Figure 2. IMF1 is the original signal and IMF10 is the residual component r. The motor stator current signal contains abundant characteristic information of the motor drive system. The amplitude variation of the stator current mainly occurred in the power frequency region. Combined with the 50 Hz power frequency characteristics of the current signal, it can be considered that the IMF component that has a strong correlation with the original current signal contained a strong 50 Hz signal. The decomposed IMF components were reconstructed by the correlation coefficient and energy density to further highlight the power frequency characteristics.

Signal Reconstruction Principle Based on Energy Density and Correlation Coefficient Criterion
For signals mixed with random noise, the high-frequency IMF component is generally noise after decomposition, but low-frequency disturbance is also mixed in the low-frequency portion. Simple frequency analysis cannot eradicate noise, due to complicated work environments of the roadheader [33]. To cope with this problem, this paper selected the IMF component from the energy density and correlation coefficient angle. The steps are as follows: Step 1: Decompose the current signal of the cutting motor by EEMD to obtain N IMF components for energy analysis, and calculate the energy density of each IMF : where E k is the energy density of the k-th IMF component, K is the length of IMF, and x k (d) is the amplitude of the k-th IMF component.
Step 2: Analyze the correlation between the IMF component and the current signal of the cutting motor. The mathematical expression is: where R k denotes the mathematical expectation, µ k is the mean value of IMF components, µ is the mean value of the original current signal y, σ k is the standard deviation of the IMF component, and σ is the standard deviation of the original current signal y.
The energy density and correlation of IMF components decomposed by EEMD were analyzed. The statistical results of energy density and correlation are shown in Tables 1 and 2: According to the statistics, the energy densities of IMF2, IMF6, and IMF7 were all greater than 0.3, which significantly exceeded other-order IMFs. According to the correlation coefficient criterion [34], IMF components of order 2, 3, 4, 5, 8, and 9 with a correlation coefficient less than 0.2 are not correlated with the original current signal. IMF6 and IMF7 were selected as the target signal, and the current signal was reconstructed by superimposing them, as shown in Figure 3. Compared to the original current signal, the reconstructed current signal was smoother, which effectively filters out the high-frequency noise and low-frequency disturbances in the original signal.

Calculation of Multi-Scale Permutation Entropy
The permutation entropy algorithm can effectively enlarge the weak change in time series, and multi-scale permutation entropy can describe the characteristics of time series from multiple scales. In this paper, the multi-scale permutation entropy of stator current was used as the identification feature to accurately discern the hardness of coal rock under complex working conditions. Multi-scale permutation entropy coarse-grains the signal, calculates the permutation entropy of the coarse-grained fragments, and realizes the signal description in multiple dimensions [35]. Coarse-graining the time series X = {x e , e = 1, 2, 3 . . . , F} with the length of N to obtain the coarse-grained sequence y (s) j , the expression is: where s is the scale factor and s = 1, 2, . . . [F/s] denotes rounding F/s. When s = 1, the coarse-grained sequence is the original sequence. The multi-scale permutation entropy can be obtained by calculating the permutation entropy of each coarse-grained sequence. Figure 4 depicts the calculation process of the multi-scale permutation entropy.
where s is the scale factor and s = 1, 2, … [F/s] denotes rounding F/s. When s = 1, the coarsegrained sequence is the original sequence. The multi-scale permutation entropy can be obtained by calculating the permutation entropy of each coarse-grained sequence. Figure 4 depicts the calculation process of the multi-scale permutation entropy. The selection of scale factor is critical in coarse graining. If the scale factor s is too small, the characteristics of the signal cannot be retrieved to their full potential. If the value of s is too large, the discrepancy between the signals may be erased [36]. Therefore, the particle swarm optimization algorithm was employed in this paper to optimize the scale factor s. Skewness represents the degree of probability density a nonnormally distributed random sequence deviates from the normal distribution, which is a digital property of the degree of asymmetry of a statistical data distribution.
The permutation entropy under all scales of the time series X = {xe, e = 1, 2, 3…, F} is composed of a sequence Hp(X) = {Hp(1), Hp (2), …, Hp(s)}, and skewness Ske can be expressed as: The scale factor s = 10 of the multiscale permutation entropy algorithm is obtained by optimizing the skewness of the reconstructed current signal. The embedding dimension m is the main parameter of the permutation entropy algorithm. The embedding dimension m determines the number of states ! m , and the value of permutation entropy is highly dependent on the choice of m [37]. If the value of the embedding dimension m is too small, the mutation detection performance of the algorithm is degraded in the process of calculating the multiscale permutation entropy. If the embedding dimension m is too  The selection of scale factor is critical in coarse graining. If the scale factor s is too small, the characteristics of the signal cannot be retrieved to their full potential. If the value of s is too large, the discrepancy between the signals may be erased [36]. Therefore, the particle swarm optimization algorithm was employed in this paper to optimize the scale factor s. Skewness represents the degree of probability density a nonnormally distributed random sequence deviates from the normal distribution, which is a digital property of the degree of asymmetry of a statistical data distribution.
The permutation entropy under all scales of the time series X = {x e , e = 1, 2, 3 . . . , F} is composed of a sequence H p (X) = {H p (1), H p (2), . . . , H p (s)}, and skewness Ske can be expressed as: where H p (X) is the mean value of the sequence H p (X), and H d p (X) is the standard deviation of the sequence H p (X).
The objective function is: The scale factor s = 10 of the multiscale permutation entropy algorithm is obtained by optimizing the skewness of the reconstructed current signal. The embedding dimension m is the main parameter of the permutation entropy algorithm. The embedding dimension m determines the number of states m!, and the value of permutation entropy is highly dependent on the choice of m [37]. If the value of the embedding dimension m is too small, the mutation detection performance of the algorithm is degraded in the process of calculating the multiscale permutation entropy. If the embedding dimension m is too large, the entropy value will not reflect the subtle changes in the time series. Figure 5 indicates the variation in the MPE of different coal hardness P k with the scale factor when the embedding dimension is 3-7. large, the entropy value will not reflect the subtle changes in the time series. Figure 5 indicates the variation in the MPE of different coal hardness Pk with the scale factor when the embedding dimension is 3-7.    under various hardness values may be clearly differentiated. The MPE value decreases with the coal hardness. The delay time t = 5, m = 3, and s = 10 were determined through the aforementioned research and analysis. The sliding window method was employed for multi-scale entropy extraction to enhance the real-time performance of coal and rock recognition. The rotation of the cutting head was about 1.3 s. The sudden change in coal hardness resulted in a current change of about a 1/8 cycle. A window width of 250 ms, a window increment of 50 ms, and a window overlap rate of 80% were chosen. Table 3 indicates the average entropy values of the cutting current under different coal hardness values. Table 3 indicates that the average MPE varies greatly depending on coal hardness, and that the change trend of the MPE value with the scale factor is the same. MPE can be employed as an essential characteristic parameter of coal hardness through the above analysis.

Adaboost Improved BP Neural Network-Based Coal Hardness Estimation Algorithm
The idea of the Adaboost algorithm is to combine the outputs of multiple weak learners to generate effective predictions. The BP network is prone to falling into the local optimal solution, has a poor prediction ability, and has a limited generalization ability [38,39]. The BP neural network is regarded as a weak learner, and the Adaboost algorithm was applied to merge the output results of multiple BP weak learners to output more accurate prediction results. The process is as follows: Step 1: Data selection and network initialization. h groups of training data T = {(x 1 ,y 1 ), (x 2 ,y 2 ), . . . , (x h ,y h )} are randomly selected from the sample space to initialize the distribution weight of the test data D i = 1/h (i = 1, 2 . . . , h). The neural network structure is determined according to the sample input and output dimensions, and the BP neural network weights and thresholds are initialized.
Step 2: Weak learner prediction. The BP neural network is trained with the training data. The output of the training data is predicted to obtain the prediction error e t of the prediction sequence g(t). The calculation of the error e t is: Step 3: Calculate the weight C t of the t-th BP learner.
Step 4: Modify the weight D t+1 (i) of the next round of training samples according to the weight C t of weak learning.
where B t is the normalization factor.
Step 5: Integrate strong learners. After T iterations of training, T groups of weak learner functions f (x) and their weight vectors are generated and integrated into a strong learner function by weighted manner.
where F(x) is the strong learner function and f t (x) is the t-th weak learner function. The input layer of the BP neural network is 10 and the output layer is 1. Figure 6 shows the influence of the hidden layer, learning rate, and the number of BP weak learners in the Adaboost strong learner model on network training error. Figure 6 indicates that when the number of hidden layer nodes is 8, the training error is 0.007188; when the learning rate is 0.05, the training error is 0.0076589; when the number of weak learners is 7, the training error is 0.0036743. The number of hidden layer nodes, the learning rate, and the number of weak learners of the BP network are determined to be 8, 0.05, and 7.
Step 5: Integrate strong learners. After T iterations of training, T groups of weak learner functions f(x) and their weight vectors are generated and integrated into a strong learner function by weighted manner.
where F(x) is the strong learner function and ft(x) is the t-th weak learner function. The input layer of the BP neural network is 10 and the output layer is 1. Figure 6 shows the influence of the hidden layer, learning rate, and the number of BP weak learners in the Adaboost strong learner model on network training error. Figure 6 indicates that when the number of hidden layer nodes is 8, the training error is 0.007188; when the learning rate is 0.05, the training error is 0.0076589; when the number of weak learners is 7, the training error is 0.0036743. The number of hidden layer nodes, the learning rate, and the number of weak learners of the BP network are determined to be 8, 0.05, and 7.    Figure 7 shows the training status of the BP model and Adaboost-BP model. Figure 7a indicates that the BP model iterates 150 times with an error of 10 −2 magnitude, while the Adaboost-BP model iterates 20 steps with an error of 10 −2 magnitude. Therefore, the Adaboost-BP model outperforms the BP model in terms of convergence speed and network training error. The regression coefficient R indicates how well the network fits the data. Figure 7c,d demonstrate that the Adaboost-BP model outperforms the BP model in fitting data samples.
The relative root-mean-square error (RMSE) was applied to determine the dispersion degree of prediction results, and the coal hardness prediction model was evaluated. (22) where y (i) is the estimated hardness of coal and rock, y(i) is the actual hardness value of coal rock, and T is the sample size.
Quantitative analysis was performed on the estimation effects of the six types of coal and rock hardness. The RMSE results of the BP model and Adaboost-BP model are shown in Table 4. The RMSE value of the Adaboost-BP model was smaller than that of the BP model, and the average value decreased by 0.0633. Therefore, the estimated effect of the Adaboost-BP model was significantly higher than that of the BP model. Figure 8 indicates the stator current of the cutting motor under the change in coal rock hardness from P k = 350 to P k = 1000. In Figure 8, the hardness of coal and rock in the starting stage was 0 MPa before 1.5 s and changed from 350 MPa to 1000 MPa after 2.8 s. BP and Adaboost-BP network were used to identify the multi-scale permutation entropy of current and predict the hardness of coal. The result is shown in Figure 9.
In Figure 9, the coal hardness estimated by the Adaboost-BP model was more consistent with the changing trend of the actual coal hardness. Therefore, the coal and rock hardness estimation model based on the Adaboost-BP neural network was superior to the BP model in coal hardness identification. The relative root-mean-square error (RMSE) was applied to determine the dispersion degree of prediction results, and the coal hardness prediction model was evaluated.
where  (i) y is the estimated hardness of coal and rock, y(i) is the actual hardness value of coal rock, and T is the sample size.
Quantitative analysis was performed on the estimation effects of the six types of coal and rock hardness. The RMSE results of the BP model and Adaboost-BP model are shown in Table 4. The RMSE value of the Adaboost-BP model was smaller than that of the BP model, and the average value decreased by 0.0633. Therefore, the estimated effect of the Adaboost-BP model was significantly higher than that of the BP model.    Figure 8 indicates the stator current of the cutting motor under the change in coal rock hardness from Pk = 350 to Pk = 1000. In Figure 8, the hardness of coal and rock in the starting stage was 0 MPa before 1.5 s and changed from 350 MPa to 1000 MPa after 2.8 s. BP and Adaboost-BP network were used to identify the multi-scale permutation entropy of current and predict the hardness of coal. The result is shown in Figure 9.   In Figure 9, the coal hardness estimated by the Adaboost-BP model was more consistent with the changing trend of the actual coal hardness. Therefore, the coal and rock hardness estimation model based on the Adaboost-BP neural network was superior to the BP model in coal hardness identification.

Adaptive Speed Control Based on Coal Rock Hardness Change
The purpose of hardness identification is to realize cutting arm and cutting head speed regulation. The speed regulation method of constant power ensures that the cutting motor runs at a constant power but its cutting efficiency is not the highest. The cutting motor is often overloaded or underloaded due to the load fluctuation on the cutting head.

Adaptive Speed Control Based on Coal Rock Hardness Change
The purpose of hardness identification is to realize cutting arm and cutting head speed regulation. The speed regulation method of constant power ensures that the cutting motor runs at a constant power but its cutting efficiency is not the highest. The cutting motor is often overloaded or underloaded due to the load fluctuation on the cutting head. In order to improve the cutting performance under different coal hardness values, the fluctuation coefficient of the cutting head load K R a , K R b , K R c , and K M t and the specific energy consumption H W were used as optimization targets, the swing (v) and rotation (n) speeds were selected as the optimized variables, and the cutting motor power and speed range were used as constraints. For optimization methods and parameters, please refer to the literature [40], and the optimization results are shown in Tables 5 and 6.  It can be seen from Table 5 that the cutting performance significantly improved after optimization. In order to ensure that the roadheader was in an efficient cutting state under different coal rock hardness values, an adaptive speed control strategy based on coal rock hardness was proposed, as shown in Figure 10. EEMD and reconstruction were conducted on the real-time collection of cutting motor current. The hardness of coal was evaluated by the neural network, and the optimal cutting parameters under this hardness were matched according to the optimization results. The matching cutting parameters were employed as control signals. The swing and rotation speeds were controlled through the electrohydraulic servo valve (EHSV) and the converter, to ensure that the cutting performance of the roadheader was always in the ideal cutting state.
It can be seen from Table 5 that the cutting performance significantly improved after optimization. In order to ensure that the roadheader was in an efficient cutting state under different coal rock hardness values, an adaptive speed control strategy based on coal rock hardness was proposed, as shown in Figure 10. EEMD and reconstruction were conducted on the real-time collection of cutting motor current. The hardness of coal was evaluated by the neural network, and the optimal cutting parameters under this hardness were matched according to the optimization results. The matching cutting parameters were employed as control signals. The swing and rotation speeds were controlled through the electro-hydraulic servo valve (EHSV) and the converter, to ensure that the cutting performance of the roadheader was always in the ideal cutting state. Hardness Figure 10. Adaptive variable speed cutting control strategy.

Adaptive Speed Control Simulation
Matlab and AMESim were applied to simulate the cutting speed regulation system of the roadheader to verify the feasibility of the adaptive cutting speed regulation strategy of the roadheader based on coal rock hardness. According to Table 6, the speed of the cutting head was reduced from 48.28 r/min to 36.9 r/min, and the swing velocity of the cantilever decreased from 2.43 m/min to 1.5 m/min, when the hardness of coal changed from Pk = 350 to Pk = 1000. The speed of the cutting motor was controlled by Space Vector Pulse Width Modulation (SVPWM). The transmission ratio of the cutting head was n1 = 32. The hydraulic system adopted electro-hydraulic servo valve control. Comparing the Figure 10. Adaptive variable speed cutting control strategy.

Adaptive Speed Control Simulation
Matlab and AMESim were applied to simulate the cutting speed regulation system of the roadheader to verify the feasibility of the adaptive cutting speed regulation strategy of the roadheader based on coal rock hardness. According to Table 6, the speed of the cutting head was reduced from 48.28 r/min to 36.9 r/min, and the swing velocity of the cantilever decreased from 2.43 m/min to 1.5 m/min, when the hardness of coal changed from P k = 350 to P k = 1000. The speed of the cutting motor was controlled by Space Vector Pulse Width Modulation (SVPWM). The transmission ratio of the cutting head was n 1 = 32.
The hydraulic system adopted electro-hydraulic servo valve control. Comparing the speed control results under the Adaboost-BP and BP prediction models, the simulation results are shown in Figure 11. 32. The hydraulic system adopted electro-hydraulic servo valve control. Comparing the speed control results under the Adaboost-BP and BP prediction models, the simulation results are shown in Figure 11.  Figure 11 indicates that for the cutting motor or the hydraulic system, the speed regulation results based on the Adaboost-BP estimation model were basically consistent with  Figure 11 indicates that for the cutting motor or the hydraulic system, the speed regulation results based on the Adaboost-BP estimation model were basically consistent with the optimized results. The speed control results based on the BP estimation model fluctuated and deviated significantly from the optimization results. The adaptive speed control based on the Adaboost-BP estimation model can ensure that the roadheader is in the optimal cutting state.

Conclusions
In order to enhance the accuracy and efficiency of coal rock hardness identification, the multi-scale permutation entropy of stator current was determined as the identification trait by analyzing the relationship between stator current and load. The Adaboost-BP neural network was utilized to train the characteristics to construct the coal rock hardness estimation model. The following conclusions can be drawn from the results in this study.
The reconstruction method based on energy density and the correlation coefficient criterion could effectively filter out the noise interference in the signal. The reconstructed stator current signal was extracted by multi-scale permutation entropy after parameter optimization. The results indicated that when the delay time was 5, the embedding dimension was 3, and the scale factor was 10, the multi-scale permutation entropy could better distinguish the hardness of coal and rock.
The Adaboost-BP prediction model was applied to identify the hardness of coal and rock. The Adaboost-BP prediction model proposed in this paper outperformed the BP prediction model in terms of hardness prediction. Simultaneously, the speed regulation method based on coal rock hardness could effectively regulate the speed of the cutting head and the cutting arm.
In this paper, we proposed an efficient hardness identification approach for the detection of coal and rock hardness. However, it should be realized that using the stator current characteristics as the only input signal for identification may be insufficient. It is worth noting that multi-source data fusion technology [41] can overcome this deficiency because this technology can enrich the feature information of the identified object. In addition, the Adaboost-BP neural network has defects such as long training time and difficulty in determining the number of weak learners. The Deep Temporal Convolution Network (DTCN) has a higher efficiency in training samples and data classification [42]. Therefore, combining deep learning technology and multiple data fusion technology to identify the hardness of coal and rock will be interesting work for our future studies.