最近卡在離散拉普拉斯運算元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。