Best way to calculate the determinants of a series of matrices?

31 次查看（过去 30 天）
Christian Schröder
评论： Christian Schröder 2019-9-5
I've got a series of time-dependent matrices which I'd like to calculate the determinants of. These matrices are stored as a three-dimensional array, where the first dimension indicates the period; in other words, the matrix at time t is given by
Gt = squeeze(G(t, :, :))
I'd now like a snippet of code which, being run, will ensure that
Delta(t) = det(squeeze(G(t, :, :)))
holds for all t. Of course I could do this with a loop, but I feel that there must be a more succinct, vectorized way of doing it. Sadly, MATLAB's det function itself is of no help. Is there something else I could use, or will I have to bite the proverbial bullet and use a loop after all?

请先登录，再回答此问题。

采纳的回答

Bruno Luong
Bruno Luong 2019-9-5
编辑：Bruno Luong 2019-9-5
I reverse the order and put the page in third dimension (avoid to use squeeze).
For small size, you can save CPU time by 4 fold using MultipleQR available on FEX
A=rand(3,3,1e5);
tic
n = size(A,1);
% FEX https://fr.mathworks.com/matlabcentral/fileexchange/68976-multipleqr
[Q,R] = MultipleQR(A);
R = reshape(R,n*n,[]);
d1 = (-1)^n * prod(R(1:n+1:end,:),1);
toc % Elapsed time is 0.087167 seconds.
tic
d2 = arrayfun(@(k) det(A(:,:,k)), 1:size(A,3));
toc % Elapsed time is 0.376470 seconds.
% Check correctness
norm(d1-d2)/norm(d2) % 4.2026e-16
MultipleQR will be less efficient for large n.

更多回答（3 个）

Alex Mcaulley
Alex Mcaulley 2019-9-5
delta = arrayfun(@(t) det(squeeze(G(t,:,:))),1:size(G,1));
Fabio Freschi
Fabio Freschi 2019-9-5
Not sure if it you can speedup your code, but a single line code to do the job is
Delta = arrayfun(@(i)det(squeeze(G(i,:,:))),1:size(G,1));
Jos (10584)
Jos (10584) 2019-9-5
Elaborating on the answers using arrayfun, you can avoid the multiple squeeze operations by permuting the dimension order first:
G = permute(G,[2 3 1]) ; % last dimension is running fastest
D = arrayfun(@(k) det(G(:,:,k)), 1:size(G,3)) % per Fabio and Alex
  1 个评论
Christian Schröder
Christian Schröder 2019-9-5
Good to know, I hadn't previously been aware that squeeze() could be avoided if the singleton dimension was in fact last --- or that you could use permute to rearrange the dimensions of an array. Still very much a MATLAB newbie, so thanks to you and everyone who pointed this out!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息

标签

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by