4变异函数结构分析
4.1变异函数理论模型
4.1.1有基态值模型
球状模型
γ ( h ) = { 0 h = 0 C 0 + C ( 3 h 2 a − 1 h 3 2 a 3 ) 0 < h ≤ a C 0 + C h > a \gamma(h)=\left\{\begin{array}{lr} 0 & h=0 \\ C_{0}+C\left(\frac{3 h}{2 a}-\frac{1 h^{3}}{2 a^{3}}\right) & 0<h \leq a \\ C_{0}+C & h>a \end{array}\right. γ(h)=⎩⎪⎨⎪⎧0C0+C(2a3h−2a31h3)C0+Ch=00<h≤ah>a
C 0 , C 0 + C , C , a C_0,C_0+C,C,a C0,C0+C,C,a分别为块金常数、基台值、拱高、变程,其中球状模型是应用最广的模型,一般只用于一维,二维,和三维数据,最适合用在当空间相关性随步长增加而线性递减的情况(即曲线的斜率是线性递减的)。
指数模型:
γ ( h ) = { 0 h = 0 C 0 + C ( 1 − e − h a ) h > 0 \gamma(h)=\left\{\begin{array}{ll} 0 & h=0 \\ C_{0}+C\left(1-e^{-\frac{h}{a}}\right) & h>0 \end{array}\right. γ(h)={
0C0+C(1−e−ah)h=0h>0
C 0 , C 0 + C , C , a C_0,C_0+C,C,a C0,C0+C,C,a分别为块金常数、基台值、拱高、该模型在原点处的切线和基台值相交时的所对应的步长,其中 3 a 3a 3a约等于变程(对应于0.95*基台值)
适合一维或多维数据,由于指数模型在较大步长时也不能等于先验方差(相关系数也不能等于0),因此其变程定义为当变异函数为0.95 σ 2 σ_2 σ2( σ 2 σ_2 σ2为先验方差)所对应的步长
高斯模型
γ ( h ) = { 0 h = 0 C 0 + C ( 1 − e − h 2 a 2 ) h > 0 \gamma(h)=\left\{\begin{array}{ll} 0 & h=0 \\ C_{0}+C\left(1-e^{-\frac{h^2}{a^2}}\right) & h>0 \end{array}\right. γ(h)={
0C0+C(1−e−a2h2)h=0h>0
C 0 , C 0 + C , C , 3 a C_0,C_0+C,C,\sqrt{3}a C0,C0+C,C,3a分别为块金常数、基台值、拱高、约为变程,适合一维或多维数据,因此其变程定义为当变异函数为0.95 σ 2 σ_2 σ2( σ 2 σ_2 σ2为先验方差)所对应的步长一般适用于当数据在短距离内呈现高度连续性的情况
三种常用模型比较
对于球状模型与指数模型来说,插值结果表面比较粗糙,而对于高斯模型来说,插值结果较为平滑,是因为高斯模型的空间相关性随步长先增加后降低,连续性较强。
线性有基态值模型
γ ( h ) = { 0 h = 0 C 0 + A h 0 < h < a C 0 + C h = a \gamma(h)=\left\{\begin{array}{lc} 0 & h=0 \\ C_{0}+A h & 0<h<a \\ C_{0}+C & h=a \end{array}\right. γ(h)=⎩⎨⎧0C0+AhC0+Ch=00<h<ah=a
C 0 , C 0 + C , C , A C_0,C_0+C,C,A C0,C0+C,C,A分别为块金常数、基台值、拱高、斜率
纯块金效应模型
γ ( h ) = { 0 h = 0 C 0 h > 0 \gamma(h)=\left\{\begin{array}{ll} 0 & h=0 \\ C_{0} & h>0 \end{array}\right. γ(h)={
0C0h=0h>0
4.1.2无基态值模型
线性无基态值模型
γ ( h ) = { 0 h = 0 C 0 + A h h > 0 \gamma(h)=\left\{\begin{array}{ll} 0 & h=0 \\ C_{0}+Ah & h>0 \end{array}\right. γ(h)={
0C0+Ahh=0h>0
幂函数值模型
γ ( h ) = A h θ 1 < θ < 2 \gamma(h)=Ah^{\theta} \ 1<\theta<2 γ(h)=Ahθ 1<θ<2
4.1.3孔穴效应模型
变异函数在h大于一定距离后并非单调递增,而是以一定的周期b进行波动,此时变异函数曲线显示出一种”孔穴效应(hole effect)”.
孔穴效应模型包括有基台值模型和无基台值模型
有基台值
孔穴效应模型一般用于当数据呈现出一定周期性的情况,孔穴效应模型的变程为:当变异函数的值等于先验方差时所对应的最短步长
γ ( h ) = { 0 h = 0 C 0 + C ( 1 − s i n ( π h / r ) π h r ) h > 0 \gamma(h)=\left\{\begin{array}{ll} 0 & h=0 \\ C_0+C(1-\frac{sin(\pi h/r)}{\pi hr}) & h>0 \end{array}\right. γ(h)={
0C0+C(1−πhrsin