A Contour Integral Method for Time-Domain Response Calculations

Candy J, Futtrup C. A Contour Integral Method for Time-Domain Response Calculations. J Audio Eng Soc. 2018;66(5):360-368. doi:10.17743/jaes.2018.0017

@article{Candy2018_19566,
  author = {Candy, Jeff and Futtrup, Claus},
  title = {{A Contour Integral Method for Time-Domain Response Calculations}},
  journal = {Journal of the Audio Engineering Society},
  volume = {66},
  number = {5},
  pages = {360--368},
  year = {2018},
  month = may,
  publisher = {Audio Engineering Society},
  doi = {10.17743/jaes.2018.0017},
  url = {https://doi.org/10.17743/jaes.2018.0017}
}

TY  - JOUR
TI  - A Contour Integral Method for Time-Domain Response Calculations
AU  - Candy, Jeff
AU  - Futtrup, Claus
T2  - Journal of the Audio Engineering Society
J2  - J. Audio Eng. Soc.
VL  - 66
IS  - 5
SP  - 360
EP  - 368
PY  - 2018
DA  - 2018/05/06
DO  - 10.17743/jaes.2018.0017
UR  - https://doi.org/10.17743/jaes.2018.0017
PB  - Audio Engineering Society
LA  - en
AB  - A fundamental aspect of loudspeaker modeling is the ability to calculate the time-domain response of a driver in free air or in a specified enclosure given prior calculation of the frequency response. This report presents a numerical method to compute the time-domain response of a loudspeaker or other transducer using contour integration of its frequency response. The approach is based on Weideman’s scheme for inversion of the Laplace transform along a parabolic contour, and it is applicable to analytic functions that contain isolated singularities (poles) and branch points in the left half-plane. The new approach is motivated by the need to implement viscoelastic and semi-inductance effects into linear and nonlinear, time-dependent transducer calculations. Because the response functions that describe these phenomena contain fractional power and logarithmic singularities, solution methods based on rational function decomposition cannot be used. The new method is simple to implement, requires few function evaluations, and is remarkably accurate. For analyzing the step response in box models, the proposed contour method based on Weideman can be universally applied with only a few additional lines of new code.
ER  -