用python绘制功率谱

Question

我有一个包含301个值的数组，这些值是从一个包含
301帧。这意味着1帧中有1个值。电影剪辑时间是30分钟
fps，实际上是10秒长
现在我想得到这个“信号”的功率谱
轴）。我试过：

 X = fft(S_[:,2]);
 pl.plot(abs(X))
 pl.show()

I also tried:

 X = fft(S_[:,2]);
 pl.plot(abs(X)**2)
 pl.show()

虽然我不认为这是真正的光谱。

有人能帮忙吗？

Answer 1

Numpy有一个方便的功能np.fft.fftfreq格式计算频率
与FFT组件相关：

from __future__ import division
import numpy as np
import matplotlib.pyplot as plt

data = np.random.rand(301) - 0.5
ps = np.abs(np.fft.fft(data))**2

time_step = 1 / 30
freqs = np.fft.fftfreq(data.size, time_step)
idx = np.argsort(freqs)

plt.plot(freqs[idx], ps[idx])

请注意，您看到的最大频率不是30Hz，而是30Hz

In [7]: max(freqs)
Out[7]: 14.950166112956811

在功率谱中你永远看不到采样频率。如果你有一个
如果样本数为偶数，那么您就可以到达[Nyquist]
频率](http://en.wikipedia.org/wiki/Nyquist\u频率)，15赫兹
（尽管numpy会将其计算为-15）。