第四章-插值与拟合ppt课件.ppt

《第四章-插值与拟合ppt课件.ppt》由会员分享，可在线阅读，更多相关《第四章-插值与拟合ppt课件.ppt（146页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物插值法n插值法是一种古老的数学方法，早在一千多年前的隋唐时期定制历法时就广泛应用了二次插值。刘焯将等距节点的二次插值应用于天文计算。n插值理论却是在17世纪微积分产生后才逐步发展起来的，Newton插值公式理论是当时的重要成果。n由于计算机的使用以及航空、造船、精密仪器的加工，插值法在理论和实践上都得到进一步发展，获得了广泛的应用。我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活

2、的生物我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物4.1 引言引言 问题的提出问题的提出函数解析式未知函数解析式未知,通过实验观测得到的一组数据通过实验观测得到的一组数据, 即在即在某个区间某个区间a, b上给出一系列点的函数值上给出一系列点的函数值 yi= f(xi)或者给出函数表或者给出函数表y=f(x)y=p(x)xx0 x1x2xnyy0y1y2yn我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物插值法的基本

3、原理插值法的基本原理设函数设函数 y=f(x) 定义在区间定义在区间 a, b 上上, , 是是 a, b 上取定的上取定的n+1个互异节点个互异节点, ,且在这些点处的函数值且在这些点处的函数值 为已知为已知 , ,即即 若存在一个若存在一个f(x)的近似函数的近似函数 , ,满足满足则称则称 为为f( (x) )的一个的一个插值函数插值函数, f( (x) )为为被插函数被插函数, 点点xi为为插值节点插值节点, 称称( (4.1)4.1)式为式为插值条件插值条件, 而误差函数而误差函数R(x)= 称为称为插值余项插值余项, 区间区间 a, b 称为称为插值插值区间区间, 插值点在插值区间

4、内的称为插值点在插值区间内的称为内插内插, 否则称否则称外插外插 nxxx,10)(,),(),(10nxfxfxf)(iixfy )(x), 2 , 1()()(nixfxii)(x( (4.1)4.1)()(xxf我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物插值函数插值函数 在在n+1 个互异插值节点个互异插值节点 ( (i=0,1,n )处与处与 相等相等, ,在其它点在其它点x就用就用 的值作为的值作为f( (x) ) 的近似值。这一过程称为的近似值。这一过程称为插值插值，点，点x称为插值点。称

5、为插值点。换句话说换句话说, , 插值就是根据被插函数给出的函数表插值就是根据被插函数给出的函数表“插出插出”所要点的函数值。用所要点的函数值。用 的值作为的值作为f( (x) )的近似值的近似值, ,不仅不仅希希望望 能较好地逼近能较好地逼近 f( (x) ), ,而且还希望它计算简单而且还希望它计算简单 。由于代数多项式具有数值计算和理论分析方便的优点。由于代数多项式具有数值计算和理论分析方便的优点。所以本章主要介绍代数插值。即求一个次数不超过所以本章主要介绍代数插值。即求一个次数不超过n n次的次的多项式。多项式。 )(xix)(ixf)(x)(x)(x0111)(axaxaxaxPnn

6、nn我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物0111)(axaxaxaxPnnnn满足满足 ), 2 , 1 , 0()()(nixfxPii则称则称P(x)P(x)为为f(x)f(x)的的n n次插值多项式。这种插值法通常称次插值多项式。这种插值法通常称为代数插值法。其几何意义如下图所示为代数插值法。其几何意义如下图所示 y y=P(x) y=f(x) y1 yn x0 x1 xn x 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有

7、错：表里边有一个活的生物定理定理1 n次代数插值问题的解是存在且惟一的次代数插值问题的解是存在且惟一的 证明证明: : 设设n n次多项式次多项式 0111)(axaxaxaxPnnnn是函数是函数 在区间在区间 a, ba, b上的上的n+1n+1个互异的节点个互异的节点 ( (i=0,1,2,n )i=0,1,2,n )上的插值多项式上的插值多项式, ,则求插值多项式则求插值多项式P(x)P(x)的问题就归结为求它的系数的问题就归结为求它的系数 ( (i=0,1,2,n )i=0,1,2,n )。 )(xfy ixia由插值条件由插值条件: (: (i=0,1,2,n),i=0,1,2,n

8、),可得可得 )()(iixfxp)()()(01111011111100011010nnnnnnnnnnnnnnnnxfaxaxaxaxfaxaxaxaxfaxaxaxa我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物 这是一个关于待定参数这是一个关于待定参数 的的n+1阶线性方阶线性方程组程组, ,其系数矩阵行列式为其系数矩阵行列式为 naaa,10niijjinnnnnnxxxxxxxxxxxV110212110200)(111 称为称为Vandermonde（范德蒙）行列式，因范德蒙）行列式，因xi

9、xj（当当ij），），故故V0。根据解线性方程组的克莱姆根据解线性方程组的克莱姆（Gramer）法则，方程组的解法则，方程组的解 存在惟一，从而存在惟一，从而P(x)P(x)被惟一确定。被惟一确定。 naaa,10惟一性说明，不论用何种方法来构造，也不论用何种惟一性说明，不论用何种方法来构造，也不论用何种形式来表示插值多项式形式来表示插值多项式, ,只要满足插值条件只要满足插值条件( (4.1)4.1)其结其结果都是相互恒等的。果都是相互恒等的。 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物4.2 拉格

10、朗日（拉格朗日（Lagrange）插值插值 为了构造满足插值条件为了构造满足插值条件 (i=0,1,2,n )的便于使用的插值多项式的便于使用的插值多项式 L( (x),),先考察几种简单情形先考察几种简单情形, ,然后再推广到一般形式。（然后再推广到一般形式。（ 线性插值与抛物插值）线性插值与抛物插值）（1）线性插值）线性插值线性插值是代数插值的最简单形式。假设给定了函数线性插值是代数插值的最简单形式。假设给定了函数f(x)在两个互异的点的值，在两个互异的点的值，, ,现要求用线性函数现要求用线性函数 近似地代替近似地代替 f(x)。选选择参数择参数a和和b, 使使 。称这样的线性函数。称这

11、样的线性函数L1(x) 为为 f(x) 的线性插值函数的线性插值函数 。( )( )iiL xf x0 x1x)(),(1100 xfyxfy1( )L xax b1( )( )(0,1)iiL xf xi我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物线性插值的几何意义线性插值的几何意义: :用用通过点通过点 和和 的直线近似地代替曲线的直线近似地代替曲线 y=f(x)=f(x)由解析几何知道由解析几何知道, ,这条直线用点斜式表示为这条直线用点斜式表示为 )(,(00 xfxA)(,(11xfxB101

12、0010( )()yyL xyxxxx011010110( )xxxxL xyyxxxx01011010)(,)(xxxxxlxxxxxl0)(, 1)(1000 xlxl1)(,0)(1101xlxl1)()(10 xlxl为了便于推广，记为了便于推广，记 这是一次函这是一次函数数, ,且有性质且有性质 y=f(x) L1(x)=ax+b A(x.0,f(x.0) B(x.1,f(x.1) 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物)(0)(1)(kikixlkiik 与与 称为线性插值基函数。且有

13、称为线性插值基函数。且有 )(0 xl)(1xl1 , 0,)(10kxxxxxlkjjjkjk于是线性插值函数可以表示为与基函数的线性组合于是线性插值函数可以表示为与基函数的线性组合 10011( )( )( )L xl x yl x y例例4.14.1 已知已知 , , , , 求求 10100 11121 115y解解: : 这里这里x0=100，y0=10，x1=121，y1=11, 利用线性插值利用线性插值 1121100( )1011100121121 100 xxL x1115(115)10.714yL我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？

14、但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物（2 2）抛物插值）抛物插值 抛物插值又称二次插值，它也是常用的代数插值之一。抛物插值又称二次插值，它也是常用的代数插值之一。设已知设已知f(x)f(x)在三个互异点在三个互异点x0，x1，x2的函数值的函数值y0，y1，y2，要构造次数不超过二次的多项式要构造次数不超过二次的多项式使满足二次插值条件：使满足二次插值条件：这就是二次插值问题。其几何意义是用经过这就是二次插值问题。其几何意义是用经过3个点个点 的抛物线的抛物线 近似代替曲线近似代替曲线 , ,如下图所示。因此也称之为抛物插值。如下图所示。因此也称之为抛物插值。 22210

15、( )L xa xa xa2( )(0,1,2)iiL xyi),(),(),(221100yxyxyx2( )yL x)(xfy 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物 y y=L2(x) y0 y1 y1 y=f(x) O x0 x1 x2 x L2(x)的参数的参数直接由插值条件决定，直接由插值条件决定，即即 满足下面满足下面的代数方程组：的代数方程组： 210,aaa210,aaa222221012121100202010yxaxaayxaxaayxaxaa222211200111xxxx

16、xx该三元一该三元一次方程组次方程组的系数矩阵的系数矩阵 的行列式是范德蒙行列式，当的行列式是范德蒙行列式，当 时，时，方程组的解唯一。方程组的解唯一。 210 xxx我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物为了与下一节的为了与下一节的Lagrange插值公式比较插值公式比较, ,仿线性插值仿线性插值, ,用用基函数的方法求解方程组。先考察一个特殊的二次插值基函数的方法求解方程组。先考察一个特殊的二次插值问题：问题： 求二次式求二次式 , ,使其满足条件：使其满足条件： )(0 xl0)(,0)(,

17、 1)(201000 xlxlxl这个问题容易求解。由上式的后两个条件知这个问题容易求解。由上式的后两个条件知: : 是是 的两个零点。于是的两个零点。于是 21,xx)(0 xl)()(210 xxxxcxl再由另一条件再由另一条件 确定系数确定系数 1)(00 xl)(12010 xxxxc)()()(2010210 xxxxxxxxxl从而导出从而导出 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物类似地可以构造出满足条件：类似地可以构造出满足条件：的插值多项式的插值多项式 0)(, 0)(, 1)

18、(210111xlxlxl)()()(2101201xxxxxxxxxl及满足条件：及满足条件： 的插值多项式的插值多项式 0)(,0)(, 1)(120222xlxlxl)()()(1202102xxxxxxxxxl这样构造出来的这样构造出来的 称为抛物插值的基函数称为抛物插值的基函数 )(),(),(210 xlxlxl取已知数据取已知数据 作为线性组合系数作为线性组合系数, ,将基函数将基函数 线性组合可得线性组合可得 210,yyy)(),(),(210 xlxlxl0201122012010210122021()()()()()()( )()()()()()()x xx xx xx

19、xx xx xL xyyyxxxxxxxxxxxx容易看出容易看出, ,P(x)P(x)满足条件满足条件 2()(0,1,2)iiLxyi我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物4.2.1 4.2.1 拉格朗日插值多项式拉格朗日插值多项式 两个插值点可求出一次插值多项式两个插值点可求出一次插值多项式, ,而三个插值点可而三个插值点可求出二次插值多项式。插值点增加到求出二次插值多项式。插值点增加到n+1个时个时, ,也就是通也就是通过过n+1个不同的已知点个不同的已知点 , ,来构造一个次数来构造一个

20、次数为为n的代数多项式的代数多项式Ln(x)。与推导抛物插值的基函数类似与推导抛物插值的基函数类似, ,先构造一个特殊先构造一个特殊n次多项式次多项式 的插值问题的插值问题, ,使其在各节点使其在各节点 上满足上满足 ), 1 , 0)(,(niyxii)(xliix0)(, 0)(, 1)(, 0)(, 0)(110nkkkkkkkkxlxlxlxlxl)(0)(1)(kikixlkiik即即 由条件由条件 ( ) ( )知知, , 都是都是n n次次 的零点的零点, ,故可设故可设 0)(ikxlki nkkxxxxx,1110)(xlk我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放

21、在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物)()()()(110niiixxxxxxxxAxl 设设.为为待待定定常常数数其其中中 A可可得得由由1)( iixl)()()(1110niiiiiixxxxxxxxA 0110110()()()()( )()()()()()(0,1,)(*)()iiniiiiiiinnjjijjixxxxxxxxlxxxxxxxxxxxinxx所以称之为称之为Lagrange插值基函数插值基函数.我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一

22、个活的生物利用拉格朗日基函数利用拉格朗日基函数,可以构造多项式可以构造多项式0 01 10( )( )( )( )( ) (4.2)nnn ni iiLxy lxy lxy lxy lx.Lagrange)(.,)(,)(插插值值多多项项式式为为称称解解故故其其为为拉拉格格朗朗日日问问题题的的且且满满足足插插值值条条件件的的多多项项式式为为次次数数不不超超过过xLyxLnxLniinn插值多项式为：插值多项式为：线性插值多项式：线性插值多项式：n=1010110101xxxxyxxxxyxL)(),(,)(101iyxLii满满足足y=f (x)xyx0 x1y=L1(x).),(),()(的

23、的直直线线为为过过点点11001yxyxxLy 几何意义：几何意义：我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物抛物插值多项式：抛物插值多项式：n=2插值多项式为：插值多项式为：)()()()()()()(1202102210120120102102xxxxxxxxyxxxxxxxxyxxxxxxxxyxL),(,)(2102iyxLii满满足足xyy=L2 (x)x0y=f(x)x1x2.),(),(),()(的的一一条条抛抛物物线线和和为为过过点点2211002yxyxyxxLy 几何意义：几何意义

24、：101 ( )()()() (4.3)nnxxxxxxx引入记号)()()()( 1101nkkkkkkknxxxxxxxxx 则则得得101( ) ( ) (4.4)()()nnnkkknkxLxyxxx于是我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.2处处的的近近似似值值。在在公公式式求求，利利用用插插值值，的的值值分分别别为为：，在在设设2007408180778801086070809048370300250150100. . . . ,. ,. .)(xxexexf解解： )()()

25、()()()()()()()()()()()()()()(23130321033212023102312101320130201032103xxxxxxxxxxxxyxxxxxxxxxxxxyxxxxxxxxxxxxyxxxxxxxxxxxxyxL代代入入分分别别将将 .,.,.,.,. 3002501501002003210 xxxxx81873002003.).( L可可得得. . .52001081873080相相比比，误误差差与与准准确确结结果果e我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物定理

26、定理4.14.1nnnxaxaxaaxP2210)(niyxPiin,2 ,1 ,0)(),(jixxji若插值节点则满足插值条件则满足插值条件的插值多项式的插值多项式存在且唯一存在且唯一.反证：若不唯一，则除了反证：若不唯一，则除了Pn(x) 外还有另一外还有另一 n 阶多项式阶多项式 Ln(x) 满足满足 Ln(xi) = yi 。考察考察 则则 Qn 的阶数的阶数, )()()(xLxPxQnnn n而而 Qn 有有 个不同的根个不同的根n + 1x0 xn注：注：若不将多项式次数限制若不将多项式次数限制为为 n ，则插值多项式，则插值多项式不唯一不唯一。例如例如 也是一个插值多项也是一

27、个插值多项式，其中式，其中 可以是任意多项式。可以是任意多项式。 niinxxxpxLxP0)()()()()(xp我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.3 已知已知y=f(x)的函数表的函数表 求线性插值多项式求线性插值多项式, 并计算并计算x=1.5 的值的值X 1 3 y 1 225.1)5.1()5.1()1(2121311313)(10100101pfxxxyxxxxyxxxxxp解解: 由线性插值多项式公式得由线性插值多项式公式得我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把

28、它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.4 已知已知x=1, 4, 9 的平方根值的平方根值, 用抛物插值公式用抛物插值公式, 求求 (x0 x1)(x0 x2)(xx1)(xx2)y0+(x1x0)(x1x2)(xx0)(xx2)y1+(x2x0)(x2x1)(xx0)(xx1)y2L2(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.7L2(x) =7解解:我吓了一跳，蝎子是多么丑恶

29、和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.5 已知函数已知函数y=f(x)在节点上满足在节点上满足 x x0 x1 x2 y y0 y1 y2 求二次多项式求二次多项式 p(x) = a0 + a1x + a2x2 使之满足使之满足 p(xi) = yi i=0, 1, 2解解: 用待定系数法用待定系数法, 将各节点值依次代入所求多项式将各节点值依次代入所求多项式, 得得解上述方程解上述方程, 将求出的将求出的a0, a1, a2 代入代入p(x) = a0 + a1x + a2x2 即得所求二次多项式即得所求二次多项

30、式 201202012120122aa x a xyaa x a xyaa x a xy我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.6 求过点求过点(0,1)、(1,2)、(2,3)的三点插值多项式的三点插值多项式13) 12)(02 () 1)(0(2) 21)(01 () 2)(0(1) 20)(10 () 2)(1()(xxxxxxxxp解解:由由Lagrange 插值公式插值公式（给定的三个点在一条直线上）（给定的三个点在一条直线上）212021012101200201021)()()(

31、)()()()(yxxxxxxxxyxxxxxxxxyxxxxxxxxxP我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.7 已知已知f (x)的观测数据的观测数据 x 0 1 2 4 f (x) 1 9 23 3 构造构造Lagrange插值多项式插值多项式解解 四个点可构造三次四个点可构造三次Lagrange插值多项式插值多项式: :基函数为基函数为 1478781)40)(20)(10()4)(2)(1()(230 xxxxxxxlxxxxxxxl38231)41)(21)(01 ()4)(2

32、)(0()(231xxxxxxxl2324541)42)(12)(02()4)(1)(0()(我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物xxxxxxxl12181241)24)(14)(04()2)(1)(0()(233Lagrange插值多项式为插值多项式为 )()(303xlyxLkkk)(3)(23)(9)(3210 xlxlxlxl12144541123xxx为便于上机计算为便于上机计算, ,常将拉格朗日插值多项式常将拉格朗日插值多项式( (4.2)改写成改写成 nknkiiikiknxxxx

33、yxL00)(我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物30011223332( )( )( )( )( )43L xl x yl x yl x yl x yxx 例例4.8 已知已知f(x)的观测数据的观测数据 x 1 2 3 4f(x) 0 -5 -6 3构造插值多项式构造插值多项式 解解: 四个点可以构造三次插值多项式四个点可以构造三次插值多项式, 将数据将数据 代入插值公式，有代入插值公式，有 这个例子说明这个例子说明Ln(x)的项数不超过的项数不超过n+1项，但可以有缺项。项，但可以有缺项。

34、我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物Lagrange插值法的流程图插值法的流程图我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物x0 x1 xixi+1 xn-1 xny=f(x)y=Ln(x)ab在插值区间在插值区间 a, b 上用上用插值多项式插值多项式Ln(x)近似代替近似代替f(x), 除了除了在插值节点在插值节点xi上没有误差外，在其它点上一般是存在误上没有误差外，在其它点上一般是存在误差的。差的。

35、若记若记 R (x) = f(x) - Ln(x) 则则 R(x) 就是用就是用 Ln(x) 近似代替近似代替 f(x) 时的截断误差时的截断误差, 或称或称插值余项我们可根据后面的定理来估计它的大小。插值余项我们可根据后面的定理来估计它的大小。4.2.2 插值多项式的误差插值多项式的误差 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物定理定理4.2 设设f(x)在在 a, b 有有n+1阶导数，阶导数， x0, x1, xn 为为 a, b 上上n+1个互异的节点个互异的节点, Ln(x)为满足为满足

36、Ln(xi) = f(xi) (i=1,2, , n) 的的n 次插值多项式，那么对于任何次插值多项式，那么对于任何x a, b 有有 插值余项插值余项其中其中a b 且依赖于且依赖于x1010( )()()()(),nnniixx xx xx xx xa b(1)1( )( )( )( )( )(1)!nnnnfR xf xL xxn我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物 插值余项插值余项 /* Remainder */设节点设节点)1( nf在在a , b内存在内存在, 考察截断误差考察截断误

37、差)()()(xLxfxRnn ,baCfn bxxxan 10，且，且 f 满足条件满足条件 , Rolles Theorem: 若若 充分光滑，充分光滑， ，则，则存在存在 使得使得 。)(x 0)()(10 xx ),(10 xx 0)( 推广：推广：若若0)()()(210 xxx ),(),(211100 xxxx 使得使得0)()(10 ),(10 使得使得0)( 0)()(0 nxx 存在存在),(ba 使得使得0)()( nRn(x) 至少有至少有 个根个根n+1 niinxxxKxR0)()()(任意固定任意固定 x xi (i = 0, , n), 考察考察 niixtxK

38、tRnt0)()()()( (t)有有 n+2 个不同的根个不同的根 x0 xn x),(,0)()1(baxxn !)1()()()1(nxKRxnn 注意这里是对注意这里是对 t 求导求导 !)1)()()()1()1(nxKLfxnnxn !)1()()()1( nfxKxn niixnnxxnfxR0)1()(!)1()()( 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物注：注： 若记若记 ，则插值余项为，则插值余项为 当当 f(x) 为任一个次数为任一个次数 n 的的多项式多项式时，时， ,

39、可可知知 ，即插值多项式对于次数，即插值多项式对于次数 n 的的多项式是多项式是精精确确的。的。(若是代数插值若是代数插值,其插值函数就是其插值函数就是f(x) )0)()1( xfn0)( xRn定义定义：.max,min00为为基基本本插插值值区区间间围围称称插插值值节节点点所所界界定定的的范范iniinixx 当点当点x位于基本插值区间内时位于基本插值区间内时,插值过程称为插值过程称为内插内插,否则称为否则称为外推外推.1)1()( nnMxf)()!(xnMnn111 通常不能确定通常不能确定 , 而是估计而是估计 , x (a,b) 于是将于是将 作为误差估计上限。作为误差估计上限。

40、 niinxxx01)()( . )()!()()()(xnfxRnnn111 外推比内插效果差。外推比内插效果差。 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物对于线性插值，其误差为对于线性插值，其误差为对于抛物插值（二次插值），其误差为对于抛物插值（二次插值），其误差为011( )( )( )( )()(),2R xf xP xfx xx xa bbaxxxxxxfxPxfxR,)( )()(61)()()(210 (1)11n 1max |( )|, |( )|( )|, (1)!nna x bn

41、nfxMMR xxn 若则我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.9 已知已知 =100, =121, 用线性插值估计用线性插值估计 在在x=115时的时的截断误差截断误差xxf)(0 x1x解解: 由插值余项公式知由插值余项公式知 )()(21)(1xfxR 2341)( xxf)(81)(10231xxxxxR因为因为 )121115)(100115(81)115(231R23121,100max)121115)(100115(81)121115)(100115(1081301125.

42、010615813我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.10 已知已知x0=100, x1=121, x2=144,当用抛物插值求当用抛物插值求 在在x=115时的近似值，估计其的截断误差时的近似值，估计其的截断误差 ( )f xx解解0201122012010210122021()()()()()()( )()()()()()()115(115) 10.772756x xx xx xx xx xx xp xyyyxxxxxxxxxxxxp=0017. 010)144115)(121115

43、)(100115(161)115()144)(121)(100(161)()()()(61)(52252210) 3(2RxxxxxRxxxxxxfxR2583)( xxf我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.11 设设f(x)=x4, 用余项定理写出节点用余项定理写出节点 -1, 0, 1, 2的三次插值多项式的三次插值多项式 解解: 根据余项定理根据余项定理(4)0123432( )( )( )()()()()4!( )( 1 )( 1 )( 2)( ) 22ff x pxx x x

44、x x x x xx px xxxxpxx xx 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例4.12 已知 sin0.32=0.314567，sin0.34=0.333487，sin0.36=0.352274，分别用一、二次Lagrange插值计算 sin0.3367的值，并估计截断误差。（1）解3403203403203203403203401.sin.sin.)( xxxL3403203403203367032034032034033670336701.sin.sin.).( L3303650.

45、 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物得).)(.(sin)(34032021 xxxR 510920 .)()!()()()()()(xnfxLxfxRnnnn111 由于是|.|.|sin| ).(|34033670320336702336701 R3400000918920003300167023334870. 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物（2）3603403603403403603

46、403601.sin.sin.)( xxxL3603403603403367034036034036033670336701.sin.sin.).( L3303870. 得).)(.(sin)(36034021 xxxR )()!()()()()()(xnfxLxfxRnnnn111 由我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物于是|.|.|sin| ).(|36033670340336702336701 R700001354310023300033023522740. 700001354310023

47、300033023522740. 510351 .由此可知 稍好于 ).(336701L).(336701L我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物（3）3603403603203603403203403603403203403603203203603203403203603402.sin).)(.().)(.(.sin).)(.().)(.(.sin).)(.().)(.()( xxxxxxxL3303740336702.).( L因为).)(.)(.(!)()()(360340320332 xx

48、xfxR 我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物)(xL1)(xL1则|.|.|.|!|cos| ).(|3603403203336702 xxxR 0233000330016703320.!.cos 7101780 .oxyxsin)(xL2我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物例例4.13：已知：已知233sin,214sin,216sin 分别利用分别利用 sin x 的的1次、次、2次次 La

49、grange 插值计算插值计算 sin 50 并估计误差。并估计误差。 解：解：0 x1x2x185500 n = 1分别利用分别利用x0, x1 以及以及 x1, x2 计算计算4,610 xx利用利用21646214641/)( xxxL这里这里),(,sin)(,sin)()(462 fxxf而而)(!)()(,sin)(462222121 xxfxR00762. 0)185(01319. 01 Rsin 50 = 0.7660444)185(50sin10 L0.77614外推外推 的实际误差的实际误差 0.01010.01013,421 xx利用利用00660.018500538.0

50、1 R内插内插 的实际误差的实际误差 0.005960.00596内插通常优于外推。选择要内插通常优于外推。选择要计算的计算的 x 所在的区间的端点，所在的区间的端点，插值效果较好。插值效果较好。sin 50 = 0.76008, )(1851 L我吓了一跳，蝎子是多么丑恶和恐怖的东西，为什么把它放在这样一个美丽的世界里呢？但是我也感到愉快，证实我的猜测没有错：表里边有一个活的生物n = 22321214363463464363646342)()()()()()()( xxxxxxxL)185(50sin20 L0.76543232134632 cos;)()(!cos)(xxxxR00077