我试图绘制我的wav文件数据的频谱。我得到一个复数列表,通过用较小的数字进行测试,我知道这些值是正确的。但是,我不知道怎么策划。我只知道我应该得到一些像小条形图的东西,它们应该是对称的。但我的图表不是。下面是我的代码和当前的情节。你知道吗
from scipy.io.wavfile import read
import numpy as np
import matplotlib.pyplot as plt
import math
df = read('matches-4.wav')
data = np.array(df[1], dtype=float)[:, 1][2000:2512] # using right channel
def discrete_fourier_transformation():
c = []
for k in data:
summation = 0
for j in data:
summation += j * math.e ** (-1j * (2 * math.pi / len(data)) * j * k)
c.append(abs(1 / len(data) * summation))
return c
values = discrete_fourier_transformation()
plt.plot(values)
plt.show()
由于离散傅里叶变换的实现中存在问题,因此无法从图形中获得预期的对称性。你知道吗
更具体地说,在离散傅里叶变换实现中,复指数应该包含整数因子的乘法,我想这就是您尝试使用
j
和k
所做的。但是,您配置循环的方式j
和k
值是实际的数据值。要解决这个问题,您应该更新您的循环,使其在range(0,len(data)
范围内运行,如下所示:相关问题 更多 >
编程相关推荐