MK趋势分析

MK趋势分析

Mann-Kendall(MK)趋势分析在气象、水文学中已得到广泛应用，该方法不要求变量遵从正太分布，故可以检验水文、气象序列。

计算原理

对于一个时间序列Xt = （x1,x2,…,xn），MK趋势检验的统计量S计算公式如下：

其中，sgn为检验函数：
然后计算方差Var以及Z统计量：

利用双边检验对数据进行检验，当-Z1-a<=Zs<=Z1-a时，即表明序列具有显著变化趋势，反之趋势不显著。当Z1-a给定的显著水平a=0.05时，临界值为1.96。

工作环境

Python3.6、numpy

准备工作

首先导入numpy包，并输入数据。

计算S值及Var

注意索引值和我们平时排列计数时的区别。

s_end = []

for i in range(0,test<