ef <- efourier(bot[1], 24)
rf <- efourier(bot[1], 24)
harm_pow(ef)
harm_pow(rf)
plot(cumsum(harm_pow(ef)[-1]), type='o',
main='Cumulated harmonic power without the first harmonic',
ylab='Cumulated harmonic power', xlab='Harmonic rank')
Run the code above in your browser using DataLab