matlab离散数据求导数

需求：现实数据都是离散的，但是像希尔伯特变换求瞬时频率时，需要你对离散数据求导数。此时只能用差分近似代替求导。下面介绍matlab中实现差分(求导)的两种方法！

以一维离散的地震数据为例：
t是采样点时间序列，x是对应的振幅序列。
差分实现1(导数)：diff(x)./diff(t)
差分实现2(梯度)：gradient(x)./gradient(t)

区别：两者都能用，梯度用的是中心点差分，导数用的前后两点差分；所以从区间上看梯度用的范围比导数大一倍！所以梯度方式精度会更高一些！但是梯度法的边界可能会稍微出现点问题(不影响)。