import pyfar as pf
import numpy as np
impulses = pf.signals.impulse(
    2048, np.arange(0, 25), np.linspace(1, .5, 25))
pf.plot.phase_2d(impulses, unwrap=True, freq_scale="linear")