Modulation Analysis with Higher-Order Spectra

Abstract

A need to identify modulation of spectral power arises frequently in the analysis of physiological signals. Estimation of power in the relevant bands through filtering and envelope extraction has several limitations: the choice of filter may bias any resulting estimate, while additive Gaussian noise becomes non-Gaussian due to the nonlinearity of envelope computation. The present work considers how spectral decompositions of higher-order cumulants (higher-order spectra, HOS) avoid these limitations, with an emphasis on the use of the trispectrum to identify modulated oscillations. Specifically, it is shown: 1) The trispectrum may be interpreted as a measure of linear dependencies of power across frequencies by viewing it as the cross spectrum of the Wigner-Ville distribution. 2) A particular two-dimensional subdomain of the trispectrum is useful for identifying modulated carriers, recovering essential spectral properties of both the modulating and carrier signals while avoiding the cubic complexity of full trispectrum estimation. A representation of this subdomain, the modulogram, is demonstrated as a tool for identifying and distinguishing different forms of modulation. 3) As a cumulant-derived measure, the modulogram is not biased by additive Gaussian noise. 4) Modulogram phase retains information by which temporal patterns of modulation may be identified. 5) A recently described additive decomposition of HOS (HOSD) further aids identification when applied to the trispectrum. These developments are illustrated with the blind detection of beta bursts in rodent and human local field potential recordings. Finally, the relationship between the present approach and prior techniques of blind identification (BI) through moment maximization, including blind deconvolution and independent component analysis, is considered.