将 3D 矩阵的每个切片与 2D 矩阵相乘答案

【问题标题】：Multiplication of each slice of a 3D matrix with 2D matrix [duplicate]将 3D 矩阵的每个切片与 2D 矩阵相乘
【发布时间】：2018-07-09 05:37:54
【问题描述】：

我尝试在我的工作中使用svd，它将返回如下向量。

A=rand(10,5);
b=rand(10,5,100);
[U, S , V]=svd(A);

我要做的是将b 的每个切片与U 相乘。对于一个特定的切片，此操作是有效的：

c=U'*b(:,:,1);

但是，当我尝试使用如下矢量化方法时，它会返回数组尺寸不匹配错误。

Utran=U.';
c = cellfun(@(x) x.*b,num2cell(Utran,[1 2]),'UniformOutput',false);

第一种方法我可能会使用循环，但如果我有一个大矩阵，它就没有效率。知道我的错误是什么吗？

【问题讨论】：

【解决方案1】：

以下是矢量化解决方案。

不确定它是否比使用循环更快；即使是，请注意它使用更多内存进行中间计算：

result = permute(sum(permute(conj(U), [1 3 4 2]).*b), [4 2 3 1]);

【讨论】：