我有两个非常大的列向量,A
和 B
, 大小 ax1
和 bx1
, 分别。我想构造一个向量 C
尺寸(b*a)x1
通过计算 A(i)*B(j)
对于每个 i
和 j
.为了说明我的意思:
clear
a=10^3;
b=10^3;
A=randn(a,1);
B=randn(b,1);
Ctemp=zeros(a,b);
for i=1:a
for j=1:b
Ctemp(i,j)=A(i)*B(j);
end
end
C=reshape(Ctemp, a*b,1);
问题:有没有更高效的方式获取C
哪个避免双重循环?我的实际a
和 b
大于10^3
.
最佳答案
这是一个可以从隐式(或显式)扩展中获益的简单数组乘法示例:
% Implicit (R2016b and newer):
C = A(:) .* B(:).'; % optionally surround by reshape( ... , [], 1);
% "Explicit" (R2007a and newer):
C = bsxfun( @times, A(:), B(:).' );
从那里开始,这只是 reshape 的问题,正如您已经在做的那样(D = C(:)
或 D = C(:).'
)。
关于matlab - 将一个向量的每个元素乘以另一个向量的每个元素,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/52689323/