功能:插值求一阶导数

格式:
InterpolationDiff(x,y,h,w)
InterpolationDiff(x,y,h)
InterpolationDiff(x,y)

x:矩阵变量存储的单调递增序列
y:和x同维数的矩阵变量
h:求导时区间距离比例,默认为0.2
w:求导时两边值所占导数比重,默认为0.8

说明:
1、本函数求解在x个点处对应的导数值,所以返回的矩阵和x矩阵维度一样

原理:
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、则可以求出曲线上任意一点的值
5、最后计算x(i)两端分别相距比例为h点的斜率，最后根据权重加权w即可得到一阶导数

例子:
x=0:0.2:3;y=sin(x);dy=InterpolationDiff(x,y)