tquan.m

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

function T=Tquan(l,a,b,p,q,max)
% T权的计算
% l为预测期
% a是自回归参数
% b是滑动平均参数
% p是自回归阶数
% q是滑动平均参数
% max为要算权的个数
sum=0;
I=Iquan(a,b,p,q,max);      % I权的计算
for j=1:max
    T(j,1)=I(j);
end
if l>1
for j=1:max+1-l
    for s=2:l
        for k=1:s-1
        sum=sum+I(s-k)*T(j,k);
        end
        T(j,s)=I(j+s-1)+sum;
        sum=0;
    end
end
end