matlab求图像特征点的二阶偏导数，结构hessian矩阵

matlab求图像特征点的二阶偏导数，构造hessian矩阵
matlab怎么求图像特征点的二阶偏导数，构造hessian矩阵？我想通过hessian矩阵的特征值获得特征点的曲率，新手不知道怎么求，各位大神请留步
------解决思路----------------------
一阶导数X方向模板 [1 0 -1]，其他类似
------解决思路----------------------
在最优化的工具箱里，专门有个hessian函数。