Matlab中网格点的矢量化评估

Question

我需要以下代表性代码部分中的“for”循环尽可能高效地运行。代码中的平均函数充当我自己函数的代表占位符。

feature1

我已经对代码进行了矢量化，并通过构建一个矩阵来提高其速度约2.5倍，其中mc被复制用于每个网格点。我的实现产生了一个非常大的矩阵（3 x 22500000），填充了重复的数据。通过将矩阵转换为单精度，我减轻了这种方法的内存损失，但似乎应该有一种更有效的方法来做我想要的，避免复制如此多的数据。

Answer 1

您可以bsxfun使用reshapes -

{{1}}