ps: 这是核心的代码,混合高斯模型的概率求解公式如下:
for(k=0; k
nsubclasses; k++)
{
if(k==0) maxlike = subll[k];
if(subll[k]>maxlike) maxlike = subll[k];
}
/* Sum weighted subclass likelihoods */
subsum = 0;
for(k=0; k
nsubclasses; k++)
subsum += exp( subll[k]-maxlike )*C->SubSig[k].pi;
ll[m] = log(subsum) + maxlike;
}
http://cobweb.ecn.purdue.edu/~bouman/software/cluster/