Differential approximation of the Gaussian by short cosine sums with exponential error decay
Abstract
In this paper, we propose a method to approximate the Gaussian function on R by a short cosine sum. We generalise and extend the differential approximation method proposed in [4, 40] to approximate e-t2/2σ in the weighted space L2( R, e-t2/2) where σ, \, >0. We prove that the optimal frequency parameters λ1, … , λN for this method in the approximation problem λ1,…, λN, γ1, …, γN\|e-·2/2σ - Σj=1N γj \, eλj ·\|L2( R, e-t2/2), are zeros of a scaled Hermite polynomial. This observation leads us to a numerically stable approximation method with low computational cost of O(N3) operations. We derive a direct algorithm to solve this approximation problem based on a matrix pencil method for a special structured matrix. The entries of this matrix are determined by hypergeometric functions. For the weighted L2-norm, we prove that the approximation error decays exponentially with respect to the length N of the sum. An exponentially decaying error in the (unweighted) L2-norm is achieved using a truncated cosine sum. Our new convergence result for approximation of Gaussian functions by exponential sums of length N shows that exponential error decay rates e-cN are not only achievable for complete monotone functions.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.