我从txt文件得到了一个图。目的是过滤信号。 我绘制了一个图形以给出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))
图的图像: