def f(Y, x, N):
    total = 0
for ctr in range(len(Y)):
    total += Y[ctr] * (np.cos(x * ctr * 2 * pi / N) + 1j * np.sin(x * ctr *
        2 * pi / N))
return real(total)
tempdata = np.loadtxt('sunspots.dat')
year = tempdata[:, (0)]
wolfer = tempdata[:, (1)]
Y = fft(wolfer)
N = len(Y)
print(N)
xs = list(range(N))
gplt.plot(xs, [f(Y, x, N) for x in xs])
gplt.show()
Воссоздание данных временных рядов с использованием результатов FFT без использования ifft