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
  

各反復後に取得する共分散値は、0.75 からわずかに調整する必要がある場合、600 のオーダーです。

Answer 1

私は次のようなことを疑っています

mean(cGauss,:)

それ以外の

mean(cガウス)