jz1.m

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

function T=jz1(A,q)
% 生成矩阵T
% A为得到的迭代向量
% q为滑动平均模型的阶数
p=0;
for s=1:q+1                     %构造矩阵B
    for t=1:q+1
        if t+p>q+1
            B(s,t)=0;
        else
            B(s,t)=A(t+p);
        end
    end
    p=p+1;
end
n=0;
for s=1:q+1                       %构造矩阵C
    for t=1:q+1
        if t-n<1
            C(s,t)=0;
        else
            C(s,t)=A(t-n);
        end
    end
    n=n+1;
end
T=B+C;