4.2.5　程序设计

4.2.5.1　计算流程

根据自回归模型原理及参数估计方法，设计其参数估计流程如图4.1 所示。有了模型参数就可以根据式（4.19）进行平稳时间序列预测。

4.2.5.2　程序代码

平稳序列自回归模型程序包括参数估计子程序、预报子程序和统计检验子程序。其中χ2检验程序见3.1.4.2，不同的是统计量用式（4.17）计算。

（1）AR模型参数估计子程序。该子程序输入原始序列Y （0 ton－1）及样本容量n，输出模型阶数p、均值AVEX、方差SGM、回归系数BC（1 top）。使用时在主程序需要定义动态数组BC（）。程序代码如下：

图4.1　自回归模型参数识别流程图(https://www.daowen.com)

（2）预报子程序。程序输入原始序列Y （0 ton－1）、样本容量n、模型阶数p、均值AVEX 及回归系数BC （1 to p），输出拟合（预测）序列ZZ （p ton），其中ZZ （p to n－1）为拟合序列值，ZZ （n）为预测值。其代码如下：