import pyfar as pf
signal = pf.signals.impulse(441)
signal = pf.dsp.filter.bell(signal, 1e3, 12, 1, "III")
signal = pf.dsp.filter.bell(signal, 10e3, -60, 100, "III")
smoothed, _ = pf.dsp.smooth_fractional_octave(signal, 1)
ax = pf.plot.freq(signal, label="input")
pf.plot.freq(smoothed, label="smoothed")
ax.legend(loc=3)