import pyfar as pf
import numpy as np
import matplotlib.pyplot as plt
x = pf.signals.impulse(2**12)
y = pf.dsp.filter.reconstructing_fractional_octave_bands(x)[0]
f = pf.dsp.filter.fractional_octave_frequencies(
   frequency_range=(60, 16000))[1]
y_sum = pf.Signal(np.sum(y.time, 0), y.sampling_rate)
ax = pf.plot.time_freq(y_sum, color='k', unit='ms')
pf.plot.time(x, ax=ax[0], unit='ms')
ax[0].set_xlim(-20, 250)
ax[0].set_title("Original (blue) and reconstructed pulse (black)")
pf.plot.freq(y_sum, color='k', ax=ax[1])
pf.plot.freq(y, ax=ax[1])
ax[1].set_title(
    "Reconstructed (black) and filtered impulse (colored)")