Array of Matrices multiplication

Question

Hugo Hernández Hernández 2020-12-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/688369-array-of-matrices-multiplication

评论： Hugo Hernández Hernández 2020-12-9

Solving an exercise for orbital mechanics I found a problem while multiplying a 3x3 matrix by a 3x1 Matrix, actually these matrices have 1442 values so first one is 3x4326 and second one is 3x1442. First array is a vector of 3x3 matrices but I do not know how to separate each matrix and then multiply them by my 3x1 array of matrices. I have tried with a for loop, or by selecting the arrays or columns and then compute the operation but that did not work. The full code is very large and has other .m files involved so I am attaching the specific lines with the problem:

%% Position
R3_Thot = [cos(Tho_t), sin(Tho_t), z_0; -sin(Tho_t), cos(Tho_t), z_0; z_0, z_0, z_1]; 
%%% 3x3 array of matrices with 3x4326 length 
%%% R_p1 is an array of 3x1 matrices with a 3x1442 length
Rpos_Efix1 = R3_Thot*R_p1; %%% I want to perform this operation but I got errors like size mismatching 

I tried also with .* operator but does not work, and I am not sure how to solve or compute them using a for loop.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Hugo Hernández Hernández 2020-12-9

This question is almost the same: https://www.mathworks.com/matlabcentral/answers/496005-multiplying-a-3x3-matrix-where-values-change-with-each-iteration-by-a-3x1-matrix-with-values-that-it

请先登录，再进行评论。

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

Answer 1

Adam Danz 2020-12-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/688369-array-of-matrices-multiplication#answer_570199

编辑：Adam Danz 2020-12-9

在 MATLAB Online 中打开

One way to approach this is to break up the 3x4326 matrix into a 3x3x1442 array.

data = rand(3,4326);  % 3x4326 matrix
dataArray = reshape(data,3,3,size(data,2)/3);
size(dataArray)
ans = 1×3
           3           3        1442

Now you can multiply each 3x3 matrix by a 3x1 vector using pagemtimes() in Matlab r2020b or later.

vec = [1;2;3];
z = pagemtimes(dataArray, vec);
size(z)
ans = 1×3
           3           1        1442
●

For releases prior to Matlab r2020b, you can set up a loop

vec = [1;2;3];
z = nan(3,1,size(dataArray,3));
for i = 1:size(dataArray,3)
    z(:,:,i) = dataArray(:,:,i) * vec;
end
size(z)
ans = 1×3
           3           1        1442
●