Comparison of Matching Pursuit Algorithm with Other Signal Processing Techniques for Computation of the Time-Frequency Power Spectrum of Brain Signals

Signals recorded from the brain often show rhythmic patterns at different frequencies, which are tightly coupled to the external stimuli as well as the internal state of the subject. In addition, these signals have very transient structures related to spiking or sudden onset of a stimulus, which have durations not exceeding tens of milliseconds. Further, brain signals are highly nonstationary because both behavioral state and external stimuli can change on a short time scale. It is therefore essential to study brain signals using techniques that can represent both rhythmic and transient components of the signal, something not always possible using standard signal processing techniques such as short time fourier transform, multitaper method, wavelet transform, or Hilbert transform. In this review, we describe a multiscale decomposition technique based on an over-complete dictionary called matching pursuit (MP), and show that it is able to capture both a sharp stimulus-onset transient and a sustained gamma rhythm in local field potential recorded from the primary visual cortex. We compare the performance of MP with other techniques and discuss its advantages and limitations. Data and codes for generating all time-frequency power spectra are provided.

Universality and scaling behavior of injected power in elastic turbulence in wormlike micellar gel

We study the statistical properties of spatially averaged global injected power fluctuations for Taylor-Couette flow of a wormlike micellar gel formed by surfactant cetyltrimethylammonium tosylate. At sufficiently high Weissenberg numbers the shear rate, and hence the injected power p(t), at a constant applied stress shows large irregular fluctuations in time. The nature of the probability distribution function (PDF) of p(t) and the power-law decay of its power spectrum are very similar to that observed in recent studies of elastic turbulence for polymer solutions. Remarkably, these non-Gaussian PDFs can be well described by a universal, large deviation functional form given by the generalized Gumbel distribution observed in the context of spatially averaged global measures in diverse classes of highly correlated systems. We show by in situ rheology and polarized light scattering experiments that in the elastic turbulent regime the flow is spatially smooth but random in time, in agreement with a recent hypothesis for elastic turbulence.

Short-window spectral analysis using AMVAR and multitaper methods: a comparison

We compare two popular methods for estimating the power spectrum from short data windows, namely the adaptive multivariate autoregressive (AMVAR) method and the multitaper method. By analyzing a simulated signal (embedded in a background Ornstein-Uhlenbeck noise process) we demonstrate that the AMVAR method performs better at detecting short bursts of oscillations compared to the multitaper method. However, both methods are immune to jitter in the temporal location of the signal. We also show that coherence can still be detected in noisy bivariate time series data by the AMVAR method even if the individual power spectra fail to show any peaks. Finally, using data from two monkeys performing a visuomotor pattern discrimination task, we demonstrate that the AMVAR method is better able to determine the termination of the beta oscillations when compared to the multitaper method.

Quantum diffusion in thin disordered wires

Quantum Ohmic residual resistance of a thin disordered wire, approximated as a one-dimensional multichannel conductor, is known to scale exponentially with length. This nonadditivity is shown to imply (i) a low-frequency noise-power spectrum proportional to -ln(Ω)/Ω, and (ii) a dispersive capacitative impedance proportional to tanh(√iΩ )/ √iΩ. A deep connection to the quantum Brownian motion with linear dynamical frictional coupling to a harmonic-oscillator bath is pointed out and interpreted in physical terms.

Non-diffusive classical transport in a medium with static randomness

The probability distribution for the displacement x of a particle moving in a one-dimensional continuum is derived exactly for the general case of combined static and dynamic gaussian randomness of the applied force. The dynamics of the particle is governed by the high-friction limit of Brownian motion discussed originally by Einstein and Smoluchowski. In particular, the mean square displacement of the particle varies as t2 for t to infinity . This ballistic motion induced by the disorder does not give rise to a 1/f power spectrum, contrary to recent suggestions based on the above dynamical model.

A parameter characterizing plowing nature of surfaces close to Gaussian

A non-dimensional parameter descriptive of the plowing nature of surfaces is proposed for the case of sliding between a soft and a relatively hard metallic pair. From a set of potential parameters which can be descriptive of the phenomenon, dimensionless groups are formulated and the influence of each one of them is analyzed. A non-dimensional parameter involving the root-mean square deviation (R-q) and the centroidal frequency (F-mean) deducted from the power-spectrum is found to have a high degree of correlation (as high as 0.93) with the coefficient of friction obtained in sliding experiments under lubricated condition.

Perturbative growth of cosmological clustering .2. The two-point correlation

We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest order nonlinear correction to the two-point correlation and the pair velocity for Gaussian initial conditions in a critical density matter-dominated cosmological model. We compare our results with the results obtained using the hydrodynamic equations that neglect pressure and find that the two match, indicating that there are no effects of multistreaming at this order of perturbation. We analytically study the effect of small scales on the large scales by calculating the nonlinear correction for a Dirac delta function initial two-point correlation. We find that the induced two-point correlation has a x(-6) behavior at large separations. We have considered a class of initial conditions where the initial power spectrum at small k has the form k(n) with 0 < n less than or equal to 3 and have numerically calculated the nonlinear correction to the two-point correlation, its average over a sphere and the pair velocity over a large dynamical range. We find that at small separations the effect of the nonlinear term is to enhance the clustering, whereas at intermediate scales it can act to either increase or decrease the clustering. At large scales we find a simple formula that gives a very good fit for the nonlinear correction in terms of the initial function. This formula explicitly exhibits the influence of small scales on large scales and because of this coupling the perturbative treatment breaks down at large scales much before one would expect it to if the nonlinearity were local in real space. We physically interpret this formula in terms of a simple diffusion process. We have also investigated the case n = 0, and we find that it differs from the other cases in certain respects. We investigate a recently proposed scaling property of gravitational clustering, and we find that the lowest order nonlinear terms cause deviations from the scaling relations that are strictly valid in the linear regime. The approximate validity of these relations in the nonlinear regime in l(T)-body simulations cannot be understood at this order of evolution.

Resent radar observations of asteroid 1566 Icarus

We report Doppler-only radar observations of Icarus at Goldstone at a transmitter frequency of 8510 MHz (3.5 cm wavelength) during 8-10 June 1996, the first radar detection of the object since 1968. Optimally filtered and folded spectra achieve a maximum opposite-circular (OC) polarization signal-to-noise ratio of about 10 and help to constrain Icarus' physical properties. We obtain an OC radar cross section of 0.05 km(2) (with a 35% uncertainty), which is less than values estimated by Goldstein (1969) and by Pettengill et al. (1969), and a circular polarization (SC/OC) ratio of 0.5+/-0.2. We analyze the echo power spectrum with a model incorporating the echo bandwidth B and a spectral shape parameter it, yielding a coupled constraint between B and n. We adopt 25 Hz as the lower bound on B, which gives a lower bound on the maximum pole-on breadth of about 0.6 km and upper bounds on the radar and optical albedos that are consistent with Icarus' tentative QS classification. The observed circular polarization ratio indicates a very rough near-surface at spatial scales of the order of the radar wavelength. (C) 1999 Elsevier Science Ltd. All rights reserved.

Subdiffusion in hair bundle dynamics: The role of protein conformational fluctuations

The detection of sound signals in vertebrates involves a complex network of different mechano-sensory elements in the inner ear. An especially important element in this network is the hair bundle, an antenna-like array of stereocilia containing gated ion channels that operate under the control of one or more adaptation motors. Deflections of the hair bundle by sound vibrations or thermal fluctuations transiently open the ion channels, allowing the flow of ions through them, and producing an electrical signal in the process, eventually causing the sensation of hearing. Recent high frequency (0.1-10 kHz) measurements by Kozlov et al. Proc. Natl. Acad. Sci. U. S. A. 109, 2896 (2012)] of the power spectrum and the mean square displacement of the thermal fluctuations of the hair bundle suggest that in this regime the dynamics of the hair bundle are subdiffusive. This finding has been explained in terms of the simple Brownian motion of a filament connecting neighboring stereocilia (the tip link), which is modeled as a viscoelastic spring. In the present paper, the diffusive anomalies of the hair bundle are ascribed to tip link fluctuations that evolve by fractional Brownian motion, which originates in fractional Gaussian noise and is characterized by a power law memory. The predictions of this model for the power spectrum of the hair bundle and its mean square displacement are consistent with the experimental data and the known properties of the tip link. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4768902]

The bystander effect in optically trapped red blood cells due to Plasmodium falciparum infection

In a previous study of the properties of red blood cells (RBC) trapped in an optical tweezers trap, an increase in the spectrum of Brownian fluctuations for RBCs from a Plasmodium falciparum culture (due to increased rigidity) compared with normal RBCs was measured. A bystander effect was observed, whereby RBCs actually hosting the parasite had an effect on the physical properties of remaining non-hosting RBCs. The distribution of corner frequency (f(c)) in the power spectrum of single RBCs held in an optical tweezers trap was studied. Two tests were done to confirm the bystander effect. In the first, RBCs from an infected culture were separated into hosting and non-hosting RBCs. In the second, all RBCs were removed from the infected culture, and normal RBCs were incubated in the spent medium. The trapping environment was the same for all measurements so only changes in the properties of RBCs were measured. In the first experiment, a similar and statistically significant increase was measured both for hosting and non-hosting RBCs. In the second experiment, normal RBCs incubated in spent medium started to become rigid after a few hours and showed complete changes (comparable with RBCs from the infected culture) after 24 h. These experiments provide direct evidence of medium-induced changes in the properties of RBCs in an infected culture, regardless of whether the RBCs actually host the parasite.

Spectral analysis of Cauvery river flows

A detalied study of the maonthly Convery river flows at the krishna raja sagara (KRS) reservoir is carried out by using the techniques of spectral analysis. The correlogram and power spectrum ate platted and used to identify the peridiocities inherent in the monthly inflows. The statistical significance of these periodicities is tested. Apart from the periodiocities at 12 months and 6 months, a significant of periodicity of 4 month was also observed in the monthly inflows. The analysis prepares ground for developing an appropriate stochastic model for the item series of the monthly flows.

Structural characterization of electrodeposited boron

Structural characterization of electrodeposited boron was carried out by using transmission electron microscopy and Raman spectroscopy. Electron diffraction and phase contrast imaging were carried out by using transmission electron microscopy. Phase identification was done based on the analysis of electron diffraction patterns and the power spectrum calculated from the lattice images from thin regions of the sample. Raman spectroscopic examination was carried out to study the nature of bonding and the allotropic form of boron obtained after electrodeposition. The results obtained from transmission electron microscopy showed the presence of nanocrystallites embedded in an amorphous mass of boron. Raman microscopic studies showed that amorphous boron could be converted to its crystalline form at high temperatures.

On the duality between rate and power optimizations

Sequence design problems are considered in this paper. The problem of sum power minimization in a spread spectrum system can be reduced to the problem of sum capacity maximization, and vice versa. A solution to one of the problems yields a solution to the other. Subsequently, conceptually simple sequence design algorithms known to hold for the white-noise case are extended to the colored noise case. The algorithms yield an upper bound of 2N - L on the number of sequences where N is the processing gain and L the number of non-interfering subsets of users. If some users (at most N - 1) are allowed to signal along a limited number of multiple dimensions, then N orthogonal sequences suffice.

Cooperative Spectrum Sensing Algorithms For OFDM Systems with Frequency Selective Channels

This paper considers the problem of spectrum sensing in cognitive radio networks when the primary user employs Orthogonal Frequency Division Multiplexing (OFDM). We specifically consider the scenario when the channel between the primary and a secondary user is frequency selective. We develop cooperative sequential detection algorithms based on energy detectors. We modify the detectors to mitigate the effects of some common model uncertainties such as timing and frequency offset, IQ-imbalance and uncertainty in noise and transmit power. The performance of the proposed algorithms are studied via simulations. We show that the performance of the energy detector is not affected by the frequency selective channel. We also provide a theoretical analysis for some of our algorithms.

Bayesian Decentralized Spectrum Sensing in Cognitive Radio Networks

This paper considers the problem of spectrum sensing, i.e., the detection of whether or not a primary user is transmitting data by a cognitive radio. The Bayesian framework is adopted, with the performance measure being the probability of detection error. A decentralized setup, where N sensors use M observations each to arrive at individual decisions that are combined at a fusion center to form the overall decision is considered. The unknown fading channel between the primary sensor and the cognitive radios makes the individual decision rule computationally complex, hence, a generalized likelihood ratio test (GLRT)-based approach is adopted. Analysis of the probabilities of false alarm and miss detection of the proposed method reveals that the error exponent with respect to M is zero. Also, the fusion of N individual decisions offers a diversity advantage, similar to diversity reception in communication systems, and a tight bound on the error exponent is presented. Through an analysis in the low power regime, the number of observations needed as a function of received power, to achieve a given probability of error is determined. Monte-Carlo simulations confirm the accuracy of the analysis.