注解：
1.P5位置处的拉普拉斯算子的梯度相当于是：
P2-P5+P8-P5+P4-P5+P6-P5
最近卡在离散拉普拉斯算子del2这个函数上了，在网上查了好久，关于del2函数边缘点的处理公式都不对（通过与del2函数结果验证的），因为自己要用硬件加速算法，碰巧有拉帕拉斯算子，所以必须要知道每个点的具体运算。。。。
死磕了一个晚上，把del2函数的底层代码详细琢磨了一遍，终于搞定了。。。
矩阵中间点的计算公式，很容易，就是周边4个点的算数平均数减去这个点就ok了。。
重点是四个顶点以及边缘点的计算处理过程。。
设一个5阶矩阵为a=
A=del2(a),
左上角点A11=a11+a13+a31-(5a12+5a21+a14+a41)/4，
边缘点 A12=a32+(a11+a13-a42-5a22)/4，其他边缘角落点点类推。
可以去验证一下哦。。。。
matlab中del2这个函数结果降1/4，所以矩阵a实际的离散拉普拉斯应为4*del2(a)，即
A11=4a11+4a13+4a31-a14-a41-5a12-5a21，
A12=4a32+a11+a13-a42-5a22。
————————————————
版权声明：本文为CSDN博主「Ocean_VV」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/qq_41634276/article/details/80659218