时空序贯指示模拟和序贯高斯模拟的步骤

时空序贯指示模拟和序贯高斯模拟的步骤

一、时空序贯指示模拟
1、每个样点进行四次指示变换，小于那个阈值取1，否则取0；
2、计算四次，得到四个理论变异函数；
3、定义时空格网矩阵和随机路径；
4、随机取下一格网，对四个指示变量stok；
5、概率值内插或外推，得累积概率分布，随机取一值，作为已知点代入下一点模拟；
6、所有时空网格是否完成模拟，是，得到一次模拟结果；否，重复4-6步；
7、重复3-6步1000次。

二、序贯高斯模拟
1、确定单变量的累积分布函数；
2、将原始数据Z进行标准正态变换；
3、在某个节点进行简单克里格插值，并计算其相应的克里格方差；
4、根据克里格方差，产生一随机残差函数，该函数满足均值为0，方差等于克里格方差的正态分布；
5、将产生的残差R(x)加入到克里格估值中，Z(x)=Z*(x)+R(x)；
6、将该点产生的模拟值加入现有数据，以其作为以后模拟的条件数据；
7、采用随机顺序，逐一访问所有需要模拟的节点，重复上述计算，直到所有点都完成模拟；
8、将模拟结果进行逆高斯变换，还原为原始数据；
9、利用不同的随机“种子”数，就产生不同的模拟实现。

 
 
想了解更多，请关注微信公众号“GAO戏精光芒万丈”：