Search

We present a construction of a multi-scale Gaussian beam parametrix for the Dirichlet boundary value problem associated with the wave equation, and study its convergence rate to the true solution in the highly oscillatory regime. The construction elaborates on the wave-atom parametrix of Bao, Qian, Ying, and Zhang and extends to a multi-scale setting the technique of Gaussian beam propagation from a boundary of Katchalov, Kurylev and Lassas., Comment: 64 pages, 7 figures, minor update

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.

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 (π,Hπ) be an irreducible, square-integrable representation modulo the center Z(N) of N on a Hilbert space Hπ of formal dimension dπ. If g ∈ Hπ is an integrable vector and the set {π(λ)g : λ ∈ Λ} for a discrete subset Λ ⊆ N/Z(N) forms a frame for 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 Hπ contained in the orbit of (π,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.

We show that a real-valued function $f$ in the shift-invariant space generated by a totally positive function of Gaussian type is uniquely determined, up to a sign, by its absolute values $\{|f(\lambda)|: \lambda \in \Lambda \}$ on any set $\Lambda \subseteq \mathbb{R}$ with lower Beurling density $D^{-}(\Lambda)>2$., Comment: 7 pages

We characterize sampling and interpolating sets with derivatives in weighted Fock spaces on the complex plane in terms of their weighted Beurling densities., Comment: 27 pages. Some minor typographical errors were corrected

We analyze the problem of reconstruction of a bandlimited function $f$ from the space-time samples of its states $f_t=\phi_t\ast f$ resulting from the convolution with a kernel $\phi_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\'an type inequalities. En route, we obtain a Remez-Tur\'an inequality for prolate spheroidal wave functions. To illustrate our results, we present some numerics and explicit estimates for the heat flow problem., Comment: 29 pages

We obtain estimates for the Mean Squared Error (MSE) for the multitaper spectral estimator and certain compressive acquisition methods for multi-band signals. We confirm a fact discovered by Thomson [Spectrum estimation and harmonic analysis, Proc. IEEE, 1982]: assuming bandwidth $W$ and $N$ time domain observations, the average of the square of the first $K=2NW$ 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 measure the corresponding rate of convergence in the $L^{1}$ norm. This validates a heuristic approximation used to control the MSE of the multitaper estimator. The estimates have also consequences for the method of compressive acquisition of multi-band signals introduced by Davenport and Wakin, giving MSE approximation bounds for the dictionary formed by modulation of the critical number of prolates., Comment: 16 pages, 2 figures. (This article replaces arXiv: 1503.02991.)

We derive necessary density conditions for sampling and for interpolation in general reproducing kernel Hilbert spaces satisfying some natural conditions on the geometry of the space and the reproducing kernel. If the volume of shells is small compared to the volume of balls (weak annular decay property) and if the kernel possesses some off-diagonal decay or even some weaker form of localization, then there exists a critical density D with the following property: a set of sampling has density ⩾D, whereas a set of interpolation has density ⩽D. The main theorem unifies many known density theorems in signal processing, complex analysis, and harmonic analysis. For the special case of bandlimited function we recover Landau's fundamental density result. In complex analysis we rederive a critical density for generalized Fock spaces. In harmonic analysis we obtain the first general result about the density of coherent frames.

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., 33 pages

International audience; We present our recent work on sampling along spiral-like curves [9], and discuss the main techniques. As a first result we give a sharp density condition for sampling on spirals in terms of the separation between consecutive branches. We then further show that, below this rate, the numerical stability related to the reconstruction of compressible signals when sampled along spirals is significantly limited by the amount of undersampling.

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.

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., 24 pages

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., 31 pages, 2 figures

We present our recent work on the Weyl-Heisenberg ensemble and its statistical properties [4]. The WH ensemble is a class of determinantal point processes associated with the Schrodinger representation of the Heisenberg group. As a special example, WH ensembles include a multi-layer extension of the Ginibre ensemble modeling the distribution of electrons in higher Landau levels. We describe the hyperuniformity of WH ensembles, which characterizes a state of matter for which (scaled) density fluctuations diminish towards zero at the largest length scales. Our approach is based on methods from time-frequency analysis. We introduce the main results from [4] highlighting time-frequency techniques and connections to the theory of polyanalytic functions, and also present some small extensions.

We present our recent work on Mean Squared Error (MSE) bounds for Thomson's multitaper spectral estimator [5]. The main result provides a description of the average of the critical number of squared Slepian functions. As a further application we obtain MSE bounds for the approximation of multi-band signals using the dictionary of modulated Slepian sequences introduced by Davenport and Wakin.

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., 17 pages

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., Comment: 22 pages

We study nonuniform sampling in shift-invariant spaces and the construction of Gabor frames with respect to the class of totally positive functions whose Fourier transform factors as $ \hat g(\xi)= \prod_{j=1}^n (1+2\pi i\delta_j\xi)^{-1} \, e^{-c \xi^2}$ for $\delta_1,\ldots,\delta_n\in \mathbb{R}, c >0$ (in which case $g$ is called totally positive of Gaussian type). In analogy to Beurling's sampling theorem for the Paley-Wiener space of entire functions, we prove that every separated set with lower Beurling density $>1$ is a sampling set for the shift-invariant space generated by such a $g$. In view of the known necessary density conditions, this result is optimal and validates the heuristic reasonings in the engineering literature. Using a subtle connection between sampling in shift-invariant spaces and the theory of Gabor frames, we show that the set of phase-space shifts of $g$ with respect to a rectangular lattice $\alpha \mathbb{Z} \times \beta \mathbb{Z}$ forms a frame, if and only if $\alpha \beta, Comment: 25 pages

We investigate an inverse problem in time-frequency localization: the approximation of the symbol of a time-frequency localization operator from partial spectral information by the method of accumulated spectrograms (the sum of the spectrograms corresponding to large eigenvalues). We derive a sharp bound for the rate of convergence of the accumulated spectrogram, improving on recent results., Comment: 14 pages, 2 figures

Weyl-Heisenberg ensembles are a class of determinantal point processes associated with the Schr\"odinger representation of the Heisenberg group. Hyperuniformity characterizes a state of matter for which (scaled) density fluctuations diminish towards zero at the largest length scales. We will prove that Weyl-Heisenberg ensembles are hyperuniform. Weyl-Heisenberg ensembles include as a special case a multi-layer extension of the Ginibre ensemble modeling the distribution of electrons in higher Landau levels, which has recently been object of study in the realm of the Ginibre-type ensembles associated with polyanalytic functions. In addition, the family of Weyl-Heisenberg ensembles includes new structurally anisotropic processes, where point-statistics depend on the different spatial directions, and thus provide a first means to study directional hyperuniformity., Comment: 17 pages, 2 figures. (Typos corrected.)

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., 28 pages, 8 figures

We consider Hamiltonian deformations of Gabor systems, where the window evolves according to the action of a Schr\"odinger 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\"ostrand 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.], Comment: 11 pages. Minor revision