Academic Journal of Science and Technology ISSN: 2771-3032 | Vol. 1, No. 3, 2022 36 Estimation of EMD Envelope Based on MQ Interpolation Liming Liang1, a 1College of Mathematics and Physics, Geo-mathematics Key Laboratory of Sichuan Province, Chengdu University of Technology, Chengdu, 610059, China a1163937352@qq.com Abstract: Empirical modal decomposition (EMD) is a data-driven signal analysis method that can effectively decompose signals into components of different scales. However, there are still key issues such as endpoint effects, envelope fitting and theoretical foundation improvement that need to be solved. In this paper, we introduce Multi-Quadratic (MQ) fitting interpolation to estimate the envelope in the empirical modal decomposition, which is an efficient, high-precision and conformal fitting method without solving the system of equations. Keywords: EMD, Envelope fitting, MQ. 1. Introduction Since the EMD decomposition method was proposed by Huang et al. in 1998, it has been widely used for nonlinear signal processing in various disciplines. Although EMD has been developed for many years, there are still some problems. To solve the endpoint effect, Rilling et al. [1] proposed the mirror method, which takes the first extreme point at the ends of the signal as the boundary, maps the signal outward, and obtains an envelope curve covering the original signal by obtaining the mirror image of the original signal; Datig et al.[2] proposed the slope method, which obtained the estimate of the extension point by linear prediction of the two extreme points at the endpoints; Wu et al.[3] improved it based on Datig et al. improved the method based on Datig et al. by finding the slope from the adjacent extreme points of the signal endpoints and using the first extreme point of the endpoints and the slope to determine the linear equation where the extension point is located to determine the location of the extension point (Improved Slope Based Method (ISBM)); Xu et al.[4] proposed the Cubic Spline Based Method (CSBM), which makes use of the fact that the envelope after extension The local extrema at both ends of the signal are related to the extension points to obtain the extension points that satisfy the conditions. Regarding the modal mixing problem, many scholars have proposed improved EMD [5] and different envelope fitting algorithms [6]. In terms of screening conditions, the changes in IMF components, envelopes, or extreme value points before and after iteration [7] are mainly used to determine whether it is an eigenmodular function. In conclusion, EMD needs to be further investigated. To further solve the envelope fitting problem, a variable parameter MQ fitting interpolation algorithm is proposed. 2. Three Classical MQ Proposed Interpolation Operators The MQ function is a radial basis function. Franke compared 29 discrete point interpolation methods in terms of accuracy and stability, and found that the MQ function has the best overall performance in interpolation and is conformal. Although MQ interpolation shows good performance in scattered data interpolation, it still has some shortcomings. MQ interpolation is essentially a system of linear equations, which is very computationally intensive when the size of the scattered data points is large, and the coefficient matrix is also prone to pathology, resulting in large errors in the calculation results. Since MQ proposed interpolation not only has good approximation accuracy, but also does not need to solve the linear system of equations, so scholars have gradually started to study MQ proposed interpolation. Definition 2.1 Given a set [ , ]a b of nodes on an interval 0 1 na x x x b    � , 1 1 max( )j j j n h x x     , the MQ proposed interpolation satisfy the form 0 ( ) ( ) ( ), n j j Lf x f x x a x b     Beatson and Powell defined ( )AL f x , ( )BL f x and ( )DL f x three MQ proposed interpolation operators. The proposed ( )BL f x interpolation operator is given first, 1 0 0 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) n B j j n n j L f x f x x f x x f x x        1 0 0 1 0 1 ( ) ( ) ( ) 2 2( ) x x x x x       1 1 1 1 ( ) ( ) ( ) ( ) ( ) 2( ) 2( ) j j j j j j j j j x x x x x x x x x                1 1 1 ( ) ( ) ( ) . 2 2( ) n n n n n x x x x x         Definition of the proposed ( )AL f x interpolation operator 0 0 0 ( ) ( ) ( ) ( ) ( ) ( ) 1 2 n n A B n f x f x x x L f x L f x x x          1 0 0 0 1 0 0 ( ) ( ) ( ) ( ) ( ) + 2( ) 2( ) n n x x x x x x x x x          37 0 1 0 1 ( ) ( ) ( ) ( ) ( ) . 2( ) 2( ) n n n n n n n x x x x x x x x x             Wu et al.proposed a new MQ proposed interpolation operator 2 0 0 1 1 2 1 1 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) n D j j j n n n n L f x f x x f x x f x x f x x f x x                ， 1 0 0 1 0 1 ( ) ( ) ( ) + 2 2( ) x x x x x x      2 1 1 0 1 2 1 1 0 ( ) ( ) ( ) ( ) ( ) 2( ) 2( ) x x x x x x x x x x           1 1 2 1 1 1 2 ( ) ( ) ( ) ( ) ( ) 2( ) 2( ) n n n n n n n n n x x x x x x x x x x                  1 1 1 ( ) ( ) ( ) + . 2 2( ) n n n n n x x x x x x        3. Improved MQ Fitted Interpolation Fitting Envelope Zhang Jihong[8] first introduced MQ proposed interpolation into EMD decomposition for analyzing the similarity of DNA sequences between species, but did not analyze and discuss the effect of EMD decomposition, and chose fixed values of shape parameters. EMD decomposition is an iterative process, and the local extrema will change greatly, so it is obviously undesirable to use fixed values of shape parameters directly. Therefore, this paper proposes a new way of selecting shape parameters to make the envelope in EMD decomposition more adaptable. The improved shape parameter is defined as follows: 1 30.1c h , 1 1 max( ).j j j n h x x     Compared with the traditional interpolation envelope method, approximating the envelope by MQ proposed interpolation has the following advantages: (1) The function expression is a display expression, which is easy to implement programmatically without solving a linear system of equations and is faster to compute; (2) The MQ proposed interpolation function has good approximation ability, with convex, monotonic and linear regenerative properties, and can fit a smooth curve that matches the characteristics of the data; (3) For most of the actual signals, the triple spline interpolation has serious overshoot phenomenon, which is difficult to reflect the local characteristics of the signal and does not have the characteristics of conformality. In this case, MQ interpolation is used to estimate the envelope to avoid the above problems. 4. Experimental Results and Performance Analysis In this section, six different types of signals are designed, and the performance of Cubic Spline Interpolation (CSI), MQ proposed interpolation and the proposed method (Improved MQ, IMQ) are compared and analyzed using the orthogonal exponent IO, mean square error MSE and energy error. The parameter value of MQ proposed interpolation is chosen as the best primary parameter value for EMD decomposition results. The frequencies of the six signals are 200 Hz, 200 Hz, 100 Hz, 20 Hz, 400 Hz and 100 Hz, and the expressions are defined as follows. 0.1 1( ) sin(2 ) , [0,30]tx t t e t  1 2 3 2 2 2 2( ) ( ) ( ) ( ), [0,6]x t x t x t x t t    , 1 2 ( )=1.2sin(10 )x t t , 2 2 ( )=2cos(4 )x t t ; 3 2 ( )= cos(2 )x t t ; 1 2 3 3 3( ) ( ) ( ), [1,5]x t x t x t t   , 1 3( )=sin(200 )x t t , 2 3 ( )=sin(20 )x t t ; 1 2 4 4 4( ) ( ) ( ), [1,20]x t x t x t t   , 1 4 ( )=5sin( )x t t ; 5 ( )=(1+0.2sin(15 ))(sin(100 ) 1.5cos(60 ) sin(30 )), [0,1.5]x t t t t t t      ; 1 2 3 4 6 6 6 6 6( ) ( ) ( )+ ( )+ ( ), [1,20]x t x t x t x t x t t   , 1 6 ( )=2sin(250 t)x t  ; 2 2 6 ( )=3cos( (1.7 7.3) )x t t  , 3 0.23(1 ) 2 6 ( )=e cos( (2.58 21.95) )tx t t  , 4 6 ( )=3tx t . For the first signal, according to Figure 1(a), CSI performs the best, followed by IMQ, and MQ performs the worst. Due to the fixed parameters of MQ, the extracted IMF2, IMF3 and IMF4 components of MQ exhibit no smoothness. On the contrary, the improved IMQ envelope fitting algorithm improves most of the IMF components, and the EMD decomposition results are more satisfactory compared to MQ. From the qualitative point of view, we can also know that IMQ has better quadrature index, mean square error and energy error than MQ. However, it is slightly weaker than CSI because the signal is a damped sinusoidal signal with no local violent oscillations and a clear pattern of variation of local extrema, so the ideal envelope can be obtained directly by the third spline interpolation. In contrast, IMQ is not an interpolation method in essence, but only an approximate envelope can be obtained. The second signal is composed of high, medium and low frequency components, which also has a clear pattern, and the change of local extreme points is relatively flat. Based on the EMD decomposition results and the three metrics, we know that CSI performs best, IMQ outperforms MQ, and MQ performs the worst. For the third signal, we consider a signal consisting of high and low frequencies, which is dominated by the low frequency component. As shown in Figure 1(c), the overall performance of the EMD decomposition of the three algorithms is satisfactory, and the low-frequency components can be extracted effectively, and the IMF extracted by IMQ can be smoothed faster than CSI and MQ. The three metrics show that the performance of IMQ is better than the other two methods. It is worth mentioning that the signal appears to be a smooth and regular time series from the subjective view, but in fact, the signal amplitude fluctuates slightly due to the influence of high frequency components, which leads to severe overshoot in the third spline interpolation and increases the envelope error. In contrast, IMQ is a method of envelope estimation, which can mitigate the effects of small changes in amplitude and is strictly conformal, and the obtained envelope is more advantageous. Therefore, IMQ performs better than CSI in the signal. 38 (a) (b) (c) (d) (e) (f) Figure 1. EMD decomposition results of three envelope fitting algorithms Table 1. Comparison of IO results of three different extension methods IO CSI MQ IMQ )(1 tx )(2 tx )(3 tx )(4 tx )(5 tx )(6 tx 0.0011 0.1375 2.5964e-15 0.6055 0.0693 0.0201 0.0017 0.2023 1.7829e-15 0.2573 0.1075 0.0205 0.0015 0.1945 1.2746e-15 0.2070 0.0332 0.0091 Table 2. Comparison of MSE results of three different extension methods MSE CSI MQ IMQ ))(,( 11 txIMFMSE ))(,( 1 21 txIMFMSE ))(,( 2 31 txIMFMSE ))(,( 1 42 txIMFMSE ))(,( 51 txIMFMSE ))(,( 3 61 txIMFMSE ))(,( 2 62 txIMFMSE 0.0016 0.0790 2.5563e-26 17.1415 0.2393 12.6679 10.1776 0.0024 0.4602 2.3401e-26 8.2238 0.1583 10.4190 7.9031 0.0020 0.8458 2.2477e-26 9.4408 0.0575 7.6611 5.7826 Table 3. Comparison of  results of three different extension methods  CSI MQ IMQ )(1 tx )(2 tx )(3 tx )(4 tx )(5 tx )(6 tx 0.0011 0.1388 1.5721e-16 0.4176 0.0211 0.0139 0.0017 0.2109 9.4323e-16 0.0481 0.0748 2.1581e-04 0.0014 0.2047 9.4323e-16 0.0138 0.0109 0.0073 The signals, and are non-stationary and have no obvious variation pattern. Qualitative analysis: In Fig. 1 (d), IMF2 and IMF3 extracted by CSI-based algorithm have obvious modal blending, and IMF extracted by MQ-based algorithm has poor 39 smoothing. From Fig. 1(e), it can be found that all three methods cannot decompose the original signal into IMF components of different scales, and the extracted IMF1 all basically overlap with the original signal. For the signal, the IMF3 based on CSI and MQ screening has no clear physical meaning, and the number of IMQ screening IMF iterations is less compared to the remaining two methods. From the quantitative point of view, the overall performance of IMQ is better than CSI and MQ, and only slightly higher in energy than MQ in Table 1. Based on the analysis of the above six signals, it can be concluded that IMQ outperforms CSI and MQ for non- stationary and irregular signals; IMQ still outperforms MQ and is slightly weaker than CSI for signals with significant variation patterns. Most of the actual signals contain noise, so IMQ is more advantageous than the other two methods. 5. Conclusion Regarding the envelope fitting problem, an improved MQ fitted interpolation approximation envelope algorithm is proposed. This method defines a new way of parameter setting, and IMQ has a stronger adaptive capability than the constant-value parameter proposed by J.H. Zhang et al. (2012). The experimental results show that the overall performance of the proposed method is better than the other two methods for non-smooth and irregular real signals. References [1] Rilling G, Flandrin P, Goncalves P, On empirical mode decomposition and its algorithms. IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing, 2003, NSIP-03, 1- 5. [2] Datig M, Schlurmann T. Performance and limitations of the Hilbert-Huang transformation (HHT) with an application to irregular water waves. Ocean Engineering, 2004, 31: 1783– 1834. [3] Wu F J, Qu L S. An improved method for restraining the end effect in empirical mode decomposition and its applications to the fault diagnosis of large rotating machinery. Journal of Sound and Vibration, 2008, 314: 586-602. [4] Xu W X, Chen S H, Wang M Z, et al. Eliminating the end effect of empirical mode decomposition using a cubic spline based method. Digital Signal Processing, 2021, 110: 102936. [5] Mandic D, Rehman N, Wu Z, et al. Empirical mode decomposition-based time-frequency analysis of multivariate signals: the power of adaptive data analysis. IEEE Signal Processing Magazine, 2013, 30(6): 74-86. [6] Li H Y, Qin X Y, Zhao D, et al. An improved empirical mode decomposition method based on the cubic trigonometric B- spline interpolation algorithm. Applied Mathematics and Computation, 2018, 332: 406-419. [7] Huang N E, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non- stationary time series analysis. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1998, 454: 903-995. [8] Zhang J H. Empirical modal decomposition and some applications of radial basis functions. Dalian University of Technology, 2012.(In Chinese)