arlsq.m

应用经典prony法和AR参数建模进行模态参数识别

function a=arlsq(x,order)
n=length(x);
for i=1:n-order
    q(i,1:order)=x(order+i-1:-1:i);
end
y=x(order+1:n)';
a=(q'*q)^(-1)*(q'*y);
a=[1,-a'];