插值算法的详细描述.docVIP

PAGE  PAGE 13 重庆科技学院 数值分析上机实验报告 学 院： 数理学院 专 业： 数学与应用数学 班 级： 2011-02 学 号： 2011444839 学生姓名： 魏彩凤 指导教师： 唐利明 2013年 7月 摘要 各种岩石和矿物的磁磁性是不同的，测定地面上各部位的磁力强弱以研究地下岩石矿物的分布和地质构造，称做磁力勘探。由于地球本身就是个大磁性，所以对磁力的预测值应进行校正，求出只与岩石矿物磁性有关的磁力异常。 科学研究与工程技术中，常会遇到函数表达式过于复杂而不便于计算，且又需要计算众多点处的函数值；又或者是函数表达式未知，而仅由实验或测量得到这一函数在区间[a,b]中的n+1个节点上的值需要构造一个简单函数作为的近似表达式，使得成立，这类问题就称为插值问题，称为关于节点的插值函数。 现有新疆某石油勘探区前期地球物理勘探中的部分磁异常调查数据，其中第一列数据表示测线号，第二列号数据表示桩号，第三列数据为大地坐标系的X坐标，第四列数据为大地坐标系的Y坐标，第五列为经过基站校正的磁异常数据值。 由于地形较为恶劣，实际野外工作未完成标准网格的测量工作，请按照现有测量点的数据，自行选择两种不同数值分析插值法，在现有范围内扩展加密成标准的112*688网格数据，并用matlab绘制该区磁异常平面分布图，指出变化最快的地方。 关键字：磁异常 插值问题 matlab画图 网格数据 问题重述 图表1数据是新疆某石油勘探区前期地球物理勘探中的部分磁异常调查数据，其中第一列数据表示测线号，第二列号数据表示桩号，第三列数据为大地坐标系的X坐标，第四列数据为大地坐标系的Y坐标，第五列为经过基站校正的磁异常数据值。 （1）由于地形较为恶劣，实际野外工作未完成标准网格的测量工作，请按照现有测量点的数据，自行选择两种不同数值分析插值法，在现有范围内扩展加密成标准的112*688网格数据，并用matlab绘制该区磁异常平面分布图，指出变化最快的地方 （2）要求提交一份报告，报告内容的组成主要包括：所采用的数值分析差值法的详细介绍，附上相应的代码，该区磁异常平面分布图（在此图上说明高值区域、低值区域，并用红色粗线指出变化最快的地方, 通常情况下，低值区对应为凹陷，高值区对应为隆起，剧变带为断裂带）。 问题分析 由于条件的限制，无法测量到想要的数据，但是可以通过已测得的数据，建立一个模型，构造一个插值函数，通过这个函数，可以得到该地区磁异常的大致分布情况，由此可以通过这个插值函数，扩展成更丰富的数据情况，通过扩展可以得到想要点的磁异常情况，从而判断是否有石油。 通过插值的方法可以完成由于地形较为恶劣，实际野外工作未完成标准网格的测量工作，这只是其中的一个例子，通常类似的问题都可以用插值的问题，用可以测量到得的大量数据来推测出其该问题的扩展领域里的情况，从而得到之前无法测得的数据，在实际运用中，是非常有用的，使得插值的方法更具有实际应用的价值。 现在的主要问题有： 通过插值算法得到未知点的数据，扩展为标准的112*668的网格数据。 根据磁异常数据，运用matlab绘制磁异常平面分布图。 模型分析 模型中给出的数据量非常大，所以要运用matlab库函数xlsread调用Excel中的数据，要绘制网格图就要用到mesh函数，通常画图且是插值问题的用griddata函数。 模型建立 插值法的描述 科学研究与工程技术中，常会遇到函数表达式过于复杂而不便于计算，且又需要计算众多点处的函数值；又或者是函数表达式未知，而仅由实验或测量得到这一函数在区间[a,b]中的n+1个节点上的值需要构造一个简单函数作为的近似表达式，使得成立，这类问题就称为插值问题，称为关于节点的插值函数。 一、拉格朗日插值 1、线性插值与抛物线插值 对给定的插值点求得形如的插值多项式可以有各种不同方法，下面先讨论n=1的简单情形，假定给定的区间及端点函数值，要求线性插值多项式，使它满足。的几何意义就是通过两点与的直线，的表达式可由几何意义直接给出 由两点式看出，是由两个线性函数 线性组合得到的，其系数分别为及，即 显然，及也是线性插值多项式，在节点及上分别满足条件 称函数及为线性插值基函数。 下面讨论n=2的情况，假定插值节点为要求二次插值多项式，使它满足 我们知道在几何意义上就是通过三点的