import pyfar as pf
impulse = pf.signals.impulse(100, 10)
plots = [[pf.plot.time, pf.plot.phase],
         [pf.plot.freq, pf.plot.group_delay]]
pf.plot.custom_subplots(impulse, plots)