Source author record

José Luis Romero

José Luis Romero appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2026arXiv

Model agnostic signal encoding by leaky integrate and fire, performance and uncertainty

Integrate-and-fire is a resource efficient time-encoding mechanism that summarizes into a signed spike train those time intervals where a signal's charge exceeds a certain threshold. We analyze the IF encoder in terms of a very general notion of approximate bandwidth, which is shared by most commonly-used signal models. This complements results on exact encoding that may be overly adapted to a particular signal model. We take into account, possibly for the first time, the effect of uncertainty in the exact location of the spikes (as may arise by decimation), uncertainty of integration leakage (as may arise in realistic manufacturing), and boundary effects inherent to finite periods of exposure to the measurement device. The analysis is done by means of a concrete bandwidth-based Ansatz that can also be useful to initialize more sophisticated model specific reconstruction algorithms, and uses the earth mover's (Wasserstein) distance to measure spike discrepancy.

preprint2022arXiv

Efficient computation of the zeros of the Bargmann transform under additive white noise

We study the computation of the zero set of the Bargmann transform of a signal contaminated with complex white noise, or, equivalently, the computation of the zeros of its short-time Fourier transform with Gaussian window. We introduce the adaptive minimal grid neighbors algorithm (AMN), a variant of a method that has recently appeared in the signal processing literature, and prove that with high probability it computes the desired zero set. More precisely, given samples of the Bargmann transform of a signal on a finite grid with spacing $δ$, AMN is shown to compute the desired zero set up to a factor of $δ$ in the Wasserstein error metric, with failure probability $O(δ^4 \log^2(1/δ))$. We also provide numerical tests and comparison with other algorithms.

preprint2022arXiv

Invertibility of frame operators on Besov-type decomposition spaces

We derive an extension of the Walnut-Daubechies criterion for the invertibility of frame operators. The criterion concerns general reproducing systems and Besov-type spaces. As an application, we conclude that $L^2$ frame expansions associated with smooth and fast-decaying reproducing systems on sufficiently fine lattices extend to Besov-type spaces. This simplifies and improves recent results on the existence of atomic decompositions, which only provide a particular dual reproducing system with suitable properties. In contrast, we conclude that the $L^2$ canonical frame expansions extend to many other function spaces, and, therefore, operations such as analyzing using the frame, thresholding the resulting coefficients, and then synthesizing using the canonical dual frame are bounded on these spaces.

preprint2022arXiv

Spectral-norm risk rates for multi-taper estimation of Gaussian processes

We consider the estimation of the covariance of a stationary Gaussian process on a multi-dimensional grid from observations taken on a general acquisition domain. We derive spectral-norm risk rates for multi-taper estimators. When applied to one dimensional acquisition intervals, these show that Thomson's classical multi-taper has optimal risk rates, as they match known benchmarks. We also extend existing lower risk bounds to multi-dimensional grids and conclude that multi-taper estimators associated with certain two-dimensional acquisition domains also have almost optimal risk rates.

preprint2022arXiv

The planar low temperature Coulomb gas: separation and equidistribution

We consider planar Coulomb systems consisting of a large number $n$ of repelling point charges in the low temperature regime, where the inverse temperature $β$ grows at least logarithmically in $n$ as $n \longrightarrow \infty$, i.e., $β\gtrsim \log n$. Under suitable conditions on an external potential we prove results to the effect that the gas is with high probability uniformly separated and equidistributed with respect to the corresponding equilibrium measure (in the given external field). Our results generalize earlier results about Fekete configurations, i.e., the case $β=\infty$. There are also several auxiliary results which could be of independent interest. For example, our method of proof of equidistribution (a variant of "Landau's method") works for general families of configurations which are uniformly separated and which satisfy certain sampling and interpolation inequalities.

preprint2022arXiv

Zeros of Gaussian Weyl-Heisenberg functions and hyperuniformity of charge

We study Gaussian random functions on the complex plane whose stochastics are invariant under the Weyl-Heisenberg group (twisted stationarity). The theory is modeled on translation invariant Gaussian entire functions, but allows for non-analytic examples, in which case winding numbers can be either positive or negative. We calculate the first intensity of zero sets of such functions, both when considered as points on the plane, or as charges according to their phase winding. In the latter case, charges are shown to be in a certain average equilibrium independently of the particular covariance structure (universal screening). We investigate the corresponding fluctuations, and show that in many cases they are suppressed at large scales (hyperuniformity). This means that universal screening is empirically observable at large scales. We also derive an asymptotic expression for the charge variance. As a main application, we obtain statistics for the zero sets of the short-time Fourier transform of complex white noise with general windows, and also prove the following uncertainty principle: the expected number of zeros per unit area is minimized, among all window functions, exactly by generalized Gaussians. Further applications include poly-entire functions such as covariant derivatives of Gaussian entire functions.

preprint2021arXiv

On dual molecules and convolution-dominated operators

We show that sampling or interpolation formulas in reproducing kernel Hilbert spaces can be obtained by reproducing kernels whose dual systems form molecules, ensuring that the size profile of a function is fully reflected by the size profile of its sampled values. The main tool is a local holomorphic calculus for convolution-dominated operators, valid for groups with possibly non-polynomial growth. Applied to the matrix coefficients of a group representation, our methods improve on classical results on atomic decompositions and bridge a gap between abstract and concrete methods.

preprint2021arXiv

Sampling the flow of a bandlimited function

We analyze the problem of reconstruction of a bandlimited function $f$ from the space-time samples of its states $f_t=ϕ_t\ast f$ resulting from the convolution with a kernel $ϕ_t$. It is well-known that, in natural phenomena, uniform space-time samples of $f$ are not sufficient to reconstruct $f$ in a stable way. To enable stable reconstruction, a space-time sampling with periodic nonuniformly spaced samples must be used as was shown by Lu and Vetterli. We show that the stability of reconstruction, as measured by a condition number, controls the maximal gap between the spacial samples. We provide a quantitative statement of this result. In addition, instead of irregular space-time samples, we show that uniform dynamical samples at sub-Nyquist spatial rate allow one to stably reconstruct the function $\widehat f$ away from certain, explicitly described blind spots. We also consider several classes of finite dimensional subsets of bandlimited functions in which the stable reconstruction is possible, even inside the blind spots. We obtain quantitative estimates for it using Remez-Turán type inequalities. En route, we obtain a Remez-Turán inequality for prolate spheroidal wave functions. To illustrate our results, we present some numerics and explicit estimates for the heat flow problem.

preprint2020arXiv

Balian-Low type theorems on homogeneous groups

We prove strict necessary density conditions for coherent frames and Riesz sequences on homogeneous groups. Let $N$ be a connected, simply connected nilpotent Lie group with a dilation structure (a homogeneous group) and let $(π, \mathcal{H}_π)$ be an irreducible, square-integrable representation modulo the center $Z(N)$ of $N$ on a Hilbert space $\mathcal{H}_π$ of formal dimension $d_π$. If $g \in \mathcal{H}_π$ is an integrable vector and the set $\{ π(λ)g : λ\in Λ\}$ for a discrete subset $Λ\subseteq N / Z(N)$ forms a frame for $\mathcal{H}_π$, then its density satisfies the strict inequality $D^-(Λ)> d_π$, where $D^-(Λ)$ is the lower Beurling density. An analogous density condition $D^+(Λ) < d_π$ holds for a Riesz sequence in $\mathcal{H}_π$ contained in the orbit of $(π, \mathcal{H}_π)$. The proof is based on a deformation theorem for coherent systems, a universality result for $p$-frames and $p$-Riesz sequences, some results from Banach space theory, and tools from the analysis on homogeneous groups.

preprint2020arXiv

Filtering with Wavelet Zeros and Gaussian Analytic Functions

We present the continuous wavelet transform (WT) of white Gaussian noise and establish a connection to the theory of Gaussian analytic functions. Based on this connection, we propose a methodology that detects components of a signal in white noise based on the distribution of the zeros of its continuous WT. To illustrate that the continuous theory can be employed in a discrete setting, we establish a uniform convergence result for the discretized continuous WT and apply the proposed method to a variety of acoustic signals.

preprint2020arXiv

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.

preprint2020arXiv

The Nyquist sampling rate for spiraling curves

We consider the problem of reconstructing a compactly supported function from samples of its Fourier transform taken along a spiral. We determine the Nyquist sampling rate in terms of the density of the spiral and show that below this rate spirals suffer from an approximate form of aliasing. This sets a limit to the amount of undersampling that compressible signals admit when sampled along spirals. More precisely, we derive a lower bound on the condition number for the reconstruction of functions of bounded variation, and for functions that are sparse in the Haar wavelet basis.

preprint2019arXiv

Strict density inequalities for sampling and interpolation in weighted spaces of holomorphic functions

Answering a question of Lindholm, we prove strict density inequalities for sampling and interpolation in Fock spaces of entire functions in several complex variables defined by a plurisubharmonic weight. In particular, these spaces do not admit a set that is simultaneously sampling and interpolating. To prove optimality of the density conditions, we construct sampling sets with a density arbitrarily close to the critical density. The techniques combine methods from several complex variables (estimates for $\bar \partial$) and the theory of localized frames in general reproducing kernel Hilbert spaces (with no analyticity assumed). The abstract results on Fekete points and deformation of frames may be of independent interest.

preprint2018arXiv

Sharp results on sampling with derivatives in shift-invariant spaces and multi-window Gabor frames

We study the problem of sampling with derivatives in shift-invariant spaces generated by totally-positive functions of Gaussian type or by the hyperbolic secant. We provide sharp conditions in terms of weighted Beurling densities. As a by-product we derive new results about multi-window Gabor frames with respect to vectors of Hermite functions or totally positive functions.

preprint2017arXiv

Computing reconstructions from nonuniform Fourier samples: Universality of stability barriers and stable sampling rates

We study the problem of recovering an unknown compactly-supported multivariate function from samples of its Fourier transform that are acquired nonuniformly, i.e. not necessarily on a uniform Cartesian grid. Reconstruction problems of this kind arise in various imaging applications, where Fourier samples are taken along radial lines or spirals for example. Specifically, we consider finite-dimensional reconstructions, where a limited number of samples is available, and investigate the rate of convergence of such approximate solutions and their numerical stability. We show that the proportion of Fourier samples that allow for stable approximations of a given numerical accuracy is independent of the specific sampling geometry and is therefore universal for different sampling scenarios. This allows us to relate both sufficient and necessary conditions for different sampling setups and to exploit several results that were previously available only for very specific sampling geometries. The results are obtained by developing: (i) a transference argument for different measures of the concentration of the Fourier transform and Fourier samples; (ii) frame bounds valid up to the critical sampling density, which depend explicitly on the sampling set and the spectrum. As an application, we identify sufficient and necessary conditions for stable and accurate reconstruction of algebraic polynomials or wavelet coefficients from nonuniform Fourier data.

preprint2016arXiv

Completeness of Gabor systems

We investigate the completeness of Gabor systems with respect to several classes of window functions on rational lattices. Our main results show that the time-frequency shifts of every finite linear combination of Hermite functions with respect to a rational lattice are complete in L2(R), thus generalizing a remark of von Neumann (and proved by Bargmann, Perelomov et al.). An analogous result is proven for functions that factor into certain rational functions and the Gaussian. The results are also interesting from a conceptual point of view since they show a vast difference between the completeness and the frame property of a Gabor system. In the terminology of physics we prove new results about the completeness of coherent state subsystems.

preprint2016arXiv

Stability of Gabor frames under small time Hamiltonian evolutions

We consider Hamiltonian deformations of Gabor systems, where the window evolves according to the action of a Schrödinger propagator and the phase-space nodes evolve according to the corresponding Hamiltonian flow. We prove the stability of the frame property for small times and Hamiltonians consisting of a quadratic polynomial plus a potential in the Sjöstrand class with bounded second order derivatives. This answers a question raised in [de Gosson, M. Symplectic and Hamiltonian Deformations of Gabor Frames. Appl. Comput. Harmon. Anal. Vol. 38 No.2, (2015) p.196--221.]

preprint2016arXiv

The bias-variance trade-off in Thomson's multitaper estimator

At the heart of non-parametric spectral estimation, lies the dilemma known as the bias-variance trade-off: low biased estimators tend to have high variance and low variance estimators tend to have high bias. In 1982, Thomson introduced a multitaper method where this trade-off is made explicit by choosing a target bias resolution and obtaining a corresponding variance reduction. The method became the standard in many applications. Its favorable bias-variance trade-off is due to an empirical fact, conjectured by Thomson based on numerical evidence: assuming bandwidth W and N time domain observations, the average of the square of the first $K=\left\lfloor 2NW\right\rfloor$ Slepian functions approaches, as K grows, an ideal band-pass kernel for the interval [-W,W]. We provide an analytic proof of this fact and quantify the approximation error in the L1 norm; the approximation error is then used to control the bias of the multitaper estimator resulting from spectral leakage. This leads to new performance bounds for the method, explicit in terms of the bandwidth W and the number N of time domain observations. Our method is flexible and can be extended to higher dimensions and different geometries.

preprint2015arXiv

Deformation of Gabor systems

We introduce a new notion for the deformation of Gabor systems. Such deformations are in general nonlinear and, in particular, include the standard jitter error and linear deformations of phase space. With this new notion we prove a strong deformation result for Gabor frames and Gabor Riesz sequences that covers the known perturbation and deformation results. Our proof of the deformation theorem requires a new characterization of Gabor frames and Gabor Riesz sequences. It is in the style of Beurling's characterization of sets of sampling for bandlimited functions and extends significantly the known characterization of Gabor frames "without inequalities" from lattices to non-uniform sets.

preprint2014arXiv

On accumulated spectrograms

We study the eigenvalues and eigenfunctions of the time-frequency localization operator $H_Ω$ on a domain $Ω$ of the time-frequency plane. The eigenfunctions are the appropriate prolate spheroidal functions for an arbitrary domain $Ω$. Indeed, in analogy to the classical theory of Landau-Slepian-Pollak, the number of eigenvalues of $H_Ω$ in $[1-δ, 1]$ is equal to the measure of $Ω$ up to an error term depending on the perimeter of the boundary of $Ω$. Our main results show that the spectrograms of the eigenfunctions corresponding to the large eigenvalues (which we call the accumulated spectrogram) form an approximate partition ofunity of the given domain $Ω$. We derive both asymptotic, non-asymptotic, and weak $L^2$ error estimates for the accumulated spectrogram. As a consequence the domain $Ω$ can be approximated solely from the spectrograms of eigenfunctions without information about their phase.

preprint2014arXiv

On Minimal Trajectories for Mobile Sampling of Bandlimited Fields

We study the design of sampling trajectories for stable sampling and the reconstruction of bandlimited spatial fields using mobile sensors. The spectrum is assumed to be a symmetric convex set. As a performance metric we use the path density of the set of sampling trajectories that is defined as the total distance traveled by the moving sensors per unit spatial volume of the spatial region being monitored. Focussing first on parallel lines, we identify the set of parallel lines with minimal path density that contains a set of stable sampling for fields bandlimited to a known set. We then show that the problem becomes ill-posed when the optimization is performed over all trajectories by demonstrating a feasible trajectory set with arbitrarily low path density. However, the problem becomes well-posed if we explicitly specify the stability margins. We demonstrate this by obtaining a non-trivial lower bound on the path density of an arbitrary set of trajectories that contain a sampling set with explicitly specified stability bounds.

preprint2013arXiv

Multi-window Gabor frames in amalgam spaces

We show that multi-window Gabor frames with windows in the Wiener algebra $W(L^{\infty}, \ell^{1})$ are Banach frames for all Wiener amalgam spaces. As a byproduct of our results we positively answer an open question that was posed by [Krishtal and Okoudjou, Invertibility of the Gabor frame operator on the Wiener amalgam space, J. Approx. Theory, 153(2), 2008] and concerns the continuity of the canonical dual of a Gabor frame with a continuous generator in the Wiener algebra. The proofs are based on a recent version of Wiener's $1/f$ lemma.

preprint2012arXiv

Non-uniform painless decompositions for anisotropic Besov and Triebel-Lizorkin spaces

In this article we construct affine systems that provide a simultaneous atomic decomposition for a wide class of functional spaces including the Lebesgue spaces $L^p(\Rdst)$, $1<p<+\infty$. The novelty and difficulty of this construction is that we allow for non-lattice translations. We prove that for an arbitrary expansive matrix $A$ and any set $Λ$ - satisfying a certain spreadness condition but otherwise irregular- there exists a smooth window whose translations along the elements of $Λ$ and dilations by powers of $A$ provide an atomic decomposition for the whole range of the anisotropic Triebel-Lizorkin spaces. The generating window can be either chosen to be bandlimited or to have compact support. To derive these results we start with a known general "painless" construction that has recently appeared in the literature. We show that this construction extends to Besov and Triebel-Lizorkin spaces by providing adequate dual systems.

preprint2011arXiv

Characterization of coorbit spaces with phase-space covers

We show that coorbit spaces can be characterized in terms of arbitrary phase-space covers, which are families of phase-space multipliers associated with partitions of unity. This generalizes previously known results for time-frequency analysis to include time-scale decompositions. As a by-product, we extend the existing results for time-frequency analysis to an irregular setting.

preprint2009arXiv

Approximate reconstruction of bandlimited functions for the integrate and fire sampler

In this paper we study the reconstruction of a bandlimited signal from samples generated by the integrate and fire model. This sampler allows us to trade complexity in the reconstruction algorithms for simple hardware implementations, and is specially convenient in situations where the sampling device is limited in terms of power, area and bandwidth. Although perfect reconstruction for this sampler is impossible, we give a general approximate reconstruction procedure and bound the corresponding error. We also show the performance of the proposed algorithm through numerical simulations.