Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
  • Open Access

Nonlinear mode decomposition: A noise-robust, adaptive decomposition method

Dmytro Iatsenko, Peter V. E. McClintock, and Aneta Stefanovska
Phys. Rev. E 92, 032916 – Published 29 September 2015
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • BACKGROUND AND TERMINOLOGY
  • NONLINEAR MODE DECOMPOSITION: BASIC IDEA…
  • ADDITIONAL IMPROVEMENTS
  • METHOD PARAMETERS AND THEIR CHOICE
  • SIMULATED EXAMPLES
  • REAL LIFE EXAMPLES
  • OTHER APPLICATIONS
  • CONCLUSIONS
  • ACKNOWLEDGMENTS
  • APPENDICES
    • References

    Abstract

    The signals emanating from complex systems are usually composed of a mixture of different oscillations which, for a reliable analysis, should be separated from each other and from the inevitable background of noise. Here we introduce an adaptive decomposition tool—nonlinear mode decomposition (NMD)—which decomposes a given signal into a set of physically meaningful oscillations for any wave form, simultaneously removing the noise. NMD is based on the powerful combination of time-frequency analysis techniques—which, together with the adaptive choice of their parameters, make it extremely noise robust—and surrogate data tests used to identify interdependent oscillations and to distinguish deterministic from random activity. We illustrate the application of NMD to both simulated and real signals and demonstrate its qualitative and quantitative superiority over other approaches, such as (ensemble) empirical mode decomposition, Karhunen-Loève expansion, and independent component analysis. We point out that NMD is likely to be applicable and useful in many different areas of research, such as geophysics, finance, and the life sciences. The necessary matlab codes for running NMD are freely available for download.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    6 More
    • Received 1 April 2015

    DOI:https://doi.org/10.1103/PhysRevE.92.032916

    This article is available under the terms of the Creative Commons Attribution 3.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI.

    Published by the American Physical Society

    Authors & Affiliations

    Dmytro Iatsenko, Peter V. E. McClintock, and Aneta Stefanovska*

    • Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom

    • *aneta@lancaster.ac.uk

    Article Text

    Click to Expand

    References

    Click to Expand
    Issue

    Vol. 92, Iss. 3 — September 2015

    Reuse & Permissions
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      (a) The central 50-s section of the signal s(t) specified at the top, which is composed of three oscillations corrupted with white Gaussian noise of standard deviation equal to 0.5. The oscillations are (1) near 10 Hz, with frequency modulation (chirp); (2) near 6 Hz, with amplitude modulation; (3) near 1 Hz, with both modulations and a nonsinusoidal wave form. (b),(c) Respectively, the central parts of the WFT and WT of the signal s(t) shown in (a), where the ridge curves related to each oscillation are shown in each case by solid lines of the corresponding colors: 1, magenta (top line); 2, - green (second from the top); 3, brown (two bottom lines); the third oscillation is represented by two curves because of its complex wave form, while the apparent distortions in the ridge frequency profiles are due to noise; note also the logarithmic frequency scale for the WT, which is natural for it. The signal was sampled at 100 Hz for 100 s.

      Reuse & Permissions
    • Figure 2
      Figure 2

      Illustration of the extraction of the hth harmonic ridge curve ωp(h)(t) based on the fundamental ridge frequency ωp(1)(t). At each time, starting from the expected ridge frequency hωp(1)(t) of the harmonic (blue diamonds), one climbs (i.e., follows in the direction of TFR amplitude increase) to the nearest peak, which is then assigned to ωp(h)(t) (red circles).

      Reuse & Permissions
    • Figure 3
      Figure 3

      (a) Nonlinear mode with amplitude modulation, as specified by the equation for s(t) above the box. (b),(c) The corresponding WFT and WT amplitudes, respectively. (d),(e) The signal's FT, integrated over one frequency bin (vertical lines), with the parts responsible for different harmonics shown using different colors and markers; the shaded areas show the absolute values of the window functions ĝ(ω−ξ) or wavelet functions ψ̂*(ωψξ/ω) centered at the mean frequencies of the harmonics ω=2πh and rescaled to half of the harmonics' mean amplitudes. The signal was sampled at 50 Hz for 50 s.

      Reuse & Permissions
    • Figure 4
      Figure 4

      (a) Nonlinear mode with frequency modulation, as specified by the equation for s(t) above the box. (b),(c) The corresponding WFT and WT amplitudes, respectively. (d),(e) The signal's FT, integrated over one frequency bin (vertical lines), with the parts responsible for different harmonics shown using different colors and markers; the shaded areas show the absolute values of the window functions ĝ(ω−ξ) or wavelet functions ψ̂*(ωψξ/ω) centered at the mean frequencies of the harmonics ω=2πh and rescaled to half of the harmonics' mean amplitudes. The signal was sampled at 50 Hz for 55 s.

      Reuse & Permissions
    • Figure 5
      Figure 5

      (a) The central 50-s section of the signal s(t) specified by the equation at the top, which represents a single NM corrupted by Brownian noise of standard deviation equal to 4 [Brownian noise of unit deviation ηB(tn)∼∑m=1nηW(tm−1)tm−tm−1 is obtained as the normalized cumulative sum of the Gaussian white noise signal ηW(t)]. (b),(c) Central parts of the WFT and WT, respectively, of the signal s(t) shown in (a). The signal was sampled at 100 Hz for 100 s.

      Reuse & Permissions
    • Figure 6
      Figure 6

      The WFTs from which each true harmonic is reconstructed, with the corresponding extracted ridge curves being shown by solid red lines. The window resolution parameter f0 is adjusted individually for each harmonic. After the harmonic is reconstructed, it is subtracted from the signal, so that, e.g., the first harmonic no longer appears in the WFTs of (b)–(d). Note that the color scaling in (b)–(d) covers half the amplitude range of that in (a).

      Reuse & Permissions
    • Figure 7
      Figure 7

      The result of applying NMD to the signal shown in Fig. 5. (a) The reconstructed NM (black line) is compared with the true NM (gray background line); (b) the residual provided by NMD (black line) is compared with the actual background noise (gray background line).

      Reuse & Permissions
    • Figure 8
      Figure 8

      The result of applying EMD (left panels) and EEMD (right panels) to the signal shown in Fig. 5. Thin red lines, where present, show the real first harmonic (in C4 for EMD and C5 for EEMD), sum of the third and fourth harmonics (in C3 for EMD and C4 for EEMD), and the seventh harmonic (in C2 for EMD and C3 for EEMD). The bottom panels show the sum of components 7 to 13. For EEMD, we have used 1000 noise realizations with standard deviations ten times smaller than that of the signal.

      Reuse & Permissions
    • Figure 9
      Figure 9

      (a) The central 20-s section of the signal s(t) specified by the equation at the top, representing the sum of two NMs corrupted by Gaussian white noise of standard deviation equal to 1.725; the phases of the NMs ϕ1,2(t) were obtained as ϕ1,2(t)=∫0tν1,2(τ)dτ, with ν1,2(t)/2π=(1,2)+0.01η̃B;1,2(t), where η̃B;1,2(t) are two independent realizations of the unit deviation Brownian noise filtered in the range [0.01,0.2] Hz. (b),(c) Central parts of the WFT and WT of the signal s(t) shown in (a). The signal was sampled at 100 Hz for 100 s.

      Reuse & Permissions
    • Figure 10
      Figure 10

      The result of applying NMD to the signal shown in Fig. 9. In (a) and (b) black lines indicate the two reconstructed NMs and the gray background lines show the true NMs for comparison. (c) Similarly, the black and gray lines show the residual returned by NMD and true background noise, respectively.

      Reuse & Permissions
    • Figure 11
      Figure 11

      (a) An example of the blood flow signal, measured by LDF with the probe over the right wrist caput ulna (for more details, see [3, 53]); the signal is sampled at 40 Hz for 1800 s, and the panel shows its central 600 s part. (b) The WT of the signal shown in (a). Gray dotted lines partition the frequency axis into the regions within which the physiologically meaningful oscillations are located (according to [3, 52, 53, 54]; see the text). Bold colored lines show those extracted components which pass the surrogate test against noise, with thinner lines of the same color showing their higher harmonics. (c),(d) The reconstructed NMs corresponding to four upper magenta lines and two bottom brown lines on (b), respectively; the main graphs show the modes during 50 s and small insets during 600 s [as in (a)]. (e),(f) The wave forms of the oscillations shown in (c) and (d), respectively; the gray dashed line shows a pure sinusoidal wave form of amplitude equal to that of the first harmonic and is provided as a guide for the eye. The cardiac wave form (e) has four harmonics h=1,2,3,4, with ah=[1,0.52,0.37,0.16] and φh/π=[0,0.4,0.66,−0.97] [in the notation of (4)]; the wave form in (f) is characterized by h=1,2, ah=[1,0.08], and φh/π=[0,0.07].

      Reuse & Permissions
    • Figure 12
      Figure 12

      Same as in Fig. 11, but for the blood flow measured from a different subject. The nonlinear modes shown in (c) and (d) correspond, respectively, to the four upper lines (magenta) and the bottom (green) line in (b). The cardiac wave form (e) has four harmonics h=1,2,3,4, with ah=[1,0.32,0.15,0.04] and φh/π=[0,0.41,0.95,−0.44] [in the notation of (4)]; the respiratory oscillations have only one harmonic, so that the wave form in (f) is a pure sinusoid.

      Reuse & Permissions
    • Figure 13
      Figure 13

      (a) The EEG signal [measured using a BISTM 4-Electrode Sensor (Covidien-Medtronic, Dublin, Ireland and Fridley, Minnesota) placed on the forehead and a purpose-built signal conditioning unit (Jožef Stefan Institute, Slovenia), as described in 59]. (b) The ECG signal (three-lead, with electrodes on shoulders and the lowest left rib, see, e.g., [60]). (c),(d) The WFTs of the EEG and ECG signals shown in (a) and (b), respectively, calculated using the default resolution parameter f0=1. (e) The WFT of the EEG signal calculated using the adaptively adjusted resolution parameter f0=3.76, due to which the cardiac component becomes much more visible than in (c). (f) The cardiac artifacts extracted from the EEG, with an inset showing them over the same 500 s as presented in (a),(b); gray background lines show the ECG scaled to the dimensions of the plot. (g) The wave form of the cardiac artifacts, which has four harmonics with ah=[1,0.33,0.24,0.11] and φh/π=[0,0.49,1,−0.59] [in the notation of (4)]; the gray dashed line represents a pure sinusoid with the amplitude of the first harmonic, for comparison. The full signals are of duration 20 min and are sampled at 80 Hz; the EEG was actually measured under anaesthesia, but the same artifacts arise under all conditions.

      Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review E

    Reuse & Permissions

    It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 3.0 License. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

    ×

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×