# Frequency and damping estimation with Prony method?

I have a sinusoidal time domain signal with few damped exponential components.
i.e.

fn=[0.4 0.6 1.0];E=[0.03 0.05 0.1];A=[1.0 2.0 3.0];
wn=2*pi.*fn;a=E.*wn;wd=wn.*sqrt(1-E.^2);
x=A(1)*exp(-a(1)*t).*sin(wd(1)*t)+A(2)*exp(-a(2)*t).*...
sin(wd(2)*t)+A(3)*exp(-a(3)*t).*sin(wd(3)*t);


How can I extract the damping ratios and frequency components of this signal by Prony method.