Function 插值后的近似函数

如何近似插值点的函数？我使用自然三次样条插值点，对于n=500点，如下所示：

t=[0; 3; 6; 9]
z=[0; 6; 6; 0]
plot(t,z,'ro')
ti=linspace(0,9,500)

zn = natcubicspline(t,z,ti)
yn = line(ti,zn)

对于这些n个多点（对于大n点），是否有近似函数？或者，是否有方法将插值点视为一个函数，即求zn向量的梯度？因为zn是常数的向量，所以这不一定有用

更新：特别是，我的数据似乎形成了一个二次多项式，因此我继续使用以下Matlab函数来拟合我的数据：

p = polyfit(transpose(ti),zn,2)

它产生了一个二次多项式的系数估计。它确实适合数据，但误差值很高，我必须将这个系数向量乘以向量[1 z^2]才能得到正确的多项式。有什么办法可以简化这个过程吗