2D matrix multiply with 3D array

2 次查看（过去 30 天）
Qiu Xu
Qiu Xu 2021-10-23
评论： Qiu Xu 2021-10-24
Hi,
I have two 2D matrix u and v, and a 3D array B. I want to perform the following calculations, for example
u=rand(200,300);
v=rand(200,300);
B=rand(60,60,300)
c=0;
for j=1:300
u1=0;
v1=0;
for k=1:200
u1=u1+u(k,j)*B(:,:,k);
v1=v1+v(k,j)*B(:,:,k);
end
c=c+u1.*v1;
end
Due to two for loops, it is very slowly for the larger matrix u, v and B. Are there some good methods to compute c without using any for loop?.
Thanks a lot.
  1 个评论
Matt J
Matt J 2021-10-23
Are you sure the size of B shouldn't be 60x60x200? If not, your code example ignores the final 100 pages of B.

请先登录，再进行评论。

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

采纳的回答

Matt J
Matt J 2021-10-23
编辑：Matt J 2021-10-23
[m,n,p]=size(B);
Br=reshape(B,[],p);
Br=Br(:,1:200); %might not be necessary
U1=Br*u;
V1=Br*v;
c=reshape( sum(U1.*V1,2) , m,n);
  1 个评论
Qiu Xu
Qiu Xu 2021-10-24
Thanks, I have make a mistake that the size of B should be 60x60x200. Your method is very good.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by