BSpline-三次样条.docVIP

BSpline-三次样条 // BSpline 插值 /* 三次B－Spline 好用，可根据您自己的X 与 Y 进行任意插值。插值与matlab差距很小，保证用一次您就会喜欢！绝对，吐血推荐！ */ #includemath.h #includestdio.h #includeconio.h #define PI 3.1415926 #define total_length 50 #define orginal_length 20 int N = 4; int step = 1; int inter_points = 3; int total_number = (N/step - 1) * (inter_points + 1) + 1; // /***/***/***/ double spl(double *x,double *y,int n,double t, double eps); float Spline_Interpolation(int n,float y,float x[]); void main() { FILE *fp; double *x, *y, *output; x = new double[13]; y = new double[13]; output = new double[13]; int i, j, n; double eps, t; eps = 0.000001; printf(input the number of dots :); // eg. 0 1 2 3 we should input: 4 scanf(%d,n); for(i=0;in;i++) { x[i] = i; y[i] = sin(i); } /* x[0] = 0; y[0] = 2.3026; x[1] = 1; y[1] = 2.3979; x[2] = 2; y[2] = 2.4849; x[3] = 3; y[3] = 2.5649; */ t = x[0] ; for (j=0;jn,t n;j++) { output[j]= spl(x,y,n,t,eps); //spl(double *x,double *y,int n,double *t,double eps) t = t + 0.25; } fp=fopen(function.txt,w); // function for(j=0;jtotal_number;j++) fprintf(fp,%10.4f\n,output[j]); // output fclose(fp); } /*====================================================== // 函数名：spl // 功能描述：光滑插值 // 输入参数：x 指向存放n个结点的数据的数组的指针 // y 指向存放n个结点的函数值的数组的指针 // n 结点个数 // t 指定的插值点,eps 绝对值小于此值的数据认为是0 // 返回值： 在指定插值点的函数近似值 =========================================================*/ // #include stdio.h // #include math.h // double spl(x,y,n,t,eps) // double *x,*y,t,eps; // int n; double spl(double *x,double *y,int n,double t,double eps) { int k; double z,d[5],p[4],yy[2]; double tmp1,tmp2; if((x==NULL)||(y==NULL)) /* 检测输入指针是否为空*/ { printf(Pointer is Null\n); return(0.0); } if(n1) /* 没有提供结点，返回0.0*/ return(0.0); if(n==1) /* 只有一个结点，返回此函数值*/ return(y[0]); if(n==2)