在 octave 或 matlab 中矢量化列操作

Question

如何矢量化以下代码？它基本上是计算每一列的平均值。

mu(1) = sum(X(:,1))/C
mu(2) = sum(X(:,2))/C

还有这个（标准化每个元素，每列有不同的平均值和标准）：（X 是 47x2.mu，sigma 都是 1x2）

X_norm(:,1) = (X(:,1)-mu(1))/sigma(1)
X_norm(:,2) = (X(:,2)-mu(2))/sigma(2)

Answer 1

很简单：

mu = sum(X) ./ C

sum 默认沿第一个维度（在列上）运行。

---

编辑：

对于问题的第二部分：

X_norm = bsxfun(@rdivide, bsxfun(@minus, X, mu), sigma)

这类似于 repmat 的操作，但没有内存开销。

Answer 2

你甚至可以使用mu = mean(X)。