A Caltech Library Service

Data-driven atomic decomposition via frequency extraction of intrinsic mode functions

Chui, Charles K. and Mhaskar, H. N. and van der Walt, Maria D. (2016) Data-driven atomic decomposition via frequency extraction of intrinsic mode functions. GEM - International Journal on Geomathematics, 7 (1). pp. 117-146. ISSN 1869-2672. doi:10.1007/s13137-015-0079-3.

Full text is not posted in this repository. Consult Related URLs below.

Use this Persistent URL to link to this item:


Decomposition of functions in terms of their primary building blocks is one of the most fundamental problems in mathematical analysis and its applications. Indeed, atomic decomposition of functions in the Hardy space H^p for 0<p≤1 as infinite series of “atoms” that have the property of vanishing moments with order at least up to 1 / p has significant impacts, not only to the advances of harmonic and functional analyses, but also to the birth of wavelet analysis, which in turn allows the construction of sufficiently large dictionaries of wavelet-like basis functions for the success of atomic decomposition of more general functions or signals, by such mathematical tools as “basis pursuit” and “nonlinear basis pursuit”. However, such dictionaries are necessarily huge for atomic decomposition of real-world signals. The spirit of the present paper is to construct the atoms directly from the data, without relying on a large dictionary. Following Gabor, the starting point of this line of thought is to observe that “any” signal a can be written as a(t)=A(t)cosϕ(t) via complex extension using the Hilbert transform. Hence, if a given signal f has been decomposed by whatever available methods or schemes, as the sum of sub-signals f_k, then each sub-signal can be written as f_k(t)=Ak(t)cosϕ_k(t). Whether or not f_k is an atom of the given signal f depends on whether any of the sub-signals f_k can be further decomposed in a meaningful way. In this regard, the most popular decomposition scheme in the current literature is the sifting process of the empirical mode decomposition (EMD), where the sub-signals f_k are called intrinsic mode functions (IMF’s). The main contribution of our present paper is firstly to demonstrate that IMF’s may not be atoms, and secondly to give a computational scheme for decomposing such IMF’s into finer and meaningful signal building blocks. Our innovation is to apply the signal separation operator (SSO), introduced by the first two authors, with a clever choice of parameters, first to extract the instantaneous frequencies (IF’s) of each IMF obtained from the sifting process, and then (by using the same parameters for the SSO, with the IF’s as input) to construct finer signal building blocks of the IMF. In other words, we replace the Hilbert transform of the EMD scheme by the SSO in this present paper, first for frequency extraction, and then for constructing finer signal building blocks. As an example, we consider the problem in super-resolution of separating two Dirac delta functions that are arbitrarily close to each other. This problem is equivalent to finding the two cosine building blocks of a two-tone signal with frequencies that are arbitrarily close. While the sifting process can only yield one IMF when the frequencies are too close together, the SSO applied to this IMF extracts the two frequencies and recover the two cosine building blocks (or atoms). For this reason, we coin our scheme of sifting + SSO as “superEMD”, where “super” is used as an abbreviation of super-resolution.

Item Type:Article
Related URLs:
URLURL TypeDescription ReadCube access
Additional Information:© 2016 Springer-Verlag Berlin Heidelberg. Received: 1 December 2015; Accepted: 22 December 2015; Published online: 9 January 2016. The research of the first two authors is supported in part by Grant W911NF-15-1-0385 from the U. S. Army Research Office.
Funding AgencyGrant Number
Army Research Office (ARO)W911NF-15-1-0385
Subject Keywords:Atomic decomposition; Dictionaries for basis pursuit; Data-driven atomic decomposition; Local sifting process; Instantaneous frequencies; Empirical mode decomposition (EMD); Synchrosqueezing transform (SST); Hilbert transform; Intrinsic mode functions; Signal separation operator (SSO); Super-resolution; SuperEMD
Issue or Number:1
Classification Code:Mathematics Subject Classification: 94A12
Record Number:CaltechAUTHORS:20160707-100242858
Persistent URL:
Official Citation:Chui, C et al. Int J Geomath (2016) 7: 117. doi:10.1007/s13137-015-0079-3
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:68885
Deposited By: Tony Diaz
Deposited On:07 Jul 2016 19:31
Last Modified:11 Nov 2021 04:06

Repository Staff Only: item control page