计算声学第五章 插值法

1、第五章第五章 插值法插值法 在实际科学计算中常会出现这样的情况，由于函数的解析表达式过于复杂不便计算，但是需要计算多个点处的函数值；或者函数的解析表达式 未知，仅知道它在区间 内n+1个互异点 处对应的函数值 ，需要构造一个简单函数 作为函数 的近似表达式 ，使得这类问题称为插值问题插值问题， 称为插值函数插值函数。如果插值函数类 是代数多项式，则相应的插值问题称为代数插值代数插值；如果 是三角多项式，则相应插值问题称为三角插值三角插值。)(xf,banxxx,10nyyy,10)(xP)(xfy )()(xPxfy), 1 , 0( )()(niyxfxPiii)(xP)(xP)(xP1 拉

2、格朗日（Lagrange）插值l 代数插值代数插值定义定义：设 在区间 上有定义，且在 上的n+1个不同点 的函数值为 ，如果存在一个代数多项式 ，其中 为实数，使得成立，则称 为函数 的插值多项式插值多项式，点称为插值节点插值节点，包含插值节点的区间 称为插值区间插值区间，关系式 称为插值条件插值条件。求插值多项式的问题称为代数插值问题代数插值问题。)(xfy ,ba,babxxxan10nyyy,10nnnxaxaxaaxP2210)(ia), 2 , 1 , 0( )(niyxPiin)(xPn)(xfy nxxx,10,ba), 2 , 1 , 0( )(niyxPiin)(xPn1

3、拉格朗日（Lagrange）插值l 几何意义几何意义：1 拉格朗日（Lagrange）插值 几何意义为通过n+1个点 做一条代数曲线 ，使其近似于曲线 ，利用 上的点近似代替 上的点。（用于对函数的离散数据建立简单的数学模型）余项余项：在区间 上用 近似 ，除了在插值节点处 外，在区间其余点处一般都有误差。令 ，则 称为插值多项式的余项，它表示用 近似 时产生的截断误差截断误差。一般情况下，越小，近似程度越好。 ), 2 , 1 , 0( ),(niyxii)(xPyn)(xfy )(xPyn)(xfy ,ba)(xPn)(xfix)()(iinxfxP)()()(xPxfxRnn)(xRn)

4、(xPn)(xf)(maxxRnbxa1 拉格朗日（Lagrange）插值定理定理：在n+1个互异节点 上满足插值条件的次数不高于n次的插值多项式 存在且唯一。证明证明：如果插值多项式 的系数 可以被唯一确定，则该多项式存在并且唯一。由插值条件 ，插值多项式中的系数 满足n+1阶线性方程组 ix), 2 , 1 , 0( )(niyxPiin)(xPnnnnxaxaxaaxP2210)(naaa,10), 2 , 1 , 0( )(niyxPiinnaaa,10nnnnnnnnnnyxaxaxaayxaxaxaayxaxaxaa221011212110002020101 拉格朗日（Lagran

5、ge）插值方程组中未知量 的系数行列式为范德蒙行列式因为插值点互不相同，即 ，所以 ，方程组有唯一解 ，即插值多项式存在并且唯一。naaa,10nijjinnnnnnxxxxxxxxxxxV0212110200)(111)(jixxji0Vnaaa,101 拉格朗日（Lagrange）插值l 线性插值线性插值1 拉格朗日（Lagrange）插值l 定义定义：设函数 在区间 端点的值分别为 ，用线性函数 来近似代替 ，确定参数 ，使则称线性函数 为 的线性插值函数线性插值函数。l 几何意义几何意义：如图所示，利用通过两点 和 的直线 去近似代替曲线 。)(xfy ,10 xx)(),(1100

6、xfyxfybaxxLy)(1)(xfba,)()( ),()(111001xfxLxfxL)(1xL)(xf)(,(00 xfxA)(,(11xfxB)(1xLy )(xfy 1 拉格朗日（Lagrange）插值由直线方程的两点式方程可求得 的表达式为：记则 都为 的一次函数，并且具有下列性质我们把具有这种性质的函数 称为线性插值基函数线性插值基函数。)(1xL101001011)(yxxxxyxxxxxL01011010)( ,)(xxxxxlxxxxxl)(),(10 xlxlx) 1 , 0,( , 0 , 1)(kikikixlik)(xlk1 拉格朗日（Lagrange）插值线性插

7、值函数 用基函数可以表示为上式说明，任何一个满足插值条件的线性插值函数都可由线性插值基函数 的一个线性组合来表示。)(1xL)()()(11001xlyxlyxL)()( ),()(111001xfxLxfxL)(),(10 xlxl1 拉格朗日（Lagrange）插值定理定理：设 在区间 上连续， 在 内存在，是满足插值条件 的插值多项式，则对任何 ，插值余项（截断误差）为其中 ，且依赖于 。如果 ，则截断误差限截断误差限是)(xf ,10 xx)(xf ),(10 xx)(1xL)()( ),()(111001xfxLxfxL,10 xxx)(! 2)()()()(1011xxxxfxLx

8、fxR ),(10 xxx1)(max10Mxfxxx )(! 2)(1011xxxxMxR1 拉格朗日（Lagrange）插值l 抛物线插值抛物线插值1 拉格朗日（Lagrange）插值 设已知 在三个不同点 上的值分别为 ，做一个二次插值多项式 ，使其满足插值条件由于通过不在同一直线上的三点可做一条抛物线，所以称二次插值多项式 为 的抛抛物线插值函数物线插值函数。 设二次插值多项式为（插值基函数的线性组合）)(xfy 210,xxx210,yyy)(2xL)2 , 1 , 0( )(2iyxLii)(,(),(,(),(,(221100 xfxCxfxBxfxA)(2xL)(xf20221

9、1002 ),()()()(xxxxlyxlyxlyxL1 拉格朗日（Lagrange）插值其中 都是二次多项式，且满足已知 ，即 是 的两个零点，所以设其中 为待定常数。由 得到所以)2 , 1 , 0( )(kxlk)2 , 1 , 0,( , 0 , 1)(kikikixlik0)()(2010 xlxl21,xx)(0 xl)()(210 xxxxkxlk1)(00 xl)(11)(20102010 xxxxkxxxxk)()()(2010210 xxxxxxxxxl1 拉格朗日（Lagrange）插值同样求得所以上式又称为 的二次拉格朗日插值多项式二次拉格朗日插值多项式。)()()(

10、2101201xxxxxxxxxl)()()(1202102xxxxxxxxxl)()()()()()()(1202102210120120102102xxxxxxxxyxxxxxxxxyxxxxxxxxyxL)(20 xxx)(xf1 拉格朗日（Lagrange）插值l 截断误差截断误差与截断误差限截断误差限 如果 在区间 上连续， 在 内存在，则用 去近似 的截断误差截断误差为其中 ，并且依赖于 。如果 ，则截断误差限截断误差限为)(xf ,20 xx)(xf ),(20 xx)(2xL)(xf)()(! 3)()()()(21022xxxxxxfxLxfxR ),(20 xxx2)(ma

11、x20Mxfxxx )()(! 3)(21022xxxxxxMxR1 拉格朗日（Lagrange）插值例例：设 ，试分别应用线性插值和抛物线插值公式计算 的近似值（13.2287565553）。解解：取 ，则对应的以 为节点做线性插值以 为节点做抛物线插值xxf)()175(f196,169,144210 xxx14,13,12210yyy21,xx222222.131416919616917513196169196175)175(1L210,xxx229402.13 14)169196)(144196()169175)(144175( 13)196169)(144169()196175)(1

12、44175(12)196144)(169144()196175)(169175()175(2L1 拉格朗日（Lagrange）插值由 得到 ，所以 xxf)(2/52/383)(,41)( xxfxxf32/31101138. 016941M52/52101507. 014483M00717. 0)196175)(169175(! 2)175(11MR00098. 0)196175)(169175)(144175(! 3)175(22MR1 拉格朗日（Lagrange）插值l 拉格朗日插值多项式拉格朗日插值多项式 设函数 在节点 处的函数值为 ，做一个n次插值多项式 ，并使 在节点处满足则n次

13、插值基函数 ，就是在n+1个节点 上满足条件的n次多项式。)(xfy nxxx10), 1 , 0( )(nkxfykk)(xLn)(xLn), 2 , 1 , 0( )(nkyxLkkn), 2 , 1 , 0( )(nkxlknxxx10), 2 , 1 , 0,( , 0 , 1)(nkikikixlik1 拉格朗日（Lagrange）插值经过推导得出n次插值基函数显然 满足插值条件，所以上面插值多项式就称为n n次拉格朗日插值多项式次拉格朗日插值多项式。当 时， 分别为线性插值多项式线性插值多项式和二次插值二次插值多项式多项式。), 2 , 1 , 0( ,)()()()()()()(

14、11101110nkxxxxxxxxxxxxxxxxxxxxxlnkkkkkkknkkk)(xlknkkknxlyxL0)()(2 , 1n)(),(21xLxL1 拉格朗日（Lagrange）插值罗尔定理罗尔定理：如果函数 在 上连续， 内可导，并且 ，则至少存在一点 ，使得 。截断误差截断误差：如果 在区间 上连续， 在内存在， 是n+1个节点，则用 去近似所产生的截断误差为其中 且依赖于 ， 。)(xf,ba),(ba)()(bfaf),(ba0)(f)()(xfn,0nxx)()1(xfn),(0nxxnxxx10)(xLn)(xf)()!1()()()()(1)1(xnfxLxfxR

15、nnnn),(0nxxx)()()(101nnxxxxxxx1 拉格朗日（Lagrange）插值证明证明：由插值条件 可以得到 ，即n+1个节点是 的零点，所以设其中 是与 有关的待定函数。 为了求得 ，对区间 上异于 的任意一点 ，作辅助函数), 2 , 1 , 0( )(nkyxLkkn), 2 , 1 , 0( 0)(nkxRknnxxx10)(xRn)()()(1xxKxRnn)(xKx)(xK,0nxxkx), 2 , 1 , 0( nkxxk)()()()()(1txKtLtftFnn1 拉格朗日（Lagrange）插值将 看作是异于节点的一个固定点，则上式 满足（1） ，即 在上

16、有n+2个零点，分别为 ；（2）在 内具有n+1阶导数，并且有由罗尔定理，在 的两个零点之间至少存在一个 的零点，所以 在 内至少有n+1个互异的零点。x)(tF0)()()()(10nxFxFxFxF)(tF,0nxxnxxxx,10),(0nxx)!1()()()()1()1(nxKtftFnn)(tF)(tF)(tF),(0nxx1 拉格朗日（Lagrange）插值反复应用罗尔定理，最后可以得到 在 内至少有一个零点 ，即所以由此得到)()1(tFn),(0nxx),( , 0)!1)()()(0)1()1(nnnxxnxKfF)!1()()()1(nfxKn),( ),()!1()()

17、(101)1(xxxnfxRnnn1 拉格朗日（Lagrange）插值由于余项中含有因式如果插值点 偏离插值节点 比较远，则插值误差会比较大。如果插值点 位于插值区间内，插值过程称为内插内插，否则称为外推外推。根据余项定理，外推是不可靠的。 另外余项公式中有高阶导数项 ，就要求 足够光滑否则误差可能会比较大。 代数多项式是任意光滑的，原则上只适用于逼近光滑性好的函数。)()()(101nnxxxxxxxxxkx)()1(nf)(xf1 拉格朗日（Lagrange）插值例：例：已知 在点 的值由下表给出。试分别用线性插值与二次插值计算 的近似值，并进行误差估计。解：解：取 代入线性插值公式得xe

18、3 , 2 , 1x1 . 2exxe1230.3678794410.1353352830.0497870681 . 2, 3, 210 xxx126780462. 0 049787068. 02321 . 2135335283. 03231 . 2 ) 1 . 2(101001011yxxxxyxxxxL1 拉格朗日（Lagrange）插值取 代入二次插值公式得误差估计： 1 . 2, 3, 2, 1210 xxxx120165644. 0 049787068. 0)23)(13()21 . 2)(11 . 2( 135335283. 0)32)(12()31 . 2)(11 . 2(367

19、879441. 0)31)(21 ()31 . 2)(21 . 2( )()()()()()()(1202102210120120102102xxxxxxxxyxxxxxxxxyxxxxxxxxyxL00609009. 0)31 . 2)(21 . 2(! 2) 1 . 2(21eR00607001. 0)31 . 2)(21 . 2)(11 . 2(! 3) 1 . 2(12eR课后题课后题：1、当时 ， ，求 的二次插值多项式。2、已知函数 的观察值如下： 试求其拉格朗日插值多项式。1 拉格朗日（Lagrange）插值2 , 1, 1x 4 , 3, 0 xf xf xfy iixiy01

20、230123230-12 分段低次插值l 高次插值中的问题高次插值中的问题 一般来说，适当提高插值多项式的次数，会提高插值结果的准确程度。但是，提高插值多项式的次数，插值多项式会变得复杂，计算量加大。并且高次插值多项式往往具有数值不稳定的缺点，会产生高次插值不准确的龙格现象龙格现象。 所以当插值节点数n+1较大，特别是插值区间也较大时，通常不采用高次插值，而采用分段低次插值。常用的有分分段线性插值段线性插值和分段抛物线插值分段抛物线插值。分段低次插值的优点是公式简单，计算量小，且有较好的收敛性和稳定性，并且避免了计算机上作高次乘幂时常遇到的上溢和下溢。2 分段低次插值l 原因原因：(1)由拉格

21、朗日插值多项式余项，当差值节点增加时，的变化可能会很大，那么 可能很大；特别是当插值节点比较分散、插值区间较大时， 也比较大，这样就造成了近似时的截断误差较大；(2)当n增大时，拉格朗日插值多项式次数增加，计算量急剧增大，这样就加大了计算过程中的舍入误差。 )()1(xfn)(max)1(xfMnbxan)(1xn2 分段低次插值l 分段线性插值分段线性插值设在区间 上有节点 ，函数在上述节点处的函数值为 ，连接相邻两点 ，得到一条折线函数，如果满足：（1） 在区间 上连续；（2） ；（3） 在每个子区间 上是线性函数，则称折线函数 为分段线性插值函数分段线性插值函数。,babxxxxann1

22、10)(xf), 2 , 1 , 0( )(nixfyii) 1, 1 , 0( ),( ),(11niyxyxiiii)(x)(x)(x,ba), 2 , 1 , 0( )(niyxii,1iixx)(x2 分段低次插值 在每个子区间 上可以表示为从几何上讲，分段线性插值就是用一条过n+1个点 的折线来近似表示 。 在整个区间上用基函数来表示可以写为)(x,1iixx) 1, 2 , 1 , 0( ,)(1111niyxxxxyxxxxxiiiiiiii),( ,),(),(1100nnyxyxyx)(xf)(xbxaxlyxniii , )()(02 分段低次插值l 分段抛物线插值分段抛物

23、线插值分段抛物线插值就是把区间 分成若干个子区间，在每个子区间 上用抛物线去近似曲线，则用 表示分段抛物线插值函数 有下列性质：(1) 在区间 上是连续函数；(2) ；(3)在每个子区间上 ， 是次数不超过二次的多项式,ba) 1, 2 , 1(,11nixxii)(x11111111111111)()()()()()()(iiiiiiiiiiiiiiiiiiiiiyxxxxxxxxyxxxxxxxxyxxxxxxxxx)(x)(x)(x,ba), 2 , 1 , 0( )(niyxii,11iixx2 分段低次插值l 插值点选择插值点选择：选择插值点的原则是尽可能在插值点的邻近。公式中i的取

24、法归结为111111 , 11, 3 , 2, ,1, 3 , 2, , 1 , 1nkkkkkkkkxxnnkxxxxxxxknkxxxxxxxkxxi3 差商与牛顿（Newton）插值多项式 拉格朗日插值法有一个缺点缺点，当有了新的数据，插值节点增加时，插值多项式需要重新构造和计算，之前的计算结果无法继续利用。从构造算法的一般原则一般原则来说，应设法充分充分利用已经获得和计算的数据信息利用已经获得和计算的数据信息。为了克服拉格朗日插值法的缺点，介绍牛顿插值多项式。它使用比较灵活，增加插值节点时，只是在原来的基础上增加部分计算量，原来的计算结果仍可继续利用，节约了计算时间。3 差商与牛顿（N

25、ewton）插值多项式l 差商差商的定义 已知函数 在n+1互异节点 处的函数值分别为 ，称为 关于节点 的一阶差商一阶差商（平均变化率）。称为 关于节点 的二阶差商二阶差商。 )(xfnxxxx210)(,),(),(10nxfxfxfiiiiiixxxfxfxxf111)()(,)(xf1,iixxiiiiiiiiixxxxfxxfxxxf212121,21,iiixxx)(xf3 差商与牛顿（Newton）插值多项式 一般地，称为 关于节点 的 阶差商阶差商。当 时称 为 关于节点 的零阶差商零阶差商，记为由于所以即差商是微商的离散形式。 ikikiiikiiikiiixxxxxfxxx

26、fxxxf,11211)(xfkiiixxx,1k0k)(ixf)(xfixixfiiiixxixxxfxfxfii11)()(lim)(1,lim)(11iixxixxfxfii3 差商与牛顿（Newton）插值多项式l 差商性质：差商性质：（1）函数 关于节点 的k阶差商可以表示为函数值 的线性组合，即式中，如果 ，则其在 的导数为)(xfkxxx,10,10kxxxf)(,),(),(10kxfxfxfkjjkjkxxfxxxf0110)()(,)()()()(11101kjjjjjjjjkxxxxxxxxxxx)()()(101nnxxxxxxxixx )()()( )(lim)()(

27、lim)(11101111niiiiiiiinxxiinnxxinxxxxxxxxxxxxxxxxxxii3 差商与牛顿（Newton）插值多项式（2）差商与其所含节点的排列次序无关，即一般地，在k阶差商 中，任意调换节点的次序，其值不变。（3）设 在包含互异节点 的闭区间 上有n阶导数，则n阶差商与n阶导数之间有如下关系,11iiiixxfxxf,122121iiiiiiiiixxxfxxxfxxxf,10kxxxf)(xfnxxx,10,ba),( ,!)(,)(10banfxxxfnn3 差商与牛顿（Newton）插值多项式l 差商计算差商计算：利用插商的递推定义，差商的计算可列表计算，

28、如下表所示三阶差商二阶差商一阶差商 )( iixfx , , , )( , , )( , )( )( 321032132332102122101100 xxxxfxxxfxxfxfxxxxfxxfxfxxxfxfxxfx3 差商与牛顿（Newton）插值多项式l 牛顿插值多项式牛顿插值多项式由高等代数理论可知，任何一个不高于n次的多项式，都可以表示成函数的线性组合。所以满足插值条件的拉格朗日插值多项式又可以表示为式中 为待定系数。称这种n次插值多项式为牛顿（牛顿（NewtonNewton）插值多项式）插值多项式，记作 ，即 )()( ,),)( , 1110100nxxxxxxxxxxxx)(

29、xLn)()()()(110102010nnxxxxxxaxxxxaxxaa), 1 , 0(nkak)(xNn)()()()()(110102010nnnxxxxxxaxxxxaxxaaxN3 差商与牛顿（Newton）插值多项式由于 满足插值条件，即 ，所以由 ，得同样可以求出其他系数。)(xNn), 2 , 1 , 0( )(niyxNiin000)(yaxNn101101)()(yxxayxNn01011xxyya3 差商与牛顿（Newton）插值多项式有 是n+1个节点，对于一般情况，设 ，则由差商定义nxxx,10,baxxxi00000)()()()(,xxxfxfxxxfxfx

30、xf1100101010,xxxxfxxfxxxxfxxfxxxf221010210210210,xxxxxfxxxfxxxxxfxxxfxxxxfnnnnxxxxxfxxxfxxxxf,1010103 差商与牛顿（Newton）插值多项式得到)(,)()(000 xxxxfxfxf)(,110100 xxxxxfxxfxxf)(,221021010 xxxxxxfxxxfxxxf)(,1010110nnnnxxxxxxfxxxfxxxxf3 差商与牛顿（Newton）插值多项式将上式中第二式代入第一式式，得到式中可知 是满足插值条件的线性插值多项式线性插值多项式。而为线性插值的余项。)()(

31、 )(,)(,)()(1110100100 xRxNxxxxxxxfxxxxfxfxf)(,)()(01001xxxxfxfxN)(1xN)(,)(10101xxxxxxxfxR3 差商与牛顿（Newton）插值多项式同样，将第三式代入 得到式中是满足插值条件的二次插值多项式二次插值多项式。而为二次插值的余项。)()()(11xRxNxf)()( )()(, )(,)(,)()(22210210102100100 xRxNxxxxxxxxxxfxxxxxxxfxxxxfxfxf)(,)(,)()(1021001002xxxxxxxfxxxxfxfxN)()(,)(2102102xxxxxxxx

32、xxfxR3 差商与牛顿（Newton）插值多项式类似地将各式依次代入前式，最后可以得到其中为满足插值条件的n次插值多项式，通常称其为n n次牛顿插值次牛顿插值多项式多项式。)()( )()(, )()(, )(,)(,)()(101011010102100100 xRxNxxxxxxxxxxfxxxxxxxxxfxxxxxxxfxxxxfxfxfnnnnnn)()(, )(,)(,)()(11010102100100nnnxxxxxxxxxfxxxxxxxfxxxxfxfxN3 差商与牛顿（Newton）插值多项式与 相比较，有而为牛顿型插值余项牛顿型插值余项。 )()()()()(1101

33、02010nnnxxxxxxaxxxxaxxaaxN), 2 , 1 , 0( ,10nkxxxfakk)()(,)(1010nnnxxxxxxxxxxfxR3 差商与牛顿（Newton）插值多项式 由于满足插值条件的插值多项式存在且唯一，所以有如果 在 上有n+1阶导数，则有即)()(xLxNnn)(xf),(ba)()(xRxRnn)()!1()()(,)(1)1(10 xnfxxxxfxRnnnnn),( ,)!1()(,)1(0banfxxxfnn3 差商与牛顿（Newton）插值多项式容易看出，牛顿插值多项式具有递推性，即记 为具有节点 的牛顿插值多项式，则具有节点的牛顿插值多项式

34、为上式说明，增加一个节点 ，只要在 的基础上，增加计算即可。 )(xNk110,kxxxkxxx,10)(1xNk)()(,)()(101101kkkkxxxxxxxxxfxNxN1kx)(xNk)()(,10110kkxxxxxxxxxf3 差商与牛顿（Newton）插值多项式例例：已知一组观察数据如表，构造3次牛顿插值多项式。解解：首先计算差商iixiy012312340-5-63ixiy一阶差商二阶差商三阶差商102-5-53-6-12439513 差商与牛顿（Newton）插值多项式将计算得到的差商代入公式得到整理得到)3)(2)(1()2)(1(2) 1(50)(3xxxxxxxN3

35、4)(233xxxN625. 235 . 145 . 1)5 . 1 (233N例例：给出 的函数表求四次牛顿插值多项式，由此求 并估计误差。解解：选取最接近0.596的前5个节点，首先构造差商表3 差商与牛顿（Newton）插值多项式 xfix ixf0.400.550.650.800.901.050.41075 0.57815 0.69675 0.88811 1.02652 1.25382596. 0f3 差商与牛顿（Newton）插值多项式ix)(ixf一阶差商二阶差商三阶差商四阶差商五阶差商0.400.410750.550.578151.116000.650.696751.186000

36、.280000.800.888111.275730.358930.197330.901.026521.384100.433480.213000.031341.051.253821.515330.324930.228630.03126-0.00012 8 . 065. 055. 04 . 003134. 065. 055. 04 . 019733. 055. 04 . 028000. 04 . 011600. 141075. 04xxxxxxxxxxxN 当函数 的表达式未知或函数 的高阶导数比较复杂时，常用牛顿插值多项式余项但由于公式中的n+1阶差商 的值与 的值有关，因此不能准确计算 ，只能

37、对其做出一种估计。3 差商与牛顿（Newton）插值多项式63192. 0596. 0596. 04 Nf xf xf)()(,)(1010nnnxxxxxxxxxxfxR,10nxxxxf xf,10nxxxxf 当n+1阶差商变化不剧烈时，可用 近似代替 ，即采用此法计算 的误差，则有截断误差很小，可用忽略不计。3 差商与牛顿（Newton）插值多项式,110nnxxxxf,10nxxxxf)()(,)(10110nnnnxxxxxxxxxxfxR596. 04N9596. 0410623. 390. 055. 040. 005. 1 ,90. 0 ,80. 0 ,65. 0 ,55. 0

38、 ,40. 0596. 0 xxxxfR3 差商与牛顿（Newton）插值多项式例例：某处海洋不同深度水温如下表所示，试用牛顿插值公式求深度1000米处的水温，并估计误差。解解：计算差商x)(xf水深温度 （m）46671495014221634 （C）7.044.283.402.542.13kkx)(kxf,1kkxxf,21kkkxxxf,321kkkkxxxxf,4321kkkkkxxxxxf04667.0417144.28-0.0111329503.40-0.0037290.00001529314222.54-0.0018220.000002693-0.0000000131841634

39、2.13-0.007934-0.000000163-0.000000002750.893E-113 差商与牛顿（Newton）插值多项式用三次牛顿插值多项式 近似代替 ，得到)(3xN)(xf)950)(714)(466(10318. 1 )714)(466(10529. 1)466(10113. 104. 7)(8523xxxxxxxN)(331. 3)1000(3CN)(02878. 0)1000(1093. 8)1000(4123CR7、用拉格朗日插值和牛顿插值找经过点的三次插值多项式，并验证插值多项式的唯一性。8、利用函数表造出差商表，并利用牛顿插值公式计算 在处的近似值（计算取5位小

40、数）。3 差商与牛顿（Newton）插值多项式 10, 6,2, 3,2 , 0,1, 3x xfy 1.6151.6341.7021.8281.9212.414502.464592.652713.030353.34066 xf813. 1 ,682. 1x4 差分与等距节点插值公式 在实际应用中，常采用等距节点进行插值计算，这时插值公式可以进一步简化。由于插值节点等距分布，被插值函数的平均变化率与自变量的区间无关，差商可用差分代替。设被插值函数 在等距节点上的值 已知，其中 称为步长，则分别称为被插值函数 在 处以 为步长的向前差分向前差分和向后差分向后差分，符号 分别称为向前差分算子和向后

41、差分算子。 xfy niihxxi, 1 , 0 ,0 iixfy nxxxxhnii01,1iiiyyy1iiiyyy xfy ixh,4 差分与等距节点插值公式 高阶差分通过对低阶差分求差分来定义，如二阶差分为 阶差分为 iiiiiiiiiyyyyyyyyy111122211122iiiiiiiiiyyyyyyyyyn,111inininyyy111inininyyy4 差分与等距节点插值公式差分性质：1.差分可用函数值线性表示为式中组合表达式为2.差分与差商满足下述关系 inkinknkinninnininyyCyCyCyy112211!knknCknnkhkyxxxfkkk, 2 , 1 ,!,010nkhkyxxxfknkknnn, 2 , 1 ,!,14 差分与等距节点插值公