关于如何对离散点求导

sqqqqqq

已知excel中两列数据，分别1000+行，想通过求导的方式得出一个大致图形如下结果图形所示，我自己的程序是这样的 C=xlsread('C:\MATLAB7\work\Bdv3.xlsx');
x=C(:,1);y=C(:,2);z=C(:,3);
%%plot(x,y);
h=2.067833636*10^-15;i1=0.00357468;r1=1.1591;
m=pi*i1*i1*r1*r1;
dy=diff(y)./diff(x)
xd=x(1:max(size(x))-1);
yd=y(1:max(size(y))-1);
zd=z(1:max(size(z))-1);
I=diff(y)./diff(x);
H=-8*h*I.*zd.*xd.*yd/m;
plot(xd,H)

得出的图形却是自编图形这样的。
麻烦各位帮忙看看，谢谢！

sqqqqqq

又遇到一个问题，如何对一个响应函数求希尔伯特逆变换，请问各位希尔伯特逆变换的函数代码是什么？谢谢！