Необходимо получить частоты присутствующих в сигнале гармоник. Изначально дан файл с самим сигналом и частота дискретизации 500 МГц. Сигнал выглядит следующим образом
Для того, чтобы получить частоты написал следующий код:
data = pd.read_csv('signal_1.txt', header=None)
from scipy.fft import fft, fftfreq
fs = 500 # Переведем частоту в МГц
N = data.size
"""Прямое БПФ"""
X = fft(data)
freq = fftfreq(N, 1./fs) # Частоты полученных сигналов
f = freq[:N//2]
A = abs(X[:N//2]) # Амплитуда сигнала
# Частотный спектр
plt.subplot(4, 2, 2)
plt.plot(f, A)
plt.title('Спектр')
plt.xlabel('Частота, МГц')
plt.ylabel('Амплитуда')
plt.grid()
plt.plot(f, A)
По итогу выводится следующая картина
Однако, как мне сказали, мы должны получить всего лишь 4 частоты. Я только начинаю работать с сигналами, но просмотрев уже достаточно много примеров все равно не могу понять, где я ошибся. Был бы рад, если бы мне указали на мою ошибку