bspline-三次样条（bspline - 三次样条）.docVIP

bspline-三次样条（bspline - 三次样条） / / bspline 插值 / * 三次b－spline 好用, 可根据您自己的x 与 y 进行任意插值.插值与matlab差距很小, 保证用一次您就会喜欢! 绝对, 吐血推荐! * / # include math.h # include stdio.h # include conio.h # define pi 3.1415926 # define total _ length 50 # define original _ length 20 int n = 4; int a = 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 x, float y, []); 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; epc = 0.000001; printf (input the number of dots:); / / eg. 0 1 2 3 we should input: 4 scanf (% d, n); for (i = 0; i n; i + +) { x [i] = i; y [in] = a (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 ]； 对于（j＝0；j n，t；n；j +） { 输出= SPL（x，y，n，t，EPS）；/ SPL（双××，双×Y，int，双T，双EPS） T＝0.25； } FP = fopen（“功能。txt”，“W”）；//函数 为（J = 0；J total_number；j++） fprintf（FP，“% 10.4f \n”，输出[ J ]）；//输出 关闭文件（FP）； } / * ====================================================== / /函数名：SPL / /功能描述：光滑插值 / /输入参数：X指向存放N个结点的数据的数组的指针 / / Y指向存放N个结点的函数值的数组的指针 / / N结点个数 / / T指定的插值点，EPS绝对值小于此值的数据认为是0 / /返回值：在指定插值点的函数近似值 ========================================================= */ / / #包括“stdio. h” / / #包括“数学。” /双声压级（x，y，n，t，EPS） /双××，y，t，EPS； n； 双声压级（双××，双×Y，int N，双T，双EPS） { int k； 双Z，d [ 5 ]，p [ 4 ]，YY [ 2 ]； 双TMP1和TMP2； if (x = = null) | | (y = = null) 检测输入指针是否为空) / * * / { printf (pointer is zero \ n); return (0.0). } if (n = 1) / * 没有提供结点, 返回0.0 * / return (0.0). if (n = = 1) / * 只有一个结点, 返回此函数值 * / return (y [0]); if (n = = 1) / * 只有两个结点, 使用线性插值 * / { z = (y [0] * (t - x [1]) - y [1] * (t, x (0))) / (x [0], x [1]). return (z). } k = 0; while ((x [k] (t) (k n) 寻找合适的区间) / * * / c + +; k = k 1; if (k 1) k = 0; if (k = = (n - 1) k = n - 2. d [2] = (y [k + 1] - y [k])