3,053 research outputs found
A matrix pencil approach to the existence of compactly supported reconstruction functions in average sampling
The aim of this work is to solve a question raised for average sampling in shift-invariant spaces by using the well-known matrix pencil theory. In many common situations in sampling theory, the available data are samples of some convolution operator acting on the function itself: this leads to the problem of average sampling, also known as generalized sampling. In this paper we deal with the existence of a sampling formula involving these samples and having reconstruction functions with compact support. Thus, low computational complexity is involved and truncation errors are avoided. In practice, it is accomplished by means of a FIR filter bank. An answer is given in the light of the generalized sampling theory by using the oversampling technique: more samples than strictly necessary are used. The original problem reduces to finding a polynomial left inverse of a polynomial matrix intimately related to the sampling problem which, for a suitable choice of the sampling period, becomes a matrix pencil. This matrix pencil approach allows us to obtain a practical method for computing the compactly supported reconstruction functions for the important case where the oversampling rate is minimum. Moreover, the optimality of the obtained solution is established
Oversampling in shift-invariant spaces with a rational sampling period
8 pages, no figures.It is well known that, under appropriate hypotheses, a sampling formula allows us to recover any function in a principal shift-invariant space from its samples taken with sampling period one. Whenever the generator of the shift-invariant space satisfies the Strang-Fix conditions of order r, this formula also provides an approximation scheme of order r valid for smooth functions. In this paper we obtain sampling formulas sharing the same features by using a rational sampling period less than one. With the use of this oversampling technique, there is not one but an infinite number of sampling formulas. Whenever the generator has compact support, among these formulas it is possible to find one whose associated reconstruction functions have also compact support.This work has been supported by the Grant MTM2009-08345 from the D.G.I. of the Spanish Ministerio de Ciencia y TecnologĂa
Weighted frames of exponentials and stable recovery of multidimensional functions from nonuniform Fourier samples
In this paper, we consider the problem of recovering a compactly supported
multivariate function from a collection of pointwise samples of its Fourier
transform taken nonuniformly. We do this by using the concept of weighted
Fourier frames. A seminal result of Beurling shows that sample points give rise
to a classical Fourier frame provided they are relatively separated and of
sufficient density. However, this result does not allow for arbitrary
clustering of sample points, as is often the case in practice. Whilst keeping
the density condition sharp and dimension independent, our first result removes
the separation condition and shows that density alone suffices. However, this
result does not lead to estimates for the frame bounds. A known result of
Groechenig provides explicit estimates, but only subject to a density condition
that deteriorates linearly with dimension. In our second result we improve
these bounds by reducing the dimension dependence. In particular, we provide
explicit frame bounds which are dimensionless for functions having compact
support contained in a sphere. Next, we demonstrate how our two main results
give new insight into a reconstruction algorithm---based on the existing
generalized sampling framework---that allows for stable and quasi-optimal
reconstruction in any particular basis from a finite collection of samples.
Finally, we construct sufficiently dense sampling schemes that are often used
in practice---jittered, radial and spiral sampling schemes---and provide
several examples illustrating the effectiveness of our approach when tested on
these schemes
Designing Gabor windows using convex optimization
Redundant Gabor frames admit an infinite number of dual frames, yet only the
canonical dual Gabor system, constructed from the minimal l2-norm dual window,
is widely used. This window function however, might lack desirable properties,
e.g. good time-frequency concentration, small support or smoothness. We employ
convex optimization methods to design dual windows satisfying the Wexler-Raz
equations and optimizing various constraints. Numerical experiments suggest
that alternate dual windows with considerably improved features can be found
- …