数值分析插值法.pptVIP

Rn(x)Nn(x)ωn+1(x)多项式Nn(x)显然满足插值条件,即Nn(xj)=f(xj),(j=1,…,n),且次数不超过n,由唯一性定理它就是前述的Ln(x),其系数为Nn(x)称为牛顿均差插值多项式,它比拉格朗日插值多项式计算量省,且便于程序设计.第31页，共52页，星期日，2025年，2月5日注：?由唯一性可知Nn(x)?Ln(x)，只是算法不同，故其余项也相同，即?实际计算过程为f(x0)f(x1)f(x2)…f(xn?1)f(xn)f[x0,x1]f[x1,x2]…………f[xn?1,xn]f[x0,x1,x2]…………f[xn?2,xn?1,xn]f[x0,…,xn]f(xn+1)f[xn,xn+1]f[xn?1,xn,xn+1]f[x1,…,xn+1]f[x0,…,xn+1]均差计算可列均差表如下：第32页，共52页，星期日，2025年，2月5日例1依据如下函数值表建立不超过3次的拉格朗日插值多项式及牛顿插值多项式Nn(x),并验证插值多项式的唯一性.解:(1)拉格朗日插值多项式Ln(x).插值基函数xk0124f(xk)19233拉格朗日插值多项式为:第33页，共52页，星期日，2025年，2月5日xkf(xk)一阶均差二阶均差三阶均差0119822314343-10-8(2)牛顿插值多项式Nn(x).建立如下差商表牛顿插值多项式为:(3)唯一性验证.通过比较牛顿插值多项式和拉格朗日插值多项式,知:Nn(x)=Ln(x)这一事实与插值多项式的唯一性一致.第34页，共52页，星期日，2025年，2月5日2第35页，共52页，星期日，2025年，2月5日第36页，共52页，星期日，2025年，2月5日第37页，共52页，星期日，2025年，2月5日第1页，共52页，星期日，2025年，2月5日当精确函数y=f(x)非常复杂或未知时，在一系列节点x0…xn处测得函数值y0=f(x0),…yn=f(xn)，由此构造一个简单易算的近似函数g(x)?f(x)，满足条件g(xi)=f(xi)(i=0,…n)。这里的g(x)称为f(x)的插值函数。x0x1x2x3x4xg(x)?f(x)第2页，共52页，星期日，2025年，2月5日根据实际需要，可以用各种不同的函数来近似原来的函数。最常用的插值函数是…?多项式：代数多项式最简单，计算其值只需用到加、减乘运算，且积分和微分都很方便；所以常用它来近似表示表格函数(或复杂函数),这样的插值方法叫做代数插值法，简称插值法。第3页，共52页，星期日，2025年，2月5日§1拉格朗日多项式niyxPiin,...,0,)(==求n次多项式使得条件：无重合节点，即n=1已知x0,x1;y0,y1，求使得111001)(,)(yxPyxP==可见P1(x)是过(x0,y0)和(x1,y1)两点的直线。)(1xP101xxxx--010xxxx--=y0+y11.1线性插值两点式)()(0010101xxxxyyyxP---+=点斜式)(001010xxxxxxy---+=(())ff第4页，共52页，星期日，2025年，2月5日1.2二次插值n=2已知x0,x1,x2;y0,y1,y2,求使得002,)(yxP112)(yxP==222)(yxP=,为求P2(x),将三点代入其表达式,即可得到三个方程式,从而联立方程组解出系数a0,a1,a2即可: