结构工程插值法综述《计算方法》学习报告.docVIP

结构工程插值法综述 一、插值法的由来 1.插值法简介在许多实际问题及科学研究中，因素之间往往存在着函数关系，然而，这种关系经常很难有明显的解析表达，通常只是由观察与测试得到一些离散数值。有时，即使给出了解析表达式，却由于表达式过于复杂，不仅使用不便，而且不易于进行计算与理论分析。解决这类问题的方法有两种：一种是插值法,另一种是拟合法。插值法是一种古老的数学方法，它来自生产实践，早在一千多年前，我国科学家在研究历法上就 应用了线性插值与二次插值，但它的基本理论却是在微积分产生之后才逐渐完善的，其应用也日益增多，特别是在计算机软件中，许多库函数，如sinx ,cosx 等的计算实际上归结于它的逼近函数的计算。逼近函数一般为只含有算术运算的简单函数，如多项式、有理分式（即多项式的商）。在工程实际问题当中，我们也经常会碰到诸如此类的函数值计算问题。被计算的函数有时不容易直接计算，如表达式过于复杂或者只能通过某种手段获取该函数在某些点处的函数值信息或者导数值信息等。因此，我们希望能用一个“简单函数”逼近被计算函数，然后用该简单函数的函数值近似替代被计算函数的函数值。这种方法就叫插值逼近或者插值法。插值法要求给出函数的一个函数表，然后选定一种简单的函数形式，比如多项式、分段线性函数及三角多项式等，通过已知的函数表来确定一个简单的函数g(x)作为f(x)的近似，概括地说，就是用简单函数为离散数组建立连续模型。 二、插值法的描述 插值法又称“内插法”，是利用函数f (x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值，这种方法称为插值法。如果这特定函数是多项式，就称它为插值多项式。 插值法是函数逼近的一种重要方法，是数值计算的基本课题。本节只讨论具有唯一插值函数的多项式插值和分段多项式插值，对其中的多项式插值主要讨论n次多项式插值的方法，即给定n+1各点处的函数值后，怎样构造一个n次插值多项式的方法。虽然理论上可以用解方程组（2）（那里m=n）得到所求插值多项式，但遗憾的是方程组（2）当n较大时往往是严重是病态的。故不能用解方程组的方法获得插值多项式。本节介绍的内容有：lagrange插值，newton插值，hermite插值，分段多项式插值及样条插值。 三、插值法的相关新概念 1、多元回归法 多元回归被用来确定你的数据的大规模的趋势和图案。你可以用几个选 项来确定你需要的趋势面类型。多元回归实际上不是插值器，因为它并不试图预测未知 的 Z 值。它实际上是一个趋势面分析作图程序。 使用多元回归法时要涉及到曲面定义 和指定 XY 的最高方次设置，曲面定义是选择采用的数据的多项式类型，这些类型分别是简单平面、双线性鞍、二次曲面、三次曲面和用户定义的多项式。参数设置是指定多项式方程中 X 和 Y 组元的最高方次 2、反距离倒数乘方法 距离倒数乘方格网化方法是一个加权平均插值法，可以进行 确切的或者圆滑的方式插值。方次参数控制着权系数如何随着离开一个格网结点距离的增加而下降。对于一个较大的方次，较近的数据点被给定一个较高的权重份额，对于一个较小的方次，权重比较均匀地分配给各数据点。 计算一个格网结点时给予一个特定数据点的权值与指定方次的从结点到观测点的该结点被赋予距离倒数成比例。当计算一个格网结点时，配给的权重是一个分数，所有权重的总和等于 1.0。当一个观测点与一个格网结点重合时，该观测点被给予一个实际为 1.0 的权重，所有其它观测点被给予一个几乎为 0.0的权重。换言之，该结点被赋给与观测点一致的值。这就是一个准确插值。 距离倒数法的特征之一是要在格网区域内产生围绕观测点位置的牛眼。用距离倒数格网化时可以指定一个圆滑参数。大于零的圆滑参数保证，对于一个特定的结点，没有哪个观测点被赋予全部的权值，即使观测点与该结点重合也是如此。圆滑参数通过修匀已被插值的格网来降低牛眼影响。 3、克里金法 克里金法是一种在许多领域都很有用的地质统计格网化方法。克里金 法试图那样表示隐含在你的数据中的趋势，例如，高点会是沿一个脊连接，而不是被牛眼形等值线所孤立。 克里金法中包含了几个因子：变化图模型，漂移类型 和矿块效应。 4、最小曲率法 最小曲率法广泛用于地球科学。用最小曲率法生成的插值面类似于一个 通过各个数据值的，具有最小弯曲量的长条形薄弹性片。最小曲率法，试图在尽可能严 格地尊重数据的同时，生成尽可能圆滑的曲面。 使用最小曲率法时要涉及到两个参数： 最大残差参数和最大循环次数参数来控制最小曲率的收敛标准 5、径向基本函数法 径向基本函数法是多个数据插值方法的组合。