On the integration schemes of retrieving impulse response functions from transfer functions (Q624679)

Summary: The numerical inverse Laplace transformation (NILM) makes use of numerical integration. Generally, a high-order scheme of numerical integration renders high accuracy. However, surprisingly, this is not true for the NILM to the transfer function. Numerical examples show that the performance of higher-order schemes is no better than that of the trapezoidal scheme. In particular, the solutions from high-order scheme deviate from the exact one markedly over the rear portion of the period of interest. The underlying essence is examined. The deviation can be reduced by decreasing the frequency-sampling interval.