delta = 0.1
F = 1
X_delta = abs(freq - F) < delta
A_f_delta = abs(f[X_delta]) ** 2
maximum_ix = argrelextrema(A_f, np.greater)
E_F_delta = np.sum(A_f[maximum_ix]) / 2
E_F_delta
2453.8235776144866
Рассчитайте энергию для каждой полосы частот вокруг интересующей частоты F в Python.