33,456 research outputs found
On the numerical stability of Fourier extensions
An effective means to approximate an analytic, nonperiodic function on a
bounded interval is by using a Fourier series on a larger domain. When
constructed appropriately, this so-called Fourier extension is known to
converge geometrically fast in the truncation parameter. Unfortunately,
computing a Fourier extension requires solving an ill-conditioned linear
system, and hence one might expect such rapid convergence to be destroyed when
carrying out computations in finite precision. The purpose of this paper is to
show that this is not the case. Specifically, we show that Fourier extensions
are actually numerically stable when implemented in finite arithmetic, and
achieve a convergence rate that is at least superalgebraic. Thus, in this
instance, ill-conditioning of the linear system does not prohibit a good
approximation.
In the second part of this paper we consider the issue of computing Fourier
extensions from equispaced data. A result of Platte, Trefethen & Kuijlaars
states that no method for this problem can be both numerically stable and
exponentially convergent. We explain how Fourier extensions relate to this
theoretical barrier, and demonstrate that they are particularly well suited for
this problem: namely, they obtain at least superalgebraic convergence in a
numerically stable manner
Fast Algorithms for the computation of Fourier Extensions of arbitrary length
Fourier series of smooth, non-periodic functions on are known to
exhibit the Gibbs phenomenon, and exhibit overall slow convergence. One way of
overcoming these problems is by using a Fourier series on a larger domain, say
with , a technique called Fourier extension or Fourier
continuation. When constructed as the discrete least squares minimizer in
equidistant points, the Fourier extension has been shown shown to converge
geometrically in the truncation parameter . A fast algorithm has been described to compute Fourier extensions for the case
where , compared to for solving the dense discrete
least squares problem. We present two algorithms for
the computation of these approximations for the case of general , made
possible by exploiting the connection between Fourier extensions and Prolate
Spheroidal Wave theory. The first algorithm is based on the explicit
computation of so-called periodic discrete prolate spheroidal sequences, while
the second algorithm is purely algebraic and only implicitly based on the
theory
On the resolution power of Fourier extensions for oscillatory functions
Functions that are smooth but non-periodic on a certain interval possess
Fourier series that lack uniform convergence and suffer from the Gibbs
phenomenon. However, they can be represented accurately by a Fourier series
that is periodic on a larger interval. This is commonly called a Fourier
extension. When constructed in a particular manner, Fourier extensions share
many of the same features of a standard Fourier series. In particular, one can
compute Fourier extensions which converge spectrally fast whenever the function
is smooth, and exponentially fast if the function is analytic, much the same as
the Fourier series of a smooth/analytic and periodic function.
With this in mind, the purpose of this paper is to describe, analyze and
explain the observation that Fourier extensions, much like classical Fourier
series, also have excellent resolution properties for representing oscillatory
functions. The resolution power, or required number of degrees of freedom per
wavelength, depends on a user-controlled parameter and, as we show, it varies
between 2 and \pi. The former value is optimal and is achieved by classical
Fourier series for periodic functions, for example. The latter value is the
resolution power of algebraic polynomial approximations. Thus, Fourier
extensions with an appropriate choice of parameter are eminently suitable for
problems with moderate to high degrees of oscillation.Comment: Revised versio
- …