EM EM算法Matlab实现。最大期望（EM）算法是在概率（probabilistic）模型中寻找参数最大似然估计或者最大后验估计的算法

EM算法Matlab实现。最大期望（EM）算法是在概率（probabilistic）模型中寻找参数最大似然估计或者最大后验估计的算法，其中概率模型依赖于无法观测的隐藏变量（Latent Variable）

function [gammallm]=em1k（rk）
load （‘data.mat‘）
x=dataset1;
L=size（x）;
l=L（1）;
zc=zeros（r40）;llm=-inf;
for zzz=1:r
sigma=zeros（22k）;
mu=zeros（k2）;
pik=ones（k1）/k;
[musigma]=km1s（2）;
for zz=1:40
p=zeros（lk）;
for ii=1:k
    p（:ii）=mvnpdf（xmu（ii:）sigma（::ii））;
end
gamma=zeros（lk）;
for ii=1:l
    su=0;
    for jj=1:k
        su=pik（jj）*p（iijj）+su;
    end
    for jj=1:k
        gamma（iijj）=pik（jj）*p（iijj）/su;
    end
end
nk=zeros（k1）;
for ii=1:k
    nk（ii）=sum（gamma（:ii））;
end
pik=nk/l;
for ii=1:k
    mu（ii1）=sum（gamma（:ii）.*x（:1））/nk（ii）;
    mu（ii2）=sum（gamma（:ii）.*x（:2））/nk（ii）;
end
ssig=zeros（22l）;
for ii=1:k
    for jj=1:l
        ssig（::jj）=gamma（jjii）.*（x（jj:）-mu（ii:））‘*（x（jj:）-mu（ii:））;
    end
    sigma（11ii）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      46897  2011-04-17 10:00  DATA.mat

     文件       1456  2011-04-17 19:37  em1k.m

     文件       1629  2011-04-17 19:36  em1r.m

     文件       1419  2011-04-17 19:36  em2k.m

     文件       1592  2011-04-17 19:36  em2r.m

     文件        378  2011-08-24 21:10  example.m

----------- ---------  ---------- -----  ----

                53371                    6