matrix - 在 Octave 中求和矩阵的对角元素

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