An Adaptive Noise Reduction Method Based on Improved Dislocation Superposition Method for Abnormal Noise Fault Component of Automotive Engine

&e acoustic signal generated by mechanical motion contains the information of its motion state, but when the signal-to-noise ratio (SNR) is low, the accuracy of real-timemonitoring mechanical motion state by the acoustic signal is low.&is study proposes an adaptive noise reduction method based on the dislocation superposition method (DSM), which can realize the adaptive noise reduction and the extraction of fault a component from the automobile engine abnormal noise signal of low SNR. Firstly, the wavelet coefficients of engine abnormal noise signal are obtained by continuous wavelet transform (CWT), and the fault feature points of the abnormal noise signal in each period are extracted by setting hard threshold function, window function, and feature points extraction algorithm. &en, the signal segments containing fault components are obtained by using the position of feature points to extend the length of the fault component forward and backward, respectively, and Pearson’s correlation is calculated by traversal to determine the starting superposition point of each signal segment containing fault components. Finally, the signal segments of the odd group and even group are selected for superposition calculation. When the superposition stop condition is not satisfied, the number of superpositions increased until the stop condition is satisfied, and the superposition signal can be used as a fault component.&e experimental results show that, compared with the improved DSM, this method has a good effect on the noise reduction and extraction of fault components of automobile engine cylinder knocking fault, and the effectiveness of this method is verified. &is method is used to reduce the noise and extract the fault components of automobile engine cylinder missing fault and knock fault, and good results are obtained.


Introduction
At present, cars play an important role in human life and are a necessary condition for human travel. e engine is an important part of the automobile, and its structure is complex [1]. In the working process of the engine, due to wear, fatigue, aging, and other factors, the engine failure may cause serious economic losses and even casualties [2,3]. When the engine failure produces an abnormal noise signal, the higher SNR of the abnormal noise signal will make engine fault diagnosis easier and more accurate, so it is necessary to improve the SNR of the engine abnormal noise signal [4]. e traditional noise reduction methods mainly include Wiener filtering, spectral subtraction, and minimum mean square deviation, but the effectiveness of noise reduction for the nonstationary signal and short-term transient signal is significantly reduced [5,6]. e commonly used noise reduction methods mainly include wavelet transform (WT) [7], empirical mode decomposition (EMD) [8], local mean decomposition (LMD) [9], DSM [10], etc.
As a widely used signal processing tool, WT has strong multiresolution analysis ability in the time domain and frequency domain. e wavelet denoising generally uses setting the wavelet coefficient threshold to eliminate the noise information and then uses the inverse WT to reconstruct the signal from the threshold wavelet coefficient [11]. Li [12] used the wavelet decomposition and reconstruction algorithm to decompose and reduce noise and reconstruct and analyze the spectrum of bearing vibration signal. Experiment shows that the wavelet noise reduction method is very suitable for fault frequency detection of weak vibration signal of the rolling bearing in low SNR cases. Moshrefi et al. [13] used WT and adaptive filtering to denoise the knock signal of the internal combustion engine. is method is applied to real knock signal, which showed superiority compare to previous works and led to a 13.2% improvement in detection accuracy. Kai et al. [14] proposed a wavelet denoising method based on improved threshold function by studying soft threshold function and hard threshold function. By this method, the feature of fault signal is highlighted and the fault diagnosis effect was improved. Wen and You [15] used the wavelet decomposition and reconstruction algorithm to denoise the fault signal of highspeed rolling bearing, and the experimental data showed that most of the noise has been filtered out compared with the original data. Although wavelet denoising has good processing ability for nonstationary signals, the selection of wavelet basis function is very difficult. Different wavelet basis functions have different denoising effects. To solve this problem, Huang et al. [16] proposed the EMD, which can efficiently decompose nonlinear and nonstationary signals without any set of basis functions. Sun et al. [17] used the wavelet threshold noise reduction algorithm based on EMD to solve the problem of complex centrifugal pump vibration signals. e experimental results show that the wavelet soft threshold denoising algorithm based on EMD decomposition has a better noise reduction effect when the centrifugal pump vibration signal is used as the noise reduction object. Ren and Liu [18] proposed an adaptive reduction noise and feature extraction algorithm based on improved EMD and verified the effectiveness and feasibility of the method by simulation signals and examples. Although EMD has been successfully applied in the field of signal noise reduction, it also has its limitations, such as the frequent occurrence of modal mixing. In view of the shortcomings of EMD, Wu and Huang [19] proposed ensemble empirical mode decomposition (EEMD) in 2009. EEMD adds finite amplitude Gauss white noise to the signal, and then the signal with the white noise is decomposed as a whole, to effectively suppress the mode mixing of EMD. However, this method can only reduce modal mixing to some extent.
In recent years, scholars have also proposed some other signal noise reduction methods. Gao et al. [20] proposed a new method based on LMD and wavelet denoising to analyze the signals of bearing outer ring, inner ring, and ball. However, LMD has the phenomenon of modal mixing, which reduces the accuracy of signal decomposition and affects the accuracy of the noise reduction signal. Dayong et al. [21] proposed DSM based on random decrement technique in 2015. Compared with other methods, DSM only calculates in time domain. In the superposition process, the method does not destroy the correlation components in the mixed signal and effectively avoids the modal mixing. In 2019, Dayong et al. [10] used improved DSM to automatically extract engine fault components. Although the improved DSM can automatically extract the fault components of the automotive engine quasiperiodic signal, this method needs to use the pulse number of the encoder to determine the starting superposition point of each quasiperiodic fault signal. In practical applications, automobile engines are not suitable for installing encoders, and encoders are prone to failure when subjected to severe impact, so the improved DSM is limited in some practical applications.
To overcome the above problems, this paper proposes an adaptive noise reduction and extraction method of engine abnormal noise signal fault components based on improved DSM. e essence of this method is the superposition calculation in the time domain, which avoids the modal mixing and can better deal with the acoustic signal with low SNR. e method can adaptively select the starting superposition points, superposition length, and superposition number. By changing the superposition number, the noise reduction degree of fault components can be changed. Compared with the improved DSM, it not only reduces the use of encoders and improves the practicability of DSM but also makes it more convenient to extract fault components.

DSM Review
e mathematical expression of the traditional DSM is as follows: where S (n) is the original signal,Ŝ (n) is the signal processed by DSM (named superposition signal), K is the number of superpositions (K � 0, 1, 2, . . .), and L is the superposition step length (the period of the signal to be processed). Figure 1 shows the graphical description of DSM processing results.
In Figure 1, S is the target signal; N is the interference signal; SN is a mixed signal of S and N; SN 1 , SN 2 , and SN 3 are superposition signals obtained by superposition of 5, 15, and 20 times, respectively, using equation (1); e superposition step length L is the period of signal S. Compared with signal SN, the component proportion of target signal in signals SN 1 , SN 2 , and SN 3 increases with the increase of superposition times. On the contrary, the proportion of interference signal decreases with the increase of superposition times [22]. Generally, Pearson's correlation coefficient is used to compare the similarity between the target and the superimposed signals to test the DSM processing effect [10].

Adaptive Noise Reduction and Extraction
Method of Engine Abnormal Noise Fault Component 3.1. Influence Factors of DSM. Due to the system error of the engine, the actual speed of the engine is slightly changed, which causes the acoustic signal of the engine to be a quasiperiodic signal. Figure 2 shows the schematic of the quasi-periodic signal.
When the engine has an abnormal noise impact failure, the fault signal is quasi-periodic, resulting in different superposition step lengths. In the dislocation superposition, it may cause the phase deviation of the impact fault component and even lead to the destruction of the fault component. According to the characteristics of impact failure, it is artificially divided into the "intense change region" and the "stable region," as shown in Figure 3. When the engine impact fault acoustic signal is noise reduction by DSM, due to the less fault energy and information contained in the stable region, to improve the computational efficiency, the intense change region containing more energy and information is separated as the impact fault component [10,23]. e noise reduction degree of the fault component also depends on the number of superpositions.
eoretically, the more the number of superpositions is, the better the noise reduction effect is. However, the more the number of superpositions is, the longer the consumed time is. When the number of superpositions reaches a certain amount, the noise reduction effect shows a stable trend. erefore, the adaptive noise reduction process of engine impact fault acoustic signal is a process of automatically finding the starting superposition points of impact fault acoustic components, the length of the intense change region, and the number of superpositions.

Adaptive Noise Reduction and Extraction Method
Based on Improved DSM and CWT. Based on CWT and improved DSM, an adaptive noise reduction method for automobile engine fault acoustic signal is proposed. Figure 4 shows the flowchart of the method. e details are as follows.

Extraction of Feature Point Location of Each Period Impact Fault Component.
e "intense change region" has the characteristics of transient, periodic, and large energy which is not easily submerged by background noise [23]. e mother wavelet which is similar to the waveform of the impact fault component is selected for CWT to obtain wavelet coefficients, and the wavelet coefficients are hard thresholding processes. Finally, the window function and feature points extraction algorithm are used to extract the  (1) Load the original fault signal S (n), and select the mother wavelet similar to the intense change region of the fault signal to perform CWT on the signal S (n) to obtain the wavelet coefficients WT S . (2) Eliminate the smaller wavelet coefficients containing interference components in the fault signal through the hard threshold function and retain the larger wavelet coefficients containing the impact fault components. e mathematical model of the hard threshold function is where η H (WT S , λ) is the processed wavelet coefficient by hard threshold function, denoted as η H ; WT S are wavelet coefficients; and λ is the size of the set hard threshold. Two diagonally paired time-frequency points (t 1 , f 1 ), (t 2 , f 2 ) are selected to form a rectangular window R, and a timefrequency block is selected from the frame of the thresholding wavelet time-frequency diagram by using the rectangular window R, which is denoted as η H (a 0 , b 0 ), where Select W s range (n s ) as the reference signal Intercept W j fra (n j ) from W j erange (n j ), traverse and calculate the correlation between W s range (n s ) and W j fra (n j + τ j ) to get the optimal offset τ j opt Select number 2K of W j fra (n j + τ j opt ), and divide W j fra (n j + τ j opt ) into odd group and even group to perform superposition operations to obtain D W1K and D W2K Ifthe selected signal length is not an integer multiple of the rectangular windowlength, the maximum value of i is where [·] is rounding function and n tot is the total number of sampling points of the signal S (n). Let the wavelet coefficient corresponding to the point C i (t spe i , f spe i ) be η HC i , and two condition algorithms are set for η HC i to extract feature points.

Extraction of Impact Fault Component by
where n jj is the number of sampling points corresponding to the fault feature point, t jj is the time of the feature point of the impact fault component, and f s is the sampling frequency. e range of the impact fault component of each period is as follows: (σ −1) sampling points are extended forward and backward from the n jj th sampling point, and the signal segment is denoted as W ; σ is the length of the intense change region. Due to the interference of background noise, the feature point location of per period is different in the intense change region. erefore, the intense change regions in W range i (n j ) may have location offset of different degrees in different periods. As shown in Figure 5, there is offset τ of different degrees in the intense change regions of different periodic W range i (n j ), which will affect the dislocation superposition effect and lead to inaccurate superposition results.
To solve the problem of the offset of the intense change region, the specific method is as follows.
A segment of W e continuous intercept length of (2σ − 1) sampling points from the first sampling point in W , and record the offset signal of W fra i (n j ) as W fra i (n j + τ j ), as shown in Figure 6. e subset of τ j is denoted as L, and L is set to {0: σ − 1 + b}. Find the optimal offset τ opt j of the offset τ j in W fra j (n j + τ j ); that is, also find the optimal starting superposition point, where the superposition length is (2σ − 1) sampling points.
Traverse all the τ j values in L, and calculate Pearson's correlation coefficient ρ sj (W range s (n s ), W fra i (n j + τ j )) of the reference signal segments W range s (n s ) and W fra i (n j + τ j ), respectively; the offset corresponding to the maximum value of ρ sj (W range s (n s ), W fra i (n j + τ j )) is the optimal offset τ opt i . e equation of the proposed method is as follows: where τ opt j is the best offset; argmax [·] is arguments of the maxima; W range s (n s ) is the reference signal segment; W fra i (n j + τ j ) is the offset signal segment; and ρ sj (W range s (n s ), W fra i (n j + τ j )) is to calculate the correlation between W range s (n s ) and W fra i (n j + τ j ). 2K consecutive signals W fra i (n j + τ opt i ) are selected, which are divided into two groups according to the parity of the sequence number j, and each group is superimposed with different times according to the sequence number from low to high. e equation of the proposed method is as follows:

Shock and Vibration
where D

Parameters Setting
(1) Hard reshold λ Setting. e hard threshold λ is mainly to remove the smaller wavelet coefficients that only contain background noise and retains the larger wavelet coefficients that contain fault components. When the value of λ is too small, the smaller wavelet coefficients generated by background noise may be incompletely removed, resulting in unobvious fault periodicity. When the value of λ is appropriate, the smaller wavelet coefficients are completely removed, only the larger wavelet coefficients containing fault components are retained, and the fault periodicity is obvious. When the value of λ is too large, not only the smaller wavelet coefficients are completely removed, but also part of the larger wavelet coefficients containing fault components will be removed, resulting in a jump in the period of the fault signal. However, the location of the fault component can be found without affecting the superimposing effect. erefore, the hard threshold λ should be as large as possible.
(2) Rectangular Window R Setting. e rectangular window R is composed of two diagonally paired time-frequency points (t 1 , f 1 ), (t 2 , f 2 ); then, the length of the rectangular window R is l, l � t 2 − t 1 and the width w � f 2 − f 1 . Considering the frequency distribution of the impact fault signal and the length σ of the impact fault component, the time of the intense change region is e rectangular window starts from t � 0 s, so t 1 � 0 s. According to the characteristic frequency of the impact fault signal, which is generally distributed in the middle and low frequency, f 1 can be set between 0 Hz and 500 Hz, and f 2 ≤ 5000 Hz. If t 2 is too large, it will cause the rectangular window to select multiple fault signal periods, which is not conducive to extracting the fault signal period. If t 2 is too small, it will increase the amount of calculation. erefore, the general value range of l is 2t σ ≤ l ≤ 4t σ , which corresponds to 2σ to 4σ sampling points.
where f s is the sampling frequency of the acoustic signal.    Figure 6: W fra i (n j + τ j ) signal segment schematic diagram. e part selected by the blue frame is W fra i (n j ), and the part selected by the red frame is W fra i (n j + τ j ).
increase. When ρ generally increases to 0.9, it will tend to be stable, and there is an unobvious increase trend [10]. To achieve a better noise reduction effect for the fault component, M should be set to a value between 0.75 and 0.85, where M = 0.80. Figure 7 is the physical diagram of the automobile engine fault detection test bench, which consists of an automobile engine, sound sensor, data acquisition card, encoder, and computer. e engine model is EA211 with the detailed parameters shown in Table 1. e sound sensor with the frequency range of 20 Hz to 20 kHz is placed above the cylinder to receive the abnormal sound signal generated by the engine. e data acquisition card uses the USB-6341 produced by National Instruments Company. When collecting the sound signal, the sampling frequency is set to 44100 Hz. e encoder model is ZSP3806GC, and the resolution is 100P/R. e encoder can synchronously revolve with the crankshaft by fixing it on the front end of the crankshaft through a coupling.

Fault Component Extraction of Engine Knocking
Cylinder. Figure 8(a) shows the knocking cylinder abnormal noise signal and encoder pulse signal collected under the condition of engine running at 1800 r/min, Figure 8(b) uses db5 wavelet as the mother wavelet to perform CWT on the knocking cylinder abnormal noise signal to obtain the wavelet time-frequency diagram. e two diagonally paired time-frequency points of the rectangular window R are (0, 0), (0.0227, 4000), respectively. e hard threshold of the wavelet coefficient is λ � 0.94, and the length of the intense change region is calculated by equation (8) to obtain σ � 300, b � 299. e subset L of offset τ j is set to {0 : 598}. Table 2 is the sampling point corresponding to the fault feature point n jj of each periodic obtained by data processing using the above method. Table 3 shows the number of sampling points for each period of impact fault signal calculated by using the feature points of each periodic failure component. It can be seen from Table 3 that the number of sampling points in the 6 th -8 th periods is 8708, which is quite different from the number of sampling points in other periods.
is is due to the excessive setting of the hard threshold λ, which leads to the zeroing of the wavelet coefficients corresponding to the feature points of the fault signal components in the 7 th and 8 th periods after the hard threshold processing, so the number of sampling points in the 6 th , 7 th , and 8 th periods cannot be obtained, but it does not affect the superposition effect.
Selecting s � 1, 299 sampling points are extended forward and backward from the first feature point 3530 of the fault component to obtain the reference signal segment W range 1 (n 1 ), as shown in Figure 9. Table 4 shows the optimal offset τ opt j of each W fra j (n j + τ j ) by equation (5), and W fra j (n j + τ opt j ) is obtained according to τ opt j . e continuous 2K signal segments W fra j (n j + τ opt j ) including fault components are selected and divided into two groups according to the parity of the sequence number j for superposition operation, and the correlation coefficient ρ opt (K) of the parity group is calculated. Figure 10 shows the variation curve of ρ opt (K) and the number of superpositions.
When K � 14, ρ opt (14) � 0.8030, which is greater than the threshold M, stopping superimposing. Figure 11 shows the fault component signals D opt W114 and D opt W214 after noise reduction. Figure 12 shows the time-domain diagram of the odd group and even group superimposed signals obtained by superimposing 10 times using the improved DSM. e optimal starting point of superposition is the sampling point of the acoustic signal corresponding to the 28 th pulse of each cycle encoder signal. e superposition length is set to 600 sampling points, the odd group and even group are superimposed 10 times to obtain D opt W110 and D opt W210 , and the correlation reaches 0.8020.
According to Figures 11 and 12, the length σ � 150 of fault component is accurately extracted. Figure 13 shows the fault components W 1 , W 2 , W 3 , and W 4 extracted from D opt W114 , D opt W214 , D opt W110 , and D opt W210 , respectively. e correlation between W 1 and W 3 is 0.9029, and the correlation between W 2 and W 4 is 0.9290, which further verifies that the fault component noise reduction method achieves the same effect compared with the improved DSM. Figure 14 shows the physical picture of the first cylinder lacking of the automobile engine by pulling out the cylinder line of the first cylinder. Figure 15 shows that the engine produces abnormal noise signal for the first cylinder lacking. e db5 wavelet is selected as the mother wavelet to perform CWT on the abnormal noise signal of the lack of cylinder. e time-frequency points of the two diagonally paired of the rectangular window R are set as (0, 500), (0.0227, 4000), respectively. e hard threshold of the wavelet coefficient λ = 1, σ = 300, b = 299, and the subset L of the offset τ j is set as {0 : 598}.

Fault Component Extraction of Engine Lacking Cylinder.
14 15. Table 5 shows the sampling points corresponding to the feature points of each period of fault component obtained by the above method. Selecting s � 1, the reference signal W range 1 (n 1 ) is obtained by extending 299 sampling points forward and backward from the 1794 th sampling point, as shown in Figure 16. Table 6 shows the optimal offset τ opt j of each W fra j (n j + τ j ) signal by equation (5), and W fra j (n j + τ opt j ) is obtained according to τ opt j . Figure 17 shows the variation curve of the correlation coefficient ρ opt (K) of the superposition signals of the odd group and even group with the number of superpositions K. When K � 6, ρ opt (6) � 0.8048, which is greater than the threshold M and stopping superimposition, D opt W16 and D opt W26 of the fault component after noise reduction are got, as shown in Figure 18.

Shock and Vibration 7
(1) (3)               Figure 19 shows the engine knock signal due to the low gasoline label. e db5 wavelet is selected as the mother wavelet to perform CWT on the knocking signal. e timefrequency points of the two diagonally paired of the rectangular window R are set as (0, 0), (0.0227, 4000), respectively. e hard threshold of the wavelet coefficient λ � 2.8, σ � 300, b � 299, and the subset L of the offset τ j is set as {0 : 598}. Table 7 shows the sampling points corresponding to the fault feature points. Due to the large value of hard threshold λ of wavelet coefficient, the feature points of fault components do not show periodic feature, but the interval of sampling points between each two feature points is approximately 4000 times. Table 8 shows the optimal offset τ opt j of each W fra j (n j + τ j ), and W fra j (n j + τ opt j ) is obtained according to τ opt j . Figure 20 shows the variation curve of the    and D opt W28 of the fault components after noise reduction is got, as shown in Figure 21.

Experimental Results Analysis
In order to further verify the effectiveness of this method, the experimental results of this method and the improved DSM are compared by using the knocking cylinder experiment of automobile engine, and then the effectiveness of this method is verified. In this method, odd group signals and even group of engine knocking cylinder signals are superimposed 14 times, and the correlation between the noise reduction signals D opt W114 and D opt W214 is 0.8030. en, the noise reduction processing and extraction of engine knocking cylinder failure component are carried out by using the improved DSM. When odd group signals and even group signals are superposed 10 times, the correlation between the noise reduction signals D , where the correlation between W 1 and W 3 is 0.9029, and the correlation between W 2 and W 4 is 0.9290. e results show that this method can further replace the use of encoders and achieve the same noise reduction effect as the improved DSM.
In the engine cylinder lacking experiment, this method is used to reduce the noise of fault components. e correlation between the noise reduction signal D opt W16 and D opt W26 obtained by superposition of odd group signals and even group signals for 6 times is 0.8048. In the engine knock experiment, due to the relatively small setting of hard threshold λ, some fault feature points are not selected accurately, resulting in the correlation of superposition signals of odd and even groups after superposition for 2 and 3 times gradually decreasing, as shown in Figure 20. However, with the gradual increase in the number of superpositions, the correlation between D opt W18 and D opt W28 obtained by superposition of odd group and even group for 8 times is 0.8091.
From the analysis results, the method can replace the improved DSM, effectively reduce the noise, and extract the fault component of the engine abnormal noise signal. By increasing the threshold M and superimposition number K, the accuracy of the extracted fault components can be improved. e method can also use the extracted fault components to establish a database of automobile engine faults, which lays a foundation for the later diagnosis and classification of automobile engine faults.

Conclusion
Althoughthe improved DSM can automatically find the starting superposition point andsuperposition length of the periodic signal, it needs to use the encoder toassist in finding the starting superposition point. In practical applications, many occasions are not suitable for the use of encoders, which limits the applicability of DSM. erefore, based on the improved DSM, this paper proposes an adaptive noise reduction method of automobile engine abnormal noise signal, which cancels the use of the encoder. Using this method, the fault components of knocking cylinder, cylinder lacking, and knocking are extracted from the abnormal sound signal of automobile engine, and the effectiveness of the method is verified. e method can adaptively obtain the starting superposition point and superposition length from the engine abnormal sound signal and determine the number of superpositions, which replaces the encoder in the improved DSM to determine the starting position of the automobile engine impact fault and greatly improves the applicability of DSM. However, when the hard threshold λ is set to a small value, the selection of fault feature points will be inaccurate, and then the selected fault components will be biased, resulting in inaccurate noise reduction signals. e specific problem of hard threshold setting will be further studied. When the starting superposition point is selected, the accuracy is reduced compared with the improved DSM, resulting in more superposition times than the improved DSM. However, according to the superposition results, the increase of superposition times is within an acceptable range.
In the future, the hard threshold λ setting problem will be further improved, and the corresponding fault component database will be established.
e fault type will be determined by comparing the extracted fault component with the fault component in the database.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this study.