Мне дали сюжет из текстового файла. Цель состоит в том, чтобы отфильтровать сигнал. Я строю график, чтобы дать ck (частоты сигнала), я получил результат, которого я не ожидал.
def dft(a): # Esta função serve para calcular os ck dos correspondentes coeficientes de Fourier
N = len(a)
c = np.zeros(N//2 +1,complex)
n = np.arange(0,N)
for k in range(N//2 + 1):
c[k] = sum(a*np.exp(-1j*2*np.pik*n/N))
return c
c=dft(a)
plt.plot(abs(c))
plt.show()
print(max(c))
Изображение графика:
