Published: 26 November 2020

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
Corresponding Author:
Kai Chai
Views 219
Reads 136
Downloads 887

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.

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)=basebandcos(2πt/200)t=2000, + interrupt0.2cos2πt/12t1180,1200.

Fig. 1Original signal containing local oscillation

Original signal containing  local oscillation

Fig. 2Local 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. 3Same modal component of the decomposition result contains multiple time scale components

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

Fig. 4Gaussian white noise EMD results

Gaussian white noise EMD results

Table 1Binary 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
z1t=sin0.1πt+3sin0.08πt+0.5sin0.02πt+cos0.01πt,z2t=3sin0.08πt+0.5sin0.02πt+sin0.015πt+cos0.01πt,
t=500s, fs=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. 5Modal 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. 6EMD 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
z1t=sin0.1πt+3sin0.08πt+0.5sin0.02πt+cos0.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=xt+st,x-t=xt-st.

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

5
IMF=(IMF++IMF-)2.

Fig. 7Same 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. 8EMD 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

  • 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.
  • Xu Zheng Guang An alternative envelope approach for empirical mode decomposition. Digital Signal Processing, Vol. 20, Issue 1, 2010, p. 77-84.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.

Cited by

Decomposed intrinsic mode functions and deep learning algorithms for water quality index forecasting
Kok Poh Wai | Chai Hoon Koo | Yuk Feng Huang | Woon Chan Chong
(2024)
Multiple solutions and orbit change in energy harvesting system with a flag configuration
Jerzy Margielewicz | Damian Gąska | Sławomir Bucki | Grzegorz Litak | Sreeja Sadasivan
(2024)
Proceedings of the 2023 4th International Conference on Computer Science and Management Technology
Haojian Huang | Wenkai Yin | Chuanyu Qin | Weizhen Peng | Ziyun Huang | Gang Liu | Hongwei Mo
(2023)
Drought prediction based on an improved VMD-OS-QR-ELM model
Yang Liu | Li Hu Wang | Li Bo Yang | Xue Mei Liu | Zaher Mundher Yaseen
(2022)
Research on Dynamic Measurement Method of Flow Rate in Tea Processing
Zhangfeng Zhao | Gaohong Liu | Yueliang Wang | Jiyu Peng | Xin Qiao | Jiang Zhong
(2022)

About this article

Received
21 October 2020
Accepted
16 November 2020
Published
26 November 2020
SUBJECTS
Mathematical models in engineering
Keywords
empirical modal decomposition
modal aliasing
high frequency harmonic method
mask signal method