插值多项式.docVIP

插值与逼近 一 插值多项式 有时候我们只知道函数f(x)在区间［a,b］，而不知道它在区间［a,b］的表达式十分复杂，不便于利用函数的表达式研究问题。插值法就是构造插值函数去近似被插值函数，使之满足插值条件 。 通常我们构造插值多项式。插值多项式就是利用一些已知的函数值所做的既能反映原来函数的主要性质，又有简单形式的一种较好的替代函数。 求插值多项式的基本思想： 设函数在区间［a,b］上个互不相同的点处的值。如果多项式在点上满足 则称是函数的插值多项式。 在本章中讨论拉格朗日插值多项式、牛顿插值多项式、埃尔米特插值多项式和分段插值多项式。 1. 拉格朗日插值多项式 拉格朗日插值法是最基本、最常用的插值方法，也是其他插值方法的基础。 我们讲授的拉格朗日插值多项式包括线性插值多项式、抛物线插值多项式和n次插值多项式 拉格朗日插值多项式的公式为： 其中基函数的公式为： 余项公式为 拉格朗日插值多项式计算步骤: ⑴　准确计算插值基函数。并化简 ⑵ 代入拉格朗日插值多项式公式正确求出插值多项式。 ⑶ 求出结果后，可以用进行验算。 ⑷ 根据余项公式进行误差估计。如果要估计误差，须知道函数的表达式。 例1. 令 x0=0,x1=1，求的线性拉格朗日插值多项式，并估计误差。 【思路】参见上述拉格朗日插值多项式计算步骤 解 线性插值基函数为拉格朗日插值多项式为 拉格朗日插值多项式的误差为 ２.牛顿插值多项式 牛顿插值多项式是一种重要的计算插值多项式的方法.在学习时,应掌握节点数较少的牛顿插值多项式的计算. 牛顿插值多项式公式为 其中k阶差商的计算公式为： 牛顿插值多项式的余项公式为 牛顿插值多项式的计算步骤: （1）利用k阶差商的计算公式准确计算各阶差商。并化简 （2）代入牛顿插值多项式公式正确计算出插值多项式。 （3）求出结果后，可以用进行验算。 在计算插值多项式时, 牛顿插值多项式的计算量比拉格朗日插值多项式的计算量要小。因为应用牛顿插值多项式时，可以避免拉格朗日插值多项式在计算时,每增加或改变节点时需要重新计算插值基函数的缺陷，而只需要在已知的多项式基础上增加一项即可。(参见下例3) 差商的性质(参见课本p28-29) 例2．已知函数表 x 0 1 2 y 3 6 11 求抛物插值多项式，并计算的近似值。 【思路】1.列表计算k阶差商（根据所给互异节点的个数确定） 2.代入公式求出牛顿插值多项式 解 列差商表 将表中最上斜行中函数值与各阶差商值代入二阶牛顿插值多项式，得 一阶差商 二阶差商 0 1 2 3 6 11 3 5 1 令，得 说明：求抛物插值多项式可以用多种方法，一般采用牛顿插值多项式求解。 3.埃尔米特插值多项式 埃尔米特插值多项式又称为带导数的插值多项式,即在节点处既要求函数值和已知函数值相等又要求导数值和已知导数值相等,且埃尔米特插值多项式的精度较高. ⑴两点三次埃尔米特插值多项式公式为： ⑵一般埃尔米特插值多项式公式为： 其中和是拉格朗日插值多项式的基函数及其导数埃尔米特插值 埃尔米特插值余项公式为： 例3 求满足条件（如表所示） x 1 2 y 2 3 y/ 1 -1 试用两点三次埃尔米特插值法求埃尔米特插值多项式. 【思路】根据所给条件代入两点三次埃尔米特插值多项式公式。 解 4.分段插值多项式 （1）分段线性插值多项式公式 分段线性插值多项式余项公式 ，其中 （2）分段三次埃尔米特插值多项式 分段三次埃尔米特插值多项式公式为 分段三次埃尔米特插值多项式余项公式为 ，其中插值多项式 二、最小二乘法 构造简单函数表达式的另一种方法是曲线拟合。在自然科学和工程技术中，常通过测量和试验得到函数在一些点上的函数值，一般来说这些数据的数目都比较多，如果用插值多项式的方法构造函数表达式往往是高次插值多项式，计算量大且使用不便。另外，这些数据是通过测量和试验所得，数据本身存在着误差，而用插值多项式构造的函数表达式要求通过这些函数值，因此，保留了所有的测量误差，使计算结果与实际问题相差较大。而对数据拟合问题来说，并不要求构造的函数表达式都通过所有的已知点，而是要求与已知点的误差尽量小。最小二乘法就是数据