35 research outputs found
Complete Algebraic Reconstruction of Piecewise-Smooth Functions from Fourier Data
In this paper we provide a reconstruction algorithm for piecewise-smooth
functions with a-priori known smoothness and number of discontinuities, from
their Fourier coefficients, posessing the maximal possible asymptotic rate of
convergence -- including the positions of the discontinuities and the pointwise
values of the function. This algorithm is a modification of our earlier method,
which is in turn based on the algebraic method of K.Eckhoff proposed in the
1990s. The key ingredient of the new algorithm is to use a different set of
Eckhoff's equations for reconstructing the location of each discontinuity.
Instead of consecutive Fourier samples, we propose to use a "decimated" set
which is evenly spread throughout the spectrum
Decimated generalized Prony systems
We continue studying robustness of solving algebraic systems of Prony type
(also known as the exponential fitting systems), which appear prominently in
many areas of mathematics, in particular modern "sub-Nyquist" sampling
theories. We show that by considering these systems at arithmetic progressions
(or "decimating" them), one can achieve better performance in the presence of
noise. We also show that the corresponding lower bounds are closely related to
well-known estimates, obtained for similar problems but in different contexts
Algebraic Fourier reconstruction of piecewise smooth functions
Accurate reconstruction of piecewise-smooth functions from a finite number of
Fourier coefficients is an important problem in various applications. The
inherent inaccuracy, in particular the Gibbs phenomenon, is being intensively
investigated during the last decades. Several nonlinear reconstruction methods
have been proposed, and it is by now well-established that the "classical"
convergence order can be completely restored up to the discontinuities. Still,
the maximal accuracy of determining the positions of these discontinuities
remains an open question. In this paper we prove that the locations of the
jumps (and subsequently the pointwise values of the function) can be
reconstructed with at least "half the classical accuracy". In particular, we
develop a constructive approximation procedure which, given the first
Fourier coefficients of a piecewise- function, recovers the locations
of the jumps with accuracy , and the values of the function
between the jumps with accuracy (similar estimates are
obtained for the associated jump magnitudes). A key ingredient of the algorithm
is to start with the case of a single discontinuity, where a modified version
of one of the existing algebraic methods (due to K.Eckhoff) may be applied. It
turns out that the additional orders of smoothness produce a highly correlated
error terms in the Fourier coefficients, which eventually cancel out in the
corresponding algebraic equations. To handle more than one jump, we propose to
apply a localization procedure via a convolution in the Fourier domain
Local and global geometry of Prony systems and Fourier reconstruction of piecewise-smooth functions
Many reconstruction problems in signal processing require solution of a
certain kind of nonlinear systems of algebraic equations, which we call Prony
systems. We study these systems from a general perspective, addressing
questions of global solvability and stable inversion. Of special interest are
the so-called "near-singular" situations, such as a collision of two closely
spaced nodes.
We also discuss the problem of reconstructing piecewise-smooth functions from
their Fourier coefficients, which is easily reduced by a well-known method of
K.Eckhoff to solving a particular Prony system. As we show in the paper, it
turns out that a modification of this highly nonlinear method can reconstruct
the jump locations and magnitudes of such functions, as well as the pointwise
values between the jumps, with the maximal possible accuracy.Comment: arXiv admin note: text overlap with arXiv:1211.068