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.