牛顿插值法原理及应用汇总

牛顿插值法插值法是利用函数f(x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(x)的近似值。如果这特定函数是多项式，就称它为插值多项式。当插值节点增减时全部插值基函数均要随之变化，这在实际计算中很不方便。为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得到的一个公式：f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+...f[x0,...xn](x-x0)...(x-xn-1)+Rn(x)。插值函数插值函数的概念及相关性质白]定义：设连续函数y-f(x)在区间[a,b]上有定义，已知在n+1个互异的点x0,x1,・・・xn上取值分别为y0,y1,・・・yn(设aWx1Wx2WxnWb)。若在函数类中存在以简单函数P(x)，使得P(xi)=yi,则称P(x)为f(x)的插值函数.称x1,x2,・・・xn为插值节点，称[a,b]为插值区间。定理：n次代数插值问题的解存在且唯一。牛顿插值法C程序:鲁n育C£lII(T-trrTi岬凡b・il，・.I”*I-f母以日-T*f~『，..'.二完-*.'：/："FHjMf---程序框图#include<stdio.h>voidmain()(floatx[11],y[11][11],xx,temp,newton;inti,j,n;printf(〃Newton插值：\n请输入要运算的值：x二〃)；scanf(〃%f〃,&xx);printf("请输入插值的次数(n<11):n二〃)；scanf(〃%d〃,&n);printf("请输A%d组值：\n〃，n+1);for(i=0;i<n+1;i++)(printf(〃x%d二〃，i);scanf(〃%f〃,&x[i]);printf(〃y%d二〃，i);scanf(〃％f〃，&y[0][i]);}for(i=1;i<n+1;i++)for(j=i;j<n+1;j++)(if(i>1)y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-i]);elsey[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]-x[j-1]);printf("%f\n”,y[i][i]);}temp=1;newton=y[0][0];for(i=1;i<n+1;i++)(temp=temp*(xx-x[i-1]);newton=newton+y[i][i]*temp;}printf("求得的结果为：N(%.4f)二%9f\n",xx,newton);牛顿插值法Matlab程序functionf=Newton(x,y,x0)symst;if(length(x)==length(y))n=length(x);c(1：n)=0.0;elsedisp('x和y的维数不相等！');return;endf=y(1);y1=0;l=1;for(i=1:n-1)for(j=i+1:n)y1(j)=(y(j)-y(i))/(x(j)-x(i));endc(i)=y1(i+1);l=l*(t-x(i));f=f+c(i)*l;simplify(f);y=y1;if(i==n-1)if(nargin==3)f=subs(f,'t',x0);elsef=collect(f);%将插值多项式展开f=vpa(f,6);endend牛顿插值法摘要：值法利用函数f(x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(x)的近似值。如果这特定函数是多项式，就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化，这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。牛顿插值通过求各阶差商，递推得到的一个公式：f(x)=f[x0]+f[x0,x1](x-x0)+f[x0,x1,x2](x-x0)(x-x1)+・・・f[x0,・・.xn](x-x0)...(x-xn-1)+Rn(x)关键词：牛顿插值法流程图程序实现一、插值法的由来在许多实际问题及科学研究中，因素之间往往存在着函数关系，然而，这种关系经常很难有明显的解析表达，通常只是由观察与测试得到一些离散数值。有时，即使给出了解析表达式，却由于表达式过于复杂，不仅使用不便，而且不易于进行计算与理论分析。解决这类问题的方法有两种：一种是插值法，另一种是拟合法。插值法是一种古老的数学方法，它来自生产实践，早在一千多年前，我国科学家在研究历法上就应用了线性插值与二次插值，但它的基本理论却是在微积分产生之后才逐渐完善的，其应用也日益增多，特别是在计算机软件中，许多库函数，如等的计算实际上归结于它的逼近函数的计算。逼近函数一般为只含有算术运算的简单函数，如多项式、有理分式（即多项式的商）。在工程实际问题当中，我们也经常会碰到诸如此类的函数值计算问题。被计算的函数有时不容易直接计算，如表达式过于复杂或者只能通过某种手段获取该函数在某些点处的函数值信息或者导数值信息等。因此，我们希望能用一个“简单函数”逼近被计算函数，然后用该简单函数的函数值近似替代被计算函数的函数值。这种方法就叫插值逼近或者插值法。逐次线性插值法优点是能够最有效地计算任何给定点的函数值，而不需要写出各步用到的插值多项式的表达式。但如果解决某个问题时需要插值多项式的表达式，那么，它的这个优点就成了它的缺点了。能不能根据插值条件构造一个插值多项式，它既有具体的表达式，又很容易用它计算任何点的函数值呢？牛顿插值法能作到这一点。二、牛顿插值法的概念牛顿插值多项式的表达式设N=c+c(X-X)+c(X-X)(%-X)+...C(X-X)(%-%)...(%-X)n010201n01n-1问题是如何根据插值条件NX)=y.…ni”i=0,1,2...n来计算待定系数C0，C1，C2…气？

n(x)=y=f(x)TOC\o"1-5"\h\zn0Q0知,由知C=y=f(x)0Q知,由知N(x)=y=f(x)n111c+c(x-x)=y01101因而c=1-y0=f(x?—f气)嘴x,x]1x—xx—xJ0,11010,其中^x,x'称为函数f(x)在x，x点的一阶商。/、I011勺如J人1-^01八、、v)117。由N(x2)=y=f(x2)知C+C(x—x)+C(x—x)(x—x)=y

011022021因而y-y-f[x,x](x-x)C=01202/(xLx)(x-x)y乙f[x,x](x-x)4100__1——20_(x-x)(x-x)()2021y-y4f[x,x](x-x)-f[x,x](x-x)2101100120(x-x)(x-x)2021y-yxy+-f[x,x]x—x01,-x)20f[x,x]-f[x,x]八「2]J'1]Of[x,x,x](x-x)八0,1,2」20其中f[xxx2]称为函数f(x)在x0,x,x2点的二阶差商。实际上，它是一阶差商的差商。一般地，如果已知一阶差商f[x.-1,x]，f[x,七]，那么就可以计算阶差商f[x,x,x]=f[x，*1]-f[x1,x].-1..+1x_x.+1.-1类似于上述过程不断地推导下去，可得C=f[x，x2,¥一f[>x，x2]Of[x,x,x,x],(A—x)0123f[x,x,x:x]-f[x,x,x,x]C=1—2—34—1―2—卜Of[x,x,x,x,x],(x-x)01234■2„f[x,x.••x]-f[x,x,x.••x]C=f[1,2——n]v八v°,1,2——n^Of[x,x,x•••x],4,vv、■2n-1012n其中，f[x,x,x,x,x],f[x,x,x,x,x,x],分别称为函数01234012345f[x,x,x,x],0123f(x)在相应点处的三阶差商，四阶差商和n阶差商。实际上，xf(x)=c00xf(x)f[x,x]=c11011xf(x)f[x,x]f[x,x,x]=c22120122xf(x)f[x,x]f[x},x,x]f[x,x,x,x]=c332312301233xf(x)f[x,x]f[x,x,x]f[x},x,x,x]f[x,x,x,x,x]=c44342341234012344..♦♦♦♦..♦♦♦♦..♦♦♦♦nc0,c1,c2..（的计算可通过以下简易地构造函数的差商来完成。按上述方式构造插值多项式的方法叫做牛顿插值法。根据插值多项式的惟一性知，其截断误差与拉格朗日插值法相同，即：Rn=志"但也可以表示成差商形式。这是因为以七孔土…如|为节点的多项式N⑴=N⑴+f[x,x.••x]n3）

n+1n01n+1n+1从而f（七尸NJ%尸七（%n+1）+f[X0,x...xjnJ'/于是Nn⑴的截断误差可表为TOC\o"1-5"\h\zR(x)=f[x,x，…x,x]n(x)n01n+1n+1顺便指出，因为牛顿插值多项式具有性质:N(x)=N(x)+f[x,x.••x](x-x)(x-x).•.(x-x)nn-101n12n-1所以，类似于逐次线性插值法，也可以把上述和式中的第二项f[x,x•••x](x-x)(x-x)•••(x-x)01n12n-1看成是估计Nn-1(x)的一种实用误差估计式。与差商概念密切联系的另一个概念是差分，它是指在等距节点上函数值的差。所谓等距节点，是指对给定的常数h(称为步长)，节点x=x0+论，(.=0，1，2…n)称f(x+1)-f(x)VAfk为x.处的一阶向前差分;/(x)-/(x)◊△fii-1i为^.处的一阶向后差分;f(x+V)-f(x.l《)。8f为*・处的中心差分。一阶差分的差分称为二阶差分，Af-AfOA2f称为x处的二阶向前差分。般地，m阶向前和向后差分可定义如下:Amf=Am-1f—Am-1f,Vmf=Vm—1f—Vm—1f=2,3…三、牛顿插值法的实现1、【算法】步骤1：输入节点(xj,yj),精度&，计值点xx,f0Tp,1TT,1Ti；步骤2：对k=1,2,,i依次计算k阶均差f[xi-k,xi-k+1,…,xi]=(f[xi-k+1,…,xi]-f[xi-k,…,xi])/(xi-xi-k)步骤3：(1)、若If[x1,…,xi]-f[x0,…,xi-1]l<&，则p为最终结果Ni-1(x),余项Ri-1=f[x0,…,xi](xx-xi-1)T。(2)、否则(xx-xi-1)*TtT,p+f[x0,…,xi]*Ttp,转步骤4。步骤4：若i<n，则i+1ti，转步骤2；否则终止。2、【流程图】3、【程序清单】#include"stdio.h”#definen4//牛顿插值的次数voidmain()(floata[n+1][n+2]={0},s=0,t=1,x;inti,j;printf("请输入xi及yi的值〃要求先输入xi再输入yi然后输入下一组\n");for(i=0;i<n+1;i++)for(j=0;j<2;j++)scanf("%f",&a[i][j]);for(j=1;j<n+2;j++)//计算各阶均差for(i=j;i<n+1;i++)a[i][j+1]=(a[i][j]-a[i-1][j])/(a[i][0]-a[i-j][0]);printf(-输出xi，yi及各阶均差\n");for(i=0;i<n+1;i++){for(j=0;j<n+2;j++)printf("%6.5f",a[i][j]);printf("\n");printf(-输出牛顿插值表达式\n");printf("N%d(x)=",n);for(i=0;i<n+1;i++)(printf("%6.5f",a[i][i+1]);for(j=0;j<i;j++)printf("(x-%3.2f)",a[j][0]);if(i==n)break;printf("+");}printf("\n");printf("输入插值点x=");scanf("%f",&x);for(i=0;i<n+1;i++)//计算插值点的近似值(for(j=0;j<i;j++)t*=(x-a[j][0]);s+=a[i][i+1]*t;}printf("N%d(%4.3f)=%6.5f\n",n,x,s);}4.【程序实现】；'F^^DebugXl.eKe"国输/3及敏的道"要求先和g再输.侦然后输"一组日.4®H.550.G5^.83输出xL由一嘶丽H切.SS。藻涂).65seeJ.U13U0U"勃0®谕出牛顿插循表达式.14<«?>=0.41075+1.ilG(J0Cx-0.46>+0.2S»aa<x-0.4a><x-a.555+0.19732<x-0.40Xx-9.S5>Cx-0.65>+B.03132<x0.4O>Cx-0.55><x-6.65XxB.E0>输入插蕴点翼部■卯6|M4<3,596>=0.62998Pressanykeytocontinue0.41075B,S7S15a.G9G7E日.丽成41.UZ65Z兄及各阶均差«.41075虬讷IS6.696751-026520.丽，IE1.11&0O1.1S6»O1.27