Я делаю негибкую регистрацию с использованием B-сплайнов, и я хотел бы добавить термин регуляризации энергии изгиба в мою функцию стоимости. В основном этот член регуляризации состоит в реализации этого уравнения: 
, где T - матрица смещений, а V - размер объема. Я сделал реализацию в 2D, но это не кажется мне правильным. Так как правильно реализовать это в 2D? Спасибо за поддержку.
%Calculating the penalty terms
dx=[1;0;-1];
dy=dx.';
dIdxdy=convn( convn(grid,dx,'same') ,dy,'same');
dIdx2=convn(convn(grid, dx, 'same'), dx, 'same');
dIdy2=convn(convn(grid, dx, 'same'), dx, 'same');
C = (1/I)*sum(sum(sum(dIdx2.^2+dIdy2.^2+2*dIdxdy.^2)));