 ——估值权系数
 ——估计值
②克里格方程组：
 （6）
其中：n——用于估值的观测点的数目
 ——观测点坐标（  ）
 ——待估点坐标（x,y）
 ——两点间的协方差函数
 ——求条件机制的拉格朗日乘数
——克里格估值权系数
③点估值克里格方差：
 （7）
其中： ——变差函数基台值
将其转换为矩阵形式，令：
则普通克立格方程组为： （8）
其中：= = （9）
 代表点与点j间的距离
解方程组，可得： 其估计方差为： （10）
将克里格方程组（式6）转换成矩阵形式（式8），求方程组的解就转换为求矩阵的解。本部分用高斯消去法求矩阵的解，将矩阵的解代入式5中，进而得到待估点的克里格值和估值方差。
本程序对二维区域内的点进行克里格估值。在对选取的点进行估值时，先对周围点进行扫描，取距离待估点最近的若干个点作为信息样品点参与估值。程序计算待估点的克里格值和估值方差。
5结论
5.1本软件以地学处理软件中克里格估值功能的实现为参考，对点克里格估值的相关模型进行编程，最终实现了区域内任意点的克里格估值。
5.2实现了将空间信息点的克里格估值分解为若干个主要模块，最终在实现点克里格估值的同时，用户可根据需要在不同阶段实现对软件运行产生的数据文档进行调用。
5.3本软件可以对具有空间相关性的地质属性进行估值，实现目标区域地质储量估算。
5.4本软件具有良好的交互性，界面简洁，有针对性，可以提高用户数据处理的效率。
6 参考文献：
[1] 侯景儒，尹镇南，李维明，等.空间信息统计学.北京：地质出版社，1998，45
[2] Armstrong M.Geostatistics-Proceedings of the International GeostatisticsCongress.France:Kluwer Academic Publishers,1958,Vol,1-2
[3] Matheron G.Thetheory of Variables and its Application.France:Centre of GeostatisticsFontaine,1970
[4] Journel AG.Nonparameteric estimation of spatial distribution. Math Geol,1983,445-468
[5] 肖斌，赵大鹏，侯景儒.地质统计学新进展.地球科学进展,2000,15(3),295
[6] 孙洪泉.地质统计学及其应用. 北京：中国矿业大学出版社，1990，39-41
[7] 孙洪泉.地质统计学及其应用. 北京：中国矿业大学出版社，1990，66-69
2/2 首页 上一页 1 2 |
