关于matlab上的del2函数

Question

功能如何在边界上工作？ A[i,j] = 0.25*(L[i+1,j] + L[i,j+1] + L[i-1,j] + L[i,j-1]) - L[i,j]; - 这是内部元素，但是边界需要使用立方外推 - 我不知道如何做到这一点。

例如：

>> a = [
10 19 20 21
4 6 7 8
1 2 3 4
5 6 7 8
]
>> del2(a)
-4.25   0.5   2.5   4.5
 0.25     2   2.25  2.5
1.75      2       2     2
2.75   1.75 1.75  1.75

UPD：
当n> 3时发现算法 通过坐标行：
   g1（2：n-1）= a（3：n）-a（2：n-1）;
   g2（2：n-1）= a（2：n-1）-a（1：n-2）;
   g_x =（G1-G2）/ 2;
    g_x（1，:) = g_x（2，：）* 2 - g_x（3，:);
    g_x（n，:) = g_x（n-1，：）* 2 - g_x（n-2，:);
按坐标列：
g_y - ＆gt;像转换数组上的x坐标一样执行相同的步骤a ^ T;
之后结果= 0.5 *（g_y + g_x）！这就是它！