Analyse des signaux AM-FM basée sur une version B-splines de l'EMD-ESA

In this paper a signal analysis framework for estimating time-varying amplitude and frequency functions of multicomponent amplitude and frequency modulated (AM–FM) signals is introduced. This framework is based on local and non-linear approaches, namely Energy Separation Algorithm (ESA) and Empirical Mode Decomposition (EMD). Conjunction of Discrete ESA (DESA) and EMD is called EMD–DESA. A new modified version of EMD where smoothing instead of an interpolation to construct the upper and lower envelopes of the signal is introduced. Since extracted IMFs are represented in terms of B-spline (BS) expansions, a closed formula of ESA robust against noise is used. Instantaneous Frequency (IF) and Instantaneous Amplitude (IA) estimates of a multi- component AM–FM signal, corrupted with additive white Gaussian noise of varying SNRs, are analyzed and results compared to ESA, DESA and Hilbert transform-based algorithms. SNR and MSE are used as figures of merit. Regularized BS version of EMD– ESA performs reasonably better in separating IA and IF components compared to the other methods from low to high SNR. Overall, obtained results illustrate the effective- ness of the proposed approach in terms of accuracy and robustness against noise to track IF and IA features of a multicomponent AM–FM signal

Une nouvelle classe d'opérateurs de Teager-Kaiser multidimensionnels basée sur les dérivées directionnelles d'ordre supérieur

This work aims at introducing some energy operators linked to Teager-Kaiser energy operator and its associated higher order versions and expand them to multidimensional signals. These operators are very useful for analyzing oscillatory signals with time-varying amplitude and frequency (AM-FM). We prove that gradient tensors combined with Kronecker powers allow to express these operators by directional derivatives along any n-D vector. In particular, we show that the construction of a large class of non linear operators for AM-FM multidimensional signal demodulation is possible. Also, a new scalar function using the directional derivative along a vector giving the ”sign” of the frequency components is introduced. An application of this model to local n-D AM-FM signal is presented and related demodulation error rates estimates. To show the effectiveness and the robustness of our method in term of envelope and frequency components extraction, results obtained on synthetic and real data are compared to multi-dimensional energy separation algorithm and to our recently introduced n-D operator

Measuring Instantaneous Frequency of Local Field Potential Oscillations using the Kalman Smoother

Rhythmic local field potentials (LFPs) arise from coordinated neural activity. Inference of neural function based on the properties of brain rhythms remains a challenging data analysis problem. Algorithms that characterize non-stationary rhythms with high temporal and spectral resolution may be useful for interpreting LFP activity on the timescales in which they are generated. We propose a Kalman smoother based dynamic autoregressive model for tracking the instantaneous frequency (iFreq) and frequency modulation (FM) of noisy and non-stationary sinusoids such as those found in LFP data. We verify the performance of our algorithm using simulated data with broad spectral content, and demonstrate its application using real data recorded from behavioral learning experiments. In analyses of ripple oscillations (100–250 Hz) recorded from the rodent hippocampus, our algorithm identified novel repetitive, short timescale frequency dynamics. Our results suggest that iFreq and FM may be useful measures for the quantification of small timescale LFP dynamics.National Institutes of Health (U.S.) (NIH/NIMH R01 MH59733)National Institutes of Health (U.S.) (NIH/NIHLB R01 HL084502)Massachusetts Institute of Technology (Henry E. Singleton Presidential Graduate Fellowship Award

Multiple multidimensional morse wavelets

This paper defines a set of operators that localize a radial image in space and radial frequency simultaneously. The eigenfunctions of the operator are determined and a nonseparable orthogonal set of radial wavelet functions are found. The eigenfunctions are optimally concentrated over a given region of radial space and scale space, defined via a triplet of parameters. Analytic forms for the energy concentration of the functions over the region are given. The radial function localization operator can be generalised to an operator localizing any L-2(R-2) function. It is demonstrated that the latter operator, given an appropriate choice of localization region, approximately has the same radial eigenfunctions as the radial operator. Based on a given radial wavelet function a quaternionic wavelet is defined that can extract the local orientation of discontinuous signals as well as amplitude, orientation and phase structure of locally oscillatory signals. The full set of quaternionic wavelet functions are component by component orthogonal; their statistical properties are tractable, and forms for the variability of the estimators of the local phase and orientation are given, as well as the local energy of the image. By averaging estimators across wavelets, a substantial reduction in the variance is achieved

Rehaussement du signal de parole par EMD et opérateur de Teager-Kaiser

The authors would like to thank Professor Mohamed Bahoura from Universite de Quebec a Rimouski for fruitful discussions on time adaptive thresholdingIn this paper a speech denoising strategy based on time adaptive thresholding of intrinsic modes functions (IMFs) of the signal, extracted by empirical mode decomposition (EMD), is introduced. The denoised signal is reconstructed by the superposition of its adaptive thresholded IMFs. Adaptive thresholds are estimated using the Teager–Kaiser energy operator (TKEO) of signal IMFs. More precisely, TKEO identifies the type of frame by expanding differences between speech and non-speech frames in each IMF. Based on the EMD, the proposed speech denoising scheme is a fully data-driven approach. The method is tested on speech signals with different noise levels and the results are compared to EMD-shrinkage and wavelet transform (WT) coupled with TKEO. Speech enhancement performance is evaluated using output signal to noise ratio (SNR) and perceptual evaluation of speech quality (PESQ) measure. Based on the analyzed speech signals, the proposed enhancement scheme performs better than WT-TKEO and EMD-shrinkage approaches in terms of output SNR and PESQ. The noise is greatly reduced using time-adaptive thresholding than universal thresholding. The study is limited to signals corrupted by additive white Gaussian noise

2-D Prony-Huang Transform: A New Tool for 2-D Spectral Analysis

This work proposes an extension of the 1-D Hilbert Huang transform for the analysis of images. The proposed method consists in (i) adaptively decomposing an image into oscillating parts called intrinsic mode functions (IMFs) using a mode decomposition procedure, and (ii) providing a local spectral analysis of the obtained IMFs in order to get the local amplitudes, frequencies, and orientations. For the decomposition step, we propose two robust 2-D mode decompositions based on non-smooth convex optimization: a "Genuine 2-D" approach, that constrains the local extrema of the IMFs, and a "Pseudo 2-D" approach, which constrains separately the extrema of lines, columns, and diagonals. The spectral analysis step is based on Prony annihilation property that is applied on small square patches of the IMFs. The resulting 2-D Prony-Huang transform is validated on simulated and real data.Comment: 24 pages, 7 figure

Enhancing Missing Data Imputation of Non-stationary Signals with Harmonic Decomposition

Dealing with time series with missing values, including those afflicted by low quality or over-saturation, presents a significant signal processing challenge. The task of recovering these missing values, known as imputation, has led to the development of several algorithms. However, we have observed that the efficacy of these algorithms tends to diminish when the time series exhibit non-stationary oscillatory behavior. In this paper, we introduce a novel algorithm, coined Harmonic Level Interpolation (HaLI), which enhances the performance of existing imputation algorithms for oscillatory time series. After running any chosen imputation algorithm, HaLI leverages the harmonic decomposition based on the adaptive nonharmonic model of the initial imputation to improve the imputation accuracy for oscillatory time series. Experimental assessments conducted on synthetic and real signals consistently highlight that HaLI enhances the performance of existing imputation algorithms. The algorithm is made publicly available as a readily employable Matlab code for other researchers to use