5.4.3　程序设计

5.4.3.1　程序计算流程图

根据灾变预测的数学原理和灾变GM（1，1）模型的计算步骤，设计程序计算流程如图5.5 所示。

5.4.3.2　程序代码

灰色灾变GM（1，1）模型包括灾变映射、GM（1，1）模型参数率定、GM（1，1）模型模拟与预报及灾变时间换算等子程序。其中有关GM（1，1）模型的程序见5.3.2.3，这里仅介绍灾变映射及灾变时间换算子程序。

（1）灾变映射子程序。程序输入上灾变或下灾变的标识变量zKey （1 表示上灾变，2表示为下灾变）、灾变阈值KS，原始资料X0（m）及序列长度m；输出映射集（日期或序号集）DX（dm）及其序列长度dm。其程序代码如下：

图5.5　灾变GM（1，1）模型程序计算流程图(https://www.daowen.com)

将该子程序的输出DX（dm）作为GM（1，1）模型的输入，即可计算出GM（1，1）模型的参数，进而进行预测。

（2）灾变时间换算子程序。程序输入预报步长FT、GM （1，1）模型预测值FX（FT）及原始资料开始年份bYear，输出灾变时间Fdata（FT）。代码如下：