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

2.4.3.1　程序计算流程

根据多元门限回归的建模步骤，设计计算流程如图2.10所示。其中多元回归计算过程没有列出，其计算流程见2.2.5。

2.4.3.2　程序代码

多元门限回归模型包括门限变量选择、门限回归、多元回归、预报计算及F 检验子程序等。其中多元回归及F检验子程序见2.2节。

图2.10　多元门限回归模型计算流程图

（1）门限变量选择子程序。程序输入预报因子X（M，N）、预报对象Y（N）及因子数M和序列长度N，输出各因子的互相关系数RRK（M）及选定的门限元序号MenXianYNo。该子程序需要调用互相关系数计算子程序。其程序代码如下：

（2）互相关系数计算子程序。程序输入序列X （N）、Y （N）、序列长度N 及滞时k，输出自相关系数RRk。代码如下：

（3）多元门限回归计算子程序。程序输入门限元序号MenXianYNo，预报因子X （M，N），预报对象Y（N），统计检验的信度r0 及序列长度N 和因子数M；输出门限值TX，门限元的最大统计量计算值Fs 及查算值F0，二分割的显著性检验结果FFJG（布尔型），两段的回归系数B1 （0toM）、B2 （0toM）。该程序需要调用排序子程序PaiXu和PaiXu1、统计量F计算子程序、F检验子程序、多元回归子程序，PaiXu子程序见2.3.6，F 检验及多元回归子程序见2.2.5。

（4）根据门限元进行排序子程序。该程序根据X（iKey，N）将X（M，N）从小到大排序，使用时将所有因子和预报对象都放在二维数组X 中，将门限元从小到大排序，其他数据相应重新排序。程序输入X（M，N）及M、N、排序根据的因子序号iKey，输出排序后的序列XX（M，N）。其程序代码如下：

（5）统计量F计算子程序。程序输入预报对象值YP （N）（对应预报因子从小到大排序后的序列）、序列长度N，输出统计量F（N－1）。其程序代码如下：

该子程序需要调用计算组内方差和子程序s2。子程序s2 输入第一组数据XX1 （N1）和第二组数据XX2 （N2）、两组数据长度N1 和N2，输出组内方差和SS。其程序代码如下：

（6）预报子程序。程序输入因子数M，两段的回归系数B1（0toM）、B2（0toM），门限值TX，预报因子值ZX（1 toM），选定的门限元序号MenXianYNo；输出预报值ZY。