Multiplying all 2D square matricies in a 3D matrix
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Consider a 3D matrix 'T' made from a concatenation of N square matrices of size m. Is there a easy way to multiply them all together (m1 x m2 x m3 x m4 ... x mN) into a single 2D matrix 'total' of size m? Similar to something like;
total = prod(T(:,:,:));
%This returns a 3D matrix of each 2D matrix of size (1,m)
In a for loop it may look like:
total = eye(m);
for i = 1;N
total = total*T(:,:,i);
end
Is it also possible to stop the multiplication at somepoint along the 3D matrix at point 1 < h < N;
total = prod(T(:,:,1:h))
I should also mention in my specific application the matrices are symbolic with both variables and numbers.
2 commentaires
James Tursa
le 21 Oct 2021
Modifié(e) : James Tursa
le 21 Oct 2021
Just use the for-loop. If for some reason you were really pressed for speed you could write a mex routine that could avoid all the intermediate data copies. But even this might only save you 10% - 20% execution time.
Réponse acceptée
Paul
le 21 Oct 2021
Modifié(e) : Paul
le 21 Oct 2021
Easier than putting that loop in a function?
function mprod = matprod(T,h)
% should do some input checking to verify that 1 <= h <= size(T,3)
mprod = T(:,:,1);
for ii = 2:h
mprod = mprod*T(:,:,ii);
end
end
3 commentaires
Paul
le 21 Oct 2021
I'm not aware of one, but maybe someone else has what you're loooking for. If they do, I really want to see it.
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Logical dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!