hr function

Estimate harmonic components