向量化矩阵中不同对角线的总和

[poster]

我想向量化以下MATLAB代码。我认为它必须很简单，但是我仍然感到困惑。

r = some constant less than m or n [m,n] = size(C); S = zeros(mr,nr); for i=1:m-r+1 for j=1:n-r+1 S(i,j) = sum(diag(C(i:i+r-1,j:j+r-1))); end end 该代码从另一个得分表C计算用于动态编程算法的得分表S。
对角线求和是为所有可能的片段（大小为r）生成用于生成C的单个数据片段的分数。

预先感谢您的任何回答！抱歉，这应该很明显...


注意
事实证明，内置的conv2比convnfft更快，因为我的eye（r）很小（5