数值分析牛顿插值法

《数值分析牛顿插值法》由会员分享，可在线阅读，更多相关《数值分析牛顿插值法（27页珍藏版）》请在文档大全上搜索。

1、iiijjijiilxlbx11nnnnnnaaaaaaaaaA212222111211bAx ni, 3 ,2 2.2.2 Newton插值法插值法 2.2.3 等距节点插值公式等距节点插值公式华长生制作2)(xljnjiiijixxxx0)()(nj,2 , 1 ,0我们知道,Lagrange插值多项式的插值基函数为形式上太复杂,计算量很大,并且重复计算也很多由线性代数的知识可知,任何一个n次多项式都可以表示成, 1,0 xx ),)(10 xxxx)()(110nxxxxxx,共n+1个多项式的线性组合那么，是否可以将这n+1个多项式作为插值基函数呢？华长生制作3, 1,0 xx ),)

2、(10 xxxx)()(110nxxxxxx,显然，多项式组线性无关， 因此，可以作为插值基函数,ix设插值节点为( ),0,1,iiff xin函数值为1,2 , 1 ,0,1nixxhiiiiihhmaxnifxPii, 1 , 0,)(插值条件为)()()()()(110102010nnxxxxxxaxxxxaxxaaxP具有如下形式设插值多项式)(xP华长生制作4nifxPxPii, 1 , 0,)()(应满足插值条件000)(afxP有)()(011011xxaafxP00fa 01011xxffa)()()(12022021022xxxxaxxaafxP12010102022xxx

3、xffxxffa再继续下去待定系数的形式将更复杂为此引入差商和差分的概念华长生制作5一、差商(均差)定义1.nifxxfii, 1 , 0,)(处的函数值为在互异的节点设称)(,jixxffxxfjijiji)(,)(均差一阶差商关于节点为jixxxf)(,kjixxxxfxxfxxxfjkjikikji的二阶差商关于为kjixxxxf,)(依此类推华长生制作6,110kkiiiixxxxf阶差商的关于节点为kxxxxxfkkiiii,)(110,110kkxxxxf差商具有如下性质(请同学们自证)：且的线性组合表示可由函数值阶差商的,)(,),(),(,)()1(10110kkkxfxfxf

4、xxxxfkxf显然kkkkkiiiiiiiiixxxxxxfxxxf1210110,kkkkkxxxxxxfxxxf1210110,华长生制作7,110kkxxxxfkikiiiiiiixxxxxxxxxf0110)()()()(2) 差商具有对称性,即任意调换节点的次序,差商的值不变,210 xxxf,120 xxxf,012xxxf如,)()3(10)的区间存在时在包含节点当（kkxxxxf使得之间必存在一点在,10kxxx,10kxxxf用余项的相同证明!)()(kfk华长生制作8)()()()()()(4433221100 xfxxfxxfxxfxxfxxfxkk四阶差商三阶差商二阶

5、差商一阶差商差商的计算方法(表格法):,10 xxf,21xxf,32xxf,43xxf,210 xxxf,321xxxf,432xxxf,3210 xxxxf,4321xxxxf,410 xxxf规定函数值为零阶差商差商表Chashang.m华长生制作9xifxi fxi,xi+1fxi,xi+1,xi+2fxi,xi+1,xi+2 ,xi+200283275125621640208 1923827 493527125 9156125216 503419 10251949 14364991 105510 1261014 例1 求 f(xi)= x3在节点 x=0, 2, 3, 5, 6上的各

6、阶差商值解: 计算得如下表华长生制作10二、Newton基本插值公式)()()()()(110102010nnxxxxxxaxxxxaxxaaxP设插值多项式满足插值条件nifxPii, 1 , 0,)(则待定系数为00fa ,101xxfa ,2102xxxfa ,10nnxxxfa华长生制作11)()()()()(110102010nnnxxxxxxaxxxxaxxaaxNnkkjjkxxxxxff110100)(,称基本插值多项式次的关于节点为Newtonnxxfi)(定义3.)()()(xNxfxRnn)()!1()(1)1(xnfnn由插值多项式的唯一性,Newton基本插值公式的余

7、项为nkkkxxxxff1100)(,10)(kjjxx)(xk为k次多项式华长生制作12,10 xxxxfk,110 xxxxfk则视为一个节点若将,), 1 , 0( ,nixxi因此可得)(,)(000 xxxxffxf)(,(0110100 xxxxxxxfxxff)(,)(,10100100 xxxxxxxfxxxxffnjjnnkkjjkxxxxxxfxxxxxff010110100)(,)(,xxxxxxfxxxfkkk,11010)(,1010kkkxxxxxxfxxxf下面推导余项的另外一种形式华长生制作13)(xRn)()!1()(1)1(xnfnn)(,110 xxxxx

8、fnnnjjnnxxxxxxfxN010)(,)()()(xRxNnn因此)!1()()1(nfn,10nxxxxf!)()(kfk,10kxxxf)(xRk)(,1110 xxxxfkknk 一般Newton插值估计误差的重要公式另外华长生制作14. . 6 , 8 , 7 , 4 , 1)(,5 , 4 , 3 , 2 , 1 插值多项式求四次牛顿时设当iixfx练练习习kxkf(xk) 一阶差商 二阶差商 三阶差商 四阶差商012341234514786 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24)()4)(3)(2)(1()()3)(2)(1( 0)2)(1(3