import pyfar as pf
import numpy as np
impulses = pf.signals.impulse(
    2048, np.arange(0, 25), np.linspace(1, .5, 25))
pf.plot.group_delay_2d(impulses, unit="samples")