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

Phase Synchronization Operator for On-Chip Brain Functional Connectivity Computation

This paper presents an integer-based digital processor for the calculation of phase synchronization between two neural signals. It is based on the measurement of time periods between two consecutive minima. The simplicity of the approach allows for the use of elementary digital blocks, such as registers, counters, and adders. The processor, fabricated in a 0.18- μ m CMOS process, only occupies 0.05 mm 2 and consumes 15 nW from a 0.5 V supply voltage at a signal input rate of 1024 S/s. These low-area and low-power features make the proposed processor a valuable computing element in closed-loop neural prosthesis for the treatment of neural disorders, such as epilepsy, or for assessing the patterns of correlated activity in neural assemblies through the evaluation of functional connectivity maps.Ministerio de Economía y Competitividad TEC2016-80923-POffice of Naval Research (USA) N00014-19-1-215

Estimation de l'enveloppe et de la fréquence locales par les opérateurs de Teager-Kaiser en interférométrie en lumière blanche.

In this work, a new method for surface extraction in white light scanning interferometry (WLSI) is introduced. The proposed extraction scheme is based on the Teager-Kaiser energy operator and its extended versions. This non-linear class of operators is helpful to extract the local instantaneous envelope and frequency of any narrow band AM-FM signal. Namely, the combination of the envelope and frequency information, allows effective surface extraction by an iterative re-estimation of the phase in association with a new correlation technique, based on a recent TK crossenergy operator. Through the experiments, it is shown that the proposed method produces substantially effective results in term of surface extraction compared to the peak fringe scanning technique, the five step phase shifting algorithm and the continuous wavelet transform based method. In addition, the results obtained show the robustness of the proposed method to noise and to the fluctuations of the carrier frequency

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

Uncovering temporal regularity in atmospheric dynamics through Hilbert phase analysis

Uncovering meaningful regularities in complex oscillatory signals is a challenging problem with applications across a wide range of disciplines.Here, we present a novel approach, based on the Hilbert transform (HT).We show that temporal periodicity can be uncovered by averagingthe signal in a moving window of appropriated length,t, before applying the HT. As a case study, we investigate global gridded surface airtemperature (SAT) datasets. By analyzing the variation of the mean rotation period,T, of the Hilbert phase as a function oft, we discoverwell-de ned plateaus. In many geographical regions, the plateau corresponds to the expected 1-yr solar cycle; however, in regions where SATdynamics is highly irregular, the plateaus reveal non-trivial periodicities, which can be interpreted in terms of climatic phenomena such asEl Niño. In these regions, we also nd that Fourier analysis is unable to detect the periodicity that emerges whentincreases and graduallywashes out SAT variability. The values ofTobtained for di erentts are then given to a standard machine learning algorithm. The resultsdemonstrate that these features are informative and constitute a new approach for SAT time series classi cation. To support these results, weanalyze the synthetic time series generated with a simple model and con rm that our method extracts information that is fully consistent withour knowledge of the model that generates the data. Remarkably, the variation ofTwithtin the synthetic data is similar to that observed in thereal SAT data. This suggests that our model contains the basic mechanisms underlying the unveiled periodicities. Our results demonstrate thatHilbert analysis combined with temporal averaging is a powerful new tool for discovering hidden temporal regularity in complex oscillatorysignals.Peer ReviewedPostprint (author's final draft

Differential Phase Estimation with the SeaMARC II Bathymetric Sidescan Sonar System

A maximum-likelihood estimator is used to extract differential phase measurements from noisy seafloor echoes received at pairs of transducers mounted on either side of the SeaMARC II bathymetricsidescan sonar system. Carrier frequencies for each side are about 1 kHz apart, and echoes from a transmitted pulse 2 ms long are analyzed. For each side, phase difference sequences are derived from the full complex data consisting of base-banded and digitized quadrature components of the received echoes. With less bias and a lower variance, this method is shown to be more efficient than a uniform mean estimator. It also does not exhibit the angular or time ambiguities commonly found in the histogram method used in the SeaMARC II system. A figure for the estimation uncertainty of the phasedifference is presented, and results are obtained for both real and simulated data. Based on this error estimate and an empirical verification derived through coherent ping stacking, a single filter length of 100 ms is chosen for data processing application

Data driven optimal filtering for phase and frequency of noisy oscillations: application to vortex flowmetering

A new method for extracting the phase of oscillations from noisy time series is proposed. To obtain the phase, the signal is filtered in such a way that the filter output has minimal relative variation in the amplitude (MIRVA) over all filters with complex-valued impulse response. The argument of the filter output yields the phase. Implementation of the algorithm and interpretation of the result are discussed. We argue that the phase obtained by the proposed method has a low susceptibility to measurement noise and a low rate of artificial phase slips. The method is applied for the detection and classification of mode locking in vortex flowmeters. A novel measure for the strength of mode locking is proposed.Comment: 12 pages, 10 figure

Digital Filters for Instantaneous Frequency Estimation

This technical note is on digital filters for the high-fidelity estimation of a sinusoidal signal's frequency in the presence of additive noise. The complex noise is assumed to be white (i.e. uncorrelated) however it need not be Gaussian. The complex signal is assumed to be of (approximately) constant magnitude and (approximately) polynomial phase such as the chirps emitted by bats, whale songs, pulse-compression radars, and frequency-modulated (FM) radios, over sufficiently short timescales. Such digital signals may be found at the end of a sequence of analogue heterodyning (i.e. mixing and low-pass filtering), down to a bandwidth that is matched to an analogue-to-digital converter (ADC), followed by digital heterodyning and sample rate reduction (optional) to match the clock frequency of the processor. The spacing of the discrete frequency bins (in cycles per sample) produced by the Fast Fourier Transform (FFT) is equal to the reciprocal of the window length (in samples). However, a long FFT (for fine frequency resolution) has a high complexity and a long latency, which may be prohibitive in embedded closed-loop systems, and unnecessary when the channel only contains a single sinusoid. In such cases, and for signals of constant frequency, the conventional approach involves the (weighted) average of instantaneous phase differences. General, naive, optimal, and pragmatic (recursive), filtering solutions are discussed and analysed here using Monte-Carlo (MC) simulations.Comment: Added arXiv ID to header and fixed a few typo