期望最大化，协方差步骤

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
  

每次迭代后得到的协方差值大约是600，当它应该从0.75稍微调整一下。

Answer 1

我怀疑像

这样的东西

  

平均（cGauss，：）

而不是

  

平均（cGauss）