Abstract
A nonlinear sequence transformation is presented which is able to accelerate the convergence of Fourier series. It is tailored to be exact for a certain model sequence. As in the case of the Levin transformation and other transformations of Levin-type, in this model sequence the partial sum of the series is written as the sum of the limit (or antilimit) and a certain remainder, i.e., it is of Levin-type. The remainder is assumed to be the product of a remainder estimate and the sum of the first terms oftwo Poincaré-type expansions which are premultiplied by two different phase factors. This occurrence of two phase factors is the essential difference to the Levin transformation. The model sequence for the new transformation may also be regarded as a special case of a model sequence based on several remainder estimates leading to the generalized Richardson extrapolation process introduced by Sidi. An algorithm for the recursive computation of the new transformation is presented. This algorithm can be implemented using only two one-dimensional arrays. It is proved that the sequence transformation is exact for Fourier series of geometric type which have coefficients proportional to the powers of a numberq∈ℂ, |q|<1. It is shown that under certain conditions the algorithm indeed accelerates convergence, and the order of the convergence is estimated. Finally, numerical test data are presented which show that in many cases the new sequence transformation is more powerful than Wynn's epsilon algorithm if the remainder estimates are properly chosen. However, it should be noted that in the vicinity of singularities of the Fourier series the new sequence transformation shows a larger tendency to numerical instability than the epsilon algorithm.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
G. Baszenski and F.-J. Delvos, Error estimates for sine series expansions, Math. Nachr. 139 (1988) 155–166.
C.M. Bender and S.A. Orszag,Advanced Mathematical Methods for Scientists and Engineers (McGraw-Hill, Singapore, 1987).
K.S. Crump, Numerical inversion of Laplace transforms using a Fourier series approximation, J. ACM 23 (1976) 89–96.
W.F. Ford and A. Sidi, An algorithm for a generalization of the Richardson extrapolation process, SIAM J. Numer. Anal. 24 (1987) 1212–1232.
W.B. Jones and G. Hardy, Accelerating convergence of trigonometric expansions, Math. Comp. 24 (1970) 547–560.
J.E. Kiefer and G.H. Weiss, A comparison of two methods for accelerating the convergence of Fourier series, Comput. Math. Appl. 7 (1981) 527–535.
D. Levin, Development of non-linear transformations for improving convergence of sequences, Int. J. Comput. Math. B3 (1973) 371–388.
D. Levin and A. Sidi, Two new classes of nonlinear transformations for accelerating the convergence of infinite integrals and series, Appl. Math. Comput. 9 (1981) 175–215.
I.M. Longman, The summation of Fourier, Chebyshev, and Legendre series, Appl. Math. Comput. 23 (1987) 61–70.
J.K. Shaw, L.W. Johnson and R.D. Riess, Accelerating convergence of eigenfunction expansions, Math. Comput. 30 (1976) 469–477.
A. Sidi, Some properties of a generalization of the Richardson extrapolation process, J. Inst. Math. Appl. 24 (1979) 327–346.
D.A. Smith and W.F. Ford, Numerical comparisons of nonlinear convergence accelerators, Math. Comput. 38 (1982) 481–499.
M. Tasche, Zur Konvergenzbeschleunigung von Fourier-Reihen, Math. Nachr. 90 (1979) 123–134.
E.J. Weniger, Nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series, Comput. Phys. Rep. 10 (1989) 189–373.
P. Wynn, On a device for computing thee m(S n) transformation, Math. Tables Aids Comput. 10 (1956) 91–96 [Kiefer and Weiss [6] cite also a technical report of Wynn on transformations to accelerate the convergence of Fourier series which could not be obtained by the present author as yet].