matlab计算经验半方差(变异函数)

前言

前段时间地统计的课设代码的确有点折磨人，但是还好算法过程比较简单，没有涉及到特别复杂的步骤，matlab在矩阵运算方面有大量的库函数以及其在索引方面的优势，使得我们真正能够实现向量化(矩阵化编程)，希望大家看到这篇文章过后能够真正取思考，以前的哪些代码可以将for循环等若干繁琐步骤去掉，减少代码行数，掌握向量化编程的精髓。

1 变异函数

这里引用一下杨老师ppt上的一个简单定义，一维条件下，当空间点$x$在一维$x$轴上变化时，区域化变量$Z(x)$在点$x$和$x+h$处的值$Z(x)$与$Z(x+h)$差的方差一半定义为区域化变量Z(x)在x轴方向上的变异函数。
通常区域化变量$Z(x)$所描述的现象是二维和三维的，则变异函数可定义为在任一方向$\alpha$，相距$|h|$的两个区域化变量值$Z(x)$与$Z(x+h)$的增量的方差。
$$\begin{gathered} 2\gamma (x,h)=Var[Z(x)-Z(x+h)] \\ =E[(Z(x)-Z(x+h){)}^2]-(E[(Z(x)-Z(x+h))]{)}^2 \end{gathered}$$
其中$\gamma (x,h)$称为半变异函数，但有时为方便也称变异函数
变异函数计算公式
$$\gamma (h)=\frac{1}{N(h)}\sum _{i=1}^n[Z(x_i)-Z(x_i+h){]}^2$$
其中$h$为距离(滞后距)
N ( h ) N(h)