python - 了解 matplotlib magnitude_spectrum 输出

Question

我在理解 matplotlib 的 magnitude_spectrum 函数调用生成的输出时遇到问题。

我生成了一个频率为 50khz 的正弦波，

f_s = 488000.0 # Hz
t = np.arange(0.0, 1.0, 1/f_s)
s1 = 100*np.sin(2*np.pi*50000*t)

然后我在除以 FFT bins 的数量后绘制结果幅度谱

s1_magspec = plt.magnitude_spectrum(s1,Fs=f_s)
plt.plot(s1_magspec[0]/len(s1_magspec[0]))

结果是 50khz 的单个尖峰，但幅度为 50，而不是预期的 100。

谁能解释这是为什么？

这是一个 ipython 笔记本的链接，描述了显示上述代码和结果图:

http://nbviewer.ipython.org/gist/bkinman/22cc15d3ad3b9b2db09e

Answer 1

看起来它与使用的 FFT 窗口的默认设置有关。文档说默认是汉宁窗。如果您改为使用棚车窗:

s1_magspec = plt.magnitude_spectrum(s1, Fs=f_s, window=np.ones(s1.shape))

您将在 100 处获得峰值，就像直接使用 numpy FFT 一样。

顺便说一句，如果您编辑您的问题以放入创建 s1_magspec 的代码行而不是依赖笔记本查看器，那将会很好，指向它的链接无疑有一天会中断。