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