jz2.m

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

function D=jz2(A,E,q)
% 得到矩阵D
% A为迭代得到的矩阵
% E为各阶样本自相关函数的向量形式
% q为滑动平均阶数
sum=0;
for s=1:q+1
    for t=1:q-s+2
        sum=sum+A(t)*A(t+s-1);
    end
    D(s)=sum-E(s);
    sum=0;
end
D=D';