C语言曲线拟合.docVIP

C语言曲线拟合 #include stdio.h #include stdlib.h #include malloc.h #include math.h Smooth(double *,double *,double *,int,int, double *,double *,double *); void main() { int i,n,m; double *x,*y,*a,dt1,dt2,dt3,b; n = 20; m = 6; b = 0; /*分别为x,y,a分配存贮空间*/ x = (double *)calloc(n,sizeof(double)); if(x == NULL) { printf( 内存分配失败\n ); exit (0); } y = (double *)calloc(n,sizeof(double)); if(y == NULL) { printf( 内存分配失败\n ); exit (0); } a = (double *)calloc(n,sizeof(double)); if(a == NULL) { printf( 内存分配失败\n ); exit (0); } for(i=1;i =n;i++) { x[i-1]=b+(i-1)*0.1; /*每隔0.1取一个点，这样连续取n个点*/ y[i-1]=x[i-1]-exp(-x[i-1]); //y[i-1] = x[i-1]*x[i-1]*x[i-1]; //y[i-1] = x[i-1] ; /*计算x[i-1]点对应的y值作为拟合已知值*/ } Smooth(x,y,a,n,m,dt1,dt2,dt3); /*调用拟合函数*/ for(i=1;i =m;i++) printf( a[%d] = %.10f\n ,(i-1),a[i-1]); printf( 拟合多项式与数据点偏差的平方和为：\n ); printf( %.10e\n ,dt1); printf( 拟合多项式与数据点偏差的绝对值之和为：\n ); printf( %.10e\n ,dt2); printf( 拟合多项式与数据点偏差的绝对值最大值为：\n ); printf( %.10e\n ,dt3); free(x); /*释放存储空间*/ free(y); /*释放存储空间*/ free(a); /*释放存储空间*/ } Smooth(double *x,double *y,double *a,int n,int m,double *dt1,double *dt2,double *dt3 ) // //double *x; /*实型一维数组，输入参数，存放节点的xi值*/ // //double *y; /*实型一维数组，输入参数，存放节点的yi值*/ // //double *a;