arlsq.m

来自「应用经典prony法和AR参数建模进行模态参数识别」· M 代码 · 共 8 行

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'];