变差函数通俗理解为计算一组数据的滞后距(步长)为h时的方差的一半。公式如下
N(h)为是滞后距(步长)为h的数据对的数目, Ps(i) - Ps(i+h)为一个数据对。
以h=2为例,计算如下
# h为滞后距
# array为计算的数据集
def variogram(h):
'''变差函数'''
array = [1,3,6,7,8]
vario_1 = 0
for i in range(2,len(array),2):
g = (array[i-2]-array[i]) ** 2
vario_1 = vario_1+g
vario_2 = vario_1/(2*(len(array)//h+1))
print(vario_2)
variogram(2)