python - 使用Python对音频文件进行FFT-采样率

Question

我一直在尝试寻找一些在WAV文件中最主要的频率。我发现我可以使用python wave库以及scipy.io.wavfile来读取音频文件。

如果采样率为16000，那是否意味着我的时间步为1/16000？这意味着我的代码应如下所示查找频率:

wavData = wavfile.read(someName); # 2D array of sampling rate and values
fourier = fftpack.fft(wavData[0][1]); # FFT on the values
n = wavData[0][1].size;
timestep = 1/wavData[0]; # 1/16000
freq = np.fft.fftfreq(n, d=timestep); # getting array of sample frequencies

我正在使用 fftfreq 和 fft 的示例跟踪scipy.org的文档。

Answer 1

If the sampling rate is 16000, does that mean my timestep is 1/16000?

对，那是正确的;每秒16,000个样本意味着每个样本占1 / 16,000秒。但是，请注意，在Python 2.x中，您将需要显式使用 float来获取正确的 timestep:

timestep = 1.0 / wavData[0]

使用2.x的默认整数除法:

1 / 16000 == 0

鉴于

1.0 / 16000 == 6.25e-05