首页 理论教育 现代中长期水文预报方法及其应用

现代中长期水文预报方法及其应用

时间:2023-10-04 理论教育 版权反馈
【摘要】:4.3.3.1程序流程图根据均生函数的计算步骤,设计其流程图如图4.3 所示。互相关系数计算子程序。

现代中长期水文预报方法及其应用

4.3.3.1 程序流程图

根据均生函数的计算步骤,设计其流程图如图4.3 所示(多元回归的流程见2.2.5)。

4.3.3.2 程序代码

均生函数模型程序包含均生函数构造子程序、周期延拓子程序、影响因子选择子程序、多元回归及预报子程序等,其中多元回归程序见第2 章,其他程序介绍如下。

(1)均生函数构造子程序。程序输入原始时间序列X(N)、序列长度N,输出均生函数XJ(Lm,Lm),Lm =int(N/2)。代码如下:

(2)周期延拓子程序。程序输入均生函数XJ(LM,LM),原始序列长度N 与预报预见期k 之和NN,输出外延后的序列XE(LM,NN),其代码如下:

(3)影响因子选择子程序。程序输入延拓后的因子序列XE(Lm,N)、序列长度N、需要选择的影响因子数qYZS,输出所选因子的序号及其互相关系数R_XH(2,Lm),其中Lm =int(N/2)。该子程序需要调用相关系数计算及排序子程序,程序代码如下:

该子程序执行以后,就可以调用多元回归子程序计算多元回归参数B (0 toqYZS)。

(4)互相关系数计算子程序。程序输入序列X(N)、Y(N)、序列长度N 及滞时k,输出自相关系数RRk。代码如下:

(5)排序子程序。程序输入数组X(N)、Y(N)及其长度N,根据X 绝对值的大小将X、Y 从小到大排序。代码如下:

(6)预报(模拟)子程序。多元回归方程系数计算以后就可以调用预报子程序进行模拟或预报计算。程序输入因子数qYZS、延拓后的序列XE(Lm,N)、选择的因子序号R_XH(1,LM)及回归系数B(0 toqYZS),输出模拟序列YF(N),程序代码如下:

如果进行预报,子程序输入参数增加预见期k,子程序作如下改进:

图4.4 均生函数三角矩阵(www.xing528.com)

图4.5 均生函数延拓后形成的因子序列

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈