功能:采用分段两点插值求解(一维插值)

格式:
InterpolationPoly(x,y,x1)

x:矩阵变量存储已知的单调递增序列
y:和x同维数的矩阵变量
x1:需要插值的矩阵变量

说明:
1、本函数求解在x1对应的y1值，

原理:
1、假设在序列x、y中，每一段x(i)与x(i+1)中满足方程 y=a+b*x+c*x^2+d*x^3,
2、根据x(i)与x(i+1)中两个端点值与端点处的一阶导数，可以求出未知系数a、b、c、d
3、这样每一段x(i)与x(i+1)均有一个对应的a、b、c、d
4、则可以求出曲线上任意一点的值

例子:
x=0:0.2:3;y=sin(x);x1=0:0.02:4;y1=InterpolationPoly(x,y);plot(x,y,x1,y1);//可以绘图查看插值效果