八度中矩阵的对角元素求和

Question

给定方矩阵M,如何找到对角线上元素的总和？必须有一个比这更简单的方法:

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

Answer 1

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

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

   1   4   7
   2   5   8
   3   6   9

2> trace(M)
ans =  15

Answer 2

只是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