Ich habe mir viele Beispiele für scipy.fft und numpy.fft angesehen. Speziell dieses BeispielScipy / Numpy-FFT-Frequenzanalyse ist sehr ähnlich zu dem, was ich tun möchte. Daher habe ich die gleiche Subplot-Positionierung verwendet und alles sieht sehr ähnlich aus.

Ich möchte Daten aus einer Datei importieren, die nur eine Spalte enthält, um meinen ersten Test so einfach wie möglich zu gestalten.

Mein Code schreibt so:

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()

Das Problem ist, dass ich immer meine Spitze bei genau Null bekomme, was überhaupt nicht der Fall sein sollte. Es sollte wirklich bei ungefähr 200 Hz erscheinen.

Mit kleinerer Reichweite:

Immer noch größter Gipfel bei Null.