Jump to content

Inverse Laplace transform

From Wikipedia, the free encyclopedia

This is an old revision of this page, as edited by 205.220.129.18 (talk) at 10:19, 23 January 2022 (Post's inversion formula: Removed unnecessary and confusing argument. A function evaluated at an argument is merely a number.). The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

In mathematics, the inverse Laplace transform of a function F(s) is the piecewise-continuous and exponentially-restricted real function f(t) which has the property:

where denotes the Laplace transform.

It can be proven that, if a function F(s) has the inverse Laplace transform f(t), then f(t) is uniquely determined (considering functions which differ from each other only on a point set having Lebesgue measure zero as the same). This result was first proven by Mathias Lerch in 1903 and is known as Lerch's theorem.[1][2]

The Laplace transform and the inverse Laplace transform together have a number of properties that make them useful for analysing linear dynamical systems.

Mellin's inverse formula

An integral formula for the inverse Laplace transform, called the Mellin's inverse formula, the Bromwich integral, or the FourierMellin integral, is given by the line integral:

where the integration is done along the vertical line Re(s) = γ in the complex plane such that γ is greater than the real part of all singularities of F(s) and F(s) is bounded on the line, for example if contour path is in the region of convergence. If all singularities are in the left half-plane, or F(s) is an entire function , then γ can be set to zero and the above inverse integral formula becomes identical to the inverse Fourier transform.

In practice, computing the complex integral can be done by using the Cauchy residue theorem.

Post's inversion formula

Post's inversion formula for Laplace transforms, named after Emil Post,[3] is a simple-looking but usually impractical formula for evaluating an inverse Laplace transform.

The statement of the formula is as follows: Let f(t) be a continuous function on the interval [0, ∞) of exponential order, i.e.

for some real number b. Then for all s > b, the Laplace transform for f(t) exists and is infinitely differentiable with respect to s. Furthermore, if F(s) is the Laplace transform of f(t), then the inverse Laplace transform of F(s) is given by

for t > 0, where F(k) is the k-th derivative of F with respect to s.

As can be seen from the formula, the need to evaluate derivatives of arbitrarily high orders renders this formula impractical for most purposes.

With the advent of powerful personal computers, the main efforts to use this formula have come from dealing with approximations or asymptotic analysis of the Inverse Laplace transform, using the Grunwald–Letnikov differintegral to evaluate the derivatives.

Post's inversion has attracted interest due to the improvement in computational science and the fact that it is not necessary to know where the poles of F(s) lie, which make it possible to calculate the asymptotic behaviour for big x using inverse Mellin transforms for several arithmetical functions related to the Riemann hypothesis.

Software tools

See also

References

  1. ^ Cohen, A. M. (2007). "Inversion Formulae and Practical Results". Numerical Methods for Laplace Transform Inversion. Numerical Methods and Algorithms. Vol. 5. p. 23. doi:10.1007/978-0-387-68855-8_2. ISBN 978-0-387-28261-9.
  2. ^ Lerch, M. (1903). "Sur un point de la théorie des fonctions génératrices d'Abel". Acta Mathematica. 27: 339. doi:10.1007/BF02421315.
  3. ^ Post, Emil L. (1930). "Generalized differentiation". Transactions of the American Mathematical Society. 32 (4): 723–723. doi:10.1090/S0002-9947-1930-1501560-X. ISSN 0002-9947.
  4. ^ Abate, J.; Valkó, P. P. (2004). "Multi-precision Laplace transform inversion". International Journal for Numerical Methods in Engineering. 60 (5): 979. doi:10.1002/nme.995.

Further reading

This article incorporates material from Mellin's inverse formula on PlanetMath, which is licensed under the Creative Commons Attribution/Share-Alike License.