zhgddgj.m

来自「对平稳时间序列中的参数估计方法（包括矩估计」· M 代码 · 共 22 行

p=input('输入自回归的阶数：')
for j=1:p
    a(j)=input('输入样本自相关函数：');
end
b(1,1)=a(1);
sum1=0;
sum2=0;
for k=1:p-1
    for j=1:k
        sum1=sum1+b(k,j)*a(j);
        sum2=sum2+b(k,j)*a(k+1-j);
        b(k+1,k+1)=(a(k+1)-sum2)/(1-sum1);
    end
    for j=1:k
        b(k+1,j)=b(k,j)-b(k+1,k+1)*b(k,k+1-j);
    end
    sum1=0;
    sum2=0;
end
for j=1:p
    disp(b(p,j))
end