半方差的模型的检验
模型的检验（cross-validation，又称作jackknifing)为了检验所选模型三个参数的合理性，必须作一定的检验.但是还没有一个有效的方法检验参数的置信区间；同时，由于我们不知道半方差模型的确切形式，所选定的模型只是半方差函数的近似式，故无法以确切的函数形式对模型参数进行统计检验.交叉验证法的检验方法，一种间接的结合普通克立格的方法，为检验所选模型的参数提供了一个途径.这个方法的优点是在检验过程中对所选定的模型参数不断进行修改，直至达到一定的精度要求.交叉验证法的基本思路是：依次假设每一个实测数据点未被测定，由所选定的半方差模型，根据n-1个其它测定点数据用普通克立格估算这个点的值.设测定点的实测值为，估算值为，通过分析误差,来检验模型的合理性.
半方差的方差函数的理论模型
土壤在空间上是连续变异的，所以土壤性质的半方差函数应该是连续函数.但是，样品半方差图却是由一批间断点组成.可以用直线或曲线将这些点连接起来，用于拟合的曲线方程就称为半方差函数的理论模型.在土壤研究中常用的模型有：
式中C1/a是直线的斜率.这是一维数据拟合的最简单模型：
（（h)=C0
+C1·h/a
0在极限情况下，C1/a可以为0，这时就有纯块金效应模型：
（（h)=C0,h>0
⑷
（（0)=0
h=0
（（h)=
C0
+C1[1.5h/a-0.5(h/a)3]
0a
⑸
（（0)=0
h=0
（（h)=C0+C1[1-exp-h/a
]
h>0
⑹
（（0)=0
h=0
（（h)=C0+C1[1-exp(-h2/a2)]
h>0
⑻
（（0)=0
h=0
选定了半方差函数的拟合模型后，通常是以最小二乘法计算方程的参数，并应用Ross等的最大似然程序（MLP），得到效果最好的半方差方程.
半方差的介绍
半方差函数(Semi-variogram)及其模型，半方差函数也称为半变异函数，它是地统计学中研究土壤变异性的关键函数。如果随机函数Z(x)具有二阶平稳性，则半方差函数((h)可以用Z(x)的方差S2和空间协方差C(h)来定义：((h)= S2-C(h)。((h)反映了Z(x)中的空间相关部分，它等于所有以给定间距h相隔的样点测值之差平方的数学期望：式中N(h)是以h为间距的所有观测点的成对数目。某个特定方向半方差函数图通常是由((h)对h作图而得。通常半方差函数值都随着样点间距的增加而增大，并在一定的间距(称为变程，arrange)升大到一个基本稳定的常数(称为基台，sill)。
半方差的原理
设x0为未观测的需要估值的点，x1,x2，…，xN 为其周围的观测点，观测值相应为y(x1 ),y(x2），…，y(xN).未测点的估值记为 (x0），它由相邻观测点的已知观测值加权取和求得：此处，i为待定加权系数.和以往各种内插法不同，Kriging内插法是根据无偏估计和方差最小两项要求来确定上式中的加权系数（i的，故称为最优内插法.1. 无偏估计 设估值点的真值为y(x0).由于土壤特性空间变异性的存在，以及，y(x0）均可视为随机变量.当为无偏估计时，⑽将式⑼代入⑽式，应有⑾2. 估值和真值y(x0）之差的方差最小.即⑿利用式（3-10），经推导方差为⒀式中，（（xi,xj）表示以xi和xj两点间的距离作为间距h时参数的半方差值，（（xi,x0）则是以xi和x0两点之间的距离作为间距h时参数的半方差值.观测点和估值点的位置是已知的，相互间的距离业已知，只要有所求参数的半方差（（h）图，便可求得各个（(xi,xj）和（（xi,x0）值.因此，确定式⑼中各加权系数的问题，就是在满足式⑾的约束条件下，求目标函数以式⒀表示的方差为最小值的优化问题.求解时可采用拉格朗日法，为此构造一函数，（为待定的拉格朗日算子.由此，可导出优化问题的解应满足：i=1,2,N ⒁由式⒁和式⑾组成n+1阶线性方程组，求解此线性方程组便可得到n个加权系数（i和拉格朗日算子（.该线性方程组可用矩阵形式表示：⒂式中，（ ij为（（xi,xj）的简写.求得各（i值和（值后，由式⑼便可得出x0点的最优估值y(x0).而且还可由式⒀求出相应该估值的方差之最小值（2min.将式⒁代入式⒀，最小方差值还可由下式方便地求出：⒃上述最优化问题求解还可用其他方法，在应用Kriging内插法时还有其他方面的问题，在此都不一一列举了。
