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.     

On optimal wavelet reconstructions from Fourier samples: Linearity and universality of the stable sampling rate

In this paper we study the problem of computing wavelet coefficients of compactly supported functions from their Fourier samples. For this, we use the recently introduced framework of generalized sampling. Our first result demonstrates that using generalized sampling one obtains a stable and accurate reconstruction, provided the number of Fourier samples grows linearly in the number of wavelet coefficients recovered. For the class of Daubechies wavelets we derive the exact constant of proportionality.Our second result concerns the optimality of generalized sampling for this problem. Under some mild assumptions we show that generalized sampling cannot be outperformed in terms of approximation quality by more than a constant factor. Moreover, for the class of so-called perfect methods, any attempt to lower the sampling ratio below a certain critical threshold necessarily results in exponential ill-conditioning. Thus generalized sampling provides a nearly-optimal solution to this problem.     

Sampling and recovery of multidimensional bandlimited functions via frames

In this paper, we investigate frames for L₂d[−π,π] consisting of exponential functions in connection to oversampling and nonuniform sampling of bandlimited functions. We derive a multidimensional nonuniform oversampling formula for bandlimited functions with a fairly general frequency domain. The stability of said formula under various perturbations in the sampled data is investigated, and a computationally manageable simplification of the main oversampling theorem is given. Also, a generalization of Kadec's 1/4 theorem to higher dimensions is considered. Finally, the developed techniques are used to approximate biorthogonal functions of particular exponential Riesz bases for L₂[−π,π], and a well-known theorem of Levinson is recovered as a corollary.     

Wavelet frames and shift-invariant subspaces of periodic functions

A general approach based on polyphase splines, with analysis in the frequency domain, is developed for studying wavelet frames of periodic functions of one or higher dimensions. Characterizations of frames for shift-invariant subspaces of periodic functions and results on the structure of these subspaces are obtained. Starting from any multiresolution analysis, a constructive proof is provided for the existence of a normalized tight wavelet frame. The construction gives the minimum number of wavelets required. As an illustration of the approach developed, the one-dimensional dyadic case is further discussed in detail, concluding with a concrete example of trigonometric polynomial wavelet frames.     

On summability of Fourier coefficients of functions from Lebesgue space

In this paper we state new theorems of Hardy–Littlewood type for functions with general monotone Fourier coefficients. Sharpness of stated results is discusses.     

Fourier and Hermite series estimates of regression functions

Summary In the paper we estimate a regressionm(x)=E {Y|X=x} from a sequence of independent observations (X ₁,Y ₁),…, (X _n, Y_n) of a pair (X, Y) of random variables. We examine an estimate of a type , whereN depends onn andϕ _N is Dirichlet kernel and the kernel associated with the hermite series. Assuming, that E|Y|<∞ and |Y|≦γ≦∞, we give condition for to converge tom(x) at almost allx, provided thatX has a density. if the regression hass derivatives, then converges tom(x) as rapidly asO(nC^{−(2s−1)/4s}) in probability andO(n ^{−(2s−1)/4s} logn) almost completely.     

Fourier coefficients of Zygmund functions and analytic functions with quasiconformal deformation extensions

An open problem is to characterize the Fourier coefficients of Zygmund functions.This problem was also explicitly suggested by Nag and later by Teo and Takhtajan-Teo in the course of study of the universal Teichmu¨ller space.By a complex analysis approach,we give a characterization for the Fourier coefficients of a Zygmund function by a quadratic form.Some related topics are also discussed,including those analytic functions with quasiconformal deformation extensions.     

Convergence of double Fourier series of functions from symmetric spaces

We establish conditions for the convergence of double Fourier series in the trigonometric system of functions belonging to a symmetric space.     

Quantization of compressive samples with stable and robust recovery

In this paper we study the quantization stage that is implicit in any compressed sensing signal acquisition paradigm. We propose using Sigma–Delta (ΣΔ) quantization and a subsequent reconstruction scheme based on convex optimization. We prove that the reconstruction error due to quantization decays polynomially in the number of measurements. Our results apply to arbitrary signals, including compressible ones, and account for measurement noise. Additionally, they hold for sub-Gaussian (including Gaussian and Bernoulli) random compressed sensing measurements, as well as for both high bit-depth and coarse quantizers, and they extend to 1-bit quantization. In the noise-free case, when the signal is strictly sparse we prove that by optimizing the order of the quantization scheme one can obtain root-exponential decay in the reconstruction error due to quantization.     

Translational absolute continuity and Fourier frames on a sum of singular measures

A finite Borel measure μ in ${\mathbb{R}}^d$ is called a frame-spectral measure if it admits an exponential frame (or Fourier frame) for $L^2(\mu )$. It has been conjectured that a frame-spectral measure must be translationally absolutely continuous, which is a criterion describing the local uniformity of a measure on its support. In this paper, we show that if any measures ν and λ without atoms whose supports form a packing pair, then $\nu ?\lambda +{\delta }_t?\nu$ is translationally singular and it does not admit any Fourier frame. In particular, we show that the sum of one-fourth and one-sixteenth Cantor measure ${\mu }_4+{\mu }_{16}$ does not admit any Fourier frame. We also interpolate the mixed-type frame-spectral measures studied by Lev and the measure we studied. In doing so, we demonstrate a discontinuity behavior: For any anticlockwise rotation mapping $R_{\theta }$ with $\theta \ne \pm \pi /2$, the two-dimensional measure ${\rho }_{\theta }(?):=({\mu }_4\times {\delta }_0)(?)+({\delta }_0\times {\mu }_{16})(R_{\theta }^{?1}?)$, supported on the union of x-axis and $y=(\cot ?\theta )x$, always admit a Fourier frame. Furthermore, we can find ${\{e^{2\pi i〈\lambda ,x〉}\}}_{\lambda \in {\mathrm{\Lambda }}_{\theta }}$ such that it forms a Fourier frame for ${\rho }_{\theta }$ with frame bounds independent of θ. Nonetheless, ${\rho }_{\pm \pi /2}$ does not admit any Fourier frame.     

Norm inequalities for certain classes of functions and their Fourier transforms

We apply tools of interpolation theory and a commutative property of the Hilbert transform to prove necessary and sufficient conditions related to trigonometric series. These results extend and improve related theorems proven by several authors, summarized by Boas. In addition, we explore inequalities and operators, both connected to Hardy's inequalities, on certain classes of functions, including quasimonotone functions.     

On the recovery of a 2-D function from the modulus of its Fourier transform

A uniqueness theorem is proven for the problem of the recovery of a complex valued compactly supported 2-D function from the modulus of its Fourier transform. An application to the phase problem in optics is discussed.     

Some remarks on short‐time Fourier integral operators and classes of rapidly decaying functions

In this article, we extend the space of rapidly decaying functions to a space of rapidly decaying Boehmians. We provide convolution products, convolution theorems and generate their associated spaces of Boehmian. Then, we define the short‐time Fourier integral operator on the Boehmian spaces. Moreover, we show that the short‐time Fourier integral operator of the Boehmian is a sequentially continuous mapping that preserves certain desired properties. An inversion formula and some injections have also been obtained.     

Extensions and isomorphisms for the generalized Fourier algebras of a locally compact group

It is shown that for amenable groups, all finite-dimensional extensions of A_p(G) algebras split strongly. Furthermore, each extension of A_p(G) which splits algebraically also splits strongly. We also show that if G is an almost connected locally compact group, or a subgroup of GL_n(V) (V being a finite-dimensional vector space), and if for a fixed p∈(1,∞), all finite-dimensional singular extensions of A_p(G) split strongly, then G is amenable. Continuous order isomorphisms for the pointwise order of A_p(G) algebras, are characterized as weighted composition maps. Similarly, order isomorphisms for the pointwise order of B_p(G) algebras, are characterized as ∗-algebra isomorphisms followed by multiplication by an invertible positive multiplier. In addition, it is shown that for amenable groups, an order isomorphism for the pointwise order between A_p(G) algebras that preserve cozero sets is necessarily continuous, and hence induces an algebra isomorphism.     

A class of multiwavelets and projected frames from two-direction wavelets

This article aims at studying two-direction refinable functions and two-direction wavelets in the setting ?^s, s > 1. We give a sufficient condition for a two-direction refinable function belonging to L²(?^s). Then, two theorems are given for constructing biorthogonal (orthogonal) two-direction refinable functions in L²(?^s) and their biorthogonal (orthogonal) two-direction wavelets, respectively. From the constructed biorthogonal (orthogonal) two-direction wavelets, symmetric biorthogonal (orthogonal) multiwaveles in L²(?^s can be obtained easily. Applying the projection method to biorthogonal (orthogonal) two-direction wavelets in L²(?^s, we can get dual (tight) two-direction wavelet frames in L²(?^m, where. m ≥ s From the projected dual (tight) two-direction wavelet frames in L²(?^m, symmetric dual (tight) frames in L²(?^m can be obtained easily. In the end, an example is given to illustrate theoretical results.     

Fourier transform and distributional representation of the generalized gamma function with some applications

We present a Fourier transform representation of the generalized gamma functions, which leads to a distributional representation for them as a series of Dirac-delta functions. Applications of these representations are shown in evaluation of the integrals of products of the generalized gamma function with other functions. The results for Euler’s gamma function are deduced as special cases. The relation of the generalized gamma function with the Macdonald function is exploited to deduce new identities for it.     

Composition of Fourier Integral Operators with Fold and Blowdown Singularities

ABSTRACT

The purpose of this work is to present results about the composition of Fourier integral operators with certain singularities, for which the composition is not again a Fourier integral operator. The singularities considered here are folds and blowdowns. We prove that for such operators, the Schwartz kernel of F*F belongs to a class of distributions associated to two cleanly intersection Lagrangians. Such Fourier integral operators appear in integral geometry, inverse acoustic scattering theory and Synthetic Aperture Radar imaging, where the composition calculus can be used as a tool for finding approximate inversion formulas and for recovering images.     

Higher order Lipschitz classes of functions and absolutely convergent Fourier series

We introduce the higher order Lipschitz classes Λ _r (α) and λ _r (α) of periodic functions by means of the rth order difference operator, where r = 1, 2, ..., and 0 < α ≦ r. We study the smoothness property of a function f with absolutely convergent Fourier series and give best possible sufficient conditions in terms of its Fourier coefficients in order that f belongs to one of the above classes. This research was supported by the Hungarian National Foundation for Scientific Research under Grant T 046 192.