Computational techniques based on the Lanczos representation
HTML articles powered by AMS MathViewer
- by J. N. Lyness PDF
- Math. Comp. 28 (1974), 81-123 Request permission
Abstract:
In his book Discourse on Fourier Series, Lanczos deals in some detail with representations of $f(x)$ of the type $f(x) = {h_{p - 1}}(x) + {g_p}(x)$ where ${h_{p - 1}}(x)$ is a polynomial of degree $p - 1$ and ${g_p}(x)$ has the property that its full range Fourier coefficients converge at the rate ${r^{ - p}}$. In Part I, some properties of ${h_p}(x)$ and of the series $\{ {h_p}(x)\} _1^\infty$ are described. These properties are used here to provide criteria for the convergence or divergence of the Euler-Maclaurin series, in the case when $f(x)$ is an analytic function. The similarities and differences between this series and the Lidstone and other two-point series are briefly mentioned. In Part II, the Lanczos representation is employed to derive an approximate representation $F(x)$ for an analytic function $f(x)$ on the interval [0, 1] is derived. This has the form \[ F(x) = \sum \limits _{q = 1}^{p - 1} {{\lambda _{q - 1}}{B_q}(x)/q! + 2\sum \limits _{r = 0}^{m/2} {({\mu _r}\cos 2\pi rx + {\nu _r}\sin 2\pi rx)} } \] and requires for its determination the values of the derivatives ${f^{(q - 1)}}(1) - {f^{(q - 1)}}(0)\;(q = 1,2, \cdots p - 1)$ and the regularly spaced function values $f(j/m)\;(j = 0,1, \cdots ,m)$. It involves replacing $g_p(x)$ by a discrete Fourier expansion based on trapezoidal rule approximations to its Fourier coefficients. This representation is a powerful one. The drawback is that it requires derivatives. Most of Part II is devoted to the effect of using only approximate derivatives. It is shown that when these are successively less accurate with increasing order (the sort of behaviour encountered using finite difference formula), then the representation is still powerful and reliable. In a computational context the only penalty for using inaccurate derivatives is that a larger value of m may—or may not—be required to attain a specific accuracy.References
- Milton Abramowitz and Irene A. Stegun, Handbook of mathematical functions with formulas, graphs, and mathematical tables, National Bureau of Standards Applied Mathematics Series, No. 55, U. S. Government Printing Office, Washington, D.C., 1964. For sale by the Superintendent of Documents. MR 0167642
- Ralph Philip Boas Jr., Entire functions, Academic Press, Inc., New York, 1954. MR 0068627
- William B. Jones and G. Hardy, Accelerating convergence of trigonometric approximations, Math. Comp. 24 (1970), 547–560. MR 277086, DOI 10.1090/S0025-5718-1970-0277086-X
- Cornelius Lanczos, Discourse on Fourier series, Hafner Publishing Co., New York, 1966. MR 0199629 G. J. Lidstone, "Notes on the extension of Aitken’s theorem (for polynomial interpolation) to the Everett types," Proc. Edinburgh Math. Soc. (2), v. 2, 1929, pp. 16-19.
- D. V. Widder, Completely convex functions and Lidstone series, Trans. Amer. Math. Soc. 51 (1942), 387–398. MR 6356, DOI 10.1090/S0002-9947-1942-0006356-4 P. C. Chakravarti, Integrals and Sums, Univ. of London, The Athlone Press, 1970. W. M. Gentleman & G. Sande, Fast Fourier Transforms—For Fun and Profit, Proc. AF1PS 1966 FJCC, v. 29, Spartan Books, New York, pp. 563-578. J. N. Lyness & G. Sande, "ENTCAF and ENTCRE: Evaluation of normalised Taylor coefficients of an analytic function, algorithm 413," Comm. ACM, v. 14, 1971, pp. 669-675.
Additional Information
- © Copyright 1974 American Mathematical Society
- Journal: Math. Comp. 28 (1974), 81-123
- MSC: Primary 65D15; Secondary 42A08
- DOI: https://doi.org/10.1090/S0025-5718-1974-0334458-6
- MathSciNet review: 0334458