An improved method for EMD modal aliasing effect

Abstract. In order to eliminate the modal aliasing effect of the empirical modal decomposition (EMD) caused by signal interruption (singular event, cycling wave, local impact signal) or the interaction of each component, two elimination methods are proposed, namely, high frequency harmonic method and mask signal method. Firstly, two causes of EMD mode aliasing are analyzed: one is the local oscillation in the signal; the other is the amplitude-frequency relation of the signal. Secondly, the local impact mode aliasing is solved by the high frequency harmonic method. Finally, mask signal method is used to eliminate modal aliasing of similar frequencies. The simulation results show that the two methods of eliminating EMD modal aliasing, which are collected high frequency harmonic method and mask signal method, can quickly and effectively curb the modal aliasing phenomenon on the basis of ensuring the adaptive characteristics of EMD, and provide a novel method to expand the real-time online application of EMD.


Introduction
With the increasing demand for mechanical on-line monitoring and intelligent diagnosis, higher requirements have been put forward for the self-adaptability, computing efficiency and fault pattern recognition efficiency of signal processing technology [1]. Empirical Modal Decomposition (EMD) is a novel signal adaptive decomposition technology, which can adaptively decompose signals into several Intrinsic Mode functions (IMF) according to its own scale, and is very suitable for processing nonlinear and non-stationary signals. EMD is an adaptive decomposition process that is performed according to the time characteristic scale characteristics of the signal itself. The result of the implementation is to separate the different frequency components in the signal layer by layer according to the frequency from high to low into corresponding IMF [2]. The IMF can be intuitive understanding as composed of a set of similar frequency narrowband signal, a correct the IMF should contain only the same scale characteristics of the signal component (note: scale characteristics is not a specific value, but a scale), as the size of different characteristics of the signal components in the same the IMF, creates a modal aliasing phenomenon. Different from the phenomenon of endpoint effect, modal aliasing may occur at any time of the signal between any adjacent IMF and have a great impact on the subsequent components [3]. If it cannot be effectively eliminated, the decomposition results will lose their physical significance. Therefore, modal aliasing is another obstacle to the application of EMD [4].
Huang firstly proposed a method of interrupt detection by setting extremum spacing parameters to control the screening of extremum points in each order of IMF screening. Although this method solves such problems very well, in practical application, this time interval is difficult to be determined before the standard EMD decomposition is implemented, which damages the self-adaptability of EMD [5]. Huang's group combined with Rilling's group on the decomposition properties of noise. A set of white noise method (Ensemble Empirical Modal Decomposition, EEMD) is proposed to suppress aliasing and maintain the adaptive characteristics of EMD. However, this method requires multiple EMD operations, which results in low efficiency and limits its online application. Through a large number of studies on decomposing conditions, Riiling provided a theoretical basis for solving the problem of modal aliasing caused by signal interaction in the future. The validity of two signals decomposed by EMD is discussed and the condition of effective separation of two signals is given [6]. Derring proposed the influential mask signal method for the first time, which realized the separation of two similar frequency components by adding a set of mask signals. Although the amplitude of the added signal depends on experience, it provides a new idea for solving modal aliasing [7].
The most significant characteristic of EMD is the self-adaptability of decomposition, which is also the main reason that EMD has been widely studied and applied. How to improve the operational efficiency of eliminating modal aliasing on the basis of the self-adaptive decomposition of EMD has become an important topic in the application research.

Aliasing caused by local oscillation (interrupt signal) in a signal
For non-stationary signals, time-scale parameters are characteristic parameters based on signal characteristics. Although they have no quantitative relationship with Fourier spectrum, they can also reflect signal characteristics. According to the selected reference points, there are generally two definitions of time scale: (a) The time span of two adjacent zero crossings; (b) the time span of two adjacent extreme points.
The above two definitions of time scales can reflect the local characteristics of signals changing with time, among which the description method based on adjacent extreme points is more effective, because even if the signal does not cross zero, the modal information of the signal can also be grasped through the span of extreme points.
When such as Eq. (1) (as shown in Fig. 1) signal of z(t), the regularity of the harmonic signal in the local high frequency oscillation, the dramatic change in the local extreme value point spacing, according to the extreme value point fitting envelope local greatly distorted (Fig. 2), and eventually in the decomposition results, the phenomenon of the same IMF features include multiple scale (

Aliasing caused by signal amplitude-frequency relationship
Take a randomly generated 1000 points discrete Gaussian white noise as an example to perform EMD decomposition, and the results are shown in Fig. 4.
As can be seen from the Fig. 4, the noise signal is decomposed into 9 components according to the frequency from high to low, and the extreme value points of each component are , = 1,…,·9. Although each IMF is a narrowband signal containing a frequency band, the distribution of the number of extreme points can reflect the frequency component characteristics of each IMF. Set = ⁄ , is always around 2, i.e., on a signal cycles the IMF contains about twice the next IMF. Huang inferred from this that EMD has similar binary filtering properties. VIBROENGINEERING PROCEDIA. NOVEMBER 2020, VOLUME 35  Two simulation signals ( ) and ( ), which are composed of four single-component signals of different frequencies, are constructed to illustrate. Where, the ratio between the frequencies of each component of signal ( ) is 1.25, 4, and 2 respectively; The ratio between the frequencies of each component of signal ( ) is 4, 1.3333 and 1.5, respectively. The result of EMD is shown in Fig. 5(a, b): 1 ( ) = sin(0.1 ) + 3sin(0.08 ) + 0.5sin(0.02 ) + cos(0.01 ), 2 ( ) = 3sin(0.08 ) + 0.5sin(0.02 ) + sin(0.015 ) + cos(0.01 ), = 500 , = 4. ( As can be seen from Fig. 5, it can be seen from the waveform of ( ) and ( ) that there is no interrupt in either signal. ( ) is decomposed into three results, in which and IMF2 can basically reflect the information of two components in the original signal, while IMF1 has obvious aliasing phenomenon, because ⁄ = 1.25 < 2. As a result, EMD cannot separate these two modes, which verifies the Eq. (5). In fact, IMF1 is the superposition of the first two high-frequency components in ( ), which is regarded by IMF as an IMF with certain modulation degree without separation. ( ) is decomposed into 5 components, among which IMF1 can perfectly reflect the first component of the signal, while the subsequent components have obvious aliasing, because ⁄ = 1.33 < 2 and ⁄ = 1.5 < 2, which do not meet the frequency conditions of signal separation.

Processing method of local impact mode aliasing
According to the characteristic that most of the local impact signals are high-frequency signals, a high frequency simple harmonic signal close to the highest frequency in the signal is added to the signal to provide the scale characteristics suitable for high frequency impact. The local high frequency and high frequency simple harmonic signals are decomposed into IMF1 to reduce the impact on the subsequent components. This algorithm does not need to perform standard EMD decomposition for the signal in advance, and has a good processing effect for the modal aliasing phenomenon caused by high frequency impact, and has a small amount of computation. In this paper, the high frequency signal frequency is set near the highest analysis frequency corresponding to the sampling frequency.
The results of EMD processing by ( ) in this method are shown in Fig. 6. It can be found that although there is a small amount of residue in the impact site, the decomposition accuracy of this method is significantly. Moreover, there is no extra false component except the superposed high-frequency component.
VIBROENGINEERING PROCEDIA. NOVEMBER 2020, VOLUME 35 The two components of aliasing are 1 ( ) = sin(0.1 ) and 2 ( ) = 3sin(0.08 ). According to the analysis in Section 2.2, when the frequency ratio of two equal-amplitude signals is between 0.5 and 2, the two signals will be considered as a pattern by EMD. According to this principle, the separation of ( ) and ( ) can be realized if a signal ( ) with appropriate frequency and amplitude can be constructed to form an aliasing mode after the superposition of ( ) and ( ). In this case, = 0.05 and = 0.04, so the lower frequency = 0.05/2 = 0.025 and the upper frequency = 0.05×2 = 1. Take the frequency of = 0.18, construct ( ) = 1.5sin(0.32 ), set: EMD is applied to ( ) and ( ), and the decomposition results, IMF + and IMF -, can be obtained respectively: Taking the actual hydraulic leakage signal as an example, the decomposition results are shown in Fig. 8. It can be seen from Fig. 8(a, b) that the signal is correctly decomposed to the mask signal of corresponding scale, and the decomposition result of the original signal is obtained after taking

Conclusions
This paper analyzes the form of EMD mode aliasing, causes, influences on EMD and solutions. The principle, characteristics, advantages and disadvantages of high frequency harmonic wave method and mask signal method are analyzed, and their respective advantages and disadvantages are verified by typical simulation signals. On the basis of the EMD adaptive characteristics, the model aliasing is quickly and effectively curbed, and the effectiveness of the method is verified by the hydraulic signal simulation, which provides a new method for the expansion of the real-time online application of EMD.