Momocs (version 0.2-03)

harm.pow: 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

data(bot)
ef <- efourier(bot@coo[[1]], 24)
rf <- efourier(bot@coo[[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")