harm_pow {Momocs}

Calculates harmonic power given a list from e/t/rfourier

Description

Given a list with an, bn (and eventually cn and dn), returns the harmonic power.

Usage

harm_pow(xf)

Arguments

xf

A list with an, bn (and cn, dn) components, typically from a e/r/tfourier passed on coo_

Value

Returns a vector of harmonic power

Examples

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')