matlab - 期望最大化，协方差步骤

Question

我正在实现一个我知道是正确的算法（它已被提供），但我认为我没有正确编码它。这个方程的 Matlab 代码是什么？

在哪里：

• r 是一个 2x400 矩阵
• x 是一个 1x400 向量

• mu 是一个 1x2 向量

  for (cGauss = 1:2)
  sumR = sum(r(cGauss,:));
  cov(1,cGauss) = sum(r(cGauss,:).*((x - mean(cGauss))*(x - mean(cGauss))'))/sumR;
  end
  

每次迭代后我得到的协方差值在 60​​0 左右，此时它应该从 0.75 略微向上调整。

Answer 1

我怀疑类似

均值（cGauss,:)

代替

均值（cGauss）