將gmdistribution.fit與數據一起繪制

Question

我有一個數據（x，y）的直方圖，它是兩個高斯混合而成的（簡單的情況，即在一維中）。 為了適應它，我用

obj = gmdistribution.fit（data，2）

這給了我

obj =

高斯混合物在1維中具有2個組分的分布組分1：混合比例：0.499387均值：-0.4789

成分2：混合比例：0.500613均值：-0.4786

現在，如何將其與原始直方圖一起繪制？ 我知道這兩個鏈接，但是它們對我沒有太大幫助：

http://www.mathworks.com/matlabcentral/newsreader/view_thread/317868

http://www.mathworks.de/de/help/stats/gmdistribution.fit.html

任何明確的幫助都將受到高度贊賞！

Answer 1

我想這可能有點晚了，但是它已在matlabs 網站上進行了更新。 以防萬一

MU1 = [1 2];
SIGMA1 = [2 0; 0 .5];
MU2 = [-3 -5];
SIGMA2 = [1 0; 0 1];
X = [mvnrnd(MU1,SIGMA1,1000);mvnrnd(MU2,SIGMA2,1000)];

scatter(X(:,1),X(:,2),10,'.')
hold on

options = statset('Display','final');
obj = gmdistribution.fit(X,2,'Options',options);
10 iterations, log-likelihood = -7046.78

h = ezcontour(@(x,y)pdf(obj,[x y]),[-8 6],[-8 6]);

您可以使用scatter3對3d進行相同操作。