import pyfar as pf
impulse = pf.signals.impulse(100, 10)
pf.plot.group_delay(impulse, unit='samples')