Re: leasqr

[Gert Van den Eynde]

Hi Christian, Paul,


Another possible algorithm is Prony's algorithm. It's designed to fit
sums of exponential functions (c(i)*exp(a(i)*x)) but it also works for
imaginary a(i), and thus for sine/cosine pairs. One drawback is the
requirement that the data points need to be equidistant. My professor at
university always claimed this algorithm was more stable than the least
squares approach.... Some time ago, I implemented this for Octave but
never got around to post it to the mailing-list. So now, here it is...
Usual disclaimers applied. Hope it works. Comments and improvements
welcome.



Have a nice day,

Gert



*-*-*-*
Gert Van den Eynde
Reactor Physics & MYRRHA Dept.
SCK-CEN
Belgium
*-*-*-*

prony.m
Description: Text document