通过fft找到wav文件的频率幅度和相位

时间:2016-04-08 08:11:28

标签: python numpy fft wav spectrogram

我有一个wav文件,我想得到频率,幅度和相位。我试图这样做,但我得到的是我无法获得翻译

import math
import numpy as np
from matplotlib.pyplot import *
import scipy.io.wavfile as wave
from numpy.fft import fft

rate,data = wave.read('test.wav')
n = data.size
duree = 1.0*n/rate
print rate
spectre = np.fft.fft(data[5:10])
#freq = np.fft.fftfreq(n, 1)
print spectre

我获得了例如

[[ -9.27766766e+08+0.j  -9.27557398e+08+0.j]
[ -1.86505703e+09+0.j   2.16973235e+09+0.j]
[ -2.33588876e+08+0.j   2.33467572e+08+0.j]
[  1.76254287e+09+0.j   1.76250750e+09+0.j]
[  9.96780365e+08+0.j  -2.30269509e+09+0.j]]

1 个答案:

答案 0 :(得分:2)

您必须学习FFT算法才能理解所有内容。如您所见,它管理复杂的数字,因此需要一些添加工作来解释输出。

作为捷径,输出的一半是redundant。要查看频谱,请按以下步骤进行:

import numpy as np
import matplotlib.pyplot as plt
import scipy.io.wavfile as wave

rate,data = wave.read('57.wav')
spectre = np.fft.fft(data)
freq = np.fft.fftfreq(data.size, 1/rate)
mask=freq>0   
plt.plot(freq[mask],np.abs(spectre[mask]))

对于

spectrum