Shock Response Spectrum
From Python Wiki
Shock Response Spectrum (SRS)
Mechanical shock pulses are often analyzed in terms of shock response spectra. The shock response spectrum assumes that the shock pulse is applied as a common base input to an array of independent single-degree-of-freedom systems. The shock response spectrum gives the peak response of each system with respect to the natural frequency. Damping is typically fixed at a constant value, such as 5%, which is equivalent to an amplification factor of Q=10.
The following tutorial gives further information: 
The SRS for an an arbitrary base input can be calculate via a digital recursive filtering relationship, which models the convolution integral.
The filtering relationship can be implemented in Python by importing the lfilter function from scipy.
from scipy.signal import lfilter
A script for doing this is given at: 
The script requires the utility module: 
Home Page: []