python - Python 中的傅里叶变换

Question

我是使用 Python 进行信号处理的新手。我想找出如何将加速度计的幅度值转换为频域。我的示例代码如下:

在[44]中:

  x = np.arange(30)
  plt.plot(x, np.sin(x))
  plt.xlabel('Number of Sample')
  plt.ylabel('Magnitude Value')
  plt.show()

在这里我想将数据绘制到域频率。所需的输出可能是这样的:

Answer 1

numpy 和 scipy 有一个傅立叶变换模块 ( http://docs.scipy.org/doc/numpy/reference/routines.fft.html )。

x = np.linspace(0,5,100)
y = np.sin(2*np.pi*x)

## fourier transform
f = np.fft.fft(y)
## sample frequencies
freq = np.fft.fftfreq(len(y), d=x[1]-x[0])
plt.plot(freq, abs(f)**2) ## will show a peak at a frequency of 1 as it should.