gradient()是求数值梯度的函数。
用法
1、Fx = gradient(F,h):其中F为一维矩阵,返回F的一维数值梯度,h是F中相邻两点的间距,缺省时默认为1。
2、[Fx,Fy] = gradient(F,hx,hy):其中F为二维矩阵,返回F的二维数值梯度;hx表示x方向上相邻两点的距离,缺省时默认为1;hy表示y方向上相邻两点的距离,缺省时默认为1。
3、[Fx,Fy,Fz]=gradient(F,hx,hy,hz):其中F为三维矩阵,返回F的三维数值梯度,hx,hy,hz分别表示x、y、z三个方向上相邻两点的距离,缺省时默认为1。
理解
1、[Fx,Fy]=gradient(F,hx,hy),其中Fx为其水平方向上的梯度,Fy为其垂直方向上的梯度,Fx的第一列元素为原矩阵第二列与第一列元素之差除以hx,Fx的第二列元素为原矩阵第三列与第一列元素之差除以(2*hx),以此类推:Fx(i,j)=(F(i,j+1)-F(i,j-1))/(2*hx)。最后一列则为最后两列之差