假设我有一个
AxBxC矩阵X和一个
BxD矩阵Y
是否有一种非循环方法,我可以通过该方法将每个
C AxB矩阵乘以Y ?
回答:
您可以使用函数
NUM2CELL将矩阵X分解为一个单元格数组,并使用
CELLFUN在单元格
之间进行操作,在
一行中完成此操作:
Z = cellfun(@(x) x*Y,num2cell(X,[1 2]),'UniformOutput',false); 结果Z是一个
1×C单元阵列,其中每个单元都包含一个
A×D矩阵。如果希望Z成为
A by D D by C矩阵,则可以使用
CAT函数:
Z = cat(3,Z{:});
注意:我的旧解决方案使用
MAT2CELL而不是
NUM2CELL ,它并不那么简洁:
[A,B,C] = size(X); Z = cellfun(@(x) x*Y,mat2cell(X,A,B,ones(1,C)),'UniformOutput',false);
更多&回答...