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）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????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