数值分析-计算方法上机题集解析.docVIP

实验课程： 数值分析 专 业： 信息与计算科学 班 级： 学 号： 18 姓 名： 付强 中北大学理学院 实验一 函数差值方法 【】【】的n+1个节点值（j=0，1,…n）。试用拉格朗日公式求其差值多项式或分段二次拉格朗日插值多项式。 数据如下： x 0.4 0.55 0.65 0.80 0.95 1.05 y 0.41075 0.57815 0.69675 0.90 1.00 1.25382 球五次拉格朗日插值多项式或分段三次插值多形式，并计算的值。（提示：结果为） 【】【】 其中， 因此，实现此公式的关键步骤在于用C语言循环语句，实现此差值公式，然后运用数组输入具体实验数字，求的所需要的结果。 实验C语言代码如下： #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 #include stdlib.h #includemalloc.h float average(int n,float *x) {int i; float av; av=0; for(i=0;in;i++) av+=*(x+i); av=av/n; return(av); } //平方和 float spfh(int n,float *x) {int i; float a,b; a=0; for(i=0;in;i++) a+=(*(x+i))*(*(x+i)); return(a); } //和平方 float shpf(int n,float *x) {int i; float a,b; a=0; for(i=0;in;i++) a=a+*(x+i); b=a*