论文2__通过三维数据插值进行电阻率的拟合与估计教案分析.doc

摘 要 在实际问题中，我们无法对一个物体进行密集测量，一般等间隔的选取部分点进行测量，这时候就需要采用数学方法，根据已知位置的数据计算未知位置的数据。本文根据题目要求建立了两种合理的差值加密模型——反距离加权插值模型、克里金插值模型，对文件中给出的坐标网格大小为10m*10m*10m的三维电阻率数据进行插值，得到坐标网格大小为1m*1m*1m的数据。我们借助多种数学软件，处理大量数据并对其进行可视化处理，绘制出颜色图反映插值效果，根据插值拟合特点对插值效果给出定量指标。 针对问题一，首先对所给出的三维的电阻率数据进行分析，根据不同空间插值方法的优势，结合已知数据的网格化特点，确定插值模型——反距离加权插值模型、克里金插值模型。反距离加权插值模型算法简单易于实现，但对权重函数的选择十分敏感，而网格化的数据恰好规避了这一劣势。克里金插值是在变异函数理论及结构分析基础上，进行无偏、最优估计的一种方法，不仅考虑了观测点与待估计点的相对位置，而且考虑了各观测点之间的相对关系，插值效果较好。利用两种方法分别计算出空间某点(45.8,-32.7,68.2)处的电阻率数值。反距离加权插值法结果为Z1=196.1972234615818；克里金插值结果为Z2=194.7184939362413。 针对问题二：利用问题一中两种插值方法，分别计算出网格大小为1m*1m*1m的三维电阻率数据（由于数据量过大，具体数据见电子版附件1、2），分析两种模型的计算流程可知，克里金算法计算过程较为复杂，数据量较大，计算时间与占用内存大，而反距离加权插值计算过程较为简单，计算时间明显短于克里金方法，但计算精度较低。将原网格数据及两种方法加密网格后数据利用Excel计算出平均值与标准差，见下表。 平均值E 方差S2 原始数据 197.268937 13 反距离加权插值数据 196.549905 15 克里金差值数据 197.647521 15 针对问题三：对插值加密网格后的数据，利用颜色图直观展示。将电阻率最小值置为纯蓝色，中间值置为纯绿色，最大值置为纯红色，中间数值采用过渡的颜色，利用matlab对切片Z=0,50分别给出原数据，两种方法加密网格后数据的颜色对比图。另外分别对切片X=82，Y=47, Z=88 ，给出两种加密方法得到数据的颜色图。通过同一切面颜色图的比较，直观反映出原数据、反距离加权插值加密、与克里金插值加密的插值效果。 针对问题四：利用插值加密前后电阻率均值的偏离程度以及电阻率沿x，y，z三个轴向变化趋势的偏离程度给出定量指标，定量反应两种插值方法的差值效果。 关键词：反距离加权插值 变差函数 加权系数 克里金插值 一、问题重述 在实际问题中，由于技术与测量成本等原因，只能等间隔的选取部分点进行测量。实际中需要采用数学方法获得更多位置数据，根据已知位置的数据计算未知位置的数据。 根据坐标网格大小为10m*10m*10m的空间三维体电阻率数据，通过插值加密后获得网格大小为1m*1m*1m的三维电阻率数据，必须保证插值后的电阻率数据极值不变，且坐标位置相同。另外，必须保证插值后的三维成像结果与插值前的成像结果形态基本一致，只是前者像素更高。 （1）请给出符合条件的两种计算方法，并给出相应数学公式，证明该方法插值后的数据极值与位置不变。分别计算出空间某点(45.8,-32.7,68.2)处的电阻率数值。 (2)利用（1）中两种方法分别计算出网格大小为1m*1m*1m的三维电阻率数据，给出计算流程，并对两种方法的计算量或计算复杂性进行评估。同时给出原网格数据及两种方法加密网格后数据的平均值与标准差。 (3)对加密网格后的直观效果，可采用颜色图展示。 对每一幅需要对比显示效果的图，请将最小值置为纯蓝色（RGB为(0,0,255)），中间值置为纯绿色（RGB为(0,255,0)），最大值置为纯红色（RGB为(255,0,0)），中间数值采用过渡的颜色，可自行设计。采用这种方法，请对切片Z=0,50分别给出原数据，两种方法加密网格后数据的颜色对比图。另外分别对切片X=82，Y=47, Z=88 ，请给出两种加密方法得到数据的颜色图。 （4）对不同的插值加密方法的效果，给出定量的指标，分别计算出两种不同加密方法得到数据的指标值，并给出评价。 二、模型假设 1、基本假设： 1)假设data3D_.txt所给电阻率数据背景值真实； 2)被测电阻各向同性； 3)不考虑环境因素对电阻率的影响； 2、针对反距离加权插值法： 未知值的点受较近控制点的影响比较远控制点的影响更大； 3、针对克里金插值法： 本征假设： 设Z（x）为空间点x处的属性值即电阻率，Z(x)与Z(x + h)之间的相关性不