牛顿插值法原理及应用.pdf

实用文档 牛顿插值法 插值法 是利用函数 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) 。 插值函数 插值函数的概念及相关性质 [1] 定义： 设连续函数 y-f(x) 在区间 [a,b] 上有定义，已知在 n+1 个互异的点 x0,x1, …xn 上取值分别为 y0,y1, …yn （设a≤ x1 ≤x2 ……≤ xn ≤b) 。若在函 数类中存在以简单函数 P(x) ，使得 P(xi)=yi, 则称 P(x) 为 f(x) 的插值函数 . 称 x1,x2, …xn 为插值节点，称 [a,b] 为插值区间。 定理： n 次代数插值问题的解存在且唯一 。 实用文档 牛顿插值法 C程序 程序框图 #includestdio.h void main() { float x[11],y[11][11],xx,temp,newton; int i,j,n; printf(Newton 插值 :\n 请输入要运算的值 :x=); scanf(%f,xx); printf( 请输入插值的次数 (n11):n=); scanf(%d,n); printf( 请输入 %d组值 :\n,n+1); for(i=0;in+1;i++) { printf(x%d=,i); scanf(%f,x[i]); printf(y%d=,i); scanf(%f,y[0][i]); 实用文档 } for(i=1;in+1;i++) for(j=i;jn+1;j++) { if(i1) y[i][j]=(y[i-1][j]-y[i-1][j-1])/(x[j]- x[j-i]); else y[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;in+1;i++) { temp=temp*(xx-x[i-1]); newton=newton+y[i][i]*temp; } printf( 求得的结果为 :N(%.4f)=%9f\n,xx,newton); 牛顿插值法 Matlab 程序 fun