本文共 682 字，大约阅读时间需要 2 分钟。

k=0;  for i = 0:3    for j = 0:3 %0:L,L的值不确定         if i == 0 & j == 0            continue        elseif i == 0            ToEstMd = arima('MALags',1:j,'Constant',0); %指定模型的结构        elseif j == 0            ToEstMd = arima('ARLags',1:i,'Constant',0); %指定模型的结构        else            ToEstMd = arima('ARLags',1:i,'MALags',1:j,'Constant',0); %指定模型的结构        end        k = k + 1;        R(k) = i;        M(k) = j;        [EstMd,EstParamCov,LogL,info] = estimate(ToEstMd,da);        %模型拟合,估计模型参数        numParams = sum(any(EstParamCov));        %计算拟合参数的个数        [aic(k),bic(k)] = aicbic(LogL,numParams,n);    endendfprintf('R,M,AIC,BIC的对应值如下\n%f');%显示计算结果check  = [R',M',aic',bic'];