Построение БПФ в оцифрованных данных Python

Question

Я измеряю напряжение и пробую это. Из измеренных данных я хочу сделать БПФ, чтобы получить больше информации о гармониках напряжения. Тем не менее, результаты сбивают меня с толку... Является ли частота выборки, которую я должен здесь выбрать, той же частотой дискретизации, с которой я измеряю свое напряжение (и такой же для точек выборки)? Я приложил образец набора данных. Это должна быть синусоида с частотой 1 кГц. ссылка на данные: https://www.dropbox.com/s/bdv6uiy3lkeoomk/data.mat?dl=0

Код выглядит следующим образом:

%matplotlib inline
import numpy as np
import matplotlib.pyplot as fourier_plot
import scipy.fftpack as fft

# Number of samplepoints
N = 16384
# sample spacing
fs = 125e6/8
T = 1.0 /fs
x = np.linspace(0.0, N*T, N)
y = data
yf = scipy.fftpack.fft(y)
xf = np.linspace(0.0, 1.0/(2.0*T), N/2)


fourier_plot.plot(xf, 2.0/N * np.abs(yf[:N//2]))
fourier_plot.show()