高斯混合模型的matlab实现

原文地址：http://www.crescentmoon.info/?p=463

高斯混合函数实现部分是基本上是转载的的pluskid大神文章里的里的代码，加了一点注释，并根据他给的方法二解决
covariance 矩阵 singular 的问题。

注意：

pluskid大神这里最后返回的是px，我觉得非常奇怪，因为PRML里对点做hard assignment时是根据后验概率来判别的。于是我在大神博客上问了一下，他的解释是最大似然和最大后验的区别，前者是挑x被各个模型产生的概率最大的那个，而后者加上了先验知识，各有道理。一句话就茅塞顿开，真大神也~

然后调用该函数对数据进行聚类，要是数据是二维的或三维的，顺便画个图。

下面是对鸢尾花数据集聚类的结果。选了部分维度画图

二维图：





三维图：