Elevator Car Vibration Signal Denoising Method Based on CEEMD and Bilateral Filtering

Elevator car vibration signals are important information to monitor and diagnose the operating status of elevators, but during the process of vibration signals acquisition, vibration signals are always inevitably disturbed by noise, which affects further research. Therefore, aiming at the vibration signal with noise, we propose a new vibration signal denoising method on the basis of complementary ensemble empirical mode decomposition (CEEMD) and bilateral filtering. Firstly, the collected original vibration signals are decomposed by the CEEMD into several inherent mode functions. Then, the false components are removed by determining the correlation coefficients of mode components, and then the noise-dominate components are denoised by bilateral filtering. Finally, the processed inherent mode functions are reconstructed to require the denoised signal. We test the method through simulation and practical application. The results indicate that the proposed method can efficaciously reduce the noise in the vibration signal of an elevator car.


Introduction
With the development of the city, various high-rise buildings continue to appear. As a kind of transportation equipment, the elevator plays a vital role in these modern buildings, but it is inevitable that some faults occur in the process of putting the elevator into operation. If we can take measures to diagnose the failure of the elevator in a shorter time, we can reduce the losses. Vibration signals are often used to judge the fault type, but the vibration signals we generally collect are mixed with various noises, so we need to eliminate noise on the premise of ensuring that the pure signals are not affected.
Fourier transform is a frequently used method to denoise vibration signals. Under the influence of the Fourier transform, the vibration signal can shift from the time-domain state with time as the abscissa to the frequency-domain state with frequency as the abscissa. The process is actually the decomposition process of the signal, and the original signal finally becomes the combination of multiple sine waves. The premise of the Fourier transform application is that the signal processed is the stationary signal, but the signal we usually obtain is not a stationary signal, so it is difficult to obtain satisfactory results by using Fourier transform. Different from Fourier transform, wavelet transform can place the signal on the coordinate axis with time as the abscissa and frequency as the ordinate. Based on the wavelet transform theory of vibration signal, many denoising methods are proposed. The modulus maximum denoising method [1] calculates the modulus maximum of the signal after wavelet transformation, selects an appropriate threshold, retains the wavelet coefficient of the corresponding extreme points, and carries out signal reconstruction to achieve denoising. The spatial correlation denoising method [2] takes the scale correlation of signals as the key factor of denoising and determines whether to retain the coefficients according to the similarity between scales. The wavelet threshold denoising [3] sets the judgment threshold based on experience or other methods. By comparison and analysis, the coefficients smaller than the threshold are removed, and the coefficients larger than the threshold are retained. Finally, the purpose of denoising is achieved through wavelet inverse transform.
With the development of wavelet denoising, wavelet packet denoising methods have emerged, and there are more and more improved methods based on wavelet packet denoising. Li et al. [4] realized the signal denoising by setting the threshold to a variable value on the basis of wavelet packet decomposition of the signal. Chen et al. [5] combined sparse knowledge with wavelet packet decomposition to remove noise in the vibration signal, and the sparse knowledge used can restore wavelet packet coefficients of the pure signal to a large extent. Regardless of wavelet denoising or wavelet packet denoising, it is necessary to choose the decomposition layers and wavelet basis, but the selection of these two parameters usually needs to be based on experience once improper selection will affect the final denoising effect of the vibration signal.
In addition to wavelet analysis, Huang et al. [6] presented empirical mode decomposition (EMD) for the first time, which can represent a signal as the sum of multiple modal components and a residual component. Later, this method was applied in vibration signal denoising [7]. However, from the perspective of practical application effect, if signal is decomposed by using empirical mode decomposition, there are false components and aliasing components in the obtained components. Based on the problems existing in empirical mode decomposition, Wu et al. [8] put forward ensemble empirical mode decomposition (EEMD). The method greatly reduces the aliasing components in the decomposition process, but the operation of adding white noise makes the residual noise still exist in the reconstructed signal. Later, Yeh et al. [9] suggested complementary ensemble empirical mode decomposition (CEEMD), which greatly reduces the aliasing components. After that, a series of the improved method has been proposed to overcome the shortcomings of EMD and its improved algorithms. Dang et al. [10] introduced the method of using EMD and wavelet transform simultaneously in vibration signal denoising. Chegini et al. [11] combined EMD and wavelet analysis and applied them to denoise in bearing vibration signals. Because wavelet packet transform can achieve better analysis of medium-and high-frequency components of signals, the vibration signal denoising method combining CEEMD and wavelet packet according to the different characteristics of noise and pure signals in vibration signals was proposed in [12][13][14][15][16][17]. The above improved methods mainly combined empirical mode decomposition or its improved algorithms with wavelet decomposition or wavelet packet decomposition for denoising. The denoising process would inevitably involve the selection of wavelet base and the determination of decomposition layers. Therefore, there might be the same problems with wavelet denoising and wavelet packet denoising. In 2014, Dragomiretskiy and Zosso [18] proposed the VMD algorithm, which can decompose the signal into AM and FM component signals. Later VMD algorithm was applied to signal denoising. Long et al. [19] proposed a denoising method combining the VMD algorithm and wavelet decomposition. Yu et al. [20] applied the VMD algorithm to seismic signal denoising. The algorithm has better decomposition precision for complex data, but the number of decomposition layers and penalty factors need to be set in advance before using VMD. Once these two parameters are set improperly, the final decomposition effect will be greatly affected.
At present, Image denoising methods commonly used include gaussian filtering, median filtering, average filtering, and so on. These methods are denoising methods for two-dimensional signals, and some of them have a good filtering effect on noise points. Therefore, we can try to transform these methods with better denoising performance to make them suitable for one-dimensional vibration signals. The Gaussian filtering method has a good denoising effect for random Gaussian noise in the image, but it will affect the edge information of the image in the process of denoising so that the final image after denoising is not clear enough. The median filtering method has a good denoising effect on isolated noise and does not affect the edge information of the image. However, the method needs to determine the size of the sliding window according to the actual situation. Once the improper selection is made, the final image after denoising will become blurred. Average filtering affects the restoration of image edges while removing noise. Bilateral filtering is an evolutionary algorithm of gaussian filtering, which not only considers spatial proximity but also involves the similarity of pixel values. The method is a kind of nonlinear filtering technology that can realize edge preserving and denoising. At present, many images are denoised using a bilateral filtering algorithm [21][22][23].
The noise in the elevator car comes from the car itself and the external environment, and its frequency is mainly medium and low frequency. For the noise filtering in the elevator signal, this paper proposes a denoising method combining bilateral filtering and CEEMD. In the new method, CEEMD and bilateral filtering are combined for the first time to denoise the elevator car vibration signal. Firstly, CEEMD is used to decompose the elevator vibration signal, and the decomposition result is multiple modal components and a residual component. The paper mainly carries out further processing and analysis of the decomposed modal components. The correlation coefficient and multi-scale arrangement entropy are used in the processing, and then the components dominated by noise and the components dominated by the signal are determined. Because bilateral filtering has a good filtering effect on medium-and low-frequency signal, it is applied to the components dominated by noise to remove the noise in the elevator vibration signal. Then, the signal is reconstructed to obtain the denoised signal. Finally, the method is applied to the simulation signal and the elevator vibration signal, respectively, which verifies its superiority and shows that it has a better denoising effect.
The overall structure of this paper is as follows: Section 2 introduces empirical mode decomposition and its improved algorithm and variational mode decomposition. Section 3 introduces the method we suggest in this paper in detail. Section 4 applies the method to the simulation analysis of elevator car vibration signals. Section 5 applies the method, EMD denoising, EEMD denoising, CEEMD denoising, bilateral filtering denoising, and VMD denoising to the elevator car vibration signal and compares the final denoising effect of different methods. The last part draws the conclusion of the whole paper.

Empirical Mode Decomposition (EMD)
The formula of empirical mode decomposition is as follows: where im f i is ith inherent modal function, and r(t) is the residual component. Each modal component needs to meet the following requirements [24]: (1) The difference between the number of local extremum points and zeros crossing of the component is 0 or 1; (2) The upper envelope curve and the lower envelope curve in the component are composed of maximum points and minimum points, respectively, and one-half of their sum equals zero.
The steps of empirical mode decomposition are as follows: Step 1: Determine all maximum and minimum points of the original signal x(t).
Step 2: According to the maximum and minimum points found in Step1, selecting the cubic spline function to fit them respectively requires the upper envelope U(t) and the lower envelope D(t), and the mean value m(t) of the two envelopes can be obtained by Formula (2).
Step 3: Remove the mean envelope m(t) from the input signal x(t). The final obtained value is denoted as h 1 (t), h 1 (t) = x(t) − m(t). If h 1 (t) satisfies the condition that the modal component must have, h 1 (t) is an inherent modal function, denoted as c 1 (t), and proceed with step 4. Otherwise, let x(t) = h 1 (t), and proceed with step1-step3.
Step 4: Calculation of residual signal r 1 (t): Step 5: Let x(t) = r 1 (t), and begin execution from step1 until the next modal component c 2 (t) is found. After n times of repeated iterations, all modal components and one residual component r n (t) are finally obtained:

Ensemble Empirical Mode Decomposition (EEMD)
Ensemble empirical mode decomposition is a modified algorithm for the shortcomings of EMD. If the signal varies greatly in different time scales, the modal components obtained by empirical mode decomposition may contain more than the information of the current time scale, which is called mode aliasing. In order to reduce this phenomenon, EEMD adds white noise on the basis of unprocessed signals. Because the added signal is characterized by uniform spectrum distribution, all the signals obtained after ensemble empirical mode decomposition will be placed on the corresponding time scale. The steps of Ensemble empirical mode decomposition are shown below [8]: Step 1: Add white noise to the original signal; Step 2: The signal after processing in Step 1 is disassembled into a combination of IMF components by empirical mode decomposition (EMD); Step 3: Repeat Step1 and Step2 n times; Step 4: Take the IMF aggregate mean obtained each time as the final result.

Complementary Ensemble Empirical Mode Decomposition (CEEMD)
CEEMD improves the modal aliasing problem in the signal processing process. The difference between this method and EEMD is that the latter adds multiple white noises to the signal before decomposition, while the former adds pairs of white noise with opposite signs to the signal before decomposition, which can not only decrease the effect of noise on the decomposition result but also enhance the computational efficiency.
The main steps of CEEMD are as follows [25]: Step 1: Add a pair of noise signals to the input signal, and the added white noise has the same amplitude.
where X is the original input signal, and N i is the white noise added for the ith time.
Step 2: After adding white noise for the ith time, EMD decomposition of M i 1 and M i 2 is performed to require the corresponding components c i 1,j and c i 2,j , as well as the corresponding residuals r i 1 and r i 2 .
Step 3: Take the average values of modal components and residuals, respectively.

of 21
Step 4: Similarly, if n pairs of white noise with the opposite sign are added, the final IMF is:

Variational Mode Decomposition (VMD)
The core of variational mode decomposition is to construct and solve variational problems, which is an adaptive, completely non-recursive signal processing method. The method determines the frequency center and bandwidth of each component by iteratively searching for the optimal solution of the variational model in the process of obtaining the decomposed components so that the frequency domain division of the signal and the effective separation of each component can be adaptively realized.
The variational mode decomposition can decompose the signal into multiple amplitudefrequency modulated signals, and the decomposition formula is as follows: where f (t) is the original signal and u k (t) is the kth harmonic signal, A k (t) is the amplitude of the kth harmonic signal, ϕ k (t) is the phase of the kth component, and K is the number of modal components obtained by decomposition.
The specific decomposition steps of the variational modal decomposition algorithm are as follows [26,27]: (1) Hilbert transform for each modal signal is as follows: where δ(t) is the dirac distribution, j is the imaginary part, * is the convolution symbol.
(2) The spectrum of each mode is modulated to the corresponding fundamental frequency band.
where w k represents the center frequency of the kth mode component.
(3) Calculate the square root of the L2 norm gradient of the demodulated signal and estimate the bandwidth of each modal signal. The variational problem is expressed as follows: where ∂ t represents the partial derivative of the function with respect to time, δ(t) is the unit pulse function. The quadratic penalty factor α and the lagrangian penalty operator λ(t) are introduced to convert the constrained problem into the unconstrained problem. The extended Lagrangian expression is as follows: where α represents the bandwidth parameter. The Alternating Direction Method of Multipliers (ADMM) is used to solve the variational problem. By updating u n+1 k , w n+1 k and λ n+1 , we can find the Lagrangian saddle point, obtain the optimal solution of u n+1 k (w) and w n+1 k (w), and substitute them into ADMM to obtain a complete VMD algorithm. The specific steps are as follows: Step 1: Initialize u 1 k , w 1 k , λ 1 and initializes n to 0.
where τ represents the noise tolerance parameter. When the signal contains strong noise, we can set τ = 0.
Step 4: Repeat Step 2 and Step 3 until the following iteration constraints are satisfied.

Vibration Signal Denoising Method Based on CEEMD
In the paper, we take advantage of the CEEMD to decompose the vibration signal into several intrinsic mode functions. Then remove the false components by calculating the correlation coefficient of the IMFs. Because there is still noise in the remaining signal, the remaining signal is divided into two parts by multi-scale permutation entropy, and then the signal in the dominant position of noise is denoised by bilateral filtering. Finally, the processed signal is refactored to require the final denoising Signal. Figure 1 shows the vibration signal denoising flow chart combining CEEMD and bilateral filtering.
The main steps of the new method are as follows: (1) Let the noisy signal be y(t), and divide the signal into several intrinsic modal functions and one residual component through CEEMD. vibration signal. Therefore, considering the results of multiple denoising experiments for the elevator signal, we set the threshold of the correlation coefficient to 0.1~0.2 of the maximum correlation coefficient. Because the paper does not carry out experimental verification on other equipment except the elevator, the current setting of this threshold is only applicable to the vibration signal of the elevator car. δ in Figure 1 indicates the threshold of the correlation coefficient.
where ρ i is the correlation coefficient between the im f i and the input signal, |ρ i | ≤ 1, The main steps of the new method are as follows: (1) Let the noisy signal be ( ), and divide the signal into several intrinsic modal functions and one residual component through CEEMD.  (3) Determine the multi-scale permutation entropy of the remaining IMF components and divide the signals into the signal-dominated IMF components and the noise-dominated IMF components. Multi-scale permutation entropy is the permutation entropy of a signal at different scales. Multi-scale permutation entropy can represent the characteristics of signals at different scales, determine the complexity and randomness of signals, and highlight the small changes in signals. It is an important parameter for analyzing and processing vibration signals generated by mechanical equipment. The process is to coarse-grain the original time series to construct a multi-scale time series and then calculate the permutation entropy at each scale. The smaller the permutation entropy is, the more regular the time series is. The higher the permutation entropy, the more complex the time series. Therefore, we calculate the multi-scale permutation entropy of IMF components to determine their categories.
The calculation steps of multi-scale arrangement entropy of each modal component are as follows: Let the time series of each modal component be X = {x(i), i = 1, 2, 3, . . . . . . , n}, where n is the number of sampling points.
Step 1: The time series are coarsely granulated to obtain the processed sequences.
where s is the scale factor, and y s (j) is the time series under different scale factors.
Step 2: Reconstruct the time series y s (j).
where m stands for embedding dimension, and τ stands for delay time. If the above formula is sorted in ascending order, there are m! permutations in multi-scale order. The probabilities of time series at each scale are calculated: where K represents the number of occurrences of each type, and the permutation entropy of each scale is: Normalize the above equation: The calculated multi-scale permutation entropy of the IMF component is denoted as MPE(i), i = 1, 2, . . . , q. q represents the number of IMF components.
If MPE(i) is large, it indicates that the IMF component corresponding to this value is a noise-dominated component. If MPE(i) is small, the IMF component corresponding to this value is the signal-dominated component. The experimental object of this paper is the elevator car vibration signal. After many tests, it is more appropriate to set the threshold value of MPE between 0.7 and 0.8. ε in Figure 1 indicates the threshold of the MPE. Similar to δ, the value of ε we selected this time is only applicable to the elevator car vibration signal.
(4) Keep the signal-dominated components and remove noise from components containing a lot of noise by bilateral filter.
Bilateral filtering is a nonlinear filter, which is usually used to process image noise. When processing adjacent pixel values, it takes into account the proximity relationship in the distance and the similarity in gray scale at the same time. It realizes adaptive filtering by nonlinear combination of spatial proximity and pixel value similarity [28]. The expression of bilateral filter is as follows [29]: where g(u) is the image after bilateral filtering, M i,j indicates the pixel set of (2N + 1) × (2N + 1) spatial neighborhood centered on (i, j), f (x, y) indicates the pixel value of (x, y) in M i,j , w d (i, j) is the spatial proximity Gaussian function, w r (i, j) is the pixel value similarity Gaussian function, σ d is the spatial distance standard deviation, and σ r is the numerical similarity standard deviation. The signal studied in the paper is one dimension vibration signal. In order to illustrate how to apply bilateral filtering to vibration signal, the following example is given [30]: Given a vibration signal Y with noise, where S is the noiseless signal and N is the noise signal.
For this vibration signal, the bilateral filtering formula is changed as follows, which represents the normalized weighted average value of neighborhood points with the size of 2R + 1: where γ[ξ − x] represents the amplitude of the vibration signal at ξ − x, W[ξ, x] is the weight coefficient of bilateral filtering, which is the product of W d [ξ, x] and W r [ξ, x]. The formulas for the two parameters are as follows: (5) Finally, reconstruct the processed signal to obtain the denoising signal y (t). According to the above description, Algorithm 1 summarizes the new denoising method proposed in the paper.

Algorithm 1: The workflow of new method
Input: y(t), δ, ε Output: y (t) Procedure: Decompose y(t) to n modes by CEEMD the decomposed modes are denoted as im f i (t) (i = 1, 2, . . . . . . n) For i = 1 : n Determine the correlation coefficient between im f i (t) and the original signal the correlation coefficient is denoted as Determine the multi-scale arrangement entropy of im f i (t) the multi-scale arrangement entropy is denoted as MPE(i) If MPE(i) > ε() Bilateral filtering is needed to denoise Else Remain original im f i (t) Reconstruct the denoised signal y (t).

Simulation Signal Verification
Because the frequency of elevator vibration signal is usually 0.5~80 Hz, we set an analog signal to simulate its vibration signal and add Gaussian white noise to the analog signal as the noise in the signal. The formula for obtaining the signal is as follows: where Y is the signal with added noise, X is the original signal, N is the gaussian white noise signal, and A is a constant value. amp(i) is the corresponding amplitude of vibration signal with frequency of f (i), respectively: 13, 3, 4, 8, 5, 7, 4, and 3; f is the frequency of vibration signals, respectively: 10, 20, 30, 40, 50, 60, 70, and 80. The time domain signal diagram of the original signal and the noisy signal is shown in Figure 2, and the frequency domain signal diagram is shown in Figure 3.  EMD, EEMD, and CEEMD decomposition are performed on the noisy signal, and time domain diagram of each IMF is finally obtained, as shown in Figure 4. As can be s from Figure 4, the IMF components obtained after EMD and EEMD decomposition alw have the problem of modal aliasing, while the components obtained after CEEMD dec  EMD, EEMD, and CEEMD decomposition are performed on the noisy signal, and time domain diagram of each IMF is finally obtained, as shown in Figure 4. As can be s from Figure 4, the IMF components obtained after EMD and EEMD decomposition alw have the problem of modal aliasing, while the components obtained after CEEMD dec position only contain one time scale characteristic component from IMF4, which ef     As can be seen from Figure 5, EMD, EEMD, CEEMD, bilateral filtering, VMD, and our newly proposed method can reduce the signal amplitude of medium and high frequency, that is, reduce or eliminate the medium-and high-frequency noise, but the top five methods will greatly reduce the signal amplitude at the main frequency, resulting in the loss of useful information in the signal. The method proposed in this paper can ensure that the amplitude of the main frequency signal is kept unchanged or slightly increased to the maximum extent, which is conducive to the further analysis and processing of the signal, and with the increase in noise frequency, the noise reduction effect of vibration signal obtained by using the method proposed in this paper is more obvious. Therefore, in general, the method proposed in this paper has a better denoising effect and is suitable for elevator signal denoising.
In order to express the denoising effect of different methods more intuitively in the form of data, we employ two evaluation indicators to judge the denoising effect. The evaluation indicators involved include signal-to-noise ratio and root mean square error. The SNR is equal to the effective power of the signal divided by the effective power of the noise, and the greater the value of SNR, the better the denoising effect. The RMSE is the square root of the variance between the original signal and the denoised signal, and the smaller the value of RMSE, the stronger the denoising ability.
where N indicates the number of data, indicates pre-noised signals, and indicates post-noised signals.
The final denoising evaluation index value is shown in Table 1. It can be seen from the table that compared with the other five methods, the signal-to-noise ratio of the As can be seen from Figure 5, EMD, EEMD, CEEMD, bilateral filtering, VMD, and our newly proposed method can reduce the signal amplitude of medium and high frequency, that is, reduce or eliminate the medium-and high-frequency noise, but the top five methods will greatly reduce the signal amplitude at the main frequency, resulting in the loss of useful information in the signal. The method proposed in this paper can ensure that the amplitude of the main frequency signal is kept unchanged or slightly increased to the maximum extent, which is conducive to the further analysis and processing of the signal, and with the increase in noise frequency, the noise reduction effect of vibration signal obtained by using the method proposed in this paper is more obvious. Therefore, in general, the method proposed in this paper has a better denoising effect and is suitable for elevator signal denoising.
In order to express the denoising effect of different methods more intuitively in the form of data, we employ two evaluation indicators to judge the denoising effect. The evaluation indicators involved include signal-to-noise ratio and root mean square error. The SNR is equal to the effective power of the signal divided by the effective power of the noise, and the greater the value of SNR, the better the denoising effect. The RMSE is the square root of the variance between the original signal and the denoised signal, and the smaller the value of RMSE, the stronger the denoising ability.
where N indicates the number of data, x i indicates pre-noised signals, andx i indicates post-noised signals. The final denoising evaluation index value is shown in Table 1. It can be seen from the table that compared with the other five methods, the signal-to-noise ratio of the denoised signal obtained by the method proposed in the paper is large, and the root mean square error is small, indicating that the denoising effect of this method is better. In addition, in order to verify the denoising effect of CEEMD-BF under different SNR, we add different levels of noise to the signal to obtain the SNR and root mean square error of the denoised signal. The results are shown in Table 2. As can be seen from Table 2, as the SNR increases, the SNR of the denoised signal becomes larger and larger, and the corresponding root mean square error becomes smaller and smaller. Therefore, the larger the SNR is, the better the denoising effect will be.

Elevator Vibration Signal Denoising
In order to verify the function of the method proposed in the paper in the denoising of elevator vibration signal, we collected the vibration signal of an elevator car. The equipment used to collect data in the experiment is LE-300, as shown in Figure 6. The signal acquisition frequency is 256 Hz, and the collected vibration signal in the horizontal direction of elevator is shown in Figure 7. Because the vibration signals of the X-axis and Y-axis of elevator vibration signal have little difference, the paper firstly focuses on the comparison and analysis of the vibration signals of the X-axis elevator and then verifies the denoising effect of the new method on the Y-axis signal. The number of points collected is 1024.
We first use EMD, EEMD, CEEMD, bilateral filtering, VMD denoising, and the method proposed in the paper to denoise from the X-axis vibration signal. The time domain diagram and frequency domain diagram after denoising are shown in Figures 8 and 9.
acquisition frequency is 256 Hz, and the collected vibration signal in the horiz tion of elevator is shown in Figure 7. Because the vibration signals of the X axis of elevator vibration signal have little difference, the paper firstly focuses parison and analysis of the vibration signals of the X-axis elevator and then denoising effect of the new method on the Y-axis signal. The number of poin is 1024. We first use EMD, EEMD, CEEMD, bilateral filtering, VMD denoising, and the method proposed in the paper to denoise from the X-axis vibration signal. The time domain diagram and frequency domain diagram after denoising are shown in Figure 8 and Figure 9. The horizontal vibration frequency of the elevator is generally 0.5~80 Hz. As can be seen from Figure 9a, the collected elevator vibration signals are mainly concentrated at 0.5~70 Hz. The vibration amplitude of the elevator car signal is larger at 5.25 Hz and 38.25 Hz, and other major frequency components are 18.5 Hz, 28.25 Hz, 51.25 Hz, and 65.25 Hz. As can be seen from Figure 9b-d, the denoising effects of vibration signals obtained by using EMD, EEMD, and CEEMD have no significant difference. When the frequency is close to 40 Hz, the amplitude of the vibration signal tends to 0, while the horizontal vibration frequency of the elevator vibration signal collected in the paper is 0.5~70 Hz. Therefore, these three denoising methods will lose some vibration information of the elevator. According to the time domain diagram in Figure 8d, it seems that the denoising effect of bilateral filtering is better than that of other methods. However, combined with the frequency domain diagram in Figure 9e, when the frequency is greater than 20 Hz, the amplitude of the main frequency decreases greatly. Therefore, in a comprehensive view, bilateral filtering filters out many useful signals and cannot ensure that the information of the original signal is not affected by the denoising operation. The VMD denoising in Figure 9f can remove the noise after 80 Hz, but has almost no effect on the noise before 80 Hz. There is no obvious difference between the time-domain graphs of the vibration signal before and after denoising in Figure 8e, so the denoising effect of the VMD de-noising method is not obvious and cannot be used for the denoising of elevator car vibration signals. Therefore, comprehensively, the method proposed in this paper can ensure that the amplitude at the main frequency of the signal is kept at a large value to the maximum extent, and the noise exceeding the normal vibration frequency range of the elevator can be effectively removed. The denoising results of the method proposed in this paper are good. We first use EMD, EEMD, CEEMD, bilateral filtering, VMD denoising, and the method proposed in the paper to denoise from the X-axis vibration signal. The time domain diagram and frequency domain diagram after denoising are shown in Figure 8 and Figure 9.  (e) (f) The horizontal vibration frequency of the elevator is generally 0.5~80 Hz. As can be seen from Figure 9a  We also applied the new method to the y-axis vibration signal of the elevator. The time-domain and frequency-domain diagrams before and after denoising are shown in Figures 10 and 11. From the time-domain and frequency-domain diagrams, we can see that this method is not only applicable to the x-axis vibration signal denoising but also applicable to the y-axis vibration signal denoising. It can retain the useful information in the signal to the maximum extent, remove the noise, and obtain a good denoising effect.
Sensors 2022, 22, x FOR PEER REVIEW 20 be effectively removed. The denoising results of the method proposed in this pape good.
We also applied the new method to the y-axis vibration signal of the elevator. time-domain and frequency-domain diagrams before and after denoising are show Figures 10 and 11. From the time-domain and frequency-domain diagrams, we can that this method is not only applicable to the x-axis vibration signal denoising but applicable to the y-axis vibration signal denoising. It can retain the useful informatio the signal to the maximum extent, remove the noise, and obtain a good denoising eff

Conclusions
In the paper, an elevator car vibration signal denoising method based on CEE and bilateral filtering is presented. Firstly, the method decomposes vibration signal several inherent modal functions by CEEMD. CEEMD can better solve the proble mode aliasing in signal decomposition. Next, the method determines the false com nents on the basis of computing the correlation coefficients. Then, the method determ the partial mean value of the multi-scale permutation entropy of the remaining com nents. If the partial mean value of IMF is large, it indicates that noise is the domi component in the IMF and needs bilateral filtering to remove the noise. Finally, the v tion signal after denoising is obtained through signal reconstruction. The effective and efficiency of this method are testified by simulation and experiment. Compared EMD denoising, EEMD denoising, CEEMD denoising, bilateral filtering denoising, VMD denoising, the method can not only retain the useful signal better but also filte be effectively removed. The denoising results of the method proposed in this paper are good.
We also applied the new method to the y-axis vibration signal of the elevator. The time-domain and frequency-domain diagrams before and after denoising are shown in Figures 10 and 11. From the time-domain and frequency-domain diagrams, we can see that this method is not only applicable to the x-axis vibration signal denoising but also applicable to the y-axis vibration signal denoising. It can retain the useful information in the signal to the maximum extent, remove the noise, and obtain a good denoising effect.

Conclusions
In the paper, an elevator car vibration signal denoising method based on CEEMD and bilateral filtering is presented. Firstly, the method decomposes vibration signal into several inherent modal functions by CEEMD. CEEMD can better solve the problem of mode aliasing in signal decomposition. Next, the method determines the false components on the basis of computing the correlation coefficients. Then, the method determines the partial mean value of the multi-scale permutation entropy of the remaining components. If the partial mean value of IMF is large, it indicates that noise is the dominant component in the IMF and needs bilateral filtering to remove the noise. Finally, the vibration signal after denoising is obtained through signal reconstruction. The effectiveness and efficiency of this method are testified by simulation and experiment. Compared with EMD denoising, EEMD denoising, CEEMD denoising, bilateral filtering denoising, and VMD denoising, the method can not only retain the useful signal better but also filter out

Conclusions
In the paper, an elevator car vibration signal denoising method based on CEEMD and bilateral filtering is presented. Firstly, the method decomposes vibration signal into several inherent modal functions by CEEMD. CEEMD can better solve the problem of mode aliasing in signal decomposition. Next, the method determines the false components on the basis of computing the correlation coefficients. Then, the method determines the partial mean value of the multi-scale permutation entropy of the remaining components. If the partial mean value of IMF is large, it indicates that noise is the dominant component in the IMF and needs bilateral filtering to remove the noise. Finally, the vibration signal after denoising is obtained through signal reconstruction. The effectiveness and efficiency of this method are testified by simulation and experiment. Compared with EMD denoising, EEMD denoising, CEEMD denoising, bilateral filtering denoising, and VMD denoising, the method can not only retain the useful signal better but also filter out more noise. In the following research, we will focus on the signal processing speed of this method and how to achieve better results in the signal processing process of different specifications of elevators.