201012137040孙谦第二次作业.doc

Lagrange.cpp #includestdio.h #includeiostream.h #includemath.h #define SIZEMAX 100 float Lagrange(double x[SIZEMAX],int n,int sum,double xx); void main() { int t; double x[SIZEMAX]; double y[SIZEMAX]; double l[SIZEMAX]; printf(请问你要输入一个几个插值节点？); cint; printf(请输入X的值); for(int i=0;it;i++) { cinx[i]; } printf(请输入Y的值); for(i=0;it;i++) { ciny[i]; } printf(请输入你要算的值); double xx; cinxx; double tt=0; for(i=0;it;i++) { l[i]=Lagrange(x,i,t,xx); tt=tt+l[i]*y[i]; } printf(输出的结果为:); couttt; } float Lagrange(double x[SIZEMAX],int n,int sum,double xx) { double a; a=1; double b; b=1; int i; double m; double t; for(i=0;isum;i++) { if(i!=n) { m=xx-x[i]; a=a*m; } if(i!=n) { m=x[n]-x[i]; b=b*m; } } t=a/b; return t; } Newton #includestdio.h #includeiostream.h #define SIZEMAX 100 void main() { int t; t=0; double X[SIZEMAX],Y[SIZEMAX]; printf(请输入要输入x的个数); cint; double ff[SIZEMAX][SIZEMAX]; for(int i=0;it;i++) for(int j=0;jt;j++) { ff[i][j]=0; } printf(请分别输入x,和y的值，先输入x的值); for(i=0;it;i++) { cinX[i]; } for(i=0;it;i++) { cinY[i]; ff[i][0]=Y[i]; } for(int k=1;kt;k++) { for(int r=k;rt;r++) { ff[r][k]=(ff[r][k-1]-ff[r-1][k-1])/(X[r]-X[r-k]); } } printf(请输入要求的x的值:); double xx,nn; cinxx; nn=0; for(i=0;it;i++) { for(int j=i-1;j=0;j--) { ff[i][0]=ff[i][0]*(xx-X[j]); } nn=nn+ff[i][0]; } coutendlendl; coutnn; } Hermite 插值 #includestdio.h #includeiostream.h void main() { printf(请输入x0,x1,y0,y1,y0`,y1`的值); double x,x0,x1,y0,y1,yy0,yy1; cinx0x1y0y1yy0yy1; printf(请输入你要算的x的值); cinx; double H; H=y0*(1-2*(x-x0)/(x0-x1))*(x-x1)*(x-x1)/(x0-x1)/(x0-x1)+y1*(1-2*(x-x1)/(x1-x0))*(x-x0)*(x-x0)/(x1-x0)/(x1-x0)+yy0*(x-x0)*(x-x1)*(x-x1)/(x0-x1)/(x0-x1)+yy1*(x-x1)*(x-x0)*(x-x0)/(x1-x0)/(x1-x0); coutendl; coutH; } 最小二乘法 #includestdio.h #includeiostream.h #define SIZEMAX 100 void main() { printf(请问你要输入多少组数据); int t; cint; double x[SIZEMAX],y[