AES E-Library

AES E-Library

A Contour Integral Method for Time-Domain Response Calculations

Document Thumbnail

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.

JAES Volume 66 Issue 5 pp. 360-368; May 2018
Publication Date:

Click to purchase paper as a non-member or login as an AES member. If your company or school subscribes to the E-Library then switch to the institutional version. If you are not an AES member and would like to subscribe to the E-Library then Join the AES!

This paper costs $33 for non-members and is free for AES members and E-Library subscribers.

Learn more about the AES E-Library

E-Library Location:


Start a discussion about this paper!

AES - Audio Engineering Society