现代中长期水文预报方法程序设计

图9.5　模糊综合分析预测模型程序流程图

9.6.2.1　程序流程图

根据模糊综合分析预测的基本原理及计算步骤设计程序流程如图9.5 所示。

9.6.2.2　程序代码

模糊综合分析预测程序包括模拟子程序、预测子程序、数据标准化、初始权重计算、互相关系数计算、综合分析指标uj 计算等子程序。其VisualBasic程序代码如下：

（1）数据模拟子程序。执行此程序前需要确定预报因子的权重及参数p。当然也可以将模拟误差最小作为目标，调用初始权重计算子程序计算预报因子的初始权重，对权重和参数p 进行试算，以便求得最优权重及参数p。程序输入预报因子X（M，N）、序列长度M、因子数N、预报对象Y（M）、因子权重W（M）以及模糊综合参数p，输出线性回归系数b0，b1 及模拟过程FY（M）。

程序需要调用数据标准化子程序StandardData、模糊综合因子计算子程序UUj 及一元线性回归子程序MonoRegression。

（2）模型预测子程序。该子程序应用模拟子程序Simulate率定的参数，进行预测计算。程序输入模型率定因子资料X（M，N）、预报对象Y（M）、因子数N、序列长度M、预报因子FYZ（N）、模糊综合参数p 以及一元回归系数b0，b1，输出预报结果FResult。

（3）数据标准化子程序。程序输入预报因子X（M，N）、序列长度M、因子数N，预报对象Y（M），输出预报因子标准化数XR（M，N）。该子程序需要调用计算序列最大、最小值的子程序、互相关系数计算子程序。

（4）计算序列最大、最小值的子程序。程序输入计算序列X（N）、序列长度N，输出序列最大值xMax、最小值xMin。

（5）互相关系数计算子程序。程序输入序列X（N）、Y（N）、序列长度N、滞时k，输出自相关系数RRk。

（6）初始权重计算子程序。程序输入建模因子X（M，nc）、预报对象Y（M）、序列长度M、因子数nc，输出预报因子的初始权重W（nc）。

（7）多因子模糊综合分析合成指标计算子程序。程序输入预报因子的权重W（nc）、标准化预报因子数据XR（nr，nc）、因子数nc、序列长度nr 以及参数p，输出综合分析合成指标Uj（nr）。(www.xing528.com)

（8）一元线性回归子程序。程序输入原始序列X（N）、Y（N）、序列长度N，输出回归系数b0、b1。