matlab提取三维矩阵对角线的有效方法

Question

我想将3D矩阵的对角线（下面的Sigma）提取到另一个3D矩阵（下面的Sigma2）。

Sigma = repmat(magic(4),1,1,3);
Sigma2 = nan(1,4,3);
for i=1:3
    Sigma2(1,:,i) = diag(Sigma(:,:,i));
end

有更有效的方法吗？

Answer 1

你可以。如果您将Sigma重新整形为矩阵，则选择3D矩阵的对角线现在从矩阵中选择行。

Sigma3=reshape(Sigma,[],size(Sigma,3));
Selector=diag(true(size(Sigma,1),1));
Sigma2=Sigma3(Selector(:),:);
%Sigma2=permute(Sigma2,[3,1,2]) %optional last step to get a result with the same dimensions.