An improved method for EMD modal aliasing effect

Shao-Wei Feng1 , Kai Chai2

1, 2College of Naval Architecture and Ocean, Naval University of Engineering, Wuhan, 430033, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 35, 2020, p. 76-81. https://doi.org/10.21595/vp.2020.21778
Received 21 October 2020; received in revised form 9 November 2020; accepted 16 November 2020; published 26 November 2020

Copyright © 2020 Shao-Wei Feng, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF References
Cite this article
Views 35
Reads 18
Downloads 206
CrossRef Citations 0
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.

Keywords: empirical modal decomposition, modal aliasing, high frequency harmonic method, mask signal method.

1. 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.

2. EMD mode aliasing problem

2.1. 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 (Fig. 3):

(1)
z ( t ) = b a s e b a n d c o s ( 2 π t / 200 ) t = 2000 ,   +   i n t e r r u p t 0.2 c o s 2 π t / 12 t 1180,1200 .

Fig. 1. Original signal containing local oscillation

 Original signal containing  local oscillation

Fig. 2. Local distortion occurs in the upper and lower envelope

 Local distortion occurs  in the upper and lower envelope

2.2. 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 Pi, i= 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 λi=Pi/Pi+1, λi 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.

Fig. 3. Same modal component of the decomposition result contains multiple time scale components

 Same modal component of the decomposition result contains multiple time scale components

Fig. 4. Gaussian white noise EMD results

 Gaussian white noise EMD results

Table 1. Binary decomposition characteristics of EMD

1
2
3
4
5
6
7
8
9
Extremum points pi
714
350
178
92
44
22
12
6
4
λ i
2.040
2.040
1.934
2.090
2
1.833
2
1.500

Two simulation signals z1t and z2t, 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 z1t is 1.25, 4, and 2 respectively; The ratio between the frequencies of each component of signal z1t is 4, 1.3333 and 1.5, respectively. The result of EMD is shown in Fig. 5(a, b):

(2)
z 1 t = s i n 0.1 π t + 3 s i n 0.08 π t + 0.5 s i n 0.02 π t + c o s 0.01 π t , z 2 t = 3 s i n 0.08 π t + 0.5 s i n 0.02 π t + s i n 0.015 π t + c o s 0.01 π t ,        
t = 500 s ,           f s = 4 .

As can be seen from Fig. 5, it can be seen from the waveform of z1t and z2t that there is no interrupt in either signal. z1t is decomposed into three results, in which r2 and IMF2 can basically reflect the information of two components in the original signal, while IMF1 has obvious aliasing phenomenon, because f1/f2= 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 z1t, which is regarded by IMF as an IMF with certain modulation degree without separation. z2t 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 f2/f3= 1.33 < 2 and f3/f4= 1.5 < 2, which do not meet the frequency conditions of signal separation.

Fig. 5. Modal aliasing due to close frequencies

 Modal aliasing due to close frequencies

a)z1t

 Modal aliasing due to close frequencies

b)z2t

3. 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 zt 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.

Fig. 6. EMD results of high-frequency harmonic processing with local shock signal

 EMD results of high-frequency harmonic processing with local shock signal

4. Processing method of modal aliasing at similar frequencies

Take the aforementioned simulation signal z1t as an example:

(3)
z 1 t = s i n 0.1 π t + 3 s i n 0.08 π t + 0.5 s i n 0.02 π t + c o s 0.01 π t .

The two components of aliasing are x1(t)=sin(0.1πt) and x2(t)=3sin(0.08πt). 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 x1t and x2t can be realized if a signal st with appropriate frequency and amplitude can be constructed to form an aliasing mode after the superposition of x1t and St. In this case, f1= 0.05 and f2= 0.04, so the lower frequency fL= 0.05/2 = 0.025 and the upper frequency fL=0.05×2 = 1. Take the frequency of f¯=0.18, construct s(t)=1.5sin(0.32πt), set:

(4)
x + t = x t + s t , x - t = x t - s t .

EMD is applied to x+t and x-t, and the decomposition results, IMF+ and IMF-, can be obtained respectively:

(5)
I M F = ( I M F + + I M F - ) 2 .

Fig. 7. Same modal component of the decomposition result contains multiple time scale components

 Same modal component of the decomposition result contains multiple time scale components

a) IMF+

 Same modal component of the decomposition result contains multiple time scale components

b) IMF-

 Same modal component of the decomposition result contains multiple time scale components

c) 0.5×(IMF++IMF-)

It can be seen from Fig. 7 that in the decomposition results of x+t and x-t, x1t and auxiliary signals are formed into alias mode, and IMF2 is successfully separated. After averaging x+t and x-t, the mask component is eliminated, and the true component of xt is obtained (deformation at the end point is caused by the endpoint effect).

Fig. 8. EMD results of modal aliasing elimination method for practical hydraulic fault signals

 EMD results of modal aliasing elimination method for practical hydraulic fault signals

a) IMF+

 EMD results of modal aliasing elimination method for practical hydraulic fault signals

b) IMF-

 EMD results of modal aliasing elimination method for practical hydraulic fault signals

c) 0.5×(IMF++IMF-)

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 the mean value. Direct observation shows that there is no obvious aliasing in the decomposition results of this method. The reconstruction error can reach a level 10-15.

5. 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.

References

  1. Colominas M. A., Schloti’Hauer G., Torres M. E. Improved complete ensemble EMD: a suitable tool for biomedical signal processing. Biomedical Signal Processing and Control, Vol. 14, Issue 1, 2014, p. 19-29. [Publisher]
  2. Xu Zheng Guang An alternative envelope approach for empirical mode decomposition. Digital Signal Processing, Vol. 20, Issue 1, 2010, p. 77-84. [Publisher]
  3. Hassanpour H., Zehtabian A. Time domain signal enhancement based on an optimized singular vector denoising algorithm. Digital Signal Processing, Vol. 22, Issue 5, 2012, p. 786-794. [Publisher]
  4. He Long, Yang Lihua, Huang Daren The study of the intermittency test filtering character of Hilbert-Huang transform. Mathematics and Computers in Simulation, Vol. 70, 2005, p. 22-32. [Publisher]
  5. Wu Z. H., Huang N. E. Ensemble Empirical Mode Decomposition: a Noise Assisted Data Analysis Method. Calverton Center for Ocean-Land-Atmosphere Studies, 2005, p. 855-895. [Search CrossRef]
  6. Rilling Gabriel, Flandrin Patrick One or two frequencies? The empirical mode decomposition answers. IEEE Transactions on Signal Processing, Vol. 56, Issue 1, 2008, p. 85-96. [Publisher]
  7. Ryan Deering, Kaiser James F. The use of a masking signal to improve empirical mode decomposition. IEEE International Conference on Acoustics, Speech, and Signal Processing, 2005. [Search CrossRef]