N = 4 * fs
X = scipy.fft(x[:N])
Xdb = 20 * scipy.log10(scipy.absolute(X))
f = scipy.linspace(0, fs, N, endpoint=False)
pylab.plot(f, Xdb)
pylab.xlim(0, 5000)
Y = X * H
y = scipy.real(scipy.ifft(Y))
Python: частотный анализ звуковых файлов