手动创建样条函数分段多项式

Question

考虑在值数组中采样的函数y（x），由数组x和y表示。如果我有另一个x值x0，我可以使用样条线

来评估y（x0）

y0 = spline(x,y,x0);

现在，我也可以写

pp = spline(x,y);
y0 =  ppval(pp,x0);

我的问题：如果我已经有系数和x矩阵my_coefs（size(my_coefs) = [length(y),4]）和x，我该如何创建分段多项式{ {1}} My_pp和pp.coefs = my_coefs？

Answer 1

好的，没有＆＃34;样条对象＆＃34;，而是一个分段多项式对象。因此，如果my_coefs是通过断点my_x获得的，则所需的代码是

my_spline = mkpp(my_x,my_coefs);
y0 = ppval(my_spline, x0);

如果这里的尺寸令人眼花缭乱，那么

my_coefs是4 * n

my_x是n

y0是N

x0是N