We have a function which describes a frequency response correction to
piezo devices we use. To flatten the FFT, it is similar to:
Cdis_t = .5
N = 8192
for n in range(8192):
B3 = n * 2560 / N
Fc(n) = 1 / ((B3/((1/(Cdis_t*2*pi))**2+B3**2)**0.5)*(-0.01*log(B3)
In practice it really only matters for low frequencies.