26.如何在Matlab中计算2D矩阵的局部积分？

Question

我有一个二维矩阵Z（x，y），我需要根据它来计算方程：

我使用代码执行此操作：

Sa        = nansum(nansum(abs(Z - Zbar))) / (nX*nY);

不适用：我使用nansum是因为Z(x,y)包含一些NaN值。

现在，我发现很难计算方程式：

我试图用以下代码分解方程式：

% compute the partial sums in the equation
for jj = 1: nY-1
    for ii = 1: nX-1
        Z_abt(jj,ii) = ( Z(jj,ii)- Z(jj,ii+1) )^2 + ( Z(jj,ii)- Z(jj+1,ii) )^2;       
    end
end

% compute Sdr
Sdr = nansum( nansum( sqrt(1+Z_abt) - 1) );

Sdr值始终应该是分数，但是我得到的值却大得离谱。

我只是想知道我是否真的需要一个循环来做到这一点，以及我该怎么做？