prony2.m

来自「可以用来识别包装系统的总体参数以及原响应函数的初始相位和振幅」· M 代码 · 共 15 行

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)