IMPROVED FFT-BASED NUMERICAL INVERSION OF LAPLACE TRANSFORMS VIA FAST HARTLEY TRANSFORM ALGORITHM

被引:9
|
作者
HWANG, C [1 ]
LU, MJ [1 ]
SHIEH, LS [1 ]
机构
[1] UNIV HOUSTON,DEPT ELECT ENGN,HOUSTON,TX 77204
基金
美国国家航空航天局;
关键词
D O I
10.1016/0898-1221(91)90021-U
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The disadvantages of numerical inversion of the Laplace transform via the conventional fast Fourier transform (FFT) are identified and an improved method is presented to remedy them. The improved method is based on introducing a new integration step length DELTA-omega = pi/mT for trapezoidal-rule approximation of the Bromwich integral, in which a new parameter m, is introduced for controlling the accuracy of the numerical integration. Naturally, this method leads to multiple sets of complex FFT computations. A new inversion formula is derived such that N equally-spaced samples of the inverse Laplace transform function can be obtained by [m/2] + 1 sets of N-point complex FFT computations or by m sets of real fast Hartley transform (FHT) computations.
引用
收藏
页码:13 / 24
页数:12
相关论文
共 50 条