插值与曲线拟合实验报告

《插值与曲线拟合实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.数值计算方法插值与拟合实验报告摘要：通过实验介绍插值方法中常见的拉格朗日插值，线性分段插值和牛顿前插公式，分析计算各种方法的插值余项。在曲线拟合方面使用两种不同类型的曲线来拟合同一组数据，并计算残差向量范数，比较不同曲线拟合的效果，在此例上给出优劣的判断。关键词：拉格朗日插值；线性分段插值；牛顿前插公式；曲线拟合引言在工程和科学计算中经常碰到只知道离散的数据测量点而需要匹配其变量之间的数学函数表达式的情况，这就需要插值和拟合的数值方法来解决这些问题。插值法是在离散数据的基础上补插连续函数，使得这条连续曲

2、线通过全部给定的离散数据点，也是离散函数逼近的重要方法，利用它可通过函数在有限点处的取值状况，估算出函数在其他点处的近似值。曲线拟合则是用连续曲线近似地刻画或比拟平面上离散点所表示的坐标之间的函数关系，在各个方面也有着愈加广泛的应用。1算法介绍1.1拉格朗日插值法1.1.1算法理论对某个多项式函数，已知有给定的k+1个取值点：其中对应着自变量的位置，而对应着函数在这个位置的取值。假设任意两个不同的xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：其中每个为拉格朗日基本多项式（或称插值基

3、函数），其表达式为：...拉格朗日基本多项式的特点是在上取值为1，在其它的点上取值为0。对于给定的个点：，拉格朗日插值法的思路是找到一个在一点取值为1，而在其他点取值都是0的多项式。这样，多项式在点取值为，而在其他点取值都是0。而多项式就可以满足在其它点取值为0的多项式容易找到，例如：它在点取值为：。由于已经假定两两互不相同，因此上面的取值不等于0。于是，将多项式除以这个取值，就得到一个满足“在取值为1，而在其他点取值都是0的多项式”：这就是拉格朗日基本多项式。拉格朗日插值法的公式结构整齐紧凑，在理论分析

4、中十分方便，然而在计算中，当插值点增加或减少一个时，所对应的基本多项式就需要全部重新计算，于是整个公式都会变化，非常繁琐。此外，当插值点比较多的时候，拉格朗日插值多项式的次数可能会很高，因此具有数值不稳定的特点，也就是说尽管在已知的几个点取到给定的数值，但在附近却会和“实际上”的值之间有很大的偏差。这类现象也被称为龙格现象，解决的办法是分段用较低次数的插值多项式。1.2.2算法描述对于已给定的点和待估计的点的横坐标x，如上述理论，将其值代入计算出插值基函数的值，然后根据公式：...计算出纵坐标的估计值，由

5、此完成对该点的插值过程，其中k为该点插值的阶数。1.2线性分段插值1.2.1算法理论首先介绍线型插值，假设我们已知坐标与要得到区间内某一位置 x 在直线上的值。根据图中所示，我们得到由于 x 值已知，所以可以从公式得到y的值已知求的过程与以上过程相同，只是与要进行交换。线性插值经常用于已知函数 在两点的值要近似获得其它点数值的方法，这种近似方法的误差定义为：其中 p 表示上面定义的线性插值多项式根据罗尔定理，我们可以证明：如果 f 有二阶连续导数，那么误差范围是正如所看到的，函数上两点之间的近似随着所近似

6、的函数的二阶导数的增大而逐渐变差。从直观上来看也是这样：函数的曲率越大，简单线性插值近似的误差也越大。...而分段线型插值就是把需要插值的区间依已给定的点分为k-1段，每段利用其端点进行线型插值。1.2.2算法描述利用已给定的点对插值区间分为段，将每段的端点与作为数据点利用公式在所构成的区间进行线性插值。1.3牛顿插值法1.3.1算法理论注意在下面的讲述中着重介绍差商的概念，设：问题是如何根据插值条件来计算待定系数由知，。由因而，其中，称为函数在点的一阶差商。由知：因而...其中，称为函数在点的二阶差商。

7、实际上，它是一阶差商的差商。一般地，如果已知一阶差商，那么就可以计算二阶差商类似于上述过程不断地推导下去，可得其中，称为函数f(x)在相应点处的n阶差商。这些结果说明，要保证上面构造多项式的方法具有期望的优点，就必须要有一个好的计算系数的方法。按上述方式构造插值多项式的方法叫做牛顿插值法。根据插值多项式的惟一性知，其截断误差与拉格朗日插值法相同，即：但也可以表示成差商形式。这是因为以为节点的多项式从而于是的截断误差可表为顺便指出，因为牛顿插值多项式具有性质：所以，类似于逐次线性插值法，也可以把上述和式中的

8、第二项看成是估计的一种实用误差估计式。...与差商概念密切联系的另一个概念是差分，它是指在等距节点上函数值的差。所谓等距节点，是指对给定的常数h（称为步长），节点称为处的一阶向前差分；称为处的一阶向后差分。一阶差分的差分称为二阶差分，即称为处的二阶向前差分。一般地，m阶向前和向后差分可定义如下：差商与差分之间存在关系：其中是介于之间的某个数。将此式代入，即可得到牛顿前插公式。1.3.2算法描述利用已给定的点和待估计的点的横坐标