skip to main content
article
Free Access

An implementation of a Fourier series method for the numerical inversion of the Laplace transform

Authors Info & Claims
Published:01 September 1999Publication History
Skip Abstract Section

Abstract

Our method is based on the numerical evaluation of the integral which occurs in the Riemann Inversion formula. The trapezoidal rule approximation to this integral reduces to a Fourier series. We analyze the corresponding discretization error and demonstrate how this expression can be used in the development of an automatic routine, one in which the user needs to specify only the required accuracy.

References

  1. 1 BLANCH, G. 1964. Numerical evalution of continued fractions. SIAM Rev. 6, 4, 383-421.Google ScholarGoogle Scholar
  2. 2 CRUMP, K. S. 1976. Numerical inversion of Laplace transforms using a Fourier series approximation. J. ACM 23, 1 (Jan.), 89-96. Google ScholarGoogle Scholar
  3. 3 D'ALESSIO, A., D'AMORE, L., AND LACCETTI, G. 1994. An effective discretization error estimate of Fourier series methods for the numerical inversion of the Laplace transform. Ricerche di Matematica 43, 2, 293-307.Google ScholarGoogle Scholar
  4. 4 DAVIES,B.AND MARTIN, B. 1979. Numerical inversion of the Laplace transform: A survey and comparison of methods. J. Comput. Phys. 33, 1, 1-32.Google ScholarGoogle Scholar
  5. 5 DE HOOG,F.R.,KNIGHT,J.K.,AND STOKES, A. N. 1982. An improved method for numerical inversion of Laplace Transforms. SIAM J. Sci. Stat. Comput. 3, 3, 357-366.Google ScholarGoogle Scholar
  6. 6 DOETSCH, G. 1950. Handbuch der Laplace Transformation. Vol. 1. Birkh~user-Verlag, Basel, Switzerland.Google ScholarGoogle Scholar
  7. 7 DOETSCH, G. 1955. Handbuch der Laplace Transformation. Vol. 2. Birkh~user-Verlag, Basel, Switzerland.Google ScholarGoogle Scholar
  8. 8 DOETSCH, G. 1956. Handbuch der Laplace Transformation. Vol. 3. Birkh~user-Verlag, Basel, Switzerland.Google ScholarGoogle Scholar
  9. 9 DUBNER,H.AND ABATE, J. 1968. Numerical inversion of Laplace transforms by relating them to the finite Fourier cosine transform. J. ACM 15, 1 (Jan.), 115-123. Google ScholarGoogle Scholar
  10. 10 DUFFY, D. G. 1993. On the numerical inversion of Laplace transforms: Comparison of three new methods on characteristic problems from applications. ACM Trans. Math. Softw. 19,3 (Sept. 1993), 333-359. Google ScholarGoogle Scholar
  11. 11 DURBIN, F. 1974. Numerical inversion of Laplace Transforms: An efficient improvement to Dubner and Abate's method. Comput. J. 17, 4, 371-376.Google ScholarGoogle Scholar
  12. 12 GARBOW,B.S.,GIUNTA, G., LYNESS,J.N.,AND MURLI, A. 1988. Software for an implementation of Weeks' method for the inverse Laplace transform. ACM Trans. Math. Softw. 14,2 (June), 163-170. Google ScholarGoogle Scholar
  13. 13 HONIG,G.AND HIRDES, U. 1984. A method for numerical inversion of Laplace transforms. J. Comput. Appl. Math. 10, 1 (Feb.), 113-132.Google ScholarGoogle Scholar
  14. 14 LACCETTI, G. 1992. The incidental parameters of a numerical method for inverting a Laplace transform function. Ricerche di Matematica 41, 1, 163-184.Google ScholarGoogle Scholar
  15. 15 MURLI, A. 1970. Sull'impiego di un metodo numerico per il calcolo dell'antitrasformata di Laplace. Rend. Accad. Sci. Mat. Soc. Naz. Scienze Lettere ed Arti, Napoli 37, 89-96.Google ScholarGoogle Scholar
  16. 16 MURLI,A.AND PATRUNO, V. 1978. Un metodo per l'inversione numerica della Trasformata di Laplace. CALCOLO 15, 51-58.Google ScholarGoogle Scholar
  17. 17 PIESSENS,R.AND HUYSMANS, R. 1984. Algorithm 619: Automatic numerical inversion of the Laplace transform. ACM Trans. Math. Softw. 10, 3 (Sept.), 348-353. Google ScholarGoogle Scholar
  18. 18 RUTISHAUSER, H. 1957. der Quotienten:Differenzen:Algorithmus. Birkh~user-Verlag, Basel, Switzerland.Google ScholarGoogle Scholar

Index Terms

  1. An implementation of a Fourier series method for the numerical inversion of the Laplace transform

              Recommendations

              Reviews

              Friedemann W. Stallmann

              A software package for the inversion of the Laplace transform using the Riemann inversion formula is described. This technique assumes, of course, that the transform is known in the complex plane, information not always available in engineering applications. Step size and integration boundaries are adjusted automatically to satisfy the desired tolerances. The theory behind the error estimates is discussed in some detail, and numerous examples are tabulated for functions whose transforms are known explicitly. The performance of this method is compared with that of other software packages.

              Access critical reviews of Computing literature here

              Become a reviewer for Computing Reviews.

              Comments

              Login options

              Check if you have access through your login credentials or your institution to get full access on this article.

              Sign in

              Full Access

              PDF Format

              View or Download as a PDF file.

              PDF

              eReader

              View online with eReader.

              eReader