数值分析牛顿插值法.pptxVIP

数值分析牛顿插值法会计学第1页/共27页我们知道,Lagrange插值多项式的插值基函数为形式上太复杂,计算量很大,并且重复计算也很多由线性代数的知识可知,任何一个n次多项式都可以表示成共n+1个多项式的线性组合那么，是否可以将这n+1个多项式作为插值基函数呢？华长生制作第2页/共27页显然，多项式组因此，可以作为插值基函数线性无关，华长生制作第3页/共27页有再继续下去待定系数的形式将更复杂为此引入差商和差分的概念华长生制作第4页/共27页一、差商(均差)定义1.称依此类推华长生制作第5页/共27页显然差商具有如下性质(请同学们自证)：华长生制作第6页/共27页(2) 差商具有对称性,即任意调换节点的次序,差商的值不变如用余项的相同证明华长生制作第7页/共27页差商的计算方法(表格法):差商表Chashang.m规定函数值为零阶差商华长生制作第8页/共27页例1 求 f(xi)= x3在节点 x=0, 2, 3, 5, 6上的各阶差商值解: 计算得如下表xif[xi]f[xi,xi+1]f[xi,xi+1,xi+2]f[xi,xi+1,xi+2 ,xi+2]002832751256216华长生制作第9页/共27页二、Newton基本插值公式设插值多项式满足插值条件则待定系数为华长生制作第10页/共27页定义3.称为k次多项式由插值多项式的唯一性,Newton基本插值公式的余项为华长生制作第11页/共27页下面推导余项的另外一种形式因此可得华长生制作第12页/共27页因此一般Newton插值估计误差的重要公式另外华长生制作第13页/共27页kxkf(xk)一阶差商 二阶差商 三阶差商 四阶差商012341234514786 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24华长生制作第14页/共27页2.2.3 等距节点插值公式定义.华长生制作第15页/共27页依此类推可以证明如华长生制作第16页/共27页差分表华长生制作第17页/共27页在等距节点的前提下,差商与差分有如下关系华长生制作第18页/共27页依此类推华长生制作第19页/共27页1.Newton向前(差分)插值公式Newton插值基本公式为如果假设由差商与向前差分的关系华长生制作第20页/共27页则插值公式化为其余项化为华长生制作第21页/共27页称为Newton向前插值公式（又称为表初公式）插值余项为华长生制作第22页/共27页2.Newton向后(差分)插值公式如果假设根据向前差分和向后差分的关系可得Newton向后插值公式插值余项为华长生制作第23页/共27页 例 4 设x0=1.0,h=0.05,给出 在 处的函数值如表2-5的第3列，试用三次等距节点插值公式求f(1.01)和f(1.28)的近似值。表2-50 1.00 1.00000 0.024701 1.05 1.02470 0.02411 -0.000592 1.10 1.04881 0.02357 -0.00054 -0.000053 1.15 1.07238…… …… ……4 1.20 1.09544 0.02307 -0.00048 -0.000035 1.25 1.11803 0.02259 -0.000456 1.30 1.14017 0.02214 华长生制作 解 用Newton向前插值公式来计算f(1.01)的近似值。先构造与均差表相似的差分表，见表2-5得上半部分。由t=(x-x0)/h=0.2的得用Newton向后插值公式计算f(1.28)的近似值，可利用表2-5中的下半部分。由t=(x-x6)/h=-0.4，得　事实上，f(1.01)和f(1.28)的真值分别为11由此看出，计算结果是相当精确的。第24页/共27页 例 2.5 已知f(x)