python音频分析工具,Python的：声音文件的频率分析

I am generating some sound files that play tones at various frequencies with a certain number of harmonics.

Ultimately, these sounds will be played on a device with a small speaker.

I have the frequency response curve of the speaker and want to do the following in Python:

Plot the frequency spectrum of sound file. I need a take the FFT of the file and plot it with gnuplot

Apply a nonlinear transfer function based on the frequency response curve in the data sheet.

Plot the result after the function is applied.

Does anyone know :

What the simplest way to do this would be?

or of an Application (GNU/Linux based) that could do this for me?

解决方案

I know you didn't mention Pylab/Matplotlib, but it works. Here is an example (assumes single-channel signal):

x, fs, nbits = audiolab.wavread('schubert.wav')

audiolab.play(x, fs)

N = 4*fs # four seconds of audio

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) # view up to 5 kHz

Y = X*H

y = scipy.real(scipy.ifft(Y))