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

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