3. The method of claim 1 , wherein said at least one set of second harmonic model component values is removed in a plurality of iterations so that during each one of said plurality of iterations the following is performed until a remaining harmonic component cost function is below a threshold: analyzing new harmonic model of previous harmonic model residual to produce new set of harmonic model component values, removing said new set of harmonic component values from said previous harmonic model residual to produce a new harmonic model residual for further iterations.
3. The method of claim 1 , wherein said at least one set of second harmonic model component values is removed in a plurality of iterations so that during each one of said plurality of iterations the following is performed until a remaining harmonic component cost function is below a threshold: analyzing new harmonic model of previous harmonic model residual to produce new set of harmonic model component values, removing said new set of harmonic component values from said previous harmonic model residual to produce a new harmonic model residual for further iterations. 5. The method of claim 3 , wherein said new harmonic modeling uses at least one estimated energy envelope signal.