clear all
prony1
jie=ans(1:2);
for i=1:1001
    for j=1:2
        r(i,j)=jie(j)^(i-1);
    end
end
A=pinv(r)*x';
ampli=abs(A(1))*2
phase=atan(imag(A(1))/real(A(1)))+pi/2
pp=-log(abs(jie(1)))/(2*pi*0.001);
qq=atan(imag(jie(1))/real(jie(1)))/(2*pi*0.001);
frequency=sqrt(pp^2+qq^2)
dampratio=sqrt(1-(qq/frequency)^2)