Как функция работает на границах?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 -> Делать те же шаги, что и по координате x на преобразованном массиве a ^ T;После этого результата = 0,5 * (g_y + g_x)!ЭТО ЭТО!