import pyfar as pf
impulse = pf.signals.impulse(100, 10)
pf.plot.freq_phase(impulse, unwrap=True)