插值：三次样条插值与分段插值11组.docVIP

作业名称：插值：三次样条插值与分段插值 完成日期：2012年4月7日 【小组成员：(自评成绩选择95\90\85\80\75\70\65\60\50即可)】 班级 学号 姓名 主要工作说明 自评成绩 090711 2009071128 李涛 整理排版 95 090711 2009071101 卫胜迪 画流程图 090711 2009071115 王洪涛 编写算法程序 一、计算公式(5%) 二、算法(5%) 三、程序流程图(30%) 四、完整的程序及简要的注释(40%) 拉格朗日插值法： #include stdio.h #includeconio.h #include malloc.h float Lagrange(float *x,float *y,float xx,int n)//拉格朗日插值函算法 { int i,j; float *a,yy=0.0;//a作为临时变量，记录拉格朗日插值多项式 a=(float*)malloc(n*sizeof(float));//分配存储区 for(i=0;in;i++) { a[i]=y[i]; for(j=0;jn;j++) if(j!=i) a[i]=a[i]*((xx-x[j])/(x[i]-x[j])); yy+=a[i]; } free(a); return yy; } int main() { int i,n; float x[15],y[15],xx[29],yy[29]; printf(Intput n:);//输入已知节点个数 scanf(%d,n); for(i=0;in;i++) //输入x { printf(x[%d]:,i); scanf(%f,x[i]); } printf(\n); for(i=0;in;i++) //输入y { printf(y[%d]:,i); scanf(%f,y[i]); } printf(\n); for(i=0;i29;i++) {printf(Input xx[%d]:,i); scanf(%f,xx[i]); yy[i]=Lagrange(x,y,xx[i],n);}//调用拉格朗日函数 for(i=0;i29;i++) printf(x=%f,y=%f\n,xx[i],yy[i]); getch(); } 牛顿插值法： #include stdio.h float Newton(float x1[],float y1[],float xx,int n) //牛顿插值算法 {float x[30],y[30]; int i,j; float yy,f[30]; for(i=0;in;i++) { x[i]=x1[i];y[i]=y1[i]; } for(i=0;in;i++) { f[0]=y[i]; for(j=0;ji;j++) f[j+1]=(f[j]-y[j])/(x[i]-x[j]); y[j]=f[i]; } yy=y[n]; for(i=n-1;i=0;i--) yy=yy*(xx-x[i])+y[i]; return(yy); } int main() { float x[30],y[30]; float xx[29],yy[29]; int i,n; printf(Input n:\n); scanf(%d,n); //输入已知节点个数 printf(Input x[%d]:\n,n); for(i=0;in;i++) scanf(%f,x[i]); // 输入x printf(Input y[%d]:\n,n); for(i=0;in;i++) scanf(%f,y[i]); //输入y for(i=0;i29;i++) { printf(Input xx[%d]:,i); scanf(%f,xx[i]); yy[i]=Newton(x,y,xx[i],n);}//调用牛顿插值函数 for(i=0;i29;i++) printf(x=%lf,y=%lf\n,xx[i],yy[i]); } 三次样条插值法： #include stdio.h int main() { int i,k; floa