A Modified Complex Variational Mode Decomposition Method for Analyzing Nonstationary Signals with the Low-Frequency Trend

Complex variational mode decomposition (CVMD) has been proposed to extend the original variational mode decomposition (VMD) algorithm to analyze complex-valued data. Conventionally, CVMD divides complex-valued data into positive and negative frequency components using bandpass filters, which leads to difficulties in decomposing signals with the low-frequency trend. Moreover, both decomposition number parameters of positive and negative frequency components are required as prior knowledge in CVMD, which is difficult to satisfy in practice. This paper proposes a modified complex variational mode decomposition (MCVMD) method. First, the complex-valued data are upsampled through zero padding in the frequency domain. Second, the negative frequency component of upsampled data are shifted to be positive. Properties of analytical signals are used to get the real-valued data for standard variational mode decomposition and the complex-valued decomposition results after frequency shifting back. Compared with the conventional method, the MCVMD method gives a better decomposition of the low-frequency signal and requires less prior knowledge about the decomposition number. The equivalent filter bank structure is illustrated to analyze the behavior of MCVMD, and the MCVMD bi-directional Hilbert spectrum is provided to give the time–frequency representation. The effectiveness of the proposed algorithm is verified by both synthetic and real-world complex-valued signals.


Introduction
Nonstationary signal processing methods have attracted considerable interest over the last few decades. Typical methods such as short-time Fourier transform (STFT) [1], the continuous wavelet transform (CWT) [2], S transform (ST) [3], and Wigner-Ville distribution (WVD) [4] can provide a time-frequency or time-scale representation for analyzing nonstationary signals. However, these algorithms are based on specific bases. Their performance depends on the selection of basis and is restricted by the uncertainty principle [5]. In order to break the limitation of the uncertainty principle, the empirical mode decomposition (EMD) [6] was proposed by Huang et al. as a data-driven signal decomposition method. EMD can recursively decompose a nonstationary signal into data-dependent basis functions named the intrinsic mode function (IMF) [7]. Consequently, the instantaneous frequency characteristics of the nonstationary data can be directly obtained from the IMFs [7]. However, EMD lacks mathematical supports and has limitations such as mode mixing. Therefore, variational mode decomposition (VMD) [8] was proposed in 2014 as a newly developed method alternative to EMD, fully established on a mathematical framework [9,10]. It can adaptively decompose a multi-component signal into several sub-signals and is not restricted by the uncertainty principle. VMD has been proven to outperform EMD and many other time-frequency analysis techniques in many applications, such as detection, and analysis to verify the performance of the proposed method. At last, the conclusion is done in Section 5.

Variational Mode Decomposition
VMD was first proposed by Dragomiretskiy and Zosso in 2014, a novel signal decomposition method for time-frequency analysis [8]. It can non-recursively decompose a real-valued multi-component signal f into K number of quasi-orthogonal band-limited sub-signals u k , also called band-limited intrinsic mode functions (BLIMFs). The VMD method is realized through the constrained variational problem: where {u k (t)} = {u 1 , . . . , u K } are the BLIMFs to be solved, {ω k (t)} = {ω 1 , . . . , ω K } are the center frequencies of the BLIMFs to be solved, and K is the decomposition number that is known as an a priori parameter. In Equation (1), the δ(t) + j πt * u k (t) part represents the transformation of the real-valued BLIMF u k (t) into its analytical signal. Then the analytical signals are shifted to the zero-frequency baseband, and finally, the effective bandwidth is obtained using the L2 norm of the time derivative. The quadratic penalty and Lagrangian multiplier are introduced to address Equation (1) by the augmented Lagrange L as follows: (2) where α represents the balancing parameter of the data-fidelity constraint and λ(t) denotes the Lagrangian multipliers. Equation (2) is then solved through the alternate direction method of multipliers (ADMM) approach. Using the condition that the signal is real-valued, the BLIMFs in spectral domain can be solved bŷ where the ω k is calculated at the center of gravity of the corresponding mode's power spectrum as The BLIMFs in time domain {u k (t)} are obtained from the real part of the inverse Fourier transform of the frequency spectra.
As the standard VMD has been a mature algorithm built as a function in MATLAB since 2020, the subsequent VMD procedure is based on the MATLAB function. Without a special note, all of the parameters, such as the way of initializing the center frequencies ω k and the parameter α, are set as default.

Complex Variational Mode Decomposition
Since the VMD is designed for real-valued data, the extension method of VMD for handling complex-valued data called CVMD was proposed. CVMD divides the complexvalued signal into positive and negative frequency components and uses the relationship between the Fourier spectrum and the analytical signal to get real signals for the standard VMD algorithm. Let x(t) be the complex-valued data and X e jω , where −π ≤ ω < π is the Fourier transform of x(t). The complex-valued data are separated into positive and negative components by using a band-pass filter as The following two real signals can be generated by where * describes the complex conjugate operation, symbol R{·} defines the real part of the complex-valued data, and F −1 [·] denotes the inverse Fourier transform. According to the properties of the analytic signals, F −1 H e jω X e jω and F −1 H e jω X * e −jω are analytic signals, and their real parts contain all of the information. The reconstruction of the original complex-valued signal is determined by x + (t) and x − (t) as well as their Hilbert transform pairs, expressed as where [·] denotes the Hilbert transform operator. Since x + (n) and x − (n) are real-valued signals, the corresponding BLIMFs can be obtained using standard VMD. The VMD is separately embraced to the two real-valued signals. Suppose the decomposition number of the positive and negative frequency planes is K + and K − , respectively. The existing methods consider K + = K − = K, so the decomposition is described as where {x i (t)} K i=1 and {x i (t)} −1 i=−K denote sets of BLIMFs corresponding to x + (t) and x − (t), respectively.
Based on the Equations (8)-(10), the CVMD for a complex-valued data x(t) is eventually expressed as where z i (t) represents the i-th complex BLIMF, which can be written as

Problem Statement
The VMD mathematical theory demonstrates that the algorithm is unsuitable for complex-valued signal analysis. CVMD extent VMD to complex-valued data, but here we will illustrate its limitations in decomposing the signal with a low-frequency trend. Take the complex-valued signal s 1 (t) = 0.5e j2π(−8t) + e j2π(−t+0.6t 2 ) + 0.8e j2π(16t+0.2t 2 ) for example. The blue solid line in Figure 1 shows the frequency spectrum of s 1 (t). Three band-limited components can be seen distributed in different parts of the spectrum: the first component, x 1 = 0.5e j2π(−8t) , is mainly at negative frequency; the second component, named x 2 = e j2π(−t+0.6t 2 ) , is around zero frequency; and the third one, x 3 = 0.8e j2π(16t+0.2t 2 ) , is mainly at positive frequency. Note that the center frequency of the low-frequency signal x 2 is actually at the positive frequency, but there is energy spreading over negative frequencies. In CVMD, the positive and negative frequency components of the signal are displayed on the positive and negative frequency planes, respectively. Both of the decomposition numbers of two planes are required. The decomposition number of the negative frequency plane is considered the same as that of the positive frequency plane. However, only the whole data's decomposition number can be easily known in this example. It is unclear how many components are at the positive frequency plane or negative frequency plane without prior knowledge. Even if the correct decomposition number is known and the signals are rearranged well, the results by CVMD have energy loss in the components around zero frequency. The red dotted line in Figure 1 shows the reconstruction of the CVMD results. It can be seen clearly that signal x 2 is not recovered completely.
To further illustrate this problem, the behavior of CVMD as a filter bank is illustrated by applying the algorithm to the white Gaussian noise. Both the positive frequency and negative frequency components' decomposition number is set as 6, for example. The equivalent filter bank of CVMD is shown in Figure 2, in which the parts near zero frequency, x 1+ and x 1− , are separated artificially by the filter. This separation goes against the data-driven concept of VMD. Although the low-frequency trend can be reconstructed by combining x 1+ and x 1− components, it will not be recovered ideally.

Modified Complex Variational Mode Decomposition
CVMD cannot reconstruct the low-frequency trend completely because it separates the data into positive and negative frequency components artificially and hurts the parts around zero frequency. A modified VMD method is proposed to overcome the defect based on upsampling and frequency shifting.

MCVMD Algorithm
For real-valued data, the positive-and negative frequency components are conjugate symmetric. For complex-valued data, its positive-and negative frequency components contain different information, which is not symmetric. The idea of the proposed MCVMD method is to convert complex-valued data into real-valued data without losing information so that standard VMD can be used. This idea is similar to CVMD and CEMD, but the operation is different. The well-known property of the analytic signal utilized in the proposed method is that the analytic signal of a real-valued signal only has positive frequency spectra but contains all of the information of the origin real-valued data. Alternatively, the real part of the analytic signal has symmetrical positive and negative frequencies and contains all of the origin complex-valued analytic signal information. According to the property, if the whole spectrum of the complex-valued data are shifted to the positive frequency, then the complex-valued data become an analytic signal. The real-valued data can be obtained by taking the real part of the analytic signal and contain all of the information of the original complex-valued data. Since the signal in practical applications is almost discrete data, the proposed method is designed for discrete signals. The steps of the algorithm are listed as follows: Consider x(n), where 0 ≤ n ≤ N − 1 and n is an integer, as the discrete complexvalued data, and the sampling frequency rate is normalized as 1. Suppose the sampling rate satisfies the Nyquist sampling frequency. That is to say, the max signal frequency of the data do not exceed 1/2, and the minimum signal frequency is not less than −1/2.
First, the input data should be upsampled by the interpolation algorithm, among which the discrete sinc interpolation algorithm does not distort the signal as defined and is entirely reversible. However, it is often not used directly due to the boundary effects of sinc interpolation [31]. One of the simplest and most efficient methods for minimizing boundary effects is the mirror extension of the signal by half its length on each side. Accordingly, the mirror extension of the original input data are expressed as where x e1 denotes the refection sequence of the first half of the original data as where x e2 denotes the original sequence as where x e3 denotes the refection sequence of the last half of the original data as The extended signal can be rewritten as x e (n), where 0 ≤ n ≤ 2N − 1.
The discrete Fourier transform (DFT) is applied to the extended signal as This represents the spectral samples. Then, 2N samples valued at zero are inserted after the first N spectral samples. A one-dimensional vector with length 4N can be generated as where X e1 denotes the sequence X e (k) (0 ≤ k ≤ N − 1), Z denotes 2N samples valued at zero, and X e2 denotes the sequence X e (k) (N ≤ k ≤ 2N − 1), which can be rewritten as The corresponding time-domain data can be obtained by where F −1 [·] denotes the inverse discrete Fourier transform. The actual upsampled result is given by removing the mirror part as It is known that zero padding in the frequency domain leads to an increased sampling rate in the time domain [32]. Since the data length is doubled by zero padding, the sampling rate of x z (n) now becomes 2. In addition, the Fourier spectrum of x z (n) now has definitions from the frequency of -1 to 1 and has values from the frequency of -1/2 to 1/2.
x z (n) is then shifted by multiplying exponential factors e j2π· n 4 as follows Since the frequency of x z (n) is within −1/2 to 1/2, all of the components of x zs (n) are within 0 to 1 after frequency shifting, which means it has only positive frequency components. Thus, x zs (n) becomes an analytic signal. According to the property of the analytic signals, real parts of the analytic signal contain all of the information. A real-valued signal can be deduced as x zsr (n) = R{x zs (n)}.
Since x zsr (n) is real valued, the BLIMFs can be obtained using standard VMD. Consider that the decomposition number is K. The mode decomposition procedure can be expressed as where u i (n) denotes the i-th BLIMF corresponding to x zsr (n), and r(n) describes the residual. After VMD, the Hilbert transform is applied to generate the corresponding complex and analytic i-th BLIMF as Note the data are upsampled and modulated, so they are not the final result of MCVMD. The final complex BLIMFs can be obtained by demodulation and downsampling. The i-th upsampled and modulated complex BLIMFs y i (n) are shifted back by multiplying exponential factors e −j2π· n 4 as where z i (n) denotes the i-th upsampled complex BLIMFs, which has the double sample rate of the original data. The final i-th complex BLIMF q i (n) corresponding to the origin complex data x(n) can be obtained simply by downsampling.

MCVMD Hilbert Spectrum
The mode decomposition method is often used to get the Hilbert-Huang spectrum, which gives the time-frequency representation [33]. Specifically, mode decomposition is the first step in conventional Hilbert-Huang spectrum analysis. The Hilbert transform is then applied to the real-valued decomposed modes to get the analytic signal. Analytic signals are required to analyze each mode's instantaneous phase and amplitude [34]. Here, the complex BLIMF from MCVMD is already analytical. Thus, the instantaneous phase and instantaneous amplitude can respectively be obtained directly by the following equations where symbol I{·} defines the imaginary part of a complex function, ∅ i (n) represents the instantaneous phase, and A i (n) corresponds to the instantaneous amplitude of the i-th complex BLIMF q i (n). The instantaneous frequency can be derived as Accordingly, the Hilbert spectrum of each complex BLIMF can be obtained as The MCVMD Hilbert spectrum intuitively shows the instantaneous frequency and amplitude of each mode in the signal, which is not constrained by the uncertainty limitations. Unlike the VMD Hilbert spectrum that only contains positive frequency, the MCVMD Hilbert spectrum can represent the positive and negative frequencies. Thus, the time-varying characteristics of low-frequency trends varying from positive frequency to negative frequency or from negative to positive can also be seen.

MCVMD Equivalent Filter Bank
CVMD has been applied to white Gaussian noise to analyze its equivalent filter bank. In this subsection, the equivalent filter bank of MCVMD is also investigated. MCVMD is applied to Gaussian complex-valued data of 1024 randomly generated samples, and the decomposition number is set as 12. Thus, 12 complex BLIMFs and their spectra can be obtained. This process is repeated 2000 times independently, and then the resulting spectra are averaged. Figure 3 shows these averaged spectra for the 12 complex BLIMFs. It can be seen CVMD equivalent filter bank structures behave almost constant bandwidths. Compared to the equivalent filter bank of CVMD shown in Figure 2, the low-frequency components are no longer separated. Each resulting structure by MCVMD is natural and complete over the whole frequency range without artificial separation. The result exhibits the desired behavior of a filter bank, which confirms the proposed MCVMD being a natural extension of the standard VMD and outperforming the CVMD.

Example 1
In this example, the performance of MCVMD in decomposing the signal with a lowfrequency component is verified. Consider the signal s 1 (t) = 0.5e j2π(−8t) + e j2π(−t+0.6t 2 ) + 0.8e j2π(16t+0.2t 2 ) mentioned in Section 3. Suppose the sampling rate is 50 Hz and the number of samples is 150. The discrete signal can be written as s 1 (n) = 0.5e j2π(−8·0.2n) + e j2π(−0.2n+0.6·(0.2n) 2 ) + 0.8e j2π(16·0.2n+0.2·(0.2n) 2 ) , where n = 1, 2, . . . , 150. The mode decomposition methods aim to decompose the signal into an ensemble of sub-signals to analyze the nonstationary data. Applying the MCVMD to s(n), three complex BLIMFs can be obtained. The solid blue line in Figure 4 shows the spectrum of the original signal, and the red dotted line represents the spectrum of the reconstruction signal obtained by summing up these three sub-signals. Compared with the corresponding CVMD result in Figure 1, MCVMD better reconstructs the original data, especially the second components around zero frequency. The time-domain waveforms of the three BLIMFs obtained by CVMD and MCVMD are compared in Figure 5. Let x 1 , x 2 , and x 3 describe the three constituent modes −0.5e j2π(−8t) , e j2π(−t+0.6t 2 ) , and 0.8e j2π(16t+0.2t 2 ) , respectively. Figure 5 shows the original modes (red solid line) and the decomposition results by CVMD (green dash-dot line) and MCVMD (blue dotted line). The real part and imaginary part of each component are displayed separately. It can be seen both CVMD and MCVMD give a proper decomposition of x 1 and x 3 . However, the performance of the two methods differs in reconstructing the second mode x 2 . MCVMD gives a better decomposition and reconstruction, especially in the imaginary part of x 2 . This phenomenon is consistent with the results of the spectrum display.
The root mean square error (RMSE) of each decomposed mode is shown in Table 1, in which the real and imaginary parts are considered separately. From Figure 5, both CVMD and MCVMD yield improper values at the beginning and end of the signal. Since this problem is because of the original VMD algorithm and is not the point here, the RMSE results in this paper are given by removing the first and last five samples. Compared with CVMD, the RMSE of the decomposed modes by MCVMD is generally smaller. As can be seen, the RMSE of the imaginary part of the low-frequency mode by MCVMD is almost one-tenth of that by CVMD. Table 1. The RMSE of the decomposed modes by CVMD and MCVMD.

Example 2
To further prove the proposed method's ability to analyze the signal's time-frequency characteristics, we give another more complex example of MCVMD developing the bidirectional Hilbert spectrum. The synthetic signal has two modes as follows s 2 (t) = cos(0.2t)·e j2π(−t+0.5t 2 ) + 0.2t·e j2π(28t−sin(πt)) , where both amplitudes and frequencies are time varying. The signal length is 5 s, and the sampling rate is 100 Hz. The amplitude of the first mode varies as the cosine function, and the frequency varies linearly between -1 Hz and 4 Hz. The amplitude of the second mode increases linearly with time, and the frequency varies as the sine function. Hilbert transform is used to get the two modes' instantaneous amplitudes and frequencies. Figure 6a shows the actual Hilbert spectrum of the two simulation sub-signals, which is used for comparison later. The proposed MCVMD method is applied to s 2 (t). After that, we use Equations (27)-(30) to get the corresponding Hilbert spectrum shown in Figure 6b. The MCVMD Hilbert spectrum presents the time-frequency characteristics and shows the two modes' instantaneous amplitude and frequency as the actual Hilbert spectrum. The corresponding results by CVMD are shown in Figure 6c. As can be seen, the instantaneous amplitude of the low-frequency mode by CVMD behaves differently from the actual mode.
The RMSE of each mode's instantaneous frequency and amplitude are presented in Table 2. It shows that both the instantaneous frequency and amplitude's RMSE by MCVMD are smaller than by CVMD, especially in the low-frequency mode. It can be concluded that MCVMD not only gives a better decomposition of the multi-component signal, but also gives a better estimation for instantaneous frequency and instantaneous amplitude with the corresponding Hilbert spectrum. Table 2. The RMSE of the instantaneous frequency and amplitude by CVMD and MCVMD.

Example 3
In this subsection, we discuss the prior requirements of the proposed method on the number of decompositions. The results obtained in Example 2 rely on the precise knowledge of the decomposition number K. In the MCVMD method, the correct decomposition number is K = 2 for the entire frequency range. In the CVMD method, the appropriate decomposition number is K + = 3 for the positive frequency range and K − = 1 for the negative frequency range. Numerical experiments are presented to find out the influence of incorrect decomposition numbers.
Considering the same simulation scenario as Example 2, we set the decomposition number K from 1 to 6 in MCVMD. Figure 7a shows the second mode is lost for the case K = 1, and when K > 1, all of the modes can be obtained as in Figure 7b,c. Although one component may be divided into more than one mode when K > 2 as shown in Figure 7c, the separated modes can be summed up to recover the component without loss. This separation is not like the separation by the filter in CVMD, which is natural and data driven. The RMSEs of each decomposed mode versus decomposition number are summarized in Figure 8a, in which the RMSEs of the real and imaginary parts of each mode are averaged. The performance remains superior as long as the decomposition number exceeds two. Thus, the requirement of the decomposition number for MCVMD, in this case, is no less than two. This prior requirement is not difficult to meet in practical applications.  For CVMD, there are many more corresponding decomposition schemes to be considered for the same situation, as shown in Table 3. That is because the CVMD needs to consider the positive and negative frequencies separately. In Reference [26], K − and K + are suggested to be the same number. Here, if we set K − and K + both equal to 2, the CVMD results are obtained in Figure 9a, which can be seen the second mode is lost. The two modes can be obtained completely if we set K − = 1 and K + = 3, as shown in Figure 9b. Numerical experiments of the 15 schemes show that modes can be obtained completely only when K − ≥ 1 and K + ≥ 3. Precisely, only 3 of the 15 schemes succeed to recover the two modes, with the RMSE less than 0.11, as shown in Figure 8b. Thus, the prior requirement of CVMD is much more challenging to meet. Table 3. Decomposition schemes for CVMD.

Analysis of Float Drift Data
We apply the proposed method to a real-world complex-valued data set to demonstrate its ability further. The data were taken during the Eastern Basin experiments, which can be downloaded from the World Ocean Circulation Experiment Subsurface Float Data Assembly Center (WFDAC) http://wfdac.whoi.edu, accessed on 23 January 2022 [24]. It contains the position record of a subsurface oceanographic float deployed in the North Atlantic Ocean to track the trajectory of salty water flowing from the Mediterranean Sea [35]. The looping trajectories of the chosen float are shown in Figure 10a (blue solid line). Here, the east and north displacements can be regarded as the real and imaginary part of the complex-valued data. The application of CVMD for analyzing the float drift data are illustrated in Reference [26]. Each BLIMF associated with the real and imaginary parts indicates a different physical path during vortex evolution. However, to characterize low-frequency trends, CVMD must combine the positive frequency BLIMFs and negative frequency BLIMFs. It can be seen from Reference [26], even if the signal is rearranged well, the low-frequency physical pathway given by CVMD is quite different from the actual situation, which fails to describe the potential physical characteristic. The complex-valued input data with 548 samples are shown in Figure 11 (top row), where the real part (blue line) and the imaginary part (red line) are displayed, respectively. The proposed MCVMD is applied to the data to analyze their principal modulated oscillations. The composition number is set as 4. Thus, four modes are shown as time plots ( Figure 11) and 2D graphical plots ( Figure 10). The 2D plot represents the drifting float's modulated oscillations or rotating modes. The displacement east is the real part of the resulting BLIMF, and the displacement north is the imaginary part of the resulting BLIMF. The red solid line in Figure 10a shows the low-frequency drift trajectory of the buoy, and the other three decomposed modes in Figure 10b-d show its specific rotation characteristics. We also apply CVMD to the same dataset. The yellow dotted line in Figure 10 shows the decomposed rotating modes by CVMD. It can be seen that the low-frequency drift trajectory of the buoy in Figure 10a,e by CVMD does not fully show the true characteristics. Again, that is because CVMD uses band filters and loses the low-frequency information. This example demonstrates the ability of MCVMD to analyze physically meaningful real-world data while also verifying the successful extension of VMD to complex-valued data.

Analysis of Wind Data
In this subsection, we apply the proposed method to the real-world wind data, whose measurements consist of wind direction ϕ(t) and wind speed v(t), which can be represented as a complex variable: x(t) = v(t)e j2πϕ(t)/360 [30]. If we apply VMD to speed and direction data, respectively, the signal components' positive and negative frequencies cannot be distinguished [25,26]. Thus, the complex extension of VMD is required to analyze the complex-valued form of wind data. We apply the CVMD and MCVMD methods to the wind data recorded by the Automated Weather Observing System (AWOS), which can be downloaded from http://mesonet.agron.iastate.edu/ASOS/, accessed on 23 January 2022. Figure 12 shows the time series of the wind speed and wind direction. The decomposition number K is set as 8 in MCVMD, and the positive decomposition number K + and negative decomposition number K − in CVMD are both set as 4. Figure 13 shows the complexvalued wind data decomposition results by MCVMD and CVMD, where the real part (blue line) and the imaginary part (red line) are displayed, respectively. From top to bottom corresponds to the mode with negative to positive frequency. In CVMD, the positive and negative frequency modes are obtained separately, in which the first four are the negative frequency, and the last four are positive. Summing up all of the decomposed modes, we can reconstruct the CVMD and MCVMD results. Figure 14 compares the spectrum of the original data, the reconstruction of CVMD results, and the reconstruction of MCVMD results. It can be seen that there is energy loss in the low-frequency component by CVMD, while MCVMD gives a better reconstruction. Similarly, that is because the CVMD uses band filters to separate the original data and go against the 'data-driven' idea of VMD. Thus, the low-frequency decomposed mode by CVMD cannot represent the wind's real low-frequency trend ideally. MCVMD gives better decomposition and reconstruction of the signals with low-frequency components.

Conclusions
This paper presents a modified complex variational mode decomposition (MCVMD) method for decomposing a complex-valued signal into an ensemble of complex bandlimited intrinsic mode functions (BLIMFs). The proposed approach extends standard variational mode decomposition (VMD) to complex-valued data. In contrast to the existing complex variational mode decomposition (CVMD) method, we refrain from separating the data into positive and negative frequency planes, applying VMD twice, and rearranging the positive and negative frequency parts. Instead, we upsample and modulate the input data to have only positive frequency components to apply the standard VMD. The properties of analytical signals are used to get the final complex-valued BLIMFs. Compared with CVMD, MCVMD has the following advantages. (1) The low-frequency components can be decomposed adaptively like other frequency parts, without artificial separation. (2) The MCVMD gives better decomposition and reconstruction of complex-valued signals with lower RMSE, especially of the low-frequency component. (3) VMD is only used once, resulting in a minor computation and less prior requirement for the decomposition number. (4) It is unnecessary to rearrange the decomposition results to reconstruct the low-frequency mode. The numerical experiment illustrates that the proposed MCVMD method behaves as an ideal filter bank structure over the whole positive and negative frequency plane. Applications in synthetic and real-world complex-valued signals demonstrate the effectiveness of the proposed MCVMD method.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.