- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
克里格法
二、克里格法(Kriging)
克里格法(Kriging)是地统计学的主要内容之一,从统计意义上说,是从变量相关性和变异性出发,在有限区域内对区域化变量的取值进行无偏、最优估计的一种方法;从插值角度讲是对空间分布的数据求线性最优、无偏内插估计一种方法。克里格法的适用条件是区域化变量存在空间相关性。 克里格法,基本包括普通克里格方法(对点估计的点克里格法和对块估计的块段克里格法)、泛克里格法、协同克里格法、对数正态克里格法、指示克里格法、折取克里格法等等。随着克里格法与其它学科的渗透,形成了一些边缘学科,发展了一些新的克里金方法。如与分形的结合,发展了分形克里金法;与三角函数的结合,发展了三角克里金法;与模糊理论的结合,发展了模糊克里金法等等。 应用克里格法首先要明确三个重要的概念。一是区域化变量;二是协方差函数,三是变异函数 一、区域化变量 当一个变量呈空间分布时,就称之为区域化变量。这种变量反映了空间某种属性的分布特征。矿产、地质、海洋、土壤、气象、水文、生态、温度、浓度等领域都具有某种空间属性。区域化变量具有双重性,在观测前区域化变量Z(X)是一个随机场,观测后是一个确定的空间点函数值。 区域化变量具有两个重要的特征。一是区域化变量Z(X)是一个随机函数,它具有局部的、随机的、异常的特征;其次是区域化变量具有一般的或平均的结构性质,即变量在点X与偏离空间距离为h的点X+h处的随机量Z(X)与Z(X+h)具有某种程度的自相关,而且这种自相关性依赖于两点间的距离h与变量特征。在某种意义上说这就是区域化变量的结构性特征。 二、协方差函数 协方差又称半方差,是用来描述区域化随机变量之间的差异的参数。在概率理论中,随机向量X与Y的协方差被定义为:
区域化变量 在空间点x和x+h处的两个随机变量Z(x)和Z(x+h)的二阶混合中心矩定义为Z(x)的自协方差函数,即
区域化变量Z(x) 的自协方差函数也简称为协方差函数。一般来说,它是一个依赖于空间点x 和向量h 的函数。 设Z(x) 为区域化随机变量,并满足二阶平稳假设,即随机函数Z(x)的空间分布规律不因位移而改变, h为两样本点空间分隔距离或距离滞后,Z(xi)为Z(x)在空间位置xi处的实测值, Z(xi+h)是Z(x)在x处距离偏离h的实测值 ,根据协方差函数的定义公式,可得到协方差函数的计算公式为:
? 在上面的公式中, N(h)是分隔距离为h时的样本点对的总数, 和 分别为和的样本平均数,即
在公式中N为样本单元数。一般情况下 (特殊情况下可以认为近似相等)。若 (常数),协方差函数可改写为如下:
式中, m为样本平均数,可由一般算术平均数公式求得,即
三、变异函数 变异函数又称变差函数、变异矩,是地统计分析所特有的基本工具。在一维条件下变异函数定义为,当空间点x 在一维x 轴上变化时,区域化变量Z(x)在点x和x+h 处的值Z(x) 与Z(x+h) 差的方差的一半为区域化变量Z(x) 在 x轴方向上的变异函数,记为 ,即
在二阶平稳假设条件下,对任意的h有, 因此上式可以改写为:
从上式可知,变异函数依赖于两个自变量x 和h ,当变异函数 仅仅依赖于距离h 而与位置x 无关时,可改写成 ,即
设 Z(x)是系统某属性 Z在空间位置x处的值, Z(x)为一区域化随机变量,并满足二阶平稳假设, h为两样本点空间分隔距离,Z(xi) 和 Z(xi+h)分别是区域化变量在空间位置xi 和xi+h 处的实测值[i=1,2,...,N(h)] ,那么根据上式的定义,变异函数 的离散公式为:
变异函数揭示了在整个尺度上的空间变异格局,而且变异函数只有在最大间隔距离1/2处才有意义。 四、克里格估计量 假设x是所研究区域内任一点, Z(x)是该点的测量值,在所研究的区域内总共有n个实测点,即x1,x2,...,xn ,那么,对于任意待估点或待估块段V的实测值Zv(x) ,其估计值 是通过该待估点或待估块段影响范围内的 n个有效样本值 的线性组合来表示,即
式中, 为权重系数,是各已知样本在Z(xi) 在估计 时影响大小的系数,而估计 的好坏主要取决于怎样计算或选择权重系数 。 在求取权重系数时必须满足两个条件,一是使 的估计是无偏的,即偏差的数学期望为零;二是最优的,即使估计值 和实际值 Zv(x)之差的平方和最小,在数学上,这两个条件可表示为
五、普通克里格分析方法 设 Z(x)为区域化变量,满足二阶平稳和本征假设,其数学期望为m ,协方差
文档评论(0)