3.2.3　程序设计

3.2.3.1　计算流程

根据周期波分析的计算方法，设计其计算流程如图3.3 所示，程序设置最多能识别5 个周期波。

3.2.3.2　程序代码

周期分析叠加子程序只需要输入原始序列Y0（n）、序列长度N、信度r 值，子程序就可输出识别的周期波数k2及周期波系列值YY1（0ton＋4，k2），只要将周期波按式（3.11）线性叠加即可获得模拟系列及预测值，叠加后实测序列以后的数据为预报值（共作5步预测）。预测程序非常简单，这里没有列出。该子程序还需要调用F 检验子程序ShiSuanF，其代码见2.2.5.2。周期分析叠加子程序VisualBasic6.0程序代码如下：(https://www.daowen.com)