iquan.m

对平稳时间序列中的参数估计方法（包括矩估计

function I=Iquan(a,b,p,q,max)
% I权的计算
% a是自回归参数
% b是滑动平均参数
% p是自回归阶数
% q是滑动平均参数
% max是要算的最后一个I权
sum=0;
for j=1:max
    for h=1:j
        if h>q
            break;
        end
        if j==h
            sum=sum-b(j);
        else
            sum=sum+b(h)*I(j-h);
        end
    end
    if j>p
        I(j)=sum;
    else
        I(j)=sum+a(j);
    end
    sum=0;
end