Its popularity is due primarily to the fact that It is also a relatively straight forward approach, going back to ideas from Lanczos Weeks method is one of the most well known algorithms for the numerical inversion of a scalar Additional information about the first three numerical inversion algorithms can be found in the following links: Each has found a domain of applicationĬorresponding to the ability of the algorithm to invert certain classes of From experimentation and review, four main algorithms for numerical Laplace transform inversion have proved to be of use, the Post-Widderįormula ABATE2004, VALKO2001, VALKO2004, Fourier Series Expansion DEHOOG1982, Talbot's method TALBOT1979, and the Weeks One interpretation of the various numerical approaches to the inversion integral are that they areĭifferent regularization techniques. Īlgorithmic and finite precision errors can lead to exponential divergence of numerical solutions. One considers the need to multiply by a potentially increasing large exponential. Numerical analysis point of view, the inherent sensitivity of the inversion procedure is clear when Property can be seen as a special case of the difficulty inherent in solving an inverse problemĭescribed by a Fredholm integral equation of the first kind DAVIESB2002. From a functional analysis perspective, this can The numerical inversion of the Laplace transform Of a particular solution is not a fundamental problem. Furthermore, since inĪpplications to physical problems the time domain functions are generally continuous, the selection However, due to Lerch's theorem, one can in practice assume that the analytic inverse is well-defined. Two time domain functions whichĭiffer at a single point in time for example will have the same transform. The inversion of the Laplace transform is not a unique operation. The position along the real axis is known as the abscissa of convergence and is chosen by convention to be larger than For isolated singularities, theīromwich contour is the standard approach. In general, for a time domain scalar function withĬan be evaluated by complex integration along a contour. The latter is a well-known application of the theory of complex To understand numerical inversion of the Laplace transform, it is first necessary to understand the analytic Also provided are MATLAB scripts for computing the matrix exponential with this method.Īnalytic Inversion of the Laplace Transform For additional introductory material see recent overview talk Numerical Laplace Transform Inversion and Selected Applications. Material here comes from this published PAPERĪnd the dissertation by P. Transform to the computation of the matrix exponential. The purpose of this website is to introduce a recently studiedĪpplication of the Weeks method for the numerical inverse Laplace The numerical inverse Laplace transform is howeverĪn inherently sensitive procedure and thus requires special OneĪpproach is to compute the matrix function from a numerical inverse It can be used in many fields of mathematics, such as functional analysis, and certain areas of theoretical and applied probability.The accurate and reliable computation of the matrix exponentialįunction is a long standing problem of numerical analysis. For real-valued functions, it is the Laplace transform of a Stieltjes measure, however it is often defined for functions with values in a Banach space. The Laplace-Stieltjes transform, named for Pierre-Simon Laplace and Thomas Joannes Stieltjes, is an integral transform similar to the Laplace transform. In the past few decades, many mathematicians studied the growth and value distribution of the analytic (entire) function defined by Dirichlet series and obtained lots of interesting results (see ).Īs we know, Dirichlet series is regarded as a special example of the Laplace-Stieltjes transform. When a n, λ n, n satisfy some conditions, the series ( 1) is convergent in the whole plane or the half-plane, that is, f( s) is an analytic function or entire function in the whole plane or the half-plane. S = σ + i t ( σ, t are real variables), a n are nonzero complex numbers.
0 Comments
Leave a Reply. |