help to multiply matrices in succession

Hi,
I initialized this objects:
C = vpa(zeros(4, 4, N));
A = vpa(zeros(4,4,N));
where C is:
for i = 1:N
C(:,:,i) = [(cos(theta(i))) (-sin(theta(i))).*(cos(alpha(i))) (sin(theta(i))).*(sin(alpha(i))) (a(i)).*(cos(theta(i)));
(sin(theta(i))) (cos(theta(i))).*(cos(alpha(i))) (-cos(theta(i))).*(sin(alpha(i))) (a(i)).*(sin(theta(i)));
0 sin(alpha(i)) cos(alpha(i)) d(i);
0 0 0 1];
end
Well, I need that
A(:,:,1)=C(:,:,1);
but elements:
A(:,:,2) ....... A(:,:,N) must be obtained as
A(:,:,2) = C(:,:,1) * C(:,:,2)
A(:,:,3) = C(:,:,1) * C(:,:,2) * C(:,:,3)
etc
....
A(:,:,N) = C(:,:,1) * C(:,:,2) * C(:,:,3) * ....* C(:,:,N)
some tips ?

 Respuesta aceptada

Matt Fig
Matt Fig el 21 de Nov. de 2012
Editada: Matt Fig el 21 de Nov. de 2012
% For example:
C = randi(10,2,2,3)
A = cumprod(C,3)
Or perhaps you meant:
A(:,:,1) = C(:,:,1);
for ii = 2:size(C,3)
A(:,:,ii) = A(:,:,ii-1)*C(:,:,ii);
end

3 comentarios

Matt Fig
Matt Fig el 21 de Nov. de 2012
Marco comments:
But cumprod is not for matrix product, is it ?
Matt Fig
Matt Fig el 21 de Nov. de 2012
Hello Marco, did you see the second option?
Marco
Marco el 21 de Nov. de 2012
Yes, thanks!!!

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Loops and Conditional Statements en Centro de ayuda y File Exchange.

Productos

Preguntada:

el 21 de Nov. de 2012

Community Treasure Hunt

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

Start Hunting!

Translated by