www.gusucode.com > EM算法源程序的代码，运行的结果和实验报告 > code10/源程序代码/junzhi.m

    function U = junzhi(Z,X,m,n,k)
temp1 = 0;
temp2 = 0;
for j=1:k
  for l=1:m
        for i=1:n
            temp1 = temp1 + Z(i,j)*X(i,l);
            temp2 = temp2 + Z(i,j);   
        end
        U(j,l) = temp1/temp2;
        temp1 = 0;
        temp2 = 0;
    end
end