estep.m

一个在matlab环境下编写的采用expectation maximization方法计算高斯混合模型的程序。

function p = eStep(x, y, mu, sigma2, px)
temp = 0;
[ans number] = size(px);
for i = 1:number
    temp = temp + pyx(y, i, mu, sigma2)*px(i);
end
p = pyx(y,x,mu,sigma2)*px(x)/temp;