《数值分析--清华李庆杨五版第二章_插值法.ppt》由会员分享，可在线阅读，更多相关《数值分析--清华李庆杨五版第二章_插值法.ppt（124页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1 1引言引言 问题的提出问题的提出函数解析式未知函数解析式未知,通过实验观测得到的一组数据通过实验观测得到的一组数据,即即在某个区间在某个区间a,b上给出一系列点的函数值上给出一系列点的函数值yi=f(xi)或者给出函数表或者给出函数表y=f(x)y=p(x)xx0 x1x2xnyy0y1y2yn第二章第二章 插值法插值法插值法的基本原理插值法的基本原理设函数设函数y=y=f(x)定义在区间定义在区间 a,b 上上,是是 a,b 上上 取取 定定 的的 n+1个个 互互 异异 节节 点点,且且 在在 这这 些些 点点 处处 的的 函函 数数 值值 为为已已知知 ,即即 若若存存在在一一个个f
2、(x)的近似函数的近似函数 ,满足满足则称则称 为为f(x)的一个的一个插值函数插值函数,f(x)为为被插函数被插函数,点点xi为为插值节点插值节点,称称(2.1)2.1)式为式为插值条件插值条件,而误差函数而误差函数R(x)=称为称为插值余项插值余项,区间区间 a,b 称为称为插值插值区间区间,插值点在插值区间内的称为插值点在插值区间内的称为内插内插,否则称否则称外插外插(2.1)2.1)插值函数插值函数 在在n+1个互异插值节点个互异插值节点 (i=0,1,n)处与处与 相等相等,在其它点在其它点x就用就用 的值作为的值作为f(x)的近似值。这一过程称为的近似值。这一过程称为插值插值，点，
3、点x称为插值点。换称为插值点。换句话说句话说,插值就是根据被插函数给出的函数表插值就是根据被插函数给出的函数表“插出插出”所要点的函数值。用所要点的函数值。用 的值作为的值作为f(x)的近似值的近似值,不仅希不仅希望望 能较好地逼近能较好地逼近f(x),而且还希望它计算简单而且还希望它计算简单。由于由于代数多项式具有数值计算和理论分析方便的优点。所以代数多项式具有数值计算和理论分析方便的优点。所以本章主要介绍代数插值。即求一个次数不超过本章主要介绍代数插值。即求一个次数不超过n n次的多项次的多项式。式。满足满足则称则称P(x)P(x)为为f(x)f(x)的的n n次插值多项式。这种插值法通常
4、称为次插值多项式。这种插值法通常称为代数插值法。其几何意义如下图所示代数插值法。其几何意义如下图所示 定理定理1n次代数插值问题的解是存在且惟一的次代数插值问题的解是存在且惟一的证明证明:设设n n次多项式次多项式是函数是函数 在区间在区间 a,ba,b上的上的n+1n+1个互异的节点个互异的节点 (i=0,1,2,i=0,1,2,n),n)上的插值多项式上的插值多项式,则求插值多项式则求插值多项式P(x)P(x)的问题就归结为求它的系数的问题就归结为求它的系数 (i=0,1,2,i=0,1,2,n),n)。由插值条件由插值条件:(:(i=0,1,2,i=0,1,2,n),n),可得可得 这是
5、一个关于待定参数这是一个关于待定参数 的的n+1阶线性方阶线性方程组程组,其系数矩阵行列式为其系数矩阵行列式为 称为称为Vandermonde（范德蒙）行列式，因范德蒙）行列式，因xixj（当当ij），），故故V0。根据解线性方程组的克莱姆根据解线性方程组的克莱姆（Gramer）法则，方程组的解法则，方程组的解 存在惟一，从而存在惟一，从而P(x)P(x)被惟一确定。被惟一确定。惟一性说明，不论用何种方法来构造，也不论用何种惟一性说明，不论用何种方法来构造，也不论用何种形式来表示插值多项式形式来表示插值多项式,只要满足插值条件只要满足插值条件(2.1)2.1)其结其结果都是相互恒等的。果都是相
6、互恒等的。2拉格朗日（拉格朗日（Lagrange）插值插值 为为了了构构造造满满足足插插值值条条件件(i=0,1,2,n)的便于使用的插值多项式的便于使用的插值多项式P(x),P(x),先考察几种简单情形先考察几种简单情形,然后再推广到一般形式。（然后再推广到一般形式。（线性插值与抛物插值）线性插值与抛物插值）（1）线性插值）线性插值线性插值是代数插值的最简单形式。假设给定了函数线性插值是代数插值的最简单形式。假设给定了函数f(x)f(x)在两个互异的点的值，在两个互异的点的值，,现要求用线性函数现要求用线性函数 近似地代替近似地代替f(x)f(x)。选选择参数择参数a和和b,使使 。称这样的
7、线性函数。称这样的线性函数P(x)P(x)为为f(x)f(x)的线性插值函数的线性插值函数。线性插值的几何意义线性插值的几何意义:用用通过点通过点 和和 的直线近似地代替曲线的直线近似地代替曲线 y=f(x)=f(x)由解析几何知道由解析几何知道,这条直线用点斜式表示为这条直线用点斜式表示为为了便于推广，记为了便于推广，记 这是一次函这是一次函数数,且有性质且有性质 与与 称为线性插值基函数。且有称为线性插值基函数。且有于是线性插值函数可以表示为与基函数的线性组合于是线性插值函数可以表示为与基函数的线性组合 例例2.1 2.1 已知已知 ,求求 解解:这里这里x0=100，y0=10，x1=1
8、21，y1=11,利用线性插值利用线性插值 （2 2）抛物插值抛物插值抛抛物物插插值值又又称称二二次次插插值值，它它也也是是常常用用的的代代数数插插值值之之一一。设设已已知知f(x)f(x)在在三三个个互互异异点点x0，x1，x2的的函函数数值值y0，y1，y2，要构造次数不超过二次的多项式要构造次数不超过二次的多项式使满足二次插值条件：使满足二次插值条件：这就是二次插值问题。其几何意义是用经过这就是二次插值问题。其几何意义是用经过3个点个点 的抛物线的抛物线 近似代替曲线近似代替曲线 ,如下图所示。因此也称之为抛物插值。如下图所示。因此也称之为抛物插值。P(x)的参数的参数直接由插值条件决定
9、，直接由插值条件决定，即即 满足下面满足下面的代数方程组：的代数方程组：该三元一该三元一次方程组次方程组的系数矩阵的系数矩阵的行列式是范德蒙行列式，当的行列式是范德蒙行列式，当 时，时，方程组的解唯一。方程组的解唯一。为了与下一节的为了与下一节的Lagrange插值公式比较插值公式比较,仿线性插值仿线性插值,用用基函数的方法求解方程组。先考察一个特殊的二次插值基函数的方法求解方程组。先考察一个特殊的二次插值问题：问题：求二次式求二次式 ,使其满足条件：使其满足条件：这个问题容易求解。由上式的后两个条件知这个问题容易求解。由上式的后两个条件知:是是 的两个零点。于是的两个零点。于是 再由另一条件
10、再由另一条件 确定系数确定系数 从而导出从而导出 类似地可以构造出满足条件：类似地可以构造出满足条件：的插值多项式的插值多项式及满足条件：及满足条件：的插值多项式的插值多项式这样构造出来的这样构造出来的 称为抛物插值的基函数称为抛物插值的基函数 取已知数据取已知数据 作为线性组合系数作为线性组合系数,将基函数将基函数 线性组合可得线性组合可得 容易看出容易看出,P(x)P(x)满足条件满足条件 拉格朗日插值多项式拉格朗日插值多项式 两个插值点可求出一次插值多项式两个插值点可求出一次插值多项式,而三而三个插值点可求出二次插值多项式。插值点增加到个插值点可求出二次插值多项式。插值点增加到n+1个时
11、个时,也就是通过也就是通过n+1个不同的已知点个不同的已知点,来构造一个次数为来构造一个次数为n的代数多项式的代数多项式P(x)。与推导抛物插与推导抛物插值的基函数类似值的基函数类似,先构造一个特殊先构造一个特殊n次多项式次多项式 的插的插值问题值问题,使其在各节点使其在各节点 上满足上满足即即 由条件由条件 ()()知知,都是都是n n次次 的零点的零点,故可设故可设 其中其中 为待定常数。由条件为待定常数。由条件 ,可求得可求得于是于是 代入上式代入上式,得得称称 为关于基点为关于基点 的的n n次插值基函数次插值基函数(i=0,1,i=0,1,n),n)以以n+1个个n次基本插值多项式次
12、基本插值多项式为基础为基础,就能直接写出满足插值条件就能直接写出满足插值条件的的n次代数插值多项式。次代数插值多项式。事实上，由于每个插值基函数事实上，由于每个插值基函数都是都是n次值多项式次值多项式,所以他们的线性组合所以他们的线性组合是次数不超过是次数不超过n n次的多项式次的多项式,称形如（称形如（2.8）式的插）式的插值多项式为值多项式为n次拉格朗日插值多项式。并记为次拉格朗日插值多项式。并记为 （2.8）例例2.2 已知已知y=f(x)的函数表的函数表 求线性插值多项式求线性插值多项式,并计算并计算x=1.5 的值的值X 1 3 y 1 2解解:由线性插值多项式公式得由线性插值多项式
13、公式得例例2.3已知已知x=1,4,9的平方根值的平方根值,用抛物插值公式用抛物插值公式,求求(x0 x1)(x0 x2)(xx1)(xx2)y0+(x1x0)(x1x2)(xx0)(xx2)y1+(x2x0)(x2x1)(xx0)(xx1)y2p2(7)=x0=1,x1=4,x2=9y0=1,y1=2,y2=3(14)(19)(74)(79)*1+(41)(49)(71)(79)*2+(91)(94)(71)(74)*3=2.7p2(x)=例例2.4已知函数已知函数y=f(x)在节点上满足在节点上满足xx0 x1x2yy0y1y2求二次多项式求二次多项式p(x)=a0+a1x+a2x2使之满
14、足使之满足p(xi)=yii=0,1,2解解:用待定系数法用待定系数法,将各节点值依次代入所求多项式将各节点值依次代入所求多项式,得得解上述方程解上述方程,将求出的将求出的a0,a1,a2代入代入p(x)=a0+a1x+a2x2即得所求二次多项式即得所求二次多项式例例2.5求过点求过点(0,1)、(1,2)、(2,3)的三点插值多项的三点插值多项式式解解:由由Lagrange插值公式插值公式（给定的三个点在一条直线上）（给定的三个点在一条直线上）例例2.6已知已知f(x)的观测数据的观测数据x0124f(x)19233 构造构造Lagrange插值多项式插值多项式解解四个点可构造三次四个点可构
15、造三次Lagrange插值多项式插值多项式:基函数为基函数为Lagrange插值多项式为插值多项式为为便于上机计算为便于上机计算,常将拉格朗日插值多项式常将拉格朗日插值多项式(5.8)改写成改写成例例2.7已知已知f(x)的观测数据的观测数据x1234f(x)0-5-63构造插值多项式构造插值多项式解解:四个点可以构造三次插值多项式四个点可以构造三次插值多项式,将数据将数据代入插值公式，有代入插值公式，有这个例子说明这个例子说明p(x)的项数不超过的项数不超过n+1项，但可以有项，但可以有缺项。缺项。拉拉格格朗朗日日插插值值算算法法实实现现x0 x1xixi+1xn-1xny=f(x)y=p(
16、x)ab在插值区间在插值区间 a,b 上用上用插值多项式插值多项式p(x)近似代替近似代替f(x),除了除了在插值节点在插值节点xi上没有误差外，在其它点上一般是存在误差上没有误差外，在其它点上一般是存在误差的。的。若记若记R(x)=f(x)-p(x)则则R(x)就是用就是用p(x)近似代替近似代替f(x)时的截断误差时的截断误差,或称或称插值余项我们可根据后面的定理来估计它的大小。插值余项我们可根据后面的定理来估计它的大小。插值多项式的误差插值多项式的误差定理定理2设设f(x)在在 a,b 有有n+1阶导数，阶导数，x0,x1,xn为为 a,b 上上n+1个互异的节点个互异的节点,p(x)为
17、满足为满足 p(xi)=f(xi)(i=1,2,n)的的n 次插值多项式，那么对于任何次插值多项式，那么对于任何x a,b 有有插值余项插值余项其中其中a b且依赖于且依赖于x证明证明(略略)对于线性插值，其误差为对于线性插值，其误差为对于抛物插值（二次插值），其误差为对于抛物插值（二次插值），其误差为例例2.8已知已知 =100,=121,用线性插值估计用线性插值估计 在在x=115时的时的截断误差截断误差解解:由插值余项公式知由插值余项公式知因为因为例例2.9已知已知x0=100,x1=121,x2=144,当用抛物插值求当用抛物插值求在在x=115时的近似值，估计其的截断误差时的近似值，
18、估计其的截断误差解解=例例2.10设设f(x)=x4,用余项定理写出节点用余项定理写出节点-1,0,1,2的三次插值多项式的三次插值多项式解解:根据余项定理根据余项定理3均差与均差与牛顿插值多项式牛顿插值多项式 拉格朗日插值多项式结构对称，使用方便。拉格朗日插值多项式结构对称，使用方便。但由于是用基函数构成的插值，这样要增加一个但由于是用基函数构成的插值，这样要增加一个节点时，所有的基函数必须全部重新计算，不具节点时，所有的基函数必须全部重新计算，不具备承袭性，还造成计算量的浪费。这就启发我们备承袭性，还造成计算量的浪费。这就启发我们去构造一种具有去构造一种具有承袭性承袭性的插值多项式来克服这
19、个的插值多项式来克服这个缺点，也就是说，每增加一个节点时，只需增加缺点，也就是说，每增加一个节点时，只需增加相应的一项即可。这就是牛顿插值多项式。相应的一项即可。这就是牛顿插值多项式。由线性代数知由线性代数知,任何一个不高于任何一个不高于n次的多项式次的多项式,都可以都可以表示成函数表示成函数的线性组合的线性组合,也就是说也就是说,可以把满足插值条件可以把满足插值条件p(xi)=yi(i=0,1,n)的的n次插值多项式次插值多项式,写成如下形式写成如下形式其中其中ak(k=0,1,2,n)为待定系数为待定系数,这种形式的插值多这种形式的插值多项式称为项式称为Newton插值多项式。我们把它记为
20、插值多项式。我们把它记为Nn(x)即即（3.12）可见，牛顿插值多项式可见，牛顿插值多项式Nn(x)是是插值多项式插值多项式p(x)的另的另一种表示形式一种表示形式,与与Lagrange多项式相比它不仅克服了多项式相比它不仅克服了“增加一个节点时整个计算工作重新开始增加一个节点时整个计算工作重新开始”的缺点的缺点,且可且可以节省乘除法运算次数以节省乘除法运算次数,同时在同时在Newton插值多项式中用插值多项式中用到差分与差商等概念，又与数值计算的其他方面有密切到差分与差商等概念，又与数值计算的其他方面有密切的关系的关系.它满足它满足其中其中ak(k=0,1,2,n)为待定系数，形如（为待定系
21、数，形如（3.12）的）的插值多项式称为插值多项式称为牛顿牛顿(Newton)插值多项式插值多项式。3.1差商及其性质差商及其性质定义定义函数函数y=f(x)在区间在区间xi,xi+1上的平均变化率上的平均变化率自变量之差和因变量之差之比叫自变量之差和因变量之差之比叫差商差商称为称为f(x)关于关于xi,xi+1的一阶差商的一阶差商,并记为并记为fxi,xi+1二阶差商二阶差商m阶差商阶差商fxi,xj,xk是指是指fxi,xj,xk=fxj,xk-fxi,xjxk-xi一般的一般的,可定义区间可定义区间xi,xi+1,xi+n上的上的n阶差商为阶差商为差商及其性质差商及其性质差商表差商表xi
22、fxifxi,xi+1fxi,xi+1,xi+2fxi,xi+1,xi+2x0f(x0)x1f(x1)fx0,x1x2f(x2)fx1,x2fx0,x1,x2x3f(x3)fx2,x3fx1,x2,x3fx0,x1,x2,x3fx1,x2-fx0,x1x2x0 xifxifxi,xi+1fxi,xi+1,xi+2fxi,xi+1,xi+2,xi+2002832751256216例例2.11求求f(xi)=x3在节点在节点x=0,2,3,5,6上的各阶差商值上的各阶差商值解解:计算得如下表计算得如下表在在n+1n+1个节点处各阶差商的计算方法个节点处各阶差商的计算方法差商及其性质差商及其性质这个
23、性质可用数学归纳法证明（用这个性质可用数学归纳法证明（用Lagrange插值多项式比较最高项系数来得到插值多项式比较最高项系数来得到）性质性质1函数函数f(x)的的n 阶差商阶差商fx0,x1,xn可可由由函数值函数值f(x0),f(x1),f(xn)的线性组的线性组合表示合表示,且且差商及其性质差商及其性质fx0,x1=fx1,x0f(x1)-f(x0)x1x0f(x0)-f(x1)x0 x1=性质性质2 2 差商具有对称性差商具有对称性,即在即在k k阶差商中阶差商中 任意交换两个节点任意交换两个节点 和和 的次序的次序,其值不变。其值不变。例如例如性质性质3若若fx,x0,x1,xk是是
24、x的的m 次多项式次多项式,则则fx,x0,x1,xk,xk+1是是x 的的m-1次多项式次多项式证：由差商定义证：由差商定义右端分子为右端分子为m 次多项式次多项式,且当且当x=xk+1时时,分子分子为为0,故分子含有因子故分子含有因子xk+1x，与分母相消后，右与分母相消后，右端为端为m-1次多项式。次多项式。4.4.1 差商及其性质差商及其性质性质性质4若若f(x)是是n次多项式次多项式,则则f x,x0,x1,xn恒为恒为0证：证：f(x)是是n次多项式，则次多项式，则f x,x0是是n-1次多次多项式项式,f x,x0,x1是是n-2次多项式次多项式,依次递推依次递推，f x,x0,
25、x1,xn-1是零次多项式，所以是零次多项式，所以fx,x0,x1,xn 0性质性质5 5 k k阶差商阶差商 和和k k阶导数之间有下阶导数之间有下 列关系列关系 这这个个性性质质可可直直接接用用罗罗尔尔（RolleRolle）定定理理证证明明（或或以下方法即余项方法）以下方法即余项方法）牛顿牛顿(Newton)插值多项式插值多项式的系数的系数可根据插值条件推出可根据插值条件推出,即由即由有有这是关于这是关于 的下三角方程组的下三角方程组,可以求得可以求得一般，用数学归纳法可证明一般，用数学归纳法可证明所以所以n n次牛顿次牛顿(Newton)Newton)插值公式为插值公式为 其余项其余项
26、为为牛牛顿顿插插值值多多项项式式的的误误差差。由由插插值值多多项项式式的的存存在在惟惟一一性性定定理理知知，满满足足同同一一组组插插值值条条件件的的拉拉格格朗朗日日插插值值多多项项式式P(x)P(x)与与牛牛顿顿插插值值多多项项式式N Nn n(x)(x)实实际际上上是是同同一一个个多多项项式式，仅仅是是同同一一插插值值多多项项式式的的不不同同表表达达形形式式而而已已，因因此此得得到到牛牛顿顿插插值值多多项项式式的的误误差差与与拉拉格朗日插值多项式的误差也完全相等。故有格朗日插值多项式的误差也完全相等。故有 可以看出，牛顿插值公式计算方便，增加可以看出，牛顿插值公式计算方便，增加一个插值点，只
27、要多计算一项，而一个插值点，只要多计算一项，而Nn(x)的的各项系数恰好是各阶差商值，很有规律各项系数恰好是各阶差商值，很有规律fx0,x(x-x0)=f(x)-f(x0)f(x)+fx0,x(x-x0)=f(x0)fx1,x0,x(x-x1)=fx0,x-fx1,x0fx0,x+fx1,x0,x(x-x1)=fx1,x0f(x)+(x-x0)fx1,x0=f(x0)+(x-x0)(x-x1)fx1,x0,x牛顿插值公式牛顿插值公式(另一种推导方法）另一种推导方法）f(x)=f(x0)+(x-x0)fx1,x0+(x-x0)(x-x1)fx1,x0,xfx1,x0,x=(x-x2)fx2,x1
28、,x0,x+fx2,x1,x0f(x)=f(x0)+(x-x0)fx1,x0+(x-x0)(x-x1)fx2,x1,x0 +(x-x0)(x-x1)(x-x2)fx2,x1,x0,xNn(x)Rn(x)如当如当n=1时，时，f(x)=f(x0)+(x-x0)fx1,x0+(x-x0)(x-x1)fx1,x0,xNn(x)=f(x0)+(x-x0)fx1,x0其中其中Nn(x)称为称为牛顿插值多项式牛顿插值多项式Rn(x)称为称为牛顿插值余项牛顿插值余项xifxifxi,xi+1fxi,xi+1,xi+2fxi,xi+1,xi+2x0f(x0)x1f(x1)fx0,x1x2f(x2)fx1,x2
29、fx0,x1,x2x3f(x3)fx2,x3fx1,x2,x3fx0,x1,x2,x34.4.2牛顿插值公式牛顿插值公式xifxifxi,xi+1fxi,xi+1,xi+2114293N2(7)=1+(7-1)*0.33333+(7-1)*(7-4)*(-0.01667)=2.69992+(x-x0)(x-x1)fx1,x0,x2+(x-x0)fx1,x0=f(x0)N(x)例例2.12已知已知x=1,4,9的平方根值，求的平方根值，求解：解：牛顿插值余项牛顿插值余项由由建起了差商和导数的关系建起了差商和导数的关系用导数代替牛顿插值多项式中的差商，有用导数代替牛顿插值多项式中的差商，有差商和导
30、数的关系也可用罗尔定理证出，余项差商和导数的关系也可用罗尔定理证出，余项R(x)=f(x)-P(x)R(xi)=f(xi)-P(xi)=0i=0,1,nRn(n)(x)=f(n)(x)-Pn(n)(x)=f(n)(x)-f(x0)+(x-x0)fx0,x1+(x-x0)(x-x1)fx0,x1,x2+(x-x0)(x-x1)(x-xn-1)fx0,x1,xn(n)=f(n)(x)-n!fx0,x1,xnRn(xi)=0(i=0,1,.,n)Rn(i)=0(i=0,1,.,n-1)Rn(n)()=0(x0,x1,xn)Rn(n)()=0=f(n)()-n!fx0,x1,xn即即R(x)R(x)在
31、在xx0 0,x,xn n 有有n+1n+1个零点，根据罗尔定理个零点，根据罗尔定理R R(n)(n)(x)(x)在在xx0 0,x,xn n 有有1 1个零点，设为个零点，设为，即有即有Rn(n)()=0增加新节点增加新节点x，并且并且f(x)为为(n+1)阶可导时，有阶可导时，有(x0,x1,xn)(x0,x1,xn,x)|f(x)(n+1)|Mn+14.4.1 差商及其性质差商及其性质例例2.13已知已知x=0,2,3,5对应的函数值为对应的函数值为y=1,3,2,5,作三次作三次Newton插值多项式。插值多项式。xif(xi)一阶差商一阶差商二阶差商二阶差商三阶差商三阶差商01231
32、32-1-2/3553/25/63/10所求的三次所求的三次Newton插值多项式为插值多项式为4.4.1 差商及其性质差商及其性质例例2.14已知已知f(x)=x7+x4+3x+1求求f20,21,27及及f20,21,27,28分析：本题分析：本题f(x)是一个多项式是一个多项式,故应利用差商的性质故应利用差商的性质解解:由差商与导数之间的关系由差商与导数之间的关系例例2.15求求并估计其误差并估计其误差解：作函数解：作函数f(x)=取取x0=4,x1=9,x2=6.25,建立差商表建立差商表xf(x)f xi,xi+1,fxi,xi+1,xi+242936.25 2.5N2(7)=2+(
33、7-4)*0.2+(7-4)*(7-9)*(-0.00808)=2.64848f 3(x)=Rn(x)在区间在区间4,9上，上，余式近似余式近似0.5*10-2,N2(7)=2.64848可舍入为可舍入为2.65|f(x)(n+1)|Mn+1由由差分与等距节点插值等距节点等距节点xi+1-xi=h，函数在等距节点上的值为函数在等距节点上的值为y0,y1,yn，称，称 yi-1=yi-yi-1为函数为函数f(x)在在xi-1,xi上的上的一阶差分一阶差分。称称 2yi-1=yi-yi-1=yi+1-2yi+yi-1为函数为函数f(x)在在xi-1,xi+1上的上的二阶差分二阶差分。称称 kyi-
34、1=k-1yi-k-1yi-1为函数为函数f(x)在在xi-1,xi+k-1上的上的k阶差分阶差分。当插值节点等距分布时当插值节点等距分布时,被插值函数的变化率就可用差被插值函数的变化率就可用差分来表示分来表示,这时牛顿插值公式的形式更简单这时牛顿插值公式的形式更简单,计算量更小计算量更小xy y 2y 3y 4yx0y0 x1y1x2y2x3y3x4y4 y0=y1y0 y1=y2y1 y2=y3y2 y3=y4y3 2y0=y1-y0 2y1=y2-y1 2y2=y3-y2 3y0=2y1-2y0 3y1=2y2-2y1 4y0等距节点插值等距节点插值 y0=y1y0 y1=y2y1 y2
35、=y3y2=y22y1+y0 2y0=y1-y0 3y0=2y1-2y0=y32y2+y1(y22y1+y0)=y33y2+3y1y0 2y1=y2-y1=y32y2+y1(a-b)3=a3-3a2b+3ab2-b3(a-b)2=a2-2ab+b2 4y0=3y1-3y0=y43y3+3y2y1-(y33y2+3y1y0)=y44y3+6y24y1+y0(a-b)4=a4-4a3b+6a2b2-4ab3+b3结论：各阶差分中函数值的系数正好等于结论：各阶差分中函数值的系数正好等于 (a-b)a-b)r r展开式中的系数展开式中的系数等距节点情况下等距节点情况下xi=x0+ih，用差分表示差商：
36、用差分表示差商：=y1y0h=y01!hfx1,x2=y2y1h=y11!hfx0,x1,x2=fx1,x2-fx0,x1x2x0=y11!h y01!h2h=y1-y02h2=2y02!h2fx1,x2,x3=fx3,x2-fx2,x1x3x1=y21!h y11!h2h=y2-y12!h2=2y12!h2fx0,x1,x2,x3=2y12!h2 2y02!h23h=2y1-2y02*3h3=3y03!h3 ny0n!hn例例2.16计算计算f(x)=x3在等距节点在等距节点0，1，2，3,4上的各上的各阶差分值阶差分值xy y 2y 3y001128327464 4y17193761218
37、660牛顿前插公式牛顿前插公式取间距为取间距为h,等距节点等距节点x0 x1xn顺序建立牛顿差商公式顺序建立牛顿差商公式fx0,x1=y01!hfx0,x1,x2=2y02!h2fx0,x1,x2,x3=3y03!h3Nn(x)=y0+(x-x0)y01!h+(x-x0)(x-x1)2y02!h2+(x-x0)(x-x1)(x-xn-1)ny0n!hn牛顿前插公式牛顿前插公式Nn(x)Rn(x)因因 ,设设 ,则则xy y 2y 3y 4yx0y0 x1y1 y0 x2y2 y1 2y0 x3y3 y2 2y1 3y0 x4y4 y3 2y2 3y1 4y0向后差分向后差分函数函数y=f(x)
38、,若记若记y-1=f(x0-h),y-2=f(x0-2h),则各阶向后差分则各阶向后差分一阶一阶 y0=y0-y-1，y1=y1-y0，y2=y2-y1，二阶二阶 2y0=y0-y-1=y0-y-1-(y-1-y-2)=y0-2y-1+y-2 2y1=y1-y0=y1-y0-(y0-y-1)=y1-2y0+y-1K阶阶 ky0=k-1y0-k-1y-1 ky1=k-1y1-k-1y0同样利用向后差分可以得到牛顿向后插值公式同样利用向后差分可以得到牛顿向后插值公式其中其中 ,公式公式 称之为牛顿向后插值公式余项。称之为牛顿向后插值公式余项。x-1012y-11311解：解：建立差分表建立差分表x
39、y y 2y 3y-1-10121320211 866=-1+1+0+0.375=0.375例例5.16按下列数值表用按下列数值表用牛顿前插公式求牛顿前插公式求y(-0.5)的近似值的近似值N3(x)例例5.17估计用线性插值法计算估计用线性插值法计算lg47时的误差限时的误差限取取x0=45,x1=48,=1.671898401解：应用解：应用n=1的拉格朗日插值公式的拉格朗日插值公式x424548lgx1.62324931.65321261.6812413(45,48)误差限误差限 插值公式的唯一性及其应用插值公式的唯一性及其应用n插值公式的唯一性插值公式的唯一性若插值节点相同，则插值公式
40、是唯一的。若插值节点相同，则插值公式是唯一的。Pn(x)与与Qn(x)有相同的插值节点，有相同的插值节点，令令Rn(x)=Pn(x)-Qn(x)对于对于x=x0,x1,xn,Rn(xi)=Pn(xi)-Qn(xi)=04 4 分段线性插值分段线性插值2.4.1 2.4.1 高次插值的龙格现象高次插值的龙格现象 插值多项式余项公式说明插值节点越多，一般说插值多项式余项公式说明插值节点越多，一般说来误差越小，函数逼近越好，但这也不是绝对的，来误差越小，函数逼近越好，但这也不是绝对的，因为余项的大小既与插值节点的个数有关，也与函因为余项的大小既与插值节点的个数有关，也与函数数f(x)的高阶导数有关的
41、高阶导数有关。换句话说，适当地提高插。换句话说，适当地提高插值多项式的次数，有可能提高计算结果的准确程度值多项式的次数，有可能提高计算结果的准确程度，但并非插值多项式的次数越高越好。当插值节点，但并非插值多项式的次数越高越好。当插值节点增多时，不能保证非节点处的插值精度得到改善，增多时，不能保证非节点处的插值精度得到改善，有时反而误差更大。考察函数有时反而误差更大。考察函数考察函数考察函数右图给出了右图给出了和和 的图像的图像,当当n增大时增大时,在两端在两端会发出激烈的振荡会发出激烈的振荡,这就是所谓龙格现这就是所谓龙格现象。该现象表明象。该现象表明,在在大范围内使用高次大范围内使用高次插值
42、插值,逼近的效果往逼近的效果往往是不理想的往是不理想的 另外，从舍入误差来看，高次插值误差的传播另外，从舍入误差来看，高次插值误差的传播也较为严重，在一个节点上产生的舍入误差会在计也较为严重，在一个节点上产生的舍入误差会在计算中不断扩大，并传播到其它节点上。因此，次数算中不断扩大，并传播到其它节点上。因此，次数太高的高次插值多项式并不实用，因为节点数增加太高的高次插值多项式并不实用，因为节点数增加时，计算量增大了，但插值函数的精度并未提高。时，计算量增大了，但插值函数的精度并未提高。为克服在区间上进行高次插值所造成的龙格现象，为克服在区间上进行高次插值所造成的龙格现象，采用分段插值的方法，将插
43、值区间分成若干个小的采用分段插值的方法，将插值区间分成若干个小的区间，在每个小区间进行线性插值，然后相互连接区间，在每个小区间进行线性插值，然后相互连接，用连接相邻节点的折线逼近被插函数，这种把插，用连接相邻节点的折线逼近被插函数，这种把插值区间分段的方法就是分段线性插值法。值区间分段的方法就是分段线性插值法。2.4.2 2.4.2 分段线性插值分段线性插值 分段线性插值就是通过插值节点用折线段连接起分段线性插值就是通过插值节点用折线段连接起来逼近来逼近f(x)。设设f(x)f(x)在在n+1n+1个节点个节点 上的函数值为上的函数值为 ,在每个小区间在每个小区间 (k=0,1,k=0,1,，
44、n n）上作线性插值，得上作线性插值，得 在几何上就是用折线在几何上就是用折线替代曲线替代曲线,如右图所示如右图所示若用插值基函数表示若用插值基函数表示,则在则在 a,b 上上其中其中显然，显然，是分段线性连续函数，且是分段线性连续函数，且 称称S S(x x)为为f f(x x)的的分段线性插值函数。分段线性插值函数。由线性插值的余项估计式知由线性插值的余项估计式知,f f(x x)在每个子段在每个子段上有误差估计式上有误差估计式其中其中 例例2.19 2.19 已知已知f(x)f(x)在四个节点上的函数值如下表所示在四个节点上的函数值如下表所示304560901求求f(x)f(x)在区间在
45、区间 30,9030,90 上的上的分段连续线性插值函数分段连续线性插值函数S(x)S(x)解解 将将插值区间插值区间 30,9030,90 分成连续的三个小区间分成连续的三个小区间 30,4530,45,45,6045,60,60,9060,90 则则S(x)在区间在区间 30,4530,45 上的上的线性插值为线性插值为S(x)在区间在区间 45,6045,60 上的上的线性插值为线性插值为S(x)S(x)在区间在区间 60,9060,90 上的线性插值为上的线性插值为 将各小区间的将各小区间的线性插值函数连接在一起，得线性插值函数连接在一起，得55三次样条三次样条插值插值 我们知道我们知
46、道,给定给定n+1+1个节点上的函数值可以作个节点上的函数值可以作n次次插插值多项式值多项式,但当但当n较大时较大时,高次高次插值不仅计算复杂插值不仅计算复杂,而且可而且可能出现能出现Runge现象现象,采用分段插值虽然计算简单、也有一采用分段插值虽然计算简单、也有一致收敛性致收敛性,但不能保证整条曲线在连接点处的但不能保证整条曲线在连接点处的光滑性光滑性,如分段线性插值如分段线性插值,其图形是其图形是锯齿形的折线锯齿形的折线,虽然连续虽然连续,但但处都是处都是“尖点尖点”,因而一阶导数都不存在因而一阶导数都不存在,这在实用上这在实用上,往往不能满足某些工程技术的高精度要求。如在往往不能满足某
47、些工程技术的高精度要求。如在船体、船体、飞机等外形曲线的设计中飞机等外形曲线的设计中,不仅要求曲线连续不仅要求曲线连续,而且要有而且要有二阶光滑度二阶光滑度,即有连续的二阶导数即有连续的二阶导数。这就要求分段插值。这就要求分段插值函数在整个区间上具有连续的二阶导数。因此有必要寻函数在整个区间上具有连续的二阶导数。因此有必要寻求一种新的插值方法求一种新的插值方法,这就是样条函数插值法这就是样条函数插值法 2.5.1 2.5.1 三次样条函数三次样条函数 定义定义5.4 5.4.设函数定义在区间设函数定义在区间 a a,b b 上，给定上，给定n+1n+1个个 节点和一组与之对应的函数值，若函数节
48、点和一组与之对应的函数值，若函数 满足满足:（1 1）在每个节点上满足）在每个节点上满足 S(S(xi i)=)=f(xi i)(i=0,1,)(i=0,1,n),n)（2 2）在）在 a a,b b 上有连续的二阶导数上有连续的二阶导数 （3 3）在每个小区间）在每个小区间 xi i,xi+1i+1(i=0,1,(i=0,1,n-1),n-1)上是一个三次多项式。上是一个三次多项式。则称则称S(S(x)为为三次样条插值函数三次样条插值函数。其中四个待定系数为其中四个待定系数为 ,子区间共有子区间共有n个个所以要确定所以要确定S(S(x)需要需要4n个待定系数。个待定系数。另一方面另一方面,要
49、求分段三次多项式要求分段三次多项式S(S(x)及其导数及其导数和和 在整个插值区间在整个插值区间 a,ba,b 上连续上连续,则要求它们在则要求它们在各个子区间的连接点各个子区间的连接点 上连续，上连续，即满足条件即满足条件 三次样条插值函数三次样条插值函数S(S(x)是一个分段三次多项式是一个分段三次多项式,要求出要求出S(S(x),),在每个小区间在每个小区间 xi i,xi+1i+1 上要确定上要确定4 4个待定参数个待定参数,若用若用S Si(x)表示它在第表示它在第i个子区间个子区间 xi i,xi+1i+1 上的表达式，则上的表达式，则（1 1）插值条件插值条件 （2 2）连接条件
50、连接条件 上述二式共给出了上述二式共给出了4 4n-2n-2个条件个条件,而待定系数有而待定系数有4 4n n个个,因此还因此还需要需要2 2个条件才能确定个条件才能确定S(x),S(x),通常在区间端点上通常在区间端点上 各加一个条件各加一个条件,称为称为边界条件边界条件,常用边常用边界条件有三种类型。界条件有三种类型。第一种类型：给定两端点第一种类型：给定两端点f(x)的一阶导数值：的一阶导数值：第二种类型：给定两端点第二种类型：给定两端点f(x)的二阶导数值：的二阶导数值：作为特例作为特例,称为称为自然边界条件自然边界条件。满。满足自然边界条件的三次样条插值函数称为足自然边界条件的三次样