Multiplication of matrix and first dimension of a 3-D array without for loops

Question

Valentina Baccetti 2023-1-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1893220-multiplication-of-matrix-and-first-dimension-of-a-3-d-array-without-for-loops

回答： Matt J 2023-1-13

Hi,

I'm fairly new to Matlab but I was wondering if there is a way to multiply a 2-D matrix by the first dimension of a 3-D array without using for loops, and obtain another 3-D array.

The matrix, let's call it A, would be something like A(n,n), while the 3-D would be something like B(n,l,k). What I would like to do is somehow multiple A(n,n)*B(n,l,k) for every value of l and k, without going through two for loops.

Example below

A = rand(4,4);
B = rand(4,10,5);
for i = 1:10
    for j = 1:5
    C(:,i,j) = A(:,:)*B(:,i,j);
    end
end

Thanks

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

Rik 2023-1-13

Why exactly do you want to remove the loop? And are you aware you're doing a matrix multiplication instead of an element-wise multiplication?

请先登录，再进行评论。

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

Answer 1

Bruno Luong 2023-1-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1893220-multiplication-of-matrix-and-first-dimension-of-a-3-d-array-without-for-loops#answer_1147740

编辑：Bruno Luong 2023-1-13

在 MATLAB Online 中打开

A = rand(4,4);
B = rand(4,10,5);
% orginal code
for i = 1:10
    for j = 1:5
    C(:,i,j) = A(:,:)*B(:,i,j);
    end
end
% one loop code
Coneloop = zeros(size(A,1),size(B,2),size(B,3));
for j = 1:size(B,3)
    Coneloop(:,:,j) = A*B(:,:,j);
end
% no loop code
Cnoloop = pagemtimes(A,B);
% Check matching
norm(Coneloop(:)-C(:),'Inf')
ans = 4.4409e-16
norm(Cnoloop(:)-C(:),'Inf')
ans = 4.4409e-16
norm(Cnoloop(:)-Coneloop(:),'Inf')
ans = 0