矿床统计预测-实习4a趋势面分析.doc

实习三 趋势面分析 目的要求： 趋势面分析是用一定的函数对地质体的某种特征在空间上的分布进行分析。用函数所代表的面来逼近（或拟合）该特征的趋势变化（或区域背景）。也就是说，用数学的方法，把观测值划分为两部分：趋势部分和偏差部分。趋势部分反映了区域性的总变化，受大范围的系统性因素控制。偏差反映局部范围的变化特点。受局部因素和随机因素控制。 为适应手算，本实习将通过二元一次多项式趋势函数计算，基本掌握趋势面分析的计算原理和方法步骤。 实习资料： 某地一条含金石英脉，用钻孔揭穿得20个矿体底板高数据（表11-1），通过趋势面分析，求得含金石英脉的总体产状及局部产状变化特征。若结合厚度、品位等资料，则可进一步研究它们之间的关系。 方法步骤： 二元一次多项式趋势函数的计算： 1．整理原始观测值（数据见计算表11-1）。其中x为横坐标，y为纵坐标，（x，y为相对值），z为观测值即矿体底板标高，观测点要尽量均匀，可以是非网格分布。 2．求趋势面方程，二元一次方程为： （11-1） 其中a0a1a2为待定系数，用最小二乘法在满足观测值（zi）和趋势值（）的偏差平方和为最小的条件下，求得： 令：偏差平方和 （11-2） 把（11-1）代入（11-2）得： （11-3） 为了得到最佳的拟合趋势面，要求达到最小。为此，分别求（11-3）式中对a0、a1、a2的偏导数，并令其等于零，得： 表11-1 二维一次趋势面计算表 观测点序号10.53.54021.53.56032.54.08144.03.58550.5.2.52061.03.03572.02.53582.53.05593.52.560104.02.570110.51.515121.02.020132.01.545143.02.035153.51.545161.00.55171.51.025182.00.515193.03.025204.04.030Σ43.548801 整理后得： 将计算表11-1中所得有关计算结果代入，解联立方程，即可求得系数a0、a1、a2，联立方程为 解得： a0= a1= a2= 所以，求得的二维一次趋势面方程为： = 3．求出各点的趋势值，将计算结果填于表11-1。 4．求出各点的偏差（乘余）值。=－，即各点的观测值减去该点的趋势值，计算结果缜于表11-1。 5．求拟合系数C：由于对一组观测点做趋势面分析时，不同次数的趋势面对原始观测值的拟合（逼近）程度不同，为了衡量拟合程度，可用离差平方和的变化表示，其公式为： = （11-4） 式中： ——观测值。 ——趋势值。 ——所有观测值的算术平均值。 计算结果填于表11-1。 6．据表11-1计算的趋势值及偏差值分别作观测值等值线图（图11-1），二维一次趋势面图（图11-2）及偏差图（图11-3）。 7．对结果进行分析。 图 11-2 二次趋势面 C = % 图 11-1 观测值等值线 图 11-3 二次剩余值图 图 11-2 二次趋势面 C = %  MATLAB计算程序 导入数据文件：data.txt % 输入变量 x = data(:,2); y = data(:,3); Z = data(:,4); 调用trend23.m计算 3. 运行trendFig.m绘图