第5讲工程数学摘要.ppt

* 计算实例 例如：假设某地区降水量Z(x)（单位：mm）是二维区域化随机变量，满足二阶平稳假设，其观测值的空间正方形网格数据如图4.2.1所示（点与点之间的距离为h=1 km）。试计算其南北方向及西北和东南方向的变异函数。 计算实例 图4.2.1 空间正方形网格数据（点间距h=1 km） 从图4.2.1可以看出，空间上有些点，由于某种原因没有采集到。如果没有缺失值，可直接对正方形网格数据结构计算变异函数；在有缺失值的情况下，也可以计算变异函数。只要“跳过”缺失点位置即可（图4.2.2）。 首先计算南北方向上的变异函数值，由变异函数的计算公式可得 =385/72=5.35 图4.2.2 缺失值情况下样本数对的组成和计算过程☉为缺失值 同样计算出 最后，得到南北方向和西北—东南方向上的变异函数计算结果见下表。同样可以计算东西方向上的变异函数。 方向 南北 方向 西北—东南 h 1 2 3 4 5 h 1.41 2.82 4.24 5.65 7.07 N(h) 36 27 21 13 5 N(h) 32 21 13 8 2 5.35 9.26 17.55 25.69 22.90 7.06 12.95 30.85 58.13 50.00 变异函数的参数 变异函数有４个非常重要的参数，即基台值（sill）、变程（range）或称空间依赖范围（range of spatial dependence）、块金值（nugget）或称区域不连续性值（localized discontinuity）和分维数（fractal dimension）。 前3个参数可以直接从变异函数图中得到。它们决定变异函数的形状与结构。 变异函数的形状反映自然现象空间分布结构或空间相关的类型，同时还能给出这种空间相关的范围。 当变异函数随着间隔距离h的增大，从非零值达到一个相对稳定的常数时，该常数称为基台值C0+C。 当间隔距离h=0时，γ(0)= C0，该值称为块金值或块金方差（nugget variance）。 基台值是系统或系统属性中最大的变异，变异函数达到基台值时的间隔距离a称为变程。变程表示在h≥a以后，区域化变量Z(x)空间相关性消失。 块金值表示区域化变量在小于抽样尺度时非连续变异，由区域化变量的属性或测量误差决定。 C0 块金常数 a 变程 C 拱高 上述３个参数可从变异函数曲线图直接得到，或通过估计曲线回归参数得到。 第4个参数，即分维数用于表示变异函数的特性，由变异函数 和间隔距离h之间的关系确定 分维数D为双对数直线回归方程中的斜率，它是一个无量纲数。分维数D的大小，表示变异函数曲线的曲率，可以作为随机变异的量度。但该随机分维数D与形状分维数有本质的不同。 * 计算实例（协方差与半变异函数的实例计算） x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 z(x1)=4,z(x2)=3,z(x3)=4,z(x4)=5,z(x5)=7,z(x6)=9,z(x7)=7,z(x8)=8,z(x9)=7,z(x10)=7 * * * * * * * * * * * * * * * 计算实例（协方差与半变异函数的实例计算） * 计算实例 * 计算实例 * 计算实例 * 计算实例 * 计算实例 * 变异分析 半变异函数和协方差函数把统计相关系数的大小作为一个距离的函数，是地理学相近相似定理的定量化。下图为典型的半变异函数图和其对应的协方差函数图。 * 变异分析 * 变异分析 图1和图2显示，半变异值的变化随着距离的增大而增加，协方差随着距离的增大而减小。这主要是由于半变异函数和协方差函数都是事物空间相关系数的表现，当两事物彼此距离较小时，它们是相似的，因此协方差值较大，而半变异值较小；反之，协方差值较小，而半变异值较大。此外，协方差函数和半变异函数随着距离的加大基本呈反向变化特征 * 变异分析 半变异函数曲线和协方差函数曲线反映了一个采样点与其相邻采样点的空间关系。此外，它们对异常采样点具有很好的探测作用，在ArcGIS地统计分析模块中可以使用两者的任意一个，一般采用半变异函数。在半变异曲线图中有两个非常重要的点：间隔为0时的点和半变异函数趋近平稳时的拐点，由这两个点产生四个相应的参数：块金值（Nugget）、