.地统计学方法.ppt

* 华中农业大学 资源与环境学院 泛克里格法 漂移的概念 在普通的克里格法中，要求区域化变量Z(x)是二阶平稳的或本证的，至少是准平稳或准本证假设条件，在有限的估计领域内Z(x)的数学期望是一个常数，即E[Z(x)]=m存在。然而在许多情况下，区域化变量在研究区域内是非平稳的，其数学期望不是一个常数，即E[Z(x)]=m(x)，m(x)在地统计学上被定义为非平稳区域化变量的漂移，其表达式为 m(x)=E[Z(x)] * 华中农业大学 资源与环境学院 泛克里格法 在漂移存在的条件下就不能用普通克里格方法进行空间局部估计，而要采用泛克里格法进行估计。漂移一般采用多项式表示： 其中fi(x)为一已知多项式函数，ai为未知系数。当漂移为线性时，一维和二维条件下漂移m(x)的形式为： * 华中农业大学 资源与环境学院 泛克里格法 当漂移为二次（非线性）时，一维和二维条件下漂移m(x)的形式为 由于漂移的存在，泛克里格法在估计某一点Z(x)的估计值Z#(x)时，首先要估计该点上漂移m(x)的估计值m#(x)，这就要求在某种假设条件下确定非平稳区域化变量的协方差函数和变异函数 * 华中农业大学 资源与环境学院 泛克里格的基本假设 设Z(x)是一个非平稳的区域化变量，可表示为： 一般情况下，上述假设很难满足，因此，考虑Z(x)的增量的情况。 设Z(x)的增量具有非平稳数学期望和非平稳方差函数，并可表示为： * 华中农业大学 资源与环境学院 泛克里格的基本假设 如果非平稳区域化变量Z(x)可以分解成两部分，一部分是在较大尺度下可以观察到现象变化m(x)，另一部分是在较小尺度下的变化R(x)，即 在给定的尺度下，m(x)可以表示为一个多项式，即单项式函数fi(x)的线性组合： 其中x为领域内任一点，ai为未知系数，如果包含x0点，则ai应有n+1个 * 华中农业大学 资源与环境学院 非平稳条件下的协方差和变异函数 当Z(x)=m(x)+R(x)时，Z(x)的协方差函数CZ(x,y)为 而Z(x)的变异函数为 * 华中农业大学 资源与环境学院 即Cz(x,y)=CR(x,y)，Z(x)的协方差等于R(x)的协方差 rz(x,y)=rR(x,y)，Z(x)的变异函数等于R(x)的变异函数 如y=x+h，则rz(h)=rR(h)。只要能求出R(x)的变异函数rR(h)，就可以求得Z(x)的变异函数rz(h)。但m(x)一般为未知多项式函数，无法用R(x)=Z(x)-m(x)来计算rR(h)。为解决这个问题，地统计学泛克里格法中先对R(x)的变异函数进行估计，也就是先估计R#(x)=Z(x)-m#(x)来，即对m#(x)进行估计，根据R#(x)的变异函数rR#(h)与理论的rR(h)进行比较，当rR#(h)=rR(h)时，就可以求Z(x)的变异函数rR(h)来，因此，泛克里格方法估计有两个部分，一个是m(x)的估计，另一个是Z(x)的估计 * 华中农业大学 资源与环境学院 漂移m(x)的泛克里格法估计 在研究区域内，设Z(x)是一个非平稳区域化随机变量，并满足以下条件： * 华中农业大学 资源与环境学院 漂移m(x)的泛克里格法估计 设在研究区内有n个已知样点xa(a=1,2,…,n），其观测值为Z(xa)=Za，目的是用这些已知样点估计区域内任一固定点x处的漂移值m(x)。设m(x)的估计量为m#(x)，它可以表示这几个已知样点数据的线性组合，即 * 华中农业大学 资源与环境学院 漂移m(x)的泛克里格法估计 为使m#(x)成为m(x)的无偏最优估计，则必须满足： 1 无偏性，即 只有当 时，无偏性存在，此时L=0,1,2,…,n * 华中农业大学 资源与环境学院 漂移m(x)的泛克里格法估计 2 最优性 m#(x)估计m(x)的方差可表示为 即，上式达到最小 * 华中农业大学 资源与环境学院 漂移m(x)泛克里格线性方程组 在满足无偏性条件下，用拉格朗日乘数法求方差的最小值得到漂移m(x)的克里格线性方程组 可矩阵表达为 * 华中农业大学 资源与环境学院 其中 * 华中农业大学 资源与环境学院 这是一个具有n+K+1个未知数的n+K+1个方程组成的方程组，称为漂移m(x)泛克里格线性方程组，泛克里格估计方差为： * 华中农业大学 资源与环境学院 Z(x)的泛克里格法估计 在研究区域内设Z(x)是一个非平稳区域化变量，并满足以下条件： * 华中农业大学 资源与环境学院 Z(x)的泛克里格法估计 设在研究区内有n个已知样点xa(a=1,2,…,n），其观测值为Z(xa)=Za，目的是用这些已知样点估计区域内任一固定点x处的漂移值Z(x)。设Z(x)的估计量为Znk#(x