数值方法上机作业.doc

PAGE PAGE 1 第二章 插值法 1.通用程序编制 （1）n次拉格朗日插值计算公式 #includestdio.h void main() { int n,m,i,j,k; float a[100],b[100],c,y[100],x[100]; /*输入插值多项式次数*/ scanf(%d,n); /*输入所求未知量个数*/ scanf(%d,m); /*输入插值节点坐标*/ for(i=0;i=n;i++) { scanf(%f,a[i]); scanf(%f,b[i]); } /*输入所求点x坐标*/ for(i=0;im;i++) { scanf(%f,x[i]); } /*计算*/ for(j=0;jm;j++) { y[j]=0; for(k=0;k=n;k++) { c=1; for(i=0;i=n;i++) { if(k!=i) c=c*(x[j]-a[i])/(a[k]-a[i]); } c=c*b[k]; y[j]=y[j]+c; } } /*输出计算结果*/ for(i=0;im;i++) { printf(%f\n,y[i]); } getchar(); } （2）n次牛顿向前插值计算公式 #includestdio.h void main() { int n,m,i,j; float x,h,a,b[100],c[100],y[100],t[100]; /*输入插值多项式次数*/ scanf(%d,n); /*输入所求未知量个数*/ scanf(%d,m); /*输入步长*/ scanf(%f,h); /*输入xo点坐标*/ scanf(%f,x); /*输入插值节点坐标*/ for(i=0;i=n;i++) { scanf(%f,b[i]); } /*输入所求点t值*/ for(i=0;im;i++) { scanf(%f,t[i]); } /*造向前差分表*/ for(i=0;in;i++) { for(j=i;jn;j++) { c[j+1]=b[j+1]-b[j]; } for(j=i+1;j=n;j++) { b[j]=c[j]; } } /*计算*/ for(i=0;im;i++) { a=1; y[i]=0; for(j=0;j=n;j++) { y[i]=y[i]+a*b[j]; a=a*(t[i]-j)/(j+1); } } /*输出计算结果*/ for(i=0;im;i++) { printf(x=%f y=%f\n,x+t[i]*h,y[i]); } getchar(); } （3）n次牛顿向后插值公式 #includestdio.h void main() { int n,m,i,j; float x,h,a,b[100],c[100],y[100],t[100]; /*输入插值多项式次数*/ scanf(%d,n); /*输入所求未知量个数*/ scanf(%d,m); /*输入步长*/ scanf(%f,h); /*输入Xn点坐标*/ scanf(%f,x); /*倒着输入插值节点坐标*/ for(i=0;i=n;i++) { scanf(%f,b[i]); } /*输入所求点t值*/ for(i=0;im;i++) { scanf(%f,t[i]); } /*造向后差分表*/ for(i=0;in;i++) { for(j=i;jn;j++) { c[j+1]=b[j]-b[j+1]; } for(j=i+1;j=n;j++) { b[j]=c[j]; } } /*计算*/ for(i=0;im;i++) { a=1; y[i]=0; for(j=0;j=n;j++) { y[i]=y[i]+a*b[j]; a=a*(t[i]+j)/(j+1); } } /*输出计算结果*/ for(i=0;im;i++) { printf(x=%f y=%f\n,x+t[i]*h,y[i]); } getchar(); } 2.计