python计算变差函数

变差函数通俗理解为计算一组数据的滞后距(步长)为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)