Miré en muchos ejemplos de scipy.fft y numpy.fft. Específicamente este ejemploAnálisis de Frecuencia Scipy / Numpy FFT Es muy similar a lo que quiero hacer. Por lo tanto, utilicé el mismo posicionamiento de subtrama y todo parece muy similar.

Quiero importar datos de un archivo, que contiene una sola columna para que mi primera prueba sea lo más fácil posible.

Mi código escribe así:

import numpy as np
import scipy as sy
import scipy.fftpack as syfp
import pylab as pyl

# Read in data from file here
array = np.loadtxt("data.csv")
length = len(array)
# Create time data for x axis based on array length
x = sy.linspace(0.00001, length*0.00001, num=length)

# Do FFT analysis of array
FFT = sy.fft(array)
# Getting the related frequencies
freqs = syfp.fftfreq(array.size, d=(x[1]-x[0]))

# Create subplot windows and show plot
pyl.subplot(211)
pyl.plot(x, array)
pyl.subplot(212)
pyl.plot(freqs, sy.log10(FFT), 'x')
pyl.show()

El problema es que siempre obtendré mi pico exactamente en cero, lo que no debería ser así. Realmente debería aparecer a unos 200 Hz.

Con menor rango:

Sigue siendo el pico más grande a cero.