462 research outputs found
Multitaper estimation on arbitrary domains
Multitaper estimators have enjoyed significant success in estimating spectral
densities from finite samples using as tapers Slepian functions defined on the
acquisition domain. Unfortunately, the numerical calculation of these Slepian
tapers is only tractable for certain symmetric domains, such as rectangles or
disks. In addition, no performance bounds are currently available for the mean
squared error of the spectral density estimate. This situation is inadequate
for applications such as cryo-electron microscopy, where noise models must be
estimated from irregular domains with small sample sizes. We show that the
multitaper estimator only depends on the linear space spanned by the tapers. As
a result, Slepian tapers may be replaced by proxy tapers spanning the same
subspace (validating the common practice of using partially converged solutions
to the Slepian eigenproblem as tapers). These proxies may consequently be
calculated using standard numerical algorithms for block diagonalization. We
also prove a set of performance bounds for multitaper estimators on arbitrary
domains. The method is demonstrated on synthetic and experimental datasets from
cryo-electron microscopy, where it reduces mean squared error by a factor of
two or more compared to traditional methods.Comment: 28 pages, 11 figure
Orthogonal sets of data windows constructed from trigonometric polynomials
Suboptimal, easily computable substitutes for the discrete prolate-spheroidal windows used by Thomson for spectral estimation are given. Trigonometric coefficients and energy leakages of the window polynomials are tabulated
Compressive Sensing of Analog Signals Using Discrete Prolate Spheroidal Sequences
Compressive sensing (CS) has recently emerged as a framework for efficiently
capturing signals that are sparse or compressible in an appropriate basis.
While often motivated as an alternative to Nyquist-rate sampling, there remains
a gap between the discrete, finite-dimensional CS framework and the problem of
acquiring a continuous-time signal. In this paper, we attempt to bridge this
gap by exploiting the Discrete Prolate Spheroidal Sequences (DPSS's), a
collection of functions that trace back to the seminal work by Slepian, Landau,
and Pollack on the effects of time-limiting and bandlimiting operations. DPSS's
form a highly efficient basis for sampled bandlimited functions; by modulating
and merging DPSS bases, we obtain a dictionary that offers high-quality sparse
approximations for most sampled multiband signals. This multiband modulated
DPSS dictionary can be readily incorporated into the CS framework. We provide
theoretical guarantees and practical insight into the use of this dictionary
for recovery of sampled multiband signals from compressive measurements
Augmented Slepians: Bandlimited Functions that Counterbalance Energy in Selected Intervals
Slepian functions provide a solution to the optimization problem of joint
time-frequency localization. Here, this concept is extended by using a
generalized optimization criterion that favors energy concentration in one
interval while penalizing energy in another interval, leading to the
"augmented" Slepian functions. Mathematical foundations together with examples
are presented in order to illustrate the most interesting properties that these
generalized Slepian functions show. Also the relevance of this novel
energy-concentration criterion is discussed along with some of its
applications
- …