3.3.3　程序设计

3.3.3.1　计算流程

按照逐步回归周期分析的基本思路及计算步骤，设计计算流程如图3.6 所示。

3.3.3.2　程序代码

逐步回归周期分析包括因子生成子程序、周期分析子程序、逐步回归子程序以及检验子程序。逐步回归分析子程序见2.3.6，F 检验及t 检验程序见2.2.5。

（1）因子生成子程序。程序输入原始序列Y（0 ton－1）及序列长度n、预报步数p，输出因子数m、预报因子及预报对象数组X_Y （1 tom＋1，0ton＋p－1），其中X_Y （m＋1，0 ton－1）为预报对象序列。该子程序需要调用周期分析子程序CycleGenerating。

（2）周期分析子程序。计算各预报因子序列不同周期长度的周期波振幅（即各组组平均值）。程序输入因子数m，样本容量n，原始系列YY （0 ton），输出周期波振幅V0（m＋1，m）。使用时在主程序中将V0 定义为动态数组。其程序代码如下：(https://www.daowen.com)

接下来就可以调用逐步回归分析子程序进行计算，调用检验程序进行F检验和t检验。