以八度为单位求矩阵的对角线元素

Question

给出方阵M，如何找到对角线上元素的总和？必须有比这更简单的方法：

sum(sum(diag(diag(M), 0)))

Answer 1

sum(diag(M))似乎工作正常。

1> M = reshape(1:9, 3, 3)
M =

   1   4   7
   2   5   8
   3   6   9

2> sum(diag(M))
ans =  15

Answer 2

实际上，我正在寻找的是踪迹：

1> M = reshape(1:9, 3, 3)
M =

   1   4   7
   2   5   8
   3   6   9

2> trace(M)
ans =  15