数值分析1007024114王嘉乐.doc

实验课程： 数值分析 专 业： 信息与计算科学 班 级： 学 号： 1007024114 姓 名： 王嘉乐 中北大学理学院 实验一函数插值方法 【】【】的的节点值。试用Lagrange公式求其插值多项式或分段二次Lagrange插值多项式。 数据如下： 0.4 0.55 0.65 0.80 0.95 1.05 0.41075 0.57815 0.69675 0.90 1.00 1.25382 求五次Lagrange多项式或分段三次插值多项式，并计算，的值。（提示：结果为，） 【】【】 其中 容易求得： 程序 #include math.h #include stdio.h #include stdlib.h #includemalloc.h float lag(int n,float *x,float *y,float t) {int i,j; float p,s,la; p=1;la=0; for(i=0;in;i++) {for(j=0;jn;j++) {if(j!=i) p=p*(t-*(x+j))/(*(x+i)-*(x+j)); } s=*(y+i)*p; p=1; la=la+s; } return(la); } void main() { int m,n,k,i; float *x,*y,la,t; printf(请输入给定点的数量n); scanf(%d,n); x=(float*)calloc(n,sizeof(float)); if(x==NULL) {printf(内存分配失败); exit(1); } y=(float*)calloc(n,sizeof(float)); if(y==NULL) {printf(内存分配失败); exit(1); } x[0]=0.4;x[1]=0.55;x[2]=0.65;x[3]=0.80;x[4]=0.95;x[5]=1.05; y[0]=0.41075; y[1]=0.57815; y[2]=0.69675; y[3]=0.90; y[4]=1.00;y[5]=1.25382; for(i=0;in;i++) { printf(x[%d]=%f ,i,*(x+i)); printf(y[%d]=%f\n,i,*(y+i)); } for(i=0;;i++) {printf(请输入t的值(退出则输入0)); scanf(%f,t); if(t==0) {printf(确认退出请输入0，计算f(0)请按输入任意数字值); scanf(%d,m); if(m==0)break; } printf(请输入需要使用的点的个数); scanf(%d,k); la=lag(k,x,y,t); printf(插值结果是:\n); printf(f(%f)=%f\n,t,la); } free(x); free(y); }【】 【】【】【】与时间的拟合曲线。 0 5 10 15 20 25 30 35 40 45 50 55 0 1.27 2.16 2.86 3.44 3.87 4.15 4.37 4.51 4.58 4.02 4.64 1 1 1 1 1 1 1 1 1 1 1 【】【】 选取基函数，，则得到 ，，， ，，于是得方程组 程序实现如下： #include math.h #include stdio.h void spir(int n,int m,double x[],double y[],double a[], double dt[]) { int i,j,k; double alpha,p,q,g,w,beta,d1,d2,s[20],t[20],b[20]; for (i=0; i=m-1;i++) a[i]=0.0; if (mn) m=n; if (m20) m=20; b[0]=1.0;d1=n; alpha=0.0;q=0.0; for (i=0;i=n-1;i++) {alpha=alpha+x[i];q=q+y[i];} q=q/d1; //q0 alpha=alpha/d1; //alpha(0) a[0]=q*b[0];