第四章插值方法.pptVIP

第1页，共23页，星期日，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)§4.1多项式插值问题的一般提法第2页，共23页，星期日，2025年，2月5日§4.2拉格朗日(Lagrange)插值niyxPiin,...,0,)(==求n次多项式使得条件：无重合节点，即注:一次多项式插值---过两点直线。二次多项式插值---过三点抛物线。若不将多项式次数限制为n，则插值多项式不唯一。第3页，共23页，星期日，2025年，2月5日n=1已知x0,x1;y0,y1，求使得111001)(,)(yxPyxP==可见P1(x)是过(x0,y0)和(x1,y1)两点的直线。)()(0010101xxxxyyyxP---+=101xxxx--010xxxx--=y0+y1l0(x)l1(x)?==10)(iiiyxl称为拉格朗日插值基函数，满足条件li(xj)=?ij/*KroneckerDelta*/二.拉格朗日插值的基函数构造法第4页，共23页，星期日，2025年，2月5日n?1希望找到li(x)，i=0,…,n使得li(xj)=?ij；然后令?==niiinyxlxP0)()(，则显然有Pn(xi)=yi。拉格朗日插值多项式，常记为Ln(x)拉格朗日插值基函数与有关，而与无关节点f是n次多项式。li(x)每个li有n个根x0…xi-1…xn，第5页，共23页，星期日，2025年，2月5日三.插值余项设节点在[a,b]内存在,考察截断误差，且f满足条件,注：?通常不能确定?x,而是估计,?x?(a,b)将作为误差估计上限。?当f(x)为任一个次数?n的多项式时，,可知，即插值多项式对于次数?n的多项式是精确的。第6页，共23页，星期日，2025年，2月5日定义4.3.1差商(亦称均差)称为f关于xi和xj的1阶差商2阶差商§4.3差商与差分及其性质10111010],,...,[],...,,[],...,[+++--=kkkkkxxxxxfxxxfxxf(k+1)阶差商：第7页，共23页，星期日，2025年，2月5日事实上其中差商的值与xi的顺序无关！定义4.3.2差分当节点等距分布时:向前差分iiifff-=?+1一阶向前差分ikikikikffff1111)(-+--?-?=??=?k阶向前差分向后差分i?1iifff-=?一阶向后差分111---?-?=?ikikikfffk阶向后差分中心差分其中k阶中心差分第8页，共23页，星期日，2025年，2月5日§4.4牛顿插值公式牛顿插值公式：其中ai=f[x0,…,xi]Nn(x)Rn(x)第9页，共23页，星期日，2025年，2月5日第10页，共23页，星期日，2025年，2月5日注：?由唯一性可知Nn(x)?Ln(x)，只是算法不同，故其余项也相同，即?