请高手帮忙解释下下面程序的作用（关于ARMA）

guangdangsky

function c=mainfcoeff(ar,ma,n);
%MAINFCOEFF 求与ARMA(p,q)相对应的MA的N个系数.
%   C=MAINFCOEFF(AR,MA,N) returns N coefficients of a MA(inf) process
%   corresponding to an ARMA(p,q) process with autoregressive and
%   moving average coefficients given by vectors AR=[1,a_1,...,a_p] and
%   MA=[1,b_1,...,b_q], respectively.
% 分别确定AR、MA的阶
p = length(ar);
q = length(ma);
% 计算系数
c = [ma zeros(1,n)];
if p>1
    ar = ar(2:p);
    d = ar;
    for i = 1:n
        c(i+1:i+q) = c(i+1:i+q)-d(1)*ma;
        d = [d(2:p-1) 0]-d(1)*ar;
    end;
end;
c = c(1:n+1);
主要是棕色部分的作用，对应的公式是什么？麻烦各位了

[ 本帖最后由 zhangnan3509 于 2010-6-2 16:26 编辑 ]