Kenny S. Crump
Abstract
A method is presented for numerically inverting a Laplace transform that requires, in addition to the transform function itself, only sine, cosine, and exponential functions. The method is conceptually much like the method of Dubner and Abate, which approximates the inverse function by means of a Fourier cosine series. The method presented here, however, differs from theirs in two important respects. First of all, the Fourier series contains additional terms involving the sine function selected such that the error in the approximation is less than that of Dubner and Abate and such that the Fourier series approximates the inverse function on an interval of twice the length of the corresponding interval in Dubner and Abate's method. Second, there is incorporated into the method in this paper a transformation of the approximating series into one that converges very rapidly. In test problems using the method it has routinely been possible to evaluate inverse transforms with considerable accuracy over a wide range of values of the independent variable using a relatively few determinations of the Laplace transform itself.
- 1 DUI~NER, H , AND ABATE, J Numerical reversion of Laplace transforms by relating them to the fimte Fourier cosine transform J ACM I5, 1 (Jan 1968), 115-123. Google Scholar
Digital Library
- 2 DRANE, R E A numerical reversion techmque for Laplace transforms. Master's Th, Louisiana Tech U, Ruston, La, 1973Google Scholar
- 3 SIMON, R M, STROOT, M T, AND WEISS, C_~ H Numermal inversion of Laplace transforms with apphcatlon to percentage labeled mitoses experiments Computers and Bwmed Res. o (1972), 596-607.Google Scholar
- 4 MACDONALD, J R. Accelerated convergence, divergence, itelation, extrapolation, and curve fitting J Appl Phys 10 (1964), 3034-3041Google ScholarCross Ref
Index Terms
- Numerical Inversion of Laplace Transforms Using a Fourier Series Approximation
Recommendations
On the numerical inversion of Laplace transforms: comparison of three new methods on characteristic problems from applications
Three frequently used methods for numerically inverting Laplace transforms are tested on complicated transforms taken from the literature. The first method is a straightforward application of the trapezoidal rule to Bromwich's integral. The second ...
Numerical Inversion of Laplace Transforms by Relating Them to the Finite Fourier Cosine Transform
In this paper the problem of readily determining the inverse Laplace transform numerically by a method which meets the efficiency requirements of automatic digital computation is discussed. Because the result inverse function is given as a Fourier ...
Comments