Accurate reconstructions of functions of finite regularity from truncated Fourier series expansions

Kowledge of a truncated Fourier series expansion for a 2π-periodic function of finite regularity, which is assumed to be piecewise smooth in each period, is used to accurately reconstruct the corresponding function. An algebraic equation of degree M is constructed for the M singularity locations in each period for the function in question. The M coefficients in this algebraic equation are obtained by solving an algebraic system of M equations determined by the coefficients in the known truncated expansion. If discontinuities in the derivatives of the function are considered, in addition to discontinuities in the function itself, that algebraic system will be nonlinear with respect to the M unknown coefficients. The degree of the algebraic system will depend on the desired order of accuracy for the reconstruction, i.e., a higher degree will normally lead to a more accurate determination of the singularity locations. By solving an additional linear algebraic system for the jumps of the function and its derivatives up to the arbitrarily specified order at the calculated singularity locations, we are able to reconstruct the 2π-periodic function of finite regularity as the sum of a piecewise polynomial function and a function which is continuously differentiab1e up to the specified order

[1]  T. MacRobert.  Higher Transcendental Functions , 1955, Nature.

[2]  F. B. Hildebrand,et al.  Introduction To Numerical Analysis , 1957 .

[3]  J. C. Eilbeck.  Table errata: Higher transcendental functions. Vol. I, II (McGraw-Hill, New York, 1953) by A. Erdélyi, W. Magnus, F. Oberhettinger and F. G. Tricomi , 1971 .

[4]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[5]  Peter D. Lax,et al.  Accuracy and Resolution in the Computation of Solutions of Linear and Nonlinear Equations , 1978 .

[6]  V. Klema.  LINPACK user's guide , 1980 .

[7]  M. Powell,et al.  Approximation theory and methods , 1984 .

[8]  Steven A. Orszag,et al.  Spectral Calculations of One-Dimensional Inviscid Compressible Flows , 1981 .

[9]  Gene H. Golub,et al.  Matrix computations , 1983 .

[10]  J. N. Lyness.  The Calculation of Trigonometric Fourier Coefficients , 1984 .

[11]  Begnaud Francis Hildebrand,et al.  Introduction to numerical analysis: 2nd edition , 1987 .

[12]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[13]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[14]  Walter A. Strauss,et al.  Computation of exponentials in transient data , 1989 .

[15]  William H. Press,et al.  Numerical recipes , 1990 .

[16]  Knut S. Eckhoff.  Accurate and efficient reconstruction of discontinuous functions from truncated series expansions , 1993 .

[17]  David Gottlieb,et al.  Issues in the Application of High Order Schemes , 1993 .