训练GMM参数:
options=statset('MaxIter',1000);
gmm = gmdistribution.fit(feature,k,'Covtype','Diagonal','Regularize',1e-10,'Options',options);
k:聚类的总类别数
covtype:设置协方差矩阵的类型
regularize:防止协方差矩阵出现奇异矩阵,在协方差矩阵的对角上加一个很小的值
options:GMM模型的各种手动参数
用训练好的GMM模型实现聚类:
result = posterior(gmm, feature) # 计算后验概率,即每个特征向量属于每一类聚类的概率矩阵