第4章 插值和拟合.pptVIP

Note：由存在性证明过程，均差与点列 的顺序无关。根据唯一性，N(x)就是L(x)。 例4.2.3 用牛顿插值法构造例4.2.1中的2次插值多项式。 解用函数的3个点(0,1),(-1,5),(2,-1)作均差表 0 1 | -4 1 -1 5 | -2 2 -1 | N2(x)=1+(-4)(x?0)+1?(x?0)(x+1)=x2?3x+1 正如预料之中，结果与用Lagrange插值法相同。但却无计算Lagrange插值基函数之麻烦。第4章插值和拟合?多项式插值第30页，共68页，星期日，2025年，2月5日 Newton插值法的更方便之处是当增加新的插值节点或其他插值条件时,可以利用原有结果，构造满足新插值条件的多项式。 例4.2.4设已知f(x)的如下值： f(-1)=-2，f(0)=-1，f(1)=0，f?(0)=0 构造不超过3次的多项式p3(x)，使满足插值条件 p3(-1)=-2，p3(0)=-1，p3(1)=0，p?3(0)=0 解首先作过f(x)的3个点(-1,-2),(0,-1),(1,0)的均差表 这3个点显然在一条直线上，因此2阶均差为0。 -1 -2 | 1 0 0 -1 | 1 1 0 | N2(x)=-2+(x+1)=x?1第4章插值和拟合?多项式插值第31页，共68页，星期日，2025年，2月5日 把点的顺序变换一下，结果也一样： 1 0 | 1 0 0 -1 | 1 -1 -2 | N2(x)=0+(x?1) 令p3(x)=N2(x)+c3(x?x0)(x?x1)(x?x2)，由p?3(0)=0确定c3, p?3(x)=1+c3(x?x1)(x?x2)+c3(x?x0)(x?x2)+c3(x?x0)(x?x1) p?3(0)=1+c3x1x2+c3x0x2+c3x0x1=1+0+(-c3)+0=0，c3=1 所以 p3(x)=(x?1)+(x?x0)(x?x1)(x?x2) =(x?1)+(x+1)x(x?1)=x3?1第4章插值和拟合?多项式插值第32页，共68页，星期日，2025年，2月5日 Newton插值的余项 根据均差的定义，把x看作[a,b]上的一点，则有 f[x,x0]={f(x)?f(x0)}/(x?x0) f(x)=f(x0)+(x?x0)f[x,x0] f[x,x0,x1]={f[x,x0]?f[x0,x1]}/(x?x1) f[x,x0]=f[x0,x1]+(x?x1)f[x,x0,x1] ?? f[x,x0,?,xn?2]=f[x0,x1,?,xn?1]+(x?xn?1)f[x,x0,?,xn?1] f[x,x0,?,xn?1]=f[x0,x1,?,xn]+(x?xn)f[x,x0,?,xn] 后一式代入前一式，即得第4章插值和拟合?多项式插值第33页，共68页，星期日，2025年，2月5日 f(x)=f(x0)+(x?x0)f[x,x0] =f(x0)+(x?x0){f[x0,x1]+(x?x1)f[x,x0,x1]} =f(x0)+(x?x0)f[x0,x1]+(x?x0)(x?x1)f[x0,x1,x2]+? +(x?x0)(x?x1)?(x?xn?1)f[x0,x1,?,xn] +(x?x0)(x?x1)?(x?xn?1)(x?xn)f[x,x0,x1,?,xn] =Nn(x)+Rn(x) 其中 Rn(x)=f(x)?Nn(x)=wn+1(x)f[x,x0,x1,?,xn] (H) wn+1(x)=(x?x0)(x?x1)?(x?xn)第4章插值和拟合?多项式插值第34页，共68页，星期日，2025年，2月5日 由于 f(x)=Nn(x)+wn+1(x)f[x,x0,x1,?,xn] 根据插值多项式的唯一性，Nn(x)=Ln(x)，因此两种插值多项式的余项也相等。事实上，利用Roll’sTheorem可证明n+1阶均差与n+1阶导数的关系： Lagrange形式的插值余项要求被插函数的n+1阶导数存在，而Newton形式的插值余项不要求被插函数的导数存在，因而可用于f(x)是由函数表给出或不可导的场合。式(H)更具有一般性。第4章插值和拟合?多项式插值第35页，共68页，星期日，2025年，2月5