《地球物理反演》读书报告.docVIP

《地球物理反演》 ——读书报告 专业所在院（系、部）核工程技术学院 研究生姓名 郭猛猛 学 号 2010070807 专业名称 固体地球物理学 日 期 2011年6月30日 《地球物理反演》读书报告 拉格朗日插值法 在数值分析中，拉格朗日插值法是以法国十八世纪数学家约瑟夫·路易斯·拉格朗日命名的一种多项式插值方法。许多实际问题中都用函数来表示某种内在联系或规律，而不少函数都只能通过实验和观测来了解。如对实践中的某个物理量进行观测，在若干个不同的地方得到相应的观测值，拉格朗日插值法可以找到一个多项式，其恰好在各个观测的点取到观测到的值。这样的多项式称为拉格朗日（插值）多项式。数学上来说，拉格朗日插值法可以给出一个恰好穿过二维平面上若干个已知点的多项式函数。拉格朗日插值法最早被英国数学家爱德华·华林于1779年发现，不久后（1783年）由莱昂哈德·欧拉再次发现。1795年，拉格朗日在其著作《师范学校数学基础教程》中发表了这个插值方法，从此他的名字就和这个方法联系在一起。 定义： 对于给定的若n+1个点，对应于它们的次数不超过n的拉格朗日多项式只有一个。如果计入次数更高的多项式，则有无穷个，因为所有与相差的多项式都满足条件。 对某个多项式函数，已知有给定的k?+?1个取值点： 其中xj对应着自变量的位置，而yj对应着函数在这个位置的取值。 假设任意两个不同的xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为： 其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为： 拉格朗日基本多项式的特点是在xj 上取值为1，在其它的点 上取值为0。 证明： 存在性： 对于给定的k+1个点：，拉格朗日插值法的思路是找到一个在一点xj取值为1，而在其他点取值都是0的多项式。这样，多项式在点xj取值为yj，而在其他点取值都是0。而多项式就可以满足 在其它点取值为0的多项式容易找到，例如： 它在点xj取值为：。由于已经假定xi两两互不相同，因此上面的取值不等于0。于是，将多项式除以这个取值，就得到一个满足“在xj取值为1，而在其他点取值都是0的多项式”： 这就是拉格朗日基本多项式。 唯一性： 次数不超过k 的拉格朗日多项式至多只有一个，因为对任意两个次数不超过k 的拉格朗日多项式：P1和P2，它们的差P1 ? P2在所有k+1个点上取值都是0，因此必然是多项式的倍数。因此，如果这个差P1 ? P2不等于0，次数就一定不小于k+1。但是P1 ? P2是两个次数不超过k 的多项式之差，它的次数也不超过k。所以P1 ? P2 = 0，也就是说P1 = P2。这样就证明了唯一性。 优点与缺点： 拉格朗日插值法的公式结构整齐紧凑，在理论分析中十分方便，然而在计算中，当插值点增加或减少一个时，所对应的基本多项式就需要全部重新计算，于是整个公式都会变化，非常繁琐。这时可以用重心拉格朗日插值法或牛顿插值法来代替。此外，当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，因此具有数值不稳定的特点，也就是说尽管在已知的几个点取到给定的数值，但在附近却会和“实际上”的值之间有很大的偏差。这类现象也被称为龙格现象，解决的办法是分段用较低次数的插值多项式。 牛顿迭代法 牛顿迭代法（Newtons method）又称为牛顿-拉夫逊方法（Newton-Raphson method），它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根,此时线性收敛，但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。设r是f(x) = 0的根，选取x0作为r初始近似值，过点（x0,f(x0)）做曲线y = f(x)的切线L，L的方程为y = f(x0)+f(x0)(x-x0)，求出L与x轴交点的横坐标 x1 = x0-f(x0)/f(x0)，称x1为r的一次近似值。过点（x1,f(x1)）做曲线y = f(x)的切线，并求该切线与x轴交点的横坐标 x2 = x1-f(x1)/f(x1)，称x2为r的二次近似值。重复以上过程，得r的近似值序列，其中x(n+1)=x(n)－f(x(n))/f(x(n))，称为r的n+1次近似值，上式称为牛顿迭代公式。 解非线性方程f(x)=0的牛顿法是把非线性方程线性化的一种近似方法。把